Comparative genomics and bioenergetics

ATP synthase: Evolution, energetics, and membrane interactions

The synthesis of ATP, life’s “universal energy currency,” is the most prevalent chemical reaction in biological systems and is responsible for fueling nearly all cellular processes, from nerve impulse propagation to DNA synthesis. ATP synthases, the family of enzymes that carry out this endless task, are nearly as ubiquitous as the energy-laden molecule they are responsible for making. The F-type ATP synthase (F-ATPase) is found in every domain of life and has facilitated the survival of organisms in a wide range of habitats, ranging from the deep-sea thermal vents to the human intestine. Accordingly, there has been a large amount of work dedicated toward understanding the structural and functional details of ATP synthases in a wide range of species. Less attention, however, has been paid toward integrating these advances in ATP synthase molecular biology within the context of its evolutionary history. In this review, we present an overview of several structural and functional features of the F-type ATPases that vary across taxa and are purported to be adaptive or otherwise evolutionarily significant: ion channel selectivity, rotor ring size and stoichiometry, ATPase dimeric structure and localization in the mitochondrial inner membrane, and interactions with membrane lipids. We emphasize the importance of studying these features within the context of the enzyme’s particular lipid environment. Just as the interactions between an organism and its physical environment shape its evolutionary trajectory, ATPases are impacted by the membranes within which they reside. We argue that a comprehensive understanding of the structure, function, and evolution of membrane proteins—including ATP synthase—requires such an integrative approach.

Bacterial Oxidases of the Cytochrome bd Family: Redox Enzymes of Unique Structure, Function, and Utility As Drug Targets

Significance: Cytochrome bd is a ubiquinol:oxygen oxidoreductase of many prokaryotic respiratory chains with a unique structure and functional characteristics. Its primary role is to couple the reduction of molecular oxygen, even at submicromolar concentrations, to water with the generation of a proton motive force used for adenosine triphosphate production. Cytochrome bd is found in many bacterial pathogens and, surprisingly, in bacteria formally denoted as anaerobes. It endows bacteria with resistance to various stressors and is a potential drug target. Recent Advances: We summarize recent advances in the biochemistry, structure, and physiological functions of cytochrome bd in the light of exciting new three-dimensional structures of the oxidase. The newly discovered roles of cytochrome bd in contributing to bacterial protection against hydrogen peroxide, nitric oxide, peroxynitrite, and hydrogen sulfide are assessed. Critical Issues: Fundamental questions remain regarding the precise delineation of electron flow within this multihaem oxidase and how the extraordinarily high affinity for oxygen is accomplished, while endowing bacteria with resistance to other small ligands. Future Directions: It is clear that cytochrome bd is unique in its ability to confer resistance to toxic small molecules, a property that is significant for understanding the propensity of pathogens to possess this oxidase. Since cytochrome bd is a uniquely bacterial enzyme, future research should focus on harnessing fundamental knowledge of its structure and function to the development of novel and effective antibacterial agents.

Systems level analysis of the Chlamydomonas reinhardtii metabolic network reveals variability in evolutionary co-conservation

Metabolic networks, which are mathematical representations of organismal metabolism, are reconstructed to provide computational platforms to guide metabolic engineering experiments and explore fundamental questions on metabolism. Systems level analyses, such as interrogation of phylogenetic relationships within the network, can provide further guidance on the modification of metabolic circuitries. Chlamydomonas reinhardtii, a biofuel relevant green alga that has retained key genes with plant, animal, and protist affinities, serves as an ideal model organism to investigate the interplay between gene function and phylogenetic affinities at multiple organizational levels. Here, using detailed topological and functional analyses, coupled with transcriptomics studies on a metabolic network that we have reconstructed for C. reinhardtii, we show that network connectivity has a significant concordance with the co-conservation of genes; however, a distinction between topological and functional relationships is observable within the network. Dynamic and static modes of co-conservation were defined and observed in a subset of gene-pairs across the network topologically. In contrast, genes with predicted synthetic interactions, or genes involved in coupled reactions, show significant enrichment for both shorter and longer phylogenetic distances. Based on our results, we propose that the metabolic network of C. reinhardtii is assembled with an architecture to minimize phylogenetic profile distances topologically, while it includes an expansion of such distances for functionally interacting genes. This arrangement may increase the robustness of C. reinhardtii's network in dealing with varied environmental challenges that the species may face. The defined evolutionary constraints within the network, which identify important pairings of genes in metabolism, may offer guidance on synthetic biology approaches to optimize the production of desirable metabolites.

Physiology, phylogeny, and LUCA

Genomes record their own history. But if we want to look all the way back to life's beginnings some 4 billion years ago, the record of microbial evolution that is preserved in prokaryotic genomes is not easy to read. Microbiology has a lot in common with geology in that regard. Geologists know that plate tectonics and erosion have erased much of the geological record, with ancient rocks being truly rare. The same is true of microbes. Lateral gene transfer (LGT) and sequence divergence have erased much of the evolutionary record that was once written in genomes, and it is not obvious which genes among sequenced genomes are genuinely ancient. Which genes trace to the last universal ancestor, LUCA? The classical approach has been to look for genes that are universally distributed. Another approach is to make all trees for all genes, and sift out the trees where signals have been overwritten by LGT. What is left ought to be ancient. If we do that, what do we find?

Evolution of the F0F1 ATP synthase complex in light of the patchy distribution of different bioenergetic pathways across prokaryotes

Bacteria and archaea are characterized by an amazing metabolic diversity, which allows them to persist in diverse and often extreme habitats. Apart from oxygenic photosynthesis and oxidative phosphorylation, well-studied processes from chloroplasts and mitochondria of plants and animals, prokaryotes utilize various chemo- or lithotrophic modes, such as anoxygenic photosynthesis, iron oxidation and reduction, sulfate reduction, and methanogenesis. Most bioenergetic pathways have a similar general structure, with an electron transport chain composed of protein complexes acting as electron donors and acceptors, as well as a central cytochrome complex, mobile electron carriers, and an ATP synthase. While each pathway has been studied in considerable detail in isolation, not much is known about their relative evolutionary relationships. Wanting to address how this metabolic diversity evolved, we mapped the distribution of nine bioenergetic modes on a phylogenetic tree based on 16S rRNA sequences from 272 species representing the full diversity of prokaryotic lineages. This highlights the patchy distribution of many pathways across different lineages, and suggests either up to 26 independent origins or 17 horizontal gene transfer events. Next, we used comparative genomics and phylogenetic analysis of all subunits of the F0F1 ATP synthase, common to most bacterial lineages regardless of their bioenergetic mode. Our results indicate an ancient origin of this protein complex, and no clustering based on bioenergetic mode, which suggests that no special modifications are needed for the ATP synthase to work with different electron transport chains. Moreover, examination of the ATP synthase genetic locus indicates various gene rearrangements in the different bacterial lineages, ancient duplications of atpI and of the beta subunit of the F0 subcomplex, as well as more recent stochastic lineage-specific and species-specific duplications of all subunits. We discuss the implications of the overall pattern of conservation and flexibility of the F0F1 ATP synthase genetic locus.

Evaluation of physical and functional protein-protein interaction prediction methods for detecting biological pathways

BACKGROUND

Cellular activities are governed by the physical and the functional interactions among several proteins involved in various biological pathways. With the availability of sequenced genomes and high-throughput experimental data one can identify genome-wide protein-protein interactions using various computational techniques. Comparative assessments of these techniques in predicting protein interactions have been frequently reported in the literature but not their ability to elucidate a particular biological pathway.

METHODS

Towards the goal of understanding the prediction capabilities of interactions among the specific biological pathway proteins, we report the analyses of 14 biological pathways of Escherichia coli catalogued in KEGG database using five protein-protein functional linkage prediction methods. These methods are phylogenetic profiling, gene neighborhood, co-presence of orthologous genes in the same gene clusters, a mirrortree variant, and expression similarity.

CONCLUSIONS

Our results reveal that the prediction of metabolic pathway protein interactions continues to be a challenging task for all methods which possibly reflect flexible/independent evolutionary histories of these proteins. These methods have predicted functional associations of proteins involved in amino acids, nucleotide, glycans and vitamins & co-factors pathways slightly better than the random performance on carbohydrate, lipid and energy metabolism. We also make similar observations for interactions involved among the environmental information processing proteins. On the contrary, genetic information processing or specialized processes such as motility related protein-protein linkages that occur in the subset of organisms are predicted with comparable accuracy. Metabolic pathways are best predicted by using neighborhood of orthologous genes whereas phyletic pattern is good enough to reconstruct central dogma pathway protein interactions. We have also shown that the effective use of a particular prediction method depends on the pathway under investigation. In case one is not focused on specific pathway, gene expression similarity method is the best option.

The cytochrome bd respiratory oxygen reductases

Cytochrome bd is a respiratory quinol: O₂ oxidoreductase found in many prokaryotes, including a number of pathogens. The main bioenergetic function of the enzyme is the production of a proton motive force by the vectorial charge transfer of protons. The sequences of cytochromes bd are not homologous to those of the other respiratory oxygen reductases, i.e., the heme-copper oxygen reductases or alternative oxidases (AOX). Generally, cytochromes bd are noteworthy for their high affinity for O₂ and resistance to inhibition by cyanide. In E. coli, for example, cytochrome bd (specifically, cytochrome bd-I) is expressed under O₂-limited conditions. Among the members of the bd-family are the so-called cyanide-insensitive quinol oxidases (CIO) which often have a low content of the eponymous heme d but, instead, have heme b in place of heme d in at least a majority of the enzyme population. However, at this point, no sequence motif has been identified to distinguish cytochrome bd (with a stoichiometric complement of heme d) from an enzyme designated as CIO. Members of the bd-family can be subdivided into those which contain either a long or a short hydrophilic connection between transmembrane helices 6 and 7 in subunit I, designated as the Q-loop. However, it is not clear whether there is a functional consequence of this difference. This review summarizes current knowledge on the physiological functions, genetics, structural and catalytic properties of cytochromes bd. Included in this review are descriptions of the intermediates of the catalytic cycle, the proposed site for the reduction of O₂, evidence for a proton channel connecting this active site to the bacterial cytoplasm, and the molecular mechanism by which a membrane potential is generated.

Comparative genomics analysis of completely sequenced microbial genomes reveals the ubiquity of N-linked glycosylation in prokaryotes

Glycosylation of proteins in prokaryotes has been known for the last few decades. Glycan structures and/or the glycosylation pathways have been experimentally characterized in only a small number of prokaryotes. Even this has become possible only during the last decade or so, primarily due to technological and methodological developments. Glycosylated proteins are diverse in their function and localization. Glycosylation has been shown to be associated with a wide range of biological phenomena. Characterization of the various types of glycans and the glycosylation machinery is critical to understand such processes. Such studies can help in the identification of novel targets for designing drugs, diagnostics, and engineering of therapeutic proteins. In view of this, the experimentally characterized pgl system of Campylobacter jejuni, responsible for N-linked glycosylation, has been used in this study to identify glycosylation loci in 865 prokaryotes whose genomes have been completely sequenced. Results from the present study show that only a small number of organisms have homologs for all the pgl enzymes and a few others have homologs for none of the pgl enzymes. Most of the organisms have homologs for only a subset of the pgl enzymes. There is no specific pattern for the presence or absence of pgl homologs vis-à-vis the 16S rRNA sequence-based phylogenetic tree. This may be due to differences in the glycan structures, high sequence divergence, horizontal gene transfer or non-orthologous gene displacement. Overall, the presence of homologs for pgl enzymes in a large number of organisms irrespective of their habitat, pathogenicity, energy generation mechanism, etc., hints towards the ubiquity of N-linked glycosylation in prokaryotes.

The origin of life and the last universal common ancestor: do we need a change of perspective?

A complete tree with roots, trunk and crown remains an appropriate model to represent all steps of life's development, from the emergence of a unique genetic code up to the last universal common ancestor and its further radiation. Catalytic closure of a mixture of prebiotic polymers is a heuristic alternative to the RNA world. Conjectures about emergence of life in an infinite multiverse should not confuse probability with possibility.

A phylogenetic mixture model for gene family loss in parasitic bacteria

Gene families are frequently gained and lost from prokaryotic genomes. It is widely believed that the rate of loss was accelerated for some but not all gene families in lineages that became parasites or endosymbionts. This leads to a form of heterotachy that may be responsible for the poor performance of phylogeny estimation based on gene content. We describe a mixture model that accounts for this heterotachy. We show that this model fits data on the distribution of gene families across bacteria from the COG database much better than previous models. However, it still favors an artifactual tree topology in which parasites form a clade over the more plausible 16S topology. In contrast to a previous model of genome dynamics, our model suggests that the ancestral bacterium had a small genome. We suggest that models of gene family gain and loss are likely to be more useful for understanding genome dynamics than for estimating phylogenetic trees.

Phylogenomic dating--a method of constraining the age of microbial taxa that lack a conventional fossil record

A phylogenomic dating approach was used to identify potential age constraints for multiple archaeal groups, many of which have no fossil, isotopic, or biomarker record. First, well-resolved phylogenetic trees were inferred with the use of multiple gene sequences obtained from whole genome sequences. Next, the ability to use oxygen as a terminal electron acceptor was coded into characters, and ancestral state reconstruction was used to identify clades with taxa that metabolize oxygen and likely had an aerobic ancestor. Next, the habitat of the ancestor was inferred. If the local presence of Cyanobacteria could be excluded from the putative ancestral habitat, then these clades would have originated after the rise in atmospheric oxygen 2.32 Ga. With this method, an upper age of 2.32 Ga (an "oxygen age constraint") is proposed for four major archaeal clades: the Sulfolobales, Thermoplasmatales, Thermoproteus neutrophilus/Pyrobaculum spp., and the Thermoproteales. It was also shown that the halophilic archaea likely had an aerobic common ancestor, yet the possibility of local oxygen oases before oxygenation of the atmosphere could not be formally rejected. Thus, an oxygen age constraint was not assessed for this group. This work suggests that many archaeal groups are not as ancient as many in the research community have previously assumed, and it provides a new method for establishing upper age constraints for major microbial groups that lack a conventional fossil record.

The origin, evolution and structure of the protein world

Contemporary protein architectures can be regarded as molecular fossils, historical imprints that mark important milestones in the history of life. Whereas sequences change at a considerable pace, higher-order structures are constrained by the energetic landscape of protein folding, the exploration of sequence and structure space, and complex interactions mediated by the proteostasis and proteolytic machineries of the cell. The survey of architectures in the living world that was fuelled by recent structural genomic initiatives has been summarized in protein classification schemes, and the overall structure of fold space explored with novel bioinformatic approaches. However, metrics of general structural comparison have not yet unified architectural complexity using the 'shared and derived' tenet of evolutionary analysis. In contrast, a shift of focus from molecules to proteomes and a census of protein structure in fully sequenced genomes were able to uncover global evolutionary patterns in the structure of proteins. Timelines of discovery of architectures and functions unfolded episodes of specialization, reductive evolutionary tendencies of architectural repertoires in proteomes and the rise of modularity in the protein world. They revealed a biologically complex ancestral proteome and the early origin of the archaeal lineage. Studies also identified an origin of the protein world in enzymes of nucleotide metabolism harbouring the P-loop-containing triphosphate hydrolase fold and the explosive discovery of metabolic functions that recapitulated well-defined prebiotic shells and involved the recruitment of structures and functions. These observations have important implications for origins of modern biochemistry and diversification of life.

Transcriptional map of respiratory versatility in the hyperthermophilic crenarchaeon Pyrobaculum aerophilum

Hyperthermophilic crenarchaea in the genus Pyrobaculum are notable for respiratory versatility, but relatively little is known about the genetics or regulation of crenarchaeal respiratory pathways. We measured global gene expression in Pyrobaculum aerophilum cultured with oxygen, nitrate, arsenate and ferric iron as terminal electron acceptors to identify transcriptional patterns that differentiate these pathways. We also compared genome sequences for four closely related species with diverse respiratory characteristics (Pyrobaculum arsenaticum, Pyrobaculum calidifontis, Pyrobaculum islandicum, and Thermoproteus neutrophilus) to identify genes associated with different respiratory capabilities. Specific patterns of gene expression in P. aerophilum were associated with aerobic respiration, nitrate respiration, arsenate respiration, and anoxia. Functional predictions based on these patterns include separate cytochrome oxidases for aerobic growth and oxygen scavenging, a nitric oxide-responsive transcriptional regulator, a multicopper oxidase involved in denitrification, and an archaeal arsenate respiratory reductase. We were unable to identify specific genes for iron respiration, but P. aerophilum exhibited repressive transcriptional responses to iron remarkably similar to those controlled by the ferric uptake regulator in bacteria. Together, these analyses present a genome-scale view of crenarchaeal respiratory flexibility and support a large number of functional and regulatory predictions for further investigation. The complete gene expression data set can be viewed in genomic context with the Archaeal Genome Browser at archaea.ucsc.edu.

The last universal common ancestor: emergence, constitution and genetic legacy of an elusive forerunner

BACKGROUND

Since the reclassification of all life forms in three Domains (Archaea, Bacteria, Eukarya), the identity of their alleged forerunner (Last Universal Common Ancestor or LUCA) has been the subject of extensive controversies: progenote or already complex organism, prokaryote or protoeukaryote, thermophile or mesophile, product of a protracted progression from simple replicators to complex cells or born in the cradle of "catalytically closed" entities? We present a critical survey of the topic and suggest a scenario.

RESULTS

LUCA does not appear to have been a simple, primitive, hyperthermophilic prokaryote but rather a complex community of protoeukaryotes with a RNA genome, adapted to a broad range of moderate temperatures, genetically redundant, morphologically and metabolically diverse. LUCA's genetic redundancy predicts loss of paralogous gene copies in divergent lineages to be a significant source of phylogenetic anomalies, i.e. instances where a protein tree departs from the SSU-rRNA genealogy; consequently, horizontal gene transfer may not have the rampant character assumed by many. Examining membrane lipids suggest LUCA had sn1,2 ester fatty acid lipids from which Archaea emerged from the outset as thermophilic by "thermoreduction," with a new type of membrane, composed of sn2,3 ether isoprenoid lipids; this occurred without major enzymatic reconversion. Bacteria emerged by reductive evolution from LUCA and some lineages further acquired extreme thermophily by convergent evolution. This scenario is compatible with the hypothesis that the RNA to DNA transition resulted from different viral invasions as proposed by Forterre. Beyond the controversy opposing "replication first" to metabolism first", the predictive arguments of theories on "catalytic closure" or "compositional heredity" heavily weigh in favour of LUCA's ancestors having emerged as complex, self-replicating entities from which a genetic code arose under natural selection.

CONCLUSION

Life was born complex and the LUCA displayed that heritage. It had the "body "of a mesophilic eukaryote well before maturing by endosymbiosis into an organism adapted to an atmosphere rich in oxygen. Abundant indications suggest reductive evolution of this complex and heterogeneous entity towards the "prokaryotic" Domains Archaea and Bacteria. The word "prokaryote" should be abandoned because epistemologically unsound.

REVIEWERS

This article was reviewed by Anthony Poole, Patrick Forterre, and Nicolas Galtier.

Nitrogen metabolism in haloarchaea

The nitrogen cycle (N-cycle), principally supported by prokaryotes, involves different redox reactions mainly focused on assimilatory purposes or respiratory processes for energy conservation. As the N-cycle has important environmental implications, this biogeochemical cycle has become a major research topic during the last few years. However, although N-cycle metabolic pathways have been studied extensively in Bacteria or Eukarya, relatively little is known in the Archaea. Halophilic Archaea are the predominant microorganisms in hot and hypersaline environments such as salted lakes, hot springs or salted ponds. Consequently, the denitrifying haloarchaea that sustain the nitrogen cycle under these conditions have emerged as an important target for research aimed at understanding microbial life in these extreme environments.

The haloarchaeon Haloferax mediterranei was isolated 20 years ago from Santa Pola salted ponds (Alicante, Spain). It was described as a denitrifier and it is also able to grow using NO₃^-, NO₂^- or NH₄⁺ as inorganic nitrogen sources. This review summarizes the advances that have been made in understanding the N-cycle in halophilic archaea using Hfx mediterranei as a haloarchaeal model. The results obtained show that this microorganism could be very attractive for bioremediation applications in those areas where high salt, nitrate and nitrite concentrations are found in ground waters and soils.

The Very Early Stages of Biological Evolution and the Nature of the Last Common Ancestor of the Three Major Cell Domains

Quantitative estimates of the gene complement of the last common ancestor of all extant organisms, that is, the cenancestor, may be hindered by ancient horizontal gene transfer events and polyphyletic gene losses, as well as by biases in genome databases and methodological artifacts. Nevertheless, most reports agree that the last common ancestor resembled extant prokaryotes. A significant number of the highly conserved genes are sequences involved in the synthesis, degradation, and binding of RNA, including transcription and translation. Although the gene complement of the cenancestor includes sequences that may have originated in different epochs, the extraordinary conservation of RNA-related sequences supports the hypothesis that the last common ancestor was an evolutionary outcome of the so-called RNA/protein world. The available evidence suggests that the cenancestor was not a hyperthermophile, but it is currently not possible to assess its ecological niche or its mode of energy acquisition and carbon sources.

Protein Superfamily Evolution and the Last Universal Common Ancestor (LUCA)

By exploiting three-dimensional structure comparison, which is more sensitive than conventional sequence-based methods for detecting remote homology, we have identified a set of 140 ancestral protein domains using very restrictive criteria to minimize the potential error introduced by horizontal gene transfer. These domains are highly likely to have been present in the Last Universal Common Ancestor (LUCA) based on their universality in almost all of 114 completed prokaryotic (Bacteria and Archaea) and eukaryotic genomes. Functional analysis of these ancestral domains reveals a genetically complex LUCA with practically all the essential functional systems present in extant organisms, supporting the theory that life achieved its modern cellular status much before the main kingdom separation (Doolittle 2000). In addition, we have calculated different estimations of the genetic and functional versatility of all the superfamilies and functional groups in the prokaryote subsample. These estimations reveal that some ancestral superfamilies have been more versatile than others during evolution allowing more genetic and functional variation. Furthermore, the differences in genetic versatility between protein families are more attributable to their functional nature rather than the time that they have been evolving. These differences in tolerance to mutation suggest that some protein families have eroded their phylogenetic signal faster than others, hiding in many cases, their ancestral origin and suggesting that the calculation of 140 ancestral domains is probably an underestimate.

The Rieske Protein: A Case Study on the Pitfalls of Multiple Sequence Alignments and Phylogenetic Reconstruction

Previously published phylogenetic trees reconstructed on "Rieske protein" sequences frequently are at odds with each other, with those of other subunits of the parent enzymes and with small-subunit rRNA trees. These differences are shown to be at least partially if not completely due to problems in the reconstruction procedures. A major source of erroneous Rieske protein trees lies in the presence of a large, poorly conserved domain prone to accommodate very long insertions in well-defined structural hot spots substantially hampering multiple alignments. The remaining smaller domain, in contrast, is too conserved to allow distant phylogenies to be deduced with sufficient confidence. Three-dimensional structures of representatives from this protein family are now available from phylogenetically distant species and from diverse enzymes. Multiple alignments can thus be refined on the basis of these structures. We show that structurally guided alignments of Rieske proteins from Rieske-cytochrome b complexes and arsenite oxidases strongly reduce conflicts between resulting trees and those obtained on their companion enzyme subunits. Further problems encountered during this work, mainly consisting in database errors such as wrong annotations and frameshifts, are described. The obtained results are discussed against the background of hypotheses stipulating pervasive lateral gene transfer in prokaryotes.

A minimal estimate for the gene content of the last universal common ancestor--exobiology from a terrestrial perspective

Using an algorithm for ancestral state inference of gene content, given a large number of extant genome sequences and a phylogenetic tree, we aim to reconstruct the gene content of the last universal common ancestor (LUCA), a hypothetical life form that presumably was the progenitor of the three domains of life. The method allows for gene loss, previously found to be a major factor in shaping gene content, and thus the estimate of LUCA's gene content appears to be substantially higher than that proposed previously, with a typical number of over 1000 gene families, of which more than 90% are also functionally characterized. More precisely, when only prokaryotes are considered, the number varies between 1006 and 1189 gene families while when eukaryotes are also included, this number increases to between 1344 and 1529 families depending on the underlying phylogenetic tree. Therefore, the common belief that the hypothetical genome of LUCA should resemble those of the smallest extant genomes of obligate parasites is not supported by recent advances in computational genomics. Instead, a fairly complex genome similar to those of free-living prokaryotes, with a variety of functional capabilities including metabolic transformation, information processing, membrane/transport proteins and complex regulation, shared between the three domains of life, emerges as the most likely progenitor of life on Earth, with profound repercussions for planetary exploration and exobiology.

Genome trees from conservation profiles

The concept of the genome tree depends on the potential evolutionary significance in the clustering of species according to similarities in the gene content of their genomes. In this respect, genome trees have often been identified with species trees. With the rapid expansion of genome sequence data it becomes of increasing importance to develop accurate methods for grasping global trends for the phylogenetic signals that mutually link the various genomes. We therefore derive here the methodological concept of genome trees based on protein conservation profiles in multiple species. The basic idea in this derivation is that the multi-component “presence-absence” protein conservation profiles permit tracking of common evolutionary histories of genes across multiple genomes. We show that a significant reduction in informational redundancy is achieved by considering only the subset of distinct conservation profiles. Beyond these basic ideas, we point out various pitfalls and limitations associated with the data handling, paving the way for further improvements. As an illustration for the methods, we analyze a genome tree based on the above principles, along with a series of other trees derived from the same data and based on pair-wise comparisons (ancestral duplication-conservation and shared orthologs). In all trees we observe a sharp discrimination between the three primary domains of life: Bacteria, Archaea, and Eukarya. The new genome tree, based on conservation profiles, displays a significant correspondence with classically recognized taxonomical groupings, along with a series of departures from such conventional clusterings.

Since Darwin's Origin of Species and Haeckel's Tree of Life, systematic biology has attempted to classify species into “family trees.” Genomics has provided a new framework permitting descriptions of sibling relations between species on the basis of their complete genetic blueprints. While trees based on single genes (rRNA), or limited numbers of genes have been useful, genome trees derived from complete genome comparisons should lead to more complete pictures of phylogenetic relations between various organisms. In order to reach such a global vision, procedures to establish sibling relationships should depend on an overall comparison that captures the evolutionary fates of proteins jointly in multiple genomes. This paper aims to establish a methodological basis to use genuine multidimensional procedures in the construction of genome trees. This approach completes the derivation of trees based on more classical techniques of pair-wise comparison between species. The authors survey classification schemes emerging from this approach, which either supports traditional views, such as the separation between the three phylogenetic domains Bacteria, Archaea, and Eukarya, or challenges them by suggesting, for example, intermingled clusterings of Proteobacteria with various other bacterial species.

The last common ancestor: what's in a name?

Twenty completely sequenced cellular genomes from the three major domains were analyzed using twice one-way BLAST searches in order to define the set of the most conserved protein-encoding sequences to characterize the gene complement of the last common ancestor of extant life. The resulting set is dominated by different putative ATPases, and by molecules involved in gene expression and RNA metabolism. DEAD-type RNA helicase and enolase genes, which are known to be part of the RNA degradosome, are as conserved as many transcription and translation genes. This suggests the early evolution of a control mechanism for gene expression at the RNA level, providing additional support to the hypothesis that during early cellular evolution RNA molecules played a more prominent role. Conserved sequences related to biosynthetic pathways include those encoding putative phosphoribosyl pyrophosphate synthase and thioredoxin, which participate in nucleotide metabolism. Although the information contained in the available databases corresponds only to a minor portion of biological diversity, the sequences reported here are likely to be part of an essential and highly conserved pool of proteins domains common to all organisms.

Analysis of phenetic trees based on metabolic capabilites across the three domains of life

Here, we used data of complete genomes to study comparatively the metabolism of different species. We built phenetic trees based on the enzymatic functions present in different parts of metabolism. Seven broad metabolic classes, comprising a total of 69 metabolic pathways, were comparatively analyzed for 27 fully sequenced organisms of the domains Eukarya, Bacteria and Archaea. Phylogenetic profiles based on the presence/absence of enzymatic functions for each metabolic class were determined and distance matrices for all the organisms were then derived from the profiles. Unrooted phenetic trees based upon the matrices revealed the distribution of the organisms according to their metabolic capabilities, reflecting the ecological pressures and adaptations that those species underwent during their evolution. We found that organisms that are closely related in phylogenetic terms could be distantly related metabolically and that the opposite is also true. For example, obligate bacterial pathogens were usually grouped together in our metabolic trees, demonstrating that obligate pathogens share common metabolic features regardless of their diverse phylogenetic origins. The branching order of proteobacteria often did not match their classical phylogenetic classification and Gram-positive bacteria showed diverse metabolic affinities. Archaea were found to be metabolically as distant from free-living bacteria as from eukaryotes, and sometimes were placed close to the metabolically highly specialized group of obligate bacterial pathogens. Metabolic trees represent an integrative approach for the comparison of the evolution of the metabolism and its correlation with the evolution of the genome, helping to find new relationships in the tree of life.

Lateral gene transfer and the origins of prokaryotic groups

Lateral gene transfer (LGT) is now known to be a major force in the evolution of prokaryotic genomes. To date, most analyses have focused on either (a) verifying phylogenies of individual genes thought to have been transferred, or (b) estimating the fraction of individual genomes likely to have been introduced by transfer. Neither approach does justice to the ability of LGT to effect massive and complex transformations in basic biology. In some cases, such transformation will be manifested as the patchy distribution of a seemingly fundamental property (such as aerobiosis or nitrogen fixation) among the members of a group classically defined by the sharing of other properties (metabolic, morphological, or molecular, such as small subunit ribosomal RNA sequence). In other cases, the lineage of recipients so transformed may be seen to comprise a new group of high taxonomic rank ("class" or even "phylum"). Here we review evidence for an important role of LGT in the evolution of photosynthesis, aerobic respiration, nitrogen fixation, sulfate reduction, methylotrophy, isoprenoid biosynthesis, quorum sensing, flotation (gas vesicles), thermophily, and halophily. Sometimes transfer of complex gene clusters may have been involved, whereas other times separate exchanges of many genes must be invoked.

Common evolutionary origin of mitochondrial and rickettsial respiratory chains

Comprehensive phylogenetic analysis of the subunits of respiratory chain was carried out using a variety of mitochondrial and bacterial sequences including those from all unfinished alpha-proteobacterial genomes known to date. Maximum likelihood, neighbor-joining, and maximum parsimony consensus trees, based on four proton-translocating complexes, placed mitochondria as a sister group to the order Rickettsiales of obligate endosymbiotic bacteria to the exclusion of free-living alpha-proteobacteria. Thus, phylogenetic relationship of most eukaryotic respiratory enzymes conforms to canonical pattern of mitochondrial ancestry, prior established in analyses of ribosomal RNAs, which are encoded by residual mitochondrial genomes. These data suggest that mitochondria may have derived from a reduced intracellular bacterium and that respiration may be the only evolutionary novelty brought into eukaryotes by mitochondrial endosymbiont.

Evolutionary dynamics of the DNA-binding domains in putative R2R3-MYB genes identified from rice subspecies indica and japonica genomes

The molecular evolution of the R2R3-MYB gene family is of great interest because it is one of the most important transcription factor gene families in the plant kingdom. Comparative analyses of a gene family may reveal important adaptive changes at the protein level and thereby provide insights that relate structure to function. We have performed a range of comparative and bioinformatics analyses on R2R3-MYB genes identified from the rice (Oryza sativa subsp. japonica and indica) and Arabidopsis genome sequences. The study provides an initial framework to investigate how different evolutionary lineages in a gene family evolve new functions. Our results reveal a remarkable excess of non-synonymous substitutions, an indication of adaptive selection on protein structure that occurred during the evolution of both helix1 and helix2 of rice R2R3-MYB DNA-binding domains. These flexible alpha-helix regions associated with high frequencies of excess non-synonymous substitutions may play critical roles in the characteristic packing of R2R3-MYB DNA-binding domains and thereby modify the protein-DNA interaction process resulting in the recognition of novel DNA-binding sites. Furthermore, a co-evolutionary pattern is found between the second alpha-helix of the R2 domain and the second alpha-helix of the R3 domain by examining all the possible alpha-helix pairings in both the R2 and R3 domains. This points to the functional importance of pairing interactions between related secondary structures.

Endosymbiosis and the design of eukaryotic electron transport

The bioenergetic organelles of eukaryotic cells, mitochondria and chloroplasts, are derived from endosymbiotic bacteria. Their electron transport chains (ETCs) resemble those of free-living bacteria, but were tailored for energy transformation within the host cell. Parallel evolutionary processes in mitochondria and chloroplasts include reductive as well as expansive events: On one hand, bacterial complexes were lost in eukaryotes with a concomitant loss of metabolic flexibility. On the other hand, new subunits have been added to the remaining bacterial complexes, new complexes have been introduced, and elaborate folding patterns of the thylakoid and mitochondrial inner membranes have emerged. Some bacterial pathways were reinvented independently by eukaryotes, such as parallel routes for quinol oxidation or the use of various anaerobic electron acceptors. Multicellular organization and ontogenetic cycles in eukaryotes gave rise to further modifications of the bioenergetic organelles. Besides mitochondria and chloroplasts, eukaryotes have ETCs in other membranes, such as the plasma membrane (PM) redox system, or the cytochrome P450 (CYP) system. These systems have fewer complexes and simpler branching patterns than those in energy-transforming organelles, and they are often adapted to non-bioenergetic functions such as detoxification or cellular defense.

Diversity, taxonomy and evolution of medium-chain dehydrogenase/reductase superfamily

A comprehensive, structural and functional, in silico analysis of the medium-chain dehydrogenase/reductase (MDR) superfamily, including 583 proteins, was carried out by use of extensive database mining and the blastp program in an iterative manner to identify all known members of the superfamily. Based on phylogenetic, sequence, and functional similarities, the protein members of the MDR superfamily were classified into three different taxonomic categories: (a) subfamilies, consisting of a closed group containing a set of ideally orthologous proteins that perform the same function; (b) families, each comprising a cluster of monophyletic subfamilies that possess significant sequence identity among them and might share or not common substrates or mechanisms of reaction; and (c) macrofamilies, each comprising a cluster of monophyletic protein families with protein members from the three domains of life, which includes at least one subfamily member that displays activity related to a very ancient metabolic pathway. In this context, a superfamily is a group of homologous protein families (and/or macrofamilies) with monophyletic origin that shares at least a barely detectable sequence similarity, but showing the same 3D fold. The MDR superfamily encloses three macrofamilies, with eight families and 49 subfamilies. These subfamilies exhibit great functional diversity including noncatalytic members with different subcellular, phylogenetic, and species distributions. This results from constant enzymogenesis and proteinogenesis within each kingdom, and highlights the huge plasticity that MDR superfamily members possess. Thus, through evolution a great number of taxa-specific new functions were acquired by MDRs. The generation of new functions fulfilled by proteins, can be considered as the essence of protein evolution. The mechanisms of protein evolution inside MDR are not constrained to conserve substrate specificity and/or chemistry of catalysis. In consequence, MDR functional diversity is more complex than sequence diversity. MDR is a very ancient protein superfamily that existed in the last universal common ancestor. It had at least two (and probably three) different ancestral activities related to formaldehyde metabolism and alcoholic fermentation. Eukaryotic members of this superfamily are more related to bacterial than to archaeal members; horizontal gene transfer among the domains of life appears to be a rare event in modern organisms.

The redox protein construction kit: pre-last universal common ancestor evolution of energy-conserving enzymes

Genome analyses and the resolution of three-dimensional structures have provided evidence in recent years for hitherto unexpected family relationships between redox proteins of very diverse enzymes involved in bioenergetic electron transport. Many of these enzymes appear in fact to be constructed from only a limited set of building blocks. Phylogenetic analysis of selected units from this "redox enzyme construction kit" indicates an origin for several prominent bioenergetic enzymes that is very early, lying before the divergence of Bacteria and Archaea. Possible scenarios for the early evolution of selected complexes are proposed based on the obtained tree topologies.