BACKGROUND
photosynthetic eukaryotes are distributed across multiple branches of the eukaryotic tree of life.
currently, six putative 'super-groups' of eukaryotes have been proposed: opisthokonta, amoebozoa, rhizaria, excavata, chromalveolata, and archaeplastida  <cit> .
the origin of plastids  from a cyanobacterial endosymbiont is referred to as primary endosymbiosis.
primary plastid-containing eukaryotes, namely green algae and land plants, glaucophytes and red algae, are classified into archaeplastida, of which the monophyly is still debatable  <cit> .
subsequent to the evolution of primary plastids, two independent lineages of green algae were captured by two distinct lineages of phagotrophic protists via secondary  endosymbiosis, giving rise to the green secondary plastid-containing euglenophytes  and chlorarachniophytes .
chloroplast genome analyses suggest that the chlorarachniophyte plastid is derived from a green alga belonging to the ulvophyte-trebouxiophyte-chlorophyte group, while the ancestor of the euglenophyte plastid is related to prasinophyte green algae  <cit> .
red algae have also donated plastids to other eukaryotes by secondary endosymbiosis.
it is well established that secondary plastids in many 'chromalveolate' taxa are derived from red algal endosymbiont, but the origin and evolutionary history of 'chromalveolate' plastids are more controversial than those of green algal ancestry  <cit> .
chromalveolata is composed of four major sub-groups  and most sub-groups include non-photosynthetic members  <cit> .
recently, plastid-related genes and/or putative plastid-like organelles were found in several non-photosynthetic alveolate protists  <cit> , and a novel photosynthetic lineage chromerida was found and suggested to be a sister group of non-photosynthetic alveolates such as colpodellids  <cit> .
furthermore, phylogenomic analyses suggest that chromalveolata may be a paraphyletic super-group; two major chromalveolata lineages, stramenopiles and alveolates, are likely sister to another super-group, rhizaria, to the exclusion of cryptophytes and haptophytes  <cit> .
the recent description of roombia, a new katablepharid, has led to a proposal to establish hacrobia, a new taxonomic group that includes many protists formerly with uncertain taxonomical affiliations; katablepharids, biliphytes, centrohelids, telonemids and two traditional chromalveolata lineages, cryptophytes and haptophtyes  <cit> .
it was proposed that hacrobia included multiple lineages that secondarily lost photosynthetic ability, although the existence of cryptic non-photosynthetic plastids in some hacrobian lineages, like those seen in some apicomplexans, cannot be ruled out.
these data suggest that the history of plastid acquisition and loss in 'chromalveolates' is much more complicated than previously thought .
surprisingly, recent analyses of the nuclear genomes of the diatoms  showed that thousands of diatom genes are similar to those of prasinophyte green algae, an observation that was interpreted as evidence for the existence of a cryptic endosymbiont of green algal origin in a 'chromalveolate' ancestor prior to the secondary endosymbiosis that gave rise to the extant secondary plastid of red algal origin  <cit> .
these intriguing data still need to be investigated further, because the host nuclear phylogeny and the relationship between green algae and chromalveolata is unclear and the host components of these two groups may be specifically related to one another  <cit> .
the possible existence of multiple past endosymbioses is also debated in another super-group, the excavata.
euglenophytes are photosynthetic flagellates belonging to the euglenida   <cit> .
the euglenida include both heterotrophic and photoautotrophic protists, and share common ancestry with kinetoplastea, which include the human parasites trypanosoma and leishmania  <cit> .
morphological, biochemical and phylogenetic analyses suggest that only the last common ancestor of the extant plastid-harboring euglenophytes experienced the secondary endosymbiosis, but not the common ancestor of euglenozoa as a whole   <cit> .
however, the discovery of algal-type genes and the specific features of a mitochondrion-targeted protein in kinetoplastea showing similarity to those of euglenophytes led to the hypothesis that a plastid was present in the common ancestor of kinetoplastea, or euglenozoa  <cit> .
previously, we developed an automated pipeline for single gene phylogenetic tree construction and found a number of genes showing cyanobacterial ancestry in the amoeboflagellate naegleria .
one possible interpretation for the presence of these genes is that the primary endosymbiosis might have occurred in an ancestor of eukaryotes prior to the divergence of excavata  <cit> .
at any rate, among the extant excavata lineages, the presence of plastids is thus far only known in euglenophytes  <cit> .
a preliminary expressed sequence tag  analysis of the model euglenophyte euglena gracilis showed a complex history of nuclear genes in this organism  <cit> , but many aspects of how the e. gracilis nuclear genome integrated genes from the green algal endosymbiont via secondary endosymbiosis are unclear.
moreover, recent molecular phylogenies suggested the presence of 'red lineage' genes in the nuclear genome of e. gracilis, but their origins and evolutionary histories have not been explored in detail  <cit> .
here we provide phylogenetic evidence for the presence of a number of genes of non-green algal origins in e. gracilis through an expanded est survey using the laterally transferred gene mining pipeline  <cit> .
we also discuss the possible evolutionary origins of these genes via lateral and/or endosymbiotic gene transfer .
RESULTS
to understand the nature and extent of genome mosaicism in e. gracilis, we searched for e. gracilis genes showing specific affinity to homologues of photosynthetic eukaryotes other than green algae, the latter being the unambiguous source of the endosymbiont that gave rise to the secondary plastid in euglenophytes.
first, using the e. gracilis protein sequences generated from the est database as queries, we assembled a set of sequences showing strong similarity to green algal/plant proteins in preliminary phylogenetic trees.
we then identified the organisms to which the e. gracilis proteins showed the smallest distance to the query sequence on each tree.
<dig> and  <dig> e. gracilis queries were found to be 'closest' to the 'viridiplantae'  and the 'red lineage' , respectively .
we further checked the tree topologies of those putative 'red lineage'-like matches manually.
we then chose the trees where the e. gracilis query was nested in, or specifically associated with, the 'red lineage' clade as a monophyletic group with high support values.
finally we identified fourteen protein trees in which the e. gracilis sequence was monophyletic with the 'red lineage', not 'green' .
among them, four trees included e. gracilis sequences placed in the chromalveolata  clade  that is sister to red algal clades , another two trees showed monophyletic clades of chromalveolata sequences plus e. gracilis branching with green algae rather than red algae .
in the other eight trees, the euglenids are monophyletic with cr clades, but the sister group of the euglenids plus cr is unclear.
phylogenetic affiliations were sometimes not directly comparable due to differences in the availability of the gene/genome sequence data in each lineage.
nevertheless, the haptophytes  and stramenopiles are associated with, and sister to, e. gracilis in most of the trees.
two 'red lineage' proteins closely related to the e. gracilis counterparts were found in the plastid-lacking euglenid peranema trichophorum.
original cluster ids had "eel0000" followed by the 4-digit numbers shown.
the accession number is shown for the fatty acid desaturase gene, which was not found in the est database.
cr+red type genes
we found that the gene encoding homogentisate phytyltransferase  from e. gracilis branches within the cr sub-clade, which itself is monophyletic with the red algal clade .
these genes were included in the plastid and cyanobacteria type hpt gene family, and function in the vitamin e biosynthesis pathway  <cit> .
specific insertion/deletion sequences  shared with e. gracilis and the cr group supported the monophyly of these sequences.
e. gracilis was found to possess a gene encoding a hypothetical protein with a putative oxidoreductase domain.
this protein family is thus far found only in photosynthetic organisms and conserved in the 'red lineage' including the cryptophyte nucleomorphs, which are derived nuclei of red algal origin.
green plant sequences were distributed on a branch separated from the cr+red clade including the e. gracilis gene .
the glucokinase  gene from e. gracilis was found to belong to the prokaryote-type gene family, which was also conserved among primary and secondary algae .
interestingly, animals, fungi, land plants and the excavates trypanosoma and monocercomonoides possess a different type of genes for this glucokinase enzyme  <cit> , showing no sequence homology to the prokaryote-type genes.
although another excavate parasite giardia possesses a prokaryote-type glucokinase  <cit> , the giardia intestinalis counterpart was sister to eubacteria, separate from the cr+red clade in our preliminary analysis .
no land plant-like homologues were found in prasinophyte genomes.
in the e. gracilis est database, we found a short fragment showing similarity to the atp-dependent clp protease proteolytic subunit  .
the e. gracilis clpp was more similar to nucleomorph-encoded cryptophyte homologues and the plastid-encoded green algal counterparts than to mitochondrial-localized clpp homologues, suggesting that the e. gracilis clpp may function in the plastid.
although this fragment is too short  to allow construction of reliable phylogenetic trees, several characteristic amino acid residues are shared with cryptophyte nucleomorph-encoded proteins and those of pelagophytes .
cr+green type genes
the plastid calvin cycle enzyme phosphoribulokinase  is derived from the cyanobacterial ancestor of the organelle  <cit> .
molecular phylogenetic analysis of prk genes suggested that extant chromalveolata genes are related to their green algal counterparts, not red algae, and that these genes may not be derived from red algal endosymbionts engulfed in the ancestor of cr.
our phylogenetic analysis of prk genes recovered the previously reported kinship  <cit>  between green algal and chromalveolata genes, and confirmed that the e. gracilis prk nests within the chromalveolata clade .
folate/biopterin transporter  genes in cr have been shown to be most closely related to their counterparts in green plants  <cit> .
in the fbt protein tree , the e. gracilis sequence is monophyletic with proteins from perkinsus  and bigelowiella , and sister to the clade including the diatom thalassiosira and prasinophyte genes.
separation of the e. gracilis gene from these prasinophyte genes is supported with the high support values .
in the cr assemblage and land plants, this gene family is highly duplicated and divergent, and the phylogenetic patterns are complicated.
the tree topology suggests that both gene duplications in the ancestral lineages of cr and land plants, as well as more recent lineage specific duplications  have occurred.
however, no evidence for gene duplication was found in red algae.
e. gracilis genes with other affiliations
zeaxanthin epoxidase  is involved in the photoprotective xanthophyll cycle, catalyzing the addition of an epoxy group to zeaxanthin to form violaxanthin under low light conditions  <cit> .
zep genes are not found in cyanobacteria and red algae thus far, and molecular phylogenetic studies suggested that zep genes in chromalveolata might have been derived from prasinophytes via gene transfer events .
multiple duplicated zep genes are found in chromalveolata genomes, forming various subclades .
the e. gracilis gene was nested within a subclade of these genes, not affiliated with green plant gene clades.
fructose-bisphosphatase  is a key glycolytic enzyme in eukaryotes and eubacteria.
in photosynthetic eukaryotes, duplicated genes encoding fbp form another family functioning in carbohydrate metabolism such as the calvin cycle in the plastid  <cit> .
teich et al.
showed that both the plastid and cytosolic genes from e. gracilis were monophyletic with cr, not green plants  <cit> .
we confirmed with detailed analysis that both types of e. gracilis sequences are nested within the cr clades in the plastid and cytosolic fbp gene families .
however, the plastid fbp genes were duplicated in several lineages and the basal part of the cr fbp clade was not sufficiently resolved to verify the sister group to the e. gracilis plus cr clades.
the phylogenetic tree of 6-phosphogluconate dehydrogenase  of cyanobacterial ancestry shows that the euglenid gnds are monophyletic with the cr assemblage .
a gnd protein from the non-photosynthetic euglenid p. trichophorum is monophyletic with the e. gracilis counterpart and contained within the euglenid-cr clade.
thus the gnd gene is likely a synapomorphic character of euglenophytes and non-photosynthetic euglenids  <cit> .
we could not find a n-terminal extension on the e. gracilis homolog compared to cyanobacterial homologs , thus it is likely that this enzyme functions in cytosol in this organism.
we found a 'red lineage' gene encoding a conserved protein with an adp/atp transporter domain in e. gracilis and p. trichophorum .
this gene was found only in euglenids and chromalveolata among eukaryotes, in addition to several lineages of eubacteria.
stramenopiles possess duplicated forms of this gene.
we found another 'red lineage' gene, lepa, in e. gracilis.
the lepa protein, also known as elongation factor  <dig> , is a gtp-binding protein conserved among eubacteria and eukaryotes.
a previous study  <cit>  showed that eukaryotic lepa genes are distributed in two clades, a mitochondrial lepa clade and a plastid one, with the plastid clade genes being sister to cyanobacterial homologues, suggesting that they are likely derived from the cyanobacterial ancestor of the plastid via egt.
although escherichia coli lepa was shown to catalyze one-codon backward movement of ribosome complexes in vitro  <cit> , it has no apparent effect in vivo on the fidelity control of protein synthesis, rather presumably playing some role in protein folding  <cit> .
in this study, the e. gracilis protein was nested in the cr subclade  in the plastid clade .
methionine adenosyltransferase  catalyzes the synthesis of s-adenosylmethionine, which is the major methyl donor and used as a substrate in a variety of methylation reactions.
genes encoding a divergent form of mat, termed matx, have been found in euglenophytes and cr, and it is suggested that the matx genes were acquired via egt  <cit> .
although the e. gracilis matx is monophyletic with chromalveolata matx counterparts with high support values , the origin of the matx gene family is still unclear  <cit>  .
tripodi et al.
<cit>  showed that the gene encoding Δ <dig> fatty acid desaturase in e. gracilis was closely related to the homologues of thalassiosira and the labyrinthulid thraustochytrium .
although the sequence of the e. gracilis desaturase gene was not found in the est database, we extended the phylogenetic analysis of this gene/protein with currently available data from other taxa.
our results  are consistent with the previous study  <cit> .
examination of the protein sequence alignment revealed the presence of indels specific to the euglena/thalassiosira plus thraustochytrium clade that are not shared with other excavate proteins from trypanosoma, leishmania and naegleria .
discussion
evolutionary history of the 'red lineage' genes in euglenophytes
we have shown that a number of e. gracilis genes are monophyletic with 'red lineage' genes and that some of them are nested within the cr clade, suggesting that the e. gracilis genes were acquired from algae in this group via lgt.
in addition to the fact that no reports show a close phylogenetic relationship between the host cell components of euglenids and chromalveolata , the apparent absence of homologues of these e. gracilis genes in the available genome sequences of close relatives such as kinetoplastea and the heterolobosean amoeboflagellate n. gruberi further suggest that the gene transfer events occurred in an ancestral lineage of euglenids, rather than euglenozoans as a whole.
if the lgt scenario discussed above is true, what kind of eukaryote-to-eukaryote lgt could have occurred?
one possible and most likely scenario is lgt from prey to predator, as suggested in an est-based analysis of the mixotrophic chlorarachniophyte bigelowiella natans  <cit> .
although phagotrophy has apparently been lost in their phototrophic offspring, the ancestor of euglenophytes was likely a eukaryovorous euglenid protist  <cit> .
morphological and ultra-structural analyses have demonstrated that peranema-like eukaryovorous euglenids are similar in size to photosynthetic euglenophytes and are capable of ingesting algal prey by phagocytosis, which is presumably a necessary prerequisite for the eventual establishment of a secondary plastid  <cit> .
a study using video microscopy illustrated two types of feeding in p. trichophorum: engulfment of a prey cell whole  as well as sucking out the prey cell cytosol through the feeding apparatus   <cit> .
the presence of 'red lineage' genes in peranema  supports the idea that the lgt might have occurred prior to the divergence of euglenophytes and at least some eukaryovorous euglenids, but after the branching of eukaryovorous euglenids from bacteriovorous ancestors  <cit> .
in this scenario, the ancestor of euglenids might once have ingested algae with 'red lineage' secondary plastids as prey, and some genes might have been transferred from the prey to the nuclear genome of the ancestor of euglenids .
what kind of relationship might have existed between the above-mentioned lgt donors and recipients?
our data are consistent with the notion that these 'red lineage' genes were acquired multiple times independently via lgt from multiple sub-groups within the cr assemblage  <cit> .
under this view, some genes would have been lost while others happened to acquire a function in the host organism, perhaps ultimately contributing to the successful integration and functioning of the green algal endosymbiont that eventually became the current euglenophyte secondary plastid.
importantly, these genes were retained in the euglenophyte nuclear genome  without being replaced by the egt-derived genes from the secondary green-algal endosymbiont, which gave rise to the 'modern-day' secondary plastids in extant euglenophytes.
a previous study illustrated a remarkable example of lgt from non-green algae to euglenophytes.
triose-phosphate isomerase  genes from euglenophytes have been shown to be monophyletic with red algal homologues, and only distantly related to green algae, excavata and chromalveolata, suggesting an lgt from red algae to euglenophytes  <cit> .
interestingly, we found that the tim genes from the chlorarachniophyte bigelowiella natans, a member of another algal lineage harboring secondary plastid of green algal origin, were also nested within the red algal clade .
a unique insertion found in red algae and euglenophytes  <cit>  was also shared with the b. natans gene .
for consistency, we did not designate tim as a 'red lineage' gene due to the absence of chromalveolata genes in the euglenid plus red algae clade.
future analyses will hopefully uncover the true history of this enigmatic gene family.
endosymbiotic versus lateral gene transfer
phagocytosis of prey cells, endosymbiosis and, ultimately, establishment of an organelle is a complex process that is difficult to understand from genomic data alone  <cit> , and at the earliest stages of organellogenesis, the distinction between lgt and egt is blurred.
nevertheless, several lines of evidence from our study are consistent with the possibility that at least some of the 'red lineage' genes in euglenids are derived from egt rather than lgt.
first, some of the 'red lineage' genes appear to have plastid-related functions .
second, peranema possesses 'red lineage' genes  and the basic morphological characteristics believed to be necessary for establishment of a secondary plastid: peranema-like euglenids are larger than other bacteriovorous and dinema-like eukaryovorous euglenids, comparable in size and in the number of strips around the cell periphery  to several phototrophic euglenophytes such as euglena and eutreptiella  <cit> .
genomic data from peranema-like eukaryovores are currently very limited.
regardless, although morphological and phylogenetic analyses suggest that aplastidic  euglenids are paraphyletic in euglenid phylogenies, peranema-like eukaryovores are sister to phototrophic euglenophytes and possibly form a monophyletic group with them  <cit> .
third, although tree topologies are certainly dependent on taxon sampling , the genes we have characterized can be interpreted as having come from an ancestral lineage within the broader cr assemblage, rather than one or a few specific lineages.
on balance, it is thus possible that the incorporation  of an alga harboring a secondarily derived 'red lineage' plastid could have occurred in a peranema-like eukaryovorous ancestor .
under the egt scenario, retention of genes encoding plastid-targeted proteins such as prk, zep and plastid-type fbp  suggests that the secondary green plastid might have replaced the 'red lineage' tertiary plastid or succeeded it within a short period of time after plastid loss .
this may be comparable to the situation seen in the 'green' plastid-harboring dinoflagellate lepidodinium chlorophorum, where phylogenetic mosaicism of the nuclear-encoded plastid-targeted proteins derived from both the ancestral peridinin-type plastid, which is the most common plastid in dinoflagellates, and the 'new' secondary plastid of green algal origin, which is thought to have been acquired by plastid replacement  <cit> .
from our data alone, it is difficult to distinguish unambiguously between the lgt and egt hypotheses for any given gene.
the data do, however, represent an interesting case study for discussing which hypothesis is most likely and under what conditions.
if multiple phylogenetic trees suggest that the genes in question are most closely related to different organisms, the "independent lgt" scenario is most likely correct.
however, if the tree topologies are consistent with one another and all the genes appear to be derived from the same source, then "ancient egt" can be considered.
nevertheless, consistent tree topologies inferred from multiple genes/proteins do not necessarily prove that all the genes have a single origin, as phylogenetic artefacts can be misleading and taxon sampling is often insufficient to allow fine-scale resolution.
on a gene-by-gene basis, lgt and egt are indistinguishable on the basis of phylogenetic tree topology alone.
endosymbiosis can be considered a specific and extreme case among numerous types of prey-predator relationships  <cit> .
assuming that genes are repeatedly transferred from the same prey organism and very close relatives during the transition from endosymbiont to organelle, an "ancient egt" scenario entails stricter conditions on the nature of the organism engulfed by the host.
thus, based on parsimony, regardless of the number of transferred genes identified, "ancient egt" would seem to be less likely than multiple lgts form a single donor lineage, or "single-origin lgt".
however, this parsimony-based argument does not necessarily apply in every biological context.
for instance, when the host cell retains an endosymbiont that shares a recent common ancestor with the source organism of the transferred genes, "egt" would be regarded as the most parsimonious interpretation.
even when the host cell does not possess such an endosymbiont, if the frequency of the gene transfer correlates with the length of the period when the prey is retained in the host cell, increasingly large numbers of transferred genes would increasingly favor the "egt" scenario.
one can also argue that if the presumed functions of the transferred genes are assumed to be reflective of the relationship between the source organisms and the host cell, and if those genes possess plastid-related functions, "ancient egt" would also seem to be a likely scenario.
there are at present no known criteria with which to quantify and compare the probabilities of lgt and egt from genomic data, and further study will be necessary to verify whether the above assumptions are biologically reasonable.
alternative explanations and limitations in data interpretation
as noted above, combined with the potential to be misled by phylogenetic artifact, insufficient taxon sampling is also a concern when interpreting the phylogenies presented herein.
a previous study suggested that the plastid-bearing ancestor of e. gracilis is related to the prasinophyte genus pyramimonas   <cit> , from which a complete nuclear genome sequence is not yet available.
if so, the e. gracilis genes would be expected to be basally branching within the green plant clade, or branching between green plants and other primary plastid-containing lineages .
however, only two genera, ostreococcus and micromonas , were included in our database among prasinophytes, making interpretation difficult.
a recent genome-wide phylogenetic study demonstrated the presence of over  <dig> genes of apparent algal affinity, probably derived from lgt events, in the choanoflagellate monosiga brevicollis  <cit> .
the red algal-like glucokinase analyzed herein , which was not identified in the previous study  <cit> , may provide another example of lgt-derived algal genes in m. brevicollis.
nevertheless, given that other eukaryotes  possess a different type of enzyme for phosphorylating hexose  <cit> , our data do not rule out the possibility of differential loss of multiple gene families with overlapping functions early in eukaryotic evolution.
it is also important to recognize that the euglena and peranema est data are far from complete gene repertoires and that nuclear genome sequences, especially of basally branching green algae and euglenids, would be helpful to better resolve the early history of plastids in euglenophytes and other photosynthetic eukaryotes.
wider and richer taxon sampling will also help to reduce the impact of phylogenetic artifacts, e.g., long-branch attraction, stochastic variation or directional biases of evolutionary signals.
as discussed above, recent studies have suggested that the ancestor of chromalveolata possessed a considerable number of genes showing affinity to green plants  <cit> .
if chromalveolates are a monophyletic group, how would such 'green'-type genes be expected to behave in phylogenetic trees relative to the egt-/green algal-derived genes of euglenids and chlorarachniophytes?
'green'-type genes in euglenids and cr could be monophyletic due to phylogenetic artifacts  and we set aside trees in which their monophyly was weakly supported and/or the phylogenetic patterns were too ambiguous.
it is notable that phylogenetic patterns such as those seen in the fbt and zep trees  could be interpreted as a result of lgt from cr specifically to prasinophytes, which is in the opposite direction to the model proposed in previous studies  <cit> .
thus, it is important to recognize that the directionality of lgt events can be difficult to discern with confidence and greatly impact how we interpret global patterns of plastid gain and loss.
regardless, despite numerous uncertainties our data clearly indicate that euglenid nuclear genomes are evolutionary mosaics, the result of a complex past in which lgts from  cr to euglenids,  from green plants to cr , as well as  egts from red algae to cr and  from green algae to euglenids, appear to be overlaid upon the 'host lineage' phylogeny.
it is essential that our understanding of the evolutionary histories of these be reevaluated regularly and cautiously as more genomic data accumulate.
CONCLUSIONS
we have identified a number of 'red lineage' genes in the phototrophic euglenophyte e. gracilis, an organism that harbors a green algal-derived secondary plastid, as well as in the plastid-lacking eukaryovorous euglenid p. trichophorum.
it is likely that these genes have been acquired via eukaryote-to-eukaryote lgt, giving rise to a complex pattern of genome mosaicism in euglenids.
the possible sources of these genes are from prey organisms, and, possibly, the presence of a cryptic 'red lineage' tertiary endosymbiont in an ancestral euglenid.
such lgt- and/or egt-derived genes may have contributed to the successful integration and functioning of the green algal secondary plastid in modern-day euglenids.

 photosynthetic eukaryotes are distributed across multiple branches of the eukaryotic tree of life. currently, six putative'super-groups' of eukaryotes have been proposed: opisthokonta, amoebozoa, rhizaria, excavata, chromalveolata, and archaeplastida. the origin of plastids from a cyanobacterial endosymbiont is referred to as primary endosymbiosis. primary plastid-containing eukaryotes, namely green algae and land plants, glaucophytes and red algae, are classified into archaeplastida, of which the monophyly is still debatable. subsequent to the evolution of primary plastids, two independent lineages of green algae were captured by two distinct lineages of phagotrophic protists via secondary endosymbiosis, giving rise to the green secondary plastid-containing euglenophytes and chlorarachniophytes. chloroplast genome analyses suggest that the chlorarachniophyte plastid is derived from a green alga belonging to the ulvophyte-trebouxiophyte-chlorophyte group, while the ancestor of the euglenophyte plastid is related to prasinophyte green algae. red algae have also donated plastids to other eukaryotes by secondary endosymbiosis. it is well established that secondary plastids in many'chromalveolate' taxa are derived from red algal endosymbiont, but the origin and evolutionary history of'chromalveolate' plastids are more controversial than those of green algal ancestry. chromalveolata is composed of four major sub-groups and most sub-groups include non-photosynthetic members. recently, plastid-related genes and/or putative plastid-like organelles were found in several non-photosynthetic alveolate protists, and a novel photosynthetic lineage chromerida was found and suggested to be a sister group of non-photosynthetic alveolates such as colpodellids. furthermore, phylogenomic analyses suggest that chromalveolata may be a paraphyletic super-group; two major chromalveolata lineages, stramenopiles and alveolates, are likely sister to another super-group, rhizaria, to the exclusion of cryptophytes and haptophytes. the recent description of roombia, a new katablepharid, has led to a proposal to establish hacrobia, a new taxonomic group that includes many protists formerly with uncertain taxonomical affiliations; katablepharids, biliphytes, centrohelids, telonemids and two traditional chromalveolata lineages, cryptophytes and haptophtyes. it was proposed that hacrobia included multiple lineages that secondarily lost photosynthetic ability, although the existence of cryptic non-photosynthetic plastids in some hacrobian lineages, like those seen in some apicomplexans, cannot be ruled out. these data suggest that the history of plastid acquisition and loss in'chromalveolates' is much more complicated than previously thought. surprisingly, recent analyses of the nuclear genomes of the diatoms showed that thousands of diatom genes are similar to those of prasinophyte green algae, an observation that was interpreted as evidence for the existence of a cryptic endosymbiont of green algal origin in a'chromalveolate' ancestor prior to the secondary endosymbiosis that gave rise to the extant secondary plastid of red algal origin. these intriguing data still need to be investigated further, because the host nuclear phylogeny and the relationship between green algae and chromalveolata is unclear and the host components of these two groups may be specifically related to one another. the possible existence of multiple past endosymbioses is also debated in another super-group, the excavata. euglenophytes are photosynthetic flagellates belonging to the euglenida. the euglenida include both heterotrophic and photoautotrophic protists, and share common ancestry with kinetoplastea, which include the human parasites trypanosoma and leishmania. morphological, biochemical and phylogenetic analyses suggest that only the last common ancestor of the extant plastid-harboring euglenophytes experienced the secondary endosymbiosis, but not the common ancestor of euglenozoa as a whole. however, the discovery of algal-type genes and the specific features of a mitochondrion-targeted protein in kinetoplastea showing similarity to those of euglenophytes led to the hypothesis that a plastid was present in the common ancestor of kinetoplastea, or euglenozoa. previously, we developed an automated pipeline for single gene phylogenetic tree construction and found a number of genes showing cyanobacterial ancestry in the amoeboflagellate naegleria. one possible interpretation for the presence of these genes is that the primary endosymbiosis might have occurred in an ancestor of eukaryotes prior to the divergence of excavata. at any rate, among the extant excavata lineages, the presence of plastids is thus far only known in euglenophytes. a preliminary expressed sequence tag analysis of the model euglenophyte euglena gracilis showed a complex history of nuclear genes in this organism, but many aspects of how the e. gracilis nuclear genome integrated genes from the green algal endosymbiont via secondary endosymbiosis are unclear. moreover, recent molecular phylogenies suggested the presence of'red lineage' genes in the nuclear genome of e. gracilis, but their origins and evolutionary histories have not been explored in detail. here we provide phylogenetic evidence for the presence of a number of genes of non-green algal origins in e. gracilis through an expanded est survey using the laterally transferred gene mining pipeline. we also discuss the possible evolutionary origins of these genes via lateral and/or endosymbiotic gene transfer.  to understand the nature and extent of genome mosaicism in e. gracilis, we searched for e. gracilis genes showing specific affinity to homologues of photosynthetic eukaryotes other than green algae, the latter being the unambiguous source of the endosymbiont that gave rise to the secondary plastid in euglenophytes. first, using the e. gracilis protein sequences generated from the est database as queries, we assembled a set of sequences showing strong similarity to green algal/plant proteins in preliminary phylogenetic trees. we then identified the organisms to which the e. gracilis proteins showed the smallest distance to the query sequence on each tree. and e. gracilis queries were found to be'closest' to the'viridiplantae' and the'red lineage', respectively. we further checked the tree topologies of those putative'red lineage'-like matches manually. we then chose the trees where the e. gracilis query was nested in, or specifically associated with, the'red lineage' clade as a monophyletic group with high support values. finally we identified fourteen protein trees in which the e. gracilis sequence was monophyletic with the'red lineage', not'green'. among them, four trees included e. gracilis sequences placed in the chromalveolata clade that is sister to red algal clades, another two trees showed monophyletic clades of chromalveolata sequences plus e. gracilis branching with green algae rather than red algae. in the other eight trees, the euglenids are monophyletic with cr clades, but the sister group of the euglenids plus cr is unclear. phylogenetic affiliations were sometimes not directly comparable due to differences in the availability of the gene/genome sequence data in each lineage. nevertheless, the haptophytes and stramenopiles are associated with, and sister to, e. gracilis in most of the trees. two'red lineage' proteins closely related to the e. gracilis counterparts were found in the plastid-lacking euglenid peranema trichophorum. original cluster ids had "eel0000" followed by the 4-digit numbers shown. the accession number is shown for the fatty acid desaturase gene, which was not found in the est database. cr+red type genes we found that the gene encoding homogentisate phytyltransferase from e. gracilis branches within the cr sub-clade, which itself is monophyletic with the red algal clade. these genes were included in the plastid and cyanobacteria type hpt gene family, and function in the vitamin e biosynthesis pathway. specific insertion/deletion sequences shared with e. gracilis and the cr group supported the monophyly of these sequences. e. gracilis was found to possess a gene encoding a hypothetical protein with a putative oxidoreductase domain. this protein family is thus far found only in photosynthetic organisms and conserved in the'red lineage' including the cryptophyte nucleomorphs, which are derived nuclei of red algal origin. green plant sequences were distributed on a branch separated from the cr+red clade including the e. gracilis gene. the glucokinase gene from e. gracilis was found to belong to the prokaryote-type gene family, which was also conserved among primary and secondary algae. interestingly, animals, fungi, land plants and the excavates trypanosoma and monocercomonoides possess a different type of genes for this glucokinase enzyme, showing no sequence homology to the prokaryote-type genes. although another excavate parasite giardia possesses a prokaryote-type glucokinase, the giardia intestinalis counterpart was sister to eubacteria, separate from the cr+red clade in our preliminary analysis. no land plant-like homologues were found in prasinophyte genomes. in the e. gracilis est database, we found a short fragment showing similarity to the atp-dependent clp protease proteolytic subunit. the e. gracilis clpp was more similar to nucleomorph-encoded cryptophyte homologues and the plastid-encoded green algal counterparts than to mitochondrial-localized clpp homologues, suggesting that the e. gracilis clpp may function in the plastid. although this fragment is too short to allow construction of reliable phylogenetic trees, several characteristic amino acid residues are shared with cryptophyte nucleomorph-encoded proteins and those of pelagophytes. cr+green type genes the plastid calvin cycle enzyme phosphoribulokinase is derived from the cyanobacterial ancestor of the organelle. molecular phylogenetic analysis of prk genes suggested that extant chromalveolata genes are related to their green algal counterparts, not red algae, and that these genes may not be derived from red algal endosymbionts engulfed in the ancestor of cr. our phylogenetic analysis of prk genes recovered the previously reported kinship between green algal and chromalveolata genes, and confirmed that the e. gracilis prk nests within the chromalveolata clade. folate/biopterin transporter genes in cr have been shown to be most closely related to their counterparts in green plants. in the fbt protein tree, the e. gracilis sequence is monophyletic with proteins from perkinsus and bigelowiella, and sister to the clade including the diatom thalassiosira and prasinophyte genes. separation of the e. gracilis gene from these prasinophyte genes is supported with the high support values. in the cr assemblage and land plants, this gene family is highly duplicated and divergent, and the phylogenetic patterns are complicated. the tree topology suggests that both gene duplications in the ancestral lineages of cr and land plants, as well as more recent lineage specific duplications have occurred. however, no evidence for gene duplication was found in red algae. e. gracilis genes with other affiliations zeaxanthin epoxidase is involved in the photoprotective xanthophyll cycle, catalyzing the addition of an epoxy group to zeaxanthin to form violaxanthin under low light conditions. zep genes are not found in cyanobacteria and red algae thus far, and molecular phylogenetic studies suggested that zep genes in chromalveolata might have been derived from prasinophytes via gene transfer events. multiple duplicated zep genes are found in chromalveolata genomes, forming various subclades. the e. gracilis gene was nested within a subclade of these genes, not affiliated with green plant gene clades. fructose-bisphosphatase is a key glycolytic enzyme in eukaryotes and eubacteria. in photosynthetic eukaryotes, duplicated genes encoding fbp form another family functioning in carbohydrate metabolism such as the calvin cycle in the plastid. teich et al. showed that both the plastid and cytosolic genes from e. gracilis were monophyletic with cr, not green plants. we confirmed with detailed analysis that both types of e. gracilis sequences are nested within the cr clades in the plastid and cytosolic fbp gene families. however, the plastid fbp genes were duplicated in several lineages and the basal part of the cr fbp clade was not sufficiently resolved to verify the sister group to the e. gracilis plus cr clades. the phylogenetic tree of 6-phosphogluconate dehydrogenase of cyanobacterial ancestry shows that the euglenid gnds are monophyletic with the cr assemblage. a gnd protein from the non-photosynthetic euglenid p. trichophorum is monophyletic with the e. gracilis counterpart and contained within the euglenid-cr clade. thus the gnd gene is likely a synapomorphic character of euglenophytes and non-photosynthetic euglenids. we could not find a n-terminal extension on the e. gracilis homolog compared to cyanobacterial homologs, thus it is likely that this enzyme functions in cytosol in this organism. we found a'red lineage' gene encoding a conserved protein with an adp/atp transporter domain in e. gracilis and p. trichophorum. this gene was found only in euglenids and chromalveolata among eukaryotes, in addition to several lineages of eubacteria. stramenopiles possess duplicated forms of this gene. we found another'red lineage' gene, lepa, in e. gracilis. the lepa protein, also known as elongation factor, is a gtp-binding protein conserved among eubacteria and eukaryotes. a previous study showed that eukaryotic lepa genes are distributed in two clades, a mitochondrial lepa clade and a plastid one, with the plastid clade genes being sister to cyanobacterial homologues, suggesting that they are likely derived from the cyanobacterial ancestor of the plastid via egt. although escherichia coli lepa was shown to catalyze one-codon backward movement of ribosome complexes in vitro, it has no apparent effect in vivo on the fidelity control of protein synthesis, rather presumably playing some role in protein folding. in this study, the e. gracilis protein was nested in the cr subclade in the plastid clade. methionine adenosyltransferase catalyzes the synthesis of s-adenosylmethionine, which is the major methyl donor and used as a substrate in a variety of methylation reactions. genes encoding a divergent form of mat, termed matx, have been found in euglenophytes and cr, and it is suggested that the matx genes were acquired via egt. although the e. gracilis matx is monophyletic with chromalveolata matx counterparts with high support values, the origin of the matx gene family is still unclear. tripodi et al. showed that the gene encoding Δ fatty acid desaturase in e. gracilis was closely related to the homologues of thalassiosira and the labyrinthulid thraustochytrium. although the sequence of the e. gracilis desaturase gene was not found in the est database, we extended the phylogenetic analysis of this gene/protein with currently available data from other taxa. our results are consistent with the previous study. examination of the protein sequence alignment revealed the presence of indels specific to the euglena/thalassiosira plus thraustochytrium clade that are not shared with other excavate proteins from trypanosoma, leishmania and naegleria. discussion evolutionary history of the'red lineage' genes in euglenophytes we have shown that a number of e. gracilis genes are monophyletic with'red lineage' genes and that some of them are nested within the cr clade, suggesting that the e. gracilis genes were acquired from algae in this group via lgt. in addition to the fact that no reports show a close phylogenetic relationship between the host cell components of euglenids and chromalveolata, the apparent absence of homologues of these e. gracilis genes in the available genome sequences of close relatives such as kinetoplastea and the heterolobosean amoeboflagellate n. gruberi further suggest that the gene transfer events occurred in an ancestral lineage of euglenids, rather than euglenozoans as a whole. if the lgt scenario discussed above is true, what kind of eukaryote-to-eukaryote lgt could have occurred? one possible and most likely scenario is lgt from prey to predator, as suggested in an est-based analysis of the mixotrophic chlorarachniophyte bigelowiella natans. although phagotrophy has apparently been lost in their phototrophic offspring, the ancestor of euglenophytes was likely a eukaryovorous euglenid protist. morphological and ultra-structural analyses have demonstrated that peranema-like eukaryovorous euglenids are similar in size to photosynthetic euglenophytes and are capable of ingesting algal prey by phagocytosis, which is presumably a necessary prerequisite for the eventual establishment of a secondary plastid. a study using video microscopy illustrated two types of feeding in p. trichophorum: engulfment of a prey cell whole as well as sucking out the prey cell cytosol through the feeding apparatus. the presence of'red lineage' genes in peranema supports the idea that the lgt might have occurred prior to the divergence of euglenophytes and at least some eukaryovorous euglenids, but after the branching of eukaryovorous euglenids from bacteriovorous ancestors. in this scenario, the ancestor of euglenids might once have ingested algae with'red lineage' secondary plastids as prey, and some genes might have been transferred from the prey to the nuclear genome of the ancestor of euglenids. what kind of relationship might have existed between the above-mentioned lgt donors and recipients? our data are consistent with the notion that these'red lineage' genes were acquired multiple times independently via lgt from multiple sub-groups within the cr assemblage. under this view, some genes would have been lost while others happened to acquire a function in the host organism, perhaps ultimately contributing to the successful integration and functioning of the green algal endosymbiont that eventually became the current euglenophyte secondary plastid. importantly, these genes were retained in the euglenophyte nuclear genome without being replaced by the egt-derived genes from the secondary green-algal endosymbiont, which gave rise to the'modern-day' secondary plastids in extant euglenophytes. a previous study illustrated a remarkable example of lgt from non-green algae to euglenophytes. triose-phosphate isomerase genes from euglenophytes have been shown to be monophyletic with red algal homologues, and only distantly related to green algae, excavata and chromalveolata, suggesting an lgt from red algae to euglenophytes. interestingly, we found that the tim genes from the chlorarachniophyte bigelowiella natans, a member of another algal lineage harboring secondary plastid of green algal origin, were also nested within the red algal clade. a unique insertion found in red algae and euglenophytes was also shared with the b. natans gene. for consistency, we did not designate tim as a'red lineage' gene due to the absence of chromalveolata genes in the euglenid plus red algae clade. future analyses will hopefully uncover the true history of this enigmatic gene family. endosymbiotic versus lateral gene transfer phagocytosis of prey cells, endosymbiosis and, ultimately, establishment of an organelle is a complex process that is difficult to understand from genomic data alone, and at the earliest stages of organellogenesis, the distinction between lgt and egt is blurred. nevertheless, several lines of evidence from our study are consistent with the possibility that at least some of the'red lineage' genes in euglenids are derived from egt rather than lgt. first, some of the'red lineage' genes appear to have plastid-related functions. second, peranema possesses'red lineage' genes and the basic morphological characteristics believed to be necessary for establishment of a secondary plastid: peranema-like euglenids are larger than other bacteriovorous and dinema-like eukaryovorous euglenids, comparable in size and in the number of strips around the cell periphery to several phototrophic euglenophytes such as euglena and eutreptiella. genomic data from peranema-like eukaryovores are currently very limited. regardless, although morphological and phylogenetic analyses suggest that aplastidic euglenids are paraphyletic in euglenid phylogenies, peranema-like eukaryovores are sister to phototrophic euglenophytes and possibly form a monophyletic group with them. third, although tree topologies are certainly dependent on taxon sampling, the genes we have characterized can be interpreted as having come from an ancestral lineage within the broader cr assemblage, rather than one or a few specific lineages. on balance, it is thus possible that the incorporation of an alga harboring a secondarily derived'red lineage' plastid could have occurred in a peranema-like eukaryovorous ancestor. under the egt scenario, retention of genes encoding plastid-targeted proteins such as prk, zep and plastid-type fbp suggests that the secondary green plastid might have replaced the'red lineage' tertiary plastid or succeeded it within a short period of time after plastid loss. this may be comparable to the situation seen in the'green' plastid-harboring dinoflagellate lepidodinium chlorophorum, where phylogenetic mosaicism of the nuclear-encoded plastid-targeted proteins derived from both the ancestral peridinin-type plastid, which is the most common plastid in dinoflagellates, and the'new' secondary plastid of green algal origin, which is thought to have been acquired by plastid replacement. from our data alone, it is difficult to distinguish unambiguously between the lgt and egt hypotheses for any given gene. the data do, however, represent an interesting case study for discussing which hypothesis is most likely and under what conditions. if multiple phylogenetic trees suggest that the genes in question are most closely related to different organisms, the "independent lgt" scenario is most likely correct. however, if the tree topologies are consistent with one another and all the genes appear to be derived from the same source, then "ancient egt" can be considered. nevertheless, consistent tree topologies inferred from multiple genes/proteins do not necessarily prove that all the genes have a single origin, as phylogenetic artefacts can be misleading and taxon sampling is often insufficient to allow fine-scale resolution. on a gene-by-gene basis, lgt and egt are indistinguishable on the basis of phylogenetic tree topology alone. endosymbiosis can be considered a specific and extreme case among numerous types of prey-predator relationships. assuming that genes are repeatedly transferred from the same prey organism and very close relatives during the transition from endosymbiont to organelle, an "ancient egt" scenario entails stricter conditions on the nature of the organism engulfed by the host. thus, based on parsimony, regardless of the number of transferred genes identified, "ancient egt" would seem to be less likely than multiple lgts form a single donor lineage, or "single-origin lgt". however, this parsimony-based argument does not necessarily apply in every biological context. for instance, when the host cell retains an endosymbiont that shares a recent common ancestor with the source organism of the transferred genes, "egt" would be regarded as the most parsimonious interpretation. even when the host cell does not possess such an endosymbiont, if the frequency of the gene transfer correlates with the length of the period when the prey is retained in the host cell, increasingly large numbers of transferred genes would increasingly favor the "egt" scenario. one can also argue that if the presumed functions of the transferred genes are assumed to be reflective of the relationship between the source organisms and the host cell, and if those genes possess plastid-related functions, "ancient egt" would also seem to be a likely scenario. there are at present no known criteria with which to quantify and compare the probabilities of lgt and egt from genomic data, and further study will be necessary to verify whether the above assumptions are biologically reasonable. alternative explanations and limitations in data interpretation as noted above, combined with the potential to be misled by phylogenetic artifact, insufficient taxon sampling is also a concern when interpreting the phylogenies presented herein. a previous study suggested that the plastid-bearing ancestor of e. gracilis is related to the prasinophyte genus pyramimonas, from which a complete nuclear genome sequence is not yet available. if so, the e. gracilis genes would be expected to be basally branching within the green plant clade, or branching between green plants and other primary plastid-containing lineages. however, only two genera, ostreococcus and micromonas, were included in our database among prasinophytes, making interpretation difficult. a recent genome-wide phylogenetic study demonstrated the presence of over genes of apparent algal affinity, probably derived from lgt events, in the choanoflagellate monosiga brevicollis. the red algal-like glucokinase analyzed herein, which was not identified in the previous study, may provide another example of lgt-derived algal genes in m. brevicollis. nevertheless, given that other eukaryotes possess a different type of enzyme for phosphorylating hexose, our data do not rule out the possibility of differential loss of multiple gene families with overlapping functions early in eukaryotic evolution. it is also important to recognize that the euglena and peranema est data are far from complete gene repertoires and that nuclear genome sequences, especially of basally branching green algae and euglenids, would be helpful to better resolve the early history of plastids in euglenophytes and other photosynthetic eukaryotes. wider and richer taxon sampling will also help to reduce the impact of phylogenetic artifacts, e.g., long-branch attraction, stochastic variation or directional biases of evolutionary signals. as discussed above, recent studies have suggested that the ancestor of chromalveolata possessed a considerable number of genes showing affinity to green plants. if chromalveolates are a monophyletic group, how would such'green'-type genes be expected to behave in phylogenetic trees relative to the egt-/green algal-derived genes of euglenids and chlorarachniophytes?'green'-type genes in euglenids and cr could be monophyletic due to phylogenetic artifacts and we set aside trees in which their monophyly was weakly supported and/or the phylogenetic patterns were too ambiguous. it is notable that phylogenetic patterns such as those seen in the fbt and zep trees could be interpreted as a result of lgt from cr specifically to prasinophytes, which is in the opposite direction to the model proposed in previous studies. thus, it is important to recognize that the directionality of lgt events can be difficult to discern with confidence and greatly impact how we interpret global patterns of plastid gain and loss. regardless, despite numerous uncertainties our data clearly indicate that euglenid nuclear genomes are evolutionary mosaics, the result of a complex past in which lgts from cr to euglenids, from green plants to cr, as well as egts from red algae to cr and from green algae to euglenids, appear to be overlaid upon the'host lineage' phylogeny. it is essential that our understanding of the evolutionary histories of these be reevaluated regularly and cautiously as more genomic data accumulate.  we have identified a number of'red lineage' genes in the phototrophic euglenophyte e. gracilis, an organism that harbors a green algal-derived secondary plastid, as well as in the plastid-lacking eukaryovorous euglenid p. trichophorum. it is likely that these genes have been acquired via eukaryote-to-eukaryote lgt, giving rise to a complex pattern of genome mosaicism in euglenids. the possible sources of these genes are from prey organisms, and, possibly, the presence of a cryptic'red lineage' tertiary endosymbiont in an ancestral euglenid. such lgt- and/or egt-derived genes may have contributed to the successful integration and functioning of the green algal secondary plastid in modern-day euglenids. 
we found e. gracilis genes showing affinity not only with green algae, from which the secondary plastid in euglenophytes evolved, but also red algae and/or secondary algae containing red algal-derived plastids. euglenophytes are a group of photosynthetic flagellates possessing a plastid derived from a green algal endosymbiont, which was incorporated into an ancestral host cell via secondary endosymbiosis. phylogenetic analyses of these'red lineage' genes suggest that e. gracilis acquired at least genes via eukaryote-to-eukaryote lateral gene transfer from algal sources other than the green algal endosymbiont that gave rise to its current plastid. one possible explanation for the presence of these genes in these organisms is that some or all of them were independently acquired by lateral gene transfer and contributed to the successful integration and functioning of the green algal endosymbiont as a secondary plastid. alternative hypotheses include the presence of a phagocytosed alga as the single source of those genes, or a cryptic tertiary endosymbiont harboring secondary plastid of red algal origin, which the eukaryovorous ancestor of euglenophytes had acquired prior to the secondary endosymbiosis of a green alga. however, the impact of endosymbiosis on the euglenophyte nuclear genome is not fully understood due to its complex nature as a'hybrid' of a non-photosynthetic host cell and a secondary endosymbiont. our data show genome mosaicism in e. gracilis and p. trichophorum. we analyzed an est dataset of the model euglenophyte euglena gracilis using a gene mining program designed to detect laterally transferred genes.
 euglenophytes are a group of photosynthetic flagellates possessing a plastid derived from a green algal endosymbiont, which was incorporated into an ancestral host cell via secondary endosymbiosis. however, the impact of endosymbiosis on the euglenophyte nuclear genome is not fully understood due to its complex nature as a'hybrid' of a non-photosynthetic host cell and a secondary endosymbiont.  we analyzed an est dataset of the model euglenophyte euglena gracilis using a gene mining program designed to detect laterally transferred genes. we found e. gracilis genes showing affinity not only with green algae, from which the secondary plastid in euglenophytes evolved, but also red algae and/or secondary algae containing red algal-derived plastids. phylogenetic analyses of these'red lineage' genes suggest that e. gracilis acquired at least genes via eukaryote-to-eukaryote lateral gene transfer from algal sources other than the green algal endosymbiont that gave rise to its current plastid. we constructed an est library of the aplastidic euglenid peranema trichophorum, which is a eukaryovorous relative of euglenophytes, and also identified'red lineage' genes in its genome.  our data show genome mosaicism in e. gracilis and p. trichophorum. one possible explanation for the presence of these genes in these organisms is that some or all of them were independently acquired by lateral gene transfer and contributed to the successful integration and functioning of the green algal endosymbiont as a secondary plastid. alternative hypotheses include the presence of a phagocytosed alga as the single source of those genes, or a cryptic tertiary endosymbiont harboring secondary plastid of red algal origin, which the eukaryovorous ancestor of euglenophytes had acquired prior to the secondary endosymbiosis of a green alga. 