Research article | Open | Published:
Evolution of four gene families with patchy phylogenetic distributions: influx of genes into protist genomes
BMC Evolutionary Biologyvolume 6, Article number: 27 (2006)
Lateral gene transfer (LGT) in eukaryotes from non-organellar sources is a controversial subject in need of further study. Here we present gene distribution and phylogenetic analyses of the genes encoding the hybrid-cluster protein, A-type flavoprotein, glucosamine-6-phosphate isomerase, and alcohol dehydrogenase E. These four genes have a limited distribution among sequenced prokaryotic and eukaryotic genomes and were previously implicated in gene transfer events affecting eukaryotes. If our previous contention that these genes were introduced by LGT independently into the diplomonad and Entamoeba lineages were true, we expect that the number of putative transfers and the phylogenetic signal supporting LGT should be stable or increase, rather than decrease, when novel eukaryotic and prokaryotic homologs are added to the analyses.
The addition of homologs from phagotrophic protists, including several Entamoeba species, the pelobiont Mastigamoeba balamuthi, and the parabasalid Trichomonas vaginalis, and a large quantity of sequences from genome projects resulted in an apparent increase in the number of putative transfer events affecting all three domains of life. Some of the eukaryotic transfers affect a wide range of protists, such as three divergent lineages of Amoebozoa, represented by Entamoeba, Mastigamoeba, and Dictyostelium, while other transfers only affect a limited diversity, for example only the Entamoeba lineage. These observations are consistent with a model where these genes have been introduced into protist genomes independently from various sources over a long evolutionary time.
Phylogenetic analyses of the updated datasets using more sophisticated phylogenetic methods, in combination with the gene distribution analyses, strengthened, rather than weakened, the support for LGT as an important mechanism affecting the evolution of these gene families. Thus, gene transfer seems to be an on-going evolutionary mechanism by which genes are spread between unrelated lineages of all three domains of life, further indicating the importance of LGT from non-organellar sources into eukaryotic genomes.
During the past five years a number of reports have appeared indicating that protists acquire genes via LGT [1–4]. Recently, phylogenomic analyses of the complete genome sequences of Entamoeba histolytica and Cryptosporidium parvum indicated that several genes of these human parasites, including some key metabolic enzymes, most likely had been acquired from prokaryotes. 96 cases of relatively recent LGT from prokaryotic sources were reported for the former and 24 for the latter [5, 6]. There are reasons to believe that LGT actually does influence protist genome evolution, since foreign genetic material is constantly entering the cell via food organisms. In addition, many protists harbour prokaryotes or eukaryotes (such as those that gave rise to secondary and tertiary plastids ) as endosymbionts. As a result, the occasional incorporation of genes from engulfed cells into the nucleus may facilitate a process of directional transfer of genes from the food organisms to phagotrophic eukaryotes over evolutionary time [4, 8]. There is a growing amount of data that are consistent with this hypothesis. For instance, LGT has mostly been detected in phagotrophic lineages [4, 9]. Moreover, the introduced genes in these lineages seem to have originated from organisms sharing the same environment with the recipient organisms – the anaerobic diplomonad lineage was found to have acquired genes from anaerobic prokaryotes in most cases , 22% of candidate donors lineages in LGT cases for Entamoeba histolytica involve relatives of the Bacteroides group which are abundant in human digestive tract , while the alga Bigelowiella natans has acquired genes mostly from other algae . These observations are consistent with the idea that physical proximity in the environment of the donor and recipient lineages may greatly enhance the probability of a successful gene transfer event , a notion recently supported by phylogenetic analyses of 144 prokaryotic proteomes identifying gene pools shared between organisms (including distantly related one) occupying the same ecological niche .
Most of the claims of LGT in protists are based on unexpected phylogenetic relationships between protist and prokaryotic sequences [2, 4]. However, phylogenetic methods are susceptible to systematic error that could lead to false interpretations of transfer events . For example, a recent phylogenetic analysis indicated that the hydrogenosomal NuoF protein from Trichomonas vaginalis (a subunit of respiratory chain complex I) branched outside of a clade of mitochondrial homologs , leading the authors to propose a separate (non-mitochondrial) origin for this protein. However, these analyses failed to take into account the heterogeneity of amino acid (aa) composition displayed by sequences in this dataset . In contrast, when the dataset is analysed with methods designed to avoid this potential artefact, the T. vaginalis sequence branched within the mitochondrial cluster , in agreement with the well-supported hypothesis that Trichomonas hydrogenosomes share an evolutionary origin with mitochondria [15, 16]. Similarly, Cpn60 phylogenies with different taxonomic samplings led to important differences in the phylogenetic relationships amongst anaerobic protists including E. histolytica and two diplomonads (Giardia and Spironucleus) eliminating the possibility of an LGT event between Entamoeba and Giardia lineages[2, 17]. In both these cases, extreme divergence coupled with compositional biases in these sequences suggested, correctly, their unexpected branching patterns were due to phylogenetic artefacts. In contrast, the phylogenetic analyses of the alanyl and prolyl tRNA synthetases show the expected phylogenetic relationships amongst prokaryotes and eukaryotes with the exception that several protist sequences were found nested within Archaea as sisters to the Nanoarchaeota sequences. In this case, the observations could not be attributed to any known phylogenetic artefacts and were most easily explained in  as gene transfer events from the archaeal lineage to the protists.
Interpretations of phylogenetic analyses of proteins with a more patchy distribution in the tree of life are more challenging than the cases described above. For example, gene duplications followed by differential gene loss may also yield the unexpected phylogenetic relationships that are hallmarks of LGT. In addition, genes with a patchy distribution may only be present in one or a few lineages in each organismal group making it potentially more difficult to identify donor and recipient lineages of gene transfer events since such assignments require that recipient lineages are nested within the donor group. Fortunately, the number of complete genome sequences is steadily growing, and should clarify the patterns of gene distribution within the tree of life. In combination with thorough phylogenetic studies, analyses of the presence and absence of genes in completely sequenced genomes should be very able to differentiate putative cases of gene transfers in gene families with a patchy phylogenetic distribution from other scenarios .
To investigate whether phylogenetic artefacts, and/or unappreciated gene duplication and loss events, have influenced previous interpretations of LGT, we have broadened the taxon sampling of four gene families with patchy phylogenetic distributions, previously implicated in gene transfer events in diplomonads and E. histolytica . The updated datasets – priS (encoding a hybrid-cluster protein), fprA (A-type flavoprotein), nagB (glucosamine-6-phosphate isomerase), and adhE (alcohol dehydrogenase E) – were also analysed using more sophisticated phylogenetic methods. We have previously argued that these four genes were introduced into the genomes of diplomonads and Entamoeba from different sources based on phylogenetic analyses . If these previous observations were really indicative of LGT, increased sampling of eukaryotic and prokaryotic taxa should result in an equal or increased number of distinct eukaryotic groups in the phylogenetic analyses (i.e. eukaryotes would be polyphyletic) and stronger support for tree topologies consistent with LGT. Alternatively, a different pattern is expected if the interpretation of gene transfers were based on phylogenetic artefacts and/or differential losses. In the former case, increased taxonomic sampling should, if anything, provide evidence for a common ancestry for the diplomonad and Entamoeba sequences – as improved within-clade taxonomic sampling tends to improve phylogenetic accuracy  – reducing the number of independent eukaryote groups observed. Alternatively, if the 'polyphyletic eukaryotes' pattern was due to ancient duplications and poor paralog sampling, we would expect newly sampled sequences to cluster in the different eukaryotic clades and recover mirror eukaryotic phylogenies.
To test these alternative hypotheses, we focused our active sampling of taxa to relatives of Entamoeba; the amphizoic E. moshkovskii, the turtle parasite E. terrapinae, the snake parasite E. invadens[21, 22], and the more distantly related free-living amoeboflagellate Mastigamoeba balamuthi [23, 24], and a putative relative of diplomonads; the parabasalid Trichomonas vaginalis [18, 25–27] (the cause of trichomoniasis, a sexually transmitted disease in humans ). In addition, we updated our datasets with all currently available homologous sequences in the public databases as well as from a number of ongoing genome sequencing projects of eukaryotes. Our updated phylogenies using more sophisticated models of aa substitutions in combination with analyses of the distribution pattern of the genes indicate that gene transfer hypotheses currently best explain the data.
Results and discussion
Patchy distribution in eukaryotes
In a previous study we identified diplomonad genes potentially derived from LGT . Here we have tested if alternative hypotheses and/or phylogenetic artefacts could account for these observations. We have broadened the eukaryotic taxon sampling of four genes with a limited distribution among eukaryotes, both by cloning and sequencing new genes and mining of the available sequence databases. By using this approach we are also able to refine the timing of putative LGT events with respect to organismal divergences and to gain insights into the evolution of gene families with a patchy distribution in general. All four genes were obtained from Entamoeba invadens, Entamoeba moshkovskii and Entamoeba terrapinae, and priS, fprA and nagB (partial) sequences were obtained from Mastigamoeba balamuthi with PCR using genomic DNA as template. Two T. vaginalis cDNA clones (nagB and fprA) were also completely sequenced. A T. vaginalis priS sequence and a M. balamuthi adhE sequence had appeared in the databases since the previous analysis. Furthermore, a N. gruberi priS cDNA clone was completely sequenced (see Additional File 1 for complete listing of the datasets). To further investigate the distribution of these genes in complete or nearly complete genome sequences, we performed similarity searches against available data from ongoing eukaryotic genome projects and retrieved the significant BlastP hits. We also combined these results with the information from published genomes and mapped the occurrences of the genes onto the current hypothesis of organismal relationships among eukaryotes (Figure 1) [29–31]. The four genes show a very patchy distribution often with both presences and absences within the same eukaryotic "super group". Two extreme alternative explanations may be invoked to explain these distribution pattern within eukaryotes; (i) presence of all four genes in the last common eukaryotic ancestor followed by many differential losses within the "super groups", or (ii) absence of the genes in the ancestor followed by independent gene acquisitions in all divergent lineages that possess the genes. The duplication and gene loss scenario becomes less likely the more independent convergent gene loss events need to be postulated. Therefore, phylogenetic analyses of the individual genes should help to distinguish between these hypotheses.
In our previous study we excluded sequences that showed indications of a biased aa composition to reduce the impact of phylogenetic artefacts due to compositional heterogeneity where possible – the available methods at the time assumed aa compositional homogeneity . Here we approach this potential problem by including analyses with methods and models that are designed to mitigate the potential misleading effects of compositional heterogeneity. Each aa in the alignments was recoded to the six groups of chemically related aa that commonly replace one another [32, 33], an approach identical to the recent analyses of the NuoF protein . Previously, we were also limited by the size of the datasets , since the maximum likelihood (ML) methods were very computationally demanding at the time. The release of the PHYML software solves this problem since it is able to perform bootstrap analyses of a large number of sequences (>100) in a reasonable computational time . The recently released ModelGenerator software also ensures the usage of the optimal available model for aa substitutions in the ML analyses . These advances in the field of phylogenetics enabled us to perform more detailed analyses that include all available members of each gene family. Information about the datasets and parameters for the phylogenetic analyses are listed in Additional File 2, and the phylogenetic trees with support values from the two methods are shown in Figures 2, 3, 4, 5, 6.
All datasets in the phylogenetic analyses with grouped aa using the Metropolis-coupled Markov Chain Monte Carlo (MCMC) strategy showed convergence, indicated by good agreements between the split support values of the duplicate runs (Additional File 3). Two of the alignments (the long version of glucosamine-6-phosphate isomerase and the prismane protein) also showed a good model composition fit indicated by both posterior predictive simulations and tests for homogeneity using X2 statistics and simulations to get the null distribution (p t > 0.05 and Psim > 0.05, respectively) , while the original datasets did not (Psim < 0.05) (Additional File 2). This indicates that the recoding procedure has reduced the potential misleading effects of compositional heterogeneity in these two analyses. The other three datasets (A-type flavoprotein, the short version of glucosamine-6-phosphate isomerase, and alcohol dehydrogenase E) showed low p t and Psim values (<0.05), suggesting that compositional heterogeneity might still represent a source of artefactual results in these datasets (Additional File 2). Nevertheless, none of these grouped aa datasets failed the tests of the model composition when the χ2 curve was used to get the null distribution, while two of the these three original datasets did, suggesting that the recoding procedure had improved the model fit (Additional File 2), reducing the potential for estimation biases. At the very least, these analyses complement the more "standard" ML analyses by showing what aspects of the phylogenies are robust to aa recoding and reducing any potential effects of saturation.
The updated phylogenetic analyses show sequences highly scrambled with respect to expected organismal relationships (Figures 2, 3, 4, 5, 6), as previously observed for these genes . Thus, the earlier finding that these proteins produce phylogenetic trees that are incompatible with organismal phylogenies is robust with respect to improved taxon sampling and more detailed phylogenetic analyses – the number of eukaryotic groups (polyphyly) in the trees have increased, rather than decreased. In all analyses the eukaryotic sequences are found in at least two distinct regions of the trees nested with prokaryotic sequences (A and B boxes in Figures 2, 3, 4 &6), which are separated with strong support values. These strong separations could, in principle, be due to ancient duplication events followed by a large number of differential losses. Indeed, the presence of the same prokaryotic group in both regions in several of the phylogenetic analyses – low G+C Gram positives are for example found in both box A and B in Figures 2 and 6 – superficially supports ancient duplications. Such scenarios are expected to result in phylogenetic relationships for each paralog that mirror the organismal relationships. This is not observed in our analyses (Figures 2, 3, 4 &6). Furthermore, duplication and loss scenarios require that the gene was present in multiple copies in the last common universal ancestor and retained for a long evolutionary time. Thus, to explain the patterns we observe in the phylogenies (Figures 2, 3, 4, 5, 6) a eukaryotic ancestral genome that encoded a larger number of distantly related paralogs of the four genes than present in any of the extant eukaryote genomes would have to be inferred (Figure 1). To our knowledge, no data exist supporting a universal trend for drastic genome shrinkage in a relative recent evolutionary time. Therefore, gene duplication and differential losses alone do not seem sufficient to explain the unexpected phylogenetic relationships observed in our analyses (Figures 2, 3, 4, 5, 6).
However, the number of independent gene losses has to be weighed against the possibility of a later introduction of the genes into eukaryotes by LGT events. Yet, as none of the eukaryotic groups are found nested within a natural prokaryotic group with strong bootstrap support (Figures 2, 3, 4, 5, 6), it is difficult to identify donor and recipient lineages involved in the putative LGT events. Thus, the presence of these genes in a subset of the sampled eukaryotes are neither easily explained by vertical inheritance of the genes from the common ancestor of all eukaryotes, nor by a distinct number of easily identified gene transfer events. These phylogenies need to be carefully interpreted in combination with analysis of gene distribution patterns, as well as in the context of the biology of the available organisms.
Genes for the hybrid-cluster protein (priS) have been identified from a large number of prokaryotes, as well as several eukaryotes (Figures 1 &2). However, the cellular function of the protein is not well established; potential roles in the biological nitrogen cycle [37, 38] and the adaptive response to oxidative stress  have been suggested. Although the gene is found in all three domains of life, its distribution within the domains are patchy; for example, it is relatively widespread among proteobacterial genomes, while it has only been found in a single high G+C Gram positive species and a single cyanobacterium (Figure 2). The occurrence of the gene in a large number of unrelated lineages in combination with the absence from more closely related species is most simply explained by cross-species transmission via gene transfer. Indeed, the phylogeny of the hybrid-cluster protein strongly suggests a number of intra- and inter-domain prokaryotic LGT events, with sequences from organismal groups such as proteobacteria, low G+C Gram positives, and euryarchaeota branching in several distinct regions of the tree, often branching with unrelated lineages with strong support values (Figure 2). The eukaryotes are found within two large groups of sequences including both archaeal and bacterial homologs, separated by a long and strongly supported branch (box A and B in Figure 2). One clade contains two Trichomonas sequences that are the sole eukaryotes in one of these groups (box B in Figure 2). A prokaryote-to-eukaryote LGT event affecting the parabasalid lineage after the divergence from other eukaryotes, including diplomonads, is a more parsimonious explanation for the position of the T. vaginalis sequences than loss of this version of the gene in all other eukaryotic species, provided T. vaginalis is not basal to all the other eukaryotes included in this analysis [18, 25–27]. However, the prokaryotic donor lineage for the T. vaginalis sequences is difficult to determine from the current data and analyses.
The eukaryotic sequences found in the second group within the hybrid-cluster protein phylogeny are found in four polyphyletic groups (box A in Figure 2). However, only one of these groups is separated from the other with a significant statistical support; the Entamoeba sequences form a weakly supported group with a cyanobacterial sequence which are found as a sister group of three δ-proteobacterial sequences with a posterior probability of 1.00 in the grouped aa analysis. This is suggestive of a eubacteria-to-Entamoeba LGT event, perhaps with cyanobacteria or δ-proteobacteria as the donor lineage (Figure 2). Thus, at the very least the phylogeny of the hybrid-cluster protein suggests two transfer events from prokaryotic donors into protists. Taken at face value, the tree also supports additional transfers into various protist lineages. Indeed, the diplomonad lineages are nested within proteobacterial sequences in both the ML and grouped aa Bayesian analyses, although with weak statistical support in both cases. At any rate, this observation is suggestive of a LGT event from a proteobacterium to the diplomonad lineage. Two α-proteobacterial sequences are nested within the other eukaryotic priS sequences in box A with weak bootstrap support (Figure 2), which could indicate an origin via endosymbiotic gene transfer. However, the absence of the gene in mitochondrial genomes in combination with its absence from the nuclear genome of most eukaryotes related to pelobionts, diatoms, heterolobosea, and green algae (Figure 1), makes such an origin doubtful. Still, the weakly supported separation of diplomonad and these eukaryotic sequences may be artefactual – in reality they could represent a monophyletic group that inherited this gene from their common ancestor. If so, at least eight independent losses of priS in the apicomplexan/ciliate, oomycete, land plant, parabasalid, kinetoplastid, opisthokont, mycetozoan, and Entamoeba lineages would have to be invoked (Figure 1). Since such widespread and relatively recent independent losses appear unlikely, we favour a scenario where also the N. gruberi, M. balamuthi, T. pseudonana, and C. reinhardtii sequences have been distributed by an unknown number of gene transfer events from unsampled prokaryotic lineages or between microbial eukaryotic lineages. So far the priS gene has only been found in microbial eukaryotes. This circumstantially supports the hypothesis that the absence of a germ/soma separation in unicellular organisms increases their chance of acquiring genes by LGT . We predict that additional taxon sampling will confirm the current trend of preferential presence in unicellular eukaryotes and will further clarify the origins of the eukaryotic priS genes.
The fprA gene encodes A-type flavoprotein, a protein recently inferred to play a role in the detoxification of nitric oxide and/or oxygen in E. histolytica and was suggested to derive from a relatively recent LGT event from a prokaryotic donor . Furthermore, it has been demonstrated that T. vaginalis is able to degrade nitric oxide under microaerophilic conditions, an activity proposed to be associated with the presence of A-type flavoproteins in these parasites . Again, the gene is only found in a subset of the sequenced prokaryotes – mostly species able to grow in oxygen-poor environments (Figure 3). One exception is the widespread presence of the A-type flavoprotein within cyanobacteria, possibly indicating that the protein has evolved a different function within this group. Consistent with this hypothesis the cyanobacterial sequences are well separated from the other sequences in the tree and have a unique alignment feature; they all share a ~160 aa highly conserved C-terminal extension that is absent from all other sequences. The phylogenetic analyses of the A-type flavoprotein strongly indicate that the fprA gene has been distributed between the prokaryotic groups via LGT, rather than by vertical inheritance, since many groupings of unrelated prokaryotic taxa are observed and supported by strong support values from both analyses (Figure 3).
The eukaryotes are found in two clearly separated clusters. The two diplomonad sequences are found together with a Trichomonas sequence among a mixture of eubacterial and archaeal species, indicating a prokaryote-to-eukaryote gene transfer event to a hypothetical uniquely shared ancestor of diplomonads and parabasalids (box B in Figure 3) – unless several independent gene losses are inferred among a broad range of eukaryotic lineages (Figure 1). In the ML analysis, three additional Trichomonas homologs are found weakly associated with the strongly supported grouping of the Mastigamoeba and Entamoeba sequences (box A in Figure 3), while the Entamoeba/Mastigamoeba clade is found a sister clade to the Clostridium perfringens sequences in the grouped aa analysis with a posterior probability of 0.97 (data not shown). Thus, the relationships between amoebozoan, parabasalid and clostridial sequences are uncertain. However, Trichomonas does not share a recent common ancestor with amoebozoa (Figure 1) suggesting that the fprA gene has been acquired by separate gene transfer events in these two eukaryotic lineages. Alternatively, following a prokaryotic LGT to one of the two eukaryotic lineages, a second LGT took place between an ancestor of Entamoeba and a parabasalid. Interestingly, the three T. vaginalis sequences share a ~450 aa C-terminal extension of about 39% identity with the Clostridium perfringens 3 fprA homolog (Figure 3 and Additional File 4). This sequence, and a 433 aa long Clostridium tetani sequence (an FAD-dependent pyridine nucleotide-disulphide oxidoreductase:Rubredoxin-type 38% identical to the T. vaginalis sequences) are the most similar prokaryotic sequences in the public databases, while the most similar eukaryotic sequence, an NADH dehydrogenase from E. histolytica previously identified to be of prokaryotic origin , is only 25% identical (Additional File 4). Such a taxonomic distribution of this protein domain links the Trichomonas C-terminal extensions with the Clostridium sequences rather than to the eukaryotic sequences, suggesting that they originated via a gene transfer event from a prokaryote donor. Thus, both the N-terminal and C-terminal domains of the T. vaginalis A-type flavoproteins likely have prokaryotic origins (see Additional File 4 for discussion of plausible scenarios).
The eukaryotic lineages that encode fprA are micro-aerophilic organisms that most likely have evolved from aerobic eukaryotes , and the prokaryotes found closest to the eukaryotic sequences in the tree are found in oxygen-poor environments. These observations indicate that the transfer of the gene occurred in such an environment. The putative functional role of fprA in nitric oxide detoxification  indicate that these gene transfers might represent metabolic adaptations that allowed these different eukaryotes to better survive in anoxic environments. fprA could be part of the gene pool shared between distantly-related organisms (prokaryotic or eukaryotic) that occupy the same ecological niche.
The nagB gene encodes glucosamine-6-phosphate isomerase, an enzyme which is usually about 260 aa residues in length and is required for the biosynthesis of the cyst wall in Giardia . Apart from low G+C Gram positives and γ-proteobacteria, the nagB gene is only sparsely represented in eubacteria and not yet detected in archaea (Figure 4). It is also absent from several eukaryotic lineages (Figure 1). In the phylogenetic tree of nagB, a strongly supported group including the Entamoeba, ciliate, mycetozoa, pelobiont, parabasalid and several eubacterial sequences was detected (box B in Figure 4). All these sequences, with the exception of one of the Rhodopirellula baltica paralogs, have a roughly 500 aa residue long homologous C-terminal extension of the protein with pair-wise identities above 48%, which confirms the common ancestry of these sequences (Figure 4). To increase the resolution of this group, a separate analysis was performed which only included the sequences of the long version of the protein and therefore was based on a larger number of positions in the alignment – 560 unambiguously aligned aa residues compared to 229 (Figure 5). Interestingly, in both analyses, pelobiont and Entamoeba sequences form a group with the ciliate sequences. In the ML analyses the mycetozoan Dictyostelium discoideum is found as a sister to these sequences with a bootstrap support of 56% (Figure 5), while the Rhodopirelulla baltica 4 sequence is found as the immediate outgroup to the ciliate/Mastigamoeba/Entamoeba sequences with a posterior probability of 0.45 in the grouped aa analysis (data not shown). These weakly supported and partly incongruent phylogenies could be rationalized in the following ways: (i) a phylogenetic artefact splitting the amoebozoa sequences, in combination with differential gene loss in all sampled eukaryotic genomes with only the currently sampled ciliates and amoebozoa retaining nagB (Figure 1), (ii) inter-domain gene transfers events from closely related, but yet unsampled, prokaryotes to the amoebozoa and ciliate lineages, or (iii) the presence of the long version of nagB in the common amoebozoan ancestor followed by a transfer event to the ciliate lineage. Although none of the alternatives can be excluded, we favour the third explanation, since the expected topology within the amoebozoa is recovered, albeit with only weak bootstrap support from the ML analysis (Figure 5), if a single intra-domain LGT event is inferred. Furthermore, ciliates are known to eat other protists [43, 44], indicating that a gene transfer event from an amoebozoan to a ciliate is feasible at least in principle . Further taxonomic sampling of eukaryotic and prokaryotic genomes is obviously needed, especially within the two particular eukaryotic groups concerned, to distinguish between the different plausible scenarios. In any case, the Mastigamoeba and Entamoeba sequences form a strongly supported group that indicates the presence of the gene in their common ancestor.
In contrast, the Trichomonas homolog, that encodes the long version of the enzyme, and diplomonads, that encode the short version, have distinct origins. While the parabasalid gene likely originated via a gene transfer event, possibly from a eubacteria within the Bacteroidetes/Chlorobi group (Figure 5), the source of the diplomonad genes remains uncertain. The separation from other eukaryotes in box A appear robust with strong bootstrap support from both of the analyses with all taxa (Figure 4), as well as an additional analysis where the box B sequences and prokaryotic long branches were excluded (Additional File 5). Thus, the separation of the diplomonad sequences from the eukaryotic sequences in box A is unlikely a result of long-branch attraction; an LGT event from an unsampled eubacterial lineage seems like a more likely explanation (Figure 4).
The topology of the tree relating the opisthokont sequences to other eukaryotic lineages and prokaryotes is not easy to explain simply by vertical inheritance (box A in Figure 4); the metazoan sequences are grouped together as expected, but the fungi are split into one main group and a smaller group with three budding yeast sequences. The separation between the two fungal groups is supported by both analyses (Figure 4). In fact, the budding yeasts are found with the other fungi only in 1% of the bootstrap replicates in the ML analyses (including the analysis without long branches) and never among the 2000 sampled trees in the grouped aa analysis. Furthermore, the eukaryotes within box A are never found as a monophyletic group among the 500 bootstrap replicates in any of the ML analyses, and with a posterior probability of only 0.02 in the grouped aa analysis (data not shown). Collectively, these results indicate that the fungal nagB genes likely have separate origins; a recent introduction of a nagB gene into a common ancestor of the three budding yeast lineages Debaryomyces, Candida, and Yarrowia seems like a reasonable scenario.
The Dictyostelium sequence is found as a sister to a Fusobacterium sequence in the two ML analyses (Figure 4 and Additional File 5), while the sequence is nested between the three budding yeast sequences and the other sequences in box A in Figure 4 in the grouped aa analysis with posterior probabilities for the separations of 0.95 and 0.90, respectively (data not shown). The separation to the metazoan/fungi/euglenozoan group is strong also in the ML analysis; the Dictyostelium sequence indeed never branches with this group in any of the bootstrap replicates in the full analysis or the analysis where long branches were excluded (data not shown). Accordingly, the phylogenetic analyses indicate that a gene acquisition from a prokaryotic lineage is a plausible explanation for the origin of the Dictyostelium sequence, rather than a shared ancestry with the other eukaryotic sequences within box A.
Alcohol dehydrogenase E
Alcohol dehydrogenase E is a key enzyme in the energy metabolism of type I "amitochondriate" protists (i.e. those that lack energy-producing mitochondria or hydrogenosomes) , since it catalyzes the conversion of acetyl-CoA to ethanol in a two-step reaction which oxidizes two molecules of NADH to NAD+ . We expanded the dataset with adhE genes from three additional Entamoeba species. The failure to detect an adhE homolog in T. vaginalis in the ongoing genome project  was expected, since this organism contains hydrogenosomes (type II "amitochondrial" protist), and therefore utilizes a different set of enzymes in their energy metabolism . Interestingly, alcohol dehydrogenase E genes have been detected in the anaerobic chytrid fungus Piromyces sp. E2, which indeed does contain hydrogenosomes . However, energy metabolism of chytrids is clearly different from that of type II amitochondriate protists such as Trichomonas – chytrids exhibit a bacterial-type mixed-acid fermentation . The finding of alcohol dehydrogenase E in two green algal species, where the protein functions in aerobic mitochondria, indicates that the diversity of its functional role in eukaryotes is not fully understood .
The phylogenetic tree supports our earlier interpretations that LGT has played an important role in the evolution of this gene [10, 49], with a number of strongly supported prokaryotic relationships that most easily are explained by gene transfer events (Figure 6) – the gene is only rarely found outside low G+C Gram positives and γ-proteobacteria. The additional Entamoeba sequences form a group with the E. histolytica sequence, indicating the presence of the gene in the common ancestor of these Entamoeba species, while the Mastigamoeba sequence clearly has a distinct origin from that of its amoebozoan sisters, as observed previously . The position of the Entamoeba sequences within a eubacterial group strongly suggests a prokaryote-to-eukaryote LGT event. Sequences from two diplomonads, two green algae, a single apicomplexan, and a chytrid fungus are found in the same region of the tree as the Mastigamoeba sequence (box B in Figure 6). However, relatives of these organisms are known to lack the gene (Figure 1), arguing in favour of recent independent introductions of the gene, rather than an ancestral presence followed by differential gene loss. The green algal sequences are found as sisters to the single cyanobacterial sequence (Figure 6) with moderate to strong statistical support from both analyses, indicating a transfer event between the lineages. This seems ecologically reasonable since ancestors of these lineages could have been found in the same environment. Although the observed topology could be explained by endosymbiotic gene transfer from the plastid, the fact that land plants are lacking the gene, the absence of the gene from extant plastid genomes, and the localization of the protein in green algal mitochondria  makes a gene transfer independent of the plastid endosymbiosis somewhat more likely. Also the Mastigamoeba sequence is separated from the other eukaryotic sequences in this region with strong and moderate support from the grouped aa and ML analyses, respectively (box B in Figure 6), clearly suggesting an origin via LGT, maybe from an unsampled prokaryotic lineage.
The diplomonad, fungal, and apicomplexan alcohol dehydrogenases are found in a weakly supported cluster in both analyses. This suggests eukaryote-to-eukaryote gene transfer events, although the donor and recipient lineages are difficult to infer. Indeed, it has earlier been suggested that the green algal adhE could have been acquired by the algae from parasitizing chytrid fungi or from foraminiferan hosts to endosymbiotic algae , and similar interactions between these lineages could be invoked to explain the exchange of adhE genes. Interestingly, Cryptosporidium, Piromyces and many diplomonads are all anaerobes or microaerophiles, and many share similar, if not identical, environments; the digestive tract of various mammals. This could have facilitated gene sharing via LGT between these distantly related eukaryotic lineages, although independent acquisitions from unsampled prokaryotes cannot be excluded. Interestingly, these three distantly related eukaryotic lineages most likely have adapted to an anaerobic lifestyle independently , and the putative acquisition of adhE likely represented independent metabolic adaptations to this environment. As for the priS gene (Figure 2), our prediction is that future genome sampling will only uncover adhE genes among microbial taxa since the distribution of adhE is restricted to microbial eukaryotes.
LGT events as phylogenetic markers
LGT is usually expected to confound efforts to reconstruct organismal relationships, since it decouples the historical signals in the gene sequences from organismal lineages . However, gene transfer events can also be informative in a specific case; the shared possession of a transferred gene may indicate a phylogenetic relationship between the lineages that possess the transferred gene to the exclusion of the lineages that lack it. There certainly are limitations for such interpretations; the gene could have been lost in some of the descendants of the recipient lineage and additional transfers can complicate the correct identification of donor and recipient lineages. In any case, gene transfer events are a potentially very important source of information about organismal relationships [18, 51], especially for protists where the molecular data are scarce and phylogenetic reconstructions are difficult .
For example, the phylogenetic positions of pelobionts and Entamoeba have been difficult to resolve with molecular markers. Analyses of ribosomal RNA only weakly grouped these together , while more recently, based on a number of protein markers, it was conclusively shown that these two groups share a common ancestor to the exclusion of other eukaryotes [23, 24]. Interestingly, the Entamoeba sequences strongly group together with the Mastigamoeba sequence in two of the four analyses discussed here, fprA and nagB (Figures 3 &5). This suggests that these genes were present in the ancestor of Mastigamoeba and Entamoeba, providing further support for a specific relationship between these two eukaryotic lineages. Furthermore, the higher eukaryotic taxon Amoebozoa [23, 24] is reflected in the phylogeny of one of these genes, nagB (Figure 5), provided one accepts the recovered gene phylogeny that indicates a possible gene transfer from within this group to a ciliate lineage. If robust, this branching pattern could allow one to make inferences about the relative timing of divergences within Amoebozoa and the Ciliophora. However, improved taxonomic sampling of the nagB gene within both groups of protists will be needed to solidify such inferences. Finally, the absence of fprA in the Dictyostelium discoideum genome suggests that the presence of this eubacterial gene within various Amoebozoa lineages might be used as a synapomorphy for discerning phylogenetic relationships within the group.
Similarly, diplomonads and parabasalids have been suggested to share a common ancestor, initially mainly based on weak evidence from molecular data [53, 54]. The case for this relationship was recently strengthened by the identification of two aminoacyl-tRNA synthetase genes that appear to have been transferred to a common ancestor of the two lineages  and recent phylogenetic analyses of concatenated protein alignments [25–27]. The observation of a transfer of a gene encoding A-type flavoprotein to a uniquely shared ancestor of the two lineages (Figure 3) further supports their specific relationship. The identification of these three genes of prokaryotic origin shared between diplomonads and parabasalids in other lineages within Excavata should be useful to pinpoint relationships within this poorly resolved and diverse group of eukaryotes.
The timing of transfers relative to eukaryote diversification
The relative timing of the transfers can be addressed in more detail with our increased taxon sampling (Figure 7). For reasons outlined above, probably none of the four genes was present in the last common eukaryotic ancestor indicating that all putative transfers almost certainly happened in a more recent evolutionary time (Figure 7). However, all four genes were transferred to the diplomonad lineage before the split between Giardia and Spironucleus – they branch together in the phylogenetic reconstructions (Figures 2, 3, 4 &6). With the sampling of Trichomonas homologs for three of the genes and the absence of the fourth, we now can date the transfers of priS, nagB, and adhE to after the split between diplomonads and parabasalids [18, 25–27], but before the divergence of the two major groups of diplomonads  (Figure 7). The fourth gene (fprA) was most likely introduced into the diplomonads lineage before the split of parabasalids, but after their divergence to the other eukaryotic lineages.
Similarly, the putative transfers previously found to be affecting the Entamoeba lineage  can be dated in more detail with our updated datasets. The Entamoeba sequences branch together in the phylogenetic reconstructions for all genes, indicating that the genes were present in the common ancestor of the four species included in the analysis (Figures 2, 3, 4, 5, 6). For the other two genes, priS and adhE, the separation of the Entamoeba and Mastigamoeba sequences is strongly supported (Figures 2 &6), indicating that the transfer event to the Entamoeba lineage probably happened after the split between Entamoeba and pelobionts, but before the divergence of the Entamoeba species . A similar pattern was observed for the gene encoding alanyl-tRNA synthetase, where the ancestral eukaryotic version was replaced by a homolog from the parabasalid lineage in Entamoeba after the split from the Mastigamoeba . The timing of the transfer of the priS gene is more difficult to pinpoint since the separation of the Entamoeba and Mastigamoeba sequences is only weakly supported by the bootstrap analyses (Figure 2). As mentioned above, nagB and fprA most likely were present in the common ancestor of Mastigamoeba and Entamoeba, and the recipient of the nagB gene most likely was a common ancestor also of Dictyostelium, indicating that the transfer of these genes probably were more ancient events than the transfers of priS and adhE (Figure 7). The multiple copies found in one or more of the Entamoeba species for priS, and fprA are most likely due to recent gene duplication events within the Entamoeba lineages (Figures 2 &3), a pattern also observed from the analysis of the partial genome sequence of E. invadens. Interestingly, one of the two E. invadens priS sequences has a frameshift due to an eight nucleotide long deletion in the middle of the gene (this is unlikely to be due to a methodological artefact as several different PCR products gave identical sequences). This frameshift probably reflects the dynamics of the evolution of gene families in the Entamoeba lineage with frequent gene duplication followed by inactivation of some of the paralogs by accumulation of deleterious mutations.
Among the four sampled genes, the absence of gene transfer occurring within the Entamoeba and diplomonad groups is probably an indication that the evolutionary times since the split of these respective groups are short in comparison to the time since the last common eukaryotic ancestor, rather than an indication that the rate of inter-domain transfers have decreased in more recent evolutionary time. Indeed, one of the fifteen genes in the previous analysis, pyrG which encodes CTP synthetase, was probably introduced independently into the Giardia and Spironucleus lineages . However, the data are still scarce, and additional sampling of genes from diverse protist lineages could change the inferences of the timing of the transfers of individual genes presented here (Figure 7). Nevertheless, the data from our four genes, in combination with previously published data [10, 18], support a scenario where prokaryotic genes from various lineages have been transferred into eukaryotic lineages continuously over time (Figure 7).
A link between gene transfers and feeding habits in phagotrophic protists and their shared ecological niche?
An interesting pattern where the studied protists mostly acquire genes from prokaryotes was observed (Figure 7). This may be explained by a preference for growing in prokaryote rich environments and consuming prokaryotes by the four groups of phagotrophic protists investigated in this study – diplomonads, parabasalids, pelobionts and Entamoeba, since uptake of DNA from ingested cells is possibly an important mechanism enabling LGT in eukaryotes [4, 8]. Indeed, diplomonads generally feed on prokaryotes  and several prokaryote-to-eukaryote gene transfer events have been described for this eukaryotic group (Figure 7) [10, 41, 58, 59], while, to our knowledge, no strong case of gene transfer event from a eukaryote lineage to diplomonads has been described yet. Entamoeba, on the other hand, is able to ingest both prokaryotes and eukaryotes; it can be maintained in monoxenic cultures with bacteria as well as trypanosomatid flagellates [60, 61]. The Entamoeba lineage was recently suggested as the recipient lineage in a eukaryote-to-eukaryote gene transfer event of the alanyl-tRNA synthetase gene from the parabasalid lineage  (Figure 7), although most gene transfer events affecting Entamoeba seem to involve prokaryotic donor lineages (Figures 2, 3, 4, 5, 6) [5, 10].
In this study, the donor and recipient lineages could be inferred in one putative eukaryote-to-eukaryote gene transfer event with reasonable support; ciliates were hypothesized to have acquired a gene from an Amoebozoa lineage (Figures 5 and 7). Interestingly, ciliates were also previously shown to represent the recipient lineage in an intra-domain transfer of the alanyl-tRNA synthetase gene  (Figure 7). It is possible that the recipient lineage of these two transfers – an ancestor of Paramecium and Tetrahymena – tended to preferentially graze on eukaryotic protists rather than bacteria and was therefore exposed to eukaryotic DNA leading to LGT events; ciliates are indeed known to eat both prokaryotes and eukaryotes [43, 44]. Similarly, dinoflagellates are known to graze on eukaryotes  and have been identified as the recipient lineage in eukaryote-to-eukaryote gene transfer events [63, 64]. If this pattern holds up in light of more data, it suggests that there is a link between the genome evolution and the food content in phagotrophic protists – indicating that an understanding of eating habits is important to our understanding of gene transfer in the evolution of phagotrophic protists, as postulated by Doolittle . Global proteome phylogenies from 144 prokaryotes indicate that LGT has created pools of shared genes between distantly related prokaryotes occupying the same niche , such as mammalian mucosa. Our present analyses extend these observations to microbial eukaryotes with shared genes between microorganisms thriving on mammalian mucosa such as the trichomonads, diplomonads, apicomplexans, Piromyces and Entamoeba.
Sources of DNA
Entamoeba invadens (strain IP-1, ATCC 30994), E. moshkovskii (strain FIC, ATCC 30041), E. terrapinae (strain M, ATCC 30043) were cultured in LYI-S-2 medium at room temperature , Mastigamoeba balamuthi (ATCC 30984) were cultured in PYGC medium , and Naegleria gruberi (strain NEG-M, ATCC 30224) were cultured in modified PYNFH medium (ATCC medium 1034, cat. no. 327-X) at room temperature. The cells were harvested, followed by lysis in 0.25% SDS/0.1 M EDTA, pH 8.0, and genomic DNAs were purified using a cetyl trimethylammonium bromide extraction (CTAB) method . Trichomonas vaginalis (strain G3) nagB and fprA cDNA clones were identified in the ongoing EST project (Hirt, R.P., Embley, T.M., and Harriman, N.), and two priS Naegleria gruberi (strain NEG-M, ATCC 30224) cDNA clones were identified in an ongoing EST project (Sjögren, Å.M., Andersson, J.O., Gill, E., Roger, A.J., unpublished).
PCR and sequencing
Exact match PCR primers to obtain Entamoeba genes for independent sequencing were designed based on sequences available from the website of the ongoing genome projects  that showed similarity to the studied genes. If no such sequences were available, degenerate primers were designed against conserved regions of the alignments. Using different combinations of these primers the four genes were successfully amplified using genomic DNA in PCR reactions. The Mastigamoeba priS sequence was amplified using degenerate primers and genomic DNA, while the Mastigamoeba nagB and fprA were PCR amplified from genomic DNA using exact-match primers designed from cDNA sequences available from an in-house EST project and published  cDNA projects, respectively. The N. gruberi priS sequence was amplified from genomic DNA using exact match PCR primers based on cDNA sequences. The PCR products were purified using the Qiaquick PCR Purification Kit (Qiagen Inc., Valencia, CA) and directly sequenced using the ABI PRISM BigDye Termination Cycle Sequencing Kit (Applied Biosystems, Foster City, CA) using the primers used in the amplification as well as internal primers. Introns were identified and removed before the subsequent phylogenetic analyses from the obtained M. balamuthi nagB, priS, and fprA sequences.
Assembly of the datasets
All available homologs for the four genes were retrieved from the National Center for Biotechnology Information . In addition, similarity searches against sequences released from ongoing eukaryotic genome projects were performed to identify and retrieve additional eukaryotic homologs of the genes from ongoing genome projects at various genome sequences centres. Additional priS, and fprA T. vaginalis sequences, and nagB sequences from Tetrahymena thermophila, Trypanosoma. bruci, and Trypanosoma cruzi were retrieved from the Institute for Genomic Research . Chlamydomonas reinhardtii priS and adhE sequences, and a Thalassiosira pseudonana priS sequence were retrieved from the DOE Joint Genome Institute . Finally, a nagB sequences from Leishmania major and Paramecium tetraurelia were retrieved from the Wellcome Trust Sanger Institute , and Genoscope , respectively. None of the four genes could be detected among the available Phytophthora sojae sequences .
The aa sequence datasets were aligned using CLUSTALW , manually adjusted, and visually inspected to identify unambiguously aligned regions suitable for phylogenetic reconstructions. Only one sequence among pairs with >95% aa sequence identity within the unambiguously aligned regions were retained for further analyses. Finally, all sequences that covered less than one third of the unambiguously aligned regions were excluded. The accessions numbers and other details of the sequences within the datasets are listed in Additional File 1, and the alignments used in the study are available as Additional Files 6, 7, 8, 9. In the previously analyses data from unpublished prokaryotic genome projects were included , while only published prokaryotic sequences were included in the present analyses. Therefore, the Carboxydothermus sequence and the single Clostridium sequence, and the Carboxydothermus and Fibrobacter sequences are present in the previous priS and fprA analyses, respectively, but missing from the current datasets.
The optimal aa substitution model for each dataset was selected using the program ModelGenerator, recently developed by T. M. Keane, T. J. Naughton, and J. O. McInerney, National University of Ireland, Maynooth, Ireland . Protein maximum likelihood (ML) phylogenies were inferred using PHYML, version 2.4.4  with the optimal substitution model, and bootstrap support values were calculated based on 500 resampled datasets. Most currently available phylogenetic methods cannot deal with strong aa compositional heterogeneity in the data , which may lead to false interpretations of evolutionary events[2, 14]. Therefore, we also performed phylogenetic analyses using an approach that ameliorates (or mitigates) the compositional heterogeneity, as recently described in the analyses of the NuoF protein . The aa alignments were recoded into six categories corresponding to the PAM matrix (and most other matrices) as follows: (1) ASTGP, (2) DNEQ, (3) RKH, (4) MVIL, (5) FYW and (6) C . Bayesian phylogenetic analyses were performed on these grouped aa alignments using the program p4 . This allowed the use of a 6 × 6 general time-reversible rate matrix with free parameters rather than a fixed empirical matrix. The among-site rate variation (ASRV) was chosen by the Akaike Information Criterion (AIC) based on ML on the neighbor-joining tree. All parameters, including the composition and substitution rate matrix, were free, and the analysis used the Metropolis-coupled MCMC strategy from MrBayes . Runs were done in duplicate for 106 generations. The first halves were discarded as burn-in, while the second halves of both runs were combined for calculating the consensus trees. Convergence was assessed by plotting the split support values >0.10 in the two independent runs against each other (Additional File 3). One dataset (the short version of glucosamine-6-phosophate) did not converge well (data not shown) and was rerun in duplicate for 2 × 106 generations, and the analyses converged under these conditions (Additional File 3). The model fit of the composition was assessed using different approaches. Posterior predictive simulations were performed on the grouped aa datasets . If the tail area probabilities (p t ) is low (<0.05), the model composition does not fit the composition of the dataset. Also, tests for compositional homogeneity using X2 statistics were performed for both the original non-grouped and the grouped aa datasets using simulations to get the expected null distribution from the obtained trees and preferred substitution models used in the analyses . The more widely used X2-tests for compositional homogeneity using the χ2 curve as the null distributions are only included for comparison since they fail to take the tree-based correlation of compositions among taxa into account  (Additional File 2).
Nucleotide sequence accession numbers
The sequences reported here were deposited in GenBank under the accession numbers AJ864541-AJ864560.
Doolittle WF, Boucher Y, Nesbø CL, Douady CJ, Andersson JO, Roger AJ: How big is the iceberg of which organellar genes in nuclear genomes are but the tip?. Philos Trans R Soc Lond B Biol Sci. 2003, 358 (1429): 39-58. 10.1098/rstb.2002.1185.
Richards TA, Hirt RP, Williams BA, Embley TM: Horizontal gene transfer and the evolution of parasitic protozoa. Protist. 2003, 154 (1): 17-32. 10.1078/143446103764928468.
Gogarten JP: Gene transfer: gene swapping craze reaches eukaryotes. Curr Biol. 2003, 13 (2): R53-R54. 10.1016/S0960-9822(02)01426-4.
Andersson JO: Lateral gene transfer in eukaryotes. Cell Mol Life Sci. 2005, 62 (11): 1182-1197. 10.1007/s00018-005-4539-z.
Loftus B, Anderson I, Davies R, Alsmark UCM, Samuelson J, Amedeo P, Roncaglia P, Berriman M, Hirt RP, Mann BJ, Nozaki T, Suh B, Pop M, Duchene M, Ackers J, Tannich E, Leippe M, Hofer M, Bruchhaus I, Willhoeft U, Bhattacharya A, Chillingworth T, Churcher C, Hance Z, Harris B, Harris D, Jagels K, Moule S, Mungall K, Ormond D, Squares R, Whitehead S, Quail MA, Rabbinowitsch E, Norbertczak H, Price C, Wang Z, Guillen N, Gilchrist C, Stroup SE, Bhattacharya S, Lohia A, Foster PG, Sicheritz-Ponten T, Weber C, Singh U, Mukherjee C, El-Sayed NM, Petri WAJ, Clark CG, Embley TM, Barrell B, Fraser CM, Hall N: The genome of the protist parasite Entamoeba histolytica. Nature. 2005, 433 (7028): 865-868. 10.1038/nature03291.
Huang J, Mullapudi N, Lancto CA, Scott M, Abrahamsen MS, Kissinger JC: Phylogenomic evidence supports past endosymbiosis, intracellular and horizontal gene transfer in Cryptosporidium parvum. Genome Biol. 2004, 5 (11): R88-10.1186/gb-2004-5-11-r88.
Bhattacharya D, Yoon HS, Hackett JD: Photosynthetic eukaryotes unite: endosymbiosis connects the dots. Bioessays. 2004, 26 (1): 50-60. 10.1002/bies.10376.
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.
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 (13): 7678-7683. 10.1073/pnas.1230951100.
Andersson JO, Sjögren ÅM, Davis LAM, Embley TM, Roger AJ: Phylogenetic analyses of diplomonad genes reveal frequent lateral gene transfers affecting eukaryotes. Curr Biol. 2003, 13 (2): 94-104. 10.1016/S0960-9822(03)00003-4.
Smith MW, Feng DF, Doolittle RF: Evolution by acquisition: the case for horizontal gene transfers. Trends Biochem Sci. 1992, 17 (12): 489-493. 10.1016/0968-0004(92)90335-7.
Beiko RG, Harlow TJ, Ragan MA: Highways of gene sharing in prokaryotes. Proc Natl Acad Sci U S A. 2005, 102 (40): 14332-14337. 10.1073/pnas.0504068102.
Dyall SD, Yan W, Delgadillo-Correa MG, Lunceford A, Loo JA, Clarke CF, Johnson PJ: Non-mitochondrial complex I proteins in a hydrogenosomal oxidoreductase complex. Nature. 2004, 431 (7012): 1103-1107. 10.1038/nature02990.
Hrdy I, Hirt RP, Dolezal P, Bardonova L, Foster PG, Tachezy J, Embley TM: Trichomonas hydrogenosomes contain the NADH dehydrogenase module of mitochondrial complex I. Nature. 2004, 432 (7017): 618-622. 10.1038/nature03149.
Embley TM, van der Giezen M, Horner DS, Dyal PL, Foster P: Mitochondria and hydrogenosomes are two forms of the same fundamental organelle. Philos Trans R Soc Lond B Biol Sci. 2003, 358 (1429): 191-201. 10.1098/rstb.2002.1190.
Martin W: The missing link between hydrogenosomes and mitochondria. Trends Microbiol. 2005
Roger AJ, Svärd SG, Tovar J, Clark CG, Smith MW, Gillin FD, Sogin ML: A mitochondrial-like chaperonin 60 gene in Giardia lamblia: evidence that diplomonads once harbored an endosymbiont related to the progenitor of mitochondria. Proc Natl Acad Sci U S A. 1998, 95: 229-234. 10.1073/pnas.95.1.229.
Andersson JO, Sarchfield SW, Roger AJ: Gene transfers from Nanoarchaeota to an ancestor of diplomonads and parabasalids. Mol Biol Evol. 2005, 22 (1): 85-90. 10.1093/molbev/msh254.
Andersson JO, Roger AJ: Evolution of glutamate dehydrogenase genes: evidence for lateral gene transfer within and between prokaryotes and eukaryotes. BMC Evol Biol. 2003, 3: 14-10.1186/1471-2148-3-14.
Pollock DD, Zwickl DJ, McGuire JA, Hillis DM: Increased taxon sampling is advantageous for phylogenetic inference. Syst Biol. 2002, 51 (4): 664-671. 10.1080/10635150290102357.
Silberman JD, Clark CG, Diamond LS, Sogin ML: Phylogeny of the genera Entamoeba and Endolimax as deduced from small-subunit ribosomal RNA sequences. Mol Biol Evol. 1999, 16 (12): 1740-1751.
Clark CG, Diamond LS: Intraspecific variation and phylogenetic relationships in the genus Entamoeba as revealed by riboprinting. J Eukaryot Microbiol. 1997, 44 (2): 142-154.
Bapteste E, Brinkmann H, Lee JA, Moore DV, Sensen CW, Gordon P, Durufle L, Gaasterland T, Lopez P, Müller M, Philippe H: The analysis of 100 genes supports the grouping of three highly divergent amoebae: Dictyostelium, Entamoeba, and Mastigamoeba. Proc Natl Acad Sci U S A. 2002, 99 (3): 1414-1419. 10.1073/pnas.032662799.
Arisue N, Hashimoto T, Lee JA, Moore DV, Gordon P, Sensen CW, Gaasterland T, Hasegawa M, Müller M: The phylogenetic position of the pelobiont Mastigamoeba balamuthi based on sequences of rRNA and translation elongation factors EF-1a and EF-2. J Eukaryot Microbiol. 2002, 49: 1-10. 10.1111/j.1550-7408.2002.tb00332.x.
Arisue N, Hasegawa M, Hashimoto T: Root of the eukaryota tree as inferred from combined maximum likelihood analyses of multiple molecular sequence data. Mol Biol Evol. 2005, 22 (3): 409-420. 10.1093/molbev/msi023.
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 (12): 2508-2518. 10.1093/molbev/msi245.
Simpson AG, Inagaki Y, Roger AJ: Comprehensive multi-gene phylogenies of excavate protists reveal the evolutionary positions of 'primitive' eukaryotes. Mol Biol Evol. 2006, 23 (3): 615-625. 10.1093/molbev/msj068.
Schwebke JR, Burgess D: Trichomoniasis. Clin Microbiol Rev. 2004, 17 (4): 794-803. 10.1128/CMR.17.4.794-803.2004.
Baldauf SL: The deep roots of eukaryotes. Science. 2003, 300 (5626): 1703-1706. 10.1126/science.1085544.
Simpson AG, Roger AJ: The real 'kingdoms' of eukaryotes. Curr Biol. 2004, 14 (17): R693-6. 10.1016/j.cub.2004.08.038.
Adl SM, Simpson AG, Farmer MA, Andersen RA, Anderson OR, Barta JR, Bowser SS, Brugerolle G, Fensome RA, Fredericq S, James TY, Karpov S, Kugrens P, Krug J, Lane CE, Lewis LA, Lodge J, Lynn DH, Mann DG, McCourt RM, Mendoza L, Moestrup O, Mozley-Standridge SE, Nerad TA, Shearer CA, Smirnov AV, Spiegel FW, Taylor MF: The new higher level classification of eukaryotes with emphasis on the taxonomy of protists. J Eukaryot Microbiol. 2005, 52 (5): 399-451. 10.1111/j.1550-7408.2005.00053.x.
Dayhoff MO, Schwartz RM, Orcutt BC: A model of evolutionary change in proteins. Atlas of Protein Sequences and Structure. Edited by: Dayhoff MO. 1978, Washington DC , National Biomedical Research Foundation, 345-352.
Whelan S, Goldman N: A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol Biol Evol. 2001, 18 (5): 691-699.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.
Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McInerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol. 2006, 6: 29-10.1186/1471-2148-6-29.
Foster PG: Modeling compositional heterogeneity. Syst Biol. 2004, 53 (3): 485-495. 10.1080/10635150490445779.
Wolfe MT, Heo J, Garavelli JS, Ludden PW: Hydroxylamine reductase activity of the hybrid cluster protein from Escherichia coli. J Bacteriol. 2002, 184 (21): 5898-5902. 10.1128/JB.184.21.5898-5902.2002.
van den Berg WA, Hagen WR, van Dongen WM: The hybrid-cluster protein ('prismane protein') from Escherichia coli. Characterization of the hybrid-cluster protein, redox properties of the [2Fe-2S] and [4Fe-2S-2O] clusters and identification of an associated NADH oxidoreductase containing FAD and [2Fe-2S]. Eur J Biochem. 2000, 267 (3): 666-676. 10.1046/j.1432-1327.2000.01032.x.
Briolat V, Reysset G: Identification of the Clostridium perfringens genes involved in the adaptive response to oxidative stress. J Bacteriol. 2002, 184 (9): 2333-2343. 10.1128/JB.184.9.2333-2343.2002.
Sarti P, Fiori PL, Forte E, Rappelli P, Teixeira M, Mastronicola D, Sanciu G, Giuffre A, Brunori M: Trichomonas vaginalis degrades nitric oxide and expresses a flavorubredoxin-like protein: a new pathogenic mechanism?. Cell Mol Life Sci. 2004, 61 (5): 618-623. 10.1007/s00018-003-3413-8.
Nixon JEJ, Wang A, Field J, Morrison HG, McArthur AG, Sogin ML, Loftus BJ, Samuelson J: Evidence for lateral transfer of genes encoding ferredoxins, nitroreductases, NADH oxidase, and alcohol dehydrogenase 3 from anaerobic prokaryotes to Giardia lamblia and Entamoeba histolytica. Eukaryot Cell. 2002, 1 (2): 181-190. 10.1128/EC.1.2.181-190.2002.
Knodler LA, Svärd SG, Silberman JD, Davids BJ, Gillin FD: Developmental gene regulation in Giardia lamblia: first evidence for an encystation-specific promoter and differential 5' mRNA processing. Mol Microbiol. 1999, 34 (2): 327-340. 10.1046/j.1365-2958.1999.01602.x.
Pierce RW, Turner JP: Ecology of planktonic ciliates in marine food webs. Rev Aquat Sci. 1992, 6 (2): 139-181.
Radek R, Hausmann K: Phagotrophy of ciliates. Ciliates as organisms. Edited by: Hausmann K, Bradbury PC. 1996, Stuttgart , Gustav Fischer Verlag GmbH, 197-219.
Müller 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, Warren A. 1998, Dordrecht, the Netherlands , Kluwer, 109-131.
Sánchez LB: Aldehyde dehydrogenase (CoA-acetylating) and the mechanism of ethanol formation in the amitochondriate protist, Giardia lamblia. Arch Biochem Biophys. 1998, 354 (1): 57-64. 10.1006/abbi.1998.0664.
Lyons EJ, Carlton JM: Mind the gap: bridging the divide between clinical and molecular studies of the trichomonads. Trends Parasitol. 2004, 20 (5): 204-207. 10.1016/j.pt.2004.03.005.
Boxma B, Voncken F, Jannink S, van Alen T, Akhmanova A, van Weelden SW, van Hellemond JJ, Ricard G, Huynen M, Tielens AG, Hackstein JH: The anaerobic chytridiomycete fungus Piromyces sp. E2 produces ethanol via pyruvate:formate lyase and an alcohol dehydrogenase E. Mol Microbiol. 2004, 51 (5): 1389-1399. 10.1046/j.1365-2958.2003.03912.x.
Atteia A, van Lis R, Mendoza-Hernandez G, Henze K, Martin W, Riveros-Rosas H, Gonzalez-Halphen D: Bifunctional aldehyde/alcohol dehydrogenase (ADHE) in chlorophyte algal mitochondria. Plant Mol Biol. 2003, 53 (1-2): 175-188. 10.1023/B:PLAN.0000009274.19340.36.
Doolittle WF: Phylogenetic classification and the universal tree. Science. 1999, 284: 2124-2129. 10.1126/science.284.5423.2124.
Huang J, Xu Y, Gogarten JP: The presence of a haloarchaeal type tyrosyl-tRNA synthetase marks the opisthokonts as monophyletic. Mol Biol Evol. 2005, 22 (11): 2142-2146. 10.1093/molbev/msi221.
Philippe H, Germot A: Phylogeny of eukaryotes based on ribosomal RNA: long-branch attraction and models of sequence evolution. Mol Biol Evol. 2000, 17 (5): 830-834.
Cavalier-Smith T: The excavate protozoan phyla Metamonada Grasse emend. (Anaeromonadea, Parabasalia, Carpediemonas, Eopharyngia) and Loukozoa emend. (Jakobea, Malawimonas): their evolutionary affinities and new higher taxa. Int J Syst Evol Microbiol. 2003, 53 (Pt 6): 1741-1758. 10.1099/ijs.0.02548-0.
Embley TM, Hirt RP: Early branching eukaryotes?. Curr Opin Genet Dev. 1998, 8: 624-629. 10.1016/S0959-437X(98)80029-4.
Silberman JD, Simpson AGB, Kulda J, Cepicka I, Hampl V, Johnson PJ, Roger AJ: Retortamonad flagellates are closely related to diplomonads - implications for the history of mitochondrial function in eukaryote evolution. Mol Biol Evol. 2002, 19 (5): 777-786.
Wang Z, Samuelson J, Clark CG, Eichinger D, Paul J, Van Dellen K, Hall N, Anderson I, Loftus B: Gene discovery in the Entamoeba invadens genome. Mol Biochem Parasitol. 2003, 129 (1): 23-31. 10.1016/S0166-6851(03)00073-2.
Brugerolle G, Lee JJ: Order Diplomonadida. An Illustrated Guide to the Protozoa, 2nd edn. Edited by: Lee JJ, Leedale GF, Bradbury P. 2002, Lawrence, Kansas , Society of Protozoologists, 1125-1135.
Andersson JO, Roger AJ: Evolutionary analyses of the small subunit of glutamate synthase: gene order conservation, gene fusions and prokaryote-to-eukaryote lateral gene transfers. Eukaryot Cell. 2002, 1 (2): 304-310. 10.1128/EC.1.2.304-310.2002.
Suguri S, Henze K, Sánchez LB, Moore DV, Müller M: Archaebacterial relationships of the phosphoenolpyruvate carboxykinase gene reveal mosaicism of Giardia intestinalis core metabolism. J Eukaryot Microbiol. 2001, 48 (4): 493-497. 10.1111/j.1550-7408.2001.tb00184.x.
Diamond LS: Improved method for the monoxenic cultivation of Entamoeba histolytica Schaudinn, 1903 and E. histolytica-like amebae with trypanosomatids. J Parasitol. 1968, 54 (4): 715-719. 10.2307/3277027.
Clark CG, Diamond LS: Methods for cultivation of luminal parasitic protists of clinical importance. Clin Microbiol Rev. 2002, 15 (3): 329-341. 10.1128/CMR.15.3.329-341.2002.
Tillmann U: Interactions between planktonic microalgae and protozoan grazers. J Eukaryot Microbiol. 2004, 51 (2): 156-168. 10.1111/j.1550-7408.2004.tb00540.x.
Harper JT, Keeling PJ: Lateral gene transfer and the complex distribution of insertions in eukaryotic enolase. Gene. 2004, 340 (2): 227-235. 10.1016/j.gene.2004.06.048.
Takishita K, Ishida K, Maruyama T: An enigmatic GAPDH gene in the symbiotic dinoflagellate genus Symbiodinium and its related species (the order Suessiales): possible lateral gene transfer between two eukaryotic algae, dinoflagellate and euglenophyte. Protist. 2003, 154 (3-4): 443-454. 10.1078/143446103322454176.
Chavez LA, Balamuth W, Gong T: A light and electron microscopical study of a new, polymorphic free-living amoeba, Phreatamoeba balamuthi n. g., n. sp. J Protozool. 1986, 33 (3): 397-404.
Clark CG: DNA purification from polysaccharide-rich cells. Protocols in Protozoology, Vol 1. Edited by: Lee JJ, Soldo AT. 1992, Lawrence, Kansas , Allen Press, D-3.1-D-3.2.
Protozoan Genomes. [http://www.sanger.ac.uk/Projects/Protozoa/]
National Center for Biotechnology Information. [http://www.ncbi.nlm.nih.gov/]
The Institute for Genomic Research. [http://www.tigr.org/]
DOE Joint Genome Institute. [http://www.jgi.doe.gov/]
The Wellcome Trust Sanger Institute. [http://www.sanger.ac.uk/]
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680.
Foster PG, Hickey DA: Compositional bias may affect both DNA-based and protein-based phylogenetic reconstructions. J Mol Evol. 1999, 48: 284-290. 10.1007/PL00006471.
Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19 (12): 1572-1574. 10.1093/bioinformatics/btg180.
We thank C. Graham Clark (London School of Hygiene and Tropical Medicine) for generous gifts of cell lysates of three Entamoeba species, T. Martin Embley (University of Newcastle upon Tyne) for sharing the Trichomonas cDNA clones, Erin Gill and Lesley Davis (Dalhousie University) for a generous gift of M. balamuthi genomic DNA and a cDNA sequence, and Miklós Müller (Rockerfeller University) for the cDNA library, Stewart Sarchfield (Dalhousie University) for experimental assistance, Åsa Sjögren (Dalhousie University) for gift of N. gruberi cDNA clones and genomic DNA, experimental assistance and critical reading of the manuscript, and Iñaki Ruiz-Trillo for critical reading of the manuscript. The availability of preliminary sequences from the various sequencing centers is greatly acknowledged.
A.J.R. is supported by a fellowship from the Alfred P. Sloan Foundation, the Peter Lougheed/Canadian Institutes for Health Research New Investigator fellowship, and the Canadian Institute for Advanced Research, Program in Evolutionary Biology. R.P.H. was supported by a Wellcome Trust University Award. This work was supported by a Canadian Institutes of Health Research (CIHR) Grant (MOP-62809) awarded to A.J.R and a Swedish Research Council (VR) Grant awarded to J.O.A.
JOA carried out the molecular biology studies, most bioinformatic analyses, and ML phylogenetic analyses and drafted the manuscript. PGF performed the grouped aa Bayesian phylogenetic analyses. RPH provided two cDNA clones and carried out the analyses for additional file 4. RPH and AJR provided advice on analyses and edited the manuscript. All authors read and approved the final manuscript.