The split genes of Nanoarchaeum equitans are an ancestral character

The introns early hypothesis predicts that introns were fundamental in assembling the first genes. In Nanoarchaeum equitans some genes are split into two. If these split genes were the ancestral forms, as suggested by the introns early hypothesis, then the end-beginning of the two parts of the split protein in a multiple alignment with the orthologous proteins from the Eukarya and Arachaea domains should make a clear prediction on where the intron in the homologous eukaryotic gene should be positioned. The analysis has shown that the introns are in this position, which is therefore predictable on the basis of the split proteins of N. equitans. This corroborates the hypothesis that the split genes of N. equitans are the plesiomorphic forms of these genes. If true, this would show that the origin of genes was polyphyletic as the monophyletic origin hypothesis would deny the existence, in a real organism, of these ancestral (split) genes, which imply that they were assembled late on and after the domains of life were established. Furthermore, it would seem that hyperthermophily is also an ancestral trait because it is linked to a split gene in N. equitans.     

The non-monophyletic origin of the tRNA molecule and the origin of genes only after the evolutionary stage of the last universal common ancestor (LUCA)

A model has been proposed suggesting that the tRNA molecule must have originated by direct duplication of an RNA hairpin structure [Di Giulio, M., 1992. On the origin of the transfer RNA molecule. J. Theor. Biol. 159, 199-214]. A non-monophyletic origin of this molecule has also been theorized [Di Giulio, M., 1999. The non-monophyletic origin of tRNA molecule. J. Theor. Biol. 197, 403-414]. In other words, the tRNA genes evolved only after the evolutionary stage of the last universal common ancestor (LUCA) through the assembly of two minigenes codifying for different RNA hairpin structures, which is what the exon theory of genes suggests when it is applied to the model of tRNA origin. Recent observations strongly corroborate this theorization because it has been found that some tRNA genes are completely separate in two minigenes codifying for the 5' and 3' halves of this molecule [Randau, L., et al., 2005a. Nanoarchaeum equitans creates functional tRNAs from separate genes for their 5'- and 3'-halves. Nature 433, 537-541]. In this paper it is shown that these tRNA genes codifying for the 5' and 3' halves of this molecule are the ancestral form from which the tRNA genes continuously codifying for the complete tRNA molecule are thought to have evolved. This, together with the very existence of completely separate tRNA genes codifying for their 5' and 3' halves, proves a non-monophyletic origin for tRNA genes, as a monophyletic origin would exclude the existence of these genes which have, on the contrary, been observed. Here the polyphyletic origin of genes codifying for proteins is also suggested and discussed. Moreover, a hypothesis is advanced to suggest that the LUCA might have had a fragmented genome made up of RNA and the possibility that 'Paleokaryotes' may exist is outlined. Finally, the characteristic of the indivisibility of homology that these polyphyletic origins seem to remove at the sequence level is discussed.     

On How Many Fundamental Kinds of Cells are Present on Earth: Looking for Phylogenetic Traits that Would Allow the Identification of the Primary Lines of Descent

The phylogenetic analyses as far as the identification of the number of domains of life is concerned have not reached a clear conclusion. In the attempt to improve this circumstance, I introduce the concept that the amino acids codified in the genetic code might be of markers with outstanding phylogenetic power. In particular, I hypothesise the existence of a biosphere populated, for instance, by three groups of organisms having different genetic codes because codifying at least a different amino acid. Evidently, these amino acids would mark the proteins that are present in the three groups of organisms in an unambiguous way. Therefore, in essence, this mark would not be other than the one that we usually try to make in the phylogenetic analyses in which we transform the protein sequences in phylogenetic trees, for the purpose to identify, for example, the domains of life. Indeed, this mark would allow to classify proteins without performing phylogenetic analyses because proteins belonging to a group of organisms would be recognisable as marked in a natural way by at least a different amino acid among the diverse groups of organisms. This conceptualisation answers the question of how many fundamental kinds of cells have evolved from the Last Universal Common Ancestor (LUCA), as the genetic code has unique proprieties that make the codified amino acids excellent phylogenetic markers. The presence of the formyl-methionine only in proteins of bacteria would mark them and would identify these as domain of life. On the other hand, the presence of pyrrolysine in the genetic code of the euryarchaeota would identify them such as another fundamental kind of cell evolved from the LUCA. Overall, the phylogenetic distribution of formyl-methionine and pyrrolysine would identify at least two domains of life—Bacteria and Archaea—but their number might be actually four; that is to say, Bacteria, Euryarchaeota, archeobacteria that are not euryarchaeota and Eukarya. The usually accepted domains of life represented by Bacteria, Archaea and Eukarya are not compatible with the phylogenetic distribution of these two amino acids and therefore this last classification might be mistaken.     

GENETIC EVIDENCE FOR REPRODUCTIVE ISOLATION AND MULTIPLE ORIGINS OF SYMPATRIC TROPHIC ECOTYPES OF WHITEFISH (COREGONUS)

We assessed variation in mitochondrial DNA (mtDNA) by restriction fragment length polymorphism (RFLP) analysis and in nuclear genes by allozyme analysis among sympatric pairs of limnetic and benthic ecotypes of whitefish (Coregonus) coexisting in three lakes of southern Yukon to address three evolutionary questions regarding their origins. Are sympatric low and high gill-raker count ecotypes genetically differentiated? Are they issued from monophyletic or polyphyletic evolutionary events? If they are polyphyletic in origins, did they originate from multiple allopatric speciation events or intralacustrine radiation? Our results corroborated previous genetic and ecological studies of these ecotypes, indicating that they represent genetically distinct reproductive units, and therefore refuting the hypothesis of phenotypic polymorphism within a single population. However, the amount of gene flow between ecotypes varied among lakes, correlating with the extent of morphological differentiation and the potential for premating reproductive isolation. The results indicated a polyphyletic origin of ecotypes whereby each of them have been expressed independently more than once. In the two lakes of Squanga Creek drainage, the existence of sympatric pairs was best explained by the secondary contact of two monophyletic whitefish groups that evolved in allopatry during the last glaciation events. In Dezadeash L. of Alsek R. drainage, our results could not verify either sympatric or allopatric (or microallopatric) origin of ecotypes. Regardless of the mode of speciation involved in their origins, these sympatric whitefish populations provided further evidence that Pleistocene glaciation events created conditions favoring rapid divergence and phenotypic differentiation among northern freshwater fishes.     

The universal ancestor lived in a thermophilic or hyperthermophilic environment

Galtier et al. (Science 1999, 283, 220-221) exploit the correlation between the optimal growth temperature in prokaryotes and the G+C content of rRNAs and establish that the last universal common ancestor (LUCA) lived in a mesophilic environment. This result was achieved by estimating the G+C content of the ancestral sequences of the rRNAs of the LUCA through use of a complex Markov model. I have re-analysed their alignments of the rDNAs with maximum parsimony and I have found that their result is not robust and is, in all likelihood, incorrect. In particular, the rRNA ancestral sequences reconstructed with maximum parsimony from these rDNA alignments as well as those reconstructed after eliminating all the sites that turn out to be ambiguous to the parsimony algorithm and to a site-by-site inspection of these alignments, are such as to suggest that the LUCA lived in a thermophilic or hyperthermophilic environment. This finding is also supported by some tRNA ancestral sequences. The main conclusion of this analysis is that if the LUCA was a progenote then the origin of life might have taken place at a high temperature.     

Morphological convergence characterizes the evolution of Xanthophyceae (Heterokontophyta): evidence from nuclear SSU rDNA and plastidial rbcL genes

Xanthophyceae are a group of heterokontophyte algae. Few molecular studies have investigated the evolutionary history and phylogenetic relationships of this class. We sequenced the nuclear-encoded SSU rDNA and chloroplast-encoded rbcL genes of several xanthophycean species from different orders, families, and genera. Neither SSU rDNA nor rbcL genes show intraspecific sequence variation and are good diagnostic markers for characterization of problematic species. New sequences, combined with those previously available, were used to create different multiple alignments. Analyses included sequences from 26 species of Xanthophyceae plus three Phaeothamniophyceae and two Phaeophyceae taxa used as outgroups. Phylogenetic analyses were performed according to Bayesian inference, maximum likelihood, and maximum parsimony methods. We explored effects produced on the phylogenetic outcomes by both taxon sampling as well as selected genes. Congruent results were obtained from analyses performed on single gene multiple alignments as well as on a data set including both SSU rDNA and rbcL sequences. Trees obtained in this study show that several currently recognized xanthophycean taxa do not form monophyletic groups. The order Mischococcales is paraphyletic, while Tribonematales and Botrydiales are polyphyletic even if evidence for the second order is not conclusive. Botrydiales and Vaucheriales, both including siphonous taxa, do not form a clade. The families Botrydiopsidaceae, Botryochloridaceae, and Pleurochloridaceae as well as the genera Botrydiopsis and Chlorellidium are polyphyletic. The Centritractaceae and the genus Bumilleriopsis also appear to be polyphyletic but their monophyly cannot be completely rejected with current evidence. Our results support morphological convergence at any taxonomic rank in the evolution of the Xanthophyceae. Finally, our phylogenetic analyses exclude an origin of the Xanthophyceae from a Vaucheria-like ancestor and favor a single early origin of the coccoid cell form.     

Modern mRNA proofreading and repair: clues that the last universal common ancestor possessed an RNA genome?

RNA repair has now been demonstrated to be a genuine biological process and appears to be present in all three domains of life. In this article, we consider what this might mean for the transition from an early RNA-dominated world to modern cells possessing genetically encoded proteins and DNA. There are significant gaps in our understanding of how the modern protein-DNA world could have evolved from a simpler system, and it is currently uncertain whether DNA genomes evolved once or twice. Against this backdrop, the discovery of RNA repair in modern cells is timely food for thought and brings us conceptually one step closer to understanding how RNA genomes were replaced by DNA genomes. We have examined the available literature on multisubunit RNA polymerase structure and function and conclude that a strong case can be made that the Last Universal Common Ancestor (LUCA) possessed a repair-competent RNA polymerase, which would have been capable of acting on an RNA genome. However, while this lends credibility to the proposal that the LUCA had an RNA genome, the alternative, that LUCA had a DNA genome, cannot be completely ruled out.     

The structure of Prp40 FF1 domain and its interaction with the crn-TPR1 motif of Clf1 gives a new insight into the binding mode of FF domains

The yeast splicing factor Prp40 (pre-mRNA processing protein 40) consists of a pair of WW domains followed by several FF domains. The region comprising the FF domains has been shown to associate with the 5' end of U1 small nuclear RNA and to interact directly with two proteins, the Clf1 (Crooked neck-like factor 1) and the phosphorylated repeats of the C-terminal domain of RNA polymerase II (CTD-RNAPII). In this work we reported the solution structure of the first FF domain of Prp40 and the identification of a novel ligand-binding site in FF domains. By using chemical shift assays, we found a binding site for the N-terminal crooked neck tetratricopeptide repeat of Clf1 that is distinct and structurally separate from the previously identified CTD-RNAPII binding pocket of the FBP11 (formin-binding protein 11) FF1 domain. No interaction, however, was observed between the Prp40 FF1 domain and three different peptides derived from the CTD-RNAPII protein. Indeed, the equivalent CTD-RNAPII-binding site in the Prp40 FF1 domain is predominantly negatively charged and thus unfavorable for an interaction with phosphorylated peptide sequences. Sequence alignments and phylogenetic tree reconstructions using the FF domains of three functionally related proteins, Prp40, FBP11, and CA150, revealed that Prp40 and FBP11 are not orthologous proteins and supported the different ligand specificities shown by their respective FF1 domains. Our results also revealed that not all FF domains in Prp40 are functionally equivalent. We proposed that at least two different interaction surfaces exist in FF domains that have evolved to recognize distinct binding motifs.     

Combining phylogenetic data with co-regulated genes to identify regulatory motifs

MOTIVATION: Discovery of regulatory motifs in unaligned DNA sequences remains a fundamental problem in computational biology. Two categories of algorithms have been developed to identify common motifs from a set of DNA sequences. The first can be called a 'multiple genes, single species' approach. It proposes that a degenerate motif is embedded in some or all of the otherwise unrelated input sequences and tries to describe a consensus motif and identify its occurrences. It is often used for co-regulated genes identified through experimental approaches. The second approach can be called 'single gene, multiple species'. It requires orthologous input sequences and tries to identify unusually well conserved regions by phylogenetic footprinting. Both approaches perform well, but each has some limitations. It is tempting to combine the knowledge of co-regulation among different genes and conservation among orthologous genes to improve our ability to identify motifs. RESULTS: Based on the Consensus algorithm previously established by our group, we introduce a new algorithm called PhyloCon (Phylogenetic Consensus) that takes into account both conservation among orthologous genes and co-regulation of genes within a species. This algorithm first aligns conserved regions of orthologous sequences into multiple sequence alignments, or profiles, then compares profiles representing non-orthologous sequences. Motifs emerge as common regions in these profiles. Here we present a novel statistic to compare profiles of DNA sequences and a greedy approach to search for common subprofiles. We demonstrate that PhyloCon performs well on both synthetic and biological data. AVAILABILITY: Software available upon request from the authors. http://ural.wustl.edu/softwares.html     

The universal ancestor was a thermophile or a hyperthermophile: tests and further evidence

The existence of a correlation between the optimal growth temperature of various organisms and a thermophily index (based on the propensity of amino acids to enter more frequently into the proteins of thermophiles/hyperthermophiles) allows inferences to be made on the mesophilic or thermophilic nature of the last universal common ancestor (LUCA). By reconstructing the ancestral sequences of the various ancestors using methods based on maximum likelihood and maximum parsimony, these sequences can be attributed to the mesophiles or (hyper)thermophiles and the following conclusions can be drawn. (1) There is no evidence that the LUCA might have been a mesophile and observations seem to imply that the LUCA was a thermophile or a hyperthermophile; (2) The ancestors of the Archaea and Bacteria domains seem to be (hyper)thermophiles while that of the Eukarya domain turns out to be a mesophile. These conclusions are independent of both (i) where the root is located on the topology of the universal tree (based on that of the small subunit ribosomal RNA) and (ii) the presence of hyperthermophile bacteria near the node of the Bacteria domain ancestor. These conclusions are easier to interpret in the light of the hypotheses that see the origin of life taking place at a high temperature.     

Stochastic models for horizontal gene transfer: taking a random walk through tree space

Horizontal gene transfer (HGT) plays a critical role in evolution across all domains of life with important biological and medical implications. I propose a simple class of stochastic models to examine HGT using multiple orthologous gene alignments. The models function in a hierarchical phylogenetic framework. The top level of the hierarchy is based on a random walk process in "tree space" that allows for the development of a joint probabilistic distribution over multiple gene trees and an unknown, but estimable species tree. I consider two general forms of random walks. The first form is derived from the subtree prune and regraft (SPR) operator that mirrors the observed effects that HGT has on inferred trees. The second form is based on walks over complete graphs and offers numerically tractable solutions for an increasing number of taxa. The bottom level of the hierarchy utilizes standard phylogenetic models to reconstruct gene trees given multiple gene alignments conditional on the random walk process. I develop a well-mixing Markov chain Monte Carlo algorithm to fit the models in a Bayesian framework. I demonstrate the flexibility of these stochastic models to test competing ideas about HGT by examining the complexity hypothesis. Using 144 orthologous gene alignments from six prokaryotes previously collected and analyzed, Bayesian model selection finds support for (1) the SPR model over the alternative form, (2) the 16S rRNA reconstruction as the most likely species tree, and (3) increased HGT of operational genes compared to informational genes.     

Molecular evolution of lysin motif-type receptor-like kinases in plants

The lysin motif (LysM) domain is an ancient and ubiquitous protein module that binds peptidoglycan and structurally related molecules. A genomic survey in a large number of species spanning all kingdoms reveals that the combination of LysM and receptor kinase domains is present exclusively in plants. However, the particular biological functions and molecular evolution of this gene family remain largely unknown. We show that LysM domains in plant LysM proteins are highly diversified and that a minimum of six distinct types of LysM motifs exist in plant LysM kinase proteins and five additional types of LysM motifs exist in nonkinase plant LysM proteins. Further, motif similarities suggest that plant LysM motifs are ancient. Although phylogenetic signals are not sufficient to resolve the earliest relationships, plant LysM motifs may have arisen through common ancestry with LysM motifs in other kingdoms. Within plants, the gene family has evolved through local and segmental duplications. The family has undergone further duplication and diversification in legumes, where some LysM kinase genes function as receptors for bacterial nodulation factor. Two pairs of homeologous regions were identified in soybean (Glycine max) based on microsynteny and fluorescence in situ hybridization. Expression data show that most plant LysM kinase genes are expressed predominantly in the root and that orthologous LysM kinase genes share similar tissue expression patterns. We also examined synteny around plant LysM kinase genes to help reconstruct scenarios for the evolution of this important gene family.     

Evolution of DNA binding motifs and operators

A gene regulatory protein with helix-turn-helix (HTH) DNA-binding motif, GalS contains a functional operator within the DNA sequences encoding the HTH region (Nature 369 (1994) 314). We searched for operator-like sequences within the DNA sequences encoding the DNA binding motifs of other regulatory proteins. Five such proteins, DeoR, CytR, LRP, LuxR and PurR, were found to have actual operator or operator-like sequences in the DNA sequences encoding the DNA-binding motif. Except DeoR, all of them including GalS, are known to be auto-regulated. Auto-regulation in case of DeoR has not been investigated. Seven other proteins containing a HTH motif, do not have operator-like sequences in the DNA sequences encoding the HTH motif; none of them, except MerR, are known to be auto-regulated. The DNA binding proteins may have evolved from a common ancestor containing a DNA binding site within its gene segment that encodes the DNA-binding motif to facilitate auto-regulation. We have discussed current evidence for monophyletic or polyphyletic origin of such sequences.     

Evolution of Buchlo? dactyloides based on cloning and sequencing of matK, rbcL, and cob genes from plastid and mitochondrial genomes.

Buffalograss (Buchlo? dactyloides (Nutt.) Englem), a C4 turfgrass species, is native to the Great Plains region of North America. The evolutionary implications of buffalograss are unclear. Sequencing of rbcL and matK genes from plastid and the cob gene from mitochondrial genomes was examined to elucidate buffalo grass evolution. This study is the first to report sequencing of these genes from organelle genomes in the genus Buchlo?. Comparisons of sequence data from the mitochondrial and plastid genome revealed that all genotypes contained the same cytoplasmic origin. There were some rearrangements detected in mitochondrial genome. The buffalograss genome appears to have evolved through the rearrangements of convergent subgenomic domains. Combined analyses of plastid genes suggest that the evolutionary process in Buchlo? accessions studied was monophyletic rather than polyphyletic. However, since plastid and mitochondrial genomes are generally uniparentally inherited, the evolutionary history of these genomes may not reflect the evolutionary history of the organism, especially in a species in which out-crossing is common. The sequence information obtained from this study can be used as a genome-specific marker for investigation of the buffalograss polyploidy complex and testing of the mode of plastid and mitochondrial transmission in genus Buchlo?.     

Phylogenomics of trophically diverse cichlids disentangles processes driving adaptive radiation and repeated trophic transitions

Cichlid fishes of the tribe Tropheini are a striking case of adaptive radiation, exemplifying multiple trophic transitions between herbivory and carnivory occurring in sympatry with other established cichlid lineages. Tropheini evolved highly specialized eco‐morphologies to exploit similar trophic niches in different ways repeatedly and rapidly. To better understand the evolutionary history and trophic adaptations of this lineage, we generated a dataset of 532 targeted loci from 21 out of the 22 described Tropheini species. We resolved the Tropheini into seven monophyletic genera and discovered one to be polyphyletic. The polyphyletic genus, Petrochromis, represents three convergent origins of the algae grazing trophic specialization. This repeated evolution of grazing may have been facilitated by adaptive introgression as we found evidence for gene flow among algae grazing genera. We also found evidence of gene flow among algae browsing genera, but gene flow was restricted between herbivorous and carnivorous genera. Furthermore, we observed no evidence supporting a hybrid origin of this radiation. Our molecular evolutionary analyses suggest that opsin genes likely evolved in response to selection pressures associated with trophic ecology in the Tropheini. We found surprisingly little evidence of positive selection in coding regions of jaw‐shaping genes in this trophically diverse lineage. This suggests low degrees of freedom for further change in these genes, and possibly a larger role for regulatory variation in driving jaw adaptations. Our study emphasizes Tropheini cichlids as an important model for studying the evolution of trophic specialization and its role in speciation.     

Genomics and early cellular evolution. The origin of the DNA world.

The sequencing of several genomes from each of the three domains of life (Archaea, Bacteria and Eukarya) has provided a huge amount of data that can be used to gain insight about early cellular evolution. Some features of the universal tree of life based on rRNA polygenies have been confirmed, such as the division of the cellular living world into three domains. The monophyly of each domain is supported by comparative genomics. However, the hyperthermophilic nature of the 'last universal common ancestor' (LUCA) is not confirmed. Comparative genomics has revealed that gene transfers have been (and still are) very frequent in genome evolution. Nevertheless, a core of informational genes appears more resistant to transfer, testifying for a close relationship between archaeal and eukaryal informational processes. This observation can be explained either by a common unique history between Archaea and Eukarya or by an atypical evolution of these systems in Bacteria. At the moment, comparative genomics still does not allow to choose between a simple LUCA, possibly with an RNA genome, or a complex LUCA, with a DNA genome and informational mechanisms similar to those of Archaea and Eukarya. Further comparative studies on informational mechanisms in the three domains should help to resolve this critical question. The role of viruses in the origin and evolution of DNA genomes also appears an area worth of active investigations. I suggest here that DNA and DNA replication mechanisms appeared first in the virus world before being transferred into cellular organisms.     

The Universal Ancestor and the Ancestor of Bacteria Were Hyperthermophiles

The definition of the node of the last universal common ancestor (LUCA) is justified in a topology of the unrooted universal tree. This definition allows previous analyses based on paralogous proteins to be extended to orthologous ones. In particular, the use of a thermophily index (based on the amino acids propensity to enter the [hyper] thermophile proteins more frequently) and its correlation with the optimal growth temperature of the various organisms allow inferences to be made on the habitat in which the LUCA lived. The reconstruction of ancestral sequences by means of the maximum likelihood method and their attribution to the set of mesophilic or hyperthermophilic sequences have led to the following conclusions: the LUCA was a hyperthermophile organism, as were the ancestors of the Archaea and Bacteria domains, while the ancestor of the Eukarya domain was a mesophile. These conclusions are independent of the presence of hyperthermophile bacteria in the sample of sequences used in the analysis and are therefore independent of whether or not these are the first lines of divergence in the Bacteria domain, as observed in the topology of the universal tree of ribosomal RNA. These conclusions are thus more easily understood under the hypothesis that the origin of life took place at a high temperature.     

The origins of modern proteomes

Distributions of phylogenetically related protein domains (fold superfamilies), or FSFs, among the three Superkingdoms (trichotomy) are assessed. Very nearly 900 of the 1200 FSFs of the trichotomy are shared by two or three Superkingdoms. Parsimony analysis of FSF distributions suggests that the FSF complement of the last common ancestor to the trichotomy was more like that of modern eukaryotes than that of archaea and bacteria. Studies of length distributions among members of orthologous families of proteins present in all three Superkingdoms reveal that such lengths are significantly longer among eukaryotes than among bacteria and archaea. The data also reveal that proteins lengths of eukaryotes are more broadly distributed than they are within archaeal and bacterial members of the same orthologous families. Accordingly, selective pressure for a minimal size is significantly greater for orthologous protein lengths in archaea and bacteria than in eukaryotes. Alignments of orthologous proteins of archaea, bacteria and eukaryotes are characterized by greater sequence variation at their N-terminal and C-terminal domains, than in their central cores. Length variations tend to be localized in the terminal sequences; the conserved sequences of orthologous families are localized in a central core. These data are consistent with the interpretation that the genomes of the last common ancestor (LUCA) encoded a cohort of FSFs not very different from that of modern eukaryotes. Divergence of bacterial and archaeal genomes from that common ancestor may have been accompanied by more intensive reductive evolution of proteomes than that expressed in eukaryotes. Dollo's Law suggests that the evolution of novel FSFs is a very slow process, while laboratory experiments suggests that novel protein genesis from preexisting FSFs can be relatively rapid. Reassortment of FSFs to create novel proteins may have been mediated by genetic recombination before the advent of more efficient splicing mechanisms.     

Multiple Horizontal Gene Transfer Events and Domain Fusions Have Created Novel Regulatory and Metabolic Networks in the Oomycete Genome

Complex enzymes with multiple catalytic activities are hypothesized to have evolved from more primitive precursors. Global analysis of the Phytophthora sojae genome using conservative criteria for evaluation of complex proteins identified 273 novel multifunctional proteins that were also conserved in P. ramorum. Each of these proteins contains combinations of protein motifs that are not present in bacterial, plant, animal, or fungal genomes. A subset of these proteins were also identified in the two diatom genomes, but the majority of these proteins have formed after the split between diatoms and oomycetes. Documentation of multiple cases of domain fusions that are common to both oomycetes and diatom genomes lends additional support for the hypothesis that oomycetes and diatoms are monophyletic. Bifunctional proteins that catalyze two steps in a metabolic pathway can be used to infer the interaction of orthologous proteins that exist as separate entities in other genomes. We postulated that the novel multifunctional proteins of oomycetes could function as potential Rosetta Stones to identify interacting proteins of conserved metabolic and regulatory networks in other eukaryotic genomes. However ortholog analysis of each domain within our set of 273 multifunctional proteins against 39 sequenced bacterial and eukaryotic genomes, identified only 18 candidate Rosetta Stone proteins. Thus the majority of multifunctional proteins are not Rosetta Stones, but they may nonetheless be useful in identifying novel metabolic and regulatory networks in oomycetes. Phylogenetic analysis of all the enzymes in three pathways with one or more novel multifunctional proteins was conducted to determine the probable origins of individual enzymes. These analyses revealed multiple examples of horizontal transfer from both bacterial genomes and the photosynthetic endosymbiont in the ancestral genome of Stramenopiles. The complexity of the phylogenetic origins of these metabolic pathways and the paucity of Rosetta Stones relative to the total number of multifunctional proteins suggests that the proteome of oomycetes has few features in common with other Kingdoms.