BACKGROUND
cyanobacterial ancestors gave rise to plastids  in the ancestor of a eukaryotic lineage.
the birth of the plastid had an impact on eukaryotic genome evolution, by way of endosymbiotic gene transfer , a particular form of lateral gene transfer  from endosymbionts into the phylogenetically discontiguous host genome  <cit> .
subsequently, an algal ancestor gave rise to secondary plastids in several punctate lineages of eukaryotes.
a number of these secondarily phototrophic lineages lost their photosynthetic ability and further diverged into secondarily heterotrophic, plastid-lacking protists  <cit> .
although the position of the root of eukaryotes is still uncertain, the presence of gene fusions and insertion/deletion sequences in the marker genes have allowed us to sort eukaryotes into at least three large groups; opisthokonta, amoebozoa and bikonts   <cit>  .
most phototrophic eukaryotes harboring plastids derived from primary endosymbiosis  are classified into the super-group archaeplastida   <cit> .
although it is widely accepted that primary plastids share a single origin  and the archaeplastida are monophyletic , the evolutionary history of the primary plastids is still debatable  <cit> .
in plastid-lacking protists, 'plastid imprints' can be exemplified by genomic information, i.e.
genes with affinity to extant cyanobacterial or algal genes.
these genes were supposed to have originated from egt events, and this assumption should be affirmed by the resulting phylogenetic relationship between 'imprint' genes and the extant relatives of the putative endosymbionts.
the biggest challenge and the limitation of this 'imprint' searching process is that the inevitable incompleteness of genome information on lineages of interest and the ever-developing phylogenetic methodologies make it difficult to distinguish egt and ancient lgt  <cit> .
thus, although available eukaryotic genome data are increasingly accumulating, gene and genome phylogenies should be carefully interpreted to infer evolutionary scenarios.
chromalveolata is a large taxonomic group of eukaryotes, encompassing secondary phototrophs and secondarily heterotrophic protists  <cit> , and the 'chromalveolate hypothesis' argues that this group originated from a common ancestor harboring the chlorophyll c-containing secondary plastid derived from a red alga   <cit> .
among the secondarily heterotrophic chromalveolates, several lineages have retained remnant chloroplasts for non-photosynthetic metabolic pathways, e.g.
apicoplasts in apicomplexan parasites  <cit> .
recent genomic surveys revealed the presence of plastid-derived genes, and further suggested the presence of cryptic secondary plastids in non-photosynthetic alveolate protists  <cit> .
furthermore, re-examination of the whole genome sequences suggested the existence of algal genes in ciliates, another plastid-lacking alveolate lineage, which could support the photosynthetic ancestry of ciliates  <cit> .
oomycetes are plastid-lacking stramenopiles, or chromists, classified into chromalveolata  <cit> .
although whole genome sequence analysis showed that a number of genes with affinity to photosynthetic organisms  are encoded in the nuclear genome, most of these 'plastid imprints' candidates were only suggested by similarity search and phylogenetic analyses have not yet led to fully recovering the expected tree topology  <cit> .
considering the uncertain phylogenetic affinity of the 'best hit' in similarity search  <cit> , reassessment of the genome information is important to determine whether the evolutionary history of oomycetes is comparable to ciliates  <cit> .
one candidate of 'plastid imprints' in oomycetes has been confirmed by studies reporting the phylogeny of gnd genes, which encode 6-phosphogluconate dehydrogenase, showing that some plastid-lacking protists have plant-like, cyanobacterium-derived gnd genes  <cit> .
these analyses suggested that the gnd genes with cyanobacterial ancestry were acquired early in eukaryotic evolution, either via ancient eukaryote-to-eukaryote lgt, or primary egt that occurred earlier than had ever been thought  <cit> .
additionally, the phylogeny of gnd genes demonstrated that cyanobacterial genes are also present in several excavata protists, e.g.
the heterolobosean amoebo-flagellate naegleria gruberi.
naegleria gruberi is a non-parasitic heterotrophic species related to n. fowleri, which is the causative agent of primary amoebic meningoencephalitis in mammals  <cit> .
although the phylogenetic relationship within excavata is still unclear, heterolobosea, together with jakobida, is likely to be a sister group of euglenozoa  <cit> .
to address how many genes have cyanobacterial ancestry in plastid-lacking protists, and whether cyanobacterial ancestry is limited to this gnd gene or also found in other genes, we conducted a phylogenomic analysis using genome sequence data of a taxonomically wide range of plastid-lacking eukaryotes.
here we present a gene mining study with a novel pipeline automatically producing and summarizing one-by-one phylogenetic trees, and show phylogenetic analyses of resultant candidate genes with cyanobacterial ancestry, using the whole genome sequence data from a wide range of eukaryotic lineages.
RESULTS
to address how many genes are derived from cyanobacteria in non-photosynthetic protists, we conducted cyanobacterial gene mining using the genome sequence data of a wide range of the plastid-lacking eukaryotes .
using the whole genome data, we conducted blast searches against all 'bacteria' and selected queries showing the highest similarity to genes in the available cyanobacteria genome sequences.
we then drew the neighbor-joining  trees for genes showing homology to cyanobacterial counterparts.
after the first tree construction step, we selected the gene trees where cyanobacteria and eukaryotes formed a monophyletic group excluding other prokaryotes.
as a result, we obtained a shorter list of candidates, which we termed 'genes with cyanobacterial affinity'.
subsequently we re-analyzed the eukaryotic genes with cyanobacterial affinity by visually checking and re-drawing the bayesian and maximum likelihood  trees after manually trimming operational taxonomic units .
in total, we identified  <dig> plastid-lacking protist genes 'with cyanobacterial ancestry' in the genomes of the wide range of eukaryotes: two plastid-lacking bikonts  and three unikonts  .
these were the eukaryotic genes with cyanobacterial ancestry that shared the same origin with archaeplastida and other eukaryotes.
they were placed within a monophyletic subclade mostly composed of photosynthetic organisms  and showed an apparent cyanobacterial ancestry as far as was determined by tree topology .
we found another type of gene with cyanobacterial ancestry, which were the protist genes forming monophyletic groups mostly with genes from extant cyanobacteria .
among a number of candidate genes found through the first screening, we have presented three typical trees that were resolved with significant support values .
we postulate that these prokaryote-type genes are remnants of the bacterium-to-eukaryote lgt, which occurred 'recently' in evolution.
interestingly, while phytophthora ycf <dig> homologs, probably transferred from a relative of the extant cyanobacterial species via lgt, were placed within the cyanobacterial gene clade, the ciliate tetrahymena ycf <dig> homolog showed affinity to archaeplastida .
this gene was neither found in the paramecium genome nor in the list of the recently identified algal genes in ciliates  <cit> .
jgi gene id , gene name  and accession numbers  are listed.
we found that uroporphyrin iii methyltransferase gene homologs  consisted of two large subfamilies of genes with cyanobacterial ancestry, and that oomycete genes were included only in one of them.
given that both subfamilies include green plants, red algae, chromalveolates and cyanobacteria, it is likely that they diverged within the ancestral cyanobacteria and transferred into eukaryotic hosts via primary and secondary endosymbioses.
both of the subfamilies were concurrently present in the cyanobacterial and green algal genomes.
in land plants, red algae, diatoms, haptophytes and the plastid-lacking oomycetes, one of the subfamilies might be lost along with the loss of the plastid.
the thalassiosira homolog formed a monophyletic group with green plants, rather than red algae, suggesting that it was acquired independently of the secondary plastid of the red lineage.
in this study, the bacteriovorous choanoflagellate monosiga brevicollis gene and the proteobacterial genes  were treated as 'apparently lgt-derived genes', incongruously showing affinities to photosynthetic bikonts  <cit> .
genes encoding cobalamin-independent methionine synthase in green and red algae, diatoms, and oomycetes formed a monophyletic group with cyanobacterial homologs, while the land plants and the red alga cyanidioschyzon homologs were placed in different clades unrelated to cyanobacteria .
close association between diatom and oomycete genes suggested the deep ancestry of the genes in the chromalveolate lineage.
we failed to find the homologs in the prasinophytes ostreococcus and micromonas, suggesting that this gene family was dispensable in some plant lineages.
one of the genes with cyanobacterial ancestry found in n. gruberi is pyridoxal-dependent amino acid decarboxylase gene .
the tree indicated that green plants were split into different eukaryotic clades.
naegleria and chromalveolate genes showed robust monophyly with green plants, included in a cyanobacterial gene clade.
the tree showed that land plants possessed another subfamily, associated with red algal and fungal genes, apparently of non-cyanobacterial origin.
we also identified genes with cyanobacterial ancestry from naegleria in an oxidoreductase gene family that included genes encoding rieske iron-sulfur cluster  <dig> kda protein of chloroplast inner membrane translocon , chlorophyll a oxidase , lethal-leaf spot  <dig> ) and accelerated cell death  <dig>    <cit> .
all the members of this family in land plants were hypothesized to be located at the inner membrane of the chloroplast, and to be involved in chlorophyll metabolism  <cit> .
the phylogenetic tree of the tic55-like gene family showed intricate distribution of cyanobacterial, green plant and chromalveolate genes.
in other trees of the genes identified in this study , gene clades with cyanobacterial ancestry were mostly composed of bikonts genes, besides the choanoflagellate m. brevicollis genes .
discussion
we identified eight and seven genes with cyanobacterial ancestry in the genome sequences of the oomycete p. ramorum and the heterolobosean n. gruberi, respectively .
it was reported that the apicomplexan cryptosporidium 'recently' lost their secondary plastid, and retained two to seven putative plastid-derived genes in the genome  <cit> .
this number is comparable to our result of the gene mining study using oomycete and heterolobosean genomes.
in addition, our system resolved the hidden diversity of the gene family repertoire in eukaryotic genomes by one-by-one gene phylogenies.
secondary egt scenario
although the phylogenetic positions of cryptophyceae and haptophyta are still debatable , the chromalveolate hypothesis has been reinstated to support the evolutionary scenario that the plastid-lacking protists oomycetes and ciliates once might have had a plastid  <cit> .
according to this hypothesis, the genes with cyanobacterial ancestry found in the oomycete genomes were acquired via secondary egt in the common ancestor of chromalveolata, from the red algal ancestor of secondary plastids.
this explanation is also applicable under the alternative hypothesis for chromalveolate plastids, which proposes that a tertiary endosymbiont of the haptophyte/cryptophyceae lineage is the origin of the stramenopile/alveolate plastids  <cit> .
the phylogenetic tree of the photorespiratory glycerate kinase genes, suggesting the red algal origin of the phytophthora genes , is consistent with the chromalveolate hypothesis.
however, several other gene trees in this study showed oomycete genes with green lineage affinity, not red algae .
recently, frommolt et al.
<cit>  demonstrated that, out of  <dig> genes involved in carotenoid biosynthesis from chromalveolate algae, one third  of plastid-targeted, nuclear-encoded genes are most closely related to green algal homologs.
reyes-prieto, moustafa and bhattacharya  <cit>  identified  <dig> genes of possible algal origin in the ciliates tetrahymena thermophila and paramecium tetraurelia, and 7/ <dig> of their trees show a close relationship between green plants and chromalveolata.
frommolt et al.
<cit>  attributed the close relationships between green plants and chromalveolate genes to the secondary endosymbiosis of an ancient green plant , based on the hypothesis on the monophyly of the archaeplastida  <cit> .
this explanation might be also applicable to the plant-like genes in ciliates  <cit> .
while heterolobosea and euglenozoa are often united as the morphologically defined taxon, discicristata, within excavata  <cit> , recent morphological and molecular phylogenetic analyses suggest that the heteroloboseans  never possessed the secondary plastid of green lineage and share the same origin with euglenida  <cit> .
molecular phylogenetic analyses showed that excavata is separated from other secondary plastid-containing eukaryotes   <cit> .
therefore, it is unlikely that the genes with cyanobacterial ancestry found in the heterolobosean nuclear genomes originated from the plastid cognate with any known secondary plastids in extant photosynthetic eukaryotes.
the amino acid decarboxylase gene  and the gnd gene   <cit>  trees demonstrated the presence of genes with cyanobacterial ancestry in other heterolobosean species than n. gruberi, suggesting that the ancestor of the genus naegleria possessed this gene family.
furthermore, although ml bootstrap support or bayesian posterior probability  values were not always sufficient, the naegleria genes occupy relatively basal phylogenetic positions within the bikonts clade in all seven trees .
thus it is possible that the genes with cyanobacterial ancestry were introduced en bloc in the ancestor of heterolobosea, via a batch gene transfer, in a concerted manner.
one possible origin of such a concerted gene transfer is secondary egt from a photosynthetic eukaryote with a basal phylogenetic position within bikonts.
however, as discussed above, it is unlikely that heterolobosea experienced secondary endosymbiosis and acquired genes common to the extant secondary plastid-containing eukaryotes via secondary egt.
ancient eukaryote-to-eukaryote lgt or primary egt scenarios
alternatively, we can argue for two other explanations: a concerted eukaryote-to-eukaryote lgt scenario or a more ancient primary egt scenario.
the naegleria genes with cyanobacterial ancestry shown in table  <dig> are basally positioned within bikonts, but not intruding into any of gene clades from extant photosynthetic eukaryotes .
thus, if we assume that these genes were acquired via non-endosymbiotic lgt, they may originate from unknown ancient photosynthetic lineages basally positioned within bikonts.
meanwhile, under the primary egt scenario, in which the primary endosymbiosis occurred in the common ancestor of bikonts  , ancient primary egt occurred much earlier than the conventional hypothesis, from the cyanobacterium-like prokaryote to the common ancestor of bikonts.
primary plastids were subsequently lost in many lineages of bikonts, except for the archaeplastida lineages, but some genes originating from the cyanobacterial ancestor of the primary plastids have been retained in the nuclear genomes of the plastid-lacking lineages of bikonts .
the loss of the plastid might have triggered the loss of genes that specifically functioned within the plastid.
only a portion of the plastid-derived genes, which we can find now in the plastid-lacking protist genomes, might have escaped from or survived through eliminative pressure in a lineage-specific manner, by acquiring additional functions with other components and/or in other cellular compartments.
this might account for the observed punctate distribution of gene families among the eukaryotes  <cit> .
recently, a hypothesis for the non-monophyly of archaeplastida was proposed based on the phylogenetic analyses of slowly evolving nuclear-encoded genes  <cit> .
this non-monophyly hypothesis could be also considered within the scope of the primary egt scenario.
it is notable that a number of the trees in this study  showed intriguing topologies, depicting the split of archaeplastida and inclusion of chromalveolata and excavata genes within it, as shown in the previously reported multiple slowly-evolving gene phylogeny  <cit>  and gnd gene phylogeny  <cit> .
these results are consistent with the hypothesis for the non-monophyly of archaeplastida, and suggest that the oomycete and heterolobosean genes with cyanobacterial ancestry might reflect the host nuclear genome phylogeny.
on the other hand, the genes found in the marine choanoflagellate m. brevicollis were positioned within the bikonts clade, but not associated with the genes from other opisthokonta relatives , suggesting that the tree topologies were probably not reflective of the host phylogeny  <cit>  but eukaryote-to-eukaryote lgt .
no gene with cyanobacterial ancestry was found in d. discoideum , and only one gene in s. cerevisiae .
these results are also consistent with the ancient primary egt scenario.
a photorespiratory gene with cyanobacterial ancestry in fungi
our analysis using the genome data of the budding yeast s. cerevisiae identified one gene with cyanobacterial ancestry, encoding the glycerate kinase for photorespiration .
given that photorespiration is essential for cyanobacteria and plants, it is likely that the glycerate kinases in plants and cyanobacteria are phylogenetically and physiologically related to photorespiration  <cit> .
a previous study on glycerate kinases showed that, regardless of the complete absence of photorespiratory metabolism in fungi, the gene product from the budding yeast saccharomyces showed similar enzymatic activity and substrate specificity compared with the arabidopsis gene, suggesting that the plant and fungal genes catalyze the same reaction in different contexts of the metabolic pathway  <cit> .
another example of plant-type genes in fungi was reported in a phylogenetic study of the genes encoding high-affinity nitrate transporter nrt <dig>  which suggested that fungi probably acquired the nrt <dig> genes via lgt from one of the chromalveolate lineages  <cit> .
meanwhile, our data showed that the fungal clade was located outside the clade of plants plus oomycetes , suggesting that fungal glycerate kinase genes with cyanobacterial ancestry likely originated from an lgt event from an ancestor of cyanobacteria, or eukaryote-to-eukaryote lgt from an ancestor of archaeplastida .
one likely explanation for the presence of photorespiratory genes in oomycetes is that the ancestor of chromalveolata possessed this gene family, but some photosynthetic descendants lost this gene family or replaced it with other genes during the course of lineage-dependent customization of photorespiratory pathways , while oomycetes retained the genes without any replacement.
gene family expansion and differential reduction
another conclusion of this analysis is that rounds of gene family expansion and selective reduction are important factors in making eukaryotic genome phylogeny look like a complicated mosaic .
it is likely that the alteration of gene family repertoire contributed to the restructuring of the intracellular metabolome and a reduction of the dispensable gene families.
our data showed that all the genes identified in this study were members of multiple gene families.
algae and plastid-lacking protists retained only members of subfamilies , suggesting that the punctate distribution might be a corollary of the common mechanism by which genes with cyanobacterial ancestry were retained in their genomes.
the presence of genes from multiple subfamilies in one organism supports this idea .
discontinuous loss or gain of a metabolic pathway in a lineage might be another factor in punctate distribution; e.g.
the oxidative pentose pathway, and the cyanobacterial gnd genes functioning therein, were present in most bikonts but lost in the ciliate tetrahymena  <cit> .
a recent study on pyridoxal-dependent amino acid aminotransferase reported that, besides the ancestrally eukaryotic enzymes, land plants possess a distinct subfamily of prokaryote-type chloroplast-targeted enzymes  <cit> .
our data with richer taxon sampling identified another prokaryote-type subfamily with cyanobacterial ancestry , illustrating the hidden evolutionary diversity of protist and algal metabolomes.
future prospects
our results showed that many genes with cyanobacterial ancestry identified in this study were found only in complete genome sequences, suggesting that these genes might be difficult to discover by expressed sequence tag  library sequencing, probably due to the low-level expression of these genes.
although the whole genome data from excavate parasites  are available, they seem to be unsuited for the gene mining study because of the unusual nucleotide substitutions .
at the stage of starting the present gene mining study, n. gruberi was the only species with whole genome data released within the non-parasitic excavates, and thereby the excavate genes with cyanobacterial ancestry were mostly from n. gruberi.
more genome data from plastid-lacking protists from excavata and rhizaria as well as archaeplastida, especially red algae and glaucophytes, are needed to unravel the evolutionary history of plastids, and plastid-lacking protists.
CONCLUSIONS
the comparative analyses of the genome sequence data of the plastid-lacking eukaryotes demonstrated the potentially significant contributions of ancestral or extant cyanobacteria to the eukaryotic genomes, which probably occurred via lgt or ancient primary egt events.
furthermore, the automated phylogenetic analyses revealed the diversity and punctate distribution of gene families within the genomes in the unicellular microbes.
more genome data of the plastid-lacking excavata and rhizaria will make the evolutionary history clear and support our hypotheses.

 cyanobacterial ancestors gave rise to plastids in the ancestor of a eukaryotic lineage. the birth of the plastid had an impact on eukaryotic genome evolution, by way of endosymbiotic gene transfer, a particular form of lateral gene transfer from endosymbionts into the phylogenetically discontiguous host genome. subsequently, an algal ancestor gave rise to secondary plastids in several punctate lineages of eukaryotes. a number of these secondarily phototrophic lineages lost their photosynthetic ability and further diverged into secondarily heterotrophic, plastid-lacking protists. although the position of the root of eukaryotes is still uncertain, the presence of gene fusions and insertion/deletion sequences in the marker genes have allowed us to sort eukaryotes into at least three large groups; opisthokonta, amoebozoa and bikonts. most phototrophic eukaryotes harboring plastids derived from primary endosymbiosis are classified into the super-group archaeplastida. although it is widely accepted that primary plastids share a single origin and the archaeplastida are monophyletic, the evolutionary history of the primary plastids is still debatable. in plastid-lacking protists,'plastid imprints' can be exemplified by genomic information, i.e. genes with affinity to extant cyanobacterial or algal genes. these genes were supposed to have originated from egt events, and this assumption should be affirmed by the resulting phylogenetic relationship between'imprint' genes and the extant relatives of the putative endosymbionts. the biggest challenge and the limitation of this'imprint' searching process is that the inevitable incompleteness of genome information on lineages of interest and the ever-developing phylogenetic methodologies make it difficult to distinguish egt and ancient lgt. thus, although available eukaryotic genome data are increasingly accumulating, gene and genome phylogenies should be carefully interpreted to infer evolutionary scenarios. chromalveolata is a large taxonomic group of eukaryotes, encompassing secondary phototrophs and secondarily heterotrophic protists, and the'chromalveolate hypothesis' argues that this group originated from a common ancestor harboring the chlorophyll c-containing secondary plastid derived from a red alga. among the secondarily heterotrophic chromalveolates, several lineages have retained remnant chloroplasts for non-photosynthetic metabolic pathways, e.g. apicoplasts in apicomplexan parasites. recent genomic surveys revealed the presence of plastid-derived genes, and further suggested the presence of cryptic secondary plastids in non-photosynthetic alveolate protists. furthermore, re-examination of the whole genome sequences suggested the existence of algal genes in ciliates, another plastid-lacking alveolate lineage, which could support the photosynthetic ancestry of ciliates. oomycetes are plastid-lacking stramenopiles, or chromists, classified into chromalveolata. although whole genome sequence analysis showed that a number of genes with affinity to photosynthetic organisms are encoded in the nuclear genome, most of these'plastid imprints' candidates were only suggested by similarity search and phylogenetic analyses have not yet led to fully recovering the expected tree topology. considering the uncertain phylogenetic affinity of the'best hit' in similarity search, reassessment of the genome information is important to determine whether the evolutionary history of oomycetes is comparable to ciliates. one candidate of'plastid imprints' in oomycetes has been confirmed by studies reporting the phylogeny of gnd genes, which encode 6-phosphogluconate dehydrogenase, showing that some plastid-lacking protists have plant-like, cyanobacterium-derived gnd genes. these analyses suggested that the gnd genes with cyanobacterial ancestry were acquired early in eukaryotic evolution, either via ancient eukaryote-to-eukaryote lgt, or primary egt that occurred earlier than had ever been thought. additionally, the phylogeny of gnd genes demonstrated that cyanobacterial genes are also present in several excavata protists, e.g. the heterolobosean amoebo-flagellate naegleria gruberi. naegleria gruberi is a non-parasitic heterotrophic species related to n. fowleri, which is the causative agent of primary amoebic meningoencephalitis in mammals. although the phylogenetic relationship within excavata is still unclear, heterolobosea, together with jakobida, is likely to be a sister group of euglenozoa. to address how many genes have cyanobacterial ancestry in plastid-lacking protists, and whether cyanobacterial ancestry is limited to this gnd gene or also found in other genes, we conducted a phylogenomic analysis using genome sequence data of a taxonomically wide range of plastid-lacking eukaryotes. here we present a gene mining study with a novel pipeline automatically producing and summarizing one-by-one phylogenetic trees, and show phylogenetic analyses of resultant candidate genes with cyanobacterial ancestry, using the whole genome sequence data from a wide range of eukaryotic lineages.  to address how many genes are derived from cyanobacteria in non-photosynthetic protists, we conducted cyanobacterial gene mining using the genome sequence data of a wide range of the plastid-lacking eukaryotes. using the whole genome data, we conducted blast searches against all'bacteria' and selected queries showing the highest similarity to genes in the available cyanobacteria genome sequences. we then drew the neighbor-joining trees for genes showing homology to cyanobacterial counterparts. after the first tree construction step, we selected the gene trees where cyanobacteria and eukaryotes formed a monophyletic group excluding other prokaryotes. as a result, we obtained a shorter list of candidates, which we termed'genes with cyanobacterial affinity'. subsequently we re-analyzed the eukaryotic genes with cyanobacterial affinity by visually checking and re-drawing the bayesian and maximum likelihood trees after manually trimming operational taxonomic units. in total, we identified plastid-lacking protist genes'with cyanobacterial ancestry' in the genomes of the wide range of eukaryotes: two plastid-lacking bikonts and three unikonts. these were the eukaryotic genes with cyanobacterial ancestry that shared the same origin with archaeplastida and other eukaryotes. they were placed within a monophyletic subclade mostly composed of photosynthetic organisms and showed an apparent cyanobacterial ancestry as far as was determined by tree topology. we found another type of gene with cyanobacterial ancestry, which were the protist genes forming monophyletic groups mostly with genes from extant cyanobacteria. among a number of candidate genes found through the first screening, we have presented three typical trees that were resolved with significant support values. we postulate that these prokaryote-type genes are remnants of the bacterium-to-eukaryote lgt, which occurred'recently' in evolution. interestingly, while phytophthora ycf homologs, probably transferred from a relative of the extant cyanobacterial species via lgt, were placed within the cyanobacterial gene clade, the ciliate tetrahymena ycf homolog showed affinity to archaeplastida. this gene was neither found in the paramecium genome nor in the list of the recently identified algal genes in ciliates. jgi gene id, gene name and accession numbers are listed. we found that uroporphyrin iii methyltransferase gene homologs consisted of two large subfamilies of genes with cyanobacterial ancestry, and that oomycete genes were included only in one of them. given that both subfamilies include green plants, red algae, chromalveolates and cyanobacteria, it is likely that they diverged within the ancestral cyanobacteria and transferred into eukaryotic hosts via primary and secondary endosymbioses. both of the subfamilies were concurrently present in the cyanobacterial and green algal genomes. in land plants, red algae, diatoms, haptophytes and the plastid-lacking oomycetes, one of the subfamilies might be lost along with the loss of the plastid. the thalassiosira homolog formed a monophyletic group with green plants, rather than red algae, suggesting that it was acquired independently of the secondary plastid of the red lineage. in this study, the bacteriovorous choanoflagellate monosiga brevicollis gene and the proteobacterial genes were treated as'apparently lgt-derived genes', incongruously showing affinities to photosynthetic bikonts. genes encoding cobalamin-independent methionine synthase in green and red algae, diatoms, and oomycetes formed a monophyletic group with cyanobacterial homologs, while the land plants and the red alga cyanidioschyzon homologs were placed in different clades unrelated to cyanobacteria. close association between diatom and oomycete genes suggested the deep ancestry of the genes in the chromalveolate lineage. we failed to find the homologs in the prasinophytes ostreococcus and micromonas, suggesting that this gene family was dispensable in some plant lineages. one of the genes with cyanobacterial ancestry found in n. gruberi is pyridoxal-dependent amino acid decarboxylase gene. the tree indicated that green plants were split into different eukaryotic clades. naegleria and chromalveolate genes showed robust monophyly with green plants, included in a cyanobacterial gene clade. the tree showed that land plants possessed another subfamily, associated with red algal and fungal genes, apparently of non-cyanobacterial origin. we also identified genes with cyanobacterial ancestry from naegleria in an oxidoreductase gene family that included genes encoding rieske iron-sulfur cluster kda protein of chloroplast inner membrane translocon, chlorophyll a oxidase, lethal-leaf spot ) and accelerated cell death. all the members of this family in land plants were hypothesized to be located at the inner membrane of the chloroplast, and to be involved in chlorophyll metabolism. the phylogenetic tree of the tic55-like gene family showed intricate distribution of cyanobacterial, green plant and chromalveolate genes. in other trees of the genes identified in this study, gene clades with cyanobacterial ancestry were mostly composed of bikonts genes, besides the choanoflagellate m. brevicollis genes. discussion we identified eight and seven genes with cyanobacterial ancestry in the genome sequences of the oomycete p. ramorum and the heterolobosean n. gruberi, respectively. it was reported that the apicomplexan cryptosporidium'recently' lost their secondary plastid, and retained two to seven putative plastid-derived genes in the genome. this number is comparable to our result of the gene mining study using oomycete and heterolobosean genomes. in addition, our system resolved the hidden diversity of the gene family repertoire in eukaryotic genomes by one-by-one gene phylogenies. secondary egt scenario although the phylogenetic positions of cryptophyceae and haptophyta are still debatable, the chromalveolate hypothesis has been reinstated to support the evolutionary scenario that the plastid-lacking protists oomycetes and ciliates once might have had a plastid. according to this hypothesis, the genes with cyanobacterial ancestry found in the oomycete genomes were acquired via secondary egt in the common ancestor of chromalveolata, from the red algal ancestor of secondary plastids. this explanation is also applicable under the alternative hypothesis for chromalveolate plastids, which proposes that a tertiary endosymbiont of the haptophyte/cryptophyceae lineage is the origin of the stramenopile/alveolate plastids. the phylogenetic tree of the photorespiratory glycerate kinase genes, suggesting the red algal origin of the phytophthora genes, is consistent with the chromalveolate hypothesis. however, several other gene trees in this study showed oomycete genes with green lineage affinity, not red algae. recently, frommolt et al. demonstrated that, out of genes involved in carotenoid biosynthesis from chromalveolate algae, one third of plastid-targeted, nuclear-encoded genes are most closely related to green algal homologs. reyes-prieto, moustafa and bhattacharya identified genes of possible algal origin in the ciliates tetrahymena thermophila and paramecium tetraurelia, and 7/ of their trees show a close relationship between green plants and chromalveolata. frommolt et al. attributed the close relationships between green plants and chromalveolate genes to the secondary endosymbiosis of an ancient green plant, based on the hypothesis on the monophyly of the archaeplastida. this explanation might be also applicable to the plant-like genes in ciliates. while heterolobosea and euglenozoa are often united as the morphologically defined taxon, discicristata, within excavata, recent morphological and molecular phylogenetic analyses suggest that the heteroloboseans never possessed the secondary plastid of green lineage and share the same origin with euglenida. molecular phylogenetic analyses showed that excavata is separated from other secondary plastid-containing eukaryotes. therefore, it is unlikely that the genes with cyanobacterial ancestry found in the heterolobosean nuclear genomes originated from the plastid cognate with any known secondary plastids in extant photosynthetic eukaryotes. the amino acid decarboxylase gene and the gnd gene trees demonstrated the presence of genes with cyanobacterial ancestry in other heterolobosean species than n. gruberi, suggesting that the ancestor of the genus naegleria possessed this gene family. furthermore, although ml bootstrap support or bayesian posterior probability values were not always sufficient, the naegleria genes occupy relatively basal phylogenetic positions within the bikonts clade in all seven trees. thus it is possible that the genes with cyanobacterial ancestry were introduced en bloc in the ancestor of heterolobosea, via a batch gene transfer, in a concerted manner. one possible origin of such a concerted gene transfer is secondary egt from a photosynthetic eukaryote with a basal phylogenetic position within bikonts. however, as discussed above, it is unlikely that heterolobosea experienced secondary endosymbiosis and acquired genes common to the extant secondary plastid-containing eukaryotes via secondary egt. ancient eukaryote-to-eukaryote lgt or primary egt scenarios alternatively, we can argue for two other explanations: a concerted eukaryote-to-eukaryote lgt scenario or a more ancient primary egt scenario. the naegleria genes with cyanobacterial ancestry shown in table are basally positioned within bikonts, but not intruding into any of gene clades from extant photosynthetic eukaryotes. thus, if we assume that these genes were acquired via non-endosymbiotic lgt, they may originate from unknown ancient photosynthetic lineages basally positioned within bikonts. meanwhile, under the primary egt scenario, in which the primary endosymbiosis occurred in the common ancestor of bikonts, ancient primary egt occurred much earlier than the conventional hypothesis, from the cyanobacterium-like prokaryote to the common ancestor of bikonts. primary plastids were subsequently lost in many lineages of bikonts, except for the archaeplastida lineages, but some genes originating from the cyanobacterial ancestor of the primary plastids have been retained in the nuclear genomes of the plastid-lacking lineages of bikonts. the loss of the plastid might have triggered the loss of genes that specifically functioned within the plastid. only a portion of the plastid-derived genes, which we can find now in the plastid-lacking protist genomes, might have escaped from or survived through eliminative pressure in a lineage-specific manner, by acquiring additional functions with other components and/or in other cellular compartments. this might account for the observed punctate distribution of gene families among the eukaryotes. recently, a hypothesis for the non-monophyly of archaeplastida was proposed based on the phylogenetic analyses of slowly evolving nuclear-encoded genes. this non-monophyly hypothesis could be also considered within the scope of the primary egt scenario. it is notable that a number of the trees in this study showed intriguing topologies, depicting the split of archaeplastida and inclusion of chromalveolata and excavata genes within it, as shown in the previously reported multiple slowly-evolving gene phylogeny and gnd gene phylogeny. these results are consistent with the hypothesis for the non-monophyly of archaeplastida, and suggest that the oomycete and heterolobosean genes with cyanobacterial ancestry might reflect the host nuclear genome phylogeny. on the other hand, the genes found in the marine choanoflagellate m. brevicollis were positioned within the bikonts clade, but not associated with the genes from other opisthokonta relatives, suggesting that the tree topologies were probably not reflective of the host phylogeny but eukaryote-to-eukaryote lgt. no gene with cyanobacterial ancestry was found in d. discoideum, and only one gene in s. cerevisiae. these results are also consistent with the ancient primary egt scenario. a photorespiratory gene with cyanobacterial ancestry in fungi our analysis using the genome data of the budding yeast s. cerevisiae identified one gene with cyanobacterial ancestry, encoding the glycerate kinase for photorespiration. given that photorespiration is essential for cyanobacteria and plants, it is likely that the glycerate kinases in plants and cyanobacteria are phylogenetically and physiologically related to photorespiration. a previous study on glycerate kinases showed that, regardless of the complete absence of photorespiratory metabolism in fungi, the gene product from the budding yeast saccharomyces showed similar enzymatic activity and substrate specificity compared with the arabidopsis gene, suggesting that the plant and fungal genes catalyze the same reaction in different contexts of the metabolic pathway. another example of plant-type genes in fungi was reported in a phylogenetic study of the genes encoding high-affinity nitrate transporter nrt which suggested that fungi probably acquired the nrt genes via lgt from one of the chromalveolate lineages. meanwhile, our data showed that the fungal clade was located outside the clade of plants plus oomycetes, suggesting that fungal glycerate kinase genes with cyanobacterial ancestry likely originated from an lgt event from an ancestor of cyanobacteria, or eukaryote-to-eukaryote lgt from an ancestor of archaeplastida. one likely explanation for the presence of photorespiratory genes in oomycetes is that the ancestor of chromalveolata possessed this gene family, but some photosynthetic descendants lost this gene family or replaced it with other genes during the course of lineage-dependent customization of photorespiratory pathways, while oomycetes retained the genes without any replacement. gene family expansion and differential reduction another conclusion of this analysis is that rounds of gene family expansion and selective reduction are important factors in making eukaryotic genome phylogeny look like a complicated mosaic. it is likely that the alteration of gene family repertoire contributed to the restructuring of the intracellular metabolome and a reduction of the dispensable gene families. our data showed that all the genes identified in this study were members of multiple gene families. algae and plastid-lacking protists retained only members of subfamilies, suggesting that the punctate distribution might be a corollary of the common mechanism by which genes with cyanobacterial ancestry were retained in their genomes. the presence of genes from multiple subfamilies in one organism supports this idea. discontinuous loss or gain of a metabolic pathway in a lineage might be another factor in punctate distribution; e.g. the oxidative pentose pathway, and the cyanobacterial gnd genes functioning therein, were present in most bikonts but lost in the ciliate tetrahymena. a recent study on pyridoxal-dependent amino acid aminotransferase reported that, besides the ancestrally eukaryotic enzymes, land plants possess a distinct subfamily of prokaryote-type chloroplast-targeted enzymes. our data with richer taxon sampling identified another prokaryote-type subfamily with cyanobacterial ancestry, illustrating the hidden evolutionary diversity of protist and algal metabolomes. future prospects our results showed that many genes with cyanobacterial ancestry identified in this study were found only in complete genome sequences, suggesting that these genes might be difficult to discover by expressed sequence tag library sequencing, probably due to the low-level expression of these genes. although the whole genome data from excavate parasites are available, they seem to be unsuited for the gene mining study because of the unusual nucleotide substitutions. at the stage of starting the present gene mining study, n. gruberi was the only species with whole genome data released within the non-parasitic excavates, and thereby the excavate genes with cyanobacterial ancestry were mostly from n. gruberi. more genome data from plastid-lacking protists from excavata and rhizaria as well as archaeplastida, especially red algae and glaucophytes, are needed to unravel the evolutionary history of plastids, and plastid-lacking protists.  the comparative analyses of the genome sequence data of the plastid-lacking eukaryotes demonstrated the potentially significant contributions of ancestral or extant cyanobacteria to the eukaryotic genomes, which probably occurred via lgt or ancient primary egt events. furthermore, the automated phylogenetic analyses revealed the diversity and punctate distribution of gene families within the genomes in the unicellular microbes. more genome data of the plastid-lacking excavata and rhizaria will make the evolutionary history clear and support our hypotheses. 
eukaryotic genes with cyanobacterial ancestry in plastid-lacking protists have been regarded as important evolutionary markers implicating the presence of plastids in the early evolution of eukaryotes. our data also show that all genes identified in this study constitute multi-gene families with punctate distribution among eukaryotes, suggesting that the transferred genes could have survived through rounds of gene family expansion and differential reduction. the origins of those genes might be due to lateral gene transfer events, or an ancient primary or secondary endosymbiosis before the diversification of bikonts. although recent genomic surveys demonstrated the presence of cyanobacterial and algal ancestry genes in the genomes of plastid-lacking protists, comparative analyses on the origin and distribution of those genes are still limited. we identified gene families with cyanobacterial ancestry in the genomes of a taxonomically wide range of plastid-lacking eukaryotes using a novel phylogenetic pipeline. the eukaryotic gene clades with cyanobacterial ancestry were mostly composed of genes from bikonts. we failed to find genes with cyanobacterial ancestry in saccharomyces and dictyostelium, except for a photorespiratory enzyme conserved among fungi. meanwhile, we found several monosiga genes with cyanobacterial ancestry, which were unrelated to other opisthokonta genes. our data demonstrate that a considerable number of genes with cyanobacterial ancestry have contributed to the genome composition of the plastid-lacking protists, especially bikonts.
 eukaryotic genes with cyanobacterial ancestry in plastid-lacking protists have been regarded as important evolutionary markers implicating the presence of plastids in the early evolution of eukaryotes. although recent genomic surveys demonstrated the presence of cyanobacterial and algal ancestry genes in the genomes of plastid-lacking protists, comparative analyses on the origin and distribution of those genes are still limited.  we identified gene families with cyanobacterial ancestry in the genomes of a taxonomically wide range of plastid-lacking eukaryotes using a novel phylogenetic pipeline. the eukaryotic gene clades with cyanobacterial ancestry were mostly composed of genes from bikonts. we failed to find genes with cyanobacterial ancestry in saccharomyces and dictyostelium, except for a photorespiratory enzyme conserved among fungi. meanwhile, we found several monosiga genes with cyanobacterial ancestry, which were unrelated to other opisthokonta genes.  our data demonstrate that a considerable number of genes with cyanobacterial ancestry have contributed to the genome composition of the plastid-lacking protists, especially bikonts. the origins of those genes might be due to lateral gene transfer events, or an ancient primary or secondary endosymbiosis before the diversification of bikonts. our data also show that all genes identified in this study constitute multi-gene families with punctate distribution among eukaryotes, suggesting that the transferred genes could have survived through rounds of gene family expansion and differential reduction. 