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     

Genetic codes as evolutionary filters: subtle differences in the structure of genetic codes result in significant differences in patterns of nucleotide substitution

The codon-degeneracy model (CDM) predicts that patterns of nucleotide substitution in protein-coding genes are largely determined by the relative frequencies of four-fold (4f), two-fold, and non-degenerate sites, the attributes of which are determined by the structure of the governing genetic code. The CDM thus further predicts that genetic codes with alternative structures will "filter" molecular evolution differentially. A method, therefore, is presented by which the CDM may be applied to the unique structure of any genetic code. The mathematical relationship between the proportion of transitions at 4f degenerate nucleotide sites and the transition-to-transversion ratio is described. Predictions for five individual genetic codes, relative to the relationship between code structure and expected patterns of nucleotide substitution, are clearly defined. To test this "filter" hypothesis of genetic codes, simulated DNA sequence data sets were generated with a variety of input parameter values to estimate the relationship between patterns of nucleotide substitution and best-fit estimates of transition bias at 4f degenerate sites for both the universal genetic code and the vertebrate mitochondrial genetic code. These analyses confirm the prediction of the CDM that, all else being equal, even small differences in the structure of alternative genetic codes may result in significant shifts in the overall pattern of nucleotide substitution.     

Two computer programs for rapid entry of DNA sequence data

Two computer programs for the IBM personal computer are describedfor rapid and accurate entry of DNA sequence data. The DNA sequencefiles produced can be used directly by the DNA sequence manipulationprograms by R. Staden (the DataBase system), the Universityof Wisconsin Genetics Computer Group, DNASTAR, or D. Mount.The first program, DIGISEQ, utilizes a sonic digitizer for semi-automationof sequence entry. To enter the DNA sequence each band of agel reading is touched by the stylus of the sonic digitizer.DIGISEQ corrects for both changes in lane width and lane curvature.The algorithm is extremely efficient and rarely requires re-entenngthe centers of the lanes. The second program, TYPESEQ, usesonly the keyboard for input. The keyboard is reconfigured toplace nucleotides and ambiguity codes under the fingers of onehand, corresponding to the order of the nucleotides on the geldefined by the user Both programs produce individual tones foreach nucleotide, and certain ambiguity codes. This verifiesinput of the correct nucleotide or ambiguity code, and thuseliminates the need to visually check the screen display duringsequence entry. Received on November 16, 1986; accepted on June 16, 1987     

Characterization of the bovine endogenous retrovirus beta3 genome

We recently used degenerate PCR and locus-specific PCR methods to identify the endogenous retroviruses (ERV) in the bovine genome. Using the ovine ERV classification system, the bovine ERVs (BERVs) could be classified into four families. Here, we searched the most recently released bovine genome database with the partial nucleotide sequence of the pro/pol region of the BERV beta3 family. This allowed us to obtain and analyze the complete genome of BERV beta3. The BERV beta3 genome is 7666 nucleotides long and has the typical retroviral organization, namely, 5'-long terminal repeat (LTR)-gag-pro-pol-env-LTR-3'. The deduced open reading frames for gag, pro, pol and env of BERV Beta en- code 507, 271, 879 and 603 amino acids, respectively. BERV beta3 showed little amino acid similarity to other betaretroviruses. Phylogenetic analysis showed that it clusters with HERV-K. This is the first report describing the genetic structure and sequence of an entire BERV.     

The nucleotide sequence of the murine I-E beta b immune response gene: evidence for gene conversion events in class II genes of the major histocompatibility complex.

We have determined the DNA sequence of the murine I-E beta b immune response gene of the major histocompatibility complex (MHC) of the C57BL/10 mouse and compared it with the sequence of allelic I-E and non-allelic I-A genes from the d and k haplotypes. The polymorphic exon sequences which encode the first extracellular globular domain of the E beta domain show approximately 8% nucleotide substitutions between the E beta b and E beta d alleles compared with only approximately 2% substitutions for the intron sequences. This suggests that an active mechanism such as micro gene conversion events drive the accumulation of these mutations in the polymorphic exons. The fact that several of the nucleotide changes are clustered supports this hypothesis. The E beta b and E beta k genes show approximately 2-fold fewer nucleotide substitutions than the E beta d/E beta b pair. The A beta bm12, a mutant I-A beta b gene from the C57BL/6 mouse, has been shown to result from three nucleotide changes clustered in a short region of the beta 1 domain, which suggests that a micro gene conversion event caused this mutation. We show here that the E beta b gene is identical to the non-allelic A beta bm12 DNA sequence in the mutated region and suggest, therefore, that the E beta b gene was the donor sequence for this intergenic transfer of genetic information. Diversity in class II MHC genes appears therefore to be generated, at least in part, by the same mechanism proposed for class I genes: intergenic transfer of short DNA regions between non-allelic genes.     

Targeted isolation,sequence analysis,and physical mapping of nonTIR NBS-LRR genes in soybean

Most cloned plant disease resistance genes (R-genes) code for proteins belonging to the nucleotide binding site (NBS) leucine-rich repeat (LRR) superfamily. NBS-LRRs can be divided into two classes based on the presence of a TIR domain (Toll and interleukin receptor-like sequence) or a coiled coil motif (nonTIR) in their N-terminus. We used conserved motifs specific to nonTIR-NBS-LRR sequences in a targeted PCR approach to generate nearly 50 genomic soybean sequences with strong homology to known resistance gene analogs (RGAs) of the nonTIR class. Phylogenetic analysis classified these sequences into four main subclasses. A representative clone from each subclass was used for genetic mapping, bacterial artificial chromosome (BAC) library screening, and construction of RGA-containing BAC contigs. Of the 14 RGAs that could be mapped genetically, 12 localized to a 25-cM region of soybean linkage group F already known to contain several classical disease resistance loci. A majority of the genomic region encompassing the RGAs was physically isolated in eight BAC contigs, together spanning more than 1 Mb of genomic sequence with at least 12 RGA copies. Phylogenetic and sequence analysis, together with genetic and physical mapping, provided insights into the genome organization and evolution of this large cluster of soybean RGAs. Received: 8 May 2001 / Accepted: 30 June 2001     

DNA shape, genetic codes, and evolution

Although the three-letter genetic code that maps nucleotide sequence to protein sequence is well known, there must exist other codes that are embedded in the human genome. Recent work points to sequence-dependent variation in DNA shape as one mechanism by which regulatory and other information could be encoded in DNA. Recent advances include the discovery of shape-dependent recognition of DNA that depends on minor groove width and electrostatics, the existence of overlapping codes in protein-coding regions of the genome, and evolutionary selection for compensatory changes in nucleotide composition that facilitate nucleosome occupancy. It is becoming clear that DNA shape is important to biological function, and therefore will be subject to evolutionary constraint.     

Genetic code ambiguity confers a selective advantage on Acinetobacter baylyi

A primitive genetic code, composed of a smaller set of amino acids, may have expanded via recursive periods of genetic code ambiguity that were followed by specificity. Here we model a step in this process by showing how genetic code ambiguity could result in an enhanced growth rate in Acinetobacter baylyi.     

The genetic code in bovine mitochondria: sequence of genes for the cytochrome oxidase subunit II and two tRNAs

Bovine-heart mitochondrial DNA from a single animal was isolated and fragments representative of the entire genome cloned into multicopy plasmid vectors to facilitate determination of its complete nucleotide sequence. We present here the sequence of the region covering the gene for cytochrome oxidase subunit II. Comparison of this sequence with the amino acid sequence of the homologous beef-heart protein has enabled the determination of most of the bovine mitochondrial genetic code. The code differs from the "universal" genetic code in that UGA codes for tryptophan and not termination, and AUA codes for methionine and not isoleucine. The only codon family not represented is the AGA/AGG pair normally used for arginine; evidence from other genes suggests that these code for termination in bovine mitochondria. The sequence presented also includes the adjacent tRNAAsp and tRNALys genes. The tRNAAsp gene is separated by one nucleotide from the 5' end of the COII gene and only three bases separate the 3' end of this gene and the adjacent tRNALys gene. This highly compact gene organisation is very similar to that found in the corresponding region of the human mitochondrial genome and the gene arrangement is identical. The structure of the respective bovine and human tRNAs vary primarily the "D-" and "T psi C-loops".     

Total synthesis of a structural gene for the human peptide hormone angiotensin II.

Seven oligonucleotide chains containing between 6 and 11 nucleotide units were synthesized. The segments were phosphorylated by T4 polynucleotide 5'-hydroxyl-kinase and joined by T4 polynucleotide synthetase (ATP) to give the double-stranded DNA consisting of 33 base pairs. The DNA sequence was deduced from the known peptide sequence according to the genetic code.     

The non-standard genetic code of Candida spp.: an evolving genetic code or a novel mechanism for adaptation?

A number of yeasts of the genus Candida translate the standard leucine-CUG codon as serine. This unique genetic code change is the only known alteration to the universal genetic code in cytoplasmic mRNAs, of either eukaryotes or prokaryotes, which involves reassignment of a sense codon. Translation of CUG as serine in these species is mediated by a novel serine-tRNA (ser-tRNA_CAG), which uniquely has a guanosine at position 33, 5' to the anticodon, a position that is almost invariably occupied by a pyrimidine (uridine in general) in all other tRNAs. We propose that G-33 has two important functions: lowering the decoding efficiency of the ser-tRNA_CAG and preventing binding of the leucyl-tRNA synthetase. This implicates this nucleotide as a key player in the evolutionary reassignment of the CUG codon. In addition, the novel ser-tRNA_CAG has 1-methylguanosine (m¹G-37) at position 37, 3' to the anticodon, which is characteristic of leucine, but not serine tRNAs. Remarkably, m¹G-37 causes leucylation of the ser-tRNA_CAG both in vitro and in vivo , making the CUG codon an ambiguous codon: the polysemous codon. This indicates that some Candida species tolerate ambiguous decoding and suggests either that (i) the genetic code change has not yet been fully established and is evolving at different rates in different Candida species; or (ii) CUG ambiguity is advantageous and represents the final stage of the reassignment. We propose that such dual specificity indicates that reassignment of the CUG codon evolved through a mechanism that required codon ambiguity and that ambiguous decoding evolved to generate genetic diversity and allow for rapid adaptation to environmental challenges.     

Beta thalassemic mutations recognized by DNA mapping with Hph I and Rsa I in the Algerian population

By using Hph I and Rsa I restriction enzymes and beta globin large intervening sequence as a probe, we have investigated the DNA of 20 Algerian patients with beta(0) or beta(+) thalassemia. In any of them, we detected the nucleotide change which is known to generate an additional Hph I site at the 5' splice junction of the beta globin large intervening sequence and which yields a beta(0) phenotype. In one of them, we detected the nucleotide change which is known to generate an additional Rsa I site within the beta globin large intervening sequence and which is supposed to yield a beta(+) phenotype. These results indicate that these two types of mutation are relatively rare in the Algerian population.     

Selective advantages created by codon ambiguity allowed for the evolution of an alternative genetic code in Candida spp

Several species of the genus Candida decode the standard leucine CUG codon as serine. This and other deviations from the standard genetic code in both nuclear and mitochondrial genomes invalidate the notion that the genetic code is frozen and universal and prompt the questions ‘why alternative genetic codes evolved and, more importantly, how can an organism survive a genetic code change?’ To address these two questions, we have attempted to reconstruct the early stages of Candida albicans CUG reassignment in the closely related yeast Saccharomyces cerevisiae. These studies suggest that this genetic code change was driven by selection using a molecular mechanism that requires CUG ambiguity. Such codon ambiguity induced a significant decrease in fitness, indicating that CUG reassignment can only be selected if it introduces an evolutionary edge to counteract the negative impact of ambiguity. We have shown that CUG ambiguity induces the expression of a novel set of stress proteins and triggers the general stress response, which, in turn, creates a competitive edge under stress conditions. In addition, CUG ambiguity in S. cerevisiae induces the expression of a number of novel phenotypes that mimic the natural resistance to stress characteristic of C. albicans. The identification of an evolutionary advantage created by CUG ambiguity is the first experimental evidence for a genetic code change driven by selection and suggests a novel role for codon reassignment in the adaptation to new ecological niches.     

Pseudogene rescue: an adaptive mechanism of codon reassignment

Alterations to the genetic code – codon reassignments – have occurred many times in life’s history, despite the fact that genomes are coadapted to their genetic codes and therefore alterations are likely to be maladaptive. A potential mechanism for adaptive codon reassignment, which could trigger either a temporary period of codon ambiguity or a permanent genetic code change, is the reactivation of a pseudogene by a nonsense suppressor mutant transfer RNA. I examine the population genetics of each stage of this process and find that pseudogene rescue is plausible and also readily explains some features of extant variability in genetic codes.     

The complete nucleotide sequence of RNA beta from the type strain of barley stripe mosaic virus.

The complete nucleotide sequence of RNA beta from the type strain of barley stripe mosaic virus (BSMV) has been determined. The sequence is 3289 nucleotides in length and contains four open reading frames (ORFs) which code for proteins of Mr 22,147 (ORF1), Mr 58,098 (ORF2), Mr 17,378 (ORF3), and Mr 14,119 (ORF4). The predicted N-terminal amino acid sequence of the polypeptide encoded by the ORF nearest the 5'-end of the RNA (ORF1) is identical (after the initiator methionine) to the published N-terminal amino acid sequence of BSMV coat protein for 29 of the first 30 amino acids. ORF2 occupies the central portion of the coding region of RNA beta and ORF3 is located at the 3'-end. The ORF4 sequence overlaps the 3'-region of ORF2 and the 5'-region of ORF3 and differs in codon usage from the other three RNA beta ORFs. The coding region of RNA beta is followed by a poly(A) tract and a 238 nucleotide tRNA-like structure which are common to all three BSMV genomic RNAs.     

The distribution of the frequency of occurrence of nucleotide subsequences, based on their overlap capability

DNA's genetic code can be represented as an alphabetic sequence composed of the four letters A, C, G, and T, which represent the four types of nucleotides--adenylic, cytidylic, guanylic, and thymidylic acid--of which DNA is composed. Now that these sequences have been identified for many genes and are available in computer-readable form, scientists can analyze these data and search for patterns in an attempt to learn more about the regulatory functions of the gene. One area of study is that of the frequency of occurrence of specific nucleotide subsequences (e.g., ACAC) within part or all of a nucleotide sequence. This paper derives the probability distribution of the frequency of occurrence of a subsequence within a nucleotide sequence, under the hypothesis that the four nucleotides occur at random and with equal probability. This distribution is nontrivial because different subsequences have different "overlap capability." For example, the subsequence AAAA can occur up to 17 times in a sequence of length 20 (which would happen if the sequence were composed solely of A's), but the subsequence ACGT cannot occur more than 5 times in a sequence of length 20. Thus, the frequency distributions are different for each type of overlap capability. It is of interest to assess and compare the degree of nonrandomness for different subsequences or among different portions of a sequence; the existence and degree of nonrandomness may be related to the type and degree of functionality of a nucleotide (sub)sequence. The frequency distributions provided here can be used to perform exact significance tests of the hypothesis of randomness. An approximate test is also described for use with long sequences; this can be used to test a more general null hypothesis of nucleotides occurring with unequal probabilities.     

Complexity and sequence identification of 24 rat V beta genes

Twenty-four TCR V beta genes were cloned by anchored PCR from the Lewis rat strain and identified by nucleotide and amino acid sequence comparisons to known mouse V beta genes. Rat V beta genes exist in 17 single-member and 3 multimember subfamilies and exhibit 86 to 94 and 72 to 92% nucleotide and amino acid sequence similarities, respectively, to their mouse counterparts. A single rat gene, designated V beta 20, having no previously known mouse counterpart was identified; a closely related gene was subsequently isolated from the C57BL/6 mouse strain. Characterization of the rV beta genes provides the basis for addressing rat T cell tolerance mechanisms and V beta gene usage in immunologically mediated rat disease models.     

The sequence of a cloned cDNA for the beta subunit of bovine thyrotropin predicts a protein containing both NH2-and COOH-terminal extensions

Cloned cDNAs containing sequences coding for the beta subunit of bovine thyrotropin have been identified. The complete nucleotide sequence of the largest of the beta subunit cDNA inserts has been determined. This cDNA contains 35 nucleotides from the 5' untranslated region of thyrotropin beta subunit mRNA and 60 nucleotides coding for an NH2-terminal precursor segment. This is followed by 339 nucleotides which code for the published amino acid sequence of the thyrotropin beta subunit. Following the 339 nucleotide beta subunit coding sequence, no termination codon is encountered for another 15 nucleotides. Thus, the cDNA codes for a thyrotropin beta subunit containing an additional 5 amino acids at the COOH terminus. The cDNA also contains 82 nucleotides of 3' untranslated sequence followed by a short poly(A) segment. Comparison of the bovine cDNA sequence to the recently described mouse thyrotropin beta subunit cDNA sequence reveals considerable homology throughout the coding sequence, including the COOH-terminal extension. These findings suggest the possibility that a thyrotropin beta subunit precursor is processed at both the NH2 and COOH termini.