On the evolution of redundancy in genetic codes.

We simulate a deterministic population genetic model for the coevolution of genetic codes and protein-coding genes. We use very simple assumptions about translation, mutation, and protein fitness to calculate mutation-selection equilibria of codon frequencies and fitness in a large asexual population with a given genetic code. We then compute the fitnesses of altered genetic codes that compete to invade the population by translating its genes with higher fitness. Codes and genes coevolve in a succession of stages, alternating between genetic equilibration and code invasion, from an initial wholly ambiguous coding state to a diversified frozen coding state. Our simulations almost always resulted in partially redundant frozen genetic codes. Also, the range of simulated physicochemical properties among encoded amino acids in frozen codes was always less than maximal. These results did not require the assumption of historical constraints on the number and type of amino acids available to codes nor on the complexity of proteins, stereochemical constraints on the translational apparatus, nor mechanistic constraints on genetic code change. Both the extent and timing of amino-acid diversification in genetic codes were strongly affected by the message mutation rate and strength of missense selection. Our results suggest that various omnipresent phenomena that distribute codons over sites with different selective requirements--such as the persistence of nonsynonymous mutations at equilibrium, the positive selection of the same codon in different types of sites, and translational ambiguity--predispose the evolution of redundancy and of reduced amino acid diversity in genetic codes.

Tentative identification of a resilin gene in Drosophila melanogaster.

A search of the Drosophila genome for gene products with similarities to the amino acid sequences of three tryptic peptides from locust (Schistocerca gregaria) resilin gave two positive results: gene products CG15920 and CG9036. In both conceptual translation products a 62-residue region is present, which is identical to the resilin peptides in 29 positions. Gene product CG15920 has an amino acid composition closely resembling that of resilins from various insect species, and it has an N-terminal signal peptide sequence indicating that it is an extracellular protein. The 62-residue region shows similarity to the RR-2 sequence, which is common for a number of matrix proteins from insect solid cuticle. The N- and C-terminal regions flanking the 62-residue in CG15920 are dominated by 18 repeats of a 15-residue sequence and 11 repeats of a 13-residue sequence, respectively. The structures of the repeats predict that the peptide chain will fold in an irregular, extended beta-spiral, resembling the structures suggested for mammalian elastin and spider flagelliform silk, two materials which, like resilin, possess long-range elasticity. Accordingly, we suggest that gene product CG15920 is a Drosophila resilin precursor.

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.