The universal ancestor was a thermophile or a hyperthermophile

Revisiting the Relationships Between Genomic G + C Content, RNA Secondary Structures, and Optimal Growth Temperature

Over twenty years ago Galtier and Lobry published a manuscript entitled "Relationships between Genomic G + C Content, RNA Secondary Structure, and Optimal Growth Temperature" in the Journal of Molecular Evolution that showcased the lack of a relationship between genomic G + C content and optimal growth temperature (OGT) in a set of about 200 prokaryotes. Galtier and Lobry also assessed the relationship between RNA secondary structures (rRNA stems, tRNAs) and OGT, and in this case a clear relationship emerged. Increasing structured RNA G + C content (particularly in regions that are double-stranded) correlates with increased OGT. Both of these fundamental relationships have withstood test of many additional sequences and spawned a variety of different applications that include prediction of OGT from rRNA sequence and computational ncRNA identification approaches. In this work, I present the motivation behind Galtier and Lobry's original paper and the larger questions addressed by the work, how these questions have evolved over the last two decades, and the impact of Galtier and Lobry's manuscript in fields beyond these questions.

The Role of Meteorite Impacts in the Origin of Life

The conditions, timing, and setting for the origin of life on Earth and whether life exists elsewhere in our solar system and beyond represent some of the most fundamental scientific questions of our time. Although the bombardment of planets and satellites by asteroids and comets has long been viewed as a destructive process that would have presented a barrier to the emergence of life and frustrated or extinguished life, we provide a comprehensive synthesis of data and observations on the beneficial role of impacts in a wide range of prebiotic and biological processes. In the context of previously proposed environments for the origin of life on Earth, we discuss how meteorite impacts can generate both subaerial and submarine hydrothermal vents, abundant hydrothermal-sedimentary settings, and impact analogues for volcanic pumice rafts and splash pools. Impact events can also deliver and/or generate many of the necessary chemical ingredients for life and catalytic substrates such as clays as well. The role that impact cratering plays in fracturing planetary crusts and its effects on deep subsurface habitats for life are also discussed. In summary, we propose that meteorite impact events are a fundamental geobiological process in planetary evolution that played an important role in the origin of life on Earth. We conclude with the recommendation that impact craters should be considered prime sites in the search for evidence of past life on Mars. Furthermore, unlike other geological processes such as volcanism or plate tectonics, impact cratering is ubiquitous on planetary bodies throughout the Universe and is independent of size, composition, and distance from the host star. Impact events thus provide a mechanism with the potential to generate habitable planets, moons, and asteroids throughout the Solar System and beyond.

Computational Identification of Amino-Acid Mutations that Further Improve the Activity of a Chalcone-Flavonone Isomerase from Glycine max

Protein design for improving enzymatic activity remains a challenge in biochemistry, especially to identify target amino-acid sites for mutagenesis and to design beneficial mutations for those sites. Here, we employ a computational approach that combines multiple sequence alignment, positive selection detection, and molecular docking to identify and design beneficial amino-acid mutations that further improve the intramolecular-cyclization activity of a chalcone-flavonone isomerase from Glycine max (GmCHI). By this approach, two GmCHI mutants with higher activities were predicted and verified. The results demonstrate that this approach could determine the beneficial amino-acid mutations for improving the enzymatic activity, and may find more applications in engineering of enzymes.

An Interpretation of the Ancestral Codon from Miller's Amino Acids and Nucleotide Correlations in Modern Coding Sequences

Purine bias, which is usually referred to as an “ancestral codon”, is known to result in short-range correlations between nucleotides in coding sequences, and it is common in all species. We demonstrate that RWY is a more appropriate pattern than the classical RNY, and purine bias (Rrr) is the product of a network of nucleotide compensations induced by functional constraints on the physicochemical properties of proteins. Through deductions from universal correlation properties, we also demonstrate that amino acids from Miller’s spark discharge experiment are compatible with functional primeval proteins at the dawn of living cell radiation on earth. These amino acids match the hydropathy and secondary structures of modern proteins.

The origin of the genetic code in the ocean abysses: new comparisons confirm old observations

I have analysed the amino acid substitution pattern between two pairs of nonbarophilic-barophilic organisms in order to confirm previous results. Indeed, the pattern deriving from a different pair of such organisms led to establish that the origin of the genetic code might have occurred in the ocean abysses. The hydrostatic pressure asymmetry indices computable from these matrices of amino acid substitutions confirm the correlation previously observed, even when differences in GC content are accounted for. As the three substitution matrices are independent between them, the old conclusion that the genetic code originated in the ocean abysses is considerably strengthened.

An archaeal NADH oxidase causes damage to both proteins and nucleic acids under oxidative stress

NADH oxidases (NOXs) catalyze the two-electron reduction of oxygen to H2O2 or four-electron reduction of oxygen to H2O. In this report, we show that an NADH oxidase from Thermococcus profundus (NOXtp) displays two forms: a native dimeric protein under physiological conditions and an oxidized hexameric form under oxidative stress. Native NOXtp displays high NADH oxidase activity, and oxidized NOXtp can accelerate the aggregation of partially unfolded proteins. The aggregates formed by NOXtp have characteristics similar to beta-amyloid and Lewy bodies in neurodegenerative diseases, including an increase of beta-sheet content. Oxidized NOXtp can also bind nucleic acids and cause their degradation by oxidizing NADH to produce H2O2. Furthermore, Escherichia coli cells expressing NOXtp are less viable than cells not expressing NOXtp after treatment with H2O2. As NOXtp shares similar features with eukaryotic cell death isozymes and life may have originated from hyperthermophiles, we suggest that NOXtp may be an ancestor of cell death proteins.

Peptide synthesis in early Earth hydrothermal systems

We report here results from experiments and thermodynamic calculations that demonstrate a rapid, temperature-enhanced synthesis of oligopeptides from the condensation of aqueous glycine. Experiments were conducted in custom-made hydrothermal reactors, and organic compounds were characterized with ultraviolet-visible procedures. A comparison of peptide yields at 260 degrees C with those obtained at more moderate temperatures (160 degrees C) gives evidence of a significant (13 kJ . mol(-1)) exergonic shift. In contrast to previous hydrothermal studies, we demonstrate that peptide synthesis is favored in hydrothermal fluids and that rates of peptide hydrolysis are controlled by the stability of the parent amino acid, with a critical dependence on reactor surface composition. From our study, we predict that rapid recycling of product peptides from cool into near-supercritical fluids in mid-ocean ridge hydrothermal systems will enhance peptide chain elongation. It is anticipated that the abundant hydrothermal systems on early Earth could have provided a substantial source of biomolecules required for the origin of life.

Kinetics and activation parameter analysis for the prebiotic oligocytidylate formation on Na(+)-montmorillonite at 0-100 degrees C

The kinetic analysis of the temperature dependence of the formation of oligocytidylate (oligo(C)) from the 5'-monophosphorimidazolide moiety of cytidine (ImpC) in the presence of Na (+)-montmorillonite (Na (+)-Mont) catalyst has been carried out at 0-100 degrees C. The rate constants for the formation of oligo(C), hydrolysis of ImpC with and without Na (+)-Mont and degradation of oligo(C) were determined. The apparent activation parameters were 30.8 +/- 3.9 kJ mol (-1) ( Ea), 28.3 +/- 4.0 kJ mol (-1) (Delta H++), and -231 +/- 13 J mol (-1) K (-1) (Delta S++) for the formation of the 2-mer; 45.6 +/- 2.9 kJ mol (-1) ( Ea), 43.0 +/- 3.0 kJ mol (-1) (Delta H++), -164 +/- 10 J mol (-1) K (-1) (Delta S++) for the 3-mer; and 45.2 +/- 0.6 kJ mol (-1) ( Ea), 42.7 +/- 0.7 kJ mol (-1) (Delta H++), -159 +/- 2 J mol (-1) K (-1) (Delta S++) for the 4-mer in the presence of Na (+)-Mont. An increasing trend for the rate constants for the formation of oligo(C) in the order 2-mer << 3-mer <4-mer was observed at high temperatures, which is consistent with that observed at low temperatures. These analyses implied for the first time that the associate formation between an activated nucleotide monomer and an elongating oligonucleotide prior to the phosphodiester bond formation during the elongation of an oligonucleotide on a clay surface would be based on the interaction between the two reactants at the phosphoester and/or ribose moieties rather than at the nucleotide bases. The hydrolysis rate of ImpC at 25-100 degrees C was 5.3-10.6 times greater in the presence of Na (+)-Mont than in its absence. Although the degradation of oligo(C) in the presence of Na (+)-Mont was slower than the formation of the 3-mer and longer oligo(C) on Na (+)-Mont, its yield decreased with temperature. This is mainly because the ratios of the rate constant of the 2-mer formation to those of ImpC hydrolysis and the 3-mer and 4-mer formation decrease with an increase in temperature, which is attributed to the enthalpy and entropy changes for the formation of the 2-mer. This trend resembles the case of the template-directed formation of oligo(G) on a poly(C) template but is different from the Pb (2+)-ion-catalyzed oligo(C) formation. According to the kinetics and activation parameter analyses regarding the clay reaction and other prebiotic polymerase models, the possible pathways for the oligonucleotide formation are discussed and compared.

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.

Intense neutral drifts yield robust and evolvable consensus proteins

What changes occur when a natural protein that had been under low mutation rates is subjected to a neutral drift at high mutational loads, thus generating genetically diverse (polymorphic) gene ensembles that all maintain the protein's original function and structure? To address this question we subjected large populations of TEM-1 beta-lactamase to a prolonged neutral drift, applying high mutation rates and purifying selection to maintain TEM-1's existing penicillinase activity. Purging of deleterious mutations and enrichment of beneficial ones maintained the sequence of these ensembles closer to TEM-1's family consensus and inferred ancestor. In particular, back-to-consensus/ancestor mutations that increase TEM-1's kinetic and thermodynamic stability were enriched. These acted as global suppressors and enabled the tolerance of a broad range of deleterious mutations, thus further increasing the genetic diversity of the drifting populations. The probability of a new function emerging (cefotaxime degradation) was also substantially increased in these ensembles owing to the presence of many gene variants carrying the global suppressors. Our findings indicate the unique features of large, polymorphic neutral ensembles generated under high mutational loads and prompt the speculation that the progenitors of today's proteins may have evolved under high mutational loads. The results also suggest that predictable back-to-consensus/ancestor changes can be used in the laboratory to generate highly diverse and evolvable gene libraries.

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.

Extremely thermophilic translation system in the common ancestor commonote: ancestral mutants of Glycyl-tRNA synthetase from the extreme thermophile Thermus thermophilus

Based on phylogenetic analysis of 16 S and 18 S rRNAs, the common ancestor of all organisms (Commonote) was proposed to be hyperthermophilic. We have previously tested this hypothesis using enzymes with ancestral residues that are inferred by molecular phylogenetic analysis. The ancestral mutant enzymes involved in metabolic systems show higher thermal stability than wild-type enzymes, consistent with the hyperthermophile common ancestor hypothesis. Here, we have extended the experiments to include an enzyme of the translation system, glycyl-tRNA synthetase (GlyRS). The translation system often shows a phylogenetic tree that is similar to the rRNA tree. Thus, it is likely that the tree represents the evolutionary route of the organisms. The maximum-likelihood tree of alpha(2) type GlyRS was constructed. From this analysis the ancestral sequence of GlyRS was deduced and individual or pairs of ancestral residues were introduced into Thermus thermophilus GlyRS. The ancestral mutants were expressed in Escherichia coli, purified and activity measured. The thermostability of eight mutated proteins was evaluated by CD (circular dichroism) measurements. Six mutants showed higher thermostability than wild-type enzyme and seven mutants showed higher activity than wild-type enzyme at 70 degrees C, suggesting an extremely thermophilic translation system in the common ancestor Commonote.

The universal ancestor and the ancestors of Archaea and Bacteria were anaerobes whereas the ancestor of the Eukarya domain was an aerobe

The use of an oxyphobic index (OI) based on the propensity of amino acids to enter more frequently the proteins of anaerobes makes it possible to make inferences on the environment in which the last universal common ancestor (LUCA) lived. The reconstruction of the ancestral sequences of proteins using a method based on maximum likelihood and their attribution by means of the OI to the set of aerobe or anaerobe sequences has led to the following conclusions: the LUCA was an anaerobic 'organism', as were the ancestors of Archaea and Bacteria, whereas the ancestor of Eukarya was an aerobe. These observations seem to falsify the hypothesis that the LUCA was an aerobe and help to identify better the environment in which the first organisms lived.

Structuring of the genetic code took place at acidic pH

I have observed that in multiple regression the number of codons specifying amino acids in the genetic code is positively correlated with the isoelectric point of amino acids and their molecular weight. Therefore basic amino acids are, on average, codified in the genetic code by a larger number of codons, which seems to imply that the genetic code originated in an acidic 'intracellular' environment. Moreover, I compare the proteins from Picrophilus torridus and Thermoplasma volcanium, which have different intracellular pH and I define the ranks of acidophily for the amino acids. A simple index of acidophily (AI), which can be easily obtained from acidophily ranks, can be associated to any protein and, therefore, can also be associated to the genetic code if the number of synonymous codons attributed to the amino acids in the code is assumed to be the frequency with which the amino acids appeared in ancestral proteins. Finally, the sampling of the variable AI among organisms having an intracellular pH less than or equal to 6.6 and those having a non-acidic intracellular pH leads to the conclusion that the value of the genetic code's AI is not typical of proteins of the latter organisms. As the genetic code's AI value is also statistically not different from that of proteins of the organisms having an acidic intracellular pH, this supports the hypothesis that the structuring of the genetic code took place in acidic pH conditions.

A comparison of proteins from Pyrococcus furiosus and Pyrococcus abyssi: barophily in the physicochemical properties of amino acids and in the genetic code

A comparison is made between orthologous proteins from a nonbarophile (Pyrococcus furiosus) and a barophile (Pyrococcus abyssi) organism. The pattern of asymmetries in the amino acid substitution process identifies the amino acids arginine, serine, glycine, valine and aspartic acid as those having the most barophilic behaviour, and tyrosine and glutamine as the least barophilic. The construction of a hydrostatic pressure asymmetry index (PAI) which orders the amino acids from the most barophilic to the least barophilic makes it possible to visualise the amino acid properties that best explain barophily. The polarity of amino acids is positively correlated to the PAI values, i.e., on average, the more polar amino acids possess a higher PAI value, that is to say they are more barophilic. Moreover, the "size" of amino acids (molecular weight) is negatively correlated to the PAI value, that is to say that, on average, high PAI values are associated to "small" amino acids which are therefore more barophilic than "larger" ones. These two amino acid properties are the same ones that are known for having been important in affecting the origin of genetic code organisation. All the above, as well as the significant and positive correlation between the number of codons attributed to the amino acids in the genetic code and the PAI values, seem to favour the hypothesis that genetic code structuring took place under high hydrostatic pressure.

The ocean abysses witnessed the origin of the genetic code

The comparison of proteins from a non-barophilous and a barophilous organism makes it possible to define the barophily ranks of amino acids. The correlation of these ranks with the number of codons attributed to amino acids in the genetic code, together with another straightforward argument based on an optimisation percentage of a barophily index (BI) (easily defined by barophily ranks) which can be associated to the genetic code table, suggest that the genetic code originated under high hydrostatic pressure. Moreover, as the BI value can be calculated for the sequence of any protein, it also makes it possible to define the BI for the genetic code if the number of codons attributed to the amino acids in the code is assumed to be the frequency with which the amino acids appeared in ancestral proteins. Finally, sampling the BI variable between many non-barophile organisms and from many proteins of a single non-barophile organism leads to the conclusion that the BI value of the genetic code is not typical of these organisms. Whereas, since the genetic code BI value is statistically higher than that of these non-barophile organisms, it supports the hypothesis that genetic code structuring took place under high hydrostatic pressure.

Correlations between genomic GC levels and optimal growth temperatures are not ‘robust’

Musto et al. [FEBS Lett. 573 (2004) 73] studied the correlations between GC levels and optimal growth temperatures in 20 prokaryotic families. They reported that positive correlations are generally observed, and many of these are significant. Here, we have shown that these correlations are not "robust," i.e., correlation coefficients and/or significance of correlations can be considerably influenced by exclusion of very few (even as small as one) species from each dataset. The sensitivity of correlations is assumed as a result of high levels of bias in the family datasets. We concluded that solely based on these data, one cannot establish that GC contents of prokaryotic genomes increase as a result of growth temperature increments.

LUCA, or just a conserved Archaeon?: Comments on Xue et al. (2003)

In their recent paper, Xue et al. used an unusual technique of rooting the universal phylogenetic tree, which resulted in positioning of the last universal common ancestor within Archaea. The present paper brings some criticisms on the methods and results achieved.

The ancestor of the Bacteria domain was a hyperthermophile

Brochier and Philippe have recently re-analysed the phylogeny of ribosomal RNA using only multiple alignment positions with no phylogenetic noise. They conclude that the first branch of divergence in the Bacteria domain comprises Planctomycetales and not hyperthermophile bacteria as in classic phylogeny. In the present paper I examine the robustness of their conclusions. (1) A site-by-site reading of the RNA alignments of Brochier and Philippe seems to suggest that the number of nucleotide positions used in their analysis is not sufficiently high and their phylogenetic analysis is consequently not robust. Furthermore, (2) a different method for selecting positions with no phylogenetic noise from the rRNA alignment relocates the Aquificales and the Thermotogales as the first lines of divergence in the Bacteria domain, and sets Planctomycetales as the third branch of divergence in the phylogenetic tree built from these selected positions. These findings consolidate the hypothesis that the ancestor of the Bacteria domain was a hyperthermophile and, more generally, that the last universal common ancestor might also have been one.

Transfer RNA paralogs: evidence for genetic code-amino acid biosynthesis coevolution and an archaeal root of life

A search has been performed on 2878 tRNA sequences from 60 different genomes in order to detect the existence of closely related 'alloacceptor' tRNAs accepting dissimilar amino acids that could be paralogs generated by gene duplications. This has led to the identification of extremely conserved tRNA(Phe)-tRNA(Tyr) pairs displaying as high as 94% identity between them, and also other potentially paralogous tRNA pairs in archaeal species. These paralogous pairs are enriched for amino acid pairs belonging to the same amino acid biosynthetic family, thus providing evidence for the coevolution of genetic code and amino acid biosynthesis. Overall, the genetic distances between alloacceptor tRNAs yield estimates of how closely clustered in sequence space are the tRNAs in a genome. Among 34 Bacteria, 18 Archaea and 8 Eukarya, Methanopyrus kandleri and Aeropyrum pernix have yielded the lowest alloacceptor distances and largest number of paralogous pairs. Based on a cluster-dispersion model of tRNA evolution, such tight alloacceptor clustering is a measure of primitiveness of tRNA genotypes, and places last universal common ancestor (LUCA) between the branches leading to these two archaea in the tRNA phylogenetic tree.

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.

Early assembly of cellular life

Popular hypotheses that attempt to explain the origin of prebiotic molecules and cellular life capable of growth and division are not always agreed upon. In this manuscript, information on early bacterial life on Earth is examined using information from several disciplines. For example, knowledge can be integrated from physics, thermodynamics, planetary sciences, geology, biogeochemistry, lipid chemistry, primordial cell structures, cell and molecular biology, microbiology, metabolism and genetics. The origin of life also required a combination of elements, compounds and environmental physical-chemical conditions that allowed cells to assemble in less than a billion years. This may have been widespread in the subsurface of the early Earth located at microscopic physical domains.

Algorithms for computing parsimonious evolutionary scenarios for genome evolution, the last universal common ancestor and dominance of horizontal gene transfer in the evolution of prokaryotes

BACKGROUND

Comparative analysis of sequenced genomes reveals numerous instances of apparent horizontal gene transfer (HGT), at least in prokaryotes, and indicates that lineage-specific gene loss might have been even more common in evolution. This complicates the notion of a species tree, which needs to be re-interpreted as a prevailing evolutionary trend, rather than the full depiction of evolution, and makes reconstruction of ancestral genomes a non-trivial task.

RESULTS

We addressed the problem of constructing parsimonious scenarios for individual sets of orthologous genes given a species tree. The orthologous sets were taken from the database of Clusters of Orthologous Groups of proteins (COGs). We show that the phyletic patterns (patterns of presence-absence in completely sequenced genomes) of almost 90% of the COGs are inconsistent with the hypothetical species tree. Algorithms were developed to reconcile the phyletic patterns with the species tree by postulating gene loss, COG emergence and HGT (the latter two classes of events were collectively treated as gene gains). We prove that each of these algorithms produces a parsimonious evolutionary scenario, which can be represented as mapping of loss and gain events on the species tree. The distribution of the evolutionary events among the tree nodes substantially depends on the underlying assumptions of the reconciliation algorithm, e.g. whether or not independent gene gains (gain after loss after gain) are permitted. Biological considerations suggest that, on average, gene loss might be a more likely event than gene gain. Therefore different gain penalties were used and the resulting series of reconstructed gene sets for the last universal common ancestor (LUCA) of the extant life forms were analysed. The number of genes in the reconstructed LUCA gene sets grows as the gain penalty increases. However, qualitative examination of the LUCA versions reconstructed with different gain penalties indicates that, even with a gain penalty of 1 (equal weights assigned to a gain and a loss), the set of 572 genes assigned to LUCA might be nearly sufficient to sustain a functioning organism. Under this gain penalty value, the numbers of horizontal gene transfer and gene loss events are nearly identical. This result holds true for two alternative topologies of the species tree and even under random shuffling of the tree. Therefore, the results seem to be compatible with approximately equal likelihoods of HGT and gene loss in the evolution of prokaryotes.

CONCLUSIONS

The notion that gene loss and HGT are major aspects of prokaryotic evolution was supported by quantitative analysis of the mapping of the phyletic patterns of COGs onto a hypothetical species tree. Algorithms were developed for constructing parsimonious evolutionary scenarios, which include gene loss and gain events, for orthologous gene sets, given a species tree. This analysis shows, contrary to expectations, that the number of predicted HGT events that occurred during the evolution of prokaryotes might be approximately the same as the number of gene losses. The approach to the reconstruction of evolutionary scenarios employed here is conservative with regard to the detection of HGT because only patterns of gene presence-absence in sequenced genomes are taken into account. In reality, horizontal transfer might have contributed to the evolution of many other genes also, which makes it a dominant force in prokaryotic evolution.