BACKGROUND
phylogenomics
recent advances in sequencing technologies and reductions in the cost of sequencing have fostered an unprecedented explosion of complete genome and expressed sequence tag  data and facilitated a revolution in the field of comparative genomics.
the availability of this massive amount of genome data has been both a boon and a major challenge for biologists.
phylogenomics offers an useful avenue for dealing with genome data, allowing researchers to investigate gene phylogeny in a genome-wide context  <cit> .
this approach can be used to annotate genes of unknown function in newly sequenced genomes or to identify phylogenetic markers to infer an accurate genealogy of life, or to understand the role of endosymbiotic and horizontal gene transfer in eukaryotic evolution  <cit> .
existing phylogenomic pipelines  generate a large collection, often hundreds or thousands, of phylogenetic trees that require manual inspection to observe general patterns of genome evolution or to address specific hypotheses about gene phylogeny.
endosymbiosis
the engulfment of a free-living photosynthetic cyanobacterium by a heterotrophic protist  introduced photosynthesis into the eukaryotic domain.
the primary endosymbiosis occurred about  <dig>  billion years ago   <cit>  and was a turning point for evolution of life on our planet allowing the later development of multicellular plants and animals.
the first algae diversified over time into the three primary photosynthetic lineages, the green , red, and glaucophyte algae.
these taxa are united in the eukaryotic supergroup plantae.
some time after establishment of the endosymbiotic relationship there was selective pressure to reduce the endosymbiont genome by outright gene loss or transfer of genes to the host nuclear genome.
this latter process is termed endosymbiotic gene transfer .
determining the extent of 'primary' egt in plantae and whether only genes involved in plastid function were retained or were augmented by many other genes of non-organellar function remain controversial issues.
previously, using the flowering plant arabidopsis, it was estimated that the cyanobacterial endosymbiont contributed 18% of the total set of nuclear genes in this species  <cit> .
in another study, egt was shown to have contributed only about 4% of the nuclear genes in arabidopsis, and 12% to the reduced nuclear genome of the extremophilic red alga cyanidioschyzon merolae  <cit> .
in a more recent study using the free-living glaucophyte alga cyanophora paradoxa, we estimated that about 4% of the nuclear genome in this taxon was of cyanobacterial provenance.
however this latter study was based on incomplete est data and is a provisional result.
there exists therefore a need to apply modern methods to analyze the complete genome sequence of a mesophilic, free-living, unicellular alga to generate a robust estimate of primary egt in a relatively 'simple' ancestor of land plants  <cit> .
furthermore, once the plantae split into its constituent lineages, a red  alga was captured by the ancestor of the chromalveolate protists via secondary endosymbiosis  <cit> .
this process necessitated the movement of genes involved in plastid function from the nucleus of the algal endosymbiont to that of the chromalveolate host through 'secondary' egt.
the chromalveolate supergroup includes a broad swath of protist diversity including both photosynthetic  and plastid-lacking  lineages.
presumably the latter group have lost their plastid secondarily.
phylogenomic analysis of nuclear genes in these algae is expected therefore to return trees which show gene origin in plantae from a cyanobacterial ancestor, followed by their transfer to chromalveolates via a red or green algal secondary endosymbiosis  <cit>  and  <cit> .
given the need to better understand egt in algae using a phylogenomic approach, we developed phylosort to analyze topologies in a high-throughput fashion.
this opensource java tool inspects the topology of phylogenetic trees to address the most frequently asked question in the field: does a specific gene support the monophyly of certain operational taxonomic units ?
here, we provide an overview of phylosort and its application to the complete set of predicted proteins in the green alga chlamydomonas reinhardtii to estimate the contribution through egt of cyanobacteria to the nuclear genome of a green alga.
implementation
phylosort can be used via a graphical user interface  and a text mode command line interface.
input phylogenetic trees are read and parsed from a source folder, where trees are stored as one tree per file in newick format, which is supported and produced by many phylogenetic inference tools such as phylip  <cit> , paup*  <cit> , phyml  <cit>  and raxml  <cit> .
of the input trees, those that satisfy the search criteria are copied or moved to an output folder.
to begin the analysis, the hypothesized monophyletic taxa are selected from a pool of taxa.
this pool of taxa can be loaded as a simple list from a plain text file.
alternatively, a tree can be loaded that acts as a taxonomy reference for organizing the taxa in a phylogenetic format .
finally, if no list is loaded or no reference tree exists, the program will unite all taxa in all trees into a single non-redundant list of taxa.
in a typical phylogenomic analysis, homologs within and among different genomes give rise to multiple sequence alignments  that do not all necessarily share among them the same set of taxa.
therefore, the number of taxa that represent members of any monophyletic group  may  vary from tree to tree.
this can be explained, for example, by lineage-specific gene duplications or gene losses.
for this reason, in phylosort, the hypothesized taxa are arranged in groups such that for any number of groups to be monophyletic each group has to be represented by at least one of its constituent taxa.
in addition, because the names of the sequences used in reconstructing the trees would have different formats from one project or research group to another, "regular expressions"  <cit>  are used to extract taxa  from the sequence names.
in addition to the topological constraint of monophyly, the phylosort search can be adjusted by setting a minimum bootstrap support value associated with the monophyletic clades.
this allows the identification of trees with significantly supported  monophyletic relationships.
minimum and maximum number of taxa in a trees and average number of genes per taxon can also be chosen to sample differing levels of gene family complexity.
phylosort offers two search modes, exclusive and inclusive.
in the exclusive mode, the taxa that exist in the tree and belong to the hypothesized monophyletic taxa have to be located in a single monophyletic clade.
in the inclusive mode, any group of taxa forming a clade that matches the search criteria qualifies the trees regardless of whether other taxa belong to the hypothesized monophyletic taxa exist elsewhere in the tree.
prior to the monophyly search, each tree is searched for an outgroup and then rerooted on that taxon.
a genome-wide analysis generally produces a significant number of trees that share multiple genes due to multiple gene copies and gene families.
accordingly, to provide an estimate of the number of unique gene families, phylosort has a simple clustering functionality that can be used to merge trees by identifying overlapping genes among trees and placing the trees into 'tree clusters' representing gene families.
a minimum number of overlapping genes can be set for merging trees into clusters.
through the gui interface, a taxonomy reference tree  can be used to hierarchically collect taxa and simplify the assignment of taxa into groups, and phylogenetic trees can be visually inspected using atv  <cit> .
in addition to the gui and command line interfaces, phylosort provides a set of reusable and extendable application programming interfaces  that can be incorporated into other applications that may utilize the monophyly search or other utility components such as newick format parsing and basic phylogenetic tree manipulation.
to determine whether a set of taxa is monophyletic, there are two main steps.
first, the lowest common ancestor  is located.
second, the subtree rooted by the lca is verified to not contain outgroups.
the following subsections summarize the approaches that have been implemented to determine whether a tree contains a clade that matches the monophyletic criterion and any additional constraints in both search modes.
lca identification
to locate the lca for a set of leaves , the path from the root to each leaf in the set is determined.
next, all paths are compared to find the longest shared segment .
then, the lca is the furthest node from the root on the longest shared segment .
exclusive search
under the exclusive mode, all taxa in the tree that belong to the set of taxa being examined for monophyly are located and the lca for the entire set is determined and the clade rooted by the lca is tested for monophyly.
inclusive search
under the inclusive mode, each leaf that belongs to the groups of taxa being examined for monophyly is used as a starting point to traverse the tree in a tip to root direction.
each internal node is examined as described above in the exclusive search to test whether it contains a qualifying monophyletic clade.
if the clade rooted by the internal node has a taxon that does not belong to the hypothesized monophyletic taxa, the clade is rejected, all nodes belong to the clade are marked not to be revisited, and the next leaf is examined.
if a clade is not rejected but it does not satisfy the monophyly constraint or any additional filtering criterion, the parent of the clade is located and the clade rooted by the parent is examined.
otherwise, if the clade satisfies all criteria, the search stops returning the tree as a matching tree.
this search strategy was specifically designed to identify trees that contain paralogs, a subset of which may satisfy the monophyly criterion and therefore should be considered.
RESULTS
as an example of the application of phylosort to address biological questions, we quantified primary  egt in the green alga chlamydomonas reinhardtii.
of the  <dig>  predicted chlamydomonas nuclear-encoded proteins, we found  <dig>  that had cyanobacterial homologs with a blast e-value < 1e- <dig>  for which  <dig>  trees were inferred.
with no restriction on the bootstrap support value , there were  <dig>  trees that showed monophyly of cyanobacteria and plantae .
by enforcing the minimum bootstrap support values at ≥ 50% or ≥ 75%, we found  <dig>  and  <dig> trees, respectively, that satisfied the monophyly constraint.
clustering the 50% bootstrap category resulted into at least  <dig> unique gene families.
based on the gene ontology analysis , at least 44% of these predicted cyanobacterial genes were identified as encoding plastid-targeted proteins and at least 47% were involved in metabolic processes.
discussion
of the trees that supported the monophyly of cyanobacteria and photosynthetic eukaryotes , we present here as an example the phylogenetic tree of the thylakoid lumenal  <dig>  kda protein .
this plastid-targeted protein belongs to the "pentapeptide repeat" family of proteins of which the function has not yet been characterized in plantae.
however, it has been shown that the pentapeptide repeats are required for a proper localization of heterocyst glycolipids in cyanobacteria  <cit> .
the tree demonstrates the phylogeny predicted for plastid primary and secondary endosymbiosis  <cit> .
the nuclear-encoded plantae proteins are of cyanobacterial origin , whereas the chromalveolate sequences are rooted within red algae .
here the plantae are represented by its three member lineages, the glaucophyte cyanophora, the red alga cyanidioschyzon, the green algae ostreococcus, chlamydomonas, volvox and physcomitrella, and the plants arabidopsis and oryza, whereas the chromalveolates are represented by the diatoms phaeodactylum and thalassiosira, the haptophyte emiliania, the cryptophyte guillardia, and the dinoflagellate alexandrium).
previous analyses suggest that the chromalveolate red algal secondary endosymbiosis occurred 'soon'  after the cyanobacterial capture.
many other chlamydomonas proteins display the same topology as shown in figure  <dig> and are available as individual trees in .
specifically, of the  <dig> genes that belong to the 50% bootstrap category, about 51% of the genes were shared among red and green algae and about 12% were shared among glaucophytes, red and green algae.
the unevenness of the distribution of cyanobacterial genes among the plantae lineages can be explained by limited est data that are available from glaucophytes , the highly reduced nature of the genome of the red alga cyanidioschyzon  which likely precipitated the loss of many target genes, and the sampling bias towards the green lineages due to the use of chlamydomonas as the query for the phylogenomic analysis.
in summary, using phylosort we were able to rapidly inspect thousands of phylogenetic trees of differing complexities and number of terminal taxa and demonstrate the contribution of cyanobacterial primary egt to plantae nuclear genomes.
CONCLUSIONS
phylosort is a simple, platform independent and user-friendly tool to inspect the topology and other aspects  of large sets of phylogenetic trees in a phylogenomic context.
phylosort can be used to quantify the number of genes  that support different evolutionary hypotheses such as a taxonomic classification or to quantify endosymbiotic or horizontal gene transfer events.
in addition, phylosort can be used to search for vertically inherited phylogenetic markers, for example across the eukaryotic tree of life, by testing the monophyly in individual gene trees of members of each of the six putative eukaryotic supergroups  <cit> .
our analysis of chlamydomonas  demonstrates that cyanobacteria account for ca.
<dig> –6%  of the protein-coding genes in the nuclear genome of this taxon and that at least 63% of the transferred cyanobacterial genes are present in the red and green lineages.
a previous analysis suggested that about 90% of the gene products resulting from primary egt are targeted to the plastid  <cit> .
therefore it appears that whereas primary endosymbiosis strongly favored egt to the nucleus, selection reduced the impact of cyanobacterial gene incursions primarily to functions that are plastid localized.
furthermore, it is now clear that outright loss was the most common fate for cyanobacterial genes in algae because with ca.
<dig>  genes in the nucleus and another 150– <dig> in the plastid, this accounts for less than 50%  of the ancestral endosymbiont genome.
finally, we discuss briefly the relevance of our findings to the extent of 'real' primary egt.
our results are based on detecting phylogenetic signal in anciently diverged sequences using standard computational approaches.
given that phylogenetic signal was most certainly lost from many cyanobacterial genes over the greater than  <dig> billion years since the endosymbiotic event, it will most certainly be impossible to detect all instances of egt.
therefore, how significant is our number of  <dig> cyanobacterial-derived genes in chlamydomonas ?
in our opinion, this must clearly be viewed as an underestimate of the actual number.
many cyanobacterial genes may simply no longer be detectable as such due to highly diverged sequences.
a more accurate estimate may come from increased taxon sampling  in the phylogenomic analyses to allow these methods to potentially identify a larger proportion of highly diverged homologs among the plantae.
we also need to gain a better understanding of the process of phylogenetic signal decay among this gene set so that extrapolations to the true value are more realistic.
availability and requirements
phylosort is available as an open-source under the gnu general public license  via www at http://phylosort.sourceforge.net.
phylosort is implemented in java under compliance with compiler version  <dig> .
it has been successfully tested on machines running the operating systems sun os, red hat, windows, and mac os.

 phylogenomics recent advances in sequencing technologies and reductions in the cost of sequencing have fostered an unprecedented explosion of complete genome and expressed sequence tag data and facilitated a revolution in the field of comparative genomics. the availability of this massive amount of genome data has been both a boon and a major challenge for biologists. phylogenomics offers an useful avenue for dealing with genome data, allowing researchers to investigate gene phylogeny in a genome-wide context. this approach can be used to annotate genes of unknown function in newly sequenced genomes or to identify phylogenetic markers to infer an accurate genealogy of life, or to understand the role of endosymbiotic and horizontal gene transfer in eukaryotic evolution. existing phylogenomic pipelines generate a large collection, often hundreds or thousands, of phylogenetic trees that require manual inspection to observe general patterns of genome evolution or to address specific hypotheses about gene phylogeny. endosymbiosis the engulfment of a free-living photosynthetic cyanobacterium by a heterotrophic protist introduced photosynthesis into the eukaryotic domain. the primary endosymbiosis occurred about billion years ago and was a turning point for evolution of life on our planet allowing the later development of multicellular plants and animals. the first algae diversified over time into the three primary photosynthetic lineages, the green, red, and glaucophyte algae. these taxa are united in the eukaryotic supergroup plantae. some time after establishment of the endosymbiotic relationship there was selective pressure to reduce the endosymbiont genome by outright gene loss or transfer of genes to the host nuclear genome. this latter process is termed endosymbiotic gene transfer. determining the extent of'primary' egt in plantae and whether only genes involved in plastid function were retained or were augmented by many other genes of non-organellar function remain controversial issues. previously, using the flowering plant arabidopsis, it was estimated that the cyanobacterial endosymbiont contributed 18% of the total set of nuclear genes in this species. in another study, egt was shown to have contributed only about 4% of the nuclear genes in arabidopsis, and 12% to the reduced nuclear genome of the extremophilic red alga cyanidioschyzon merolae. in a more recent study using the free-living glaucophyte alga cyanophora paradoxa, we estimated that about 4% of the nuclear genome in this taxon was of cyanobacterial provenance. however this latter study was based on incomplete est data and is a provisional result. there exists therefore a need to apply modern methods to analyze the complete genome sequence of a mesophilic, free-living, unicellular alga to generate a robust estimate of primary egt in a relatively'simple' ancestor of land plants. furthermore, once the plantae split into its constituent lineages, a red alga was captured by the ancestor of the chromalveolate protists via secondary endosymbiosis. this process necessitated the movement of genes involved in plastid function from the nucleus of the algal endosymbiont to that of the chromalveolate host through'secondary' egt. the chromalveolate supergroup includes a broad swath of protist diversity including both photosynthetic and plastid-lacking lineages. presumably the latter group have lost their plastid secondarily. phylogenomic analysis of nuclear genes in these algae is expected therefore to return trees which show gene origin in plantae from a cyanobacterial ancestor, followed by their transfer to chromalveolates via a red or green algal secondary endosymbiosis and. given the need to better understand egt in algae using a phylogenomic approach, we developed phylosort to analyze topologies in a high-throughput fashion. this opensource java tool inspects the topology of phylogenetic trees to address the most frequently asked question in the field: does a specific gene support the monophyly of certain operational taxonomic units? here, we provide an overview of phylosort and its application to the complete set of predicted proteins in the green alga chlamydomonas reinhardtii to estimate the contribution through egt of cyanobacteria to the nuclear genome of a green alga. implementation phylosort can be used via a graphical user interface and a text mode command line interface. input phylogenetic trees are read and parsed from a source folder, where trees are stored as one tree per file in newick format, which is supported and produced by many phylogenetic inference tools such as phylip, paup*, phyml and raxml. of the input trees, those that satisfy the search criteria are copied or moved to an output folder. to begin the analysis, the hypothesized monophyletic taxa are selected from a pool of taxa. this pool of taxa can be loaded as a simple list from a plain text file. alternatively, a tree can be loaded that acts as a taxonomy reference for organizing the taxa in a phylogenetic format. finally, if no list is loaded or no reference tree exists, the program will unite all taxa in all trees into a single non-redundant list of taxa. in a typical phylogenomic analysis, homologs within and among different genomes give rise to multiple sequence alignments that do not all necessarily share among them the same set of taxa. therefore, the number of taxa that represent members of any monophyletic group may vary from tree to tree. this can be explained, for example, by lineage-specific gene duplications or gene losses. for this reason, in phylosort, the hypothesized taxa are arranged in groups such that for any number of groups to be monophyletic each group has to be represented by at least one of its constituent taxa. in addition, because the names of the sequences used in reconstructing the trees would have different formats from one project or research group to another, "regular expressions" are used to extract taxa from the sequence names. in addition to the topological constraint of monophyly, the phylosort search can be adjusted by setting a minimum bootstrap support value associated with the monophyletic clades. this allows the identification of trees with significantly supported monophyletic relationships. minimum and maximum number of taxa in a trees and average number of genes per taxon can also be chosen to sample differing levels of gene family complexity. phylosort offers two search modes, exclusive and inclusive. in the exclusive mode, the taxa that exist in the tree and belong to the hypothesized monophyletic taxa have to be located in a single monophyletic clade. in the inclusive mode, any group of taxa forming a clade that matches the search criteria qualifies the trees regardless of whether other taxa belong to the hypothesized monophyletic taxa exist elsewhere in the tree. prior to the monophyly search, each tree is searched for an outgroup and then rerooted on that taxon. a genome-wide analysis generally produces a significant number of trees that share multiple genes due to multiple gene copies and gene families. accordingly, to provide an estimate of the number of unique gene families, phylosort has a simple clustering functionality that can be used to merge trees by identifying overlapping genes among trees and placing the trees into'tree clusters' representing gene families. a minimum number of overlapping genes can be set for merging trees into clusters. through the gui interface, a taxonomy reference tree can be used to hierarchically collect taxa and simplify the assignment of taxa into groups, and phylogenetic trees can be visually inspected using atv. in addition to the gui and command line interfaces, phylosort provides a set of reusable and extendable application programming interfaces that can be incorporated into other applications that may utilize the monophyly search or other utility components such as newick format parsing and basic phylogenetic tree manipulation. to determine whether a set of taxa is monophyletic, there are two main steps. first, the lowest common ancestor is located. second, the subtree rooted by the lca is verified to not contain outgroups. the following subsections summarize the approaches that have been implemented to determine whether a tree contains a clade that matches the monophyletic criterion and any additional constraints in both search modes. lca identification to locate the lca for a set of leaves, the path from the root to each leaf in the set is determined. next, all paths are compared to find the longest shared segment. then, the lca is the furthest node from the root on the longest shared segment. exclusive search under the exclusive mode, all taxa in the tree that belong to the set of taxa being examined for monophyly are located and the lca for the entire set is determined and the clade rooted by the lca is tested for monophyly. inclusive search under the inclusive mode, each leaf that belongs to the groups of taxa being examined for monophyly is used as a starting point to traverse the tree in a tip to root direction. each internal node is examined as described above in the exclusive search to test whether it contains a qualifying monophyletic clade. if the clade rooted by the internal node has a taxon that does not belong to the hypothesized monophyletic taxa, the clade is rejected, all nodes belong to the clade are marked not to be revisited, and the next leaf is examined. if a clade is not rejected but it does not satisfy the monophyly constraint or any additional filtering criterion, the parent of the clade is located and the clade rooted by the parent is examined. otherwise, if the clade satisfies all criteria, the search stops returning the tree as a matching tree. this search strategy was specifically designed to identify trees that contain paralogs, a subset of which may satisfy the monophyly criterion and therefore should be considered.  as an example of the application of phylosort to address biological questions, we quantified primary egt in the green alga chlamydomonas reinhardtii. of the predicted chlamydomonas nuclear-encoded proteins, we found that had cyanobacterial homologs with a blast e-value < 1e- for which trees were inferred. with no restriction on the bootstrap support value, there were trees that showed monophyly of cyanobacteria and plantae. by enforcing the minimum bootstrap support values at ≥ 50% or ≥ 75%, we found and trees, respectively, that satisfied the monophyly constraint. clustering the 50% bootstrap category resulted into at least unique gene families. based on the gene ontology analysis, at least 44% of these predicted cyanobacterial genes were identified as encoding plastid-targeted proteins and at least 47% were involved in metabolic processes. discussion of the trees that supported the monophyly of cyanobacteria and photosynthetic eukaryotes, we present here as an example the phylogenetic tree of the thylakoid lumenal kda protein. this plastid-targeted protein belongs to the "pentapeptide repeat" family of proteins of which the function has not yet been characterized in plantae. however, it has been shown that the pentapeptide repeats are required for a proper localization of heterocyst glycolipids in cyanobacteria. the tree demonstrates the phylogeny predicted for plastid primary and secondary endosymbiosis. the nuclear-encoded plantae proteins are of cyanobacterial origin, whereas the chromalveolate sequences are rooted within red algae. here the plantae are represented by its three member lineages, the glaucophyte cyanophora, the red alga cyanidioschyzon, the green algae ostreococcus, chlamydomonas, volvox and physcomitrella, and the plants arabidopsis and oryza, whereas the chromalveolates are represented by the diatoms phaeodactylum and thalassiosira, the haptophyte emiliania, the cryptophyte guillardia, and the dinoflagellate alexandrium). previous analyses suggest that the chromalveolate red algal secondary endosymbiosis occurred'soon' after the cyanobacterial capture. many other chlamydomonas proteins display the same topology as shown in figure and are available as individual trees in. specifically, of the genes that belong to the 50% bootstrap category, about 51% of the genes were shared among red and green algae and about 12% were shared among glaucophytes, red and green algae. the unevenness of the distribution of cyanobacterial genes among the plantae lineages can be explained by limited est data that are available from glaucophytes, the highly reduced nature of the genome of the red alga cyanidioschyzon which likely precipitated the loss of many target genes, and the sampling bias towards the green lineages due to the use of chlamydomonas as the query for the phylogenomic analysis. in summary, using phylosort we were able to rapidly inspect thousands of phylogenetic trees of differing complexities and number of terminal taxa and demonstrate the contribution of cyanobacterial primary egt to plantae nuclear genomes.  phylosort is a simple, platform independent and user-friendly tool to inspect the topology and other aspects of large sets of phylogenetic trees in a phylogenomic context. phylosort can be used to quantify the number of genes that support different evolutionary hypotheses such as a taxonomic classification or to quantify endosymbiotic or horizontal gene transfer events. in addition, phylosort can be used to search for vertically inherited phylogenetic markers, for example across the eukaryotic tree of life, by testing the monophyly in individual gene trees of members of each of the six putative eukaryotic supergroups. our analysis of chlamydomonas demonstrates that cyanobacteria account for ca. –6% of the protein-coding genes in the nuclear genome of this taxon and that at least 63% of the transferred cyanobacterial genes are present in the red and green lineages. a previous analysis suggested that about 90% of the gene products resulting from primary egt are targeted to the plastid. therefore it appears that whereas primary endosymbiosis strongly favored egt to the nucleus, selection reduced the impact of cyanobacterial gene incursions primarily to functions that are plastid localized. furthermore, it is now clear that outright loss was the most common fate for cyanobacterial genes in algae because with ca. genes in the nucleus and another 150– in the plastid, this accounts for less than 50% of the ancestral endosymbiont genome. finally, we discuss briefly the relevance of our findings to the extent of'real' primary egt. our results are based on detecting phylogenetic signal in anciently diverged sequences using standard computational approaches. given that phylogenetic signal was most certainly lost from many cyanobacterial genes over the greater than billion years since the endosymbiotic event, it will most certainly be impossible to detect all instances of egt. therefore, how significant is our number of cyanobacterial-derived genes in chlamydomonas? in our opinion, this must clearly be viewed as an underestimate of the actual number. many cyanobacterial genes may simply no longer be detectable as such due to highly diverged sequences. a more accurate estimate may come from increased taxon sampling in the phylogenomic analyses to allow these methods to potentially identify a larger proportion of highly diverged homologs among the plantae. we also need to gain a better understanding of the process of phylogenetic signal decay among this gene set so that extrapolations to the true value are more realistic. availability and requirements phylosort is available as an open-source under the gnu general public license via www at http://phylosort.sourceforge.net. phylosort is implemented in java under compliance with compiler version. it has been successfully tested on machines running the operating systems sun os, red hat, windows, and mac os. 
phylosort can be applied to quantify the number of genes that support different evolutionary hypotheses such as a taxonomic classification or endosymbiotic or horizontal gene transfer events. phylosort is an open-source tool to sort phylogenetic trees by searching for user specified subtrees that contain a monophyletic group of interest defined by operational taxonomic units in a phylogenomic context. plastid establishment precipitated extensive gene transfer from the endosymbiont to the nuclear genome of the'host'. in our application, we demonstrate that cyanobacteria account for –6% of the protein-coding genes in the nuclear genome of chlamydomonas. phylogenomic pipelines generate a large collection of phylogenetic trees that require manual inspection to answer questions about gene or genome evolution. using phylosort, we identified chlamydomonas genes of putative cyanobacterial origin, of which had bootstrap support values ≥ 50% for the grouping of the algal and cyanobacterial homologs. in the case of primary endosymbiosis, a heterotrophic protist engulfed a cyanobacterium, giving rise to the first photosynthetic eukaryote. estimating the magnitude of this endosymbiotic gene transfer and determining the functions of the prokaryotic genes remain controversial issues. to facilitate this procedure, we developed phylosort to rapidly search large collection of trees for monophyletic relationships.
 phylogenomic pipelines generate a large collection of phylogenetic trees that require manual inspection to answer questions about gene or genome evolution. a notable application of phylogenomics is to photosynthetic organelle endosymbiosis. in the case of primary endosymbiosis, a heterotrophic protist engulfed a cyanobacterium, giving rise to the first photosynthetic eukaryote. plastid establishment precipitated extensive gene transfer from the endosymbiont to the nuclear genome of the'host'. estimating the magnitude of this endosymbiotic gene transfer and determining the functions of the prokaryotic genes remain controversial issues. we used phylogenomics to study egt in the model green alga chlamydomonas reinhardtii. to facilitate this procedure, we developed phylosort to rapidly search large collection of trees for monophyletic relationships. here we present phylosort and its application to estimating egt in chlamydomonas.  phylosort is an open-source tool to sort phylogenetic trees by searching for user specified subtrees that contain a monophyletic group of interest defined by operational taxonomic units in a phylogenomic context. using phylosort, we identified chlamydomonas genes of putative cyanobacterial origin, of which had bootstrap support values ≥ 50% for the grouping of the algal and cyanobacterial homologs.  phylosort can be applied to quantify the number of genes that support different evolutionary hypotheses such as a taxonomic classification or endosymbiotic or horizontal gene transfer events. in our application, we demonstrate that cyanobacteria account for –6% of the protein-coding genes in the nuclear genome of chlamydomonas. 