Spectinomycin operon of Micrococcus luteus: evolutionary implications of organization and novel codon usage

Singular features of trypanosomatids' phosphotransferases involved in cell energy management

Trypanosomatids are responsible for economically important veterinary affections and severe human diseases. In Africa, Trypanosoma brucei causes sleeping sickness or African trypanosomiasis, while in America, Trypanosoma cruzi is the etiological agent of Chagas disease. These parasites have complex life cycles which involve a wide variety of environments with very different compositions, physicochemical properties, and availability of metabolites. As the environment changes there is a need to maintain the nucleoside homeostasis, requiring a quick and regulated response. Most of the enzymes required for energy management are phosphotransferases. These enzymes present a nitrogenous group or a phosphate as acceptors, and the most clear examples are arginine kinase, nucleoside diphosphate kinase, and adenylate kinase. Trypanosoma and Leishmania have the largest number of phosphotransferase isoforms ever found in a single cell; some of them are absent in mammals, suggesting that these enzymes are required in many cellular compartments associated to different biological processes. The presence of such number of phosphotransferases support the hypothesis of the existence of an intracellular enzymatic phosphotransfer network that communicates the spatially separated intracellular ATP consumption and production processes. All these unique features make phosphotransferases a promising start point for rational drug design for the treatment of human trypanosomiasis.

Structural peculiarities of linear megaplasmid, pLMA1, from Micrococcus luteus interfere with pyrosequencing reads assembly

Different strains of Micrococcus luteus, isolated from high-altitude Argentinean wetlands, were recently reported to harbour the linear plasmids pLMA1, pLMH5 and pLMV7, all of which with 5′-covalently attached terminal proteins. The link between pLMA1 and the host’s erythromycin resistance as well as further presumptive qualities prompted us to perform a detailed characterization. When the 454 technology was applied for direct sequencing of gel-purified pLMA1, assembly of the reads was impossible. However, combined Sanger/454 sequencing of cloned pLMA1 fragments, covering altogether 23 kb of the 110-kb spanning plasmid, allowed numerous sequence repeats of varying in lengths to be identified thus rendering an explanation for the above 454 assembly failure. A large number of putative transposase genes were identified as well. Furthermore, a region with five putative iteron sequences is possibly involved in pLMA1 replication.

Comprehensive analysis of prokaryotic mechanosensation genes: their characteristics in codon usage

In the present study, we examined GC nucleotide composition, relative synonymous codon usage (RSCU), effective number of codons (ENC), codon adaptation index (CAI) and gene length for 308 prokaryotic mechanosensitive ion channel (MSC) genes from six evolutionary groups: Euryarchaeota, Actinobacteria, Alphaproteobacteria, Betaproteobacteria, Firmicutes, and Gammaproteobacteria. Results showed that: (1) a wide variation of overrepresentation of nucleotides exists in the MSC genes; (2) codon usage bias varies considerably among the MSC genes; (3) both nucleotide constraint and gene length play an important role in shaping codon usage of the bacterial MSC genes; and (4) synonymous codon usage of prokaryotic MSC genes is phylogenetically conserved. Knowledge of codon usage in prokaryotic MSC genes may benefit from the study of the MSC genes in eukaryotes in which few MSC genes have been identified and functionally analysed.

A 50-kb plasmid rich in mobile gene sequences isolated from a marine micrococcus

A 50,709-bp cryptic plasmid isolated from a marine Micrococcus has been sequenced and found to contain a number of putative mobile genetic elements. The coding regions for 11 putative transposases comprise approximately 17% of the total plasmid sequence. The majority of these transposases are located within a 13-kb cluster which includes a 1553-bp direct repeat consisting of a duplicated pair of transposase genes. The remaining putative ORFs showed similarity to a variety of proteins, the most notable being spider silk.

Translational selection shapes codon usage in the GC-rich genome of Chlamydomonas reinhardtii

In unicellular species codon usage is determined by mutational biases and natural selection. Among prokaryotes, the influence of these factors is different if the genome is skewed towards AT or GC, since in AT-rich organisms translational selection is absent. On the other hand, in AT-rich unicellular eukaryotes the two factors are present. In order to understand if GC-rich genomes display a similar behavior, the case of Chlamydomonas reinhardtii was studied. Since we found that translational selection strongly influences codon usage in this species, we conclude that there is not a common pattern among unicellular organisms.

Cloning and heterologous expression of Entamoeba histolytica adenylate kinase and uridylate/cytidylate kinase

We have isolated two cDNA clones encoding Entamoeba histolytica nucleotide kinases, EhAK and EhUK, expressed them in E. coli and performed functional studies of the recombinant enzymes. Nucleotide sequence analysis showed that EhAK and EhUK genes exhibited the features characteristic of E. histolytica genes, such as transcripts with relatively short 5' and 3' untranslated flanking regions containing the conserved E. histolytica transcription promoter elements located 5' to the initiation codon and a polyadenylation signal in the 3' UTR, a distinctive codon usage bias for A or T in the third position and an AT bias greater than 75% in the flanking regions of the transcripts. At the protein level, both enzymes belong to the short variant nucleoside monophosphate (NMP) kinases, which lack a 29amino acid LID region present in the long variant isoenzymes. EhAK was 30-38% identical to the members of the adenylate kinase (AK) family while EhUK was more similar (48-49% identity) to UMP/CMP kinases. Both enzymes used ATP as preferred phosphate-group donor but each one exhibited strict specificity for the acceptor NMP, EhAK for AMP and EhUK for the pyrimidine nucleoside monophosphates UMP and CMP. Biochemical characterization of the enzymes and phylogenetic reconstruction showed that EhUK is an authentic and well conserved member of the UMP/CMP kinase group while EhAK is the most divergent member known of the AK1 isoenzymes.

Organization of a large gene cluster encoding ribosomal proteins in the cyanobacterium Synechococcus sp. strain PCC 6301: comparison of gene clusters among cyanobacteria, eubacteria and chloroplast genomes

The structure of a large gene cluster containing 22 ribosomal protein (r-protein) genes of the cyanobacterium Synechococcus sp. strain PCC6301 is presented. Based on DNA and protein sequence analyses, genes encoding r-proteins L3, L4, L23, L2, S19, L22, S3, L16, L29, S17, L14, L24, L5, S8, L6, L18, S5, L15, L36, S13, S11, L17, SecY, adenylate kinase (AK) and the alpha subunit of RNA polymerase were identified. The gene order is similar to that of the E. coli S10, spc and alpha operons. Unlike the corresponding E. coli operons, the genes for r-proteins S4, S10, S14 and L30 are not present in this cluster. The organization of Synechococcus r-protein genes also resembles that of chloroplast (cp) r-protein genes of red and brown algal species. This strongly supports the endosymbiotic theory that the cp genome evolved from an ancient photosynthetic bacterium.

Sequence, overproduction and purification of Vibrio proteolyticus ribosomal protein L18 for in vitro and in vivo studies

A strategy suggested by comparative genomic studies was used to amplify the entire Vibrio proteolyticus (Vp) gene for ribosomal protein L18. Vp L18 and its flanking regions were sequenced and compared with the deduced amino acid (aa) sequences of other known L18 proteins. A 26-aa residue segment at the carboxy terminus contains many strongly conserved residues and may be critical for the L18 interaction with 5S rRNA. This approach should allow rapid characterization of L18 from large numbers of bacteria. Both Vp L18 and Escherichia coli (Ec) L18 were overproduced and purified using a T7 expression vector which fuses an N-terminal peptide segment (His-tag) containing 6 histidine residues to the recombinant protein. The purified fusion proteins, Vp His::L18 and Ec His::L18, were both found to bind to either the Vp 5S or Ec 5S rRNAs in vitro. Vp His::L18 protein was also shown to incorporate into Ec ribosomes in vivo. This His-tag strategy likely will have general applicability for the study of ribosomal proteins in vitro and in vivo.

Cloning and sequencing of the secY homolog from Streptomyces lividans 1326

Two conserved regions of SecY proteins from six Gram+ bacteria were exploited in a PCR-based strategy for isolating a secY homolog from Streptomyces lividans (Sl). The nucleotide sequence of part of a 3.8-kb fragment showed that the secY homolog is flanked, at the 5' end, by the gene encoding ribosomal protein L15 and, at the 3' end, by an adenylate kinase-encoding gene. The deduced gene product of secY would have 437 amino acids (aa) and an M(r) of 47,200. Sl SecY shows 89.5, 56.1, 42 and 40% identity to its homologs from Streptomyces scabies, Brevibacterium flavum, Bacillus subtilis and Escherichia coli, respectively. Promoterprobe analyses indicated that the secY gene probably contains its own promoter.

Cloning and characterization of the gene encoding Halobacterium halobium adenylate kinase

The gene (AK) encoding adenylate kinase (AK) of Halobacterium halobium was cloned. AK consisted of 648 bp and coded for 216 amino acids (aa). S1 mapping and primer extension experiments indicated that the transcription start point (tsp) was located immediately upstream from the start codon. The TAT-like promoter sequence was found at a position 20-24 bp upstream from tsp. The most striking property of the enzyme was a putative Zn finger-like structure with four cysteines. It might contribute to the structural stability of the molecule in high-salt conditions. Phylogenetic analysis indicated two lineages of the AK family, the short and long types which diverged a long time ago, possibly before the separation of prokaryotes and eukaryotes. Although the H. halobium AK belongs to the long-type AK lineage, it is located in an intermediary position between the two lineages of the phylogenetic tree, indicating early divergence of the gene along the long-type lineage.

Structural evidence for specific S8-RNA and S8-protein interactions within the 30S ribosomal subunit: ribosomal protein S8 from Bacillus stearothermophilus at 1.9 A resolution

BACKGROUND

Prokaryotic ribosomal protein S8 is an important RNA-binding protein that occupies a central position within the small ribosomal subunit. It interacts extensively with 16S rRNA and is crucial for the correct folding of the central domain of the rRNA. S8 also controls the synthesis of several ribosomal proteins by binding to mRNA. It binds specifically to very similar sites in the two RNA molecules.

RESULTS

S8 is divided into two tightly associated domains and contains three regions that are proposed to interact with other ribosomal components: two potential RNA-binding sites, and a hydrophobic patch that may interact with a complementary hydrophobic region of S5. The N-terminal domain fold is found in several proteins including two that bind double-stranded DNA.

CONCLUSIONS

These multiple RNA-binding sites are consistent with the role of S8 in organizing the central domain and agree with the latest models of the 16S RNA which show that the S8 location coincides with a region of complicated nucleic-acid structure. The presence in a wide variety of proteins of a region homologous to the N-terminal domain supports the idea that ribosomal proteins must represent some of the earliest protein molecules.

Ancient divergence of long and short isoforms of adenylate kinase: molecular evolution of the nucleoside monophosphate kinase family

Adenylate kinases (AK) from vertebrates are separated into three isoforms, AK1, AK2 and AK3, based on structure, subcellular localization and substrate specificity. AK1 is the short type with the amino acid sequence being 27 residues shorter than sequences of the long types, AK2 and AK3. A phylogenetic tree prepared for the AK isozymes and other members of the nucleoside monophosphate (NMP) kinase family shows that the divergence of long and short types occurred first and then differentiation in subcellular localization or substrate specificity took place. The first step involved a drastic change in the three-dimensional structure of the LID domain. The second step was caused mainly by smaller changes in amino acid sequences.

Purification of a cytochrome bd terminal oxidase encoded by the Escherichia coli app locus from a delta cyo delta cyd strain complemented by genes from Bacillus firmus OF4

Escherichia coli GK100, with deletions in the operons encoding its two terminal oxidases, cytochrome bo and ctyochrome bd, was complemented for growth on succinate by a recombinant plasmid (pMS100) containing a 3.4-kb region of DNA from alkaliphilic Bacillus firmus OF4. The complementing DNA was predicted to encode five proteins, but neither sequence analysis nor complementation experiments with subclones provided insight into the basis for the complementation. Cytochrome difference spectra of everted membrane vesicles from the transformed strain had characteristics of a cytochrome bd spectrum but with features different from those observed for alkaliphile membranes. To determine the bacterial source and identity of the structural genes for the cytochrome bd in the transformed mutant, the complex was extracted and partially purified. On sodium dodecyl sulfate-polyacrylamide gels, two polypeptides were resolved from the preparation, 43 (subunit I) and 27 (subunit II) kDa. An internal peptide from subunit I was sequenced, and it yielded the same primary sequence as is found in positions 496 to 510 of E. coli appC. Consistent with the microsequencing results pMS100 failed to complement a triple mutant of E. coli carrying a deletion in appB as well as in the cyo and cyd loci. The deduced sequence of AppBC had been predicted to be very similar to the sequence of CydAB (J. Dassa et al., Mol. Gen. Genet. 229:341-352, 1991) but this is the first demonstration that the former is indeed a cytochrome bd terminal oxidase. The enzyme catalyzed oxygen uptake coupled to quinol or N,N,N',N'-tetramethyl-p-phenylenediamine oxidation, and the activity was sensitive to cyanide. No cross-reactivity to subunit-specific polyclonal antibodies directed against the two individual subunits of cyd-encoded cytochrome bd was detected. Since this is the second cytochrome bd discovered in E. coli, it is proposed that the two complexes be designated cytochrome bd-I (cydAB-encoded enzyme) and cytochrome bd-II (appBC-encoded enzyme). In addition, cbdAB is suggested as a more appropriate gene designation for cytochrome bd than either appBC or cyxAB.

Genetic and transcriptional organization of the Bacillus subtilis spc-alpha region

We used chromosomal walking methods to isolate a 10.8-kb region from the major ribosomal protein (r-protein) gene cluster of Bacillus subtilis (Bs). The gene order in this region, given by gene product, was r-proteins L16-L29-S17-L14-L24-L5-S14-S8-L6-L18-S5-L30-L15-SecY-adenylate kinase (Adk)-methionine aminopeptidase (Map)-initiation factor 1 (IF1)-L36-S13-S11-alpha subunit of RNA polymerase-L17. The region cloned, therefore, contains the homologues for the last three genes of the Escherichia coli (Ec) S10 operon, together with entire spc and alpha operons. This Bs organization differs from the corresponding region in Ec by the inclusion of the genes encoding Adk, Map and IF1 between the genes encoding SecY and L36. Plasmid integration experiments indicated that all 22 genes comprise a single large transcriptional unit controlled from a major promoter which lies upstream from the gene encoding r-protein L16. Promoter probe experiments located lesser activities internal to this large transcriptional unit, the secY and map promoters. The secY promoter region (psecY) contained two activities, each principally functioning in the stationary growth phase when high protein export is required. Thus, the Bs S10-spc-alpha region differs from its Ec counterpart in both genetic and transcriptional organization. Given this difference in transcriptional organization, the mechanisms coordinating expression of the translational apparatus are also likely to differ between Ec and Bs.

Characterization of an unusual Rho factor from the high G + C gram-positive bacterium Micrococcus luteus

A transcription termination factor (Rho) was purified from the Gram-positive bacterium Micrococcus luteus, and the complete gene sequence was determined. The M. luteus Rho polypeptide has 690 residues, which is 271 residues more than its homolog from Escherichia coli. Most of the additional residues compose a highly charged, hydrophilic segment that is inserted in a non-conserved region between two conserved regions of the RNA-binding domain of the known Rho homolog proteins. This segment extends from residues 49 to 311 and includes a stretch of 238 residues that contain no hydrophobic side chains. Biochemical studies indicate that the M. luteus protein is very similar to E. coli Rho in terms of its RNA-dependent NTPase activity and its sensitivity to the Rho-specific inhibitor bicyclomycin. However, the M. luteus protein has a less stringent RNA cofactor specificity. It also acts to terminate RNA transcription with E. coli RNA polymerase on the lambda cro DNA template, but at much earlier termination stop points than those recognized by E. coli Rho. Thus, the M. luteus protein functions as a true Rho factor, but with a different specificity than that of E. coli Rho. We propose that this altered specificity is consistent with its need to function on transcripts that have a high content of G + C residues.

The thiostrepton-resistance-encoding gene in Streptomyces laurentii is located within a cluster of ribosomal protein operons

A common approach to identify and clone biosynthetic gene from an antibiotic-producing streptomycete is to clone the resistance gene for the antibiotic of interest and then use that gene to clone DNA that is linked to it. As a first step toward cloning the genes responsible for the biosynthesis of thiostrepton (Th) in Streptomyces laurentii (Sl), the Th resistance-encoding gene (tsnR) was cloned as a 1.5-kb BamHI-PvuII fragment in Escherichia coli (Ec), and shown to confer Th resistance when introduced into S. lividans TK24. The tsnR-containing DNA fragment was used as a probe to isolate clones from cosmid libraries of DNA in the Ec cosmid vector SuperCos, and pOJ446 (an Ec/streptomycete) cosmid vector. Sequence and genetic analysis of the DNA flanking the tsnR indicates that the Sl tsnR is not closely linked to biosynthetic genes. Instead it is located within a cluster of ribosomal protein operons.

Purification and cloning of Micrococcus luteus ultraviolet endonuclease, an N-glycosylase/abasic lyase that proceeds via an imino enzyme-DNA intermediate

Although Micrococcus luteus UV endonuclease has been reported to be an 18-kDa enzyme with possible homology to the 16-kDa endonuclease V from bacteriophage T4 (Gordon, L. K., and Haseltine, W. A. (1980) J. Biol. Chem. 255, 12047-12050; Grafstrom, R. H., Park, L., and Grossman, L. (1982) J. Biol. Chem. 257, 13465-13474), this study describes three independent purification schemes in which M. luteus UV damage-specific or pyrimidine dimer-specific nicking activity was associated with two proteins of apparent molecular masses of 31 and 32 kDa. An 18-kDa contaminant copurified with the doublet through many of the chromatographic steps, but it was determined to be a homolog of Escherichia coli ribosomal protein L6. Edman degradation analyses of the active proteins yielded identical NH2-terminal amino acid sequences. The corresponding gene (pdg, pyrimidine dimer glycosylase) was cloned. The protein bears strong sequence similarities to the E. coli repair proteins endonuclease III and MutY. Nonetheless, traditionally purified M. luteus protein acted exclusively on cis-syn thymine dimers; it was unable to cleave site-specific oligonucleotide substrates containing a trans-syn -I, (6-4), or Dewar thymine dimer, a 5,6-dihydrouracil lesion, or an A:G or A:C mismatch. The UV endonuclease incised cis-syn dimer-containing DNA in a dose-dependent manner and exhibited linear kinetics within that dose range. Enzyme activity was inhibited by the presence of NaCN or NaBH4 with NaBH4 additionally being able to trap a covalent enzyme-substrate product. These last findings confirm that the catalytic mechanism of M. luteus UV endonuclease, like those of other glycosylase/AP lyases, involves an imino intermediate.

Cloning and sequencing of a secY homolog from Streptomyces scabies

The complete DNA sequence of the Streptomyces scabies (Ss) secY homolog and partial sequences of adjacent upstream and downstream open reading frames (ORFs) have been determined. The nucleotide sequence of a 2-kb region predicts a polypeptide of 437 amino acids in length with homology to the SecY protein family. The Ss secY homolog lies upstream from a sequence that has homology to the adenylate kinase gene (adk) family. The translational stop codon of the putative SecY ORF overlaps the predicted start codon for the Adk ORF. Another ORF that lies upstream from the secY homolog has sequence similarity to the genes that code for the L15 r-protein. Within the 243-bp intergenic region between the L15 and SecY coding sequences, the presence of a streptomycete-like promoter sequence and an 18-bp inverted repeat suggests that the secY homolog and the adjacent downstream sequences may be transcribed independently of the L15 coding sequence. Transcript analysis indicates that the secY homolog is expressed in both Ss and Streptomyces lividans. The proposed gene and transcript organization of the L15-SecY-Adk coding regions in the Ss clone resembles that of Micrococcus luteus which, like the streptomycetes, has a G+C-rich genome.

Novel protein families in archaean genomes

In a quest for novel functions in archaea, all archaean hypothetical open reading frames (ORFs), as annotated in the Swiss-Prot protein sequence database, were used to search the latest databases for the identification of characterized homologues. Of the 95 hypothetical archaean ORFs, 25 were found to be homologous to another hypothetical archaean ORF, while 36 were homologous to non-archaean proteins, of which as many as 30 were homologous to a characterized protein family. Thus the level of sequence similarity in this set reaches 64%, while the level of function assignment is only 32%. Of the ORFs with predicted functions, 12 homologies are reported here for the first time and represent nine new functions and one gene duplication at an acetyl-coA synthetase locus. The novel functions include components of the transcriptional and translational apparatus, such as ribosomal proteins, modification enzymes and a translation initiation factor. In addition, new enzymes are identified in archaea, such as cobyric acid synthase, dCTP deaminase and the first archaean homologues of a new subclass of ATP binding proteins found in fungi. Finally, it is shown that the putative laminin receptor family of eukaryotes and an archaean homologue belong to the previously characterized ribosomal protein family S2 from eubacteria. From the present and previous work, the major implication is that archaea seem to have a mode of expression of genetic information rather similar to eukaryotes, while eubacteria may have proceeded into unique ways of transcription and translation. In addition, with the detection of proteins in various metabolic and genetic processes in archaea, we can further predict the presence of additional proteins involved in these processes.

Sequencing analysis reveals a unique gene organization in the gyrB region of Mycoplasma hominis

The homolog of the gyrB gene, which has been reported to be present in the vicinity of the initiation site of replication in bacteria, was mapped on the Mycoplasma hominis genome, and the region was subsequently sequenced. Five open reading frames were identified flanking the gyrB gene, one of which showed similarity to that which encodes the LicA protein of Haemophilus influenzae. The organization of the genes in the region showed no resemblance to that in the corresponding regions of other bacteria sequenced so far. The gyrA gene was mapped 35 kb downstream from the gyrB gene.

Analysis of directional mutation pressure and nucleotide content in mitochondrial cytochrome b genes

We present a new approach for analyzing directional mutation pressure and nucleotide content in protein-coding genes. Directional mutation pressure, the heterogenicity in the likelihood of different nucleotide substitutions, is used to explain the increasing or decreasing guanine-cytosine content (GC%) in DNA and is represented by microD, in agreement with Sueoka (1962, Proc Natl Acad Sci USA 48:582-592). The new method uses simulation to facilitate identification of significant A+T or G+C pressure as well as the comparison of directional mutation pressure among genes, even when they are translated by different genetic codes. We use the method to analyze the evolution of directional mutation pressure and nucleotide content of mitochondrial cytochrome b genes. Results from a survey of 110 taxa indicate that the cytochrome b genes of most taxa are subjected to significant directional mutation pressure and that the gene is subject to A+T pressure in most cases. Only in the anseriform bird Cairina moschata is the cytochrome b gene subject to significant G+C pressure. The GC% at nonsynonymous codon sites decreases proportionately with increasing A+T pressure, and with a slope less than one, indicating a presence of selective constraints. The cytochrome b genes of insects, nematodes, and eumycotes are subject to extreme A+T pressures (microD = 0.123, 0.224, and 0.130) and, in parallel, the GC% of the nonsynonymous codon sites has decreased from about 0.44 in organisms that are not subjected to A+T or G+C pressure to about 0.332, 0.323, and 0.367, respectively. The distribution of taxa according to the GC% at nonsynonymous codon sites and directional mutation pressure supports the notion that variation in these parameters is a phylogenetic component.

Primary structure of the hydrogenosomal adenylate kinase of Trichomonas vaginalis and its phylogenetic relationships

Hydrogenosomal adenylate kinase of the amitochondriate protist, Trichomonas vaginalis, has been purified and the sequence of its 39 amino-terminal residues established. Based on this sequence and a conserved internal region of the enzyme, a probe was obtained by DNA polymerase chain reaction and used to isolate a genomic DNA clone containing the gene of this enzyme. This gene exists probably as a single copy in T. vaginalis and is not interrupted by introns. The open reading frame obtained codes for a large type adenylate kinase with a mature molecular mass of 24.5 kDa. The T. vaginalis enzyme is homologous with adenylate kinases of other eukaryotes and eubacteria. Strongly conserved parts and residues of the molecule are conserved also in this enzyme. Phylogenetic trees obtained with various methods placed the T. vaginalis adenylate kinase close to the point where the different subfamilies of this enzyme branch from each other, indicating that the T. vaginalis enzyme has no close relationship to any of these subfamilies and that it separated early from other adenylate kinases. The conceptual translation predicts the existence of an amino-terminal nonapeptide absent from the protein purified from hydrogenosomes, similar to the processed amino-terminal extensions of other hydrogenosomal proteins. These extensions have been considered as putative targeting and import signals.

Phylogenetic analysis of sequences from diverse bacteria with homology to the Escherichia coli rho gene

Genes from Pseudomonas fluorescens, Chromatium vinosum, Micrococcus luteus, Deinococcus radiodurans, and Thermotoga maritima with homology to the Escherichia coli rho gene were cloned and sequenced, and their sequences were compared with other available sequences. The species for all of the compared sequences are members of five bacterial phyla, including Thermotogales, the most deeply diverged phylum. This suggests that a rho-like gene is ubiquitous in the Bacteria and was present in their common ancestor. The comparative analysis revealed that the Rho homologs are highly conserved, exhibiting a minimum identity of 50% of their amino acid residues in pairwise comparisons. The ATP-binding domain had a particularly high degree of conservation, consisting of some blocks with sequences of residues that are very similar to segments of the alpha and beta subunits of F1-ATPase and of other blocks with sequences that are unique to Rho. The RNA-binding domain is more diverged than the ATP-binding domain. However, one of its most highly conserved segments includes a RNP1-like sequence, which is known to be involved in RNA binding. Overall, the degree of similarity is lowest in the first 50 residues (the first half of the RNA-binding domain), in the putative connector region between the RNA-binding and the ATP-binding domains, and in the last 50 residues of the polypeptide. Since functionally defective mutants for E. coli Rho exist in all three of these segments, they represent important parts of Rho that have undergone adaptive evolution.

Evolutionary changes in the genetic code

1. The genetic code was thought to be identical ("universal") in all biological systems until 1981, when it was discovered that the coding system in mammalian mitochondria differed from the universal code in the use of codons AUA, UGA, AGA and AGG. 2. Many other differences have since been discovered, some in mitochondria of various phyla, others in bacteria, ciliated protozoa, algae and yeasts. 3. The original thesis that the code was universal and "frozen" depended on the precept that any mutational change in the code would be lethal, because it would produce widespread alterations in the amino acid sequences of proteins. Such changes would destroy protein function, and hence would be intolerable. 4. The objection was "by-passed" by nature. It is possible for a codon to disappear from mRNA molecules, often as a result of directional mutation pressure in DNA: thus all UGA stop codons can be replaced by UAA. 5. The missing UGA codon can then reappear when some UGG tryptophan codons mutate to UGA. The new UGA codons will be translated as tryptophan, as is the case in non-plant mitochondria and Mycoplasma. Therefore, no changes have taken place in the amino acid sequences of proteins. 6. Variations of this procedure have occurred, affecting various codons, and discoveries are still being made. The findings illustrate the evolutionary interplay between tRNA, release factors and codon-anticodon pairing.

Cloning and molecular characterization of the secY genes from Bacillus licheniformis and Staphylococcus carnosus: comparative analysis of nine members of the SecY family

SecY is a central component of the export machinery that mediates the translocation of secretory proteins across the plasma membrane of Escherichia coli. We have cloned and sequenced the secY genes from Bacillus licheniformis and Staphylococcus carnosus. The deduced amino acid sequences are highly homologous to those of other known SecY polypeptides, all having the potential to form 10 transmembrane segments. Comparative analysis of 9 SecY polypeptides, derived from different bacteria, revealed that 14 amino acid positions (2.7%) are identical in all SecY proteins and 89 (16.9%) show conservative changes. Clusters of conserved amino acid residues were found in 4 of the 10 transmembrane segments and 2 of the 6 cytoplasmic domains. It is suggested that the conserved regions might be involved in the translocation activity of SecY or might be required for the correct interaction of SecY with other components of the secretion apparatus.

Protein secretion in Gram-positive bacteria

Gram-positive bacteria often secrete large amounts of proteins into the surrounding medium. This feature makes them attractive as hosts for the industrial production of extracellular enzymes. Compared to Escherichia coli, relatively little is known about the mechanism of protein secretion in these organisms. However, the recent identification of Bacillus subtilis genes whose gene products are highly homologous to some of the Sec (secretion) proteins of E. coli strongly suggests that important principles of protein translocation across the plasma membrane might be highly conserved. In contrast, the steps following the actual translocation event might be different in Gram-positive and Gram-negative bacteria. The scope of this review is to outline the recent progress that has been made in the elucidation of the secretion pathway in Gram-positive bacteria and to discuss potential applications in strain improvement for the industrial production of extracellular proteins.

Codon usage in the G+C-rich Streptomyces genome

The codon usage (CU) patterns of 64 genes from the Gram+ prokaryotic genus Streptomyces were analysed. Despite the extremely high overall G+C content of the Streptomyces genome (estimated at 0.74), individual genes varied in G+C content from 0.610 to 0.797, and had third codon position G+C contents (GC3s) that varied from 0.764 to 0.983. The variation in GC3s explains a significant proportion of the variation in CU patterns. This is consistent with an evolutionary model of the Streptomyces genome where biased mutation pressure has led to a high average G+C content with random variation about the mean, although the variation observed is greater than that expected from a simple binomial model. The only gene in the sample that can be confidently predicted to be highly expressed, EF-Tu of Streptomyces coelicolor A3(2) (GC3s = 0.927), shows a preference for a third position C in several of the four codon families, and for CGY and GGY for Arg and Gly codons, respectively (Y = pyrimidine); similar CU patterns are found in highly expressed genes of the G+C-rich Micrococcus luteus genome. It thus appears that codon usage in Streptomyces is determined predominantly by mutation bias, with weak translational selection operating only in highly expressed genes. We discuss the possible consequences of the extreme codon bias of Streptomyces and consider how it may have evolved. A set of CU tables is provided for use with computer programs that locate protein-coding regions.

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.

Mutagenesis of ribosomal protein S8 from Escherichia coli: defects in regulation of the spc operon

The structural features of Escherichia coli ribosomal protein S8 that are involved in translational regulation of spc operon expression and, therefore, in its interaction with RNA have been investigated by use of a genetic approach. The rpsH gene, which encodes protein S8, was first inserted into an expression vector under the control of the lac promoter and subsequently mutagenized with methoxylamine or nitrous acid. A screening procedure based on the regulatory role of S8 was used to identify mutants that were potentially defective in their ability to associate with spc operon mRNA and, by inference, 16S mRNA. In this way, we isolated 39 variants of the S8 gene containing alterations at 34 different sites, including 37 that led to single amino acid substitutions and 2 that generated premature termination codons. As the mutations were distributed throughout the polypeptide chain, our results indicate that amino acid residues important for the structural integrity of the RNA-binding domain are not localized to a single segment. Nonetheless, the majority were located within three short sequences at the N terminus, middle, and C terminus that are phylogenetically conserved among all known eubacterial and chloroplast versions of this protein. We conclude that these sites encompass the main structural determinants required for the interaction of protein S8 with RNA.

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.

Presence of a gene in the archaebacterium Methanococcus vannielii homologous to secY of eubacteria

The nucleotide sequence of a gene located at the promoter-distal side of the 'spectinomycin-operon' homologue of the archaebacterium Methanococcus vannielii was determined. Its derived amino acid sequence displayed 20% (identical positions) or 52% (including conservative exchanges) similarity, respectively, to SECY from E coli. An alignment of the Methanococcus SECY with eubacterial SECY sequences showed the existence of 10 membrane-associated primary structure domains in equivalent positions. The 5' and 3' ends of the secY transcript were mapped and the gene was expressed in the T7 promoter/polymerase system in E col. The temperature-sensitive growth of the E coli mutant IQ292 which harbours a secYts mutation could be complemented by the secY gene from Methanococcus. This indicates that a protein integral to an archaebacterial ether-lipid membrane can be inserted into a eubacterial phospholipid membrane without apparent loss of function.

Protein synthesis in vitro by Micrococcus luteus

Bacillus subtilis and related gram-positive bacteria which have low to moderate genomic G + C contents are unable to efficiently translate mRNA derived from gram-negative bacteria, whereas Escherichia coli and other gram-negative bacteria are able to translate mRNA from both types of organisms. This phenomenon has been termed translational species specificity. Ribosomes from the low-G + C-content group (low-G + C group) of gram-positive organisms (B. subtilis and relatives) lack an equivalent to Escherichia ribosomal protein S1. The requirement for S1 for translation in E. coli (G. van Dieijen, P. H. van Knippenberg, J. van Duin, B. Koekman, and P. H. Pouwels, Mol. Gen. Genet. 153:75-80, 1977) and its specific role (A.R. Subramanian, Trends Biochem. Sci. 9:491-494, 1984) have been proposed. The group of gram-positive bacteria characterized by high genomic G + C content (formerly Actinomyces species and relatives) contain S1, in contrast to the low-G + C group (K. Mikulik, J. Smardova, A. Jiranova, and P. Branny, Eur. J. Biochem. 155:557-563, 1986). It is not known whether members of the high-G + C group are translationally specific, although there is evidence that one genus, Streptomyces, can express Escherichia genes in vivo (M. J. Bibb and S. N. Cohen, Mol. Gen. Genet. 187:265-277, 1985; J. L. Schottel, M. J. Bibb, and S. N. Cohen, J. Bacteriol. 146:360-368, 1981). In order to determine whether the organisms of this group are translationally specific, we examined the in vitro translational characteristics of a member of the high-G + C group, Micrococcus luteus, whose genomic G + C content is 73%. A semipurified coupled transcription-translation system of M. luteus translates Escherichia mRNA as well as Bacillus and Micrococcus mRNA. Therefore, M. luteus is translationally nonspecific and resembles E. coli rather than B. subtilis in its translational characteristics.

Structure and evolution of the Tetrahymena thermophila gene encoding ribosomal protein L21

The single-copy gene encoding ribosomal protein (r-protein) L21 in the macronucleus of the ciliate Tetrahymena thermophila was cloned and characterized. Sequencing of the L21 gene and a corresponding cDNA clone showed the gene to contain three introns, all located in the 3' half of the coding region. Primer extension and nuclease protection analyses revealed five transcription start points (tsp) 56-73 nucleotides upstream from the start codon. The uppermost tsp mapped to the first T in a sequence, TATAA, often found at tsp in T. thermophila. A comparison of amino acid sequences of r-proteins revealed that T. thermophila L21 belongs to a family of r-proteins that have been conserved in eubacteria, archaebacteria and eukaryotes. On the basis of structural and functional considerations, we suggest that the eukaryotic, like the prokaryotic, members of this protein family interact with the 5S RNA complex in ribosomes.

Prokaryotic genetic code

The prokaryotic genetic code has been influenced by directional mutation pressure (GC/AT pressure) that has been exerted on the entire genome. This pressure affects the synonymous codon choice, the amino acid composition of proteins and tRNA anticodons. Unassigned codons would have been produced in bacteria with extremely high GC or AT genomes by deleting certain codons and the corresponding tRNAs. A high AT pressure together with genomic economization led to a change in assignment of the UGA codon, from stop to tryptophan, in Mycoplasma.

Nucleotide sequence of four genes encoding ribosomal proteins from the 'S10 and spectinomycin' operon equivalent region in the archaebacterium Halobacterium marismortui

Four genes encoding ribosomal proteins HmaS17, HmaL14, HmaL24 and HS3, have been identified in the lambda EMBL3 clone PP*7 from a genomic library of the archaebacterium Halobacterium marismortui. The clone contains genes from the 'S10 and spectinomycin' operon equivalent region. Three of the deduced proteins are homologous to the corresponding Escherichia coli and Methancoccus vannielii S17, L14 and L24 proteins, as well as to eukaryotic proteins from rat or yeast. HS3 was identified as an extra protein corresponding to the gene product for orfc in M. vannielii and the eukaryotic ribosomal protein RS4 from rat. The equivalence of HmaL24 (HL16) and E. coli L24, which share only 28% identical amino acid residues, could now be shown by localizing the HmaL24 gene at the same position in the cluster.

Role of GC-biased mutation pressure on synonymous codon choice in Micrococcus luteus, a bacterium with a high genomic GC-content

The GC (G + C, or G or C)-contents of codon silent positions in all two-codon sets and three codons AUY/A (IIe), and in most of the family boxes of Micrococcus luteus (genomic GC-content: 74%) are 95% to 100% in both the highly and weakly expressed genes. In some family boxes, there is a decrease in NNC codons and an increase in NNG codons from the highly expressed to weakly expressed genes without apparent involvement of NNU and NNA codons. From these observations, we conclude that the selective use of synonymous codons in M. luteus may be largely determined by GC-biased mutation pressure and that in the highly expressed genes tRNAs would act as a weak selection pressure in some family boxes. Available data suggest that the effect of selection pressure by tRNAs on the synonymous codon choice becomes more apparent in the highly expressed genes in eubacteria with intermediate GC-contents such as Escherichia coli and Bacillus subtilis, and that the U/C ratio of the codon third positions in NNU/C-type two-codon sets in the weakly expressed genes would represent the approximate magnitude of directional mutation pressure throughout eubacteria.