The glycolytic pathway of Trimastix pyriformis is an evolutionary mosaic
© Stechmann et al; licensee BioMed Central Ltd. 2006
Received: 09 August 2006
Accepted: 23 November 2006
Published: 23 November 2006
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.
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.
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.
Eukaryotes catabolize glucose to pyruvate via the Embden-Meyerhof-Parnas (EMP) 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 [2–4]. Although often referred to as 'amitochondriates', most of these protist taxa have been shown to harbor organelles of mitochondrial descent [5–7]. 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 [1, 8] and it is likely that this variability was achieved, at least in part, by lateral gene transfer (LGT) [9–15]. 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 [11, 13]. Similarly, amitochondriate protists such as parabasalids, diplomonads, pelobionts and entamoebids cleave fructose-1,6-bisphosphate with a class II fructose-bisphosphate aldolase (FBA), while most other eukaryotes possess a class I FBA [14, 16–18]. These two enzymes are non-homologous and belong to different superfamilies [1, 16]. 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 (GPI) are more closely related to the GPIs of cyanobacteria and chloroplasts than to cytosolic GPIs of other eukaryotes . The glyceraldehyde-3-phosphate dehydrogenase (GAPDH) 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 [10, 19, 20]. Similarly, the evolution of phosphofructokinase (PFK) 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 [9, 12, 15].
In some protist taxa, additional enzymes involved in glycolytic reactions are found that are rare amongst most eukaryotes. For instance, pyruvate phosphate dikinase (PPDK) is an analogous enzyme to pyruvate kinase (PK), 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 [21, 22], 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 [25, 26] 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 (e.g. kinetoplastids) and jakobid flagellates as well as anaerobic lineages such as diplomonads (Giardia), Carpediemonas and parabasalids (e.g. Trichomonas). Although, ultrastructural evidence for the Excavata is convincing , molecular phylogenetic support remains weak [30, 31].
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
ORFs coding for glycolytic enzymes of Trimastix pyriformis: Number of nucleotides, G+C content, and accession number
length of ORF in nucleotides
G+C content in %
glucose phosphate isomerase
fructose bisphosphate aldolase
triose phosphate isomerase
glyceraldehyde phosphate dehydrogenase
The evolution of hexokinases is characterized by several independent gene duplications which gave rise to numerous isoenzymes within different phyla: plants have two isoenzymes (hxk1, hxk2), vertebrates have four (hexokinase A-D) and yeasts have two (PI and PII) ; 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 [11, 13]. It is uncertain if bacteria use a real hexokinase, however search of several databases revealed three prokaryotes (Bacteroides and 2 spirochaetes) 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 (see Additional file 1). 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 (GPI)
The GPI gene family can be divided into three major types of GPIs which show only a low degree of sequence similarity [11, 33]: 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 (see Additional file 2). 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 (Additional file 2): 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 (Proteobacteria and Cytophaga). 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 (Additional file 2) where it groups with the alpha-proteobacterium Agrobacterium tumefaciens (ML analysis only) 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 (parabasalids and diplomonads).
Fructose-1,6-bisphosphate aldolase (FBA)
Triosephosphate isomerase (TPI)
The TPIs of all eukaryotes including amitochondriates are monophyletic (see Additional file 3). In our analysis, the monophyly of major eukaryotic phyla (with representatives from more than one taxon) 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 (GAPDH)
Phosphoglycerate kinase (PGK)
Phosphoglycerate mutase (PGAM)
The enolases of eukaryotes are more closely related to archaebacteria than to bacteria (see Additional file 4), [39–41] and it was concluded that the enolase genes were transmitted vertically from the archaebacterial ancestors to eukaryotes [39, 41]. 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 (not shown) . 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 (Additional file 4, ), 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 [40, 43]. 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 (Additional file 4). Given this phylogeny, it is simplest to assume that Trimastix inherited its enolase gene by vertical descent.
Pyruvate kinase (PK)
Pyruvate-phosphate dikinase (PPDK)
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 (see Additional file 5, [45, 46]). 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 (high GC) 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.
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 (FBA, GAPDH, PGK, and PK), have been acquired from a bacterium by LGT and the phylogenies of two more enzymes (PGAM and PPDK) 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 (e. g. FBA, Fig. 1A) 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 (e. g. introns/exons) . 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 (e. g. FBA Fig. 1A). In other cases, such as GAPDH (Fig. 2A), 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 (Figs. 1A, 2A, 5) 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 [48–51].
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 (Additional Files 3 and 4). This indicates that the selection of enzymes in glycolysis that are replaced by LGTs in amitochondriate (and other) 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 (Figs 1A, 2A, 3A, Additional Files 2 and 5). 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: (1) diplomonads plus parabasalids (and other taxa collectively known as Metamonada); (2) Trimastix and oxymonads (Preaxostyla); and (3) Euglenoids and kinetoplastids (Euglenozoa) . 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 (or possibly two) separate event from a eubacterial donor (Additional file 2). 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 (Fig. 1A). 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 (like other eukaryotes), 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 (Fig. 2A). 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 (Additional file 5). 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 10 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 [47, 52].
Overview of the datasets used in phylogenetic reconstructions
number of sequences
number of aa
sequences with divergent aa composition2
glucose phosphate isomerase
fructose bisphosphate aldolase
triose phosphate isomerase
glyceraldehyde phosphate dehydrogenase
pyruvate phosphate dikinase
Not only is there conspicuous lack of evidence for 'phylogenetic artefacts' accounting for the aberrant branching patterns we describe, in a number of cases (Figs 1, 2, 3) 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 [54–56]. Pyrophosphate (PPi) 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 [54, 55]. 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 (PEP) to pyruvate (PPDK can also work in the other direction towards gluconeogenesis). 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 (very similar) PPDK genes (Fig. 5). 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 30 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 (two copies in fact), 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 CO2 . 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 [19, 58]. 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 1,6-bisphosphate [57, 58]. 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 [11, 47, 53]. 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 [52, 62]: 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 (the reduced nucleus of the eukaryotic endosymbiont) to the nucleus of the mixotrophic alga Bigelowiella natans . The extent of LGT from non-organellar (in contrast to endosymbiotic organelle) donors to eukaryotic genomes seems to vary from lineage to lineage. For rumen ciliates, diplomonads, Bigelowiella natans, and Trimastix pyriformis (this study), it was shown that they acquired a significant number of genes by LGT from different bacterial, and even eukaryotic, lineages [64–66]. More sporadic, but steadily increasing in number, are documented LGTs to other protist lineages [11, 13, 48, 50, 51, 65, 66]. 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 [64, 66]. 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.
Source of sequences and alignments
All sequences obtained in our EST project from Trimastix pyriformis were used for searches against the non-redundant NCBI database (BLAST, ) and cDNA clones with high similarity to 10 glycolytic enzymes were found. Sequencing of all clones was completed using the primer walking method. Truncated 5' ends of cDNAs were amplified using the GeneRacer Kit (Invitrogen). Length and G+C contents of the coding regions plus accession numbers of the deduced protein sequences are summarized in table 1.
Sequences homologous to the glycolytic enzymes of Trimastix were retrieved from the NCBI non-redundant protein database. The PGK, TPI, PK, and PPDK sequences for Trichomonas vaginalis were downloaded from the TIGR database . The TPI dataset was complemented with the TPI sequences of Entamoeba histolytica, Spironucleus barkhanus, and Acrasis rosea (accession numbers EF064144–EF064146), which were obtained by PCR using degenerate primers directed against the N-terminal motifs VGGNWK (TF-1) and VGGNFK (TF-2) and the C-terminal VGGASL (TR-1). PCRs were carried out with genomic DNA using standard methods. The PGK sequence of Naegleria gruberi was obtained from another EST project (accession number EF064143).
The individual datasets were aligned with ClustalW using default settings  and subsequently adjusted manually. An alignment of enolases, to which we manually added the Trimastix sequence, was kindly provided by Patrick Keeling (University of British Columbia, Vancouver). Regions of ambiguous alignment were excluded from further analyses. Care was taken that each dataset contained a reasonable number of sequences representing all major taxonomic groups and sequence clades. The number of sequences and amino acid positions used for phylogenetic analyses of each dataset are listed in table 2.
Maximum likelihood (ML) phylogenies were inferred using PMBML , with the JTT substitution model and a mixed four category discrete gamma-model of among-site variation. The gamma shape parameter alpha was estimated using TREE-PUZZLE  (see table 2 for alpha values). TREE-PUZZLE also provided the distance matrix (using the same models) which was used to construct ML distance trees using the Fitch-Margoliash algorithm as implemented in FITCH from the Phylip package . In both types of analyses, ten random additions with global rearrangements were used to find the optimal tree. ML distance bootstrap values for bipartitions were calculated by analyses of 100 resampled data sets using PUZZLEBOOT ( distributed by A. J. Roger and M. E. Holder). The alpha parameter for each resampled dataset was calculated separately using the same models as described above. Only one random sequence addition was done. ML bootstrap values were calculated with the program phyml , using a discrete gamma model with four rate categories and 100 resamplings; the gamma parameter alpha was calculated individually for each resampled dataset. All trees were caluclated unrooted but are shown with outgroups for display purposes.
expressed sequence tags
lateral gene transfer
glucose phosphate isomerase
pyruvate phosphate dikinase
AS is supported by a postdoctoral fellowship from the Canadian Institutes of Health Research (CIHR grant MFE-76510). JDS and MB were supported by the Genome Canada/Genome Atlantic Large Scale Project: "The Protist EST Program (PEP)". The EST data from Trimastix was also generated as part of PEP. This work was supported by the operating grant MOP-62809 from the Canadian Institutes for Health Research awarded to AJR. We also thank two anonymous reviewers for their comments.
- Fothergill-Gilmore LA, Michels PA: Evolution of glycolysis. Prog Biophys Mol Biol. 1993, 59: 105-235. 10.1016/0079-6107(93)90001-Z.View ArticlePubMedGoogle Scholar
- Muller M: Energy metabolism of protozoa without mitochondria. Annu Rev Microbiol. 1988, 42: 465-488. 10.1146/annurev.mi.42.100188.002341.View ArticlePubMedGoogle Scholar
- Muller M: The hydrogenosome. J Gen Microbiol. 1993, 139: 2879-2889.View ArticlePubMedGoogle Scholar
- Muller M: Enzymes and compartmentation of core energy metabolism of anaerobic protists: a special case in eukaryotic evolution. Evolutionary relationships among protozoa. Edited by: Coombs, GH., Vickerman, K., Sleigh, MA. and Warren A. 1998, London, Chapman and Hall, 109-127.Google Scholar
- Tovar J, Fischer A, Clark CG: The mitosome, a novel organelle related to mitochondria in the amitochondrial parasite Entamoeba histolytica. Mol Microbiol. 1999, 32: 1013-1021. 10.1046/j.1365-2958.1999.01414.x.View ArticlePubMedGoogle Scholar
- Williams BA, Hirt RP, Lucocq JM, Embley TM: A mitochondrial remnant in the microsporidian Trachipleistophora hominis. Nature. 2002, 418: 865-869. 10.1038/nature00949.View ArticlePubMedGoogle Scholar
- Tovar J, Leon-Avila G, Sanchez LB, Sutak R, Tachezy J, van der Giezen M, Hernandez M, Muller M, Lucocq JM: Mitochondrial remnant organelles of Giardia function in iron-sulphur protein maturation. Nature. 2003, 426: 172-176. 10.1038/nature01945.View ArticlePubMedGoogle Scholar
- Muller M: Enzymes and compartmentation of core energy metabolism of anaerobic protists - a special case inm eukaryotic evolution?. Evolutionary relationships among protozoa. Edited by: G. H. Coombs KVMASAW. 1998, London, Chapman & HallGoogle Scholar
- Mertens E, Ladror US, Lee JA, Miretsky A, Morris A, Rozario C, Kemp RG, Muller M: The pyrophosphate-dependent phosphofructokinase of the protist, Trichomonas vaginalis, and the evolutionary relationships of protist phosphofructokinases. J Mol Evol. 1998, 47: 739-750. 10.1007/PL00006433.View ArticlePubMedGoogle Scholar
- Viscogliosi E, Muller M: Phylogenetic relationships of the glycolytic enzyme, glyceraldehyde-3-phosphate dehydrogenase, from parabasalid flagellates. J Mol Evol. 1998, 47: 190-199. 10.1007/PL00006376.View ArticlePubMedGoogle Scholar
- Henze K, Horner DS, Suguri S, Moore DV, Sanchez LB, Muller M, Embley TM: Unique phylogenetic relationships of glucokinase and glucosephosphate isomerase of the amitochondriate eukaryotes Giardia intestinalis, Spironucleus barkhanus and Trichomonas vaginalis. Gene. 2001, 281: 123-131. 10.1016/S0378-1119(01)00773-9.View ArticlePubMedGoogle Scholar
- Muller M, Lee JA, Gordon P, Gaasterland T, Sensen CW: Presence of prokaryotic and eukaryotic species in all subgroups of the PP(i)-dependent group II phosphofructokinase protein family. J Bacteriol. 2001, 183: 6714-6716. 10.1128/JB.183.22.6714-6716.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu G, Henze K, Muller M: Evolutionary relationships of the glucokinase from the amitochondriate protist, Trichomonas vaginalis. Gene. 2001, 264: 265-271. 10.1016/S0378-1119(01)00329-8.View ArticlePubMedGoogle Scholar
- Sanchez L, Horner D, Moore D, Henze K, Embley T, Muller M: Fructose-1,6-bisphosphate aldolases in amitochondriate protists constitute a single protein subfamily with eubacterial relationships. Gene. 2002, 295: 51-59. 10.1016/S0378-1119(02)00804-1.View ArticlePubMedGoogle Scholar
- Bapteste E, Moreira D, Philippe H: Rampant horizontal gene transfer and phospho-donor change in the evolution of the phosphofructokinase. Gene. 2003, 318: 185-191. 10.1016/S0378-1119(03)00797-2.View ArticlePubMedGoogle Scholar
- Marsh JJ, Lebherz HG: Fructose-bisphosphate aldolases: an evolutionary history. Trends Biochem Sci. 1992, 17: 110-113. 10.1016/0968-0004(92)90247-7.View ArticlePubMedGoogle Scholar
- Henze K, Morrison HG, Sogin ML, MullerM: Sequence and phylogenetic position of a class II aldolase gene in the amitochondriate protist, Giardia lamblia. Gene. 1998, 222: 163-168. 10.1016/S0378-1119(98)00499-5.View ArticlePubMedGoogle Scholar
- Patron NJ, Rogers MB, Keeling PJ: Gene replacement of fructose-1,6-bisphosphate aldolase supports the hypothesis of a single photosynthetic ancestor of chromalveolates. Eukaryot Cell. 2004, 3: 1169-1175. 10.1128/EC.3.5.1169-1175.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Markos A, Miretsky A, Muller M: A glyceraldehyde-3-phosphate dehydrogenase with eubacterial features in the amitochondriate eukaryote, Trichomonas vaginalis. J Mol Evol. 1993, 37: 631-643.View ArticlePubMedGoogle Scholar
- Figge RM, Cerff R: GAPDH gene diversity in spirochetes: a paradigm for genetic promiscuity. Mol Biol Evol. 2001, 18: 2240-2249.View ArticlePubMedGoogle Scholar
- Hrdy I, Mertens E, Nohynkova E: Giardia intestinalis: detection and characterization of a pyruvate phosphate dikinase. Exp Parasitol. 1993, 76: 438-441. 10.1006/expr.1993.1052.View ArticlePubMedGoogle Scholar
- Reeves RE: A new enzyme with the glycolytic function of pyruvate kinase. J Biol Chem. 1968, 243: 3202-3204.PubMedGoogle Scholar
- Slamovits CH, Keeling PJ: Pyruvate-phosphate dikinase of oxymonads and parabasalia and the evolution of pyrophosphate-dependent glycolysis in anaerobic eukaryotes. Eukaryot Cell. 2006, 5: 148-154. 10.1128/EC.5.1.148-154.2006.PubMed CentralView ArticlePubMedGoogle Scholar
- O'Kelly CJ, Farmer MA, Nerad TA: Ultrastructure of Trimastix pyriformis (Klebs) Bernard et al.: similarities of Trimastix species with retortamonad and jakobid flagellates. Protist. 1999, 150: 149-162.View ArticlePubMedGoogle Scholar
- Bernard C: Some free-living flagellates (Protista) from anoxic habitats. Ophelia. 2000, 52: 113-142.View ArticleGoogle Scholar
- Brugerolle G: Ultrastructure of Trimastix convexa Hollande, an amitochondriate anaerobic flagellate with a previously undescribed organization. European Journal of Protistology. 1997, 33: 121-130.View ArticleGoogle Scholar
- Simpson AG, Roger AJ: The real 'kingdoms' of eukaryotes. Curr Biol. 2004, 14: R693-6. 10.1016/j.cub.2004.08.038.View ArticlePubMedGoogle Scholar
- Dacks JB, Silberman JD, Simpson AG, Moriya S, Kudo T, Ohkuma M, Redfield RJ: Oxymonads are closely related to the excavate taxon Trimastix. Mol Biol Evol. 2001, 18: 1034-1044.View ArticlePubMedGoogle Scholar
- Simpson AG: Cytoskeletal organization, phylogenetic affinities and systematics in the contentious taxon Excavata (Eukaryota). Int J Syst Evol Microbiol. 2003, 53: 1759-1777. 10.1099/ijs.0.02578-0.View ArticlePubMedGoogle Scholar
- Simpson AG, Inagaki Y, Roger AJ: Comprehensive multigene phylogenies of excavate protists reveal the evolutionary positions of "primitive" eukaryotes. Mol Biol Evol. 2006, 23: 615-625. 10.1093/molbev/msj068.View ArticlePubMedGoogle Scholar
- Hampl V, Horner DS, Dyal P, Kulda J, Flegr J, Foster PG, Embley TM: Inference of the phylogenetic position of oxymonads based on nine genes: support for metamonada and excavata. Mol Biol Evol. 2005, 22: 2508-2518. 10.1093/molbev/msi245.View ArticlePubMedGoogle Scholar
- Cardenas ML, Cornish-Bowden A, Ureta T: Evolution and regulatory role of the hexokinases. Biochim Biophys Acta. 1998, 1401: 242-264. 10.1016/S0167-4889(97)00150-X.View ArticlePubMedGoogle Scholar
- Nowitzki U, Flechner A, Kellermann J, Hasegawa M, Schnarrenberger C, Martin W: Eubacterial origin of nuclear genes for chloroplast and cytosolic glucose-6-phosphate isomerase from spinach: sampling eubacterial gene diversity in eukaryotic chromosomes through symbiosis. Gene. 1998, 214: 205-213. 10.1016/S0378-1119(98)00229-7.View ArticlePubMedGoogle Scholar
- Qian Q, Keeling PJ: Diplonemid glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and prokaryote-to-eukaryote lateral gene transfer. Protist. 2001, 152: 193-201. 10.1078/1434-4610-00059.View ArticlePubMedGoogle Scholar
- Adje CA, Opperdoes FR, Michels PA: Molecular analysis of phosphoglycerate kinase in Trypanoplasma borreli and the evolution of this enzyme in kinetoplastida. Gene. 1998, 217: 91-99. 10.1016/S0378-1119(98)00356-4.View ArticlePubMedGoogle Scholar
- Galperin MY, Bairoch A, Koonin EV: A superfamily of metalloenzymes unifies phosphopentomutase and cofactor-independent phosphoglycerate mutase with alkaline phosphatases and sulfatases. Protein Sci. 1998, 7: 1829-1835.PubMed CentralView ArticlePubMedGoogle Scholar
- Fraser HI, Kvaratskhelia M, White MF: The two analogous phosphoglycerate mutases of Escherichia coli. FEBS Lett. 1999, 455: 344-348. 10.1016/S0014-5793(99)00910-2.View ArticlePubMedGoogle Scholar
- van der Oost J, Huynen MA, Verhees CH: Molecular characterization of phosphoglycerate mutase in archaea. FEMS Microbiol Lett. 2002, 212: 111-120.View ArticlePubMedGoogle Scholar
- Hannaert V, Brinkmann H, Nowitzki U, Lee JA, Albert MA, Sensen CW, Gaasterland T, Muller M, Michels P, Martin W: Enolase from Trypanosoma brucei, from the amitochondriate protist Mastigamoeba balamuthi, and from the chloroplast and cytosol of Euglena gracilis: pieces in the evolutionary puzzle of the eukaryotic glycolytic pathway. Mol Biol Evol. 2000, 17: 989-1000.View ArticlePubMedGoogle Scholar
- Keeling PJ, Palmer JD: Parabasalian flagellates are ancient eukaryotes. Nature. 2000, 405: 635-637. 10.1038/35015167.View ArticlePubMedGoogle Scholar
- Tracy MR, Hedges SB: Evolutionary history of the enolase gene family. Gene. 2000, 259: 129-138. 10.1016/S0378-1119(00)00439-X.View ArticlePubMedGoogle Scholar
- Keeling PJ, Palmer JD: Lateral transfer at the gene and subgenic levels in the evolution of eukaryotic enolase. Proc Natl Acad Sci U S A. 2001, 98: 10745-10750. 10.1073/pnas.191337098.PubMed CentralView ArticlePubMedGoogle Scholar
- Keeling PJ: Polymorphic insertions and deletions in parabasalian enolase genes. J Mol Evol. 2004, 58: 550-556. 10.1007/s00239-003-2577-4.View ArticlePubMedGoogle Scholar
- Schramm A, Siebers B, Tjaden B, Brinkmann H, Hensel R: Pyruvate kinase of the hyperthermophilic crenarchaeote Thermoproteus tenax: physiological role and phylogenetic aspects. J Bacteriol. 2000, 182: 2001-2009. 10.1128/JB.182.7.2001-2009.2000.PubMed CentralView ArticlePubMedGoogle Scholar
- Nevalainen L, Hrdy I, Muller M: Sequence of a Giardia lamblia gene coding for the glycolytic enzyme, pyruvate,phosphate dikinase. Mol Biochem Parasitol. 1996, 77: 217-223. 10.1016/0166-6851(96)02604-7.View ArticlePubMedGoogle Scholar
- Maldonado RA, Fairlamb AH: Cloning of a pyruvate phosphate dikinase from Trypanosoma cruzi. Mol Biochem Parasitol. 2001, 112: 183-191. 10.1016/S0166-6851(00)00362-5.View ArticlePubMedGoogle Scholar
- Andersson JO: Lateral gene transfer in eukaryotes. Cell Mol Life Sci. 2005, 62: 1182-1197. 10.1007/s00018-005-4539-z.View ArticlePubMedGoogle Scholar
- Field J, Rosenthal B, Samuelson J: Early lateral transfer of genes encoding malic enzyme, acetyl-CoA synthetase and alcohol dehydrogenases from anaerobic prokaryotes to Entamoeba histolytica. Mol Microbiol. 2000, 38: 446-455. 10.1046/j.1365-2958.2000.02143.x.View ArticlePubMedGoogle Scholar
- Sanchez LB, Galperin MY, Muller M: Acetyl-CoA synthetase from the amitochondriate eukaryote Giardia lamblia belongs to the newly recognized superfamily of acyl-CoA synthetases (Nucleoside diphosphate-forming). J Biol Chem. 2000, 275: 5794-5803. 10.1074/jbc.275.8.5794.View ArticlePubMedGoogle Scholar
- Suguri S, Henze K, Sanchez LB, Moore DV, Muller M: Archaebacterial relationships of the phosphoenolpyruvate carboxykinase gene reveal mosaicism of Giardia intestinalis core metabolism. J Eukaryot Microbiol. 2001, 48: 493-497. 10.1111/j.1550-7408.2001.tb00184.x.View ArticlePubMedGoogle Scholar
- Andersson JO, Hirt RP, Foster PG, Roger AJ: Evolution of four gene families with patchy phylogenetic distributions: influx of genes into protist genomes. BMC Evol Biol. 2006, 6: 27-10.1186/1471-2148-6-27.PubMed CentralView ArticlePubMedGoogle Scholar
- Schnarrenberger C, Martin W: Evolution of the enzymes of the citric acid cycle and the glyoxylate cycle of higher plants. A case study of endosymbiotic gene transfer. Eur J Biochem. 2002, 269: 868-883. 10.1046/j.0014-2956.2001.02722.x.View ArticlePubMedGoogle Scholar
- Doolittle WF: You are what you eat: a gene transfer ratchet could account for bacterial genes in eukaryotic nuclear genomes. Trends Genet. 1998, 14: 307-311. 10.1016/S0168-9525(98)01494-2.View ArticlePubMedGoogle Scholar
- Reeves RE: How useful is the energy in inorganic pyrophosphate?. Trends Biochem Sci. 1976, 1: 53-55.Google Scholar
- Mertens E: ATP versus pyrophosphate: glycolysis revisited in parasitic protists. Parasitol Today. 1993, 9: 122-126. 10.1016/0169-4758(93)90169-G.View ArticlePubMedGoogle Scholar
- Mertens E: Pyrophosphate-dependent phosphofructokinase, an anaerobic glycolytic enzyme?. FEBS Lett. 1991, 285: 1-5. 10.1016/0014-5793(91)80711-B.View ArticlePubMedGoogle Scholar
- Park JH, Schofield PJ, Edwards MR: Pyruvate kinase is present in Giardia intestinalis. Exp Parasitol. 1997, 87: 153-156. 10.1006/expr.1997.4206.View ArticlePubMedGoogle Scholar
- Mertens E, Van Schaftingen E, Muller M: Pyruvate kinase from Trichomonas vaginalis, an allosteric enzyme stimulated by ribose 5-phosphate and glycerate 3-phosphate. Mol Biochem Parasitol. 1992, 54: 13-20. 10.1016/0166-6851(92)90090-7.View ArticlePubMedGoogle Scholar
- Acosta H, Dubourdieu M, Quinones W, Caceres A, Bringaud F, Concepcion JL: Pyruvate phosphate dikinase and pyrophosphate metabolism in the glycosome of Trypanosoma cruzi epimastigotes. Comp Biochem Physiol B Biochem Mol Biol. 2004, 138: 347-356. 10.1016/j.cbpc.2004.04.017.View ArticlePubMedGoogle Scholar
- Matsuoka M: The gene for pyruvate, orthophosphate dikinase in C4 plants: structure, regulation and evolution. Plant Cell Physiol. 1995, 36: 937-943.PubMedGoogle Scholar
- Canback B, Andersson SG, Kurland CG: The global phylogeny of glycolytic enzymes. Proc Natl Acad Sci U S A. 2002, 99: 6097-6102. 10.1073/pnas.082112499.PubMed CentralView ArticlePubMedGoogle Scholar
- Esser C, Ahmadinejad N, Wiegand C, Rotte C, Sebastiani F, Gelius-Dietrich G, Henze K, Kretschmann E, Richly E, Leister D, Bryant D, Steel MA, Lockhart PJ, Penny D, Martin W: A genome phylogeny for mitochondria among alpha-proteobacteria and a predominantly eubacterial ancestry of yeast nuclear genes. Mol Biol Evol. 2004, 21: 1643-1660. 10.1093/molbev/msh160.View ArticlePubMedGoogle Scholar
- Martin W, Rujan T, Richly E, Hansen A, Cornelsen S, Lins T, Leister D, Stoebe B, Hasegawa M, Penny D: Evolutionary analysis of Arabidopsis, cyanobacterial, and chloroplast genomes reveals plastid phylogeny and thousands of cyanobacterial genes in the nucleus. Proc Natl Acad Sci U S A. 2002, 99: 12246-12251. 10.1073/pnas.182432999.PubMed CentralView ArticlePubMedGoogle Scholar
- Archibald JM, Rogers MB, Toop M, Ishida K, Keeling PJ: Lateral gene transfer and the evolution of plastid-targeted proteins in the secondary plastid-containing alga Bigelowiella natans. Proc Natl Acad Sci U S A. 2003, 100: 7678-7683. 10.1073/pnas.1230951100.PubMed CentralView ArticlePubMedGoogle Scholar
- Andersson JO, Roger AJ: A cyanobacterial gene in nonphotosynthetic protists--an early chloroplast acquisition in eukaryotes?. Curr Biol. 2002, 12: 115-119. 10.1016/S0960-9822(01)00649-2.View ArticlePubMedGoogle Scholar
- Ricard G, McEwan NR, Dutilh BE, Jouany JP, Macheboeuf D, Mitsumori M, McIntosh FM, Michalowski T, Nagamine T, Nelson N, Newbold CJ, Nsabimana E, Takenaka A, Thomas NA, Ushida K, Hackstein JH, Huynen MA: Horizontal gene transfer from Bacteria to rumen Ciliates indicates adaptation to their anaerobic, carbohydrates-rich environment. BMC Genomics. 2006, 7: 22-10.1186/1471-2164-7-22.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- TIGR: http://www.tigr.org/tdb/e2k1/tvg/.
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Veerassamy S, Smith A, Tillier ER: A transition probability model for amino acid substitutions from blocks. J Comput Biol. 2003, 10: 997-1010. 10.1089/106652703322756195.View ArticlePubMedGoogle Scholar
- Schmidt HA, Strimmer K, Vingron M, von Haeseler A: TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002, 18: 502-504. 10.1093/bioinformatics/18.3.502.View ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP (Phylogeny inference package), version 3.6. 2004, , Department of Genome Sciences, University of WashingtonGoogle Scholar
- PUZZLEBOOT: http://rogerlab.biochemistryandmolecularbiology.dal.ca/Software/Software.htm.
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.