Unusually high evolutionary rate of the elongation factor 1 alpha genes from the Ciliophora and its impact on the phylogeny of eukaryotes

The elongation factor 1 alpha (EF-1 alpha) has become widely employed as a phylogenetic marker for studying eukaryotic evolution. However, a disturbing problem, the artifactual polyphyly of ciliates, is always observed. It has been suggested that the addition of new sequences will help to circumvent this problem. Thus, we have determined 15 new ciliate EF-1 alpha sequences, providing for a more comprehensive taxonomic sampling of this phylum. These sequences have been analyzed together with a representation of eukaryotic sequences using distance-, parsimony-, and likelihood-based phylogenetic methods. Such analyses again failed to recover the monophyly of Ciliophora. A study of the substitution rate showed that ciliate EF-1 alpha genes exhibit a high evolutionary rate, produced in part by an increased number of variable positions. This acceleration could be related to alterations of the accessory functions acquired by this protein, likely to those involving interactions with the cytoskeleton, which is very modified in the Ciliophora. The high evolutionary rate of these sequences leads to an artificial basal emergence of some ciliates in the eukaryotic tree by effecting a long-branch attraction artifact that produces an asymmetric topology for the basal region of the tree. The use of a maximum-likelihood phylogenetic method (which is less sensitive to long-branch attraction) and the addition of sequences to break long branches allow retrieval of more symmetric topologies, which suggests that the asymmetric part of the tree is most likely artifactual. Therefore, the sole reliable part of the tree appears to correspond to the apical symmetric region. These kinds of observations suggest that the general eukaryotic evolution might have consisted of a massive radiation followed by an increase in the evolutionary rates of certain groups that emerge artificially as early branches in the asymmetric base of the tree. Ciliates in the case of the EF-1 alpha genes would offer clear evidence for this hypothesis.     

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.     

Congruent evidence from alpha-tubulin and beta-tubulin gene phylogenies for a zygomycete origin of microsporidia

The origin of microsporidia and the evolutionary relationships among the major lineages of fungi have been examined by molecular phylogeny using alpha-tubulin and beta-tubulin. Chytrids, basidiomycetes, ascomycetes, and microsporidia were all recovered with high support, and the zygomycetes were consistently paraphyletic. The microsporidia were found to branch within zygomycetes, and showed relationships with members of the Entomophthorales and Zoopagales. This provides support for the microsporidia having evolved from within the fungi, however, the tubulin genes are difficult to interpret unambiguously since fungal and microsporidian tubulins are very divergent. Rapid evolutionary rates a characteristic of practically all microsporidian genes studied, so determining their evolutionary history will never be completely free of such difficulties. While the tubulin phylogenies do not provide a decisive conclusion, they do further narrow the probable origin of microsporidia to a zygomycete-like ancestor.     

Evidence from beta-tubulin phylogeny that microsporidia evolved from within the fungi

Microsporidia are obligate intracellular parasites that were thought to be an ancient eukaryotic lineage based on molecular phylogenies using ribosomal RNA and translation elongation factors. However, this ancient origin of microsporidia has been contested recently, as several other molecular phylogenies suggest that microsporidia are closely related to fungi. Most of the protein trees that place microsporidia with fungi are not well sampled, however, and it is impossible to resolve whether microsporidia evolved from a fungus or from a protistan relative of fungi. We have sequenced beta-tubulins from 3 microsporidia, 4 chytrid fungi, and 12 zygomycete fungi, expanding the representation of beta-tubulin to include all four fungal divisions and a wide diversity of microsporidia. In phylogenetic trees including these new sequences, the overall topology of the fungal beta-tubulins generally matched the expected relationships among the four fungal divisions, although the zygomycetes were polyphyletic in some analyses. The microsporidia consistently fell within this fungal diversification, and not as a sister group to fungi. Overall, beta-tubulin phylogeny suggests that microsporidia evolved from a fungus sometime after the divergence of chytrids. We also found that chytrid alpha- and beta-tubulins are much less divergent than are tubulins from other fungi or microsporidia. In trees in which the only fungal representatives were the chytrids, microsporidia still branched with fungi (i.e., with chytrids), suggesting that the affiliation between microsporidian and fungal tubulins is not an artifact of long-branch attraction.     

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.     

Interaction of process partitions in phylogenetic analysis: an example from the swallowtail butterfly genus Papilio

In this study, we explored how the concept of the process partition may be applied to phylogenetic analysis. Sequence data were gathered from 23 species and subspecies of the swallowtail butterfly genus Papilio, as well as from two outgroup species from the genera Eurytides and Pachliopta. Sequence data consisted of 1,010 bp of the nuclear protein-coding gene elongation factor-1 alpha (EF-1 alpha) as well as the entire sequences (a total of 2,211 bp) of the mitochondrial protein-coding genes cytochrome oxidase I and cytochrome oxidase II (COI and COII). In order to examine the interaction between the nuclear and mitochondrial partitions in a combined analysis, we used a method of visualizing branch support as a function of partition weight ratios. We demonstrated how this method may be used to diagnose error at different levels of a tree in a combined maximum-parsimony analysis. Further, we assessed patterns of evolution within and between subsets of the data by implementing a multipartition maximum-likelihood model to estimate evolutionary parameters for various putative process partitions. COI third positions have an estimated average substitution rate more than 15 times that of EF-1 alpha, while COII third positions have an estimated average substitution rate more than 22 times that of EF-1 alpha. Ultimately, we found that although the mitochondrial and nuclear data were not significantly incongruent, homoplasy in the fast-evolving mitochondrial data confounded the resolution of basal relationships in the combined unweighted parsimony analysis despite the fact that there was relatively strong support for the relationships in the nuclear data. We conclude that there may be shortcomings to the methods of "total evidence" and "conditional combination" because they may fail to detect or accommodate the type of confounding bias we found in our data.     

Phylogenetic Analysis of the Complete Genome Sequence of <Emphasis Type="Italic">Encephalitozoon cuniculi</Emphasis> Supports the Fungal Origin of Microsporidia and Reveals a High Frequency of Fast-Evolving Genes

Microsporidia are unicellular eukaryotes living as obligate intracellular parasites. Lacking mitochondria, they were initially considered as having diverged before the endosymbiosis at the origin of mitochondria. That microsporidia were primitively amitochondriate was first questioned by the discovery of microsporidial sequences homologous to genes encoding mitochondrial proteins and then refuted by the identification of remnants of mitochondria in their cytoplasm. Various molecular phylogenies also cast doubt on the early divergence of microsporidia, these organisms forming a monophyletic group with or within the fungi. The 2001 proteins putatively encoded by the complete genome of Encephalitozoon cuniculi provided powerful data to test this hypothesis. Phylogenetic analysis of 99 proteins selected as adequate phylogenetic markers indicated that the E. cuniculi sequences having the lowest evolutionary rates preferentially clustered with fungal sequences or, more rarely, with both animal and fungal sequences. Because sequences with low evolutionary rates are less sensitive to the long-branch attraction artifact, we concluded that microsporidia are evolutionarily related to fungi. This analysis also allowed comparing the accuracy of several phylogenetic algorithms for a fast-evolving lineage with real rather than simulated sequences.This article contains online supplementary material.Reviewing Editor: Dr. Wen-Hsiung LiSupplementary material is available at     

Phylogenetic Relationships of Fungi, Plantae, and Animalia Inferred from Homologous Comparison of Ribosomal Proteins

The complete set of available ribosomal proteins was utilized, at both the peptidic and the nucleotidic level, to establish that plants and metazoans form two sister clades relative to fungi. Different phylogenetic inference methods are applied to the sequence data, using archeans as the outgroup. The evolutionary length of the internal branch within the eukaryotic crown trichotomy is demonstrated to be, at most, one-tenth of the evolutionary length of the branch leading to the cenancester of these three kingdoms. Received: 1 November 1997 / Accepted: 7 January 1998     

α‐ and β‐Tubulin Phylogenies Support a Close Relationship Between the Microsporidia Brachiola algerae and Antonospora locustae

ABSTRACT. Microsporidia are a large and diverse group of intracellular parasites related to fungi. Much of our understanding of the relationships between microsporidia comes from phylogenies based on a single gene, the small subunit (SSU) rRNA, because only this gene has been sampled from diverse microsporidia. However, SSUrRNA trees are limited in their ability to resolve basal branches and some microsporidian affiliations are inconsistent between different analyses. Protein phylogenies have provided insight into relationships within specific groups of microsporidia, but have rarely been applied to the group as a whole. We have sequenced α‐ and β‐tubulins from microsporidia from three different subgroups, including representatives from what have previously been inferred to be the basal branches, allowing the broadest sampled protein‐based phylogenetic analysis to date. Although some relationships remain unresolved, many nodes uniting subgroups are strongly supported and consistent in both individual trees as well as a concatenate of both tubulins. One such relationship that was previously unclear is between Brachiola algerae and Antonospora locustae, and their close association with Encephalitozoon and Nosema. Also, an uncultivated microsporidian that infects cyclopoid copepods is shown to be related to Edhazardia aedis.     

An empirical assessment of long-branch attraction artefacts in deep eukaryotic phylogenomics

In the context of exponential growing molecular databases, it becomes increasingly easy to assemble large multigene data sets for phylogenomic studies. The expected increase of resolution due to the reduction of the sampling (stochastic) error is becoming a reality. However, the impact of systematic biases will also become more apparent or even dominant. We have chosen to study the case of the long-branch attraction artefact (LBA) using real instead of simulated sequences. Two fast-evolving eukaryotic lineages, whose evolutionary positions are well established, microsporidia and the nucleomorph of cryptophytes, were chosen as model species. A large data set was assembled (44 species, 133 genes, and 24,294 amino acid positions) and the resulting rooted eukaryotic phylogeny (using a distant archaeal outgroup) is positively misled by an LBA artefact despite the use of a maximum likelihood-based tree reconstruction method with a complex model of sequence evolution. When the fastest evolving proteins from the fast lineages are progressively removed (up to 90%), the bootstrap support for the apparently artefactual basal placement decreases to virtually 0%, and conversely only the expected placement, among all the possible locations of the fast-evolving species, receives increasing support that eventually converges to 100%. The percentage of removal of the fastest evolving proteins constitutes a reliable estimate of the sensitivity of phylogenetic inference to LBA. This protocol confirms that both a rich species sampling (especially the presence of a species that is closely related to the fast-evolving lineage) and a probabilistic method with a complex model are important to overcome the LBA artefact. Finally, we observed that phylogenetic inference methods perform strikingly better with simulated as opposed to real data, and suggest that testing the reliability of phylogenetic inference methods with simulated data leads to overconfidence in their performance. Although phylogenomic studies can be affected by systematic biases, the possibility of discarding a large amount of data containing most of the nonphylogenetic signal allows recovering a phylogeny that is less affected by systematic biases, while maintaining a high statistical support.     

Phylogenetic positions of several amitochondriate protozoa-Evidence from phylogenetic analysis of DNA topoisomerase II

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, Entamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group―Archezoa. The main evidence for this is their 'lacking mitochondria' and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of 'long-branch attraction' appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplomonads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be polyphyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.     

Phylogenetic positions of several amitochondriate protozoa

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, Entamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group―Archezoa. The main evidence for this is their ‘lacking mitochondria’ and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of ‘long-branch attraction’ appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplo-monads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be poly-phyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.     

Microsporidia: accumulating molecular evidence that a group of amitochondriate and suspectedly primitive eukaryotes are just curious fungi

Microsporidia are obligate intracellular parasites that have long been considered to be primitive eukaryotes, both on the basis of morphological features and on the basis of molecular, mainly ribosomal RNA-based, phylogenies. However, accumulating sequence data and the use of more sophisticated tree construction methods now seem to suggest that microsporidia share a common origin with fungi and are therefore most probably just curious fungi. In this paper, we describe the current views on the phylogenetic position of the microsporidia and present additional evidence for a close relationship between fungi and microsporidia on the basis of reanalyzed ribosomal RNA data. In this respect, the importance of incorporating detailed knowledge of the substitution pattern of sequences into phylogenetic methods is discussed.     

Phylogenetic analysis of 70 kD heat shock protein sequences suggests a chimeric origin for the eukaryotic cell nucleus

Background: The evolutionary relationships between archaebacteria, eubacteria and eukaryotic cells are of central importance in biology. The current view is that each of these three groups of organisms constitutes a monophyletic domain, and that eukaryotic cells have evolved from an archaebacterial ancestor. Recent studies on a number of highly conserved protein sequences do not, however, support this view and raise important questions concerning the evolutionary relationships between all extant organisms, particularly regarding the origin of eukaryotic cells.Results We have used sequences of 70 kD heat shock protein (hsp70) — the most conserved protein found to date in all species — to examine the evolutionary relationship between various species. We have obtained two new archaebacterial hsp70 sequences from the species, Thermoplasma acidophilum and Halobacterium cutirubrum. A global comparison of hsp70 sequences, including our two new sequences, shows that all known archaebacterial homologs share a number of sequence signatures with the Gram-positive group of bacteria that are not found in any other prokaryotic or eukaryotic species. In contrast, the eukaryotic homologs are shown to share a number of unique sequence features with the Gram-negative bacteria that are not present in any archaebacteria. Detailed phylogenetic analyses of hsp70 sequences strongly support a specific evolutionary relationship between archaebacteria and Gram-positive bacteria on the one hand, and Gram-negative bacteria and eukaryotes on the other. The phylogenetic analyses also indicate a polyphyletic branching of archaebacteria within the Gram-positive species. The possibility that the observed relationships are due to horizontal gene transfers can be excluded on the basis of sequence characteristics of different groups of homologs.Conclusion Our results do not support the view that archaebacteria constitute a monophyletic domain, but instead suggest a close evolutionary linkage between archaebacteria and Gram-positive bacteria. Furthermore, in contrast to the presently accepted view, eukaryotic hsp70s show a close and specific relationship to those from Gram-negative species. To explain the phylogenies based on different gene sequences, a chimeric model for the origin of the eukaryotic cell nucleus involving fusion between an archaebacterium and a Gram-negative eubacterium is proposed. Several predictions from the chimeric model are discussed.     

Compositional statistics: An improvement of evolutionary parsimony and its application to deep branches in the tree of life

Summary We present compositional statistics, a new method of phylogenetic inference, which is an extension of evolutionary parsimony. Compositional statistics takes account of the base composition of the compared sequences by using nucleotide positions that evolutionary parsimony ignores. It shares with evolutionary parsimony the features of rate invariance and the fundamental distinction between transitions and transversions. Of the presently available methods of phylogenetic inference, compositional statistics is based on the fewest and mildest assumptions about the mode of DNA sequence evolution. It is therefore applicable to phylogenetic studies of the most distantly related organisms or molecules. This was illustrated by analyzing conservative positions in the DNA sequences of the large subunit of RNA polymerase from three archaebacterial groups, a eubacterium, a chloroplast, and the three eukaryotic polymerases. Internally consistent results, which are in accord with our knowledge of organelle origin and archaebacterial physiology, were achieved.     

Multiple losses of sex within a single genus of Microsporidia

Background  

Most asexual eukaryotic lineages have arisen recently from sexual ancestors and contain few ecologically distinct species, providing evidence for long-term advantages of sex. Ancient asexual lineages provide rare exceptions to this rule and so can yield valuable information relating to the evolutionary forces underlying the maintenance of sex. Microsporidia are parasitic, unicellular fungi. They include many asexual species which have traditionally been grouped together into large, presumably ancient taxonomic groups. However, these putative ancient asexual lineages have been identified on the basis of morphology, life cycles and small subunit ribosomal RNA (16S rRNA) gene sequences, all of which hold questionable value in accurately inferring phylogenetic relationships among microsporidia.     

Are Ichthyosporea animals or fungi? Bayesian phylogenetic analysis of elongation factor 1alpha of Ichthyophonus irregularis

Ichthyosporea is a recently recognized group of morphologically simple eukaryotes, many of which cause disease in aquatic organisms. Ribosomal RNA sequence analyses place Ichthyosporea near the divergence of the animal and fungal lineages, but do not allow resolution of its exact phylogenetic position. Some of the best evidence for a specific grouping of animals and fungi (Opisthokonta) has come from elongation factor 1alpha, not only phylogenetic analysis of sequences but also the presence or absence of short insertions and deletions. We sequenced the EF-1alpha gene from the ichthyosporean parasite Ichthyophonus irregularis and determined its phylogenetic position using neighbor-joining, parsimony and Bayesian methods. We also sequenced EF-1alpha genes from four chytrids to provide broader representation within fungi. Sequence analyses and the presence of a characteristic 12 amino acid insertion strongly indicate that I. irregularis is a member of Opisthokonta, but do not resolve whether I. irregularis is a specific relative of animals or of fungi. However, the EF-1alpha of I. irregularis exhibits a two amino acid deletion heretofore reported only among fungi.     

Molecular phylogenies based on ribosomal protein L11, L1, L10, and L12 sequences

Summary Available sequences that correspond to the E. coli ribosomal proteins L11, L1, L10, and L12 from eubacteria, archaebacteria, and eukaryotes have been aligned. The alignments were analyzed qualitatively for shared structural features and for conservation of deletions or insertions. The alignments were further subjected to quantitative phylogenetic analysis, and the amino acid identity between selected pairs of sequences was calculated. In general, eubacteria, archaebacteria, and eukaryotes each form coherent and well-resolved nonoverlapping phylogenetic domains. The degree of diversity of the four proteins between the three groups is not uniform. For L11, the eubacterial and archaebacterial proteins are very similar whereas the eukaryotic L11 is clearly less similar. In contrast, in the case of the L12 proteins and to a lesser extent the L10 proteins, the archaebacterial and eukaryotic proteins are similar whereas the eubacterial proteins are different. The eukaryotic L1 equivalent protein has yet to be identified. If the root of the universal tree is near or within the eubacterial domain, our ribosomal protein-based phylogenies indicate that archaebacteria are monophyletic. The eukaryotic lineage appears to originate either near or within the archaebacterial domain. Correspondence to: P. Dennis     

The phylogenetic position of red algae revealed by multiple nuclear genes from mitochondria-containing eukaryotes and an alternative hypothesis on the origin of plastids

Abstract Red algae are one of the main photosynthetic eukaryotic lineages and are characterized by primitive features, such as a lack of flagella and the presence of phycobiliproteins in the chloroplast. Recent molecular phylogenetic studies using nuclear gene sequences suggest two conflicting hypotheses (monophyly versus non-monophyly) regarding the relationships between red algae and green plants. Although kingdom-level phylogenetic analyses using multiple nuclear genes from a wide-range of eukaryotic lineages were very recently carried out, they used highly divergent gene sequences of the cryptomonad nucleomorph (as the red algal taxon) or incomplete red algal gene sequences. In addition, previous eukaryotic phylogenies based on nuclear genes generally included very distant archaebacterial sequences (designated as the outgroup) and/or amitochondrial organisms, which may carry unusual gene substitutions due to parasitism or the absence of mitochondria. Here, we carried out phylogenetic analyses of various lineages of mitochondria-containing eukaryotic organisms using nuclear multigene sequences, including the complete sequences from the primitive red alga Cyanidioschyzon merolae. Amino acid sequence data for two concatenated paralogous genes (α- and β-tubulin) from mitochondria-containing organisms robustly resolved the basal position of the cellular slime molds, which were designated as the outgroup in our phylogenetic analyses. Phylogenetic analyses of 53 operational taxonomic units (OTUs) based on a 1525-amino-acid sequence of four concatenated nuclear genes (actin, elongation factor-1α, α-tubulin, and β-tubulin) reliably resolved the phylogeny only in the maximum parsimonious (MP) analysis, which indicated the presence of two large robust monophyletic groups (Groups A and B) and the basal eukaryotic lineages (red algae, true slime molds, and amoebae). Group A corresponded to the Opisthokonta (Metazoa and Fungi), whereas Group B included various primary and secondary plastid-containing lineages (green plants, glaucophytes, euglenoids, heterokonts, and apicomplexans), Ciliophora, Kinetoplastida, and Heterolobosea. The red algae represented the sister lineage to Group B. Using 34 OTUs for which essentially the entire amino acid sequences of the four genes are known, MP, distance, quartet puzzling, and two types of maximum likelihood (ML) calculations all robustly resolved the monophyly of Group B, as well as the basal position of red algae within eukaryotic organisms. In addition, phylogenetic analyses of a concatenated 4639-amino-acid sequence for 12 nuclear genes (excluding the EF-2 gene) of 12 mitochondria-containing OTUs (including C. merolae) resolved a robust non-sister relationship between green plants and red algae within a robust monophyletic group composed of red algae and the eukaryotic organisms belonging to Group B. A new scenario for the origin and evolution of plastids is suggested, based on the basal phylogenetic position of the red algae within the large clade (Group B plus red algae). The primary plastid endosymbiosis likely occurred once in the common ancestor of this large clade, and the primary plastids were subsequently lost in the ancestor(s) of the Discicristata (euglenoids, Kinetoplastida, and Heterolobosea), Heterokontophyta, and Alveolata (apicomplexans and Ciliophora). In addition, a new concept of “Plantae” is proposed for phototrophic and nonphototrophic organisms belonging to Group B and red algae, on the basis of the common history of the primary plastid endosymbiosis. The Plantae include primary plastid-containing phototrophs and nonphototrophic eukaryotes that possibly contain genes of cyanobacterial origin acquired in the primary endosymbiosis.