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 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     

Physicochemical Optimization in the Genetic Code Origin as the Number of Codified Amino Acids Increases

We have assumed that the coevolution theory of genetic code origin (Wong JT, Proc Natl Acad Sci USA 72:1909–1912, 1975) is essentially correct. This theory makes it possible to identify at least 10 evolutionary stages through which genetic code organization might have passed prior to reaching its current form. The calculation of the minimization level of all these evolutionary stages leads to the following conclusions. (1) The minimization percentages increased linearly with the number of amino acids codified in the codes of the various evolutionary stages when only the sense changes are considered in the analysis. This seems to favor the physicochemical theory of genetic code origin even if, as discussed in the paper, this observation is also compatible with the coevolution theory. (2) For the first seven evolutionary stages of the genetic code, this trend is less clear and indeed is inverted when we consider the global optimisation of the codes due to both sense changes and synonymous changes. This inverse correlation between minimization percentages and the number of amino acids codified in the codes of the intermediate stages seems to favor neither the physicochemical nor the stereochemical theories of genetic code origin, as it is in the early and intermediate stages of code development that these theories would expect minimization to have played a crucial role, and this does not seem to be the case. However, these results are in agreement with the coevolution theory, which attributes a role to the physicochemical properties of amino acids that, while important, is nevertheless subordinate to the mechanism which concedes codons from the precursor amino acids to the product amino acids as the primary factor determining the evolutionary structuring of the genetic code. The results are therefore discussed in the context of the various theories proposed to explain genetic code origin. Received: 25 October 1998 / Accepted: 19 February 1999     

An Analysis of the Metabolic Theory of the Origin of the Genetic Code

A computer program was used to test Wong's coevolution theory of the genetic code. The codon correlations between the codons of biosynthetically related amino acids in the universal genetic code and in randomly generated genetic codes were compared. It was determined that many codon correlations are also present within random genetic codes and that among the random codes there are always several which have many more correlations than that found in the universal code. Although the number of correlations depends on the choice of biosynthetically related amino acids, the probability of choosing a random genetic code with the same or greater number of codon correlations as the universal genetic code was found to vary from 0.1% to 34% (with respect to a fairly complete listing of related amino acids). Thus, Wong's theory that the genetic code arose by coevolution with the biosynthetic pathways of amino acids, based on codon correlations between biosynthetically related amino acids, is statistical in nature. Received: 8 August 1996 / Accepted: 26 December 1996     

A Novel Theory on the Origin of the Genetic Code: A GNC-SNS Hypothesis

We have previously proposed an SNS hypothesis on the origin of the genetic code (Ikehara and Yoshida 1998). The hypothesis predicts that the universal genetic code originated from the SNS code composed of 16 codons and 10 amino acids (S and N mean G or C and either of four bases, respectively). But, it must have been very difficult to create the SNS code at one stroke in the beginning. Therefore, we searched for a simpler code than the SNS code, which could still encode water-soluble globular proteins with appropriate three-dimensional structures at a high probability using four conditions for globular protein formation (hydropathy, α-helix, β-sheet, and β-turn formations). Four amino acids (Gly [G], Ala [A], Asp [D], and Val [V]) encoded by the GNC code satisfied the four structural conditions well, but other codes in rows and columns in the universal genetic code table do not, except for the GNG code, a slightly modified form of the GNC code. Three three-amino acid systems ([D], Leu and Tyr; [D], Tyr and Met; Glu, Pro and Ile) also satisfied the above four conditions. But, some amino acids in the three systems are far more complex than those encoded by the GNC code. In addition, the amino acids in the three-amino acid systems are scattered in the universal genetic code table. Thus, we concluded that the universal genetic code originated not from a three-amino acid system but from a four-amino acid system, the GNC code encoding [GADV]-proteins, as the most primitive genetic code. Received: 11 June 2001 / Accepted: 11 October 2001     

On Error Minimization in a Sequential Origin of the Standard Genetic Code

Distances between amino acids were derived from the polar requirement measure of amino acid polarity and Benner and co-workers' (1994) 74-100 PAM matrix. These distances were used to examine the average effects of amino acid substitutions due to single-base errors in the standard genetic code and equally degenerate randomized variants of the standard code. Second-position transitions conserved all distances on average, an order of magnitude more than did second-position transversions. In contrast, first-position transitions and transversions were about equally conservative. In comparison with randomized codes, second-position transitions in the standard code significantly conserved mean square differences in polar requirement and mean Benner matrix-based distances, but mean absolute value differences in polar requirement were not significantly conserved. The discrepancy suggests that these commonly used distance measures may be insufficient for strict hypothesis testing without more information. The translational consequences of single-base errors were then examined in different codon contexts, and similarities between these contexts explored with a hierarchical cluster analysis. In one cluster of codon contexts corresponding to the RNY and GNR codons, second-position transversions between C and G and transitions between C and U were most conservative of both polar requirement and the matrix-based distance. In another cluster of codon contexts, second-position transitions between A and G were most conservative. Despite the claims of previous authors to the contrary, it is shown theoretically that the standard code may have been shaped by position-invariant forces such as mutation and base content. These forces may have left heterogeneous signatures in the code because of differences in translational fidelity by codon position. A scenario for the origin of the code is presented wherein selection for error minimization could have occurred multiple times in disjoint parts of the code through a phyletic process of competition between lineages. This process permits error minimization without the disruption of previously useful messages, and does not predict that the code is optimally error-minimizing with respect to modern error. Instead, the code may be a record of genetic process and patterns of mutation before the radiation of modern organisms and organelles. Received: 28 July 1997 / Accepted: 23 January 1998     

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     

Codon Usage Decreases the Error Minimization Within the Genetic Code

The genetic code is not random but instead is organized in such a way that single nucleotide substitutions are more likely to result in changes between similar amino acids. This fidelity, or error minimization, has been proposed to be an adaptation within the genetic code. Many models have been proposed to measure this adaptation within the genetic code. However, we find that none of these consider codon usage differences between species. Furthermore, use of different indices of amino acid physicochemical characteristics leads to different estimations of this adaptation within the code. In this study, we try to establish a more accurate model to address this problem. In our model, a weighting scheme is established for mistranslation biases of the three different codon positions, transition/transversion biases, and codon usage. Different indices of amino acids physicochemical characteristics are also considered. In contrast to pervious work, our results show that the natural genetic code is not fully optimized for error minimization. The genetic code, therefore, is not the most optimized one for error minimization, but one that balances between flexibility and fidelity for different species.     

The Impact of Message Mutation on the Fitness of a Genetic Code

The Standard Genetic Code is organized such that similar codons encode similar amino acids. One explanation suggested that the Standard Code is the result of natural selection to reduce the fitness ``load' that derives from the mutation and mistranslation of protein-coding genes. We review the arguments against the mutational load-minimizing hypothesis and argue that they need to be reassessed. We review recent analyses of the organization of the Standard Code and conclude that under cautious interpretation they support the mutational load-minimizing hypothesis. We then present a deterministic asexual model with which we study the mode of selection for load minimization. In this model, individual fitness is determined by a protein phenotype resulting from the translation of a mutable set of protein-coding genes. We show that an equilibrium fitness may be associated with a population with the same genetic code and that genetic codes that assign similar codons to similar amino acids have a higher fitness. We also show that the number of mutant codons in each individual at equilibrium, which determines the strength of selection for load minimization, reflects a long-term evolutionary balance between mutations in messages and selection on proteins, rather than the number of mutations that occur in a single generation, as has been assumed by previous authors. We thereby establish that selection for mutational load minimization acts at the level of an individual in a single generation. We conclude with comments on the shortcomings and advantages of load minimization over other hypotheses for the origin of the Standard Code. Received: 4 April 2001 / Accepted: 22 October 2001     

Relationships Among Isoacceptor tRNAs Seems to Support the Coevolution Theory of the Origin of the Genetic Code

A new method for looking at relationships between nucleotide sequences has been used to analyze divergence both within and between the families of isoaccepting tRNA sets. A dendrogram of the relationships between 21 tRNA sets with different amino acid specificities is presented as the result of the analysis. Methionine initiator tRNAs are included as a separate set. The dendrogram has been interpreted with respect to the final stage of the evolutionary pathway with the development of highly specific tRNAs from ambiguous molecular adaptors. The location of the sets on the dendrogram was therefore analyzed in relation to hypotheses on the origin of the genetic code: the coevolution theory, the physicochemical hypothesis, and the hypothesis of ambiguity reduction of the genetic code. Pairs of 16 sets of isoacceptor tRNAs, whose amino acids are in biosynthetic relationships, occupied contiguous positions on the dendrogram, thus supporting the coevolution theory of the genetic code. Received: 4 May 1998 / Accepted: 11 July 1998     

How Mitochondria Redefine the Code

Annotated, complete DNA sequences are available for 213 mitochondrial genomes from 132 species. These provide an extensive sample of evolutionary adjustment of codon usage and meaning spanning the history of this organelle. Because most known coding changes are mitochondrial, such data bear on the general mechanism of codon reassignment. Coding changes have been attributed variously to loss of codons due to changes in directional mutation affecting the genome GC content (Osawa and Jukes 1988), to pressure to reduce the number of mitochondrial tRNAs to minimize the genome size (Anderson and Kurland 1991), and to the existence of transitional coding mechanisms in which translation is ambiguous (Schultz and Yarus 1994a). We find that a succession of such steps explains existing reassignments well. In particular, (1) Genomic variation in the prevalence of a codon's third-position nucleotide predicts relative mitochondrial codon usage well, though GC content does not. This is because A and T, and G and C, are uncorrelated in mitochondrial genomes. (2) Codons predicted to reach zero usage (disappear) do so more often than expected by chance, and codons that do disappear are disproportionately likely to be reassigned. However, codons predicted to disappear are not significantly more likely to be reassigned. Therefore, low codon frequencies can be related to codon reassignment, but appear to be neither necessary nor sufficient for reassignment. (3) Changes in the genetic code are not more likely to accompany smaller numbers of tRNA genes and are not more frequent in smaller genomes. Thus, mitochondrial codons are not reassigned during demonstrable selection for decreased genome size. Instead, the data suggest that both codon disappearance and codon reassignment depend on at least one other event. This mitochondrial event (leading to reassignment) occurs more frequently when a codon has disappeared, and produces only a small subset of possible reassignments. We suggest that coding ambiguity, the extension of a tRNA's decoding capacity beyond its original set of codons, is the second event. Ambiguity can act alone but often acts in concert with codon disappearance, which promotes codon reassignment. Received: 26 October 2000 / Accepted: 19 January 2001     

Obcells as Proto-Organisms: Membrane Heredity, Lithophosphorylation, and the Origins of the Genetic Code, the First Cells, and Photosynthesis

I attempt to sketch a unified picture of the origin of living organisms in their genetic, bioenergetic, and structural aspects. Only selection at a higher level than for individual selfish genes could power the cooperative macromolecular coevolution required for evolving the genetic code. The protein synthesis machinery is too complex to have evolved before membranes. Therefore a symbiosis of membranes, replicators, and catalysts probably mediated the origin of the code and the transition from a nucleic acid world of independent molecular replicators to a nucleic acid/protein/lipid world of reproducing organisms. Membranes initially functioned as supramolecular structures to which different replicators attached and were selected as a higher-level reproductive unit: the proto-organism. I discuss the roles of stereochemistry, gene divergence, codon capture, and selection in the code's origin. I argue that proteins were primarily structural not enzymatic and that the first biological membranes consisted of amphipathic peptidyl-tRNAs and prebiotic mixed lipids. The peptidyl-tRNAs functioned as genetically-specified lipid analogues with hydrophobic tails (ancestral signal peptides) and hydrophilic polynucleotide heads. Protoribosomes arose from two cooperating RNAs: peptidyl transferase (large subunit) and mRNA-binder (small subunit). Early proteins had a second key role: coupling energy flow to the phosphorylation of gene and peptide precursors, probably by lithophosphorylation by membrane-anchored kinases scavenging geothermal polyphosphate stocks. These key evolutionary steps probably occurred on the outer surface of an `inside out-cell' or obcell, which evolved an unambiguous hydrophobic code with four prebiotic amino acids and proline, and initiation by isoleucine anticodon CAU; early proteins and nucleozymes were all membrane-attached. To improve replication, translation, and lithophosphorylation, hydrophilic substrate-binding and catalytic domains were later added to signal peptides, yielding a ten-acid doublet code. A primitive proto-ecology of molecular scavenging, parasitism, and predation evolved among obcells. I propose a new theory for the origin of the first cell: fusion of two cup-shaped obcells, or hemicells, to make a protocell with double envelope, internal genome and ribosomes, protocytosol, and periplasm. Only then did water-soluble enzymes, amino acid biosynthesis, and intermediary metabolism evolve in a concentrated autocatalytic internal cytosolic soup, causing 12 new amino acid assignments, termination, and rapid freezing of the 22-acid code. Anticodons were recruited sequentially: GNN, CNN, INN, and *UNN. CO₂ fixation, photoreduction, and lipid synthesis probably evolved in the protocell before photophosphorylation. Signal recognition particles, chaperones, compartmented proteases, and peptidoglycan arose prior to the last common ancestor of life, a complex autotrophic, anaerobic green bacterium. Received: 19 February 2001 / Accepted: 9 April 2001     

Genetic Code Variations in Mitochondria: tRNA as a Major Determinant of Genetic Code Plasticity

Characteristic features of tRNA such as the anticodon sequence and modified nucleotides in the anticodon loop are thought to be crucial effectors for promoting or restricting codon reassignment. Our recent findings on basepairing rules between anticodon and codon in various metazoan mitochondria suggest that the complete loss of a codon is not necessarily essential for codon reassignment to take place. We postulate that a possible competition between two tRNAs with cognate anticodon sequences towards the relevant codon to be varied has a potential role in codon reassignment. Our proposition can be viewed as an expanded version of the codon capture theory proposed by Osawa and Jukes (J Mol Evol 28: 271–278, 1989). Received: 28 December 2000 / Accepted: 12 March 2001     

Directionally Evolving Genetic Code: The UGA Codon from Stop to Tryptophan in Mitochondria

For the comprehensive analyses of deviant codes in protistan mitochondria (mt), we sequenced about a 1.1-kb region of a mitochondrial (mt) gene, the cytochrome c oxidase subunit I (coxI) in two chlorarachniophytes, the filose amoeba Euglypha rotunda, the cryptomonad Cryptomonas ovata, the prymnesiophyte (haptophyte) Diacronema vlkianum (Pavlovales), and the diatom Melosira ambigua. As a result of this analysis, we noticed that the UGA codon is assigned to tryptophan (Trp) instead of being a signal for translational termination in two chlorarachniophytes and in E. rotunda. The same type of deviant code was reported previously in animals, fungi, ciliates, kinetoplastids, Chondrus crispus (a red alga), Acanthamoeba castellanii (an amoeboid protozoon), and three of the four prymnesiophyte orders with the exception of the Pavlovales. A phylogenetic analysis based on the COXI sequences of 56 eukaryotes indicated that the organisms bearing the modified code, UGA for Trp, are not monophyletic. Based on these studies, we propose that the ancestral mitochondrion was bearing the universal genetic code and subsequently reassigned the codon to Trp independently, at least in the lineage of ciliates, kinetoplastids, rhodophytes, prymnesiophytes, and fungi. We also discuss how this codon was directionally captured by Trp tRNA. Received: 26 January 1998 / Accepted: 24 April 1998     

Synonymous and nonsynonymous rate variation in nuclear genes of mammals

A maximum likelihood approach was used to estimate the synonymous and nonsynonymous substitution rates in 48 nuclear genes from primates, artiodactyls, and rodents. A codon-substitution model was assumed, which accounts for the genetic code structure, transition/transversion bias, and base frequency biases at codon positions. Likelihood ratio tests were applied to test the constancy of nonsynonymous to synonymous rate ratios among branches (evolutionary lineages). It is found that at 22 of the 48 nuclear loci examined, the nonsynonymous/synonymous rate ratio varies significantly across branches of the tree. The result provides strong evidence against a strictly neutral model of molecular evolution. Our likelihood estimates of synonymous and nonsynonymous rates differ considerably from previous results obtained from approximate pairwise sequence comparisons. The differences between the methods are explored by detailed analyses of data from several genes. Transition/transversion rate bias and codon frequency biases are found to have significant effects on the estimation of synonymous and nonsynonymous rates, and approximate methods do not adequately account for those factors. The likelihood approach is preferable, even for pairwise sequence comparison, because more-realistic models about the mutation and substitution processes can be incorporated in the analysis. Received: 17 May 1997 / Accepted: 28 September 1997     

On the Crucial Stages in the Origin of Animate Matter

Theories of the origin of life have proposed hypotheses to link inanimate to animate matter. The theory proposed here derived the crucial stages in the origin of animate matter directly from the basic properties of inanimate matter. It asked what were the general characteristics of the link, rather than what might have been its chemical details. Life and its origin are shown to be one continuous physicochemical process of replication, random variation, and natural selection. Since life exists here and now, animate properties must have been initiated in the past somewhere. According to the theory, life originated from an as yet unknown elementary autocatalyst which occurred spontaneously, then replicated autocatalytically. As it multiplied to macroscopic abundance, its replicas gradually exhausted their reactants. Random chemical drift initiated diversity among autocatalysts. Diversity led to competition. Competition and depletion of reactants slowed down the rates of net replication of the autocatalysts. Some reached negative rates and became extinct, while those which stayed positive ``survived.' Thus chemical natural selection appeared, the first step in the transition from inanimate to animate matter. It initiated the first animate property, fitness, i.e., the capacity to adapt to the environment and to survive. As the environment was depleted of reactants, it was enriched with sequels—namely, with decomposition products and all other products which accompany autocatalysis. The changing environment exerted a selective pressure on autocatalysts to replace dwindling reactants by accumulating sequels. Sequels that were incorporated into the autocatalytic process became internal components of complex autocatalytic systems. Primitive forms of metabolism and organization were thus initiated. They evolved further by the same mechanism to ever higher levels of complexity, such as homochirality (handedness) and membranal enclosure. Subsequent evolution by the same mechanism generated cellular metabolism, cell division, information carriers, and a genetic code. Theories of self-organization without natural selection are refuted. Received: 29 March 1996 / Accepted: 30 May 1996     

A Sequential Scenario for the Origin of Biological Chirality

A sequential model is proposed regarding the origin of biological chirality. Three major stages are presumed: a symmetry breaking (prebiotic chiral disruption in enantiomeric mixtures of monomers), a chiral amplification (prebiotic increase of the chiral character of the monomers affected first by the symmetry breaking), and a chiral expansion (proto biological increase of the chiral character and spread of the chirality to molecules which were less affected by prebiotic chiralizations). As a symmetry-breaking mechanism, the model proposed by Deutsch (1991) is used, which involves a dissymmetric exposure of amino acids (AA) to ultraviolet circularly polarized light (UV-CPL) on evaporative seashores. It is presumed that the chiral amplification, up to a protobiologic significance, was influenced by a periodic overlapping of two abiotic events, a synchronization between tidal-based hydrous–anhydrous cycles, and littoral asymmetric photolysis cycles. This long-term astronomic asymmetry acted around 3.8–4.2 billion years ago and was unique to the Earth in our solar system. It is also presumed that the abiotic symmetry breaking is heterogenous, that only a few l-AAs were used in the beginning, and that the chirality expanded later to all 20 AAs based on a coevolutionary strategy of the genetic code and on a physiological relationship between AAs. In this scenario the d-chirality of pentoses in polynucleotides was attributed to both d-pentose/l-AA relationships and to a structural evolution. Received: 10 May 1996 / Accepted: 13 August 1996     

Amino Acids as RNA Ligands: A Direct-RNA-Template Theory for the Code's Origin

Numerous RNA binding sites for specific amino acids are now known, coming predominantly from selection-amplification experiments. These sites are chemically discriminating despite being predominantly small, simple RNA structures: internal and bulge loops. Recent studies of sites for hydrophobic side chains suggest that there are other generalizable structural features which recur in hydrophobic RNA sites. Further, sites for hydrophobic side chains can contain codons for the bound amino acid, as has also long been known for the polar amino acid arginine. Such findings are comprehensively reviewed, and the implications for the origin of coded peptide synthesis are considered. An origins hypothesis which accommodates all the data, DRT (direct RNA templating), is formulated. Received: 22 December 1997 / Accepted: 13 February 1998     

Estimation of the Transition/Transversion Rate Bias and Species Sampling

The transition/transversion (ti/tv) rate ratios are estimated by pairwise sequence comparison and joint likelihood analysis using mitochondrial cytochrome b genes of 28 primate species, representing both the Strepsirrhini (lemurs and lories) and the Anthropoidea (monkeys, apes, and humans). Pairwise comparison reveals a strong negative correlation between estimates of the ti/tv ratio and the sequence distance, even when both are corrected for multiple substitutions. The maximum-likelihood estimate of the ti/tv ratio changes with the species included in the analysis. The ti/tv bias within the lemuriform taxa is found to be as strong as in the anthropoids, in contradiction to an earlier study which sampled only one lemuriform. Simulations show the surprising result that both the pairwise correction method and the joint likelihood analysis tend to overcorrect for multiple substitutions and overestimate the ti/tv ratio, especially at low sequence divergence. The bias, however, is not large enough to account for the observed patterns. Nucleotide frequency biases, variation of substitution rates among sites, and different evolutionary dynamics at the three codon positions can be ruled out as possible causes. The likelihood-ratio test suggests that the ti/tv rate ratios may be variable among evolutionary lineages. Without any biological evidence for such a variation, however, we are left with no plausible explanations for the observed patterns other than a possible saturation effect due to the unrealistic nature of the model assumed. Received: 1 October 1997 / Accepted: 29 September 1998