BACKGROUND
eukaryotes catabolize glucose to pyruvate via the embden-meyerhof-parnas  glycolytic pathway  <cit> .
a number of protists are known that lack classical mitochondrial structures and electron-transport-linked atp production, instead relying on substrate-level phosphorylation.
in these organisms, glycolysis and subsequent fermentation reactions play a central role in their energy metabolism  <cit> .
although often referred to as 'amitochondriates', most of these protist taxa have been shown to harbor organelles of mitochondrial descent  <cit> .
hydrogenosomes for example metabolize pyruvate via a fermentative pathway that yields atp and hydrogen  <cit> .
other mitochondrion-derived organelles, such as the 'mitosomes' of the parasites entamoeba histolytica and giardia lamblia, seem to have lost all core metabolic functions  <cit> .
in these organisms, pyruvate is catabolized by fermentation in the cytosol.
although the core emp pathway appears to be almost universal amongst eukaryotes, amitochondriate protists have been shown to be surprisingly flexible in terms of the specific enzyme families they utilize to carry out glycolytic reactions.
more interestingly these enzymes appear to have rather diverse origins  <cit>  and it is likely that this variability was achieved, at least in part, by lateral gene transfer   <cit> .
for instance, in the first step of glycolysis, diplomonads and parabasalids use glucokinase, rather than hexokinase which carries out this reaction in most other eukaryotes  <cit> .
similarly, amitochondriate protists such as parabasalids, diplomonads, pelobionts and entamoebids cleave fructose- <dig> -bisphosphate with a class ii fructose-bisphosphate aldolase , while most other eukaryotes possess a class i fba  <cit> .
these two enzymes are non-homologous and belong to different superfamilies  <cit> .
in many cases, the glycolytic enzymes that are replaced by lgt appear to come from a bacterial donor lineage.
for instance, diplomonad and parabasalid glucose phosphate isomerases  are more closely related to the gpis of cyanobacteria and chloroplasts than to cytosolic gpis of other eukaryotes  <cit> .
the glyceraldehyde-3-phosphate dehydrogenase  of parabasalids is not part of the gapc clade which encompasses most eukaryotic gapdhs including those of most amitochondriates.
instead, it clusters most closely with the sequences of the spirochaete genus borrelia within the bacterial gapab clade  <cit> .
similarly, the evolution of phosphofructokinase  shows phospho-donor changes and frequent lateral gene transfer  <cit> .
in amitochondriates, pyrophosphate-dependent, atp-dependent or both types of pfk are present – most of them acquired by lgt  <cit> .
in some protist taxa, additional enzymes involved in glycolytic reactions are found that are rare amongst most eukaryotes.
for instance, pyruvate phosphate dikinase  is an analogous enzyme to pyruvate kinase , with a key difference being that ppdk uses pyrophosphate instead of atp to catalyze the conversion of phosphoenolpyruvate into pyruvate and generates atp in the process.
ppdk is present in a number of eukaryotes including the anaerobes entamoeba histolytica, giardia lamblia  <cit> , trichomonas vaginalis and streblomastix strix  <cit>  and the phylogenetic affinites of the eukaryote enzymes are diverse  <cit> .
these examples demonstrate the diversity in the equipment of glycolytic enzymes and in their evolutionary history in protists in general, and anaerobic species in particular.
however, our knowledge of the biochemistry of diverse protist lineages still remains severely limited to taxa of biomedical or agricultural importance.
here we expand this knowledge to include the anaerobic flagellate trimastix pyriformis  <cit> .
species of the genus trimastix are free-living nanoflagellates which thrive in anoxic environments, where they feed on bacteria  <cit> .
they harbor small double membrane-bound organelles of unknown function which structurally resemble hydrogenosomes  <cit>  but are probably derived from mitochondria  <cit> .
phylogenetic analyses indicate that trimastix species form a sister group to the oxymonads, another goup of 'amitochondriate' protist that are symbionts of termites and cockroaches  <cit> .
the placement of these organisms in the eukaryotic tree remains contentious, but on the basis of ultrastructural comparisons they have been suggested to be part of a eukaryote super-group called the excavata  <cit> .
the excavata is comprised of a number of aerobic protist lineages including the heterolobosea, euglenozoa  and jakobid flagellates as well as anaerobic lineages such as diplomonads , carpediemonas and parabasalids .
although, ultrastructural evidence for the excavata is convincing  <cit> , molecular phylogenetic support remains weak  <cit> .
to improve our knowledge of the biochemical diversity of anaerobic protists we initiated an est project of trimastix pyriformis.
we identified a large number of ests coding for enzymes involved in glycolysis and subsequent fermentation.
our phylogenetic analyses reveal that, while some glycolytic enzymes in this organism were inherited by vertical descent, about half of the glycolytic enzymes of trimastix were acquired by lateral gene transfer events, likely from different bacterial donors.
although the glycolytic pathway is universally conserved in eukaryotes, its enzyme components are evolutionarily labile and have been repeatedly replaced in separate eukaryotic lineages by gene transfer from diverse eubacterial donors.
RESULTS
phylogenies of glycolytic enzymes of trimastix pyriformis
from our est data we identified  <dig> homologs of glycolytic enzymes typically found in other eukaryotes plus a second enzyme for the final glycolytic step ).
in order to completely sequence each homolog, 5'-race and primer walking sequencing were performed.
putative enzyme activities, lengths, and gc contents of the encoded proteins are shown in table  <dig>  in the following sections we describe the properties and phylogenetic affinities of each enzyme in turn.
hexokinase
the evolution of hexokinases is characterized by several independent gene duplications which gave rise to numerous isoenzymes within different phyla: plants have two isoenzymes , vertebrates have four  and yeasts have two   <cit> ; the vertebrate hexokinase d is often referred to as glucokinase.
like most aerobic eukaryotes and the amitochondriate parasite entamoeba histolytica, trimastix appears to use hexokinase for the first enzymatic step in glycolysis whereas the amitochondriate 'excavates' giardia intestinalis, spironucleus barkhanus and trichomonas vaginalis all appear to use glucokinase  <cit> .
it is uncertain if bacteria use a real hexokinase, however search of several databases revealed three prokaryotes  with hexokinases that were readily aligned with the eukaryotic hexokinase genes.
these three prokaryotic sequences form a distinct clade that is clearly separated from the eukaryotic hexokinase sequences .
the eukaryotic part of the tree shows very little resolution between major groups with fungi appearing to be paraphyletic with respect to bacteria and other eukaryotes, a feature that is not supported by bootstrap analysis.
the trimastix homolog emerges firmly within the eukaryotes, but is not strongly allied with any other lineages.
thus the trimastix hexokinase gene is most probably of eukaryotic origin via vertical descent.
glucose phosphate isomerase
the gpi gene family can be divided into three major types of gpis which show only a low degree of sequence similarity  <cit> : type i gpis are found in the eukaryotic cytosol and in many bacteria.
type ii are typical for cyanobacteria and chloroplasts but were recently also detected in the amitochondriate excavates giardia intestinalis, spironucleus barkhanus, and trichomonas vaginalis  <cit> .
finally, the gpis of the archaeon methanococcus jannaschii and some gram positive bacteria are assigned to type iii  <cit> .
for tree reconstruction, we used only gpis of type i and ii .
the gpi of trimastix pyriformis is of type i and not similar to the type ii gpis of diplomonads and parabasalids.
in the type i subtree, several bacterial lineages and two clades of eukaryotic gpis can be distinguished : one of the clades consists of cytosolic gpis from alveolates, stramenopiles, and green plants.
the second comprises cytosolic gpis of fungi and animals, glycosomal gpis of kinetoplastids, and a clade of bacterial gpis .
likelihood comparisons of different possible trees neither supported nor rejected this seemingly robust grouping of animals and bacteria  <cit>  and it was concluded that the branching pattern within the type i subtree is very unstable and currently unresolved.
except for the most basal branches, this also describes our results.
the trimastix gpi falls in the unresolved part of the type i subtree  where it groups with the alpha-proteobacterium agrobacterium tumefaciens  but with no statistical support.
given the poor resolution of the type i gpi subtree, we cannot trace the origin of the trimastix gpi other than it is quite distinct from the type ii forms of the other amitochondriate flagellates .
fructose- <dig> -bisphosphate aldolase
a phylogeny of type b class ii fba enzymes is shown in fig 1a.
trimastix pyriformis possesses a type b class ii fba like the other amitochondriates  <cit>  and branches robustly with diplomonad fba sequences in our phylogenetic analysis .
the other amitochondriates form a clearly separated clade branching within diverse bacteria.
an alignment of class ii type b fba protein sequences supports the phylogenetic separation of these groups of amitochondriate taxa; the fbas of the diplomonads and trimastix share some specific indels to the exclusion of the other amitochondriates .
however, none of those indels is exclusive to the diplomonad and trimastix sequences but is also found in related fbas of bacteria .
our tree shows some similarity to an earlier study, where mastigamoeba and entamoeba also group together  <cit> ; in that study trichomonas groups with these two amitochondriates and this monophyly was not rejected by a shimodaira-hasegawa test   <cit> .
in our study however, trichomonas groups with treponema, albeit bootstrap support is not very high.
lateral transfer scenarios that account for this particular phylogenetic distribution in eukaryotes are complex, especially when relationships amongst the excavate taxa are concerned.
regardless, our phylogeny and indel data suggest that amitochondrial protists have acquired their fbas from at least two different sources implying at least two individual lgt events.
triosephosphate isomerase
the tpis of all eukaryotes including amitochondriates are monophyletic .
in our analysis, the monophyly of major eukaryotic phyla  like animals, fungi, green plants, kinetoplastids, stramenopiles, and diplomonads are recovered although monophyly of fungi is only poorly supported.
the relationships between the eukaryotic lineages are not resolved.
the amitochondriates trimastix pyriformis, trichomonas vaginalis, the diplomonads, and entamoeba histolytica each branch individually in different parts of the tree.
since there is no clear affiliation between different phyla nor is one nested in another, the tpi tree does not support, nor conflict with, the commonly accepted higher-order eukaryotic relationships  <cit> .
therefore, it seems likely that the tpi genes in all eukaryotes studied here including trimastix were inherited vertically from the last common ancestor of eukaryotes.
glyceraldehyde-3-phosphate dehydrogenase
the phylogenetic relationships of gapdhs have been studied in detail in a wide variety of organisms and the evolutionary history of this enzyme has turned out to be very complex  <cit> .
for simplicity, the gapdhs may be divided into two general groups, the eukaryote cytosolic gapc group and the bacterial gapab group .
green plants contain in addition to their cytosolic gapc homolog, gapa and gapb duplicates which they acquired from a single ancestral gene in the ancestral chloroplast endosymbiont and which is used in the calvin cycle  <cit> .
in diplonemids  an obscure additional gapab gene of a different origin was detected  <cit> .
in parabasalids, no gapc has ever been found.
they instead possess bacterial-like gapab homologs which are located in the cytosol and are active in glycolysis  <cit> .
curiously, in our phylogenetic analyses, the gapdh sequence of trimastix pyriformis is the sister of the gapdh of parabasalids  supported by high bootstrap values.
the tight parabasalid/trimastix clade is loosely associated with the sequence of the spirochaete borrelia burgdorferi  similar to earlier reports of the affinities of the parabasalid homologs  <cit> .
despite the poor bootstrap support for this grouping, a common origin for the parabasalid, trimastix and spirochaete homologs is confirmed by a unique homologous insertion, which is  <dig> amino acids long in parabasalids, but shorter in trimastix and borrelia .
another characteristic found in the trimastix gapdh is a eubacterial-like sequence in the s-loop of the enzyme, which is also seen in trichomonas  <cit>  and the other parabasalids .
phosphoglycerate kinase
pgks of archaebacteria, bacteria and eukaryotes are homologous  <cit> .
with few exceptions, eukaryotic pgk genes constitute a monophyletic group nested within the bacterial clade .
the structure of the bacterial subtree is not resolved and only few bacterial groups are recovered.
in an earlier study, pgks of kinetoplastids were located basal in the eukaryotic clade to the exclusion of chloroplast and bacterial pgks  <cit> .
however, with a broader sampling of bacterial pgks, the homologs in kinetoplastids are separated from the major eukaryote clade by the sequence of aquifex  and the chloroplast targeted pgks branch with those from cyanobacteria.
the pgk of trimastix pyriformis also falls within the bacterial pgks but its sister lineage can not be determined exactly.
a relationship to clostridium perfringens  or to deinococcus radiodurans is suggested  but without any statistical support.
a possible monophyly of the trimastix and kinetoplastids pgks can not be excluded entirely.
in contrast, the pgks of the amitochondriates giardia lamblia and trichomonas vaginalis are part of the eukaryotic clade.
pgks of the eukaryotic type are supposed to have a typical surface loop, which is shortened in bacterial pgks  <cit> .
in accordance with its phylogenetic position, the loop is missing completely in the pgk of trimastix .
it is present but shorter in the pgks of kinetoplastids and the bacterium staphylococcus aureus.
on the other hand, the loop is also shorter in the pgks of the protists naegleria gruberi and giardia lambia which branch in the eukaryotic clade.
a one amino acid insertion is typical for pgks of eukaryotes and not found in pgk of bacteria, kinetoplastids, chloroplasts and trimastix .
in summary, the pgk of trimastix appears to be of different origin than pgks of most eukaryotes and likely derives from an lgt event, although a bacterial donor lineage can not be determined.
phosphoglycerate mutase
two types of pgam can be distinguished that seem to be evolutionarily unrelated: a cofactor-dependent form  found in vertebrates and yeasts and a cofactor-independent form , mainly found in plants, algae and invertebrates  <cit> .
several eubacteria possess both forms and recently it has been shown that archaebacteria also have both forms  <cit> .
trimastix harbors a cofactor-independent pgam, which groups together with the two kinetoplastid and the plant sequences  in a cluster that also contains the amitochondriate giardia lamblia, a spirochaete , and a proteobacterium .
giardia branches together with dechloromonas, whereas leptospira branches at the base of the trimastix/plant/kinetoplastid clade.
although bootstrap support is very high for the whole cluster, it is not clear if both bacterial sequences are involved in lgt events or only one of them.
the main cluster of the tree in fig.
<dig> comprises prokaryotes and several eukaryotic sequences .
within that cluster, the two red algal chloroplast sequences branch with high bootstrap support with cyanobacteria as expected given the endosymbiotic origin of chloroplasts.
the remaining eukaryotes from this part of the tree  branch in a cluster with a proteobacterium and cytophaga, but bootstrap support is very low.
our tree suggests that ipgam was involved in numerous lateral gene transfer events within all major groups of organisms including archaebacteria, eubacteria, and the eukaryotes.
eukaryotes seem to have acquired ipgam several times independently, however weak or missing bootstrap support leaves it open which sequences were derived from which donor group.
enolase
the enolases of eukaryotes are more closely related to archaebacteria than to bacteria ,  <cit>  and it was concluded that the enolase genes were transmitted vertically from the archaebacterial ancestors to eukaryotes  <cit> .
the only exceptions are the cytosolic and probably plastid-targeted enolases of euglena gracilis, which are missing a eukaryote-specific indel, and cluster among bacterial homologs with the spirochaete treponema pallidum   <cit> .
within the eukaryotic clade, intracellular transfer of enolase genes from the eukaryotic endosymbiont to the nucleus of the host cell was detected in chlorarchniophytes and cryptophytes  <cit> .
on the subgenic level, insertions were probably transmitted between alveolates and land plants, two distantly related lineages, by lateral transfer and fine-scale recombination, resulting in a mosaic gene  <cit> .
parabasalids are the deepest branch in the eukaryotic tree , but the discovery that a two amino acid deletion which was believed to be unique for parabasalia is in fact a polymorphic character for this group does not support the hypothesis that they are the earliest diverging eukaryotes  <cit> .
the enolase of trimastix pyriformis is recovered as related to kinetoplastids with moderate bootstrap support, however the relationship of this clade to other eukaryotic taxa is not well resolved despite the fact that the eukaryote grouping itself is strongly supported .
given this phylogeny, it is simplest to assume that trimastix inherited its enolase gene by vertical descent.
pyruvate kinase
phylogenetic trees of pks divide the bacterial sequences into two main clusters .
this split topology may be explained by an ancient gene duplication  <cit>  or a complex lgt scenario.
in any case, cytoplasmic pks of fungi, kinetoplastids, animals, and plants form a tight clade, which is embedded in one of the bacterial clusters.
by contrast, the pks of the amitochondriate protists have a different origin: the pk of mastigamoeba balamuthi is in the same subtree as the cytoplasmatic pks, but it is closely related with an enzyme of borrelia burgdorferi, a relationship that is supported by 100% bootstrap values.
the pks of trimastix pyriformis, giardia lamblia, and trichomonas vaginalis fall within the second subtree, however, the branching order within this subtree is not well resolved.
the enzyme from giardia, although highly divergent in its amino acid sequence, very strongly groups with proteobacterial homologs to the exclusion of the other eukaryotes.
both the trichomonas and the trimastix sequences emerge in separate parts of this subtree, but the internal branching order depends on methods of analysis and is poorly supported in all cases.
thus although we cannot reconstruct the history of the trimastix and trichomonas genes with any precision, it does seem likely that they acquired their pk homologs in at least one event of lgt from a eubacterium, probably separately from giardia.
pyruvate-phosphate dikinase
in our phylogenetic analyses of diverse ppdks sequences, the chloroplast targeted ppdks from plants and a ppdk of a stramenopile form a tight clade to the exclusion of ppdks of bacteria, archaebacteria, kinetoplastids, and amitochondriate protists .
the ppkds of protists are not monophyletic; the two possible ppdk genes of trichomonas vaginalis are very similar and branch as sistergroup to the ppdks of the gram positive  bacteria.
ppdks of the amoeboid amitochondriates entamoeba histolytica and mastigamoeba balamuthi form a well-supported clade, which surprisingly branches as sistergroup with the ppdk of the proteobacterium methylococcus capsulatus which in turn branches with the hyperthermophilic archaeon pyrobaculum aerophilum, a grouping that was recently shown also by slamovits and keeling  <cit> .
while the relationship of the kinetoplastid ppdks within the tree is not fully resolved, they branch as the sister group to proteobacteria and chlorobium plus bacteroides, but receive no statistical support.
trimastix pyriformis and the two streblomastix sequences form a well supported clade which in turn branches together with giardia lamblia, with low to moderate support.
overall our tree does not support a common origin of all amitochondriate ppdks as the parabasalid and the amoebozoan homologs are each strongly allied with a distinct prokaryotic group.
however, the placement of the trimastix ppdk with oxymonads and the diplomonad in an unresolved portion of the tree prevents definitive conclusions to be made regarding its origins but is suggestive of an lgt event to a common ancestral excavate.
discussion
phylogenetic patterns in glycolytic enzymes of amitochondriates
the extent of lateral gene transfer among glycolytic enzymes in trimastix pyriformis is comparable with that in giardia lamblia, trichomonas vaginalis, and entamoeba histolytica.
at least four of the ten glycolytic enzymes that we identified in the ests of trimastix , have been acquired from a bacterium by lgt and the phylogenies of two more enzymes  suggest lgts as well.
although we propose that these events happened in the direction of prokaryote to eukaryote we cannot completely exclude the possibility that the opposite has happened, where a eukaryote donated a gene to a prokaryote.
for instance, for some of these enzymes  the backbone phylogeny of prokaryotes is scrambled relative to accepted taxonomic relationships implying that gene transfers must have been occurring between these lineages as well.
however, there are several reasons why we suspect most of the examples we discuss are more likely prokaryote-to-eukaryote transfers.
the existing literature on lgt indicates that well substantiated cases of eukaryote to prokaryote transfer are rare, possibly because of difficulties in prokaryotes accommodating eukaryotic gene structure   <cit> .
furthermore, in some of the cases that we discuss, the number of eukaryotes that possess homologs of an enzyme sub-type is very small relative to the number of prokaryotes and the eukaryotes are often distantly separated on the tree .
in other cases, such as gapdh , the putative eukaryotic recipients of the transfers are distantly related to major eukaryotic clades nested deep within a group of prokaryotic and organellar homologs.
for the glycolytic enzymes we have examined, the number of inter-domain lgt events required to explain the data are far fewer if eukaryotes are the recipients and prokaryotes the donor.
although the trees supported separation of the eukaryote lineages and their emergence within predominantly prokaryotic clades, we could not determine the bacterial donor lineage implicated because of the poorly resolved branching patterns in the trees.
taking also the other amitochondriates into account, there are only few examples where we can trace back the donor lineages.
spirochaetes seem to be involved quite often in lgt events  as was previously noted  <cit>  and there was a conspicuous lack of putative lgt events between archaebacteria and eukaryotes.
this is curious because lgts between archaebacteria and amitochondriates are not exceptional and have been reported for enzymes catalyzing the reactions downstream from the glycolytic pathway involved in the fermentation of pyruvate  <cit> .
another curious pattern is that two enzymes in trimastix that were clearly of vertical origin – tpi and enolase – were also apparently inherited vertically in the other amitochondriate eukaryotes .
this indicates that the selection of enzymes in glycolysis that are replaced by lgts in amitochondriate  eukaryotes is not random.
selection for particular properties of these enzymes in amitochondriate protists must play some kind of role in this process, although precisely what these properties are remains a complete mystery.
alliances between 'excavate' homologs
in our trees, glycolytic enzyme homologs from excavates tend to cluster together to some extent .
however, the relationships between excavate homologs in the different enzyme trees are mutually contradictory.
reconciling these groupings with known excavate relationships is a complex undertaking, especially given the uncertainty in the validity of an 'excavata' clade  <cit> .
nevertheless, multiple gene phylogenies support several robust groupings of excavates, of which three are relevant here:  diplomonads plus parabasalids ;  trimastix and oxymonads ; and  euglenoids and kinetoplastids   <cit> .
if we assume these groupings are correct, then we can suggest plausible gene transfer scenarios.
for instance, the type ii gpi of diplomonads and parabasalids is explained by a gene transfer event in their common 'metamonad' ancestor, while all other eukaryotes acquired their type i gpi homologs in one  separate event from a eubacterial donor .
on the other hand, the fba tree suggests a common ancestor of diplomonads and trimastix inherited an lgt of the class ii type b enzyme to the exclusion of parabasalids .
in this case, a common ancestor of the preaxostyla and metamonada may have acquired the fba enzyme by lgt from eubacteria, with a second separate event replacing the enzyme later in parabasalids.
in this scenario the kinetoplastids may retain the original 'excavate' class i fba , while a separate lgt must be postulated to explain why euglena and yeasts possess the type a version of the class ii enzyme  <cit> .
the gapdh tree indicates that trimastix and parabasalids share a 'eubacterial' type enzyme whereas diplomonads and euglenozoa have the canonical 'eukaryotic' versions of the enzyme .
in this case it seems likely that either trimastix or an ancestral parabasalid acquired the enzyme from a spirochaete-like donor and 'passed' this version to the other eukaryote by lgt.
evidence for such eukaryote-eukaryote lgt events is beginning to accumulate  <cit>  and makes ecological sense if ancestors of these lineages coexisted in anoxic environments.
finally, in the ppdk tree trimastix groups with the oxymonad streblomastix and the diplomonad giardia, while the parabasalids branch independently among the high-gc gram positive bacteria .
as both of these clades emerge from within the bacteria it seems possible that either an ancestral 'excavate' lineage acquired the ppdk from a eubacterial donor with the parabasalids replacing their version later by an additional lgt event from eubacteria or trimastix or the diplomonads acquired ppdk from the other by lgt.
in either scenario, two lgt events are required.
the foregoing interpretations, while reasonable, are ad hoc explanations of observed branching patterns and insertion/deletion characters.
alternative lgt scenarios that account for the data may be equally plausible.
however, to test these scenarios more information of two sorts is sorely needed.
first, full genome sequences from the taxa represented in the tree are needed to determine if the homologs we have analyzed are the only members of the enzyme family present in the genome; if multiple copies exist in these genomes, then the phylogenetic histories of the copies may provide clues as to which ones are 'ancestral' and which are newly acquired by lgt.
second, a much better sampling of protistan and eubacterial genomes related to those implicated in the gene transfer events could provide key information regarding when such gene transfers took place, and allow the direction of transfer to be discerned.
indeed, a recent investigation aimed at broadening taxonomic sampling was able to provide much more precise information regarding the timing of lgt events in eukaryotic genomes  <cit> .
can phylogenetic error be the explanation for the aberrant branching patterns?
as in every phylogenetic analysis, we may encounter random and systematic error.
the lack of resolution in our trees can not be avoided due to the limited amount of data.
bootstrap values allow us at least to judge the statistical significance of bipartitions and determine which nodes are unresolved.
we did not compare different tree topologies by likelihood ratio tests since this is too time consuming to complete for multiple topologies and the  <dig> enzyme families we examined.
it was more important for us to determine if a gene from an amitochondriate branches within the main eukaryotic cluster than to trace back the actual bacterial donor lineage.
for rather ancient lgts, identifying the donor lineage might be impossible given the high background of bacterial lgts, as pointed out by others  <cit> .
we included some sequences in our analyses which did not pass the chi-squared tests for homogeneity of amino acid frequencies as implemented in tree-puzzle  and are therefore prone to phylogenetic artefacts.
in most cases, these sequences branched with closely-related organisms in regions of the tree that were irrelevant to our conclusions.
one exception is the pk of giardia lamblia that is indeed represented by a very long branch and has significantly deviant amino acid composition .
however, it clusters with the short branches of proteobacterial pks  that are themselves not compositionally biased, suggesting that this grouping is likely historical rather than artefactual.
dictyostelium discoideum
gracilaria tenuistipitata
porphyra purpurea
encephalitozoon cuniculi
<dig> the gamma shape parameter of the gamma distribution estimated using tree-puzzle
<dig> based on the chi square tests for deviation of amino acid frequencies implemented in tree-puzzle
not only is there conspicuous lack of evidence for 'phylogenetic artefacts' accounting for the aberrant branching patterns we describe, in a number of cases  we have also found characteristic insertion and deletion patterns in the enzyme homologs that support aspects of the phylogeny and the inference of lgt.
as these regions were always removed from the data prior to the phylogenetic analyses and are not subject to phylogenetic 'artefacts' of the same sort as the trees, we suggest that the strongly supported phylogenetic patterns that we report likely represent the true history of the protein families.
what about ancient paralogy and differential loss?
provided that the phylogenies were correctly inferred, therefore, the unexpected relationships can be explained generally either by proposing lgt events or by ancient gene duplications followed by selective gene loss events in independent lineages  <cit> .
however, the hallmarks of such 'paralogy scenarios' are multiple copies of enzymes in many taxa that generate 'mirror' organismal trees with random taxonomic gaps reflecting differential gene loss.
our phylogenies do not show this pattern.
thus if ancient paralogy were invoked to explain our data, it would require positing hypothetical ancestral organisms that retained progressively more and more paralogs of the enzymes in question as one moves deeper in the tree of life.
furthermore, cataclysmic numbers of differential losses of these paralogs would have to occur on independent lineages to generate the extant pattern.
this scenario seems unlikely at best.
however, although we do not see any evidence for paralogy, it is possible that some mixture of paralogy/gene loss and lgt events have generated the observed distribution of enzyme types in eukaryotes.
once again, full genome sequences from a much greater diversity of organisms are required to evaluate the relative likelihood of these alternative scenarios.
CONCLUSIONS
replacement of glycolytic enzymes by lgt: selection or neutral evolution?
do trimastix pyriformis and other amitochondriates gain any advantages by replacing the canonical eukaryotic versions of the glycolytic enzymes or is it just a random process  <cit> ?
the use of pyrophosphate-linked instead of atp-dependent enzymes can increase the overall efficiency of glycolysis, which is especially important if atp production is solely carried out by glycolysis  <cit> .
pyrophosphate  is created during biosynthetic polymerization reactions and, in most organisms is hydrolyzed by inorganic pyrophosphatase in order to thermodynamically favor the anabolic reactions.
organisms like giardia lamblia, which presumably lack inorganic pyrophosphatase, can use pyrophosphate instead of atp as the phosphor-donor in some reactions  <cit> .
in trimastix we found in addition to an atp-dependent pk, the pyrophosphate-linked ppdk as a second enzyme that catalyzes the conversion of phosphoenolpyruvate  to pyruvate .
for giardia, in vitro tests suggested that both pk and ppdk play a role in glycolysis  <cit> .
for trichomonas, on the other hand, only pk activity could be detected  <cit>  but the preliminary tigr database of the trichomonas genome contains two  ppdk genes .
like in trimastix, the mere presence of these two genes allows no conclusions to be made regarding their metabolic function.
however, one could speculate that only one of the two is active in glycolysis, whereas the other one is still present in the genome but on its way to losing its function.
it is noteworthy that we found over  <dig> times more ests for ppdk than for pk in the trimastix library suggesting that the expression of ppdk is significantly higher in trimastix than that of pk.
the closely related oxymonad streblomastix strix has also been shown to harbour ppdk , but an atp-dependent pk has not yet been found  <cit> .
ppdk is also found in kinetoplastids, where it is involved in pyrophosphate recycling  <cit> , and in plants, where it is involved in fixation of co <dig>  <cit> .
our tree suggests that the ppdks of all amitochondriates sampled here have been acquired by lateral gene transfer events.
it is tempting to assume that the acquisition of this ppi-linked enzyme is the result of an adaptation to the amitochondriate nature of these organisms, thus enabling them to increase atp yield during glycolysis.
interestingly, laterally acquired glycolytic enzymes of amitochondriates were shown to have biochemical properties more similar to the bacterial homologs than to the canonical eukaryotic enzymes consistent with their phylogenetic position  <cit> .
in trichomonas pk is stimulated by ribose 5-phosphate and glycerate 3-phosphate, as is the case for many bacteria, while most eukaryotic pks are allosterically activated by fructose  <dig> -bisphosphate  <cit> .
similarly, the trichomonas gapdh exhibits sensitivity to the inhibitor koningic acid at comparable levels to eubacterial homologs and two orders of magnitude less than eukaryotic gapdhs.
this sensitivity is directly linked to the sequence of the s-loop domain of the enzyme which has a eubacterial signature in trichomonas  <cit>  and trimastix.
a selective advantage may be hidden behind these differences but it is not yet recognized.
the same is true for fbas of amitochondriates, which are all of class ii and have been acquired in several independent lgt events.
thus, it seems that laterally transferred genes often retain some of their 'bacterial properties', although the biochemical properties of these enzymes need to be studied in much more detail.
the most abundant pattern of lgt in eukaryotes is the transfer of genes from the bacterial ancestors of the modern mitochondrion and plastid to the cell nucleus during the process of organelle genome reduction  <cit> .
several proteins involved in the energy metabolism of eukaryotes are sister to or nested within alpha-proteobacteria in phylogenetic inferences which is seen as straightforward evidence for mitochondrion-to-host gene transfer  <cit> .
some go further and regard every affiliation of a eukaryotic protein to any bacterial phylum as probable endosymbiotic gene transfer  <cit> : first, given the extent of lgt and gene loss among bacteria, the respective homolog may no longer be found in extant alpha proteobacteria; or second, during the course of the gene transfer from endosymbiont to the host genome, the gene had no functional constraints and may have acquired numerous mutations which now mask its origin.
plants recruited a huge number of their nuclear genes from the plastid genome as revealed by a comparison of the genomes of arabidopsis thaliana, cyanobacteria, and chloroplasts  <cit> .
massive endosymbiotic gene transfers from eukaryote to eukaryote were observed from the nucleomorph  to the nucleus of the mixotrophic alga bigelowiella natans  <cit> .
the extent of lgt from non-organellar  donors to eukaryotic genomes seems to vary from lineage to lineage.
for rumen ciliates, diplomonads, bigelowiella natans, and trimastix pyriformis , it was shown that they acquired a significant number of genes by lgt from different bacterial, and even eukaryotic, lineages  <cit> .
more sporadic, but steadily increasing in number, are documented lgts to other protist lineages  <cit> .
the pattern that seems to emerge is that lgt is not infrequent in phagotrophic protists, while its impact on the evolution of other eukaryotic lineages is either minor or simply not yet known  <cit> .
a mechanism has been proposed, where small pieces of dna of the engulfed bacteria escape digestion and are incorporated into the protist genome replacing the ancestral eukaryotic genes over time in a rather random manner  <cit> .
as already discussed, the high frequency of lgt in amitochondriates may be partially but not entirely caused by the fact that these protists live phagotrophically.
selection as a result of adaptation to certain environments might favor the uptake of specific proteins more suitable to the host organism, making lgt a non-random process  <cit> .
yet, the relative importance of random neutral evolution versus the effect of selection in causing lgt in eukaryotes is something we can not answer until we understand more about the general frequency of lgt in eukaryotic microbes and the functional properties of the enzymes involved.

 eukaryotes catabolize glucose to pyruvate via the embden-meyerhof-parnas glycolytic pathway. a number of protists are known that lack classical mitochondrial structures and electron-transport-linked atp production, instead relying on substrate-level phosphorylation. in these organisms, glycolysis and subsequent fermentation reactions play a central role in their energy metabolism. although often referred to as'amitochondriates', most of these protist taxa have been shown to harbor organelles of mitochondrial descent. hydrogenosomes for example metabolize pyruvate via a fermentative pathway that yields atp and hydrogen. other mitochondrion-derived organelles, such as the'mitosomes' of the parasites entamoeba histolytica and giardia lamblia, seem to have lost all core metabolic functions. in these organisms, pyruvate is catabolized by fermentation in the cytosol. although the core emp pathway appears to be almost universal amongst eukaryotes, amitochondriate protists have been shown to be surprisingly flexible in terms of the specific enzyme families they utilize to carry out glycolytic reactions. more interestingly these enzymes appear to have rather diverse origins and it is likely that this variability was achieved, at least in part, by lateral gene transfer. for instance, in the first step of glycolysis, diplomonads and parabasalids use glucokinase, rather than hexokinase which carries out this reaction in most other eukaryotes. similarly, amitochondriate protists such as parabasalids, diplomonads, pelobionts and entamoebids cleave fructose- -bisphosphate with a class ii fructose-bisphosphate aldolase, while most other eukaryotes possess a class i fba. these two enzymes are non-homologous and belong to different superfamilies. in many cases, the glycolytic enzymes that are replaced by lgt appear to come from a bacterial donor lineage. for instance, diplomonad and parabasalid glucose phosphate isomerases are more closely related to the gpis of cyanobacteria and chloroplasts than to cytosolic gpis of other eukaryotes. the glyceraldehyde-3-phosphate dehydrogenase of parabasalids is not part of the gapc clade which encompasses most eukaryotic gapdhs including those of most amitochondriates. instead, it clusters most closely with the sequences of the spirochaete genus borrelia within the bacterial gapab clade. similarly, the evolution of phosphofructokinase shows phospho-donor changes and frequent lateral gene transfer. in amitochondriates, pyrophosphate-dependent, atp-dependent or both types of pfk are present – most of them acquired by lgt. in some protist taxa, additional enzymes involved in glycolytic reactions are found that are rare amongst most eukaryotes. for instance, pyruvate phosphate dikinase is an analogous enzyme to pyruvate kinase, with a key difference being that ppdk uses pyrophosphate instead of atp to catalyze the conversion of phosphoenolpyruvate into pyruvate and generates atp in the process. ppdk is present in a number of eukaryotes including the anaerobes entamoeba histolytica, giardia lamblia, trichomonas vaginalis and streblomastix strix and the phylogenetic affinites of the eukaryote enzymes are diverse. these examples demonstrate the diversity in the equipment of glycolytic enzymes and in their evolutionary history in protists in general, and anaerobic species in particular. however, our knowledge of the biochemistry of diverse protist lineages still remains severely limited to taxa of biomedical or agricultural importance. here we expand this knowledge to include the anaerobic flagellate trimastix pyriformis. species of the genus trimastix are free-living nanoflagellates which thrive in anoxic environments, where they feed on bacteria. they harbor small double membrane-bound organelles of unknown function which structurally resemble hydrogenosomes but are probably derived from mitochondria. phylogenetic analyses indicate that trimastix species form a sister group to the oxymonads, another goup of'amitochondriate' protist that are symbionts of termites and cockroaches. the placement of these organisms in the eukaryotic tree remains contentious, but on the basis of ultrastructural comparisons they have been suggested to be part of a eukaryote super-group called the excavata. the excavata is comprised of a number of aerobic protist lineages including the heterolobosea, euglenozoa and jakobid flagellates as well as anaerobic lineages such as diplomonads, carpediemonas and parabasalids. although, ultrastructural evidence for the excavata is convincing, molecular phylogenetic support remains weak. to improve our knowledge of the biochemical diversity of anaerobic protists we initiated an est project of trimastix pyriformis. we identified a large number of ests coding for enzymes involved in glycolysis and subsequent fermentation. our phylogenetic analyses reveal that, while some glycolytic enzymes in this organism were inherited by vertical descent, about half of the glycolytic enzymes of trimastix were acquired by lateral gene transfer events, likely from different bacterial donors. although the glycolytic pathway is universally conserved in eukaryotes, its enzyme components are evolutionarily labile and have been repeatedly replaced in separate eukaryotic lineages by gene transfer from diverse eubacterial donors.  phylogenies of glycolytic enzymes of trimastix pyriformis from our est data we identified homologs of glycolytic enzymes typically found in other eukaryotes plus a second enzyme for the final glycolytic step ). in order to completely sequence each homolog, 5'-race and primer walking sequencing were performed. putative enzyme activities, lengths, and gc contents of the encoded proteins are shown in table in the following sections we describe the properties and phylogenetic affinities of each enzyme in turn. hexokinase the evolution of hexokinases is characterized by several independent gene duplications which gave rise to numerous isoenzymes within different phyla: plants have two isoenzymes, vertebrates have four and yeasts have two; the vertebrate hexokinase d is often referred to as glucokinase. like most aerobic eukaryotes and the amitochondriate parasite entamoeba histolytica, trimastix appears to use hexokinase for the first enzymatic step in glycolysis whereas the amitochondriate'excavates' giardia intestinalis, spironucleus barkhanus and trichomonas vaginalis all appear to use glucokinase. it is uncertain if bacteria use a real hexokinase, however search of several databases revealed three prokaryotes with hexokinases that were readily aligned with the eukaryotic hexokinase genes. these three prokaryotic sequences form a distinct clade that is clearly separated from the eukaryotic hexokinase sequences. the eukaryotic part of the tree shows very little resolution between major groups with fungi appearing to be paraphyletic with respect to bacteria and other eukaryotes, a feature that is not supported by bootstrap analysis. the trimastix homolog emerges firmly within the eukaryotes, but is not strongly allied with any other lineages. thus the trimastix hexokinase gene is most probably of eukaryotic origin via vertical descent. glucose phosphate isomerase the gpi gene family can be divided into three major types of gpis which show only a low degree of sequence similarity: type i gpis are found in the eukaryotic cytosol and in many bacteria. type ii are typical for cyanobacteria and chloroplasts but were recently also detected in the amitochondriate excavates giardia intestinalis, spironucleus barkhanus, and trichomonas vaginalis. finally, the gpis of the archaeon methanococcus jannaschii and some gram positive bacteria are assigned to type iii. for tree reconstruction, we used only gpis of type i and ii. the gpi of trimastix pyriformis is of type i and not similar to the type ii gpis of diplomonads and parabasalids. in the type i subtree, several bacterial lineages and two clades of eukaryotic gpis can be distinguished: one of the clades consists of cytosolic gpis from alveolates, stramenopiles, and green plants. the second comprises cytosolic gpis of fungi and animals, glycosomal gpis of kinetoplastids, and a clade of bacterial gpis. likelihood comparisons of different possible trees neither supported nor rejected this seemingly robust grouping of animals and bacteria and it was concluded that the branching pattern within the type i subtree is very unstable and currently unresolved. except for the most basal branches, this also describes our results. the trimastix gpi falls in the unresolved part of the type i subtree where it groups with the alpha-proteobacterium agrobacterium tumefaciens but with no statistical support. given the poor resolution of the type i gpi subtree, we cannot trace the origin of the trimastix gpi other than it is quite distinct from the type ii forms of the other amitochondriate flagellates. fructose- -bisphosphate aldolase a phylogeny of type b class ii fba enzymes is shown in fig 1a. trimastix pyriformis possesses a type b class ii fba like the other amitochondriates and branches robustly with diplomonad fba sequences in our phylogenetic analysis. the other amitochondriates form a clearly separated clade branching within diverse bacteria. an alignment of class ii type b fba protein sequences supports the phylogenetic separation of these groups of amitochondriate taxa; the fbas of the diplomonads and trimastix share some specific indels to the exclusion of the other amitochondriates. however, none of those indels is exclusive to the diplomonad and trimastix sequences but is also found in related fbas of bacteria. our tree shows some similarity to an earlier study, where mastigamoeba and entamoeba also group together; in that study trichomonas groups with these two amitochondriates and this monophyly was not rejected by a shimodaira-hasegawa test. in our study however, trichomonas groups with treponema, albeit bootstrap support is not very high. lateral transfer scenarios that account for this particular phylogenetic distribution in eukaryotes are complex, especially when relationships amongst the excavate taxa are concerned. regardless, our phylogeny and indel data suggest that amitochondrial protists have acquired their fbas from at least two different sources implying at least two individual lgt events. triosephosphate isomerase the tpis of all eukaryotes including amitochondriates are monophyletic. in our analysis, the monophyly of major eukaryotic phyla like animals, fungi, green plants, kinetoplastids, stramenopiles, and diplomonads are recovered although monophyly of fungi is only poorly supported. the relationships between the eukaryotic lineages are not resolved. the amitochondriates trimastix pyriformis, trichomonas vaginalis, the diplomonads, and entamoeba histolytica each branch individually in different parts of the tree. since there is no clear affiliation between different phyla nor is one nested in another, the tpi tree does not support, nor conflict with, the commonly accepted higher-order eukaryotic relationships. therefore, it seems likely that the tpi genes in all eukaryotes studied here including trimastix were inherited vertically from the last common ancestor of eukaryotes. glyceraldehyde-3-phosphate dehydrogenase the phylogenetic relationships of gapdhs have been studied in detail in a wide variety of organisms and the evolutionary history of this enzyme has turned out to be very complex. for simplicity, the gapdhs may be divided into two general groups, the eukaryote cytosolic gapc group and the bacterial gapab group. green plants contain in addition to their cytosolic gapc homolog, gapa and gapb duplicates which they acquired from a single ancestral gene in the ancestral chloroplast endosymbiont and which is used in the calvin cycle. in diplonemids an obscure additional gapab gene of a different origin was detected. in parabasalids, no gapc has ever been found. they instead possess bacterial-like gapab homologs which are located in the cytosol and are active in glycolysis. curiously, in our phylogenetic analyses, the gapdh sequence of trimastix pyriformis is the sister of the gapdh of parabasalids supported by high bootstrap values. the tight parabasalid/trimastix clade is loosely associated with the sequence of the spirochaete borrelia burgdorferi similar to earlier reports of the affinities of the parabasalid homologs. despite the poor bootstrap support for this grouping, a common origin for the parabasalid, trimastix and spirochaete homologs is confirmed by a unique homologous insertion, which is amino acids long in parabasalids, but shorter in trimastix and borrelia. another characteristic found in the trimastix gapdh is a eubacterial-like sequence in the s-loop of the enzyme, which is also seen in trichomonas and the other parabasalids. phosphoglycerate kinase pgks of archaebacteria, bacteria and eukaryotes are homologous. with few exceptions, eukaryotic pgk genes constitute a monophyletic group nested within the bacterial clade. the structure of the bacterial subtree is not resolved and only few bacterial groups are recovered. in an earlier study, pgks of kinetoplastids were located basal in the eukaryotic clade to the exclusion of chloroplast and bacterial pgks. however, with a broader sampling of bacterial pgks, the homologs in kinetoplastids are separated from the major eukaryote clade by the sequence of aquifex and the chloroplast targeted pgks branch with those from cyanobacteria. the pgk of trimastix pyriformis also falls within the bacterial pgks but its sister lineage can not be determined exactly. a relationship to clostridium perfringens or to deinococcus radiodurans is suggested but without any statistical support. a possible monophyly of the trimastix and kinetoplastids pgks can not be excluded entirely. in contrast, the pgks of the amitochondriates giardia lamblia and trichomonas vaginalis are part of the eukaryotic clade. pgks of the eukaryotic type are supposed to have a typical surface loop, which is shortened in bacterial pgks. in accordance with its phylogenetic position, the loop is missing completely in the pgk of trimastix. it is present but shorter in the pgks of kinetoplastids and the bacterium staphylococcus aureus. on the other hand, the loop is also shorter in the pgks of the protists naegleria gruberi and giardia lambia which branch in the eukaryotic clade. a one amino acid insertion is typical for pgks of eukaryotes and not found in pgk of bacteria, kinetoplastids, chloroplasts and trimastix. in summary, the pgk of trimastix appears to be of different origin than pgks of most eukaryotes and likely derives from an lgt event, although a bacterial donor lineage can not be determined. phosphoglycerate mutase two types of pgam can be distinguished that seem to be evolutionarily unrelated: a cofactor-dependent form found in vertebrates and yeasts and a cofactor-independent form, mainly found in plants, algae and invertebrates. several eubacteria possess both forms and recently it has been shown that archaebacteria also have both forms. trimastix harbors a cofactor-independent pgam, which groups together with the two kinetoplastid and the plant sequences in a cluster that also contains the amitochondriate giardia lamblia, a spirochaete, and a proteobacterium. giardia branches together with dechloromonas, whereas leptospira branches at the base of the trimastix/plant/kinetoplastid clade. although bootstrap support is very high for the whole cluster, it is not clear if both bacterial sequences are involved in lgt events or only one of them. the main cluster of the tree in fig. comprises prokaryotes and several eukaryotic sequences. within that cluster, the two red algal chloroplast sequences branch with high bootstrap support with cyanobacteria as expected given the endosymbiotic origin of chloroplasts. the remaining eukaryotes from this part of the tree branch in a cluster with a proteobacterium and cytophaga, but bootstrap support is very low. our tree suggests that ipgam was involved in numerous lateral gene transfer events within all major groups of organisms including archaebacteria, eubacteria, and the eukaryotes. eukaryotes seem to have acquired ipgam several times independently, however weak or missing bootstrap support leaves it open which sequences were derived from which donor group. enolase the enolases of eukaryotes are more closely related to archaebacteria than to bacteria, and it was concluded that the enolase genes were transmitted vertically from the archaebacterial ancestors to eukaryotes. the only exceptions are the cytosolic and probably plastid-targeted enolases of euglena gracilis, which are missing a eukaryote-specific indel, and cluster among bacterial homologs with the spirochaete treponema pallidum. within the eukaryotic clade, intracellular transfer of enolase genes from the eukaryotic endosymbiont to the nucleus of the host cell was detected in chlorarchniophytes and cryptophytes. on the subgenic level, insertions were probably transmitted between alveolates and land plants, two distantly related lineages, by lateral transfer and fine-scale recombination, resulting in a mosaic gene. parabasalids are the deepest branch in the eukaryotic tree, but the discovery that a two amino acid deletion which was believed to be unique for parabasalia is in fact a polymorphic character for this group does not support the hypothesis that they are the earliest diverging eukaryotes. the enolase of trimastix pyriformis is recovered as related to kinetoplastids with moderate bootstrap support, however the relationship of this clade to other eukaryotic taxa is not well resolved despite the fact that the eukaryote grouping itself is strongly supported. given this phylogeny, it is simplest to assume that trimastix inherited its enolase gene by vertical descent. pyruvate kinase phylogenetic trees of pks divide the bacterial sequences into two main clusters. this split topology may be explained by an ancient gene duplication or a complex lgt scenario. in any case, cytoplasmic pks of fungi, kinetoplastids, animals, and plants form a tight clade, which is embedded in one of the bacterial clusters. by contrast, the pks of the amitochondriate protists have a different origin: the pk of mastigamoeba balamuthi is in the same subtree as the cytoplasmatic pks, but it is closely related with an enzyme of borrelia burgdorferi, a relationship that is supported by 100% bootstrap values. the pks of trimastix pyriformis, giardia lamblia, and trichomonas vaginalis fall within the second subtree, however, the branching order within this subtree is not well resolved. the enzyme from giardia, although highly divergent in its amino acid sequence, very strongly groups with proteobacterial homologs to the exclusion of the other eukaryotes. both the trichomonas and the trimastix sequences emerge in separate parts of this subtree, but the internal branching order depends on methods of analysis and is poorly supported in all cases. thus although we cannot reconstruct the history of the trimastix and trichomonas genes with any precision, it does seem likely that they acquired their pk homologs in at least one event of lgt from a eubacterium, probably separately from giardia. pyruvate-phosphate dikinase in our phylogenetic analyses of diverse ppdks sequences, the chloroplast targeted ppdks from plants and a ppdk of a stramenopile form a tight clade to the exclusion of ppdks of bacteria, archaebacteria, kinetoplastids, and amitochondriate protists. the ppkds of protists are not monophyletic; the two possible ppdk genes of trichomonas vaginalis are very similar and branch as sistergroup to the ppdks of the gram positive bacteria. ppdks of the amoeboid amitochondriates entamoeba histolytica and mastigamoeba balamuthi form a well-supported clade, which surprisingly branches as sistergroup with the ppdk of the proteobacterium methylococcus capsulatus which in turn branches with the hyperthermophilic archaeon pyrobaculum aerophilum, a grouping that was recently shown also by slamovits and keeling. while the relationship of the kinetoplastid ppdks within the tree is not fully resolved, they branch as the sister group to proteobacteria and chlorobium plus bacteroides, but receive no statistical support. trimastix pyriformis and the two streblomastix sequences form a well supported clade which in turn branches together with giardia lamblia, with low to moderate support. overall our tree does not support a common origin of all amitochondriate ppdks as the parabasalid and the amoebozoan homologs are each strongly allied with a distinct prokaryotic group. however, the placement of the trimastix ppdk with oxymonads and the diplomonad in an unresolved portion of the tree prevents definitive conclusions to be made regarding its origins but is suggestive of an lgt event to a common ancestral excavate. discussion phylogenetic patterns in glycolytic enzymes of amitochondriates the extent of lateral gene transfer among glycolytic enzymes in trimastix pyriformis is comparable with that in giardia lamblia, trichomonas vaginalis, and entamoeba histolytica. at least four of the ten glycolytic enzymes that we identified in the ests of trimastix, have been acquired from a bacterium by lgt and the phylogenies of two more enzymes suggest lgts as well. although we propose that these events happened in the direction of prokaryote to eukaryote we cannot completely exclude the possibility that the opposite has happened, where a eukaryote donated a gene to a prokaryote. for instance, for some of these enzymes the backbone phylogeny of prokaryotes is scrambled relative to accepted taxonomic relationships implying that gene transfers must have been occurring between these lineages as well. however, there are several reasons why we suspect most of the examples we discuss are more likely prokaryote-to-eukaryote transfers. the existing literature on lgt indicates that well substantiated cases of eukaryote to prokaryote transfer are rare, possibly because of difficulties in prokaryotes accommodating eukaryotic gene structure. furthermore, in some of the cases that we discuss, the number of eukaryotes that possess homologs of an enzyme sub-type is very small relative to the number of prokaryotes and the eukaryotes are often distantly separated on the tree. in other cases, such as gapdh, the putative eukaryotic recipients of the transfers are distantly related to major eukaryotic clades nested deep within a group of prokaryotic and organellar homologs. for the glycolytic enzymes we have examined, the number of inter-domain lgt events required to explain the data are far fewer if eukaryotes are the recipients and prokaryotes the donor. although the trees supported separation of the eukaryote lineages and their emergence within predominantly prokaryotic clades, we could not determine the bacterial donor lineage implicated because of the poorly resolved branching patterns in the trees. taking also the other amitochondriates into account, there are only few examples where we can trace back the donor lineages. spirochaetes seem to be involved quite often in lgt events as was previously noted and there was a conspicuous lack of putative lgt events between archaebacteria and eukaryotes. this is curious because lgts between archaebacteria and amitochondriates are not exceptional and have been reported for enzymes catalyzing the reactions downstream from the glycolytic pathway involved in the fermentation of pyruvate. another curious pattern is that two enzymes in trimastix that were clearly of vertical origin – tpi and enolase – were also apparently inherited vertically in the other amitochondriate eukaryotes. this indicates that the selection of enzymes in glycolysis that are replaced by lgts in amitochondriate eukaryotes is not random. selection for particular properties of these enzymes in amitochondriate protists must play some kind of role in this process, although precisely what these properties are remains a complete mystery. alliances between'excavate' homologs in our trees, glycolytic enzyme homologs from excavates tend to cluster together to some extent. however, the relationships between excavate homologs in the different enzyme trees are mutually contradictory. reconciling these groupings with known excavate relationships is a complex undertaking, especially given the uncertainty in the validity of an'excavata' clade. nevertheless, multiple gene phylogenies support several robust groupings of excavates, of which three are relevant here: diplomonads plus parabasalids; trimastix and oxymonads; and euglenoids and kinetoplastids. if we assume these groupings are correct, then we can suggest plausible gene transfer scenarios. for instance, the type ii gpi of diplomonads and parabasalids is explained by a gene transfer event in their common'metamonad' ancestor, while all other eukaryotes acquired their type i gpi homologs in one separate event from a eubacterial donor. on the other hand, the fba tree suggests a common ancestor of diplomonads and trimastix inherited an lgt of the class ii type b enzyme to the exclusion of parabasalids. in this case, a common ancestor of the preaxostyla and metamonada may have acquired the fba enzyme by lgt from eubacteria, with a second separate event replacing the enzyme later in parabasalids. in this scenario the kinetoplastids may retain the original'excavate' class i fba, while a separate lgt must be postulated to explain why euglena and yeasts possess the type a version of the class ii enzyme. the gapdh tree indicates that trimastix and parabasalids share a'eubacterial' type enzyme whereas diplomonads and euglenozoa have the canonical'eukaryotic' versions of the enzyme. in this case it seems likely that either trimastix or an ancestral parabasalid acquired the enzyme from a spirochaete-like donor and'passed' this version to the other eukaryote by lgt. evidence for such eukaryote-eukaryote lgt events is beginning to accumulate and makes ecological sense if ancestors of these lineages coexisted in anoxic environments. finally, in the ppdk tree trimastix groups with the oxymonad streblomastix and the diplomonad giardia, while the parabasalids branch independently among the high-gc gram positive bacteria. as both of these clades emerge from within the bacteria it seems possible that either an ancestral'excavate' lineage acquired the ppdk from a eubacterial donor with the parabasalids replacing their version later by an additional lgt event from eubacteria or trimastix or the diplomonads acquired ppdk from the other by lgt. in either scenario, two lgt events are required. the foregoing interpretations, while reasonable, are ad hoc explanations of observed branching patterns and insertion/deletion characters. alternative lgt scenarios that account for the data may be equally plausible. however, to test these scenarios more information of two sorts is sorely needed. first, full genome sequences from the taxa represented in the tree are needed to determine if the homologs we have analyzed are the only members of the enzyme family present in the genome; if multiple copies exist in these genomes, then the phylogenetic histories of the copies may provide clues as to which ones are'ancestral' and which are newly acquired by lgt. second, a much better sampling of protistan and eubacterial genomes related to those implicated in the gene transfer events could provide key information regarding when such gene transfers took place, and allow the direction of transfer to be discerned. indeed, a recent investigation aimed at broadening taxonomic sampling was able to provide much more precise information regarding the timing of lgt events in eukaryotic genomes. can phylogenetic error be the explanation for the aberrant branching patterns? as in every phylogenetic analysis, we may encounter random and systematic error. the lack of resolution in our trees can not be avoided due to the limited amount of data. bootstrap values allow us at least to judge the statistical significance of bipartitions and determine which nodes are unresolved. we did not compare different tree topologies by likelihood ratio tests since this is too time consuming to complete for multiple topologies and the enzyme families we examined. it was more important for us to determine if a gene from an amitochondriate branches within the main eukaryotic cluster than to trace back the actual bacterial donor lineage. for rather ancient lgts, identifying the donor lineage might be impossible given the high background of bacterial lgts, as pointed out by others. we included some sequences in our analyses which did not pass the chi-squared tests for homogeneity of amino acid frequencies as implemented in tree-puzzle and are therefore prone to phylogenetic artefacts. in most cases, these sequences branched with closely-related organisms in regions of the tree that were irrelevant to our conclusions. one exception is the pk of giardia lamblia that is indeed represented by a very long branch and has significantly deviant amino acid composition. however, it clusters with the short branches of proteobacterial pks that are themselves not compositionally biased, suggesting that this grouping is likely historical rather than artefactual. dictyostelium discoideum gracilaria tenuistipitata porphyra purpurea encephalitozoon cuniculi the gamma shape parameter of the gamma distribution estimated using tree-puzzle based on the chi square tests for deviation of amino acid frequencies implemented in tree-puzzle not only is there conspicuous lack of evidence for'phylogenetic artefacts' accounting for the aberrant branching patterns we describe, in a number of cases we have also found characteristic insertion and deletion patterns in the enzyme homologs that support aspects of the phylogeny and the inference of lgt. as these regions were always removed from the data prior to the phylogenetic analyses and are not subject to phylogenetic'artefacts' of the same sort as the trees, we suggest that the strongly supported phylogenetic patterns that we report likely represent the true history of the protein families. what about ancient paralogy and differential loss? provided that the phylogenies were correctly inferred, therefore, the unexpected relationships can be explained generally either by proposing lgt events or by ancient gene duplications followed by selective gene loss events in independent lineages. however, the hallmarks of such'paralogy scenarios' are multiple copies of enzymes in many taxa that generate'mirror' organismal trees with random taxonomic gaps reflecting differential gene loss. our phylogenies do not show this pattern. thus if ancient paralogy were invoked to explain our data, it would require positing hypothetical ancestral organisms that retained progressively more and more paralogs of the enzymes in question as one moves deeper in the tree of life. furthermore, cataclysmic numbers of differential losses of these paralogs would have to occur on independent lineages to generate the extant pattern. this scenario seems unlikely at best. however, although we do not see any evidence for paralogy, it is possible that some mixture of paralogy/gene loss and lgt events have generated the observed distribution of enzyme types in eukaryotes. once again, full genome sequences from a much greater diversity of organisms are required to evaluate the relative likelihood of these alternative scenarios.  replacement of glycolytic enzymes by lgt: selection or neutral evolution? do trimastix pyriformis and other amitochondriates gain any advantages by replacing the canonical eukaryotic versions of the glycolytic enzymes or is it just a random process? the use of pyrophosphate-linked instead of atp-dependent enzymes can increase the overall efficiency of glycolysis, which is especially important if atp production is solely carried out by glycolysis. pyrophosphate is created during biosynthetic polymerization reactions and, in most organisms is hydrolyzed by inorganic pyrophosphatase in order to thermodynamically favor the anabolic reactions. organisms like giardia lamblia, which presumably lack inorganic pyrophosphatase, can use pyrophosphate instead of atp as the phosphor-donor in some reactions. in trimastix we found in addition to an atp-dependent pk, the pyrophosphate-linked ppdk as a second enzyme that catalyzes the conversion of phosphoenolpyruvate to pyruvate. for giardia, in vitro tests suggested that both pk and ppdk play a role in glycolysis. for trichomonas, on the other hand, only pk activity could be detected but the preliminary tigr database of the trichomonas genome contains two ppdk genes. like in trimastix, the mere presence of these two genes allows no conclusions to be made regarding their metabolic function. however, one could speculate that only one of the two is active in glycolysis, whereas the other one is still present in the genome but on its way to losing its function. it is noteworthy that we found over times more ests for ppdk than for pk in the trimastix library suggesting that the expression of ppdk is significantly higher in trimastix than that of pk. the closely related oxymonad streblomastix strix has also been shown to harbour ppdk, but an atp-dependent pk has not yet been found. ppdk is also found in kinetoplastids, where it is involved in pyrophosphate recycling, and in plants, where it is involved in fixation of co. our tree suggests that the ppdks of all amitochondriates sampled here have been acquired by lateral gene transfer events. it is tempting to assume that the acquisition of this ppi-linked enzyme is the result of an adaptation to the amitochondriate nature of these organisms, thus enabling them to increase atp yield during glycolysis. interestingly, laterally acquired glycolytic enzymes of amitochondriates were shown to have biochemical properties more similar to the bacterial homologs than to the canonical eukaryotic enzymes consistent with their phylogenetic position. in trichomonas pk is stimulated by ribose 5-phosphate and glycerate 3-phosphate, as is the case for many bacteria, while most eukaryotic pks are allosterically activated by fructose -bisphosphate. similarly, the trichomonas gapdh exhibits sensitivity to the inhibitor koningic acid at comparable levels to eubacterial homologs and two orders of magnitude less than eukaryotic gapdhs. this sensitivity is directly linked to the sequence of the s-loop domain of the enzyme which has a eubacterial signature in trichomonas and trimastix. a selective advantage may be hidden behind these differences but it is not yet recognized. the same is true for fbas of amitochondriates, which are all of class ii and have been acquired in several independent lgt events. thus, it seems that laterally transferred genes often retain some of their'bacterial properties', although the biochemical properties of these enzymes need to be studied in much more detail. the most abundant pattern of lgt in eukaryotes is the transfer of genes from the bacterial ancestors of the modern mitochondrion and plastid to the cell nucleus during the process of organelle genome reduction. several proteins involved in the energy metabolism of eukaryotes are sister to or nested within alpha-proteobacteria in phylogenetic inferences which is seen as straightforward evidence for mitochondrion-to-host gene transfer. some go further and regard every affiliation of a eukaryotic protein to any bacterial phylum as probable endosymbiotic gene transfer: first, given the extent of lgt and gene loss among bacteria, the respective homolog may no longer be found in extant alpha proteobacteria; or second, during the course of the gene transfer from endosymbiont to the host genome, the gene had no functional constraints and may have acquired numerous mutations which now mask its origin. plants recruited a huge number of their nuclear genes from the plastid genome as revealed by a comparison of the genomes of arabidopsis thaliana, cyanobacteria, and chloroplasts. massive endosymbiotic gene transfers from eukaryote to eukaryote were observed from the nucleomorph to the nucleus of the mixotrophic alga bigelowiella natans. the extent of lgt from non-organellar donors to eukaryotic genomes seems to vary from lineage to lineage. for rumen ciliates, diplomonads, bigelowiella natans, and trimastix pyriformis, it was shown that they acquired a significant number of genes by lgt from different bacterial, and even eukaryotic, lineages. more sporadic, but steadily increasing in number, are documented lgts to other protist lineages. the pattern that seems to emerge is that lgt is not infrequent in phagotrophic protists, while its impact on the evolution of other eukaryotic lineages is either minor or simply not yet known. a mechanism has been proposed, where small pieces of dna of the engulfed bacteria escape digestion and are incorporated into the protist genome replacing the ancestral eukaryotic genes over time in a rather random manner. as already discussed, the high frequency of lgt in amitochondriates may be partially but not entirely caused by the fact that these protists live phagotrophically. selection as a result of adaptation to certain environments might favor the uptake of specific proteins more suitable to the host organism, making lgt a non-random process. yet, the relative importance of random neutral evolution versus the effect of selection in causing lgt in eukaryotes is something we can not answer until we understand more about the general frequency of lgt in eukaryotic microbes and the functional properties of the enzymes involved. 
the glycolytic pathway consists of ten enzymatic steps which appear to be universal amongst eukaryotes. 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. however, it has been shown that the origins of these enzymes in specific eukaryote lineages can differ, and sometimes involve lateral gene transfer events. 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. 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. we have conducted an expressed sequence tag survey of the anaerobic flagellate trimastix pyriformis to investigate the nature of the evolutionary origins of the glycolytic enzymes in this relatively unstudied organism. a number of the glycolytic enzymes of trimastix have been transferred laterally from bacteria instead of being inherited from the last common eukaryotic ancestor. glycolysis and subsequent fermentation is the main energy source for many anaerobic organisms.
 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 survey of the anaerobic flagellate trimastix pyriformis to investigate the nature of the evolutionary origins of the glycolytic enzymes in this relatively unstudied organism.  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 to pyruvate 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.  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. 