Inference of the Phylogenetic Position of Oxymonads Based on Nine Genes: Support for Metamonada and Excavata

Phylogenetic and morphological diversity of free-living diplomonads

Diplomonadida is a lineage of anaerobic protists belonging to Fornicata, Metamonada. Most diplomonads are endobiotic or parasitic, such as Giardia intestinalis, which is a famous human pathogen, but several free-living species exist as well. Although it has been proposed that the free-living diplomonads are descendants of endobiotic organisms and thus interesting from the evolutionary point of view, they have been largely neglected. We obtained 58 cultures of free-living diplomonads belonging to four genera (Hexamita, Trepomonas, Gyromonas, and Trimitus) and six strains of endobiotic diplomonads and analyzed their SSU rRNA gene sequences. We also studied light-microscopic morphology of selected strains and the ultrastructure of Trepomonas rotans for the first time. Our phylogenetic analysis showed that the genus Hexamita, and, possibly, also the genus Trepomonas, are polyphyletic. Trepomonas rotans, which may represent a novel genus, is unique among Diplomonadida by having the cell covered in scales. Our results suggest that the evolution of the endobiotic life style and cell organization in diplomonads is more complicated than previously thought.

A mitochondrion-free eukaryote contains proteins capable of import into an exogenous mitochondrion-related organelle

The endobiotic flagellate Monocercomonoides exilis is the only known eukaryote to have lost mitochondria and all its associated proteins in its evolutionary past. This final stage of the mitochondrial evolutionary pathway may serve as a model to explain events at their very beginning such as the initiation of protein import. We have assessed the capability of proteins from this eukaryote to enter emerging mitochondria using a specifically designed in vitro assay. Hydrogenosomes (reduced mitochondria) of Trichomonas vaginalis were incubated with a soluble protein pool derived from a cytosolic fraction of M. exilis, and proteins entering hydrogenosomes were subsequently detected by mass spectrometry. The assay detected 19 specifically and reproducibly imported proteins, and in 14 cases the import was confirmed by the overexpression of their tagged version in T. vaginalis. In most cases, only a small portion of the signal reached the hydrogenosomes, suggesting specific but inefficient transport. Most of these proteins represent enzymes of carbon metabolism, and none exhibited clear signatures of proteins targeted to hydrogenosomes or mitochondria, which is consistent with their inefficient import. The observed phenomenon may resemble a primaeval type of protein import which might play a role in the establishment of the organelle and shaping of its proteome in the initial stages of endosymbiosis.

High quality genome assembly of the amitochondriate eukaryote Monocercomonoides exilis

Monocercomonoides exilis is considered the first known eukaryote to completely lack mitochondria. This conclusion is based primarily on a genomic and transcriptomic study which failed to identify any mitochondrial hallmark proteins. However, the available genome assembly has limited contiguity and around 1.5 % of the genome sequence is represented by unknown bases. To improve the contiguity, we re-sequenced the genome and transcriptome of M. exilis using Oxford Nanopore Technology (ONT). The resulting draft genome is assembled in 101 contigs with an N50 value of 1.38 Mbp, almost 20 times higher than the previously published assembly. Using a newly generated ONT transcriptome, we further improve the gene prediction and add high quality untranslated region (UTR) annotations, in which we identify two putative polyadenylation signals present in the 3′UTR regions and characterise the Kozak sequence in the 5′UTR regions. All these improvements are reflected by higher BUSCO genome completeness values. Regardless of an overall more complete genome assembly without missing bases and a better gene prediction, we still failed to identify any mitochondrial hallmark genes, thus further supporting the hypothesis on the absence of mitochondrion.

Mechanistic concepts of iron-sulfur protein biogenesis in Biology

Iron-sulfur (Fe/S) proteins are present in virtually all living organisms and are involved in numerous cellular processes such as respiration, photosynthesis, metabolic reactions, nitrogen fixation, radical biochemistry, protein synthesis, antiviral defense, and genome maintenance. Their versatile functions may go back to the proposed role of their Fe/S cofactors in the origin of life as efficient catalysts and electron carriers. More than two decades ago, it was discovered that the in vivo synthesis of cellular Fe/S clusters and their integration into polypeptide chains requires assistance by complex proteinaceous machineries, despite the fact that Fe/S proteins can be assembled chemically in vitro. In prokaryotes, three Fe/S protein biogenesis systems are known; ISC, SUF, and the more specialized NIF. The former two systems have been transferred by endosymbiosis from bacteria to mitochondria and plastids, respectively, of eukaryotes. In their cytosol, eukaryotes use the CIA machinery for the biogenesis of cytosolic and nuclear Fe/S proteins. Despite the structural diversity of the protein constituents of these four machineries, general mechanistic concepts underlie the complex process of Fe/S protein biogenesis. This review provides a comprehensive and comparative overview of the various known biogenesis systems in Biology, and summarizes their common or diverging molecular mechanisms, thereby illustrating both the conservation and diverse adaptions of these four machineries during evolution and under different lifestyles. Knowledge of these fundamental biochemical pathways is not only of basic scientific interest, but is important for the understanding of human 'Fe/S diseases' and can be used in biotechnology.

Phylogenetic Reconstruction Shows Independent Evolutionary Origins of Mitochondrial Transcription Factors from an Ancient Family of RNA Methyltransferase Proteins

Here, we generate a robust phylogenetic framework for the rRNA adenine N(6)-methyltransferase (RAMTase) protein family that shows a more ancient and complex evolutionary history within the family than previously reported. RAMTases occur universally by descent across the three domains of life, and typical orthologs within the family perform methylation of the small subunits of ribosomal RNA (rRNA). However, within the RAMTase family, two different groups of mitochondrial transcription factors, mtTFB1 and mtTFB2, have evolved in eukaryotes through neofunctionalization. Previous phylogenetic analyses have suggested that mtTFB1 and mtTFB2 comprise sister clades that arose via gene duplication, which occurred sometime following the endosymbiosis event that produced the mitochondrion. Through dense and taxonomically broad sampling of RAMTase family members especially within bacteria, we found that these eukaryotic mitochondrial transcription factors, mtTFB1 and mtTFB2, have independent origins in phylogenetically distant clades such that their divergence most likely predates the last universal common ancestor of life. The clade of mtTFB2s comprises orthologs in Opisthokonts and the clade of mtTFB1s includes orthologs in Amoebozoa and Metazoa. Thus, we clearly demonstrate that the neofunctionalization producing the transcription factor function evolved twice independently within the RAMTase family. These results are consistent with and help to elucidate outcomes from prior experimental studies, which found that some members of mtTFB1 still perform the ancestral rRNA methylation function, and the results have broader implications for understanding the evolution of new protein functions. Our phylogenetic reconstruction is also in agreement with prior studies showing two independent origins of plastid RAMTases in Viridiplantae and other photosynthetic autotrophs. We believe that this updated phylogeny of RAMTases should provide a robust evolutionary framework for ongoing studies to identify and characterize the functions of these proteins within diverse organisms.

Arginine deiminase pathway enzymes: evolutionary history in metamonads and other eukaryotes

BACKGROUND

Multiple prokaryotic lineages use the arginine deiminase (ADI) pathway for anaerobic energy production by arginine degradation. The distribution of this pathway among eukaryotes has been thought to be very limited, with only two specialized groups living in low oxygen environments (Parabasalia and Diplomonadida) known to possess the complete set of all three enzymes. We have performed an extensive survey of available sequence data in order to map the distribution of these enzymes among eukaryotes and to reconstruct their phylogenies.

RESULTS

We have found genes for the complete pathway in almost all examined representatives of Metamonada, the anaerobic protist group that includes parabasalids and diplomonads. Phylogenetic analyses indicate the presence of the complete pathway in the last common ancestor of metamonads and heterologous transformation experiments suggest its cytosolic localization in the metamonad ancestor. Outside Metamonada, the complete pathway occurs rarely, nevertheless, it was found in representatives of most major eukaryotic clades.

CONCLUSIONS

Phylogenetic relationships of complete pathways are consistent with the presence of the Archaea-derived ADI pathway in the last common ancestor of all eukaryotes, although other evolutionary scenarios remain possible. The presence of the incomplete set of enzymes is relatively common among eukaryotes and it may be related to the fact that these enzymes are involved in other cellular processes, such as the ornithine-urea cycle. Single protein phylogenies suggest that the evolutionary history of all three enzymes has been shaped by frequent gene losses and horizontal transfers, which may sometimes be connected with their diverse roles in cellular metabolism.

Marine Isolates of Trimastix marina Form a Plesiomorphic Deep-branching Lineage within Preaxostyla, Separate from Other Known Trimastigids (Paratrimastix n. gen.)

Trimastigids are free-living, anaerobic protists that are closely related to the symbiotic oxymonads, forming together the taxon Preaxostyla (Excavata: Metamonada). We isolated fourteen new strains morphologically corresponding to two species assigned to Trimastix (until now the only genus of trimastigids), Trimastix marina and Trimastix pyriformis. Unexpectedly, marine strains of Trimastix marina branch separately from freshwater strains of this morphospecies in SSU rRNA gene trees, and instead form the sister group of all other Preaxostyla. This position is confirmed by three-gene phylogenies. Ultrastructural examination of a marine isolate of Trimastix marina demonstrates a combination of trimastigid-like features (e.g. preaxostyle-like I fibre) and ancestral characters (e.g. absence of thickened flagellar vane margins), consistent with inclusion of marine T. marina within Preaxostyla, but also supporting its distinctiveness from 'freshwater T. marina' and its deep-branching position within Preaxostyla. Since these results indicate paraphyly of Trimastix as currently understood, we transfer the other better-studied trimastigids to Paratrimastix n. gen. and Paratrimastigidae n. fam. The freshwater form previously identified as T. marina is described as Paratrimastix eleionoma n. sp., and Trimastix pyriformis becomes Paratrimastix pyriformis n. comb. Because of its phylogenetic position, 'true' Trimastix is potentially important for understanding the evolution of mitochondrion-related organelles in metamonads.

The role of host phylogeny varies in shaping microbial diversity in the hindguts of lower termites

The hindguts of lower termites and Cryptocercus cockroaches are home to a distinct community of archaea, bacteria, and protists (primarily parabasalids and some oxymonads). Within a host species, the composition of these hindgut communities appears relatively stable, but the evolutionary and ecological factors structuring community composition and stability are poorly understood, as are differential impacts of these factors on protists, bacteria, and archaea. We analyzed the microbial composition of parabasalids and bacteria in the hindguts of Cryptocercus punctulatus and 23 species spanning 4 families of lower termites by pyrosequencing variable regions of the small-subunit rRNA gene. Especially for the parabasalids, these data revealed undiscovered taxa and provided a phylogenetic basis for a more accurate understanding of diversity, diversification, and community composition. The composition of the parabasalid communities was found to be strongly structured by the phylogeny of their hosts, indicating the importance of historical effects, although exceptions were also identified. Particularly, spirotrichonymphids and trichonymphids likely were transferred between host lineages. In contrast, host phylogeny was not sufficient to explain the majority of bacterial community composition, but the compositions of the Bacteroidetes, Elusimicrobia, Tenericutes, Spirochaetes, and Synergistes were structured by host phylogeny perhaps due to their symbiotic associations with protists. All together, historical effects probably resulting from vertical inheritance have had a prominent role in structuring the hindgut communities, especially of the parabasalids, but dispersal and environmental acquisition have played a larger role in community composition than previously expected.

Multigene phylogenies of diverse Carpediemonas-like organisms identify the closest relatives of 'amitochondriate' diplomonads and retortamonads

Diplomonads, retortamonads, and "Carpediemonas-like" organisms (CLOs) are a monophyletic group of protists that are microaerophilic/anaerobic and lack typical mitochondria. Most diplomonads and retortamonads are parasites, and the pathogen Giardia intestinalis is known to possess reduced mitochondrion-related organelles (mitosomes) that do not synthesize ATP. By contrast, free-living CLOs have larger organelles that superficially resemble some hydrogenosomes, organelles that in other protists are known to synthesize ATP anaerobically. This group represents an excellent system for studying the evolution of parasitism and anaerobic, mitochondrion-related organelles. Understanding these evolutionary transitions requires a well-resolved phylogeny of diplomonads, retortamonads and CLOs. Unfortunately, until now the deep relationships amongst these taxa were unresolved due to limited data for almost all of the CLO lineages. To address this, we assembled a dataset of up to six protein-coding genes that includes representatives from all six CLO lineages, and complements existing rRNA datasets. Multigene phylogenetic analyses place CLOs as well as the retortamonad Chilomastix as a paraphyletic basal assemblage to the lineage comprising diplomonads and the retortamonad Retortamonas. In particular, the CLO Dysnectes was shown to be the closest relative of the diplomonads + Retortamonas clade, with strong support. This phylogeny is consistent with a drastic degeneration of mitochondrion-related organelles during the evolution from a free-living organism resembling extant CLOs to a probable parasite/commensal common ancestor of diplomonads and Retortamonas.

Molecular phylogeny and evolution of parabasalia with improved taxon sampling and new protein markers of actin and elongation factor-1α

Background

Inferring the evolutionary history of phylogenetically isolated, deep-branching groups of taxa—in particular determining the root—is often extraordinarily difficult because their close relatives are unavailable as suitable outgroups. One of these taxonomic groups is the phylum Parabasalia, which comprises morphologically diverse species of flagellated protists of ecological, medical, and evolutionary significance. Indeed, previous molecular phylogenetic analyses of members of this phylum have yielded conflicting and possibly erroneous inferences. Furthermore, many species of Parabasalia are symbionts in the gut of termites and cockroaches or parasites and therefore formidably difficult to cultivate, rendering available data insufficient. Increasing the numbers of examined taxa and informative characters (e.g., genes) is likely to produce more reliable inferences.

Principal Findings

Actin and elongation factor-1α genes were identified newly from 22 species of termite-gut symbionts through careful manipulations and seven cultured species, which covered major lineages of Parabasalia. Their protein sequences were concatenated and analyzed with sequences of previously and newly identified glyceraldehyde-3-phosphate dehydrogenase and the small-subunit rRNA gene. This concatenated dataset provided more robust phylogenetic relationships among major groups of Parabasalia and a more plausible new root position than those previously reported.

Conclusions/Significance

We conclude that increasing the number of sampled taxa as well as the addition of new sequences greatly improves the accuracy and robustness of the phylogenetic inference. A morphologically simple cell is likely the ancient form in Parabasalia as opposed to a cell with elaborate flagellar and cytoskeletal structures, which was defined as most basal in previous inferences. Nevertheless, the evolution of Parabasalia is complex owing to several independent multiplication and simplification events in these structures. Therefore, systematics based solely on morphology does not reflect the evolutionary history of parabasalids.

Phylogeny of parasitic parabasalia and free-living relatives inferred from conventional markers vs. Rpb1, a single-copy gene

BACKGROUND

Parabasalia are single-celled eukaryotes (protists) that are mainly comprised of endosymbionts of termites and wood roaches, intestinal commensals, human or veterinary parasites, and free-living species. Phylogenetic comparisons of parabasalids are typically based upon morphological characters and 18S ribosomal RNA gene sequence data (rDNA), while biochemical or molecular studies of parabasalids are limited to a few axenically cultivable parasites. These previous analyses and other studies based on PCR amplification of duplicated protein-coding genes are unable to fully resolve the evolutionary relationships of parabasalids. As a result, genetic studies of Parabasalia lag behind other organisms.

PRINCIPAL FINDINGS

Comparing parabasalid EF1α, α-tubulin, enolase and MDH protein-coding genes with information from the Trichomonas vaginalis genome reveals difficulty in resolving the history of species or isolates apart from duplicated genes. A conserved single-copy gene encodes the largest subunit of RNA polymerase II (Rpb1) in T. vaginalis and other eukaryotes. Here we directly sequenced Rpb1 degenerate PCR products from 10 parabasalid genera, including several T. vaginalis isolates and avian isolates, and compared these data by phylogenetic analyses. Rpb1 genes from parabasalids, diplomonads, Parabodo, Diplonema and Percolomonas were all intronless, unlike intron-rich homologs in Naegleria, Jakoba and Malawimonas.

CONCLUSIONS/SIGNIFICANCE

The phylogeny of Rpb1 from parasitic and free-living parabasalids, and conserved Rpb1 insertions, support Trichomonadea, Tritrichomonadea, and Hypotrichomonadea as monophyletic groups. These results are consistent with prior analyses of rDNA and GAPDH sequences and ultrastructural data. The Rpb1 phylogenetic tree also resolves species- and isolate-level relationships. These findings, together with the relative ease of Rpb1 isolation, make it an attractive tool for evaluating more extensive relationships within Parabasalia.

Patterns of kinesin evolution reveal a complex ancestral eukaryote with a multifunctional cytoskeleton

BACKGROUND

The genesis of the eukaryotes was a pivotal event in evolution and was accompanied by the acquisition of numerous new cellular features including compartmentalization by cytoplasmic organelles, mitosis and meiosis, and ciliary motility. Essential for the development of these features was the tubulin cytoskeleton and associated motors. It is therefore possible to map ancient cell evolution by reconstructing the evolutionary history of motor proteins. Here, we have used the kinesin motor repertoire of 45 extant eukaryotes to infer the ancestral state of this superfamily in the last common eukaryotic ancestor (LCEA).

RESULTS

We bioinformatically identified 1624 putative kinesin proteins, determined their protein domain architectures and calculated a comprehensive Bayesian phylogeny for the kinesin superfamily with statistical support. These data enabled us to define 51 anciently-derived kinesin paralogs (including three new kinesin families) and 105 domain architectures. We then mapped these characters across eukaryotes, accounting for secondary loss within established eukaryotic groupings, and alternative tree topologies.

CONCLUSIONS

We show that a minimum of 11 kinesin families and 3 protein domain architectures were present in the LCEA. This demonstrates that the microtubule-based cytoskeleton of the LCEA was surprisingly highly developed in terms of kinesin motor types, but that domain architectures have been extensively modified during the diversification of the eukaryotes. Our analysis provides molecular evidence for the existence of several key cellular functions in the LCEA, and shows that a large proportion of motor family diversity and cellular complexity had already arisen in this ancient cell.

Evolution and diversity of the Golgi body

Often considered a defining eukaryotic feature, the Golgi body is one of the most recognizable and functionally integrated cellular organelles. It is therefore surprising that some unicellular eukaryotes do not, at first glance, appear to possess Golgi stacks. Here we review the molecular evolutionary, genomic and cell biological evidence for Golgi bodies in these organisms, with the organelle likely present in some form in all cases. This, along with the overwhelming prevalence of stacked cisternae in most eukaryotes, implies that the ancestral eukaryote possessed a stacked Golgi body, with at least eight independent instances of Golgi unstacking in our cellular history.

Phylogenomic analyses support the monophyly of Excavata and resolve relationships among eukaryotic "supergroups"

Nearly all of eukaryotic diversity has been classified into 6 suprakingdom-level groups (supergroups) based on molecular and morphological/cell-biological evidence; these are Opisthokonta, Amoebozoa, Archaeplastida, Rhizaria, Chromalveolata, and Excavata. However, molecular phylogeny has not provided clear evidence that either Chromalveolata or Excavata is monophyletic, nor has it resolved the relationships among the supergroups. To establish the affinities of Excavata, which contains parasites of global importance and organisms regarded previously as primitive eukaryotes, we conducted a phylogenomic analysis of a dataset of 143 proteins and 48 taxa, including 19 excavates. Previous phylogenomic studies have not included all major subgroups of Excavata, and thus have not definitively addressed their interrelationships. The enigmatic flagellate Andalucia is sister to typical jakobids. Jakobids (including Andalucia), Euglenozoa and Heterolobosea form a major clade that we name Discoba. Analyses of the complete dataset group Discoba with the mitochondrion-lacking excavates or "metamonads" (diplomonads, parabasalids, and Preaxostyla), but not with the final excavate group, Malawimonas. This separation likely results from a long-branch attraction artifact. Gradual removal of rapidly-evolving taxa from the dataset leads to moderate bootstrap support (69%) for the monophyly of all Excavata, and 90% support once all metamonads are removed. Most importantly, Excavata robustly emerges between unikonts (Amoebozoa + Opisthokonta) and "megagrouping" of Archaeplastida, Rhizaria, and chromalveolates. Our analyses indicate that Excavata forms a monophyletic suprakingdom-level group that is one of the 3 primary divisions within eukaryotes, along with unikonts and a megagroup of Archaeplastida, Rhizaria, and the chromalveolate lineages.

Flavodiiron protein from Trichomonas vaginalis hydrogenosomes: the terminal oxygen reductase

Trichomonas vaginalis is one of a few eukaryotes that have been found to encode several homologues of flavodiiron proteins (FDPs). Widespread among anaerobic prokaryotes, these proteins are believed to function as oxygen and/or nitric oxide reductases to provide protection against oxidative/nitrosative stresses and host immune responses. One of the T. vaginalis FDP homologues is equipped with a hydrogenosomal targeting sequence and is expressed in the hydrogenosomes, oxygen-sensitive organelles that participate in carbohydrate metabolism and assemble iron-sulfur clusters. The bacterial homologues characterized thus far have been dimers or tetramers; the trichomonad protein is a dimer of identical 45-kDa subunits, each noncovalently binding one flavin mononucleotide. The protein reduces dioxygen to water but is unable to utilize nitric oxide as a substrate, similarly to its closest homologue from another human parasite Giardia intestinalis and related archaebacterial proteins. T. vaginalis FDP is able to accept electrons derived from pyruvate or NADH via ferredoxin and is proposed to play a role in the protection of hydrogenosomes against oxygen.

Environmental PCR survey to determine the distribution of a non-canonical genetic code in uncultivable oxymonads

The universal genetic code is conserved throughout most living systems, but a non-canonical code where TAA and TAG encode glutamine has evolved in several eukaryotes, including oxymonad protists. Most oxymonads are uncultivable, so environmental RT-PCR and PCR was used to examine the distribution of this rare character. A total of 253 unique isolates of four protein-coding genes were sampled from the hindgut community of the cockroach, Cryptocercus punctulatus, an environment rich in diversity from two of the five subgroups of oxymonad, saccinobaculids and polymastigids. Four alpha-tubulins were found with non-canonical glutamine codons. Environmental RACE confirmed that these and related genes used only TGA as stop codons, as expected for the non-canonical code, whereas other genes used TAA or TAG as stop codons, as expected for the universal code. We characterized alpha-tubulin from manually isolated Saccinobaculus ambloaxostylus, confirming it uses the universal code and suggesting, by elimination, that the non-canonical code is used by a polymastigid. HSP90 and EF-1alpha phylogenies also showed environmental sequences falling into two distinct groups, and are generally consistent with previous hypotheses that polymastigids and Streblomastix are closely related. Overall, we propose that the non-canonical genetic code arose once in a common ancestor of Streblomastix and a subgroup of polymastigids.

Ancestral roles of eukaryotic frataxin: mitochondrial frataxin function and heterologous expression of hydrogenosomal Trichomonas homologues in trypanosomes

Frataxin is a small conserved mitochondrial protein; in humans, mutations affecting frataxin expression or function result in Friedreich's ataxia. Much of the current understanding of frataxin function comes from informative studies with yeast models, but considerable debates remain with regard to the primary functions of this ubiquitous protein. We exploit the tractable reverse genetics of Trypanosoma brucei in order to specifically consider the importance of frataxin in an early branching lineage. Using inducible RNAi, we show that frataxin is essential in T. brucei and that its loss results in reduced activity of the marker Fe-S cluster-containing enzyme aconitase in both the mitochondrion and cytosol. Activities of mitochondrial succinate dehydrogenase and fumarase also decreased, but the concentration of reactive oxygen species increased. Trypanosomes lacking frataxin also exhibited a low mitochondrial membrane potential and reduced oxygen consumption. Crucially, however, iron did not accumulate in frataxin-depleted mitochondria, and as T. brucei frataxin does not form large complexes, it suggests that it plays no role in iron storage. Interestingly, RNAi phenotypes were ameliorated by expression of frataxin homologues from hydrogenosomes of another divergent protist Trichomonas vaginalis. Collectively, the data suggest trypanosome frataxin functions primarily only in Fe-S cluster biogenesis and protection from reactive oxygen species.

Order within a mosaic distribution of mitochondrial c-type cytochrome biogenesis systems?

Mitochondrial cytochromes c and c(1) are present in all eukaryotes that use oxygen as the terminal electron acceptor in the respiratory chain. Maturation of c-type cytochromes requires covalent attachment of the heme cofactor to the protein, and there are at least five distinct biogenesis systems that catalyze this post-translational modification in different organisms and organelles. In this study, we use biochemical data, comparative genomic and structural bioinformatics investigations to provide a holistic view of mitochondrial c-type cytochrome biogenesis and its evolution. There are three pathways for mitochondrial c-type cytochrome maturation, only one of which is present in prokaryotes. We analyze the evolutionary distribution of these biogenesis systems, which include the Ccm system (System I) and the enzyme heme lyase (System III). We conclude that heme lyase evolved once and, in many lineages, replaced the multicomponent Ccm system (present in the proto-mitochondrial endosymbiont), probably as a consequence of lateral gene transfer. We find no evidence of a System III precursor in prokaryotes, and argue that System III is incompatible with multi-heme cytochromes common to bacteria, but absent from eukaryotes. The evolution of the eukaryotic-specific protein heme lyase is strikingly unusual, given that this protein provides a function (thioether bond formation) that is also ubiquitous in prokaryotes. The absence of any known c-type cytochrome biogenesis system from the sequenced genomes of various trypanosome species indicates the presence of a third distinct mitochondrial pathway. Interestingly, this system attaches heme to mitochondrial cytochromes c that contain only one cysteine residue, rather than the usual two, within the heme-binding motif. The isolation of single-cysteine-containing mitochondrial cytochromes c from free-living kinetoplastids, Euglena and the marine flagellate Diplonema papillatum suggests that this unique form of heme attachment is restricted to, but conserved throughout, the protist phylum Euglenozoa.

Genetic evidence for a mitochondriate ancestry in the 'amitochondriate' flagellate Trimastix pyriformis

Most modern eukaryotes diverged from a common ancestor that contained the α-proteobacterial endosymbiont that gave rise to mitochondria. The ‘amitochondriate’ anaerobic protist parasites that have been studied to date, such as Giardia and Trichomonas harbor mitochondrion-related organelles, such as mitosomes or hydrogenosomes. Yet there is one remaining group of mitochondrion-lacking flagellates known as the Preaxostyla that could represent a primitive ‘pre-mitochondrial’ lineage of eukaryotes. To test this hypothesis, we conducted an expressed sequence tag (EST) survey on the preaxostylid flagellate Trimastix pyriformis, a poorly-studied free-living anaerobe. Among the ESTs we detected 19 proteins that, in other eukaryotes, typically function in mitochondria, hydrogenosomes or mitosomes, 12 of which are found exclusively within these organelles. Interestingly, one of the proteins, aconitase, functions in the tricarboxylic acid cycle typical of aerobic mitochondria, whereas others, such as pyruvate:ferredoxin oxidoreductase and [FeFe] hydrogenase, are characteristic of anaerobic hydrogenosomes. Since Trimastix retains genetic evidence of a mitochondriate ancestry, we can now say definitively that all known living eukaryote lineages descend from a common ancestor that had mitochondria.

Implications of the new eukaryotic systematics for parasitologists

An accurate understanding of evolutionary relationships is central in biology. For parasitologists, understanding the relationships among eukaryotic organisms allows the prediction of virulence mechanisms, reconstruction of metabolic pathways, identification of potential drug targets, elucidation of parasite-specific cellular processes and understanding of interactions with the host or vector. Here we consider the impact of major recent revisions of eukaryotic systematics and taxonomy on parasitology. The previous, ladder-like model placed some protists as early diverging, with the remaining eukaryotes "progressing" towards a "crown radiation" of animals, plants, Fungi and some additional protistan lineages. This model has been robustly disproven. The new model is based on vastly increased amounts of molecular sequence data, integration with morphological information and the rigorous application of phylogenetic methods to those data. It now divides eukaryotes into six major supergroups; the relationships between those groups and the order of branching remain unknown. This new eukaryotic phylogeny emphasizes that organisms including Giardia, Trypanosoma and Trichomonas are not primitive, but instead highly evolved and specialised for their specific environments. The wealth of newly available comparative genomic data has also allowed the reconstruction of ancient suites of characteristics and mapping of character evolution in diverse parasites. For example, the last common eukaryotic ancestor was apparently complex, suggesting that lineage-specific adaptations and secondary losses have been important in the evolution of protistan parasites. Referring to the best evidence-based models for eukaryotic evolution will allow parasitologists to make more accurate and reliable inferences about pathogens that cause significant morbidity and mortality.

An expanded inventory of conserved meiotic genes provides evidence for sex in Trichomonas vaginalis

Meiosis is a defining feature of eukaryotes but its phylogenetic distribution has not been broadly determined, especially among eukaryotic microorganisms (i.e. protists)-which represent the majority of eukaryotic 'supergroups'. We surveyed genomes of animals, fungi, plants and protists for meiotic genes, focusing on the evolutionarily divergent parasitic protist Trichomonas vaginalis. We identified homologs of 29 components of the meiotic recombination machinery, as well as the synaptonemal and meiotic sister chromatid cohesion complexes. T. vaginalis has orthologs of 27 of 29 meiotic genes, including eight of nine genes that encode meiosis-specific proteins in model organisms. Although meiosis has not been observed in T. vaginalis, our findings suggest it is either currently sexual or a recent asexual, consistent with observed, albeit unusual, sexual cycles in their distant parabasalid relatives, the hypermastigotes. T. vaginalis may use meiotic gene homologs to mediate homologous recombination and genetic exchange. Overall, this expanded inventory of meiotic genes forms a useful "meiosis detection toolkit". Our analyses indicate that these meiotic genes arose, or were already present, early in eukaryotic evolution; thus, the eukaryotic cenancestor contained most or all components of this set and was likely capable of performing meiotic recombination using near-universal meiotic machinery.

Where do animal α-amylases come from? An interkingdom trip

Alpha-amylases are widely found in eukaryotes and prokaryotes. Few amino acids are conserved among these organisms, but at an intra-kingdom level, conserved protein domains exist. In animals, numerous conserved stretches are considered as typical of animal alpha-amylases. Searching databases, we found no animal-type alpha-amylases outside the Bilateria. Instead, we found in the sponge Reniera sp. and in the sea anemone Nematostella vectensis, alpha-amylases whose most similar cognate was that of the amoeba Dictyostelium discoideum. We found that this "Dictyo-type" alpha-amylase was shared not only by these non-Bilaterian animals, but also by other Amoebozoa, Choanoflagellates, and Fungi. This suggested that the Dictyo-type alpha-amylase was present in the last common ancestor of Unikonts. The additional presence of the Dictyo-type in some Ciliates and Excavates, suggests that horizontal gene transfers may have occurred among Eukaryotes. We have also detected putative interkingdom transfers of amylase genes, which obscured the historical reconstitution. Several alternative scenarii are discussed.

A genomic survey of the fish parasite Spironucleus salmonicida indicates genomic plasticity among diplomonads and significant lateral gene transfer in eukaryote genome evolution

BACKGROUND

Comparative genomic studies of the mitochondrion-lacking protist group Diplomonadida (diplomonads) has been lacking, although Giardia lamblia has been intensively studied. We have performed a sequence survey project resulting in 2341 expressed sequence tags (EST) corresponding to 853 unique clones, 5275 genome survey sequences (GSS), and eleven finished contigs from the diplomonad fish parasite Spironucleus salmonicida (previously described as S. barkhanus).

RESULTS

The analyses revealed a compact genome with few, if any, introns and very short 3' untranslated regions. Strikingly different patterns of codon usage were observed in genes corresponding to frequently sampled ESTs versus genes poorly sampled, indicating that translational selection is influencing the codon usage of highly expressed genes. Rigorous phylogenomic analyses identified 84 genes--mostly encoding metabolic proteins--that have been acquired by diplomonads or their relatively close ancestors via lateral gene transfer (LGT). Although most acquisitions were from prokaryotes, more than a dozen represent likely transfers of genes between eukaryotic lineages. Many genes that provide novel insights into the genetic basis of the biology and pathogenicity of this parasitic protist were identified including 149 that putatively encode variant-surface cysteine-rich proteins which are candidate virulence factors. A number of genomic properties that distinguish S. salmonicida from its human parasitic relative G. lamblia were identified such as nineteen putative lineage-specific gene acquisitions, distinct mutational biases and codon usage and distinct polyadenylation signals.

CONCLUSION

Our results highlight the power of comparative genomic studies to yield insights into the biology of parasitic protists and the evolution of their genomes, and suggest that genetic exchange between distantly-related protist lineages may be occurring at an appreciable rate in eukaryote genome evolution.

The glycolytic pathway of Trimastix pyriformis is an evolutionary mosaic

Background

Glycolysis and subsequent fermentation is the main energy source for many anaerobic organisms. The glycolytic pathway consists of ten enzymatic steps which appear to be universal amongst eukaryotes. However, it has been shown that the origins of these enzymes in specific eukaryote lineages can differ, and sometimes involve lateral gene transfer events. We have conducted an expressed sequence tag (EST) survey of the anaerobic flagellate Trimastix pyriformis to investigate the nature of the evolutionary origins of the glycolytic enzymes in this relatively unstudied organism.

Results

We have found genes in the Trimastix EST data that encode enzymes potentially catalyzing nine of the ten steps of the glycolytic conversion of glucose to pyruvate. Furthermore, we have found two different enzymes that in principle could catalyze the conversion of phosphoenol pyruvate (PEP) to pyruvate (or the reverse reaction) as part of the last step in glycolysis. Our phylogenetic analyses of all of these enzymes revealed at least four cases where the relationship of the Trimastix genes to homologs from other species is at odds with accepted organismal relationships. Although lateral gene transfer events likely account for these anomalies, with the data at hand we were not able to establish with confidence the bacterial donor lineage that gave rise to the respective Trimastix enzymes.

Conclusion

A number of the glycolytic enzymes of Trimastix have been transferred laterally from bacteria instead of being inherited from the last common eukaryotic ancestor. Thus, despite widespread conservation of the glycolytic biochemical pathway across eukaryote diversity, in a number of protist lineages the enzymatic components of the pathway have been replaced by lateral gene transfer from disparate evolutionary sources. It remains unclear if these replacements result from selectively advantageous properties of the introduced enzymes or if they are neutral outcomes of a gene transfer 'ratchet' from food or endosymbiotic organisms or a combination of both processes.

Evaluating support for the current classification of eukaryotic diversity

Perspectives on the classification of eukaryotic diversity have changed rapidly in recent years, as the four eukaryotic groups within the five-kingdom classification—plants, animals, fungi, and protists—have been transformed through numerous permutations into the current system of six “supergroups.” The intent of the supergroup classification system is to unite microbial and macroscopic eukaryotes based on phylogenetic inference. This supergroup approach is increasing in popularity in the literature and is appearing in introductory biology textbooks. We evaluate the stability and support for the current six-supergroup classification of eukaryotes based on molecular genealogies. We assess three aspects of each supergroup: (1) the stability of its taxonomy, (2) the support for monophyly (single evolutionary origin) in molecular analyses targeting a supergroup, and (3) the support for monophyly when a supergroup is included as an out-group in phylogenetic studies targeting other taxa. Our analysis demonstrates that supergroup taxonomies are unstable and that support for groups varies tremendously, indicating that the current classification scheme of eukaryotes is likely premature. We highlight several trends contributing to the instability and discuss the requirements for establishing robust clades within the eukaryotic tree of life.

Evolutionary perspectives, including the classification of living organisms, provide the unifying scaffold on which biological knowledge is assembled. Researchers in many areas of biology use evolutionary classifications (taxonomy) in many ways, including as a means for interpreting the origin of evolutionary innovations, as a framework for comparative genetics/genomics, and as the basis for drawing broad conclusions about the diversity of living organisms. Thus, it is essential that taxonomy be robust. Here the authors evaluate the stability of and support for the current classification system of eukaryotic cells (cells with nuclei) in which eukaryotes are divided into six kingdom level categories, or supergroups. These six supergroups unite diverse microbial and macrobial eukaryotic lineages, including the well-known groups of plants, animals, and fungi. The authors assess the stability of supergroup classifications through time and reveal a rapidly changing taxonomic landscape that is difficult to navigate for the specialist and generalist alike. Additionally, the authors find variable support for each of the supergroups in published analyses based on DNA sequence variation. The support for supergroups differs according to the taxonomic area under study and the origin of the genes (e.g., nuclear, plastid) used in the analysis. Encouragingly, combining a conservative approach to taxonomy with increased sampling of microbial eukaryotes and the use of multiple types of data is likely to produce a robust scaffold for the eukaryotic tree of life.

Reconstructing the mosaic glycolytic pathway of the anaerobic eukaryote Monocercomonoides

All eukaryotes carry out glycolysis, interestingly, not all using the same enzymes. Anaerobic eukaryotes face the challenge of fewer molecules of ATP extracted per molecule of glucose due to their lack of a complete tricarboxylic acid cycle. This may have pressured anaerobic eukaryotes to acquire the more ATP-efficient alternative glycolytic enzymes, such as pyrophosphate-fructose 6-phosphate phosphotransferase and pyruvate orthophosphate dikinase, through lateral gene transfers from bacteria and other eukaryotes. Most studies of these enzymes in eukaryotes involve pathogenic anaerobes; Monocercomonoides, an oxymonad belonging to the eukaryotic supergroup Excavata, is a nonpathogenic anaerobe representing an evolutionarily and ecologically distinct sampling of an anaerobic glycolytic pathway. We sequenced cDNA encoding glycolytic enzymes from a previously established cDNA library of Monocercomonoides and analyzed the relationships of these enzymes to those from other organisms spanning the major groups of Eukaryota, Bacteria, and Archaea. We established that, firstly, Monocercomonoides possesses alternative versions of glycolytic enzymes: fructose-6-phosphate phosphotransferase, both pyruvate kinase and pyruvate orthophosphate dikinase, cofactor-independent phosphoglycerate mutase, and fructose-bisphosphate aldolase (class II, type B). Secondly, we found evidence for the monophyly of oxymonads, kinetoplastids, diplomonads, and parabasalids, the major representatives of the Excavata. We also found several prokaryote-to-eukaryote as well as eukaryote-to-eukaryote lateral gene transfers involving glycolytic enzymes from anaerobic eukaryotes, further suggesting that lateral gene transfer was an important factor in the evolution of this pathway for denizens of this environment.

Evolutionary Relationships of Apusomonads Inferred from Taxon-Rich Analyses of 6 Nuclear Encoded Genes

The phylogenetic relationships of the biflagellate protist group Apusomonadidae have been unclear despite the availability of some molecular data. We analyzed sequences from 6 nuclear encoded genes-small-subunit rRNA, large-subunit rRNA, alpha-tubulin, beta-tubulin, actin, and heat shock protein 90-to infer the phylogenetic position of Apusomonas proboscidea Aléxéieff 1924. To increase the taxon richness of the study, we also obtained new sequences from representatives of several other major eukaryotic groups: Chrysochromulina sp. National Institute for Environmental Studies 1333 (Haptophyta), Cyanophora paradoxa (Glaucophyta), Goniomonas truncata (Cryptophyceae), Leucocryptos marina (Kathablepharidae), Mesostigma viride (Streptophyta, Viridiplantae), Peridinium limbatum (Alveolata), Pterosperma cristatum (Prasinophytae, Viridiplantae), Synura sphagnicola (Stramenopiles), and Thaumatomonas sp. (Rhizaria). In most individual gene phylogenies, Apusomonas branched close to either of the 2 related taxa-Opisthokonta (including animals, fungi, and choanoflagellates) or Amoebozoa. Combined analyses of all 4 protein-coding genes or all 6 studied genes strongly supported the hypothesis that Apusomonadidae is closely related to Opisthokonta (or to all other eukaryotic groups except Opisthokonta, depending on the position of the eukaryotic root). Alternative hypotheses were rejected in approximately unbiased tests at the 5% level. However, the strong phylogenetic signal supporting a specific affiliation between Apusomonadidae and Opisthokonta largely originated from the alpha-tubulin data. If alpha-tubulin is not considered, topologies in which Apusomonadidae is sister to Opisthokonta or is sister to Amoebozoa were more or less equally supported. One current model for deep eukaryotic evolution holds that eukaryotes are divided into primary "unikont" and "bikont" clades and are descended from a "uniflagellate" common ancestor. Together with other information, our data suggest instead that unikonts (=Opisthokonta and Amoebozoa) are not strictly monophyletic and are descended from biflagellate ancestors.

The origin and diversification of eukaryotes: problems with molecular phylogenetics and molecular clock estimation

Determining the relationships among and divergence times for the major eukaryotic lineages remains one of the most important and controversial outstanding problems in evolutionary biology. The sequencing and phylogenetic analyses of ribosomal RNA (rRNA) genes led to the first nearly comprehensive phylogenies of eukaryotes in the late 1980s, and supported a view where cellular complexity was acquired during the divergence of extant unicellular eukaryote lineages. More recently, however, refinements in analytical methods coupled with the availability of many additional genes for phylogenetic analysis showed that much of the deep structure of early rRNA trees was artefactual. Recent phylogenetic analyses of a multiple genes and the discovery of important molecular and ultrastructural phylogenetic characters have resolved eukaryotic diversity into six major hypothetical groups. Yet relationships among these groups remain poorly understood because of saturation of sequence changes on the billion-year time-scale, possible rapid radiations of major lineages, phylogenetic artefacts and endosymbiotic or lateral gene transfer among eukaryotes. Estimating the divergence dates between the major eukaryote lineages using molecular analyses is even more difficult than phylogenetic estimation. Error in such analyses comes from a myriad of sources including: (i) calibration fossil dates, (ii) the assumed phylogenetic tree, (iii) the nucleotide or amino acid substitution model, (iv) substitution number (branch length) estimates, (v) the model of how rates of evolution change over the tree, (vi) error inherent in the time estimates for a given model and (vii) how multiple gene data are treated. By reanalysing datasets from recently published molecular clock studies, we show that when errors from these various sources are properly accounted for, the confidence intervals on inferred dates can be very large. Furthermore, estimated dates of divergence vary hugely depending on the methods used and their assumptions. Accurate dating of divergence times among the major eukaryote lineages will require a robust tree of eukaryotes, a much richer Proterozoic fossil record of microbial eukaryotes assignable to extant groups for calibration, more sophisticated relaxed molecular clock methods and many more genes sampled from the full diversity of microbial eukaryotes.

Telonemia, a new protist phylum with affinity to chromist lineages

Recent molecular investigations of marine samples taken from different environments, including tropical, temperate and polar areas, as well as deep thermal vents, have revealed an unexpectedly high diversity of protists, some of them forming deep-branching clades within important lineages, such as the alveolates and heterokonts. Using the same approach on coastal samples, we have identified a novel group of protist small subunit (SSU) rDNA sequences that do not correspond to any phylogenetic group previously identified. Comparison with other sequences obtained from cultures of heterotrophic protists showed that the environmental sequences grouped together with Telonema, a genus known since 1913 but of uncertain taxonomic affinity. Phylogenetic analyses using four genes (SSU, Hsp90, alpha-tubulin and beta-tubulin), and accounting for gamma- and covarion-distributed substitution rates, revealed Telonema as a distinct group of species branching off close to chromist lineages. Consistent with these gene trees, Telonema possesses ultrastructures revealing both the distinctness of the group and the evolutionary affinity to chromist groups. Altogether, the data suggest that Telonema constitutes a new eukaryotic phylum, here defined as Telonemia, possibly representing a key clade for the understanding of the early evolution of bikont protist groups, such as the proposed chromalveolate supergroup.

The phylogenetic position of the oxymonad Saccinobaculus based on SSU rRNA

The oxymonads are a group of structurally complex anaerobic flagellates about which we know very little. They are found in association with complex microbial communities in the guts of animals. There are five recognized families of oxymonads; molecular data have been acquired for four of these. Here, we describe the first molecular data from the last remaining group, represented by Saccinobaculus, an organism that is found exclusively in the hindgut of the wood-eating cockroach Cryptocercus. We sequenced small subunit ribosomal RNA (SSU rRNA) from total gut DNA to describe Saccinobaculus SSU rRNA diversity. We also sequenced SSU rRNA from manually isolated cells of the two most abundant and readily identifiable species: the type species Saccinobaculus ambloaxostylus and the taxonomically contentious Saccinobaculus doroaxostylus. We inferred phylogenetic trees including all five known oxymonad subgroups in order to elucidate the internal phylogeny of this poorly studied group, to resolve some outstanding issues of the taxonomy and identification of certain Saccinobaculus species, and to investigate the evolution of character states within it. Our analysis recovered strong support for the existence of the five subgroups of oxymonads, and consistently united the subgroups containing Monocercomonoides and Streblomastix, but was unable to resolve any further higher-order branching patterns.

Monophyly of Rhizaria and multigene phylogeny of unicellular bikonts

Reconstructing a global phylogeny of eukaryotes is an ongoing challenge of molecular phylogenetics. The availability of genomic data from a broad range of eukaryotic phyla helped in resolving the eukaryotic tree into a topology with a rather small number of large assemblages, but the relationships between these "supergroups" are yet to be confirmed. Rhizaria is the most recently recognized "supergroup," but, in spite of this important position within the tree of life, their representatives are still missing in global phylogenies of eukaryotes. Here, we report the first large-scale analysis of eukaryote phylogeny including data for 2 rhizarian species, the foraminiferan Reticulomyxa filosa and the chlorarachniophyte Bigelowiella natans. Our results confirm the monophyly of Rhizaria (Foraminifera + Cercozoa), with very high bootstrap supports in all analyses. The overall topology of our trees is in agreement with the current view of eukaryote phylogeny with basal division into "unikonts" (Opisthokonts and Ameobozoa) and "bikonts" (Plantae, alveolates, stramenopiles, and excavates). As expected, Rhizaria branch among bikonts; however, their phylogenetic position is uncertain. Depending on the data set and the type of analysis, Rhizaria branch as sister group to either stramenopiles or excavates. Overall, the relationships between the major groups of unicellular bikonts are poorly resolved, despite the use of 85 proteins and the largest taxonomic sampling for this part of the tree available to date. This may be due to an acceleration of evolutionary rates in some bikont phyla or be related to their rapid diversification in the early evolution of eukaryotes.

A high density of ancient spliceosomal introns in oxymonad excavates

Background

Certain eukaryotic genomes, such as those of the amitochondriate parasites Giardia and Trichomonas, have very low intron densities, so low that canonical spliceosomal introns have only recently been discovered through genome sequencing. These organisms were formerly thought to be ancient eukaryotes that diverged before introns originated, or at least became common. Now however, they are thought to be members of a supergroup known as excavates, whose members generally appear to have low densities of canonical introns. Here we have used environmental expressed sequence tag (EST) sequencing to identify 17 genes from the uncultivable oxymonad Streblomastix strix, to survey intron densities in this most poorly studied excavate group.

Results

We find that Streblomastix genes contain an unexpectedly high intron density of about 1.1 introns per gene. Moreover, over 50% of these are at positions shared between a broad spectrum of eukaryotes, suggesting theyare very ancient introns, potentially present in the last common ancestor of eukaryotes.

Conclusion

The Streblomastix data show that the genome of the ancestor of excavates likely contained many introns and the subsequent evolution of introns has proceeded very differently in different excavate lineages: in Streblomastix there has been much stasis while in Trichomonas and Giardia most introns have been lost.

Evolution of four gene families with patchy phylogenetic distributions: influx of genes into protist genomes

BACKGROUND

Lateral gene transfer (LGT) in eukaryotes from non-organellar sources is a controversial subject in need of further study. Here we present gene distribution and phylogenetic analyses of the genes encoding the hybrid-cluster protein, A-type flavoprotein, glucosamine-6-phosphate isomerase, and alcohol dehydrogenase E. These four genes have a limited distribution among sequenced prokaryotic and eukaryotic genomes and were previously implicated in gene transfer events affecting eukaryotes. If our previous contention that these genes were introduced by LGT independently into the diplomonad and Entamoeba lineages were true, we expect that the number of putative transfers and the phylogenetic signal supporting LGT should be stable or increase, rather than decrease, when novel eukaryotic and prokaryotic homologs are added to the analyses.

RESULTS

The addition of homologs from phagotrophic protists, including several Entamoeba species, the pelobiont Mastigamoeba balamuthi, and the parabasalid Trichomonas vaginalis, and a large quantity of sequences from genome projects resulted in an apparent increase in the number of putative transfer events affecting all three domains of life. Some of the eukaryotic transfers affect a wide range of protists, such as three divergent lineages of Amoebozoa, represented by Entamoeba, Mastigamoeba, and Dictyostelium, while other transfers only affect a limited diversity, for example only the Entamoeba lineage. These observations are consistent with a model where these genes have been introduced into protist genomes independently from various sources over a long evolutionary time.

CONCLUSION

Phylogenetic analyses of the updated datasets using more sophisticated phylogenetic methods, in combination with the gene distribution analyses, strengthened, rather than weakened, the support for LGT as an important mechanism affecting the evolution of these gene families. Thus, gene transfer seems to be an on-going evolutionary mechanism by which genes are spread between unrelated lineages of all three domains of life, further indicating the importance of LGT from non-organellar sources into eukaryotic genomes.

Pyruvate-phosphate dikinase of oxymonads and parabasalia and the evolution of pyrophosphate-dependent glycolysis in anaerobic eukaryotes

In pyrophosphate-dependent glycolysis, the ATP/ADP-dependent enzymes phosphofructokinase (PFK) and pyruvate kinase are replaced by the pyrophosphate-dependent PFK and pyruvate phosphate dikinase (PPDK), respectively. This variant of glycolysis is widespread among bacteria, but it also occurs in a few parasitic anaerobic eukaryotes such as Giardia and Entamoeba spp. We sequenced two genes for PPDK from the amitochondriate oxymonad Streblomastix strix and found evidence for PPDK in Trichomonas vaginalis and other parabasalia, where this enzyme was thought to be absent. The Streblomastix and Giardia genes may be related to one another, but those of Entamoeba and perhaps Trichomonas are distinct and more closely related to bacterial homologues. These findings suggest that pyrophosphate-dependent glycolysis is more widespread in eukaryotes than previously thought, enzymes from the pathway coexists with ATP-dependent more often than previously thought and may be spread by lateral transfer of genes for pyrophosphate-dependent enzymes from bacteria.