Organization and codon usage of the streptomycin operon in Micrococcus luteus, a bacterium with a high genomic G + C content.

The DNA sequence of the Micrococcus luteus str operon, which includes genes for ribosomal proteins S12 (str or rpsL) and S7 (rpsG) and elongation factors (EF) G (fus) and Tu (tuf), has been determined and compared with the corresponding sequence of Escherichia coli to estimate the effect of high genomic G + C content (74%) of M. luteus on the codon usage pattern. The gene organization in this operon and the deduced amino acid sequence of each corresponding protein are well conserved between the two species. The mean G + C content of the M. luteus str operon is 67%, which is much higher than that of E. coli (51%). The codon usage pattern of M. luteus is very different from that of E. coli and extremely biased to the use of G and C in silent positions. About 95% (1,309 of 1,382) of codons have G or C at the third position. Codon GUG is used for initiation of S12, EF-G, and EF-Tu, and AUG is used only in S7, whereas GUG initiates only one of the EF-Tu's in E. coli. UGA is the predominant termination codon in M. luteus, in contrast to UAA in E. coli.     

Analysis of the nucleotide sequence of the Mycoplasma hominis tuf gene and its flanking region

A gene from Mycoplasma hominis PG21 similar to the tuf gene encoding the elongation factor Tu (EF-Tu) of Escherichia coli was cloned and sequenced. The 1193-bp open reading frame flanked by a putative promoter and a potential stem-and-loop structure encoded a 44-kDa polypeptide. The tuf gene of M. hominis PG21 has the lowest G + C content seen in prokaryotes (38.2%). A gene (mhlmp1) encoding a variable surface exposed membrane protein (LMP1) was found downstream the 3' end of the tuf gene. It was found that the highly conserved tuf gene was linked to the highly variable mhlmp1 gene in 26 different M. hominis strains.     

Evidence for horizontal gene transfer in evolution of elongation factor Tu in enterococci

The elongation factor Tu, encoded by tuf genes, is a GTP binding protein that plays a central role in protein synthesis. One to three tuf genes per genome are present, depending on the bacterial species. Most low-G+C-content gram-positive bacteria carry only one tuf gene. We have designed degenerate PCR primers derived from consensus sequences of the tuf gene to amplify partial tuf sequences from 17 enterococcal species and other phylogenetically related species. The amplified DNA fragments were sequenced either by direct sequencing or by sequencing cloned inserts containing putative amplicons. Two different tuf genes (tufA and tufB) were found in 11 enterococcal species, including Enterococcus avium, Enterococcus casseliflavus, Enterococcus dispar, Enterococcus durans, Enterococcus faecium, Enterococcus gallinarum, Enterococcus hirae, Enterococcus malodoratus, Enterococcus mundtii, Enterococcus pseudoavium, and Enterococcus raffinosus. For the other six enterococcal species (Enterococcus cecorum, Enterococcus columbae, Enterococcus faecalis, Enterococcus sulfureus, Enterococcus saccharolyticus, and Enterococcus solitarius), only the tufA gene was present. Based on 16S rRNA gene sequence analysis, the 11 species having two tuf genes all have a common ancestor, while the six species having only one copy diverged from the enterococcal lineage before that common ancestor. The presence of one or two copies of the tuf gene in enterococci was confirmed by Southern hybridization. Phylogenetic analysis of tuf sequences demonstrated that the enterococcal tufA gene branches with the Bacillus, Listeria, and Staphylococcus genera, while the enterococcal tufB gene clusters with the genera Streptococcus and Lactococcus. Primary structure analysis showed that four amino acid residues encoded within the sequenced regions are conserved and unique to the enterococcal tufB genes and the tuf genes of streptococci and Lactococcus lactis. The data suggest that an ancestral streptococcus or a streptococcus-related species may have horizontally transferred a tuf gene to the common ancestor of the 11 enterococcal species which now carry two tuf genes.     

Elongation factor (EF-Tu) gene probe detects polymorphism in Mycoplasma strains

Abstract Polymorphism in mycoplasma strains was observed by Southern blot hybridization of the digested mycoplasma DNAs with the elongation factor (EF-Tu) gene tuf of Escherichia coli . The hybridization patterns revealed genotypic heterogeneity among Mycoplasma gallisepticum strains and a remarkable degree of homogeneity among Mucoplasma pneumoniae strains isolated from pneumonia patients. The distinction among M. gallisepticum strain clusters achieved by the tuf gene probe corresponded exactly with that obtained with the rRNA gene probe pMC5. The tuf gene probe may thus be added as another effective tool in the taxonomy of Mollicutes and in epidemiological surveys of mycoplasma infections.     

Molecular cloning and sequence determination of the tuf gene coding for the elongation factor Tu of Thermus thermophilus HB8

The tuf gene, which encodes the elongation factor Tu (EF-Tu) of Thermus thermophilus HB8, and its flanking regions were cloned and sequenced. The gene encoding EF-G was found upstream of the 5' end of the tuf gene. The tuf gene of T. thermophilus HB8 had a very high G + C content and 84.5% of the third base in codon usage was either G or C. The deduced primary structure of the EF-Tu was composed of 405 amino acid residues with a Mr = 44658. A comparison of the amino acid sequence of EF-Tu from T. thermophilus HB8 with those of Escherichia coli and Saccharomyces cerevisiae mitochondria showed a very high sequence homology (65-70%). Two Cys residues out of the three found in E. coli EF-Tu had been replaced with Val in T. thermophilus HB8 EF-Tu. An extra amino acid sequence of ten residues, consisting predominantly of basic amino acids (Met-182-Gly-191), which does not occur in EF-Tu of E. coli, was found in T. thermophilus HB8.     

Essential factors determining codon usage in ubiquitin genes

Summary Ubiquitin is ubiquitous in all eukaryotes and its amino acid sequence shows extreme conservation. Ubiquitin genes comprise direct repeats of the ubiquitin coding unit with no spacers. The nucleotide sequences coding for 13 ubiquitin genes from 11 species reported so far have been compiled and analyzed. The G+C content of codon third base reveals a positive linear correlation with the genome G+C content of the corresponding species. The slope strongly suggests that the overall G+C content of codons of polyubiquitin genes clearly reflects the genome G+C content by AT/GC substitutions at the codon third position. The G+C content of ubiquitin codon third base also shows a positive linear correlation with the overall G+C content of coding regions of compiled genes, indicating the codon choices among synonymous codons reflect the average codon usage pattern of corresponding species. On the other hand, the monoubiquitin gene, which is different from the polyubiquitin gene in gene organization, gene expression, and function of the encoding protein, shows a different codon usage pattern compared with that of the polyubiquitin gene. From comparisons of the levels of synonymous substitutions among ubiquitin repeats and the homology of the amino acid sequence of the tail of monomeric ubiquitin genes, we propose that the molecular evolution of ubiquitin genes occurred as follows: Plural primitive ubiquitin sequences were dispersed on genome in ancestral eukaryotes. Some of them situated in a particular environment fused with the tail sequence to produce monomeric ubiquitin genes that were maintained across species. After divergence of species, polyubiquitin genes were formed by duplication of the other primitive ubiquitin sequences on different chromosomes. Differences in the environments in which ubiquitin genes are embedded reflect the differences in codon choice and in gene expression pattern between poly- and monomeric ubiquitin genes.     

Proliferation and deterioration of Rickettsia palindromic elements

It has been suggested that Rickettsia Palindromic Elements (RPEs) have evolved as selfish DNA that mediate protein sequence evolution by being targeted to genes that code for RNA and proteins. Here, we have examined the phylogenetic depth of two RPEs that are located close to the genes encoding elongation factors Tu (tuf) and G (fus) in Rickettsia. An exceptional organization of the elongation factor genes was found in all 11 species examined, with complete or partial RPEs identified downstream of the tuf gene (RPE-tuf) in six species and of the fus gene (RPE-fus) in 10 species. A phylogenetic reconstruction shows that both RPE-tuf and RPE-fus have evolved in a manner that is consistent with the expected species divergence. The analysis provides evidence for independent loss of RPE-tuf in several species, possibly mediated by short repetitive sequences flanking the site of excision. The remaining RPE-tuf sequences evolve as neutral sequences in different stages of deterioration. Likewise, highly fragmented remnants of the RPE-fus sequence were identified in two species. This suggests that genome-specific differences in the content of RPEs are the result of recent loss rather than recent proliferation.     

Divergence among genes encoding the elongation factor Tu of Yersinia Species

Elongation factor Tu (EF-Tu), encoded by tuf genes, carries aminoacyl-tRNA to the ribosome during protein synthesis. Duplicated tuf genes (tufA and tufB), which are commonly found in enterobacterial species, usually coevolve via gene conversion and are very similar to one another. However, sequence analysis of tuf genes in our laboratory has revealed highly divergent copies in 72 strains spanning the genus Yersinia (representing 12 Yersinia species). The levels of intragenomic divergence between tufA and tufB sequences ranged from 8.3 to 16.2% for the genus Yersinia, which is significantly greater than the 0.0 to 3.6% divergence observed for other enterobacterial genera. We further explored tuf gene evolution in Yersinia and other Enterobacteriaceae by performing directed sequencing and phylogenetic analyses. Phylogenetic trees constructed using concatenated tufA and tufB sequences revealed a monophyletic genus Yersinia in the family Enterobacteriaceae. Moreover, Yersinia strains form clades within the genus that mostly correlate with their phenotypic and genetic classifications. These genetic analyses revealed an unusual divergence between Yersinia tufA and tufB sequences, a feature unique among sequenced Enterobacteriaceae and indicative of a genus-wide loss of gene conversion. Furthermore, they provided valuable phylogenetic information for possible reclassification and identification of Yersinia species.     

Mycoplasma cavipharyngis and Mycoplasma fastidiosum, the closest relatives to Eperythrozoon spp. and Haemobartonella spp.

The 16S rRNA gene sequences of Mycoplasma cavipharyngis and Mycoplasma fastidiosum have been determined. Phylogenetic analysis showed that these species formed a new cluster within the so-called pneumoniae group of the mollicutes (class Mollicutes). This cluster will be referred to as the M. fastidiosum cluster. Interestingly, the M. fastidiosum cluster formed a sister lineage to the haemotrophic bacteria. Eperythrozoon spp. and Haemobartonella spp. The two latter genera, formerly classified as rickettsias, formed a stable phylogenetic entity in the tree as judged from branch lengths, bootstrap values and sequence signatures. Thus, the members of the M. fastidiosum cluster are the closest known relatives to the haemotrophic bacteria. Our data strongly support that the haemotrophic bacteria should be reclassified to reflect their actual phylogenetic affiliation.     

Processes of genome evolution reflected by base frequency differences among Serratia marcescens genes

The G + C content of silent sites in codons varies greatly among Serratia marcescens genes; the value in any one gene seems to reflect a balance between mutation pressure towards high G + C content and natural selection constraining choice among synonymous codons. Interestingly, non-coding sequences have substantially lower G + C content than silent sites thought to be under little selective constraint.     

Gram-positive bacteria with a high DNA G+C content are characterized by a common insertion within their 23S rRNA genes.

An insertion of about 100 bases within the central part of the 23S rRNA genes was found to be a phylogenetic marker for the bacterial line of descent of Gram-positive bacteria with a high DNA G + C content. The insertion was present in 23S rRNA genes of 64 strains representing the major phylogenetic groups of Gram-positive bacteria with a high DNA G+C content, whereas it was not found in 23S rRNA genes of 55 (eu)bacteria representing Gram-positive bacteria with a low DNA G + C content and all other known (eu)bacterial phyla. The presence of the insertion could be easily demonstrated by comparative gel electrophoretic analysis of in vitro-amplified 23S rDNA fragments, which contained the insertion. The nucleotide sequences of the amplified fragments were determined and sequence similarities of at least 44% were found. The overall similarity values are lower than those of 16S and 23S rRNA sequences of the particular organism. Northern hybridization experiments indicated the presence of the insertion within the mature 23S rRNA of Corynebacterium glutamicum.     

Selection on GGU and CGU Codons in the High Expression Genes in Bacteria

The fourfold degenerate site (FDS) in coding sequences is important for studying the effect of any selection pressure on codon usage bias (CUB) because nucleotide substitution per se is not under any such pressure at the site due to the unaltered amino acid sequence in a protein. We estimated the frequency variation of nucleotides at the FDS across the eight family boxes (FBs) defined as Um(g), the unevenness measure of a gene g. The study was made in 545 species of bacteria. In many bacteria, the Um(g) correlated strongly with Nc′—a measure of the CUB. Analysis of the strongly correlated bacteria revealed that the U-ending codons (GGU, CGU) were preferred to the G-ending codons (GGG, CGG) in Gly and Arg FBs even in the genomes with G+C % higher than 65.0. Further evidence suggested that these codons can be used as a good indicator of selection pressure on CUB in genomes with higher G+C %.     

Nucleotide sequence and characteristics of the gene for L-lactate dehydrogenase of Thermus caldophilus GK24 and the deduced amino-acid sequence of the enzyme

The gene for L-lactate dehydrogenase (LDH) (EC 1.1.1.27) of Thermus caldophilus GK24 was cloned in Escherichia coli using synthetic oligonucleotides as hybridization probes. The nucleotide sequence of the cloned DNA was determined. The primary structure of the LDH was deduced from the nucleotide sequence. The deduced amino acid sequence agreed with the NH2-terminal and COOH-terminal sequences previously reported and the determined amino acid sequences of the peptides obtained from trypsin-digested T. caldophilus LDH. The LDH comprised 310 amino acid residues and its molecular mass was determined to be 32,808. On alignment of the whole amino acid sequences, the T. caldophilus LDH showed about 40% identity with the Bacillus stearothermophilus, Lactobacillus casei and dogfish muscle LDHs. The T. caldophilus LDH gene was expressed with the E. coli lac promoter in E. coli, which resulted in the production of the thermophilic LDH. The gene for the T. caldophilus LDH showed more than 40% identity with those for the human and mouse muscle LDHs on alignment of the whole nucleotide sequences. The G + C content of the coding region for the T. caldophilus LDH was 74.1%, which was higher than that of the chromosomal DNA (67.2%). The G + C contents in the first, second and third positions of the codons used were 77.7%, 48.1% and 95.5% respectively. The high G + C content in the third base caused extremely non-random codon usage in the LDH gene. About half (48.7%) the codons in the LDH gene started with G, and hence there were relatively high contents of Val, Ala, Glu and Gly in the LDH. The contents of Pro, Arg, Ala and Gly, which have high G + C contents in their codons, were also high. Rare codons with U or A as the third base were sometimes used to avoid the TCGA sequence, the recognition site for the restriction endonuclease, TaqI. Two TCGA sequences were found only in the sequence of CTCGAG (XhoI site) in the sequenced region of the T. caldophilus DNA. There were three segments with similar sequences in the two 5' non-coding regions, probably the promoter and ribosome-binding regions, of the genes for the T. caldophilus LDH and the Thermus thermophilus 3-isopropylmalate dehydrogenase.     

Characteristics of Human and Mouse Orthologous Protein-Coding Nucleotide Sequences with Large G+C Content Variations

Characteristics of human and mouse orthologous gene sequences which have large G+C content variations were investigated in this study. The orthologous gene pairs were classified into two groups according to the deviation between human and mouse G+C content at the third codon position (GC3) and were subsequently analyzed. In one group, mouse genes had higher GC3 than the corresponding human genes and in another group, human genes had higher GC3 than mouse. Furthermore, the orthologous pairs were separated based on the deviation between human or mouse GC3 and the G+C content at the third codon position of identical codons (IC3), to examine the effect of increased or decreased G+C content in human or mouse sequences. The nucleotide substitution patterns between human and mouse sequences in the two groups were remarkably distinct, and consistent with the state of G+C-rich or G+C-poor sequences. The effect of increase or decrease of G+C content in human or mouse sequences was not clear in the nucleotide substitution patterns. The chromosomal locations of human and mouse orthologous gene pairs were different between the two groups. The genes located on an identical syntenic segment showed the trend of having similar G+C content. Moreover, the same gene order of some genes on different chromosomes of both species demonstrated the gene rearrangements between human and mouse. Our study indicated that the chromosomal locations and rearrangements are associated with the GC3 variation between human and mouse sequences.Key Words: Human mouse orthologs, G+C content variation, nucleotide substitution, gene location, gene rearrangement.     

Comparison of the conserved region in the dnaA gene from three mollicute species

Abstract Polymerase chain reaction was carried out to amplify the conserved region (789 bp in the case of Mycoplasma capricolum ) of the dnaA gene (1350 bp in the case of M. capricolum ) of 15 representatives of the class Mollicutes using degenerate oligonucleotide primers. The dnaA gene fragments were amplified from M. mycoides subsp. capri, Spiroplasma apis and S. citri . The amino acid sequences deduced from the nucleotide sequences of the amplified fragments showed very low similarities to those of the corresponding regions of four walled bacteria. The values of similarity between any two of the three mollicute species were lower than those between any two of the four walled bacteria.     

The conservation of DNA sequences over very long periods of evolutionary time. Evidence against intergeneric chromosomal transfer as an explanation for the presence of Escherichia coli tuf gene sequences in taxonomically-unrelated prokaryotes

In the present study we tried to determine whether the presence of DNA sequences homologous to the Escherichia coli tuf gene (encodes peptide chain elongation factor Tu) in many taxonomically-unrelated prokaryotes is due to selective pressure for these sequences or due to the transfer of chromosomal material subsequent to the divergence of the genera from their progenitors. We found that the degree of sequence homology to the DNA immediately adjacent to the E. coli tuf A gene is either nonexistent or much less than that found for the tuf gene. Furthermore, the tuf-homologous sequences present in one prokaryote were found to be in large part the same as or a subset of those present in others. That is, various prokaryotes share a common subset of tuf-homologous sequences. These findings suggest that strong selective pressure and not recent intergeneric chromosomal transfer is responsible for the ubiquitous presence of certain tuf-homologous sequences. Because the genetic code is degenerate, DNA sequence need not be conserved to conserve protein sequence. Therefore, if the only function of these sequences is to encode protein, their persistence must mean that in some instances codon sequence is selected for.     

Molecular cloning of HSP70 in Mycoplasma ovipneumoniae and comparison with that of other mycoplasmas

Mycoplasma ovipneumoniae, a bacterial species that specifically affects ovine and goat, is the cause of ovine infectious pleuropneumonia. We cloned, sequenced and analyzed heat shock protein 70 (HSP70) (dnaK) gene of M. ovipneumoniae. The full length open reading frame of the M. ovipneumoniae HSP70 gene consists of 1812 nucleotides, with a G+C content of 34.16%, encoding 604 amino acids. Comparative analysis with the HSP70 sequences of 15 Mycoplasma species revealed 59 to 87% DNA sequence identity, with an amino acid sequence identity range of 58 to 94%. M. ovipneumoniae and M. hyopneumoniae shared the highest DNA and amino acid sequence identity (87 and 94%, respectively). Based on phylogenetic analysis, both the DNA and amino acid identities of M. ovipneumoniae with other mycoplasmal HSP70 were correlated with the degree of relationship between the species. The C-terminus of the HSP70 was cloned into a bacterial expression vector and expressed in Escherichia coli cells. The recombinant C-terminal portion of HSP70 protein strongly reacted with convalescent sera from M. ovipneumoniae-infected sheep, based on an immunoblotting assay. This indicates that HSP70 is immunogenic in a natural M. ovipneumoniae infection and may be a relevant antigen for vaccine development.     

Molecular characterization of the P1-like adhesin gene from Mycoplasma pirum.

A DNA fragment has been isolated from the genome of Mycoplasma pirum by use of a genetic probe derived from the conserved region within the genes for the major adhesins of Mycoplasma genitalium and Mycoplasma pneumoniae. A gene encoding an adhesin-like polypeptide was localized, and sequence analysis indicated a G + C content of only 28%, with A- and T-rich codons being preferentially used. A total of 91% of positions 3 were either A or T. The deduced polypeptide is 1,144 amino acids long (126 kDa) and shows 26% identity with the adhesins of M. genitalium and M. pneumoniae. Other features in common with these two membrane proteins include a similar hydropathic profile and a proline-rich C terminus. Antibodies were prepared by using as an immunogen a peptide derived from the C terminus of the M. pirum adhesin-like polypeptide and were found to recognize on immunoblots a 126-kDa polypeptide from an M. pirum cellular extract. The characterization of the adhesin-like gene is a first step toward a better understanding of the mechanisms allowing this human mycoplasma to attach to host cells.     

Evidence for recent intergeneric transfer of a new tetracycline resistance gene, tet(W), isolated from Butyrivibrio fibrisolvens, and the occurrence of tet(O) in ruminal bacteria

We have previously reported high-frequency transfer of tetracycline resistance between strains of the rumen anaerobic bacterium Butyrivibrio fibrisolvens . Donor strains were postulated to carry two Tc^R genes, one of which is transferred on a novel chromosomal element. It is shown here that coding sequences within the non-transmissible gene in B. fibrisolvens 1.230 are identical to those of the Streptococcus pneumoniae tet (O) gene. This provides the first evidence for genetic exchange between facultatively anaerobic bacteria and rumen obligate anaerobes. In contrast, the product of the transmissible Tc^R gene shares only 68% amino acid sequence identity with the TetO and TetM proteins and represents a new class of ribosome protection tetracycline resistance determinant, designated Tet W. The tet (W) coding region shows a higher DNA G + C content (53%) than other B. fibrisolvens genes or other ribosome protection-type tet genes, suggesting recent acquisition from a high G + C content genome. Tet (W) genes with almost identical sequences are also shown to be present in Tc^R strains of B. fibrisolvens from Australian sheep and in Tc^R strains of two other genera of rumen obligate anaerobes, Selenomonas ruminantium and Mitsuokella multiacidus . This provides compelling evidence for recent intergeneric transfer of resistance genes between ruminal bacteria. Tet (W) is not restricted to ruminal bacteria, as it was also present in a porcine strain of M. multiacidus .