Methylation at nucleotide G745 or G748 in 23S rRNA distinguishes Gram-negative from Gram-positive bacteria

Bacteria tune the function of their ribosomes by methylating specific rRNA nucleotides. Nucleotide G745 in Escherichia coli 23S rRNA is methylated by the methyltransferase enzyme RrmA, whereas in Streptomyces fradiae, the neighbouring nucleotide G748 is methylated by the enzyme TlrB. Both nucleotides line the peptide exit channel of the ribosome at the binding site of macrolide, lincosamide and streptogramin B antibiotics. Resistance to the macrolide tylosin, which is produced by S. fradiae, is conferred by methylation of G748. RrmA and TlrB are homologues (29% identical), and a database search against all presently available sequences revealed a further two dozen homologues from a wide variety of Bacteria. No homologues were found among the Archaea or Eukarya. The bacterial sequences adhere to the species phylogeny and segregate into two groups, in which the Gram-negative sequences align with RrmA and the Gram-positives with TlrB. Consistently, in more than 20 species tested, the distribution of methylation in the Gram-negative rRNAs (methylated at G745) and the Gram-positives (methylated at G748) perfectly matches the bacterial phylogeny. Cloning and expression of representative methyltransferase genes showed that this specificity of methylation is determined solely by the methyltransferase enzyme and is independent of the origin of the rRNA substrate. This is the first case in which the position of an RNA methylation defines a sharp division between the Gram-negative and Gram-positive bacteria. Given the specificities and distribution of these methyltransferases, we propose a change in the nomenclature of RrmA to RlmAI (rRNA large subunit methyltransferase) and of TlrB to RlmAII.     

Teleost introns are characterized by a high A+T content

We previously observed that Antarctic fish genes contain intron sequences of high A+T content (60-70% average A+T) which are in stark contrast with adjacent protein coding-sequences. Here, we report that this disparity in intron/exon base composition is a common feature among teleosts. We analyzed 483 teleost genomic DNA sequences, containing 2583 introns, from 80 teleost genera that populate polar, temperate, or tropical habitats. Eighty-nine percent of teleost introns display an A+T content between 50-84% A+T with a mean of 60% A+T. In contrast, only 37% of teleost exons have an A+T content greater-than 50% with a mean of 48% A+T. A comparison to homologous mammalian genes showed a striking difference; in this case, introns and exons have similar base compositions, averaging 45-47% A+T. This indicates that most teleost genes exhibit a large difference in base composition between their introns and exons. There was no correlation of teleost intron A+T content to intron length or habitat temperature range. Thus, teleost intron sequences tend to show the common feature of being much higher in A+T content then neighboring exons.     

tRNA genes are found between 16S and 23S rRNA genes in Bacillus subtilis.

There are at least nine, and probably ten, ribosomal RNA gene sets in the genome of Bacillus subtilis. Each gene set contains sequences complementary to 16S, 23S and 5S rRNAs. We have determined the nucleotide sequences of two DNA fragments which each contain 165 base pairs of the 16S rRNA gene, 191 base pairs of the 23S rRNA gene, and the spacer region between them. The smaller space region is 164 base pairs in length and the larger one includes an additional 180 base pairs. The extra nucleotides could be transcribed in tRNAIIe and tRNA Ala sequences. Evidence is also presented for the existence of a second spacer region which also contains tRNAIIe and tRNA Ala sequences. No other tRNAs appear to be encoded in the spacer regions between the 16S and 23S rRNA genes. Whereas the nucleotide sequences corresponding to the 16S rRNA, 23S rRNA and the spacer tRNAs are very similar to those of E. coli, the sequences between these structural genes are very different.     

The ESAT-6 gene cluster of Mycobacterium tuberculosis and other high G+C Gram-positive bacteria

Background  

The genome of Mycobacterium tuberculosis H37Rv has five copies of a cluster of genes known as the ESAT-6 loci. These clusters contain members of the CFP-10 (lhp) and ESAT-6 (esat-6) gene families (encoding secreted T-cell antigens that lack detectable secretion signals) as well as genes encoding secreted, cell-wall-associated subtilisin-like serine proteases, putative ABC transporters, ATP-binding proteins and other membrane-associated proteins. These membrane-associated and energy-providing proteins may function to secrete members of the ESAT-6 and CFP-10 protein families, and the proteases may be involved in processing the secreted peptide.     

Transfer RNA genes associated with the 16S and 23S rRNA genes of Euglena chloroplast DNA

Hybridization studies of Euglena chloroplast ¹²⁵I-labeled tRNAs to restriction fragments of Euglena chloroplast DNA have shown that the spacer between the 16S and 23S rRNA genes, in two and possibly all three of the ribosomal DNA units, contains genes for tRNA^Ile and tRNA^Ala, whereas a tRNA gene (for either tRNA^Trp or tRNA^Glu) is located before probably all four 16S rRNA genes present on the chloroplast DNA molecule.     

Characterization of the 23S and 5S rRNA genes of Coxiella burnetii and identification of an intervening sequence within the 23S rRNA gene.

Characterization of the rRNA operon from the obligate intracellular bacterium Coxiella burnetii has determined the order of the rRNA genes to be 16S-23S-5S. A 444-bp intervening sequence (IVS) was identified to interrupt the 23S rRNA gene beginning at position 1176. The IVS is predicted to form a stem-loop structure formed by flanking inverted repeats, and the absence of intact 23S rRNA molecules suggests that the loop is removed. An open reading frame in the IVS has been identified that shows 70% similarity at the amino acid level to IVS open reading frames characterized from four species of Leptospira.     

How high G+C Gram-positive bacteria and in particular bifidobacteria cope with heat stress: protein players and regulators

The Actinobacteridae group of bacteria includes pathogens, plant commensals, endosymbionts as well as inhabitants of the gastrointestinal tract. For various reasons, these microorganisms represent a growing area of interest with respect to genomics, molecular biology and genetics. This review will discuss the current knowledge on the molecular players that allow actinobacteria to contend with heat stress, with an emphasis on bifidobacteria. We describe the principal molecular chaperones involved in heat stress. Temporal expression of heat-shock genes based on functional genomics in members of the Actinobacteridae group is also discussed, as well as the emerging molecular mechanisms controlling the heat-stress response.     

Codon usage and G+C content in Bradyrhizobium japonicum genes are not uniform

To date, the sequences of 45 Bradyrhizobium japonicum genes are known. This provides sufficient information to determine their codon usage and G+C content. Surprisingly, B. japonicum nodulation and NifA-regulated genes were found to have a less biased codon usage and a lower G+C content than genes not belonging to these two groups. Thus, the coding regions of nodulation genes and NifA-regulated genes could hardly be identified in codon preference plots whereas this was not difficult with other genes. The codon frequency table of the highly biased genes was used in a codon preference plot to analyze the RSRj9 sequence which is an insertion sequence (IS)-like element. The plot helped identify a new open reading frame (ORF355) that escaped previous detection because of two sequencing errors. These were now corrected. The deduced gene product of ORF355 in RSRj9 showed extensive similarity to a putative protein encoded by an ORF in the T-DNA of Agrobacterium rhizogenes. The DNA sequences bordering both ORFs showed inverted repeats and potential target site duplications which supported the assumption that they were IS-like elements.     

Specific oligonucleotide probes for in situ detection of a major group of gram-positive bacteria with low DNA G + C content.

Almost one thousand 16S rRNA sequences of Gram-positive bacteria with a low DNA G + C content from public databases were analyzed using the ARB software package. A signature region was identified between positions 354 and 371 (E. coli numbering) for the Bacillus sub-branch of the Gram-positive bacteria with a low DNA G + C content, the former orders Bacillales and Lactobacillales. Three oligonucleotide probes, namely LGC354A, LGC354B, and LGC354C, were designed to target this diagnostic site. Their fluorescent derivatives were suitable for whole cell detection by fluorescence in situ hybridization (FISH). Hybridization conditions were adjusted for differentiation of target and related non-target reference species. When applying FISH to whole bacterial cells in a sample of activated sludge from a communal wastewater treatment plant, members of the Bacillus sub-branch were detected at levels from 0.01% of cells in samples fixed with paraformaldehyde to over 8 percent in the same samples fixed with ethanol and treated with lysozyme. The problems of quantitative in situ analysis of Gram-positive bacteria with a low DNA G + C content in biofilm flocs are discussed and recommendations made. Members of the Bacillus sub-branch were detected in different abundances in activated sludge samples from different wastewater plants.     

Identification of a structural motif of 23S rRNA interacting with 5S rRNA.

To identify RNA motifs interacting with 5S rRNA, a systematic evolution of ligands by exponential enrichment experiment was applied. Some of the resulting RNA aptamers contained a consensus sequence similar to the sequence in the loop region of helix 89 of 23S rRNA. We show that the synthetic helix 89 RNA motif indeed interacted with 5S rRNA and that the region around loop B of 5S rRNA was involved in this interaction. These results suggest the presence of a novel RNA-RNA interaction between 23S rRNA and 5S rRNA which may play an important role in the ribosome function.     

Increased G + C content of DNA stabilizes methyl CpG dinucleotides.

The vertebrate genome is a mosaic of regions differing dramatically in their G + C content. Those regions with a high G + C content contain the expected number of CpG dinucleotides and we propose that following methylation these have been protected from deamination by the increased stability of the surrounding DNA duplex. This argument applies both to the microenvironment of the CpG dinucleotide and to whole gene regions.     

Affinity purification of ribosomes with a lethal G2655C mutation in 23 S rRNA that affects the translocation

A method for preparation of Escherichia coli ribosomes carrying lethal mutations in 23 S rRNA was developed. The method is based on the site-directed incorporation of a streptavidin binding tag into functionally neutral sites of the 23 S rRNA and subsequent affinity chromatography. It was tested with ribosomes mutated at the 23 S rRNA position 2655 (the elongation factor (EF)-G binding site). Ribosomes carrying the lethal G2655C mutation were purified and studied in vitro. It was found in particular that this mutation confers strong inhibition of the translocation process but only moderately affects GTPase activity and binding of EF-G.     

DNA G+C content of the third codon position and codon usage biases of human genes

The human genome, as in other eukaryotes, has a wide heterogeneity in the DNA base composition. The evolutionary basis for this heterogeneity has been unknown. A previous study of the human genome (846 genes analyzed) has shown that, in the major range of the G+C content in the third codon position (0.25-0.75), biases from the Parity Rule 2 (PR2) among the synonymous codons of the four-codon amino acids are similar except in the highest G+C range (Sueoka, N., 1999. Translation-coupled violation of Parity Rule 2 in human genes is not the cause of heterogeneity of the DNA G+C content of third codon position. Gene 238, 53-58.). PR2 is an intra-strand rule where A=T and G=C are expected when there are no biases between the two complementary strands of DNA in mutation and selection rates (substitution rates). In this study, 14,026 human genes were analyzed. In addition, the third codon positions of two-codon amino acids were analyzed. New results show the following: (a) The G+C contents of the third codon position of human genes are scattered in the G+C range of 0.22-0.96 in the third codon position. (b) The PR2 biases are similar in the range of 0.25-0.75, whereas, in the high G+C range (0.75-0.96; 13% of the genes), the PR2-bias fingerprints are different from those of the major range. (c) Unlike the PR2 biases, the G+C contents of the third codon position for both four-codon and two-codon amino acids are all correlated almost perfectly with the G+C content of the third codon position over the total G+C ranges. These results support the notion that the directional mutation pressure, rather than the directional selection pressure, is mainly responsible for the heterogeneity of the G+C content of the third codon position.     

FramePlot: a new implementation of the Frame analysis for predicting protein-coding regions in bacterial DNA with a high G+C content

FramePlot is a web-based tool for predicting protein-coding regions in bacterial DNA with a high G + C content, such as Streptomyces. The graphical output provides for easy distinction of protein-coding regions from non-coding regions. The plot is a clickable map. Clicking on an ORF provides not only the nucleotide sequence but also its deduced amino acid sequence. These sequences can then be compared to the NCBI sequence database over the Internet. The program is freely available for academic purposes at http://www.nih.go.jp/jun/cgi-bin/frameplot.pl.     

Two-dimensional display and whole genome comparison of bacterial pathogen genomes of high G+C DNA content

Three highly mutable loci of the wall-less pathogens Mycoplasma bovis, Mycoplasma pulmonis and Mycoplasma agalactiae undergo high-frequency genomic rearrangements and generate extensive antigenic variation of major surface lipoproteins. Adjacent to each locus, an open reading frame exists as a single chromosomal copy and is predicted to encode a site-specific DNA recombinase exhibiting high homology to the recombinases XerD of Escherichia coli and CodV of Bacillus subtilis. Each of the mycoplasmal proteins are members of the lambda integrase family of tyrosine site-specific recombinases and likely mediates site-specific DNA inversions observed within the adjacent, variable loci.     

DNA probes with different specificities from a cloned 23S rRNA gene of Micrococcus luteus

A 7500 bp PstI restriction fragment of chromosomal DNA from Micrococcus luteus containing a 23S rRNA gene was cloned in vector pHE3 in E. coli RR 28 (the recombinant plasmid was designated pAR1). A recombinant phage (pAR5) hybridizing to all eubacteria tested was constructed by shotgun subcloning of the PstI fragment in phage M13mp8. Further subcloning of the fragments of the 23S rRNA gene in the vectors pTZ18R and pTZ19R using selected restriction sites of the gene enabled us to select cloned fragments of the 23S rRNA gene representing different specificities. Probes specific for Micrococcus luteus-Micrococcus lylae (pAR28), for the Arthrobacter-Micrococcus group (pAR27), for eubacteria (pAR5), and for the detection of eu- and archaebacteria (the so-called universal probe pAR17) were constructed. The specificity of each probe was analysed by dot hybridization to the chromosomal DNAs of representatives of most of the main phyla of eu- and archaebacteria.     

A novel method to calculate the G+C content of genomic DNA sequences.

The base composition of a DNA fragment or genome is usually measured by the proportion of A+T or G+C in the sequence. The G+C content along genomic sequences is usually calculated using an overlapping or non-overlapping sliding window method. The result and accuracy of such an approach depends on the size of the window and the moving distance adopted. In this paper, a novel windowless technique to calculate the G+C content of genomic sequences is proposed. By this method, the G+C content can be calculated at different "resolution". In an extreme case, the G+C content may be computed at a specific point, rather than in a window of finite size. This is particularly useful to analyze the fine variation of base composition along genomic sequences. As the first example, the variation of G+C content along each of 16 yeast chromosomes is analyzed. The G+C-rich regions with length larger than 5 kb sequences are detected and listed in details. It is found that each chromosome consists of several G+C-rich and G+C-poor regions alternatively, i.e., a mosaic structure. Another example is to analyze the G+C content for each of the two chromosomes of the Vibrio cholerae genome. Based on the variations of the G+C content in each chromosome, it is shown that some fragments in the Vibrio cholerae genome may have been transferred from other species. Especially, the position and size of the large integron island on the smaller chromosome was precisely predicted. This method would be a useful tool for analyzing genomic sequences.     

Estimating evolutionary distance from restriction maps of mitochondrial DNA with arbitrary G+C content

Summary We develop a mathematical model for estimating evolutionary distance from restriction enzyme maps, which incorporates non-uniformity of the rate of base substitution into the theory and allows for an arbitrary G+C content at equilibrium. When the G+C content differs significantly from 1/2, the traditional model of base changes can introduce a systematic bias which depends upon the base composition of the restriction site. In addition, the accuracy of estimated evolutionary distance depends heavily upon the choice of restriction enzyme in that the expected number of sites is also affected. Monte Carlo experiments are conducted to check the validity of the present theoretical treatment and from which we draw several cautionary notes on estimation. An application is made to the available data on restriction enzyme maps of human mitochondrial DNA where the G+C content is approximately 1/3.Contribution No. 1372 from the National Institute of Genetics, Mishima, 411 Japan