The organization and evolution of transfer RNA genes in Mycoplasma capricolum

The proteomic constraint and its role in molecular evolution

Recently, the concept of a "Proteomic Constraint" was introduced to explain the frequency of genetic code deviations in mitochondrial genomes. The Proteomic Constraint was proposed to be proportional to the size of the mitochondrially encoded proteome, hence small proteomes are expected to experience smaller total numbers of errors resulting from genetic code deviations, leading to less likelihood of causing lethality. The concept is now extended to encompass several other aspects of the genetic information system. When the Proteomic Constraint is small, it is proposed that there is little selective pressure to evolve or maintain error correction mechanisms, as a result of the smaller total number of errors that accumulate. Conversely, a large Proteomic Constraint is proposed to result in a correspondingly large selective pressure to evolve or maintain error correction mechanisms. Differences in the size of the Proteomic Constraint can help to explain differences in replicational, transcriptional, and translational fidelities between genomes. A key piece of evidence is the existence of negative power law relationships between proteome size and error rates; these are demonstrated to be diagnostic of the action of the Proteomic Constraint. The Proteomic Constraint is argued to be a major factor determining mutation rates in a diverse range of DNA genomes, implying that mutation rates are clock like. A small Proteomic Constraint partly explains why RNA viruses possess high mutation rates. A reduced Proteomic Constraint in intracellular pathogenic bacteria predicts a drift upwards in mutation rates. Differences in the Proteomic Constraint also appear to be linked to differences in recombination rates between eukaryotes. In addition, a reduced Proteomic Constraint may explain features of resident genomes, such as loss of DNA repair pathways, increased substitution rates, and AT biases, in addition to the occurrence of genetic code deviations. Thus, it is argued that the Proteomic Constraint is a universal factor that influences a wide range of properties of the genetic information system.

Archaeology and evolution of transfer RNA genes in the Escherichia coli genome

Transfer RNA genes tend to be presented in multiple copies in the genomes of most organisms, from bacteria to eukaryotes. The evolution and genomic structure of tRNA genes has been a somewhat neglected area of molecular evolution. Escherichia coli, the first phylogenetic species for which more than two different strains have been sequenced, provides an invaluable framework to study the evolution of tRNA genes. In this work, a detailed analysis of the tRNA structure of the genomes of Escherichia coli strains K12, CFT073, and O157:H7, Shigella flexneri 2a 301, and Salmonella typhimurium LT2 was carried out. A phylogenetic analysis of these organisms was completed, and an archaeological map depicting the main events in the evolution of tRNA genes was drawn. It is shown that duplications, deletions, and horizontal gene transfers are the main factors driving tRNA evolution in these genomes. On average, 0.64 tRNA insertions/duplications occur every million years (Myr) per genome per lineage, while deletions occur at the slower rate of 0.30 per million years per genome per lineage. This work provides a first genomic glance at the problem of tRNA evolution as a repetitive process, and the relationship of this mechanism to genome evolution and codon usage is discussed.

Prediction of signal recognition particle RNA genes

We describe a method for prediction of genes that encode the RNA component of the signal recognition particle (SRP). A heuristic search for the strongly conserved helix 8 motif of SRP RNA is combined with covariance models that are based on previously known SRP RNA sequences. By screening available genomic sequences we have identified a large number of novel SRP RNA genes and we can account for at least one gene in every genome that has been completely sequenced. Novel bacterial RNAs include that of Thermotoga maritima, which, unlike all other non-gram-positive eubacteria, is predicted to have an Alu domain. We have also found the RNAs of Lactococcus lactis and Staphylococcus to have an unusual UGAC tetraloop in helix 8 instead of the normal GNRA sequence. An investigation of yeast RNAs reveals conserved sequence elements of the Alu domain that aid in the analysis of these RNAs. Analysis of the human genome reveals only two likely genes, both on chromosome 14. Our method for SRP RNA gene prediction is the first convenient tool for this task and should be useful in genome annotation.

Investigating the relationship between genome structure, composition, and ecology in prokaryotes

Our thesis is that the DNA composition and structure of genomes are selected in part by mutation bias (GC pressure) and in part by ecology. To illustrate this point, we compare and contrast the oligonucleotide composition and the mosaic structure in 36 complete genomes and in 27 long genomic sequences from archaea and eubacteria. We report the following findings (1) High-GC-content genomes show a large underrepresentation of short distances between G(n) and C(n) homopolymers with respect to distances between A(n) and T(n) homopolymers; we discuss selection versus mutation bias hypotheses. (2) The oligonucleotide compositions of the genomes of Neisseria (meningitidis and gonorrhoea), Helicobacter pylori and Rhodobacter capsulatus are more biased than the other sequenced genomes. (3) The genomes of free-living species or nonchronic pathogens show more mosaic-like structure than genomes of chronic pathogens or intracellular symbionts. (4) Genome mosaicity of intracellular parasites has a maximum corresponding to the average gene length; in the genomes of free-living and nonchronic pathogens the maximum occurs at larger length scales. This suggests that free-living species can incorporate large pieces of DNA from the environment, whereas for intracellular parasites there are recombination events between homologous genes. We discuss the consequences in terms of evolution of genome size. (5) Intracellular symbionts and obligate pathogens show small, but not zero, amount of chromosome mosaicity, suggesting that recombination events occur in these species.

Molecular biology and pathogenicity of mycoplasmas

The recent sequencing of the entire genomes of Mycoplasma genitalium and M. pneumoniae has attracted considerable attention to the molecular biology of mycoplasmas, the smallest self-replicating organisms. It appears that we are now much closer to the goal of defining, in molecular terms, the entire machinery of a self-replicating cell. Comparative genomics based on comparison of the genomic makeup of mycoplasmal genomes with those of other bacteria, has opened new ways of looking at the evolutionary history of the mycoplasmas. There is now solid genetic support for the hypothesis that mycoplasmas have evolved as a branch of gram-positive bacteria by a process of reductive evolution. During this process, the mycoplasmas lost considerable portions of their ancestors' chromosomes but retained the genes essential for life. Thus, the mycoplasmal genomes carry a high percentage of conserved genes, greatly facilitating gene annotation. The significant genome compaction that occurred in mycoplasmas was made possible by adopting a parasitic mode of life. The supply of nutrients from their hosts apparently enabled mycoplasmas to lose, during evolution, the genes for many assimilative processes. During their evolution and adaptation to a parasitic mode of life, the mycoplasmas have developed various genetic systems providing a highly plastic set of variable surface proteins to evade the host immune system. The uniqueness of the mycoplasmal systems is manifested by the presence of highly mutable modules combined with an ability to expand the antigenic repertoire by generating structural alternatives, all compressed into limited genomic sequences. In the absence of a cell wall and a periplasmic space, the majority of surface variable antigens in mycoplasmas are lipoproteins. Apart from providing specific antimycoplasmal defense, the host immune system is also involved in the development of pathogenic lesions and exacerbation of mycoplasma induced diseases. Mycoplasmas are able to stimulate as well as suppress lymphocytes in a nonspecific, polyclonal manner, both in vitro and in vivo. As well as to affecting various subsets of lymphocytes, mycoplasmas and mycoplasma-derived cell components modulate the activities of monocytes/macrophages and NK cells and trigger the production of a wide variety of up-regulating and down-regulating cytokines and chemokines. Mycoplasma-mediated secretion of proinflammatory cytokines, such as tumor necrosis factor alpha, interleukin-1 (IL-1), and IL-6, by macrophages and of up-regulating cytokines by mitogenically stimulated lymphocytes plays a major role in mycoplasma-induced immune system modulation and inflammatory responses.

RNase P RNA structure and cleavage reflect the primary structure of tRNA genes

The function of RNase P RNA depends on its folding in space. A majority of RNase P RNAs from various bacteria show a similar secondary structure to that of Escherichia coli (M1 RNA). However, there are exceptions as exemplified by the RNase P RNA derived from the low GC-content Gram-positive bacteria Bacillus subtilis and Mycoplasma hyopneumoniae (Hyo P RNA). Previous studies using M1 RNA and Hyo P RNA suggest differences both with respect to the kinetics of cleavage as well as to cleavage site recognition. Here we have studied cleavage by these two structurally different RNase P RNAs as a function of changes in the 5' leader and the 3'-terminal CCA motif in the substrate. Our data suggest that the nucleotide at the -2 position in the 5' leader plays a role both for cleavage site recognition and for the rate of cleavage. However, depending on the identity of the -2 residue differences in the cleavage pattern comparing these two types of RNase P RNAs were observed. The results also suggest that the identity of the -1/+73 base-pair in the substrate influences the cleavage site recognition process. These findings will be related to differences in structure comparing these types of RNase P RNAs and the "RCCA-RNase P RNA" interaction. In addition, our findings will be discussed with respect to the primary structure of the tRNA genes in different bacteria.

The comparative metabolism of the mollicutes (Mycoplasmas): the utility for taxonomic classification and the relationship of putative gene annotation and phylogeny to enzymatic function in the smallest free-living cells

Mollicutes or mycoplasmas are a class of wall-less bacteria descended from low G + C% Gram-positive bacteria. Some are exceedingly small, about 0.2 micron in diameter, and are examples of the smallest free-living cells known. Their genomes are equally small; the smallest in Mycoplasma genitalium is sequenced and is 0.58 mb with 475 ORFs, compared with 4.639 mb and 4288 ORFs for Escherichia coli. Because of their size and apparently limited metabolic potential, Mollicutes are models for describing the minimal metabolism necessary to sustain independent life. Mollicutes have no cytochromes or the TCA cycle except for malate dehydrogenase activity. Some uniquely require cholesterol for growth, some require urea and some are anaerobic. They fix CO2 in anaplerotic or replenishing reactions. Some require pyrophosphate not ATP as an energy source for reactions, including the rate-limiting step of glycolysis: 6-phosphofructokinase. They scavenge for nucleic acid precursors and apparently do not synthesize pyrimidines or purines de novo. Some genera uniquely lack dUTPase activity and some species also lack uracil-DNA glycosylase. The absence of the latter two reactions that limit the incorporation of uracil or remove it from DNA may be related to the marked mutability of the Mollicutes and their tachytelic or rapid evolution. Approximately 150 cytoplasmic activities have been identified in these organisms, 225 to 250 are presumed to be present. About 100 of the core reactions are graphically linked in a metabolic map, including glycolysis, pentose phosphate pathway, arginine dihydrolase pathway, transamination, and purine, pyrimidine, and lipid metabolism. Reaction sequences or loci of particular importance are also described: phosphofructokinases, NADH oxidase, thioredoxin complex, deoxyribose-5-phosphate aldolase, and lactate, malate, and glutamate dehydrogenases. Enzymatic activities of the Mollicutes are grouped according to metabolic similarities that are taxonomically discriminating. The arrangements attempt to follow phylogenetic relationships. The relationships of putative gene assignments and enzymatic function in My. genitalium, My. pneumoniae, and My. capricolum subsp. capricolum are specially analyzed. The data are arranged in four tables. One associates gene annotations with congruent reports of the enzymatic activity in these same Mollicutes, and hence confirms the annotations. Another associates putative annotations with reports of the enzyme activity but from different Mollicutes. A third identifies the discrepancies represented by those enzymatic activities found in Mollicutes with sequenced genomes but without any similarly annotated ORF. This suggests that the gene sequence is significantly different from those already deposited in the databanks and putatively annotated with the same function. Another comparison lists those enzymatic activities that are both undetected in Mollicutes and not associated with any ORF. Evidence is presented supporting the theory that there are relatively small gene sequences that code for functional centers of multiple enzymatic activity. This property is seemingly advantageous for an organism with a small genome and perhaps under some coding restraint. The data suggest that a concept of "remnant" or "useless genes" or "useless enzymes" should be considered when examining the relationship of gene annotation and enzymatic function. It also suggests that genes in addition to representing what cells are doing or what they may do, may also identify what they once might have done and may never do again.

A bacterial RNA that functions as both a tRNA and an mRNA

Bacterial tmRNA (transfer-messenger RNA, also known as 10Sa RNA) contains a tRNA-like structure in the 5'- and 3'-end sequences and an internal reading frame encoding a 'tag' peptide. The dual function of this molecule as both a tRNA and an mRNA facilitates a trans-translation reaction, in which a ribosome can switch between translation of a truncated mRNA and the tmRNA's tag sequence. The result is a chimeric protein with the tag peptide attached to the C-terminus of the truncated peptide.

Exploiting unassigned codons in Micrococcus luteus for tRNA-based amino acid mutagenesis

An alternative to suppression of stop codons for the biosynthetic insertion of non-natural amino acids has been developed. Micrococcus luteus , a Gram-positive bacterium, is incapable of translating at least two codons. One of these unused codons was inserted in a gene to act as a nonsense site. An aminoacylated tRNA was synthesized which was complementary to this codon. The gene containing the missing codon was expressed in vitro in a M.luteus transcription/translation system. Read-through of the missing codon occurred only when the complementary tRNA was included. The results demonstrate that M.luteus can be used for incorporation of amino acids via synthetically prepared aminoacylated tRNAs. The use of a M. luteus translation system provides a method for incorporation of non-natural amino acids which avoids the use of stop codons.

Identification of N2-fixing plant- and fungus-associated Azoarcus species by PCR-based genomic fingerprints

Most species of the diazotrophic Proteobacteria Azoarcus spp. occur in association with grass roots, while A. tolulyticus and A. evansii are soil bacteria not associated with a plant host. To facilitate species identification and strain comparison, we developed a protocol for PCR-generated genomic fingerprints, using an automated sequencer for fragment analysis. Commonly used primers targeted to REP (repetitive extragenic palindromic) and ERIC (enterobacterial repetitive intergenic consensus) sequence elements failed to amplify fragments from the two species tested. In contrast, the BOX-PCR assay (targeted to repetitive intergenic sequence elements of Streptococcus) yielded species-specific genomic fingerprints with some strain-specific differences. PCR profiles of an additional PCR assay using primers targeted to tRNA genes (tDNA-PCR, for tRNA(IIe)) were more discriminative, allowing differentiation at species-specific (for two species) or infraspecies-specific level. Our protocol of several consecutive PCR assays consisted of 16S ribosomal DNA (rDNA)-targeted, genus-specific PCR followed by BOX- and tDNA-PCR; it enabled us to assign new diazotrophic isolates originating from fungal resting stages (sclerotia) to known species of Azoarcus. The assignment was confirmed by phylogenetic analysis of 16S rDNA sequences. Additionally, the phylogenetic distances and the lack of monophyly suggested emendment of the genus Azoarcus: the unnamed species Azoarcus groups C and D and a new group (E) of Azoarcus, which was detected in association with fungi, are likely to have the taxonomic rank of three different genera. According to its small subunit rRNA, the sclerotium-forming basidiomycete was related to the Ustilagomycetes, facultatively biotrophic parasites of plants. Since they occurred in a field which was under cultivation with rice and wheat, these fungi might serve as a niche for survival for Azoarcus in the soil and as a source for reinfection of plants.

Codon usage and base composition in Rickettsia prowazekii

Codon usage and base composition in sequences from the A + T-rich genome of Rickettsia prowazekii, a member of the alpha Proteobacteria, have been investigated. Synonymous codon usage patterns are roughly similar among genes, even though the data set includes genes expected to be expressed at very different levels, indicating that translational selection has been ineffective in this species. However, multivariate statistical analysis differentiates genes according to their G + C contents at the first two codon positions. To study this variation, we have compared the amino acid composition patterns of 21 R. prowazekii proteins with that of a homologous set of proteins from Escherichia coli. The analysis shows that individual genes have been affected by biased mutation rates to very different extents: genes encoding proteins highly conserved among other species being the least affected. Overall, protein coding and intergenic spacer regions have G + C content values of 32.5% and 21.4%, respectively. Extrapolation from these values suggests that R. prowazekii has around 800 genes and that 60-70% of the genome may be coding.

RNase P from bacteria. Substrate recognition and function of the protein subunit

RNase P recognizes many different precursor tRNAs as well as other substrates and cleaves all of them accurately at the expected position. RNase P recognizes the tRNA structure of the precursor tRNA by a set of interactions between the catalytic RNA subunit and the T- and acceptor-stems mainly, although residues in the 5'-leader sequence as well as the 3'-terminal CCA are important. These conclusions have been reached by several studies on mutant precursor tRNAs as well as cross-linking studies between RNase P RNA and precursor tRNAs. The protein subunit of RNase P seems also to affect the way that the substrate is recognized as well as the range of substrates that can be used by RNase P, although the protein does not seem to interact directly with the substrates. The interaction between the protein and RNA subunits of RNase P has been extensively studied in vitro. The protein subunit sequence is not highly conserved among bacteria, however different proteins are functionally equivalent as heterologous reconstitution of the RNase P holoenzyme can be achieved in many cases.

Transfer RNA gene organization and RNase P

Mature tRNAs are remarkably similar in all cells. However, the primary transcripts from tRNA genes can vary considerably due to differences in gene organization. RNase P must be able to recognize the elements that are common to all tRNA precursors to accurately remove the 5'-leader sequences.

RNase P RNA of Mycoplasma capricolum

There are at least six small stable RNAs in Mycoplasma capricolum cells besides tRNAs and rRNAs. One of them, MCS5 RNA, is a homolog of RNase P RNA. The predicted secondary structure of this RNA is essentially the same as that of other eubacterial RNase P RNAs. MCS5 RNA is more similar to the RNase P RNA of B. Subtilis than to that of E. coli. This is consistent with previous conclusions that mycoplasmas are phylogenetically related to the low G + C Gram-positive bacterial group. The major substrates for MCS5 RNA must be the precursors of tRNAs. The precursor of MCS6 RNA, which is a homolog of the E. coli 10Sa RNA, may also be a substrate for the MCS5 RNA because this RNA has a tRNA-like structure at its 5' and 3' ends.

Characterization of IS1221 from Mycoplasma hyorhinis: expression of its putative transposase in Escherichia coli incorporates a ribosomal frameshift mechanism

Seven complete and two partial copies of IS1221 variants from Mycoplasma hyorhinis and Mycoplasma hyopneumoniae characterized to date have established a consensus IS1221 as a 1513 bp element with unique structural characteristics resembling the IS3 family of bacterial insertion sequences. Each IS1221 copy contains highly conserved 28 bp imperfect terminal inverted repeats and three distinctive internal inverted repeats (LIR, RIR and IIR). IIR is located within the coding region of the element and it is proposed that it plays a critical role in the regulation of putative transposase expression. Consensus IS1221 and one particular copy, G1135.2, contain a single long open reading frame (ORF). Two potential initiation codons are present at nucleotide 46 (AUG46) and nucleotide 397 (AUG397) and both are preceded by strong ribosome-binding sites. Both initiation codons can be used efficiently in an Escherichia coli T7 expression system. The LIR has a negative regulatory effect on translation initiation from AUG46. A-1 translational frameshift event is shown to be involved in expression of the IS1221 ORF and results in the production of 20 kDa and 6 kDa truncated proteins from the respective upstream initiation codons of the IS1221 ORF. Base substitution and deletion mutations in sequences resembling characterized motifs in documented examples of translational frameshifting resulted in a significant increase in the full-length products and a corresponding decrease in the truncated products from the IS1221 ORF. In contrast to the usual -1 frameshift regulatory event in the IS3 family, which produces a transframe fusion product as the active transposase, IS1221 may have evolved a high-frequency -1 frameshift mechanism that produces a truncated product from the upstream coding domain and thereby results in the regulated low-level production of the full-length presumptive transposase.

tRNA-like structures in 10Sa RNAs of Mycoplasma capricolum and Bacillus subtilis

The stable RNAs, whose sequences are homologous to 10Sa RNA of Escherichia coli, have been isolated from Mycoplasma capricolum and Bacillus subtilis, both belonging to the Gram-positive bacterial group. The total nucleotide sequences of the RNAs have been determined by partial RNA sequencing and DNA sequencing of their genes. A comparison of the sequences, together with those of other bacterial 10Sa RNAs so far known, has shown that the 5'- and 3'-end sequences are well conserved among species, while the central parts reveal little homologies. Unexpectedly, the conserved 5'- and 3'-regions can be folded in a common tRNA-like structure containing an amino acid-acceptor stem and a T phi C-stem/loop. The 3'-terminal CCA sequence of B.subtilis 10Sa RNA is not encoded on the DNA, but is added after transcription. Furthermore, the RNA is aminoacylatable with alanine in vitro, and binds to the 70S ribosome in vivo.

Codon reading scheme in Mycoplasma pneumoniae revealed by the analysis of the complete set of tRNA genes

The 33 genes encoding the complete set of tRNA species in Mycoplasma pneumoniae have been cloned and sequenced. They are organized into 5 clusters in addition to 9 single genes. No redundant gene was found, indicating that 33 tRNAs correspond to 32 different anticodons and decode all 62 codons used in this organism. There is only one single tRNA for each of the Ala, Leu, Pro, and Val family boxes. Therefore, a simplified decoding system resembling that recently described for Mycoplasma capricolum (1) has to also exist in M.pneumoniae. However, analysis of the anticodon set and codon usage revealed features characteristic of the latter: (i) there is no obvious preference toward AT rich synonymous codons, (ii) CGG codons are assigned for arginine and are translated by tRNA Arg(UCG), and (iii) CNN or GNN anticodons are encountered in the Ser, Thr, Arg, and Gly family boxes. We thus propose that this codon-anticodon recognition pattern has emerged in the 'M.pneumoniae cluster' under a genomic economization strategy but without the influence of AT pressure.

Directional mutation pressure, mutator mutations, and dynamics of molecular evolution

Using a general form of the directional mutation theory, this paper analyzes the effect of mutations in mutator genes on the G+C content of DNA, the frequency of substitution mutations, and evolutionary changes (cumulative mutations) under various degrees of selective constraints. Directional mutation theory predicts that when the mutational bias between A/T and G/C nucleotide pairs is equilibrated with the base composition of a neutral set of DNA nucleotides, the mutation frequency per gene will be much lower than the frequency immediately after the mutator mutation takes place. This prediction explains the wide variation of the DNA G+C content among unicellular organisms and possibly also the wide intragenomic heterogeneity of third codon positions for the genes of multicellular eukaryotes. The present analyses lead to several predictions that are not consistent with a number of the frequently held assumptions in the field of molecular evolution, including belief in a constant rate of evolution, symmetric branching of phylogenetic trees, the generality of higher mutation frequency for neutral sets of nucleotides, the notion that mutator mutations are generally deleterious because of their high mutation rates, and teleological explanations of DNA base composition.

Staphylococcus aureus has clustered tRNA genes

The polymerase chain reaction (PCR) was used to detect large tRNA gene clusters in Bacillus subtilis, Bacillus badius, Bacillus megaterium, Lactobacillus brevis, Lactobacillus casei, and Staphylococcus aureus. The primers were based on conserved sequences of known gram-positive bacterial tRNA(Arg) and tRNA(Phe) genes. This PCR procedure detected an unusually large tRNA gene cluster in S. aureus. PCR-generated probes were used to identify a 4.5-kb EcoRI fragment that contained 27 tRNA genes immediately 3' to an rRNA operon. Some of these 27 tRNA genes are very similar, but only 1 is exactly repeated in the cluster. The 5' end of this cluster has a gene order similar to that found in the 9- and 21-tRNA gene clusters of B. subtilis. The 3' end of this S. aureus cluster exhibits more similarity to the 16-tRNA gene cluster of B. subtilis. The 24th, 25th, and 26th tRNA genes of this S. aureus tRNA gene cluster code for three similar, unusual Gly-tRNAs that may be used in the synthesis of the peptidoglycan in the cell wall but not in protein synthesis. Southern analysis of restriction digests of S. aureus DNA indicate that there are five to six rRNA operons in this bacterium's genome and that most or all may have large tRNA gene clusters at the 3' end.

The genes aroA and trnQ are located upstream of psbO in the chromosome of Synechocystis 6803

We have identified the existence of two genes, trnQ and aroA, located upstream of the psbO gene in Synechocystis sp. PCC 6803. The trnQ gene encodes a glutamine-specific transfer RNA (tRNA(Gln)) and the sequence given is the first reported for any cyanobacterium. The gene seems to exist as a single copy since its deletion results in non-viable mutation. The aroA gene encodes for 5-enolpyruvylshikimate 3-phosphate synthase and its discovery in the genome of Synechocystis 6803 is the first genetic evidence for the existence of the shikimate biosynthetic pathway in cyanobacteria. Interestingly, the partial sequence shares close homologies with the sequences of aroA from Gram-positive bacteria.

A small RNA of Mycoplasma capricolum that resembles eukaryotic U6 small nuclear RNA

Mycoplasma capricolum, a parasitic prokaryote, contains several small stable RNAs, besides rRNAs and tRNAs. One of them, designated MCS4 RNA (125 nucleotides in length), has been isolated and sequenced. This RNA is abundant in the cell, and is encoded by two genes. Unexpectedly, MCS4 RNA has been found to reveal extensive sequence similarity to eukaryotic U6 snRNAs. This finding suggests that MCS4 and U6 snRNAs are derived from a common ancestral RNA that has existed before the divergence of prokaryotes and eukaryotes.

Spiroplasma citri UGG and UGA tryptophan codons: sequence of the two tryptophanyl-tRNAs and organization of the corresponding genes

From the total tRNAs of Spiroplasma citri, we isolated and purified two tRNA(Trp) species by using chromatography on an RPC-5 column followed by denaturing polyacrylamide gel electrophoresis. The sequence of the two tRNAs, as well as the sequences of the corresponding genes, were determined. One of the two tRNA(Trp) species has a CCA anticodon and is able to pair with the universal UGG tryptophan codon, while the second has a U*CA (U* is a modified uridine) anticodon and is able to pair with UGA but also with UGG in accordance with the "U:N wobble" rule. Thus, in S. citri, UGA is not a stop codon but codes for tryptophan. The two tRNA(Trp) genes, together with a third tRNA gene, tRNA(Ser) (CGA), belong to a single transcription unit. The nucleotide sequences of the two tRNA(Trp) species show 82.9% similarity. The two spiroplasmal tRNA(Trp) species can be aminoacylated by using an aminoacyl-tRNA synthetase fraction from S. citri. In contrast, the enzyme fraction from Escherichia coli aminoacylates tRNA(Trp) (CCA) but not tRNA(Trp) (U*CA).

A physical and genetic map of the Spiroplasma citri genome

A physical and genetic map of the Spiroplasma citri genome has been constructed using several restriction enzymes and pulsed field gel electrophoresis. A number of genes were subsequently localized on the map by the use of appropriate probes. The genome size of the spiroplasma estimated from restriction fragments is close to 1780 kbp, the largest of all Mollicutes studied so far. It contains multisite insertions of Spiroplasma virus 1 (SpV1) sequences. The physical and genetic map of the S. citri genome shares several features with that of other Mollicutes, especially those in the Mycoplasma mycoides cluster. This supports the finding that S. citri and these Mycoplasma spp. are phylogenetically related.

Length polymorphisms in tRNA intergenic spacers detected by using the polymerase chain reaction can distinguish streptococcal strains and species

Intergenic tRNA spacers from strains of streptococcal groups A, B, and G were amplified by using the polymerase chain reaction (PCR) at low stringency with consensus tRNA gene primers. Cloning and sequencing showed that many of the homologous intergenic spacers differed in length between species. The sequences of the tRNA genes that flank these polymorphic spacers were determined and used to synthesize fully complementary primers. With these primers at high stringency, PCR products which varied in lengths from 53 to 71 bp, depending on the species or strain, were obtained from streptococcal DNAs, even in the presence of a 1,000-fold mass excess of human DNA. PCR products, the lengths of which could also be used for classification, were obtained at high stringency from a few genera closely related to Streptococcus. No products were obtained from genomic DNAs from more distantly related genera. Production of species- or strain-specific tRNA intergenic length polymorphisms with primers that generate characteristic products from a variety of species within the same genus should be applicable to many organisms, including those that would otherwise be difficult to culture or identify.

PCR-amplified length polymorphisms in tRNA intergenic spacers for categorizing staphylococci

The intergenic spacers between some adjacent tRNA genes were shown to be polymorphic in length when closely related Staphylococcus species were compared. A simple procedure was developed to detect and sequence these tRNA intergenic length polymorphisms (tRNA-ILPs). A comparison of homologous tRNA gene sequences flanking these ILPs in three Staphylococcus species was used to derive primers for high-stringency amplification of the ILPs by the polymerase chain reaction (PCR). The detection of tRNA-ILPs by PCR allowed the classification of virtually all strains from the five species of Staphylococcus that were examined. The procedure used to identify, sequence and derive primers for PCR detection of tRNA-ILPs in Staphylococcus should be applicable to many other genera of eubacteria. These primers could be used on uncultured material such as clinical samples.

A cluster of nine tRNA genes between ribosomal gene operons in Bacillus subtilis

A cluster of nine tRNA genes located in the 1-kb region between ribosomal operons rrnJ and rrnW in Bacillus subtilis has been cloned and sequenced. This cluster contains the genes for tRNA(UACVal), tRNA(UGUThr), tRNA(UUULys), tRNA(UAGLeu). tRNA(GCCGly), tRNA(UAALeu), tRNA(ACGArg), tRNA(UGGPro), and tRNA(UGCAla). The newly discovered tRNA gene cluster combines features of the 3'-end of trnI, a cluster of 6 tRNA genes between ribosomal operons rrnI and rrnH, and of the 5'-end of trnB, a cluster of 21 tRNA genes found immediately 3' to rrnB. Neither the tRNA(UAGLeu) gene nor its product has been found previously in B. subtilis. With the discovery of this new set of tRNA genes, a total of 60 such genes have now been found in B. subtilis. These known genes account for almost all of the tRNA hybridizing restriction fragments of the B. subtilis genome. The 60 known tRNA genes of B. subtilis code for only 28 different anticodons, compared with a total of 41 different anticodons for 78 tRNA genes in Escherichia coli. This may indicate that B. subtilis does not need as many anticodons because of more flexible translation rules, similar to the situation in Mycoplasma capricolum.

The acceptor stem in pre-tRNAs determines the cleavage specificity of RNase P

As the result of an unusual RNase P specificity, some special, mature tRNAs have acceptor stems with eight instead of the common seven base pairs. The data from numerous studies suggest that some features in the tRNA domain of pre-tRNAs are important for this behaviour. Here, we show that only five base pairs in the acceptor stem of bacterial histidine tRNAs are required to obtain the changed cleavage site in an unrelated eukaryotic serine tRNA.

Recent evidence for evolution of the genetic code

The genetic code, formerly thought to be frozen, is now known to be in a state of evolution. This was first shown in 1979 by Barrell et al. (G. Barrell, A. T. Bankier, and J. Drouin, Nature [London] 282:189-194, 1979), who found that the universal codons AUA (isoleucine) and UGA (stop) coded for methionine and tryptophan, respectively, in human mitochondria. Subsequent studies have shown that UGA codes for tryptophan in Mycoplasma spp. and in all nonplant mitochondria that have been examined. Universal stop codons UAA and UAG code for glutamine in ciliated protozoa (except Euplotes octacarinatus) and in a green alga, Acetabularia. E. octacarinatus uses UAA for stop and UGA for cysteine. Candida species, which are yeasts, use CUG (leucine) for serine. Other departures from the universal code, all in nonplant mitochondria, are CUN (leucine) for threonine (in yeasts), AAA (lysine) for asparagine (in platyhelminths and echinoderms), UAA (stop) for tyrosine (in planaria), and AGR (arginine) for serine (in several animal orders) and for stop (in vertebrates). We propose that the changes are typically preceded by loss of a codon from all coding sequences in an organism or organelle, often as a result of directional mutation pressure, accompanied by loss of the tRNA that translates the codon. The codon reappears later by conversion of another codon and emergence of a tRNA that translates the reappeared codon with a different assignment. Changes in release factors also contribute to these revised assignments. We also discuss the use of UGA (stop) as a selenocysteine codon and the early history of the code.

Evolution of tRNAs and tRNA genes in Acholeplasma laidlawii

The genes for 22 tRNA species from Acholeplasma laidawii, belonging to the class Mollicutes (Mycoplasmas), have been cloned and sequenced. Sixteen genes are organized in 3 clusters consisting of eleven, three and two tRNA genes, respectively, and the other 6 genes exist as a single gene. The arrangement of tRNA genes in the 11-gene, the 3-gene and the 2-gene clusters reveals extensive similarity to several parts of the 21-tRNA or 16-tRNA gene cluster in Bacillus subtilis. The 11-gene cluster is also similar to the tRNA gene clusters found in other mycoplasma species, the 9-tRNA gene cluster in M.capricolum and in M.mycoides, and the 10-tRNA gene cluster in Spiroplasma meliferm. The results suggest that the tRNA genes in mycoplasmas have evolved from large tRNA gene clusters in the ancestral Gram-positive bacterial genome common to mycoplasmas and B.subtilis. The anticodon sequences including base modifications of 15 tRNA species from A.laidlawii were determined. The anticodon composition and codon-recognition patterns of A.laidlawii resemble those of Bacillus subtilis rather than those of other mycoplasma species.

Levels of tRNAs in bacterial cells as affected by amino acid usage in proteins

Transfer RNAs of Mycoplasma capricolum were separated by two-dimensional polyacrylamide gel electrophoresis, and the relative abundance of each of the 28 known tRNA species was measured. There existed a correlation between the relative amount of isoacceptor tRNAs and the frequency in choosing synonymous codons that could be translated by the isoacceptors. Furthermore, it was observed that the total amount of tRNAs for a particular amino acid was paralleled by the composition of the amino acid in ribosomal proteins. A similar relationship was obtained from reexamination of the previous data on Escherichia coli tRNAs, suggesting that the amount of tRNAs for an amino acid is affected by the usage of the amino acid in proteins.

A random sequencing approach for placing markers on the physical map of Mycoplasma genitalium

A physical map of the Mycoplasma genitalium genome has been prepared using pulsed-field gel electrophoresis. This report details recent efforts made to add markers or specific loci to this map in the absence of any mutants or system of genetic exchange. A total of 44 random clones were partially sequenced. Computer analysis was performed in an attempt to identify homologies with genes already recorded in the DNA sequence database. Clones with a large extent of homology to genes from other microorganisms have been assigned to specific loci on the M. genitalium map by hybridization to selected restriction digests. The additional data has facilitated an updated version of the physical map, and verified this random sequencing method as a useful mapping procedure as well as offering new insight into the physiological processes of this fastidious organism.

Novel anticodon composition of transfer RNAs in Micrococcus luteus, a bacterium with a high genomic G + C content. Correlation with codon usage

The number and relative amount of isoacceptor tRNAs for each amino acid in Micrococcus luteus, a Gram-positive bacterium with high genomic G + C content, have been determined by sequencing their anticodon loop and its adjacent regions and by selective labelling of tRNAs. Thirty-one tRNA species with 29 different anticodon sequences have been detected. All the tRNAs have G or C at the anticodon first position except for tRNA(ICGArg) and tRNA(NGASer), in response to the abundant usage of NNC and NNG codons. No tRNA with the anticodon UNN capable of translating codon NNA has been detected, in accordance with a very low or zero usage of NNA codons. The relative amount of isoacceptor tRNAs for an amino acid determined by selective labelling strongly correlates with usage of the corresponding codons. On the basis of these and other observations in this and other eubacterial species, we conclude that the relative amount and anticodon composition of isoacceptor tRNA species are flexible, and their changes are mainly adaptive phenomena that have been primarily affected by codon usage, which in turn is affected by directional mutation pressure.

CGG: an unassigned or nonsense codon in Mycoplasma capricolum

CGG is an arginine codon in the universal genetic code. We previously reported that in Mycoplasma capricolum, a relative of Gram-positive eubacteria, codon CGG did not appear in coding frames, including termination sites, and tRNA(ArgCCG) pairing with codon CGG, was not detected. These facts suggest that CGG is a nonsense (unassigned and untranslatable) codon--i.e., not assigned to arginine or to any other amino acid. We have investigated whether CGG is really an unassigned codon by using a cell-free translation system prepared from M. capricolum. Translation of synthetic mRNA containing in-frame CGG codons does not result in "read-through" to codons beyond the CGG codons--i.e., translation ceases just before CGG. Sucrose-gradient centrifugation profiles of the reaction mixture have shown that the bulk of peptide that has been synthesized is attached to 70S ribosomes and is released upon further incubation with puromycin. The result suggests that the peptide is in the P site of ribosome in the form of peptidyl-tRNA, leaving the A site empty. When in-frame CGG codons are replaced by UAA codons in mRNA, no read-through occurs beyond UAA, just as in the case of CGG. However, the synthesized peptide is released from 70S ribosomes, presumably by release factor 1. These data suggest strongly that CGG is an unassigned codon and differs from UAA in that CGG is not used for termination.