On the Origin of Metabolic Pathways

The heterotrophic theory of the origin of life is the only proposal available with experimental support. This comes from the ease of prebiotic synthesis under strongly reducing conditions. The prebiotic synthesis of organic compounds by reduction of CO₂ to monomers used by the first organisms would also be considered an heterotrophic origin. Autotrophy means that the first organisms biosynthesized their cell constituents as well as assembling them. Prebiotic synthetic pathways are all different from the biosynthetic pathways of the last common ancestor (LCA). The steps leading to the origin of the metabolic pathways are closer to prebiotic chemistry than to those in the LCA. There may have been different biosynthetic routes between the prebiotic and the LCAs that played an early role in metabolism but have disappeared from extant organisms. The semienzymatic theory of the origin of metabolism proposed here is similar to the Horowitz hypothesis but includes the use of compounds leaking from preexisting pathways as well as prebiotic compounds from the environment.     

The Path from the RNA World

We describe a sequential (step by step) Darwinian model for the evolution of life from the late stages of the RNA world through to the emergence of eukaryotes and prokaryotes. The starting point is our model, derived from current RNA activity, of the RNA world just prior to the advent of genetically-encoded protein synthesis. By focusing on the function of the protoribosome we develop a plausible model for the evolution of a protein-synthesizing ribosome from a high-fidelity RNA polymerase that incorporated triplets of oligonucleotides. With the standard assumption that during the evolution of enzymatic activity, catalysis is transferred from RNA → RNP → protein, the first proteins in the ``breakthrough organism' (the first to have encoded protein synthesis) would be nonspecific chaperone-like proteins rather than catalytic. Moreover, because some RNA molecules that pre-date protein synthesis under this model now occur as introns in some of the very earliest proteins, the model predicts these particular introns are older than the exons surrounding them, the ``introns-first' theory. Many features of the model for the genome organization in the final RNA world ribo-organism are more prevalent in the eukaryotic genome and we suggest that the prokaryotic genome organization (a single, circular genome with one center of replication) was derived from a ``eukaryotic-like' genome organization (a fragmented linear genome with multiple centers of replication). The steps from the proposed ribo-organism RNA genome → eukaryotic-like DNA genome → prokaryotic-like DNA genome are all relatively straightforward, whereas the transition prokaryotic-like genome → eukaryotic-like genome appears impossible under a Darwinian mechanism of evolution, given the assumption of the transition RNA → RNP → protein. A likely molecular mechanism, ``plasmid transfer,' is available for the origin of prokaryotic-type genomes from an eukaryotic-like architecture. Under this model prokaryotes are considered specialized and derived with reduced dependence on ssRNA biochemistry. A functional explanation is that prokaryote ancestors underwent selection for thermophily (high temperature) and/or for rapid reproduction (r selection) at least once in their history. Received: 14 January 1997 / Accepted: 19 May 1997     

The Prebiotic Synthesis of Modified Purines and Their Potential Role in the RNA World

Modified purines are found in all organisms in the tRNA, rRNA, and even DNA, raising the possibility of an early role for these compounds in the evolution of life. These include N ⁶-methyladenine, 1-methyladenine, N ⁶,N ⁶-dimethyladenine, 1-methylhypoxanthine, 1-methylguanine, and N ²-methylguanine. We find that these bases as well as a number of nonbiological modified purines can be synthesized from adenine and guanine by the simple reaction of an amine or an amino group with adenine and guanine under the concentrated conditions of the drying-lagoon or drying-beach model of prebiotic synthesis with yields as high as 50%. These compounds are therefore as prebiotic as adenine and guanine and could have played an important role in the RNA world by providing additional functional groups in ribozymes, especially for the construction of hydrophobic binding pockets. Received: 7 August 1998 / Accepted: 31 December 1998     

Relics from the RNA World

An RNA world is widely accepted as a probable stage in the early evolution of life. Two implications are that proteins have gradually replaced RNA as the main biological catalysts and that RNA has not taken on any major de novo catalytic function after the evolution of protein synthesis, that is, there is an essentially irreversible series of steps RNA → RNP → protein. This transition, as expected from a consideration of catalytic perfection, is essentially complete for reactions when the substrates are small molecules. Based on these principles we derive criteria for identifying RNAs in modern organisms that are relics from the RNA world and then examine the function and phylogenetic distribution of RNA for such remnants of the RNA world. This allows an estimate of the minimum complexity of the last ribo-organism—the stage just preceding the advent of genetically encoded protein synthesis. Despite the constraints placed on its size by a low fidelity of replication (the Eigen limit), we conclude that the genome of this organism reached a considerable level of complexity that included several RNA-processing steps. It would include a large protoribosome with many smaller RNAs involved in its assembly, pre-tRNAs and tRNA processing, an ability for recombination of RNA, some RNA editing, an ability to copy to the end of each RNA strand, and some transport functions. It is harder to recognize specific metabolic reactions that must have existed but synthetic and bio-energetic functions would be necessary. Overall, this requires that such an organism maintained a multiple copy, double-stranded linear RNA genome capable of recombination and splicing. The genome was most likely fragmented, allowing each ``chromosome' to be replicated with minimum error, that is, within the Eigen limit. The model as developed serves as an outgroup to root the tree of life and is an alternative to using sequence data for inferring properties of the earliest cells. Received: 14 January 1997 / Accepted: 19 May 1997     

The Genetic Code Is One in a Million

Statistical and biochemical studies of the genetic code have found evidence of nonrandom patterns in the distribution of codon assignments. It has, for example, been shown that the code minimizes the effects of point mutation or mistranslation: erroneous codons are either synonymous or code for an amino acid with chemical properties very similar to those of the one that would have been present had the error not occurred. This work has suggested that the second base of codons is less efficient in this respect, by about three orders of magnitude, than the first and third bases. These results are based on the assumption that all forms of error at all bases are equally likely. We extend this work to investigate (1) the effect of weighting transition errors differently from transversion errors and (2) the effect of weighting each base differently, depending on reported mistranslation biases. We find that if the bias affects all codon positions equally, as might be expected were the code adapted to a mutational environment with transition/transversion bias, then any reasonable transition/transversion bias increases the relative efficiency of the second base by an order of magnitude. In addition, if we employ weightings to allow for biases in translation, then only 1 in every million random alternative codes generated is more efficient than the natural code. We thus conclude not only that the natural genetic code is extremely efficient at minimizing the effects of errors, but also that its structure reflects biases in these errors, as might be expected were the code the product of selection. Received: 25 July 1997 / Accepted: 9 January 1998     

Evolutionary History of 4.5SI RNA and Indication That It Is Functional

To date, the small nuclear 4.5S_I RNA has only been studied in the rat (Rattus norvegicus). Combining PCR and hybridization analyses, we have revealed 4.5S_I RNA homologues sequences in the genomes of four myomorph rodent families (Muridae, Cricetidae, Spalicidae, and Rhizomyidae), and not in other myomorph families (Dipodidae, Zapodidae, Geomyidae, and Heteromyidae) or sciuromorph and caviomorph rodents. By Northern-hybridization, 4.5S_I RNA has been detected in the common rat (R. norvegicus, Muridae), golden hamster (Mesocricetus auratus, Cricetidae), and Russian mole rat (Spalax microphthalmus, Spalacidae), but not in the related great jerboa (Allactaga jaculus, Dipodidae) or in four non-myomorph rodent species tested. cDNA derived from 4.5S_I RNA of M. auratus and S. microphthalmus has been cloned and sequenced. The hamster RNA is found to differ from rat 4.5S_I RNA by only one nucleotide substitution. For the mole rat, two variants of 4.5S_I RNA are detected: short (S) and long (L) with length 101 and 108 nt, respectively. The L variant differs from the S variant as well as from murid and cricetid 4.5S_I RNAs by both a 7 nt insertion and a varying number of nucleotide substitutions. The sequence similarity between the spalacid S-variant and murid/crecitid variants of 4.5S_I RNA is 90%. Judging from species distribution, 4.5S_I RNA genes emerged during the same period of time as the related short interspersed element B2 arose. This occurred after the divergence of Dipodidae lineage but before the branching of Spalicidae/Rhizomyidae lineage from a common myomorph rodent stem. S variant genes seemed to emerge in a common ancestor of spalacids and rhizomyds whereas L variant genes formed in spalacids following the divergence of these two families. The low rate of evolutionary changes of 4.5S_I RNA, at least, in murids and cricetids (6 × 10⁻⁴ substitutions per site per million years), suggests that this RNA is under selection constraint and have a function. This is a remarkable fact if the recent origin and narrow species distribution range of 4.5S_I RNA genes is taken into account. Genes with narrow species distribution are proposed to be referred to as stenogenes. Received: 11 December 2000 / Accepted: 27 August 2001     

The Length Distribution of Perfect Dimer Repetitive DNA Is Consistent with Its Evolution by an Unbiased Single-Step Mutation Process

We have examined the length distribution of perfect dimer repeats, where perfect means uninterrupted by any other base, using data from GenBank on primates and rodents. Virtually no lengths greater than 30 repeats are found, except for rodent AG repeats, which extend to 35. Comparable numbers of long AC and AG repeats suggest that they have not been selected for special functions or DNA structures. We have compared the data with predictions of two models: (1) a Bernoulli Model in which bases are assumed equally likely and distributed at random and (2) an Unbiased Random Walk Model (URWM) in which repeats are permitted to change length by plus or minus one unit, with equal probabilities, and in which base substitutions are allowed to destroy long perfect repeats, producing two shorter perfect repeats. The source of repeats is assumed to be from single base substutions from neighboring sequences, i.e., those differing from the perfect repeat by a single base. Mutation rates either independent of repeat length or proportional to length were considered. An upper limit to the lengths L≈ 30 is assumed and isolated dimers are assumed unable to expand, so that there are absorbing barriers to the random walk at lengths 1 and L+ 1, and a steady state of lengths is reached. With these assumptions and estimated values for the rates of length mutation and base substitution, reasonable agreement is found with the data for lengths > 5 repeats. Shorter repeats, of lengths ≤ 3 are in general agreement with the Bernoulli Model. By reducing the rate of length mutations for n≤ 5, it is possible to obtain reasonable agreement with the full range of data. For these reduced rates, the times between length mutations become comparable to those suggested for a bottleneck in the evolution of Homo sapiens, which may be the reason for low heterozygosity of short repeats.     

The Historical Factor: The Biosynthetic Relationships Between Amino Acids and Their Physicochemical Properties in the Origin of the Genetic Code

Two forces are in general, hypothesized to have influenced the origin of the organization of the genetic code: the physicochemical properties of amino acids and their biosynthetic relationships. In view of this, we have considered a model incorporating these two forces. In particular, we have studied the optimization level of the physicochemical properties of amino acids in the set of amino acid permutation codes that respects the biosynthetic relationships between amino acids. Where the properties of amino acids are represented by polarity and molecular volume we obtain indetermination percentages in the organization of the genetic code of approximately 40%. This indicates that the contingent factor played a significant role in structuring the genetic code. Furthermore, this result is in agreement with the genetic code coevolution hypothesis, which attributes a merely ancillary role to the properties of amino acids while it suggests that it was their biosynthetic relationships that organized the code. Furthermore, this result does not favor the stereochemical models proposed to explain the origin of the genetic code. On the other hand, where the properties of amino acids are represented by polarity alone, we obtain an indetermination percentage of at least 21.5%. This might suggest that the polarity distances played an important role and would therefore provide evidence in favor of the physicochemical hypothesis of genetic code origin. Although, overall, the analysis might have given stronger support to the latter hypothesis, this did not actually occur. The results are therefore discussed in the context of the different theories proposed to explain the origin of the genetic code. Received: 10 September 1996 / Accepted: 3 March 1997     

A Blind Empiricism Against the Coevolution Theory of the Origin of the Genetic Code

Ronneberg et al. (Proc Natl Acad Sci USA 97:13690–13695, 2000) recently suggested abandoning the coevolution theory of genetic code origin on the basis of two pieces of evidence. They (1) criticize the use of several pairs of amino acids in a precursor–product relationship to support this theory and (2) suggest a new set of codes in which to investigate the statistical bases of the coevolution theory, reaching the conclusion that this theory is not statistically validated in this set. In this paper I critically analyze the robustness of these conclusions. Observations and arguments lead to the belief that the pairs of amino acids in a precursor–product relationship originally used by the coevolution theory are such, or may at least be interpreted as such, and are therefore a manifestation of this theory. Furthermore, the new set of codes that Ronneberg et al. suggest is open to criticism and is thus substituted by the set of amino acid permutation codes, in which even the pairs of amino acids they favor end up by supporting the coevolution theory. Overall, the analysis seems to show that the paper by Ronneberg et al. is of minor scientific value while the coevolution theory seems to be one of the best theories at our disposal for explaining the evolutionary organisation of the genetic code and is, contrary to their claims, statistically well validated. Received: 21 February 2001 / Accepted: 22 May 2001     

Molecular Evolution of Nitrogen Fixation: The Evolutionary History of the nifD, nifK, nifE, and nifN Genes

The pairs of nitrogen fixation genes nifDK and nifEN encode for the α and β subunits of nitrogenase and for the two subunits of the NifNE protein complex, involved in the biosynthesis of the FeMo cofactor, respectively. Comparative analysis of the amino acid sequences of the four NifD, NifK, NifE, and NifN in several archaeal and bacterial diazotrophs showed extensive sequence similarity between them, suggesting that their encoding genes constitute a novel paralogous gene family. We propose a two-step model to reconstruct the possible evolutionary history of the four genes. Accordingly, an ancestor gene gave rise, by an in-tandem paralogous duplication event followed by divergence, to an ancestral bicistronic operon; the latter, in turn, underwent a paralogous operon duplication event followed by evolutionary divergence leading to the ancestors of the present-day nifDK and nifEN operons. Both these paralogous duplication events very likely predated the appearance of the last universal common ancestor. The possible role of the ancestral gene and operon in nitrogen fixation is also discussed. Received: 21 June 1999 / Accepted: 1 March 2000     

The puzzle of the Krebs citric acid cycle: Assembling the pieces of chemically feasible reactions,and opportunism in the design of metabolic pathways during evolution

The evolutionary origin of the Krebs citric acid cycle has been for a long time a model case in the understanding of the origin and evolution of metabolic pathways: How can the emergence of such a complex pathway be explained? A number of speculative studies have been carried out that have reached the conclusion that the Krebs cycle evolved from pathways for amino acid biosynthesis, but many important questions remain open: Why and how did the full pathway emerge from there? Are other alternative routes for the same purpose possible? Are they better or worse? Have they had any opportunity to be developed in cellular metabolism evolution? We have analyzed the Krebs cycle as a problem of chemical design to oxidize acetate yielding reduction equivalents to the respiratory chain to make ATP. Our analysis demonstrates that although there are several different chemical solutions to this problem, the design of this metabolic pathway as it occurs in living cells is the best chemical solution: It has the least possible number of steps and it also has the greatest ATP yielding. Study of the evolutionary possibilities of each one-taking the available material to build new pathways-demonstrates that the emergence of the Krebs cycle has been a typical case of opportunism in molecular evolution. Our analysis proves, therefore, that the role of opportunism in evolution has converted a problem of several possible chemical solutions into asingle-solution problem, with the actual Krebs cycle demonstrated to be the best possible chemical design. Our results also allow us to derive the rules under which metabolic pathways emerged during the origin of life.     

The Rooting of the Universal Tree of Life Is Not Reliable

Several composite universal trees connected by an ancestral gene duplication have been used to root the universal tree of life. In all cases, this root turned out to be in the eubacterial branch. However, the validity of results obtained from comparative sequence analysis has recently been questioned, in particular, in the case of ancient phylogenies. For example, it has been shown that several eukaryotic groups are misplaced in ribosomal RNA or elongation factor trees because of unequal rates of evolution and mutational saturation. Furthermore, the addition of new sequences to data sets has often turned apparently reasonable phylogenies into confused ones. We have thus revisited all composite protein trees that have been used to root the universal tree of life up to now (elongation factors, ATPases, tRNA synthetases, carbamoyl phosphate synthetases, signal recognition particle proteins) with updated data sets. In general, the two prokaryotic domains were not monophyletic with several aberrant groupings at different levels of the tree. Furthermore, the respective phylogenies contradicted each others, so that various ad hoc scenarios (paralogy or lateral gene transfer) must be proposed in order to obtain the traditional Archaebacteria–Eukaryota sisterhood. More importantly, all of the markers are heavily saturated with respect to amino acid substitutions. As phylogenies inferred from saturated data sets are extremely sensitive to differences in evolutionary rates, present phylogenies used to root the universal tree of life could be biased by the phenomenon of long branch attraction. Since the eubacterial branch was always the longest one, the eubacterial rooting could be explained by an attraction between this branch and the long branch of the outgroup. Finally, we suggested that an eukaryotic rooting could be a more fruitful working hypothesis, as it provides, for example, a simple explanation to the high genetic similarity of Archaebacteria and Eubacteria inferred from complete genome analysis.     

Taxonomy of 5 S Ribosomal RNA by the Linguistic Technique: Probing with Mitochondrial and Mammalian Sequences

Linguistic similarities and dissimilarities between 5 S rRNA sequences allowed taxonomical separation of species and classes. Comparisons with the molecule from mammals distinguished fungi and plants from protists and animals. Similarities to mammalians progressively increased from protists to invertebrates and to somatic-type molecules of the vertebrates lineage. In this, deviations were detected in avian, oocyte type, and pseudogene sequences. Among bacteria, actinobacteria were most similar to the mammalians, which could be related to the high frequency of associations among members of these groups. Some archaebacterial species most similar to the mammalians belonged to the Thermoproteales and Halobacteria groups. Comparisons with the soybean mitochondrial molecule revealed high internal homogeneity among plant mitochondria. The eubacterial groups most similar to it were Thermus and Rhodobacteria γ-1 and α-2. Other procedures have already indicated similarities of Rhodobacteria α to mitochondria but the linguistic similarities were on the average higher with the first two groups. Received: 5 August 1996 / Accepted: 9 April 1997     

RNA Folding Argues Against a Hot-Start Origin of Life

Opinion is strongly divided on whether life arose on earth under hot or cold conditions, the hot-start and cold-start scenarios, respectively. The origin of life close to deep thermal vents appears as the majority opinion among biologists, but there is considerable biochemical evidence that high temperatures are incompatible with an RNA world. To be functional, RNA has to fold into a three-dimensional structure. We report both theoretical and experimental results on RNA folding and show that (as expected) hot conditions strongly reduce RNA folding. The theoretical results come from energy-minimization calculations of the average extent of folding of RNA, mainly from 0–90°C, for both random sequences and tRNA sequences. The experimental results are from circular-dichroism measurements of tRNA over a similar range of temperatures. The quantitative agreement between calculations and experiment is remarkable, even to the shape of the curves indicating the cooperative nature of RNA folding and unfolding. These results provide additional evidence for a lower temperature stage being necessary in the origin of life. Received: 1 March 2000 / Accepted: 14 June 2000     

An Alphoid-Like Satellite DNA Sequence Is Present in the Genome of a Lacertid Lizard

A PstI DNA family was isolated from the genome of a lacertid, Lacerta graeca. The 185-bp monomeric unit (pGPS) was cloned and hybridized to DNAs and chromosomes of several lacertid species. The data showed that pGPS hybridizes to the (1) centromeric or pericentromeric heterochromatin of almost all the chromosomes of L. graeca and (2) genomic DNA of species phylogenetically related and unrelated to L. graeca. The presence of pGPS even in species immunologically apart more than 30 million years suggests that this repeated family might be either very ancient or have been conserved during evolution due to its functional role. The latter hypothesis might be supported by the results of sequence analysis which showed some homology with both several alphoid sequences of primates and the CDEIII centromeric sequence of yeast. Segments of the satellite sequence are similar to the mammalian CENP-B box. These observations suggest that pGPS might have a role in determining the centromeric function in lacertid lizards. Received: 6 February 1997 / Accepted: 14 May 1997     

Contribution of Mutation and RNA Recombination to the Evolution of a Plant Pathogenic RNA

The nucleotide sequence of 17 variants of the satellite RNA of cucumber mosaic virus (CMV-satRNA) isolated from field-infected tomato plants in the springs of 1989, 1990, and 1991 was determined. The sequence of each of the 17 satRNAs was unique and was between 334 and 340 nucleotides in length; 57 positions were polymorphic. There was much genetic divergence, ranging from 0.006 to 0.141 nucleotide substitutions per site for pairwise comparisons, and averaging 0.074 for any pair. When the polymorphic positions were analyzed relative to a secondary structure model proposed for CMV-satRNAs, it was found that there were significantly different numbers of changes in base-paired and non–base-paired positions, and that mutations that did not disrupt base pairing were preferred at the putatively paired sites. This supports the concept that the need to maintain a functional structure may limit genetic divergence of CMV-satRNA. Phylogenetic analyses showed that the 17 CMV-satRNA variants clustered into two subgroups, I and II, and evolutionary lines proceeding by the sequential accumulation of mutations were apparent. Three satRNA variants were outliers for these two phylogenetic groups. They were shown to be recombinants of subgroup I and II satRNAs by calculating phylogenies for different molecular regions and by using Sawyer's test for gene conversion. At least two recombination events were required to produce these three recombinant satRNAs. Thus, recombinants were found to be frequent (∼17%) in natural populations of CMV-satRNA, and recombination may make an important contribution to the generation of new variants. To our knowledge this is the first report of data allowing the frequency of recombinant isolates in natural populations of an RNA replicon to be estimated. Received: 14 May 1996 / Accepted: 17 July 1996     

The Robust Statistical Bases of the Coevolution Theory of Genetic Code Origin

A paper (Amirnovin R, J Mol Evol 44:473–476, 1997) seems to undermine the validity of the coevolution theory of genetic code origin by shedding doubt on the connection between the biosynthetic relationships between amino acids and the organization of the genetic code, at a time when the literature on the topic takes this for granted. However, as a few papers cite this paper as evidence against the coevolution theory, and to cast aside all doubt on the subject, we have decided to reanalyze the statistical bases on which this theory is founded. We come to the following conclusions: (1) the methods used in the above referred paper contain certain mistakes, and (2) the statistical foundations on which the coevolution theory is based are extremely robust. We have done this by critically appraising Amirnovin's paper and suggesting an alternative method based on the generation of random codes which, along with the method reported in the literature, allows us to evaluate the significance, in the genetic code, of different sets of amino acid pairs in biosynthetic relationships. In particular, by using this method and after building up a certain set of amino acid pairs reflecting the expectations of the coevolution theory, we show that the presence of this set in the genetic code would be obtained, purely by chance, with a probability of 6 × 10⁻⁵. This observation seems to provide particularly strong support to the coevolution theory. Received: 28 June 1999 / Accepted: 23 October 1999     

The MADS-box floral homeotic gene lineages predate the origin of seed plants: Phylogenetic and molecular clock estimates

Flower development in angiosperms is controlled in part by floral homeotic genes, many of which are members of the plant MADS-box regulatory gene family. The evolutionary history of these developmental genes was reconstructed using 74 loci from 15 dicot, three monocot, and one conifer species. Molecular clock estimates suggest that the different floral homeotic gene lineages began to diverge from one another about 450–500 mya, around the time of the origin of land plants themselves. Received: 31 January 1997 / Accepted: 9 April 1997