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.     

Heterogeneity of nucleotide sequences of 16S ribosomal RNA genes from the Desulfotomaculum kuznetsovii type strain]

Two copies of the 16S rRNA gene, rrnA and rrnB, of the type strain 17T of the thermophilic sulfate-reducing bacterium Desulfotomaculum kuznetsovii were cloned and completely sequenced. The comparison of the determined sequences revealed considerable heterogeneity (8.3%) of the two genes, rrnA and rrnB. The main differences were associated with superlong inserts located at the variable 5'- and 3'-terminal regions of the 16S rRNA genes. Comparative analysis that involved analogous genes from the phylogenetically closest representatives of the genus Desulfotomaculum showed that disregard of the heterogeneity of the two gene copies distorts the position of the bacterium studied in the phylogenetic tree.     

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.     

Structure of the promoter region for the rrnB gene in Escherichia coli.

The nucleotide sequence of the promoter region for the rrnB gene of E. coli had been determined by the Maxam-Gilbert technique. The 700 bp long sequence had been compared with the published sequences of four other rRNA promoter regions. The rrnB sequence was found to be homologous with the rrnA promoter sequence till the 370th base upstream from the coding region of mature 16S rRNA. The significance of this homology is discussed and a tentative model is proposed to account for the unusual properties of the rRNA promoters.     

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.     

Precise characterization of rDNA genes by intraspecies and inter-loci comparison of rDNA sequences and biochemical analysis of ribosomal RNA molecules in Agrobacterium tumefaciens

Annotation of rRNA genes has been incomplete in Agrobacterium species although a number of Agrobacterial rDNA fragments have been sequenced. In this study, precise characterization of rRNA operons (rrn) was carried out in two biovar 1 strains, C58 and MAFF301001. Complete DNA sequencing of four rrns in MAFF301001 indicated that each operon codes for 16S, 23S and 5S rRNA as well as three tRNAs, trn(Ile), trn(Ala) and trn(Met). The genes and 16S-23S ITS of a given locus were exactly identical with those in the other three loci, except for a T-base loss in the 23S rRNA gene of rrnA and in the 5S rRNA gene of rrnB. Comparison with the four C58 rDNAs available in the DNA database indicated extensive sequence and size variations in the 23S rRNA gene, suggesting the presence of an intervening sequence (IVS). Biochemical RNA analysis, including Northern hybridization and 5' end mapping, in MAFF301001 revealed 2886-base and 2571-base precursors, two 1.3-kb major fragments, a 150-base fragment and removal of an IVS for 23S rRNA. We confirmed similar biochemical characteristics in the C58 strain. The features of rDNA detected here enable correction of previously reported information about Agrobacterial rRNAs and rRNA genes and should be useful for phylogenetic considerations.     

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.     

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.     

A comparative study of the ribosomal RNA operons of Streptomyces coelicolor A3(2) and sequence analysis of rrnA.

S. coelicolor A3(2) contains six ribosomal RNA operons. Here we describe the cloning of rrnA, rrnC and rrnE, thereby completing the cloning of all operons. Southern hybridisation of genomic DNA with a heterologous probe from the E.coli rrnB 16S rRNA gene showed differences in hybridisation among the six rRNA operon-containing bands. The nucleotide sequence of the 16S rRNA gene and the upstream region of rrnA was determined and compared with the corresponding sequence of rrnD, showing that the 16S rRNA genes are 99% identical. Substantial differences were found, however, in the upstream regions corresponding to the P1 and P2 promoters of rrnD. Southern analysis showed that some of the other rRNA operons of S.coelicolor A3(2) also differed in this part of the upstream region.     

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.     

Characterization of Streptomyces venezuelae ATCC 10595 rRNA gene clusters and cloning of rrnA.

Streptomyces venezuelae ATCC 10595 harbors seven rRNA gene clusters which can be distinguished by BglII digestion. The three rRNA genes present in each set are closely linked with the general structure 16S-23S-5S. We cloned rrnA and sequenced the 16S-23S spacer region and the region downstream of the 5S rRNA gene. No tRNA gene was found in these regions.     

Sequence Heterogeneity of the Ten rRNA Operons in Clostridium perfringens

We have cloned and sequenced rRNA operons of Clostridium perfringens strain 13 and analyzed the sequence structure in view of the phylogenesis. The organism had ten copies of rRNA operons all of that comprised of 16S, 23S and 5S rDNAs except for one operon. The operons clustered around the origin of replication, ranging within one-third of the whole genome sequence as it is arranged in a circle. Seven operons were transcribed in clockwise direction, and the remaining three were transcribed in counter clockwise direction assuming that the gyrA was transcribed in clockwise direction. Two of the counter clockwise operons contained tRNA^Ile genes between the 16S and 23S rDNAs, and the other had a tRNA^Ile genes between the 16S and 23S rDNAs and a tRNA^Asn gene in the place of the 5S rDNA. Microheterogeneity was found within the rRNA structural genes and spacer regions. The length of each 16S, 23S and 5S rDNA were almost identical among the ten operons, however, the intergenic spacer region of 16S-23S and 23S-5S were variable in the length depending on loci of the rRNA operons on the chromosome. Nucleotide sequences of the helix 19, helix 19a, helix 20 and helix 21 of 23S rDNA were divergent and the diversity appeared to be correlated with the loci of the rRNA operons on the chromosome.     

Diversity and structure of the archaeal community in the leachate of a full-scale recirculating landfill as examined by direct 16S rRNA gene sequence retrieval

The diversity and structure of the archaeal community in the effluent leachate from a full-scale recirculating landfill was characterized by direct 16S rRNA gene (16S rDNA) retrieval. Total-community DNA was extracted from the microbial assemblages in the landfill leachate, and archaeal 16S rDNAs were amplified with a universally conserved primer and an Archaea-specific primer. The amplification product was then used to construct a 16S rDNA clone library, and 70 randomly selected archaeal clones in the library were grouped by restriction fragment length polymorphism (RFLP) analysis. Sequencing and phylogenetic analysis of representatives from each unique RFLP type showed that the archaeal library was dominated by methanogen-like rDNAs. Represented in the kingdom of Euryarchaeota were phylotypes highly similar to the methanogenic genera Methanoculleus, Methanosarcina, Methanocorpusculum, Methanospirillum and Methanogenium, where the clone distribution was 48, 11, 3, 1 and 1, respectively. No sequences related to known Methanosaeta spp. were retrieved. Four rDNA clones were not affiliated with the known methanogenic Archaea, but instead, they were clustered with the uncultured archaeal sequences recently recovered from anaerobic habitats. Two chimeric sequences were identified among the clones analyzed.