Identification of characteristic oligonucleotides in the bacterial 16S ribosomal RNA sequence dataset

MOTIVATION: The phylogenetic structure of the bacterial world has been intensively studied by comparing sequences of 16S ribosomal RNA (16S rRNA). This database of sequences is now widely used to design probes for the detection of specific bacteria or groups of bacteria one at a time. The success of such methods reflects the fact that there are local sequence segments that are highly characteristic of particular organisms or groups of organisms. It is not clear, however, the extent to which such signature sequences exist in the 16S rRNA dataset. A better understanding of the numbers and distribution of highly informative oligonucleotide sequences may facilitate the design of hybridization arrays that can characterize the phylogenetic position of an unknown organism or serve as the basis for the development of novel approaches for use in bacterial identification. RESULTS: A computer-based algorithm that characterizes the extent to which any individual oligonucleotide sequence in 16S rRNA is characteristic of any particular bacterial grouping was developed. A measure of signature quality, Q(s), was formulated and subsequently calculated for every individual oligonucleotide sequence in the size range of 5-11 nucleotides and for 15mers with reference to each cluster and subcluster in a 929 organism representative phylogenetic tree. Subsequently, the perfect signature sequences were compared to the full set of 7322 sequences to see how common false positives were. The work completed here establishes beyond any doubt that highly characteristic oligonucleotides exist in the bacterial 16S rRNA sequence dataset in large numbers. Over 16,000 15mers were identified that might be useful as signatures. Signature oligonucleotides are available for over 80% of the nodes in the representative tree.     

A possible role of the 5'' terminal sequence of 16S ribosomal RNA in the recognition of initiation sequences for protein synthesis.

Extensive complementarity is found between the 5' end of 16S ribosomal RNA and protein synthesis initiation sites of bacteriophage RNA. Hybrids can be constructed from base sequences of 16S-RNA and two initiation regions on phage RNA. A model is proposed for the involvement of 16S-RNA in the unfolding of hairpin loops containing the initiation codon AUG.     

Human 28S ribosomal RNA sequence heterogeneity.

DNA sequencing of several cloned human 28S ribosomal RNA gene fragments has revealed sequence heterogeneity (1) but it was not clear whether these are inactive pseudogenes or are active genes that are transcribed and represented in ribosomes. S1 nuclease analysis allowed us to examine the population of ribosomal RNA molecules of a cell, and we found that 28S rRNA is a heterogeneous assortment of molecules in both mono- and polysomal preparations. Sequence variation, although largely concentrated in variable regions of the molecule, apparently also occurs in the conserved regions.     

Identification of probiotic lactobacilli used for animal feeds on the basis of 16S ribosomal RNA gene sequence

The use of probiotics such as Lactobacillus in animal feeds has gained popularity in recent years. In this study the 16S rRNA gene sequence of L. acidophilus in two commercial agents which have been used in animal feeds, LAB‐MOS and Ghenisson 22, was determined. Phylogenetic tree analysis revealed that the two agents, strain MNFLM01 in LAB‐MOS and strain GAL‐2 in Ghenisson 22, belonged to L. rhamnosus (a member of the L. casei group) and L. johnsonii (a member of the L. acidophilus group), respectively. Biochemical tests assigned the two as L. rhamnosus and ambiguously as L. acidophilus. The data suggest that 16S rRNA gene sequence analysis provides more accurate identification of Lactobacillus species than biochemical tests and would allow quality assurance of relevant commercial products. The 16S rRNA gene sequences of strains MNFLM01 and GAL‐2 determined in this study have been submitted to the DDBJ/EMBL/GenBank accession numbers under accession numbers AB288235 and AB295648, respectively.     

Identification and characterization of an intervening sequence within the 23S ribosomal RNA genes of Campylobacter jejuni

Campylobacter jejuni is a significant cause of bacterial enteritis in humans. Three of seven C. jejuni isolates examined were found to contain fragmented 23S rRNA. The occurrence of fragmented 23S rRNA correlated with the presence of an intervening sequence (IVS) within the 23S rRNA genes. The IVS is 157 nucleotides in length and replaces an eight nucleotide sequence in the 23S rRNA genes of C. jejuni isolates that contain intact 23S rRNA. The two ends of the IVS share 31 bases of complementarity that could form a stem-loop structure. Fragmentation of the 23S ribosomal RNA results from the excision of the IVS from the transcribed RNA; the 3’ cleavage site maps within the putative stem-loop formed by the IVS. Southern hybridization analysis revealed that the IVS is not present in the genomes of isolates that contain intact 23S rRNA, suggesting that the IVS is not derived from Campylobacter chromosomal sequences. The C. jejuni IVS is located at a position analogous to that of the IVSs found in both Salmonella and Yersinia spp.     

Identification of a site of psoralen crosslinking in E. coli 16S ribosomal RNA

We have developed a 2-dimensional gel method for identification of RNA sequences crosslinked by the intercalative drug 4'-hydroxymethyl-4,5',8-trimethylpsoralen (HMT). This method is being used to localize such sites in E. coli ribosomal RNA. We report here the identification of a site for HMT crosslinking within positions 434 and 497 of 16S rRNA. We suggest a likely site for HMT intercalation, in which residues U548 and U473 become crosslinked via the drug.     

Identification and characterization of an intervening sequence within the 23S ribosomal RNA genes of Edwardsiella ictaluri

Comparison of the 23S rRNA gene sequences of Edwardsiella tarda and Edw. ictaluri confirmed a close phylogenetic relationship between these two fish pathogen species and a distant relation with the 'core' members of the Enterobacteriaceae family. Analysis of the rrl gene for 23S rRNA in Edw. ictaluri revealed the presence of an intervening sequence (IVS) in helix-45. This new 98bp IVS shared 97% nucleotide identity with Salmonella typhimurium helix-45 IVS. Edw. ictaluri helix-45 IVS was present in all Edw. ictaluri strains analyzed and in at least six rrl operons within each cell. Fragmentation of 23S rRNA due to IVS excision by RNase III was observed by methylene blue staining of ribosomal RNA extracted from Edw. ictaluri isolates. This is the first report of an IVS in the 23S rRNA gene of the genus Edwardsiella.     

Identification of the initiation site of poliovirus polyprotein synthesis.

The complete nucleotide sequence of poliovirus RNA has a long open reading frame capable of encoding the precursor polyprotein NCVP00. The first AUG codon in this reading frame is located 743 nucleotides from the 5' end of the RNA and is preceded by eight AUG codons in all three reading frames. Because all proteins that map at the amino terminus of the polyprotein (P1-1a, VP0, and VP4) are blocked at their amino termini and previous studies of ribosome binding have been inconclusive, direct identification of the initiation site of protein synthesis was difficult. We separated and identified all of the tryptic peptides of capsid protein VP4 and correlated these peptides with the amino acid sequence predicted to follow the AUG codon at nucleotide 743. Our data indicate that VP4 begins with a blocked glycine that is encoded immediately after the AUG codon at nucleotide 743. An S1 nuclease analysis of poliovirus mRNA failed to reveal a splice in the 5' region. We concluded that synthesis of the poliovirus polyprotein is initiated at nucleotide 743, the first AUG codon in the long open reading frame.     

The sequence of the nucleotides at the alpha-sarcin cleavage site in rat 28 S ribosomal ribonucleic acid

The sequence of the 521 nucleotides at the 3' end of a rat 28 S rRNA gene was determined. The region encompasses the site of cleavage of 28 S rRNA by the cytotoxin alpha-sarcin. The toxin hydrolyzes a phosphodiester bond on the 3' side of a guanine residue 393 nucleotides from the 3' end. The alpha-sarcin domain is composed of a purine-rich sequence of 14 highly conserved nucleotides.     

Evidence for maternal imprinting of 45S ribosomal RNA genes in Xenopus hybrids

We discovered that gene clusters of 45S ribosomal RNA in Xenopus hybrid frogs are maternally imprinted, similar to X chromosome inactivation in marsupial females. Paternal expression was partly restored after chemical inhibition of histone deacetylation during larval stages. This provides a new spectacular example of epigenetic silencing and first evidence of genomic imprinting in amphibians.     

The ribosomal binding site for eukaryotic elongation factor EF-2 contains 5 S ribosomal RNA

The possible location of RNA in the ribosomal attachment site for the eukaryotic elongation factor EF-2 was analysed. Stable EF-2 · ribosome complexes formed in the presence of the non-hydrolysable GTP analogue GuoPP[CH₂]P were cross-linked with the short (4 Å between the reactive groups) bifunctional reagent, diepoxybutane. Non-cross-linked EF-2 was removed and the covalent factor-ribosome complex isolated. No interaction between EF-2 and 18 S or 28 S rRNA could be demonstrated. However, density gradient centrifugation of the cross-linked ribosomal complexes showed an increased density (1.25 g/cm³) of the factor, as expected from a covalent complex between EF-2 and a low-molecular-weight RNA species. Treatment of the covalent ribosome-factor complexes with EDTA released approx 50% of the cross-linked EF-2 from the ribosome together with the 5 S rRNA · protein L5 complex. Furthermore, the complex co-migrated with the 5S rRNA · L5 particle in sucrose gradients. Polyacrylamide gel electrophoresis showed that EF-2 was directly linked to 5 S rRNA in the 5 S rRNA · L5 complex, as well as in the complexes isolated by density gradient centrifugation. No traces of 5.8 S rRNA or tRNA could be demonstrated. The data indicate that the ribosomal binding domain for EF-2 contains the 5 S rRNA · protein L5 particle and that EF-2 is located in close proximity to 5 S rRNA within the EF-2 · GuoPP[CH₂]P · ribosome complex.     

Site-directed mutagenesis of the binding site for ribosomal protein S8 within 16S ribosomal RNA from Escherichia coli.

Twelve specific alterations have been introduced into the binding site for ribosomal protein S8 in Escherichia coli 16S rRNA. Appropriate rDNA segments were first cloned into bacteriophage M13 vectors and subjected to bisulfite and oligonucleotide-directed mutagenesis in vitro. Subsequently, the mutagenized sequences were placed within the rrnB operon of plasmid pNO1301 and the mutant plasmids were used to transform E. coli recipients. The growth rates of cells containing the mutant plasmids were determined and compared with that of cells containing the wild-type plasmid. Only those mutations which occurred at highly conserved positions, or were expected to disrupt the secondary structure of the binding site, increased the doubling time appreciably. The most striking changes in growth rate resulted from mutations that altered a small internal loop within the S8 binding site. This structure is phylogenetically conserved in prokaryotic 16S rRNAs and may play a direct role in S8-16S rRNA recognition and interaction.