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.     

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.     

Distinct Types of rRNA Operons Exist in the Genome of the Actinomycete Thermomonospora chromogena and Evidence for Horizontal Transfer of an Entire rRNA Operon

We describe here the presence of two distinct types of rRNA operons in the genome of a thermophilic actinomycete Thermomonospora chromogena. The genome of T. chromogena contains six rRNA operons (rrn), of which four complete and two incomplete ones were cloned and sequenced. Comparative analysis revealed that the operon rrnB exhibits high levels of sequence variations to the other five nearly identical ones throughout the entire length of the operon. The coding sequences for the 16S and 23S rRNA genes differ by approximately 6 and 10%, respectively, between the two types of operons. Normal functionality of rrnB is concluded on the basis of the nonrandom distribution of nucleotide substitutions, the presence of compensating nucleotide covariations, the preservation of secondary and tertiary rRNA structures, and the detection of correctly processed rRNAs in the cell. Comparative sequence analysis also revealed a close evolutionary relationship between rrnB operon of T. chromogena and rrnA operon of another thermophilic actinomycete Thermobispora bispora. We propose that T. chromogena acquired rrnB operon from T. bispora or a related organism via horizontal gene transfer.     

Comparative and functional analysis of the rRNA-operons and their tRNA gene complement in different lactic acid bacteria

The complete genome sequences of the lactic acid bacteria (LAB), Lactobacillus plantarum, Lactococcus lactis, and Lactobacillus johnsonii were used to compare location, sequence, organisation, and regulation of the ribosomal RNA (rrn) operons. All rrn operons of the examined LAB diverge from the origin of replication, which is compatible with their efficient expression. All operons show a common organisation of 5'-16S-23S-5S-3' structure, but differ in the number, location and specificity of the tRNA genes. In the 16S-23S intergenic spacer region, two of the five rrn operons of Lb. plantarum and three of the six of Lb. johnsonii contain tRNA-ala and tRNA-ile genes, while L. lactis has a tRNA-ala gene in all six operons. The number of tRNA genes following the 5S rRNA gene ranges up to 14, 16, and 21 for L. lactis, Lb. johnsonii and Lb. plantarum, respectively. The tRNA gene complements are similar to each other and to those of other bacteria. Micro-heterogeneity was found within the rRNA structural genes and spacer regions of each strain. In the rrn operon promoter regions of Lb. plantarum and L. lactis marked differences were found, while the promoter regions of Lb. johnsonii showed a similar tandem promoter structure in all operons. The rrn promoters of L. lactis show either a single or a tandem promoter structure. All promoters of Lb. plantarum contain two or three -10 and -35 regions, of which either zero to two were followed by an UP-element. The Lb. plantarum rrnA, rrnB, and rrnC promoter regions display similarity to the rrn promoter structure of Esherichia coli. Differences in regulation between the five Lb. plantarum promoters were studied using a low copy promoter-probe plasmid. Taking copy number and growth rate into account, a differential expression over time was shown. Although all five Lb. plantarum rrn promoters are significantly different, this study shows that their activity was very similar under the circumstances tested. An active promoter was also identified within the Lb. plantarum rrnC operon preceding a cluster of 17 tRNA genes.     

Physical mapping of stable RNA genes in Bacillus subtilis using polymerase chain reaction amplification from a yeast artificial chromosome library.

A new approach for mapping genes which utilizes yeast artificial chromosome clones carrying parts of the Bacillus subtilis genome and the polymerase chain reaction technique is described. This approach was used to physically map stable RNA genes of B. subtilis. Results from over 400 polymerase chain reactions carried out with the yeast artificial chromosome clone library, using primers specific for the genes of interest and designed from published sequences, were collected. The locations of 10 known rRNA gene regions (rrnO, rrnA, rrnE, rrnD, rrnB, rrnJ-rrnW, and rrnI-rrnH-rrnG) have been determined by this method, and these results correlate with those observed by standard genetic mapping. All rRNA operons, except rrnB, are found between 0 and 90 degrees, while rrnB has been placed in the area of 270 degrees on the chromosome map. Also localized were the tRNA gene clusters associated with the following ribosomal operons: rrnB (21 tRNAs), rrnJ (9 tRNAs), rrnD (16 tRNAs), and rrnO and rrnA (2 internal tRNAs). A previously unmapped four-tRNA gene cluster, trnY, a tRNA gene region that is not associated with a ribosomal operon, was found near the origin of replication. The P-RNA gene, important for processing of tRNAs, was found between map locations 197 and 204 degrees.     

Structure, expression and products of the ribosomal RNA operons of Rhodopseudomonas palustris No. 7

Rhodopseudomonas palustris strains carry one or two ribosomal rRNA operons, and those with duplicated rrn operons grow faster. The two rrn operons in R. palustris No. 7 are virtually identical over a 54,70-bp stretch containing the genes for 16S rRNA, tRNAile, tRNAala, 23S rRNA and 5S rRNA, as well as the intergenic spacers and part of the extragenic spacer. In R. palustris, unlike most bacteria with multiple rrn operons, the putative promoter sequences of the two operons are highly diverged, suggesting possible functional differentiation. By simultaneous primer-extension analysis of both pre-rRNAs, we detected a two-fold higher level of expression from rrnA under photoautotrophic conditions. Alteration of the conditions of growth leads to changes in the relative levels of expression of the two operons. Within the 5,470-bp segment, only two sequence differences are found between the 23S rRNA genes; one is at the center of the 23S rRNA molecule and affects a site of unknown function, and the other is within or immediately adjacent to sequences involved in processing of the 5' 23S rRNA IVS. In vitro processing of 5' IVS-containing 23S rRNA precursors from each operon does not reveal any detectable difference between them. The 5' ends of the mature 16S, 23S, and 5S rRNAs were determined by primer-extension analysis, and the 3' end of 23S rRNA was determined by RNA linker ligation-mediated cDNA cloning. The 5' and 3' ends of the R. palustris 23S rRNA molecule are extensively processed, suggesting that, unlike the situation in the established eubacterial model, these ends cannot basepair.     

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.     

Tn10-mediated inversions fuse uridine phosphorylase (udp) and rRNA genes of Escherichia coli.

Two strains carrying metE::Tn10 insertions (upstream of the udp gene) were used to isolate mutants of Escherichia coli overexpressing udp. These strains differ in their gene order; one contains an inversion between the rrnD and rrnE rRNA operons. Selection was based on the ability of overexpressed Udp to complement thymine auxotrophy. Chromosomal rearrangements that connect the udp gene and promoters of different rrn operons were obtained by this selection. Seven of 14 independent mutants selected in one of the initial strains contained similar inversions of the metE-rrnD segment of the chromosome (about 12% of its length). Another mutant contained traces of a more complicated event, inversion between rrnB and rrnG operons, which was followed by reinversion of the segment between metE and the hybrid rrnG/B operon. Similar inversions (udp-rrn) in a strain already carrying an rrnE-rrnD inversion flip the chromosomal segment between metE and rrnD/E in the opposite direction. In this case, inversions are also accompanied by duplications of the chromosomal region between the rrnA and hybrid udp-rrnD/E operons. PCR amplification with a set of oligonucleotides from the rrn, Tn5, and met genes was used for more detailed mapping. Amplified fragments of the rearranged chromosomes connecting rrnD sequences and insertion elements were sequenced, and inversion endpoints were established.     

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.     

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 analysis of 16S rRNA genes amplified from two ribosomal RNA gene clusters of Bifidobacterium bifidum

Two rRNA gene clusters were detected in the genome of Bifidobacterium bifidum KCTC 3202^T using Southern blot analysis. To analyse the sequences of the 16S rRNA genes from rrnA and rrnB, 16S rDNAs were amplified by PCR using DNA fragments purified from gel slices containing each of the rRNA gene clusters. The amplified 16S rDNAs from rrnA and rrnB were cloned into vectors and three clones of each gene sequenced. The resultant sequences were confirmed by direct sequencing of the 16S rDNAs from rrnA and rrnB. Sequence differences were not found between rrnA and rrnB in 1488 bp of the 16S rRNA genes.