Evolution of eukaryotic translation elongation and termination factors: variations of evolutionary rate and genetic code deviations.

Translation is carried out by the ribosome and several associated protein factors through three consecutive steps: initiation, elongation, and termination. Termination remains the least understood of them, partly because of the nonuniversality of the factors involved. To get some insights on the evolution of eukaryotic translation termination, we have compared the phylogeny of the release factors eRF1 and eRF3 to that of the elongation factors EF-1alpha and EF-2, with special focus on ciliates. Our results show that these four translation proteins have experienced different modes of evolution. This is especially evident for the EF-1alpha, EF-2, and eRF1 ciliate sequences. Ciliates appear as monophyletic in the EF-2 phylogenetic tree but not in the EF-1alpha and eRF1 phylogenetic trees. This seems to be mainly because of phylogeny reconstruction artifacts (the long-branch attraction) produced by the acceleration of evolutionary rate of ciliate EF-1alpha and eRF1 sequences. Interaction with the highly divergent actin found in ciliates, or on the contrary, loss of interaction, could explain the acceleration of the evolutionary rate of the EF-1alpha sequences. In the case of ciliate eRF1 sequences, their unusually high evolutionary rate may be related to the deviations in the genetic code usage found in diverse ciliates. These deviations involve a relaxation (or even abolition) of the recognition of one or two stop codons by eRF1. To achieve this, structural changes in eRF1 are needed, and this may affect its evolutionary rate. Eukaryotic translation seems to have followed a mosaic evolution, with its different elements governed by different selective pressures. However, a correlation analysis shows that, beneath the disagreement shown by the different translation proteins, their concerted evolution can still be made apparent when they are compared with other proteins that are not involved in translation.     

Covarion shifts cause a long-branch attraction artifact that unites microsporidia and archaebacteria in EF-1alpha phylogenies

Microsporidia branch at the base of eukaryotic phylogenies inferred from translation elongation factor 1alpha (EF-1alpha) sequences. Because these parasitic eukaryotes are fungi (or close relatives of fungi), it is widely accepted that fast-evolving microsporidian sequences are artifactually "attracted" to the long branch leading to the archaebacterial (outgroup) sequences ("long-branch attraction," or "LBA"). However, no previous studies have explicitly determined the reason(s) why the artifactual allegiance of microsporidia and archaebacteria ("M + A") is recovered by all phylogenetic methods, including maximum likelihood, a method that is supposed to be resistant to classical LBA. Here we show that the M + A affinity can be attributed to those alignment sites associated with large differences in evolutionary site rates between the eukaryotic and archaebacterial subtrees. Therefore, failure to model the significant evolutionary rate distribution differences (covarion shifts) between the ingroup and outgroup sequences is apparently responsible for the artifactual basal position of microsporidia in phylogenetic analyses of EF-1alpha sequences. Currently, no evolutionary model that accounts for discrete changes in the site rate distribution on particular branches is available for either protein or nucleotide level phylogenetic analysis, so the same artifacts may affect many other "deep" phylogenies. Furthermore, given the relative similarity of the site rate patterns of microsporidian and archaebacterial EF-1alpha proteins ("parallel site rate variation"), we suggest that the microsporidian orthologs may have lost some eukaryotic EF-1alpha-specific nontranslational functions, exemplifying the extreme degree of reduction in this parasitic lineage.     

An evaluation of elongation factor 1 alpha as a phylogenetic marker for eukaryotes

Elongation factor 1 alpha (EF-1 alpha) is a highly conserved ubiquitous protein involved in translation that has been suggested to have desirable properties for phylogenetic inference. To examine the utility of EF-1 alpha as a phylogenetic marker for eukaryotes, we studied three properties of EF-1 alpha trees: congruency with other phyogenetic markers, the impact of species sampling, and the degree of substitutional saturation occurring between taxa. Our analyses indicate that the EF-1 alpha tree is congruent with some other molecular phylogenies in identifying both the deepest branches and some recent relationships in the eukaryotic line of descent. However, the topology of the intermediate portion of the EF-1 alpha tree, occupied by most of the protist lineages, differs for different phylogenetic methods, and bootstrap values for branches are low. Most problematic in this region is the failure of all phylogenetic methods to resolve the monophyly of two higher-order protistan taxa, the Ciliophora and the Alveolata. JACKMONO analyses indicated that the impact of species sampling on bootstrap support for most internal nodes of the eukaryotic EF-1 alpha tree is extreme. Furthermore, a comparison of observed versus inferred numbers of substitutions indicates that multiple overlapping substitutions have occurred, especially on the branch separating the Eukaryota from the Archaebacteria, suggesting that the rooting of the eukaryotic tree on the diplomonad lineage should be treated with caution. Overall, these results suggest that the phylogenies obtained from EF-1 alpha are congruent with other molecular phylogenies in recovering the monophyly of groups such as the Metazoa, Fungi, Magnoliophyta, and Euglenozoa. However, the interrelationships between these and other protist lineages are not well resolved. This lack of resolution may result from the combined effects of poor taxonomic sampling, relatively few informative positions, large numbers of overlapping substitutions that obscure phylogenetic signal, and lineage-specific rate increases in the EF-1 alpha data set. It is also consistent with the nearly simultaneous diversification of major eukaryotic lineages implied by the "big-bang" hypothesis of eukaryote evolution.     

Molecular evolution of ciliates (Ciliophora) and some related groups of protozoans

The review summarizes current evidence, including the findings related to molecular phylogeny of ciliates (type Ciliophora) and some related groups of protozoans. Based on comparison of the sequences of genes encoding various ribosomal RNAs (rRNAs), the phylogenetic relationships in seven out of eight known classes of ciliates are discussed. The events related to early branching of the eukaryotic tree are briefly presented. The evolutionary history of amitochondrial protists ids considered with regard to reductionistic evolution and archeozoic hypothesis. The phylogenetic relationships among ciliates and sister groups of apicomplexans and dinoflagellates are considered.     

Insights into the Evolution of Nuclear Dualism in the Ciliates Revealed by Phylogenetic Analysis of rRNA Sequences

The small subunit rRNA gene sequences of the karyorelictean ciliates, Loxodes striatus and Protocruzia sp., and the heterotrichian ciliates, Climacostomum virens and Eufolliculina uhligi , were used to test the evolution of nuclear dualism in the Phylum Ciliophora. Phylogenies derived using a least squares distance method, neighbour joining, and maximum parsimony demonstrate that the karyorelictean ciliates sensu Small and Lynn, 1985 do not form a monophyletic group. However, Loxodes and the heterotrich ciliates form the first branch in the ciliate lineage, and Protocruzia branches, in distance methods, basal to the spirotrich lineage. It is proposed that Protocruzia be removed from the Class Karyorelictea, and placed in closer taxonomic association with the spirotrich lineage. The distribution of nuclear division types along the phylogenetic tree is consistent with the notion that macronuclei incapable of division represent a derived rather than a primitive or "karyorelictid" character trait.     

Long-branch attraction and the rDNA model of early eukaryotic evolution

Phylogenetic analyses of ribosomal RNA genes have become widely accepted as a framework for understanding broad-scale eukaryotic evolution. Nevertheless, conflicts exist between the phylogenetic placement of certain taxa in rDNA trees and their expected position based on fossils, cytology, or protein-encoding gene sequences. For example, pelobiont amoebae appear to be an ancient group based on cytologic features, but they are not among the early eukaryotic brances in rDNA analyses. In this report, the derived position of pelobionts in rDNA trees is shown to be unreliable and likely due to long-branch attraction among more deeply branching sequences. All sequences that branch near the base of the tree suffer from relatively high apparent substitution rates and exhibit greater variation in ssu rDNA sequence length. Moreover, the order of the branches leading from the root of the eukaryotic tree to the base of the so-called "crown taxa" is consistent with a sequential attachment, due to "long-branch" effects, of sequences with increasing rates of evolution. These results suggest that the basal eurkaryotic topology drawn from rDNA analyses may be, in reality, an artifact of variation in the rate of molecular evolution among eukaryotic taxa.     

Is the Felsenstein zone a fly trap?

Although long-branch attraction, the incorrect grouping of long lineages in a phylogeny because of systematic error, has been identified as a potential source of error in phylogenetic analysis for almost two decades, no empirical examples of the phenomenon exist. Here, I outline several criteria for identifying long-branch attraction and apply these criteria to 18S ribosomal DNA (rDNA) sequence data for 13 insects. Parsimony and minimum evolution with p distances group the two longest branches together (those leading to Strepsiptera and Diptera). Simulation studies show that the long branches are long enough to attract. When a tree is assumed in which Strepsiptera and Diptera are separated and many data sets are simulated for that tree (using the parameter estimates for that tree for the original data), parsimony analysis of the simulated data consistently groups Strepsiptera and Diptera. Analyses of the 18S rDNA sequences using methods that are less sensitive to the problem of long-branch attraction estimate trees in which the long branches are separate.     

Escaping from the Felsenstein Zone by Detecting Long Branches in Phylogenetic Data

Long branches in a true phylogeny tend to disrupt hierarchical character covariation (phylogenetic signal) in the distribution of traits among organisms. The distortion of hierarchical structure in character-state matrices can lead to errors in the estimation of phylogenetic relationships and inconsistency of methods of phylogenetic inference. Examination of trees distorted by long-branch attraction will not reveal the identities of problematic taxa, in part because the distortion can mask long branches by reducing inferred branch lengths and through errors in branching order. Here we present a simple method for the detection of taxa whose placement in evolutionary trees is made difficult by the effects of long-branch attraction. The method is an extension of a tree-independent conceptual framework of phylogenetic data exploration (RASA). Taxa that are likely to attract are revealed because long branches leave distinct footprints in the distribution of character states among taxa, and these traces can be directly observed in the error structure of the RASA regression. Problematic taxa are identified using a new diagnostic plot called the taxon variance plot, in which the apparent cladistic and phenetic variances contributed by individual taxa are compared. The procedure for identifying long edges employs algorithms solved in polynomial time and can be applied to morphological, molecular, and mixed characters. The efficacy of the method is demonstrated using simulated evolution and empirical evidence of long branches in a set of recently published sequences. We show that the accuracy of evolutionary trees can be improved by detecting and combating the potentially misleading influences of long-branch taxa.     

Rooting the eukaryotic tree with mitochondrial and bacterial proteins

By exploiting the large body of genome data and the considerable progress in phylogenetic methodology, recent phylogenomic studies have provided new insights into the relationships among major eukaryotic groups. However, confident placement of the eukaryotic root remains a major challenge. This is due to the large evolutionary distance separating eukaryotes from their closest relatives, the Archaea, implying a weak phylogenetic signal and strong long-branch attraction artifacts. Here, we apply a new approach to the rooting of the eukaryotic tree by using a subset of genomic information with more recent evolutionary origin-mitochondrial sequences, whose closest relatives are α-Proteobacteria. For this, we identified and assembled a data set of 42 mitochondrial proteins (mainly encoded by the nuclear genome) and performed Bayesian and maximum likelihood analyses. Taxon sampling includes the recently sequenced Thecamonas trahens, a member of the phylogenetically elusive Apusozoa. This data set confirms the relationships of several eukaryotic supergroups seen before and places the eukaryotic root between the monophyletic "unikonts" and "bikonts." We further show that T. trahens branches sister to Opisthokonta with significant statistical support and question the bikont/excavate affiliation of Malawimonas species. The mitochondrial data set developed here (to be expanded in the future) constitutes a unique alternative means in resolving deep eukaryotic relationships.     

Phylogenetic Relationships of the Nassulida Within the Phylum Ciliophora Inferred from the Complete Small Subunit rRNA Gene Sequences of Furgasonia blochmanni, Obertrumia georgiana, and Pseudomicrothorax dubius

ABSTRACT. Using comparisons of complete small subunit rRNA sequences from the ciliated protozoans Furgasonia blochmanni, Obertrumia georgiana , and Pseudomicrothorax dubius we inferred the phylogenetic position of the Nassulida (Class Nassophorea) within the Ciliophora. In distance matrix analyses the Nassulida share a common ancestry with the colpodean ciliate Colpoda inflata. Distance matrix and parsimony methods convincingly demonstrate that the Nassulida plus Colpodida are members of a complex ciliate assemblage that also includes the oligohymenophorans and phyllopharyngeans. These phylogenetic inferences are largely congruent with recent analyses of 23S-like rRNA gene sequences and morphogenetic features. Groups traditionally thought to represent ancestral lineages now appear as highly derived ciliates. In contrast, heterotrichs which were considered to represent a highly evolved group, diverge at the base of the ciliates.     

LDL binds to surface-expressed human T-cadherin in transfected HEK293 cells and influences homophilic adhesive interactions

The interaction between elongation factor 1alpha (EF-1alpha) and alpha/beta-tubulins has been analyzed in vivo and in vitro. An association of both alpha- and beta-tubulins with EF-1alpha in the lysate of Tetrahymena pyriformis was detected by co-immunoprecipitation analysis. In contrast, in vitro biomolecular interaction analysis with glutathione S-transferase (GST) fusion proteins revealed that GST-beta-tubulin, but not GST-alpha-tubulin, can bind to GST-EF-1alpha. Two beta-tubulin binding sites have been identified to reside in the domains I and III of EF-1alpha. In addition, beta-tubulin itself seems to have two distinct interaction sites for each of the domains. Since domain II of EF-1alpha did not interact with beta-tubulin, we have re-evaluated the phylogenetic status of ciliates using EF-1alpha sequences devoid of domain II. The phylogenetic tree thus obtained was significantly different from that inferred from the whole sequence of EF-1alpha, suggesting the presence of functional constraints on the molecular evolution of EF-1alpha.     

Nuclear and Nucleomorph SSU rDNA Phylogeny in the Cryptophyta and the Evolution of Cryptophyte Diversity

The plastid-bearing members of the Cryptophyta contain two functional eukaryotic genomes of different phylogenetic origin, residing in the nucleus and in the nucleomorph, respectively. These widespread and diverse protists thus offer a unique opportunity to study the coevolution of two different eukaryotic genomes within one group of organisms. In this study, the SSU rRNA genes of both genomes were PCR-amplified with specific primers and phylogenetic analyses were performed on different data sets using different evolutionary models. The results show that the composition of the principal clades obtained from the phylogenetic analyses of both genes was largely congruent, but striking differences in evolutionary rates were observed. These affected the topologies of the nuclear and nucleomorph phylogenies differently, resulting in long-branch attraction artifacts when simple evolutionary models were applied. Deletion of long-branch taxa stabilized the internal branching order in both phylogenies and resulted in a completely resolved topology in the nucleomorph phylogeny. A comparison of the tree topologies derived from SSU rDNA sequences with characters previously used in cryptophyte systematics revealed that the biliprotein type was congruent, but the type of inner periplast component incongruent, with the molecular trees. The latter is indicative of a hidden cellular dimorphism (cells with two periplast types present in a single clonal strain) of presumably widespread occurrence throughout cryptophyte diversity, which, in consequence, has far-reaching implications for cryptophyte systematics as it is practiced today.     

Assessing whether alpha-tubulin sequences are suitable for phylogenetic reconstruction of Ciliophora with insights into its evolution in euplotids

The current understanding of ciliate phylogeny is mainly based on analyses of a single gene, the small subunit ribosomal RNA (SSU-rDNA). However, phylogenetic trees based on single gene sequence are not reliable estimators of species trees, and SSU-rDNA genealogies are not useful for resolution of some branches within Ciliophora. Since congruence between multiple loci is the best tool to determine evolutionary history, we assessed the usefulness of alpha-tubulin gene, a protein-coding gene that is frequently sequenced, for ciliate phylogeny. Here, we generate alpha-tubulin gene sequences of 12 genera and 30 species within the order Euplotida, one of the most frequently encountered ciliate clades with numerous apparently cosmopolitan species, as well as four genera within its putative sister order Discocephalida. Analyses of the resulting data reveal that: 1) the alpha-tubulin gene is suitable phylogenetic marker for euplotids at the family level, since both nucleotide and amino acid phylogenies recover all monophyletic euplotid families as defined by both morphological criteria and SSU-rDNA trees; however, alpha-tubulin gene is not a good marker for defining species, order and subclass; 2) for seven out of nine euplotid species for which paralogs are detected, gene duplication appears recent as paralogs are monophyletic; 3) the order Euplotida is non-monophyletic, and the family Uronychiidae with sequences from four genera, is non-monophyletic; and 4) there is more genetic diversity within the family Euplotidae than is evident from dargyrome (geometrical pattern of dorsal "silverline system" in ciliates) patterns, habit and SSU-rDNA phylogeny, which indicates the urgent need for taxonomic revision in this area.     

Comparative analysis of the mitochondrial cytochrome c oxidase subunit I (COI) gene in ciliates (Alveolata,Ciliophora) and evaluation of its suitability as a biodiversity marker

The mitochondrial cytochrome c oxidase subunit 1 (COI) gene of ciliates was first successfully sequenced in species of the genera Tetrahymena and Paramecium (Class Oligohymenophorea). The sequence of the COI gene is extremely divergent from other eukaryotes and includes an insert, which is over 300 nucleotides long. In this study, we designed a primer pair that successfully amplified the COI gene of ciliates from five different classes: Heterotrichea, Spirotrichea, Oligohymenophorea, Nassophorea and Colpodea. These classes represent the diversity of the phylum Ciliophora very well, since they are widely distributed on the ciliate small subunit rRNA tree. The amplified region is approximately 850 nucleotides long and corresponds to the general barcoding region; it also includes the insert region. In this study, 58 new COI sequences from over 38 species in 13 orders are analysed and compared, and distance trees are constructed. While the COI gene shows high divergence within ciliates, the insert region, which is present in all classes, is even more divergent. Genetic distances calculated with and without the insert region remain in the same range at the intraspecific level, but they differ considerably at or above genus level. This suggests that the entire barcoding region is under similar selective constraints and that the evolutionary rate of the ciliate COI is extremely high and shows unequal rate variation. Although many problems still remain regarding standardization of barcoding methods in ciliates, the development of a universal or almost universal primer combination for the Phylum Ciliophora represents important progress. As shown in four examples, the resolution of COI at the intraspecific level is much greater than that of any nuclear genes and shows great potential to (1) identify species based on molecular data if a reliable database exists, and (2) resolve the relationships of closely related ciliate taxa and uncover cryptic species.     

Phylogenetic position of Licnophora,Lechriopyla, and Schizocaryum,three unusual ciliates (phylum Ciliophora) endosymbiotic in echinoderms (phylum Echinodermata)

Various echinoderms are colonized by species from several classes of the Phylum Ciliophora, indicating that the echinoderm "habitat" has been invaded independently on numerous occasions throughout evolutionary history. Two "echinoderm" ciliates whose phylogenetic positions have been problematic are Licnophora macfarlandi Stevens, 1901 and Schizocaryum dogieli Poljansky and Golikova, 1957. Licnophora macfarlandi is an endosymbiont of the respiratory trees of holothuroids, and S. dogieli is found in the esophagus of echinoids. A third species, Lechriopyla mystax Lynch, 1930, is a plagiopylid ciliate found in the intestine of echinoids. Host echinoderms were collected near the Friday Harbor Laboratories, San Juan Island, WA. Specimens of S. dogieli and L. mystax were obtained from the esophagus and intestine, respectively, of the sea urchin Strongylocentrotus pallidus. Specimens of L. macfarlandi were collected from the fluid obtained from the respiratory trees of Parastichopus californicus. Using small subunit ribosomal RNA (SSrRNA) sequences of these three ciliates and a global alignment of SSrRNA sequences of other ciliates, we established the following. 1) Licnophora is a spirotrich ciliate, clearly related to the hypotrichs and stichotrichs; this is corroborated by its possession of macronuclear replication bands. 2) Lechriopyla is the sister genus to Plagiopyla and is a member of the Class Plagiopylea, which was predicted based on its cytology. 3) Schizocaryum clusters in the Class Oligohymenophorea and is most closely related to the scuticociliates; there are currently no morphological features known to relate Schizocaryum to the scuticociliates.     

Evolution of Histone H4 and H3 Genes in Different Ciliate Lineages

The histones H4 are known as highly conserved proteins. However, in ciliates a high degree of variation was found compared both to other eukaryotes and between the ciliate species. To date, only H4 histones of species belonging to two distantly related classes have been investigated. In order to obtain more detailed information on histone H4 variation in ciliates we undertook a comprehensive sequence analysis of PCR-amplified internal H4 fragments from 12 species belonging to seven out of the nine currently recognized ciliate classes. In addition, we used PCR primers to amplify longer fragments of H3 and H4 genes including the intergenic region. The encoded amino acid sequences reveal a high number of differences when compared with those of other eukaryotes and the ciliate species investigated. Furthermore, in some species H4 gene variants were detected, which result in amino acid differences. The greatest number of substitutions and insertions found was in the amino terminal region of the H4 histones. However, all sequences possess a conserved region corresponding to those of all other eukaryotic H4 histones. The histone gene variations were used to reconstruct phylogenetic relationships. The tree from our data matches perfectly with the ribosomal RNA data: The heterotrichs, which were considered as a late branching lineage, diverge at the base of the ciliate tree and groups formerly thought to represent ancestral lineages now appear as highly derived ciliates. Received: 4 April 1997 / Accepted: 1 August 1997     

Can sensitivity analysis help to detect long-branch attraction?

The behavior of nodal support and stability in the presence of long branches were examined under simulations and an analysis of real data. Relatively short branches were typically correctly resolved, received high bootstrap support, and were stable in sensitivity analyses. Longer branches received lower support and stability measures, and were often incorrectly resolved due to the long-branch attraction. Support and stability does not always correlate, and in the case of mammalian mitochondrial tree, well supported but unstable nodes were typically associated with long-branch attraction. Very long branches, on the other hand, may be incorrectly resolved with high support and stability indices. These patterns were observed both in simulations, and in the real data. The results indicate that sensitivity analysis may help to reveal phylogenetic uncertainty hidden behind artificially high support.     

The Root of the Tree of Life in the Light of the Covarion Model

A few duplicated genes have been found useful to root the universal tree of life. Despite controversial results, the consensus led to locate the root in the eubacterial branch. However, we demonstrated (Philippe and Forterre 1999) that all these markers were in fact unsuitable for any firm conclusion, mainly because of their high level of mutational saturation, which masks a major part of the phylogenetic signal. But then, the very persistence of signal for events as early as the separation of the three domains becomes puzzling. This paradox was studied here for translation elongation factor proteins, EF-1α and EF-2, which appeared to be one of the least confusing markers. We showed that these proteins do not conform to a classical rate-across-sites pattern, as those modeled by a gamma law, but rather to a covarion-based model, because the evolutionary rate of a given position often changes between taxonomic groups. Conservation of the very ancient signal can thus be better explained by the covarion model: a substitution can occur in deep branches, and the position remains constant afterward, as ``fossilized' by a change of covation. As no reconstruction method has up to now taken into account this complex model, we devised a simple method for extracting the phylogenetic signal, by considering the variability of sequence positions within predefined phylogenetic groups. We showed that noise quantitatively prevailed upon signal. Parsimony will produce erroneous topologies, because it has to minimize primarily the number of steps of the noise. In contrast, our method effectively concentrated the signal and was more suitable for inferring ancient events. We consequently found the eubacterial rooting to be presumably due to a long branch attraction artifact, because of the higher evolutionary rate of Eubacteria for these proteins. Among the two other rooting possibilities, the eukaryotic rooting appeared to be more supported, although not enough to be conclusive.     

Resting Cysts in the Ciliate Class Polyhymenophorea: Phylogenetic Implications1

Distinctive organic-walled resting cysts of at least three different types with a highly conservative morphology appear to characterize specific orders or groups of genera within the Class Polyhymenophorea (Protozoa, Ciliophora), contrasting markedly with the great diversity of form seen in trophic stages. Polyhymenophorean ciliates have been considered in the past to form a cohesive class within the Phylum Ciliophora and, possibly, to represent the pinnacle of ciliate evolution. Evidence from cysts challenges the cohesive nature of the class, suggesting that the hypotrichs should be subdivided and that they have a different phylogenetic origin from the heterotrichs, tintinnids, and oligotrichs.     

Evolutionary history of "early-diverging" eukaryotes: the excavate taxon Carpediemonas is a close relative of Giardia

Diplomonads, such as Giardia, and their close relatives retortamonads have been proposed as early-branching eukaryotes that diverged before the acquisition-retention of mitochondria, and they have become key organisms in attempts to understand the evolution of eukaryotic cells. In this phylogenetic study we focus on a series of eukaryotes suggested to be relatives of diplomonads on morphological grounds, the "excavate taxa". Phylogenies of small subunit ribosomal RNA (SSU rRNA) genes, alpha-tubulin, beta-tubulin, and combined alpha- + beta-tubulin all scatter the various excavate taxa across the diversity of eukaryotes. But all phylogenies place the excavate taxon Carpediemonas as the closest relative of diplomonads (and, where data are available, retortamonads). This novel relationship is recovered across phylogenetic methods and across various taxon-deletion experiments. Statistical support is strongest under maximum-likelihood (ML) (when among-site rate variation is modeled) and when the most divergent diplomonad sequences are excluded, suggesting a true relationship rather than an artifact of long-branch attraction. When all diplomonads are excluded, our ML SSU rRNA tree actually places retortamonads and Carpediemonas away from the base of the eukaryotes. The branches separating excavate taxa are mostly not well supported (especially in analyses of SSU rRNA data). Statistical tests of the SSU rRNA data, including an "expected likelihood weights" approach, do not reject trees where excavate taxa are constrained to be a clade (with or without parabasalids and Euglenozoa). Although diplomonads and retortamonads lack any mitochondria-like organelle, Carpediemonas contains double membrane-bounded structures physically resembling hydrogenosomes. The phylogenetic position of Carpediemonas suggests that it will be valuable in interpreting the evolutionary significance of many molecular and cellular peculiarities of diplomonads.