The organization of two rRNA (rrn) operons of the slow-growing pathogen Mycobacterium celatum provides key insights into mycobacterial evolution

The slow-growing Mycobacterium celatum is known to have two different 16S rRNA gene sequences. This study confirms the presence of two rrn operons and describes their organization. One operon (rrnA) was found to be located downstream from murA and the other (rrnB) was found downstream from tyrS. The promoter regions were sequenced, and also the intergenic transcribed spacer (ITS1 and ITS2) regions separating the 16S rRNA, 23S rRNA and 5S rRNA gene coding regions. Analysis of the RNA fraction revealed that rrnA is regulated by two (P1 and PCL1) promoters and rrnB is regulated by one (P1). These data show that the two rrn operons of M. celatum are organized in the same way as the two rrn operons of classical fast-growing mycobacteria. This information was incorporated into a phylogenetic analysis of the genus based on both 16S rRNA gene sequences and (where possible) the number of rrn operons per genome. The results suggest that the ancestral Mycobacterium possessed two (rrnA and rrnB) operons per genome and that subsequently, on two separate occasions, an operon (rrnB) was lost, leading to two clusters of species having a single operon (rrnA); one cluster includes the classical pathogens and the other includes Mycobacterium abscessus and Mycobacterium chelonae.     

Characterization of an rRNA operon (rrnB) of Mycobacterium fortuitum and other mycobacterial species: implications for the classification of mycobacteria.

Mycobacteria are thought to have either one or two rRNA operons per genome. All mycobacteria investigated to date have an operon, designated rrnA, located downstream from the murA gene. We report that Mycobacteriun fortuitum has a second rrn operon, designated rrnB, which is located downstream from the tyrS gene; tyrS is very close to the 3' end of a gene (3-mag) coding for 3-methylpurine-DNA-glycosylase. The second rrn operon of Mycobacterium smegmatis was shown to have a similar organization, namely, 5' 3-mag-tyrS-rrnB 3'. The rrnB operon of M. fortuitum was found to have a single dedicated promoter. During exponential growth in a rich medium, the rrnB and rrnA operons were the major and minor contributors, respectively, to pre-rRNA synthesis. Genomic DNA was isolated from eight other fast-growing mycobacterial species. Samples were investigated by Southern blot analysis using probes for murA, tyrS, and 16S rRNA sequences. The results revealed that both rrnA and rrnB operons were present in each species. The results form the basis for a proposed new scheme for the classification of mycobacteria. The approach, which is phylogenetic in concept, is based on particular properties of the rrn operons of a cell, namely, the number per genome and a feature of 16S rRNA gene sequences.     

Molecular cloning and comparative sequence analyses of rRNA operons in Streptomyces nodosus ATCC 14899.

The genome of Streptomyces nodosus contains six ribosomal RNA (rRNA) operons. Four of the rRNA operons; rrnB, rrnD, rrnE and rrnF were cloned. We have completely sequenced all four operons, including a region 750 base pairs (bp) upstream of the 16S rRNA gene. The three rRNA genes present in each operon were closely linked in the order 16S-23S-5S. A sequence comparison of the four operons showed more than 99% sequence similarity between the corresponding 16S and 23S rRNA genes, and more than 97% similarity between 5S rRNA genes. The sequence differences observed between 23S rRNA genes appeared to be localized in two specific regions. Substantial sequence differences were found in the region upstream of the 16S rRNA gene as well as in the internal transcribed spacers. No tRNA gene was found in the 16S-23S spacer regions.     

First complete nucleotide sequence and heterologous gene organization of the two rRNA operons in the phytoplasma genome

Phytoplasmas are cell-wallless Gram-positive low G + C bacteria belonging to the Mollicutes that inhabit the cytoplasm of plants and insects. Although phytoplasmas possess two ribosomal RNA (rrn) operons, only one has been fully sequenced. Here, we determined the complete nucleotide sequence of both rrn operons (designated rrnA and rrnB) of onion yellows (OY) phytoplasma. Both operons have rRNA genes organized as 5'-16S-23S-5S-3' with very highly conserved sequences; the 16S, 23S, and 5S rRNA genes are 99.9, 99.8, and 99.1% identical between the two operons. However, the organization of tRNA genes in the upstream region from 16S rRNA gene and in the downstream region from 5S rRNA gene differs markedly. Several promoter candidates were detected upstream from both operons, which suggests that both operons are functional. Interestingly, both have a tRNA(Ile) gene in the 16S-23S spacer region, while the reported rrnB operon of loofah witches' broom phytoplasma does not, indicating heterogenous gene organization of rrnB within phytoplasmas. The phytoplasma tRNA gene organization is similar to that of acholeplasmas, a closely related mollicute, and different from that of mycoplasmas, another mollicute. Moreover, the organization suggests that the rrn operons were derived from that of a related nonmollicute bacterium, Bacillus subtilis. This data should shed light on the evolutionary relationships and phylogeny of the mollicutes.     

Rates and consequences of recombination between rRNA operons

A mutant strain of Escherichia coli was created by inserting a cassette encoding sucrose sensitivity and neomycin resistance (sacB-neo) into the small-subunit rRNA-encoding gene rrs in the rrnB operon. During growth in a complex medium, the cassette was lost from the population, and a complete rrs gene was restored at a rate of 5 x 10(-9) per cell division. Repair of this lesion required flanking regions of DNA that were similar to the six remaining intact rRNA operons and reestablished the full complement of seven rRNA operons. The relative fitness of strains with restored rrnB operons was 1 to 2% higher than that of the mutant strain. The rrnB operon normally contains a spacer region between the 16S and 23S rRNA-encoding genes that is similar in length and tRNA gene content to the spacer in rrnC, -E, and -G. In 2 of the 14 strains in which rrnB was restored, the spacer region had the same length as the spacer region in rrnA, -D, and -H. The requirement for flanking regions of nearly identical DNA and the replication of the spacer region from other rRNA operons during the repair of rrnB suggest that the restoration was accomplished via gene conversion. The rate of gene conversion was 10-fold less than the fixation of point mutations in the same region of the chromosome but was apparently sufficient to homogenize the sequences of rRNA genes in E. coli. These findings are discussed in the context of a conceptual model describing the presence of sequence heterogeneity in coevolving rRNA genes.     

Restriction fragment length polymorphism of rRNA operons for discrimination and intergenic spacer sequences for cataloging of Bacillus subtilis sub-groups

Restriction fragment length polymorphism of rRNA operons (RFLP) and 16S-23S rRNA intergenic region (ISR) sequences of Bacillus subtilis subsp. subtilis, B. subtilis subsp. spizizenii, and B. atrophaeus were compared. ISR sequences of the B. subtilis subspecies were extremely similar (W23 versus 168 rrn H, J, G,W; 96.8%; rrn D, E; 98.4%; rrnB; 97.9%) and, therefore, not useful for their differentiation. However, RFLP of rRNA operons of the B. subtilis subspecies were distinct in terms of numbers and organization within the genome (e.g. the 168 sub-group generally contained 8.3- and 8.0-kb fragments absent in the W23 sub-group). The more distantly related B. atrophaeus was distinct from both B. subtilis subspecies in terms of ISR sequence and rRNA operon number and organization. RFLP of rRNA operons discriminates the two sub-groups of Bacillus subtilis that are indistinguishable by ISR sequence. However, ISR sequence defines the relatedness of B. subtilis to other species (e.g. B. atrophaeus) within the genus Bacillus.     

Mycoplasmas (Mollicutes) have a low number of rRNA genes.

DNA from Mycoplasma, Ureaplasma, Acholeplasma, and Spiroplasma species digested by restriction endonucleases was hybridized with probes consisting of portions of the rrnB rRNA operon of Escherichia coli and the rRNA operon of Mycoplasma capricolum. The results indicate the presence of only one or two sets of rRNA genes in the genome of Mollicutes linked in the procaryotic fashion, 16S-23S-5S.     

Analysis of operons encoding 23S rRNA of Clostridium botulinum type A.

Southern hybridization analysis of Clostridium botulinum type A chromosomal DNA indicated the presence of six copies of the 23S rRNA gene. Fragments of DNA encoding 23S rRNA were amplified by polymerase chain reaction and cloned in Escherichia coli. Three clones examined by restriction enzyme and sequence analysis were found to be derived from different operons. Sequence determination of the entire insert of two clones revealed nine nucleotide changes in the genes coding for 23S rRNA (99.7% sequence identity) between operons encoded on the same chromosome, showing microheterogeneity in the rRNA operons of this organism.     

Divergence and redundancy of 16S rRNA sequences in genomes with multiple rrn operons

The level of sequence heterogeneity among rrn operons within genomes determines the accuracy of diversity estimation by 16S rRNA-based methods. Furthermore, the occurrence of widespread horizontal gene transfer (HGT) between distantly related rrn operons casts doubt on reconstructions of phylogenetic relationships. For this study, patterns of distribution of rrn copy numbers, interoperonic divergence, and redundancy of 16S rRNA sequences were evaluated. Bacterial genomes display up to 15 operons and operon numbers up to 7 are commonly found, but ~40% of the organisms analyzed have either one or two operons. Among the Archaea, a single operon appears to dominate and the highest number of operons is five. About 40% of sequences among 380 operons in 76 bacterial genomes with multiple operons were identical to at least one other 16S rRNA sequence in the same genome, and in 38% of the genomes all 16S rRNAs were invariant. For Archaea, the number of identical operons was only 25%, but only five genomes with 21 operons are currently available. These considerations suggest an upper bound of roughly threefold overestimation of bacterial diversity resulting from cloning and sequencing of 16S rRNA genes from the environment; however, the inclusion of genomes with a single rrn operon may lower this correction factor to ~2.5. Divergence among operons appears to be small overall for both Bacteria and Archaea, with the vast majority of 16S rRNA sequences showing <1% nucleotide differences. Only five genomes with operons with a higher level of nucleotide divergence were detected, and Thermoanaerobacter tengcongensis exhibited the highest level of divergence (11.6%) noted to date. Overall, four of the five extreme cases of operon differences occurred among thermophilic bacteria, suggesting a much higher incidence of HGT in these bacteria than in other groups.     

Construction and Initial Characterization of Escherichia coli Strains with Few or No Intact Chromosomal rRNA Operons

The Escherichia coli genome carries seven rRNA (rrn) operons, each containing three rRNA genes. The presence of multiple operons has been an obstacle to many studies of rRNA because the effect of mutations in one operon is diluted by the six remaining wild-type copies. To create a tool useful for manipulating rRNA, we sequentially inactivated from one to all seven of these operons with deletions spanning the 16S and 23S rRNA genes. In the final strain, carrying no intact rRNA operon on the chromosome, rRNA molecules were expressed from a multicopy plasmid containing a single rRNA operon (prrn). Characterization of these rrn deletion strains revealed that deletion of two operons was required to observe a reduction in the growth rate and rRNA/protein ratio. When the number of deletions was extended from three to six, the decrease in the growth rate was slightly more than the decrease in the rRNA/protein ratio, suggesting that ribosome efficiency was reduced. This reduction was most pronounced in the Delta7 prrn strain, in which the growth rate, unlike the rRNA/protein ratio, was not completely restored to wild-type levels by a cloned rRNA operon. The decreases in growth rate and rRNA/protein ratio were surprisingly moderate in the rrn deletion strains; the presence of even a single operon on the chromosome was able to produce as much as 56% of wild-type levels of rRNA. We discuss possible applications of these strains in rRNA studies.     

Variation in 16S-23S rRNA intergenic spacer regions among Bacillus subtilis 168 isolates

The genome of the Bacillus subtilis 168-type strain contains 10 ribosomal RNA (rRNA) operons. In the intergenic spacer region (ISR) between the 16S and 23S rRNA genes, five rRNA operons, rrnI-H-G and rrnJ-W, lack a trinucleotide signature region. Precise determination of molecular weight (MW), using electrospray mass spectrometry (MS), of the polymerase chain reaction (PCR) products from a segment of the ISR from the 168-type strain and B. subtilis 168-like strain 23071 demonstrated 114 and 111 basepair (bp) PCR products (due to the presence or absence of the insert in the operons) as predicted from sequence. However, PCR of the ISR segment for five other B. subtilis 168 isolates generated only a 114 bp PCR product, suggesting the presence of the trinucleotide signature region in all rRNA operons for these strains. Additional genetic variability between the seven B. subtilis 168 isolates was demonstrated by restriction fragment length polymorphism (RFLP) of the rRNA operons, with three distinct patterns found upon Southern blot analysis. The 168-type strain and three others (23066, 23067, and 23071) exhibited the same Southern pattern. Thus, operon deletion is not responsible for the absence of a 111 bp product on MS analysis for strains 23066 and 23067. Restriction analysis confirmed the presence of the trinucleotide signature region in the ISR of all rRNA operons for five B. subtilis 168 isolates; sequencing of rrnW/H from a representative strain also upheld this finding. These results help provide a better understanding of variations in sequence, operon number and chromosomal organization, both within a genome and among isolates of B. subtilis subgroup 168. It is also hypothesized that the presence of the trinucleotide insert in certain rRNA operons may play a role in rRNA maturation and protein synthesis.     

Loss of the spacer loop sequence from the rrnB operon in the Escherichia coli K-12 subline that bears the relA1 mutation.

A polymorphism affecting the spacer region of the rrnB rRNA operon is described. Strains from a major Escherichia coli K-12 subbranch are missing a 106-nucleotide portion of the rrnB 16S-to-23S spacer, and a 20-nucleotide sequence is found in its place. We have called this mutant operon rrnB2. The rrnB2 spacer was most probably derived from either rrnC or rrnE. This alteration of rrnB may have occurred by a recombinational exchange or by gene conversion. In the genealogy of E. coli K-12 strains, the appearance of rrnB2 is associated with the spontaneous occurrence of the first relaxed mutation, but attempts to show a selective relationship between the two mutational events have had negative results. The sequences of the rrnG and rrnC 16S-to-23S spacers have also been determined and their comparisons to the other rrn operons encoding tRNAGlu2 are presented.     

Sequence Heterogeneity between the Two Genes Encoding 16s Rrna from the Halophilic Archaebacterium Haloarcula Marismortui

The halophilic archaebacterium, Haloarcula marismortui, contains two nonadjacent ribosomal RNA operons, designated rrnA and rrnB, in its genome. The 16S rRNA genes within these operons are 1472 nucleotides in length and differ by nucleotide substitutions at 74 positions. The substitutions are not uniformly distributed but rather are localized within three domains of 16S rRNA; more than two-thirds of the differences occur within the domain bounded by nucleotides 508 and 823. This domain is known to be important for P site binding of aminoacylated tRNA and for 30-50S subunit association. Using S1 nuclease protection, it has been shown that the 16S rRNAs transcribed from both operons are equally represented in the functional 70S ribosome population. Comparison of these two H. marismortui sequences to the 16S gene sequences from related halophilic genera suggests that (i) in diverging genera, mutational differences in 16S gene sequences are not clustered but rather are more generally distributed throughout the length of the 16S sequence, and (ii) the rrnB sequence, particularly within the 508-823 domain, is more different from the out group sequences than is the rrnA sequence. Several possible explanations for the evolutionary origin and maintenance of this sequence heterogeneity within 16S rRNA of H. marismortui are discussed.