Harnessing the intragenomic variability of rRNA operons to improve differentiation of Vibrio species

Although the 16S rRNA gene is frequently used as a phylogenetic marker in analysis of environmental DNA, this marker often fails to distinguish closely related species, including those in the genus Vibrio. Here, we investigate whether inclusion and analysis of 23S rRNA sequence can help overcome the intrinsic weaknesses of 16S rRNA analyses for the differentiation of Vibrio species. We construct a maximum likelihood 16S rRNA gene tree to assess the use of this gene to identify clades of Vibrio species. Within the 16S rRNA tree, we identify the putative informative bases responsible for polyphyly, and demonstrate the association of these positions with tree topology. We demonstrate that concatenation of 16S and 23S rRNA genes increases the number of informative nucleotide positions, thereby overcoming ambiguities in 16S rRNA-based phylogenetic reconstructions. Finally, we experimentally demonstrate that this approach considerably improves the differentiation and identification of Vibrio species in environmental samples.

Bacteria can maintain rRNA operons solely on plasmids for hundreds of millions of years

It is generally assumed that all bacteria must have at least one rRNA operon (rrn operon) on the chromosome, but some strains of the genera Aureimonas and Oecophyllibacter carry their sole rrn operon on a plasmid. However, other related strains and species have chromosomal rrn loci, suggesting that the exclusive presence of rrn operons on a plasmid is rare and unlikely to be stably maintained over long evolutionary periods. Here, we report the results of a systematic search for additional bacteria without chromosomal rrn operons. We find that at least four bacterial clades in the phyla Bacteroidota, Spirochaetota, and Pseudomonadota (Proteobacteria) lost chromosomal rrn operons independently. Remarkably, Persicobacteraceae have apparently maintained this peculiar genome organization for hundreds of millions of years. In our study, all the rrn-carrying plasmids in bacteria lacking chromosomal rrn loci possess replication initiator genes of the Rep_3 family. Furthermore, the lack of chromosomal rrn operons is associated with differences in copy numbers of rrn operons, plasmids, and chromosomal tRNA genes. Thus, our findings indicate that the absence of rrn loci in bacterial chromosomes can be stably maintained over long evolutionary periods.

RNA polymerase redistribution supports growth in E. coli strains with a minimal number of rRNA operons

Bacterial transcription by RNA polymerase (RNAP) is spatially organized. RNAPs transcribing highly expressed genes locate in the nucleoid periphery, and form clusters in rich medium, with several studies linking RNAP clustering and transcription of rRNA (rrn). However, the nature of RNAP clusters and their association with rrn transcription remains unclear. Here we address these questions by using single-molecule tracking to monitor the subcellular distribution of mobile and immobile RNAP in strains with a heavily reduced number of chromosomal rrn operons (Δrrn strains). Strikingly, we find that the fraction of chromosome-associated RNAP (which is mainly engaged in transcription) is robust to deleting five or six of the seven chromosomal rrn operons. Spatial analysis in Δrrn strains showed substantial RNAP redistribution during moderate growth, with clustering increasing at cell endcaps, where the remaining rrn operons reside. These results support a model where RNAPs in Δrrn strains relocate to copies of the remaining rrn operons. In rich medium, Δrrn strains redistribute RNAP to minimize growth defects due to rrn deletions, with very high RNAP densities on rrn genes leading to genomic instability. Our study links RNAP clusters and rrn transcription, and offers insight into how bacteria maintain growth in the presence of only 1-2 rrn operons.

Microbial Identification Using rRNA Operon Region: Database and Tool for Metataxonomics with Long-Read Sequence

Recent development of long-read sequencing platforms has enabled researchers to explore bacterial community structure through analysis of full-length 16S rRNA gene (∼1,500 bp) or 16S-ITS-23S rRNA operon region (∼4,300 bp), resulting in higher taxonomic resolution than short-read sequencing platforms. Despite the potential of long-read sequencing in metagenomics, resources and protocols for this technology are scarce. Here, we describe MIrROR, the database and analysis tool for metataxonomics using the bacterial 16S-ITS-23S rRNA operon region. We collected 16S-ITS-23S rRNA operon sequences extracted from bacterial genomes from NCBI GenBank and performed curation. A total of 97,781 16S-ITS-23S rRNA operon sequences covering 9,485 species from 43,653 genomes were obtained. For user convenience, we provide an analysis tool based on a mapping strategy that can be used for taxonomic profiling with MIrROR database. To benchmark MIrROR, we compared performance against publicly available databases and tool with mock communities and simulated data sets. Our platform showed promising results in terms of the number of species covered and the accuracy of classification. To encourage active 16S-ITS-23S rRNA operon analysis in the field, BLAST function and taxonomic profiling results with 16S-ITS-23S rRNA operon studies, which have been reported as BioProject on NCBI are provided. MIrROR (http://mirror.egnome.co.kr/) will be a useful platform for researchers who want to perform high-resolution metagenome analysis with a cost-effective sequencer such as MinION from Oxford Nanopore Technologies. IMPORTANCE Metabarcoding is a powerful tool to investigate community diversity in an economic and efficient way by amplifying a specific gene marker region. With the advancement of long-read sequencing technologies, the field of metabarcoding has entered a new phase. The technologies have brought a need for development in several areas, including new markers that long-read can cover, database for the markers, tools that reflect long-read characteristics, and compatibility with downstream analysis tools. By constructing MIrROR, we met the need for a database and tools for the 16S-ITS-23S rRNA operon region, which has recently been shown to have sufficient resolution at the species level. Bacterial community analysis using the 16S-ITS-23S rRNA operon region with MIrROR will provide new insights from various research fields.

Structural Basis for Transcript Elongation Control by NusG Family Universal Regulators

NusG/RfaH/Spt5 transcription elongation factors are the only transcription regulators conserved across all life. Bacterial NusG regulates RNA polymerase (RNAP) elongation complexes (ECs) across most genes, enhancing elongation by suppressing RNAP backtracking and coordinating ρ-dependent termination and translation. The NusG paralog RfaH engages the EC only at operon polarity suppressor (ops) sites and suppresses both backtrack and hairpin-stabilized pausing. We used single-particle cryoelectron microscopy (cryo-EM) to determine structures of ECs at ops with NusG or RfaH. Both factors chaperone base-pairing of the upstream duplex DNA to suppress backtracking, explaining stimulation of elongation genome-wide. The RfaH-opsEC structure reveals how RfaH confers operon specificity through specific recognition of an ops hairpin in the single-stranded nontemplate DNA and tighter binding to the EC to exclude NusG. Tight EC binding by RfaH sterically blocks the swiveled RNAP conformation necessary for hairpin-stabilized pausing. The universal conservation of NusG/RfaH/Spt5 suggests that the molecular mechanisms uncovered here are widespread.

Differential Regulation of rRNA and tRNA Transcription from the rRNA-tRNA Composite Operon in Escherichia coli

Escherichia coli contains seven rRNA operons, each consisting of the genes for three rRNAs (16S, 23S and 5S rRNA in this order) and one or two tRNA genes in the spacer between 16S and 23S rRNA genes and one or two tRNA genes in the 3’ proximal region. All of these rRNA and tRNA genes are transcribed from two promoters, P1 and P2, into single large precursors that are afterward processed to individual rRNAs and tRNAs by a set of RNases. In the course of Genomic SELEX screening of promoters recognized by RNA polymerase (RNAP) holoenzyme containing RpoD sigma, a strong binding site was identified within 16S rRNA gene in each of all seven rRNA operons. The binding in vitro of RNAP RpoD holoenzyme to an internal promoter, referred to the promoter of riRNA (an internal RNA of the rRNA operon), within each 16S rRNA gene was confirmed by gel shift assay and AFM observation. Using this riRNA promoter within the rrnD operon as a representative, transcription in vitro was detected with use of the purified RpoD holoenzyme, confirming the presence of a constitutive promoter in this region. LacZ reporter assay indicated that this riRNA promoter is functional in vivo. The location of riRNA promoter in vivo as identified using a set of reporter plasmids agrees well with that identified in vitro. Based on transcription profile in vitro and Northern blot analysis in vivo, the majority of transcript initiated from this riRNA promoter was estimated to terminate near the beginning of 23S rRNA gene, indicating that riRNA leads to produce the spacer-coded tRNA. Under starved conditions, transcription of the rRNA operon is markedly repressed to reduce the intracellular level of ribosomes, but the levels of both riRNA and its processed tRNA^Glu stayed unaffected, implying that riRNA plays a role in the continued steady-state synthesis of tRNAs from the spacers of rRNA operons. We then propose that the tRNA genes organized within the spacers of rRNA-tRNA composite operons are expressed independent of rRNA synthesis under specific conditions where further synthesis of ribosomes is not needed.

Open complex scrunching before nucleotide addition accounts for the unusual transcription start site of E. coli ribosomal RNA promoters

Most Escherichia coli promoters initiate transcription with a purine 7 or 8 nt downstream from the -10 hexamer, but some promoters, including the ribosomal RNA promoter rrnB P1, start 9 nt from the -10 element. We identified promoter and RNA polymerase determinants of this noncanonical rrnB P1 start site using biochemical and genetic approaches including mutational analysis of the promoter, Fe(2+) cleavage assays to monitor template strand positions near the active-site, and Bpa cross-linking to map the path of open complex DNA at amino acid and nucleotide resolution. We find that mutations in several promoter regions affect transcription start site (TSS) selection. In particular, we show that the absence of strong interactions between the discriminator region and σ region 1.2 and between the extended -10 element and σ region 3.0, identified previously as a determinant of proper regulation of rRNA promoters, is also required for the unusual TSS. We find that the DNA in the single-stranded transcription bubble of the rrnB P1 promoter complex expands and is "scrunched" into the active site channel of RNA polymerase, similar to the situation in initial transcribing complexes. However, in the rrnB P1 open complex, scrunching occurs before RNA synthesis begins. We find that the scrunched open complex exhibits reduced abortive product synthesis, suggesting that scrunching and unusual TSS selection contribute to the extraordinary transcriptional activity of rRNA promoters by increasing promoter escape, helping to offset the reduction in promoter activity that would result from the weak interactions with σ.

Ribosomal multi-operon diversity: an original perspective on the genus Aeromonas

16S rRNA gene (rrs) is considered of low taxonomic interest in the genus Aeromonas. Here, 195 Aeromonas strains belonging to populations structured by multilocus phylogeny were studied using an original approach that considered Ribosomal Multi-Operon Diversity. This approach associated pulsed-field gel electrophoresis (PFGE) to assess rrn operon number and distribution across the chromosome and PCR-temporal temperature gel electrophoresis (TTGE) to assess rrs V3 region heterogeneity. Aeromonads harbored 8 to 11 rrn operons, 10 operons being observed in more than 92% of the strains. Intraspecific variability was low or nul except for A. salmonicida and A. aquariorum suggesting that large chromosomic rearrangements might occur in these two species while being extremely rarely encountered in the evolution of other taxa. rrn operon number at 8 as well as PFGE patterns were shown valuable for taxonomic purpose allowing resolution of species complexes. PCR-TTGE revealed a high rate of strains (41.5%) displaying intragenomic rrs heterogeneity. Strains isolated from human samples more frequently displayed intragenomic heterogeneity than strains recovered from non-human and environmental specimens. Intraspecific variability ranged from 0 to 76.5% of the strains. The observation of species-specific TTGE bands, the recovery of identical V3 regions in different species and the variability of intragenomic heterogeneity (1-13 divergent nucleotides) supported the occurrence of mutations and horizontal transfer in aeromonad rrs evolution. Altogether, the presence of a high number of rrn operon, the high proportion of strains harboring divergent rrs V3 region and the previously demonstrated high level of genetic diversity argued in favor of highly adaptative capabilities of aeromonads. Outstanding features observed for A. caviae supported the ongoing process of adaptation to a specialized niche represented by the gut, previously hypothesized. 16S rRNA gene is an informative marker in the genus Aeromonas for both evolutionary and polyphasic taxonomic studies provided that multi-operon fingerprinting approaches are used.

Tidying up international nucleotide sequence databases: ecological, geographical and sequence quality annotation of its sequences of mycorrhizal fungi

Sequence analysis of the ribosomal RNA operon, particularly the internal transcribed spacer (ITS) region, provides a powerful tool for identification of mycorrhizal fungi. The sequence data deposited in the International Nucleotide Sequence Databases (INSD) are, however, unfiltered for quality and are often poorly annotated with metadata. To detect chimeric and low-quality sequences and assign the ectomycorrhizal fungi to phylogenetic lineages, fungal ITS sequences were downloaded from INSD, aligned within family-level groups, and examined through phylogenetic analyses and BLAST searches. By combining the fungal sequence database UNITE and the annotation and search tool PlutoF, we also added metadata from the literature to these accessions. Altogether 35,632 sequences belonged to mycorrhizal fungi or originated from ericoid and orchid mycorrhizal roots. Of these sequences, 677 were considered chimeric and 2,174 of low read quality. Information detailing country of collection, geographical coordinates, interacting taxon and isolation source were supplemented to cover 78.0%, 33.0%, 41.7% and 96.4% of the sequences, respectively. These annotated sequences are publicly available via UNITE (http://unite.ut.ee/) for downstream biogeographic, ecological and taxonomic analyses. In European Nucleotide Archive (ENA; http://www.ebi.ac.uk/ena/), the annotated sequences have a special link-out to UNITE. We intend to expand the data annotation to additional genes and all taxonomic groups and functional guilds of fungi.

A coarse-grained biophysical model of E. coli and its application to perturbation of the rRNA operon copy number

We propose a biophysical model of Escherichia coli that predicts growth rate and an effective cellular composition from an effective, coarse-grained representation of its genome. We assume that E. coli is in a state of balanced exponential steady-state growth, growing in a temporally and spatially constant environment, rich in resources. We apply this model to a series of past measurements, where the growth rate and rRNA-to-protein ratio have been measured for seven E. coli strains with an rRNA operon copy number ranging from one to seven (the wild-type copy number). These experiments show that growth rate markedly decreases for strains with fewer than six copies. Using the model, we were able to reproduce these measurements. We show that the model that best fits these data suggests that the volume fraction of macromolecules inside E. coli is not fixed when the rRNA operon copy number is varied. Moreover, the model predicts that increasing the copy number beyond seven results in a cytoplasm densely packed with ribosomes and proteins. Assuming that under such overcrowded conditions prolonged diffusion times tend to weaken binding affinities, the model predicts that growth rate will not increase substantially beyond the wild-type growth rate, as indicated by other experiments. Our model therefore suggests that changing the rRNA operon copy number of wild-type E. coli cells growing in a constant rich environment does not substantially increase their growth rate. Other observations regarding strains with an altered rRNA operon copy number, such as nucleoid compaction and the rRNA operon feedback response, appear to be qualitatively consistent with this model. In addition, we discuss possible design principles suggested by the model and propose further experiments to test its validity.

A bacterium like E. coli can be thought of as a self-replicating factory, where inventory synthesis, degradation, and management is concerted according to a well-defined set of rules encoded in the organism's genome. Since the organism's survival depends on this set of rules, these rules were most likely optimized by evolution. Therefore, by writing down these rules, what could one learn about Escherichia coli? We examined E. coli growing in the simplest imaginable environment, one constant in space and time and rich in resources, and attempted to identify the rules that relate the genome to the cell composition and self-replication time. With more than 4,400 genes, a full-scale model would be prohibitively complicated, and therefore we “coarse-grained” E. coli by lumping together genes and proteins of similar function. We used this model to examine measurements of strains with reduced copy number of ribosomal-RNA genes, and also to show that increasing this copy number overcrowds the cell with ribosomes and proteins. As a result, there appears to be an optimum copy number with respect to the wild-type genome, in agreement with observation. We hope that this model will improve and further challenge our understanding of bacterial physiology, also in more complicated environments.

Structural organization of the transfer RNA operon I of Vibrio cholerae: differences between classical and El Tor strains

Nine major transfer RNA (tRNA) gene clusters were analysed in various Vibrio cholerae strains. Of these, only the tRNA operon I was found to differ significantly in V. cholerae classical (sixth pandemic) and El Tor (seventh pandemic) strains. Amongst the sixteen tRNA genes contained in this operon, genes for tRNA Gln3 (CAA) and tRNA Leu6 (CUA) were absent in classical strains as compared to El Tor strains. The observation strongly supported the view that the above two pandemic strains constitute two different clones.

Phylogenetic identity, growth-response time and rRNA operon copy number of soil bacteria indicate different stages of community succession

Bacterial community succession corresponds to changes in the phylogenetic identity, growth-response time and rRNA operon (rrn) copy number of culturable populations. To test this hypothesis, we compared the bacterial fractions culturable from the oxic zone of flooded, unplanted paddy soil microcosms after 1-day (early succession) and 70-day (late succession) incubation periods. The proportion of bacteria that was cultivable on solid media corresponded for early and late succession to 37-40% and 31-35% of total DAPI cell counts, which were 7.40 (+/-0.36) x 10(8) and 5.54 (+/-0.28) x 10(8) cells per gram of dry soil, respectively. In colony-forming curve analysis, late successional bacteria showed a significant delay in their growth response compared with those from early succession. A total of 59 early successional isolates grouped into 16 species-level clusters (SLC) plus three Bacilli-like SLC, while 66 late successional isolates formed 25 SLC plus five Bacilli-like SLC. Except Bacilli-like spp., isolates from early succession always belonged to different SLC than those from late succession. Betaproteobacteria and Gammaproteobacteria were typical of the early stage, while Alphaproteobacteria and Actinobacteria prevailed in late succession. Considering all SLC except those assigned to Bacilli, growth-response time and rrn copy number were significantly correlated with successional stage. Isolates of most early successional SLC (14 of 16) formed visible colonies within 1 (11 SLC) or 2 days (three SLC) and contained >or= 4 rrn copies. In contrast, isolates of late successional SLC (23 of 25) formed visible colonies within 2 days (four SLC) or, in most cases, only within 3-15 days (19 SLC) and contained <or= 2 rrn copies. Regardless of whether isolated from early or late succession, Bacilli-like isolates always showed a colony-forming time of 2 days and had 9-11 rrn copies.

Heterogeneous rRNAs are differentially expressed during the morphological development ofStreptomyces coelicolor

It is generally assumed that all mature rRNA molecules assembled into ribosomes within a single cell are identical. However, sequence analysis of Streptomyces coelicolor genome revealed that it harbors six copies of divergent rRNA operons that may express and constitute three and five different kinds of small subunit (SSU) and large subunit (LSU) rRNA molecules, respectively, in a single cell. Phylogenetic analyses of the LSU rRNA genes and the internal transcribed spacer between SSU and LSU genes indicated that the LSU gene of rrnA and rrnE operons might be the result of interspecies recombination between rRNA genes in closely related streptomycetes. Profiling of rRNA species using primer extension analysis showed that heterogeneous rRNA transcripts are expressed and assembled into ribosomes in the cell. As the cells developed from germination to sporulation, the relative amount of LSU rRNA molecules derived from three rRNA operons (rrnA, D, and E) gradually decreased from approximately 85% to approximately 60%, whereas the distribution of LSU rRNA molecules from two other operons (rrnB and F) and rrnC operon gradually increased from approximately 10% to approximately 20% of the total LSU rRNA. These findings indicate that heterogeneous rRNA molecules are differentially expressed during the life cycle of this developmentally complex microorganism.

Intragenomic variation and evolution of the internal transcribed spacer of the rRNA operon in bacteria

Variation in the internal transcribed spacer (ITS) of the rRNA (rrn) operon is increasingly used to infer population-level diversity in bacterial communities. However, intragenomic ITS variation may skew diversity estimates that do not correct for multiple rrn operons within a genome. This study characterizes variation in ITS length, tRNA composition, and intragenomic nucleotide divergence across 155 Bacteria genomes. On average, these genomes encode 4.8 rrn operons (range: 2-15) and contain 2.4 unique ITS length variants (range: 1-12) and 2.8 unique sequence variants (range: 1-12). ITS variation stems primarily from differences in tRNA gene composition, with ITS regions containing tRNA-Ala + tRNA-Ile (48% of sequences), tRNA-Ala or tRNA-Ile (10%), tRNA-Glu (11%), other tRNAs (3%), or no tRNA genes (27%). Intragenomic divergence among paralogous ITS sequences grouped by tRNA composition ranges from 0% to 12.11% (mean: 0.94%). Low divergence values indicate extensive homogenization among ITS copies. In 78% of alignments, divergence is <1%, with 54% showing zero variation and 81% containing at least two identical sequences. ITS homogenization occurs over relatively long sequence tracts, frequently spanning the entire ITS, and is largely independent of the distance (basepairs) between operons. This study underscores the potential contribution of interoperon ITS variation to bacterial microdiversity studies, as well as unequivocally demonstrates the pervasiveness of concerted evolution in the rrn gene family.

Variation in the 16S-23S rRNA intergenic spacer regions in Vibrio parahaemolyticus strains are due to indels nearby their tRNAGlu

Vibrio parahaemolyticus contains 11 rRNA operons each including one of six 16S-23S rRNA gene intergenic spacer classes differing in size and nucleotide sequence. Some of the spacer classes may differ between isolates. We observed that the differences in the spacers between isolates are generally in two spacer classes present in single copies in the genome, one class containing tRNA(Ala) and tRNA(Glu) and the other tRNA(Glu) exclusively. Moreover, these differences are due to indels located nearby their tRNA(Glu). Comparison of the nucleotide sequence between spacer classes suggests that intragenomic nonreciprocal recombination causes the size variations observed in the spacer regions of V. parahaemolyticus strains.

The role of LRP and H-NS in transcription regulation: involvement of synergism, allostery and macromolecular crowding

LRP has recently been shown to interact with the regulatory regions of bacterial ribosomal RNA promoters. Here we study details of the LRP-rDNA interaction by gel retardation and high-resolution footprinting techniques. We show that a second regulator for rRNA transcription, H-NS, facilitates the formation of a higher-order LRP-nucleoprotein complex, probably acting transiently as a DNA chaperone. The macromolecular crowding substance ectoine stabilizes the formation of this dynamic complex, while the amino acid leucine, as a metabolic effector, has the opposite effect. DNase I and hydroxyl radical footprint experiments with LRP-DNA complexes reveal a periodic change of the target DNA structure, which implies extensive DNA wrapping reaching into the promoter core region. We show furthermore that LRP binding is able to constrain supercoils, providing a link between DNA topology and regulation. The results support the conclusion that the bacterial DNA-binding protein LRP, assisted by H-NS, forms a repressive nucleoprotein structure involved in regulation of rRNA transcription. The formation of this regulatory structure appears to be directly affected by environmental changes.

The mosaic nature of intergenic 16S-23S rRNA spacer regions suggests rRNA operon copy number variation in Clostridium difficile strains

Clostridium difficile is a major spore-forming environmental pathogen that causes serious health problems in patients undergoing antibiotic therapy. Consequently, reliable and sensitive methods for typing individual strains are required for epidemiological and environmental studies. Ribotyping is generally considered the best method, but it fails to account for sequence diversity which might exist in intergenic 16S-23S rRNA spacer regions (ISRs) within and among strains of this organism. Therefore, this study was undertaken to compare the sequence of each individual ISR in five strains of C. difficile to explore the extent of this diversity and see whether such information might provide the basis for more sensitive and discriminatory strain typing methods. After targeted PCR amplification, cloning, and sequencing, the diversity of the ISRs was used as a measure of rRNA operon copy number. In C. difficile strains 630, ATCC 43593, A, and B, 11, 11, 7, and 8 ISR length variants, respectively, were found (containing different combinations of sequence groups [i to xiii]), suggesting 11, 11, 7, and 8 rrn copies in the respective strains. Many ISRs of the same length differed markedly in their sequences, and some of these were restricted in occurrence to a single strain. Most of these ISRs did not contain any tRNA genes, and only single copies of the tRNA(Ala) gene were found in those that did. The presence of ISR sequence groups (i to xiii) varied between strains, with some found in one, two, three, four, or all five strains. We conclude that the intergenic 16S-23S rRNA spacer regions showed a high degree of diversity, not only among the rrn operons in different strains and different rrn copies in a single strain but also among ISRs of the same length. It appears that C. difficile ISRs vary more at the inter- and intragenic levels than those of other species as determined by empirical comparison of sequences. The precise characterization of these sequences has demonstrated a high level of mosaic sequence block rearrangements that are present or absent in multiple strain-variable rrn copies within and between five different strains of C. difficile.

A cyanobacterial strain with all chromosomal rRNA operons inactivated: a single nucleotide mutation of 23S rRNA confers temperature-sensitive phenotypes

The presence of a multicopy chromosome, with each copy containing two rRNA operons (rrnA and rrnB), has been an obstacle to analysing mutated rRNA in Synechococcus PCC 7942. To create a system for expressing homogeneous mutated rRNA, the chromosomal rrn operons were sequentially inactivated and a final strain was successfully obtained with all the chromosomal rrn operons inactivated but carrying a replaceable multicopy plasmid containing a single rrn operon. The lag time required for growth response on dark/light shift of mutant strains with chromosomal rrnA or rrnB inactivated was increased 50 % over that of the wild-type strain; however, the presence of the plasmid-borne rrn operon restored the lag time. The doubling time of mutant strains carrying only a functional rrnB operon, but not strains carrying only a functional rrnA operon, was significantly longer than that of the wild-type strain. A strain in which essentially all the cellular 23S rRNA contained the mutation C2588A was temperature sensitive at 16 degrees C and 45 degrees C. Position C2588 is equivalent to C2611 of the peptidyltransferase centre in domain V of Escherichia coli 23S rRNA.

Active transcription of rRNA operons is a driving force for the distribution of RNA polymerase in bacteria: effect of extrachromosomal copies of rrnB on the in vivo localization of RNA polymerase

In contrast to eukaryotes, bacteria such as Escherichia coli contain only one form of RNA polymerase (RNAP), which is responsible for all cellular transcription. Using an RNAP-green fluorescent protein fusion protein, we showed previously that E. coli RNAP is partitioned exclusively in the nucleoid and that stable RNA synthesis, particularly rRNA transcription, is critical for concentrating a significant fraction of RNAP in transcription foci during exponential growth. The extent of focus formation varies under different physiological conditions, supporting the proposition that RNAP redistribution is an important element for global gene regulation. Here we show that extra, plasmid-borne copies of an rRNA operon recruit RNAP from the nucleoid into the cytoplasmic space and that this is accompanied by a reduction in the growth rate. Transcription of an intact rRNA operon is not necessary, although a minimal transcript length is required for this phenotype. Replacement of the ribosomal promoters with another strong promoter, Ptac, abolished the effect. These results demonstrate that active synthesis from rRNA promoters is a major driving force for the distribution of RNAP in bacteria. The implications of our results for the regulation of rRNA synthesis and cell growth are discussed.

Evolution of paralogous genes: Reconstruction of genome rearrangements through comparison of multiple genomes within Staphylococcus aureus

Analysis of evolution of paralogous genes in a genome is central to our understanding of genome evolution. Comparison of closely related bacterial genomes, which has provided clues as to how genome sequences evolve under natural conditions, would help in such an analysis. With species Staphylococcus aureus, whole-genome sequences have been decoded for seven strains. We compared their DNA sequences to detect large genome polymorphisms and to deduce mechanisms of genome rearrangements that have formed each of them. We first compared strains N315 and Mu50, which make one of the most closely related strain pairs, at the single-nucleotide resolution to catalogue all the middle-sized (more than 10 bp) to large genome polymorphisms such as indels and substitutions. These polymorphisms include two paralogous gene sets, one in a tandem paralogue gene cluster for toxins in a genomic island and the other in a ribosomal RNA operon. We also focused on two other tandem paralogue gene clusters and type I restriction-modification (RM) genes on the genomic islands. Then we reconstructed rearrangement events responsible for these polymorphisms, in the paralogous genes and the others, with reference to the other five genomes. For the tandem paralogue gene clusters, we were able to infer sequences for homologous recombination generating the change in the repeat number. These sequences were conserved among the repeated paralogous units likely because of their functional importance. The sequence specificity (S) subunit of type I RM systems showed recombination, likely at the homology of a conserved region, between the two variable regions for sequence specificity. We also noticed novel alleles in the ribosomal RNA operons and suggested a role for illegitimate recombination in their formation. These results revealed importance of recombination involving long conserved sequence in the evolution of paralogous genes in the genome.

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.

Combining two genomes in one cell: stable cloning of the Synechocystis PCC6803 genome in the Bacillus subtilis 168 genome

Cloning the whole 3.5-megabase (Mb) genome of the photosynthetic bacterium Synechocystis PCC6803 into the 4.2-Mb genome of the mesophilic bacterium Bacillus subtilis 168 resulted in a 7.7-Mb composite genome. We succeeded in such unprecedented large-size cloning by progressively assembling and editing contiguous DNA regions that cover the entire Synechocystis genome. The strain containing the two sets of genome grew only in the B. subtilis culture medium where all of the cloning procedures were carried out. The high structural stability of the cloned Synechocystis genome was closely associated with the symmetry of the bacterial genome structure of the DNA replication origin (oriC) and its termination (terC) and the exclusivity of Synechocystis ribosomal RNA operon genes (rrnA and rrnB). Given the significant diversity in genome structure observed upon horizontal DNA transfer in nature, our stable laboratory-generated composite genome raised fundamental questions concerning two complete genomes in one cell. Our megasize DNA cloning method, designated megacloning, may be generally applicable to other genomes or genome loci of free-living organisms.

Diversity of genome structure in Salmonella enterica serovar Typhi populations

The genomes of most strains of Salmonella and Escherichia coli are highly conserved. In contrast, all 136 wild-type strains of Salmonella enterica serovar Typhi analyzed by partial digestion with I-CeuI (an endonuclease which cuts within the rrn operons) and pulsed-field gel electrophoresis and by PCR have rearrangements due to homologous recombination between the rrn operons leading to inversions and translocations. Recombination between rrn operons in culture is known to be equally frequent in S. enterica serovar Typhi and S. enterica serovar Typhimurium; thus, the recombinants in S. enterica serovar Typhi, but not those in S. enterica serovar Typhimurium, are able to survive in nature. However, even in S. enterica serovar Typhi the need for genome balance and the need for gene dosage impose limits on rearrangements. Of 100 strains of genome types 1 to 6, 72 were only 25.5 kb off genome balance (the relative lengths of the replichores during bidirectional replication from oriC to the termination of replication [Ter]), while 28 strains were less balanced (41 kb off balance), indicating that the survival of the best-balanced strains was greater. In addition, the need for appropriate gene dosage apparently selected against rearrangements which moved genes from their accustomed distance from oriC. Although rearrangements involving the seven rrn operons are very common in S. enterica serovar Typhi, other duplicated regions, such as the 25 IS200 elements, are very rarely involved in rearrangements. Large deletions and insertions in the genome are uncommon, except for deletions of Salmonella pathogenicity island 7 (usually 134 kb) from fragment I-CeuI-G and 40-kb insertions, possibly a prophage, in fragment I-CeuI-E. The phage types were determined, and the origins of the phage types appeared to be independent of the origins of the genome types.

[Adaptive reactions of mycoplasmas in vitro: "viable but unculturable forms" and nanocells of Acholeplasma laidlawii]

The adaptation of Acholeplasma laidlawii to conditions unfavorable for growth has been found to be accompanied by cell transformation into special morphological structures known as ultramicroforms (nanocells). The ratio of the cells of the two morphological types in the population depended on the growth conditions. Nanocells retained viability for a long time under conditions unfavorable for growth and showed resistance to stressors. Reduction in the cell size occurred due to unequal division, which involved the loss of cytoplasmic material. A. laidlawii ultramicroforms (nanocells) were able to restore proliferative activity and to revert to their initial vegetative form; they measured less than 0.2 microm and are the smallest cells known at present. Nanocells formed in vitro under exposure to abiogenic stressors may correspond to the A. laidlawii minibodies observed in infected plants upon exposure to biogenic stressors. The transformation of A. laidlawii cells into ultramicroforms was accompanied by condensation of the nucleoid, a change in the polypeptide spectrum, and a change in the availability of rRNA operons for in vitro amplification. All these changes are indicative of reorganization of the genetic and metabolic systems of mycoplasmas.

Construction of a combined physical and genetic map of the chromosome of Lactobacillus acidophilus ATCC 4356 and characterization of the rRNA operons

The combination of PFGE and hybridization approaches was used to study the genome of Lactobacillus acidophilus neotype strain ATCC 4356. PFGE analysis of chromosomal DNA after digestion with each of the rare-cutting restriction enzymes I-CeuI, NotI, CspI, SmaI, ApaI and SgrAI allowed the size of the circular chromosome of L. acidophilus to be estimated at 2.061 Mbp. The physical map contained 86 restriction sites for the six enzymes employed, with intervals between the sites varying from 1 to 88 kbp (approximately 0.05-4.3 % of the chromosome). Based on the physical map, a genetic map was constructed via Southern blot analyses of L. acidophilus DNA using specific gene probes. A total of 73 probes representing key genes, including 12 rRNA (rrn) genes, were positioned on the latter map. Mapping analysis also indicated the presence of four rrn operons (rrnA-D) on the chromosome, each containing a single copy of each of the three rrn genes 16S (rrl), 23S (rrs) and 5S (rrf). Operon rrnD was inverted in orientation with respect to the others and contained a long 16S-23S intergenic spacer region with tRNAIle and tRNAAla genes, whereas the other operons contained a short spacer lacking any tRNA genes. The high-resolution physical/genetic map constructed in this study provides a platform for genomic and genetic studies of Lactobacillus species and for improving industrial and probiotic strains.

Characterization of Paenibacillus popilliae rRNA operons

The terminal 39 nucleotides on the 3' end of the 16S rRNA gene, along with the complete DNA sequences of the 5S rRNA, 23S rRNA, tRNA(Ile), and tRNA(Ala) genes were determined for Paenibacillus popilliae using strains NRRL B-2309 and Dutky 1. Southern hybridization analysis with a 16S rDNA hybridization probe and restriction-digested genomic DNA demonstrated 8 copies of the 16S rRNA gene in P. popilliae strains KLN 3 and Dutky 1. Additionally, the 23S rRNA gene in P. popilliae strains NRRL B-2309, KLN 3, and Dutky 1 was shown by I-CeuI digestion and pulsed-field gel electrophoresis of genomic DNA to occur as 8 copies. It was concluded that these 3 P. popilliae strains contained 8 rrn operons. The 8 operon copies were preferentially located on approximately one-half of the chromosome and were organized into 3 different patterns of genes, as follows: 16S-23S-5S, 16S-ala-23S-5S, and 16S-5S-ile-ala-23S-5S. This is the first report to identify a 5S rRNA gene between the 16S and 23S rRNA genes of a bacterial rrn operon. Comparative analysis of the nucleotides on the 3' end of the 16S rRNA gene suggests that translation of P. popilliae mRNA may occur in Bacillus subtilis and Escherichia coli.

The mechanism of upstream activation in the rrnB operon of Mycobacterium smegmatis is different from the Escherichia coli paradigm

Mycobacteria are slow-growing bacteria with a generation time of from 2-3 h up to several weeks. Consistent with the low growth rate, mycobacterial species have a maximum of two rRNA operons, rrnA and rrnB. The rrnA operon is present in all mycobacteria and has between two and five promoters, depending on species, whereas the rrnB operon, with a single promoter, is only found in some of the faster-growing species. The promoter region of the rrnB operon of a typical fast grower, Mycobacterium smegmatis, was investigated. By using lacZ reporter gene fusions it was demonstrated that the rrnB operon contains a highly activating region upstream of the core promoter, comparable to other bacterial rrn operons. However, the results suggest that, unlike the situation in, for example, Escherichia coli, the activating mechanism is solely factor dependent, and that no UP element is involved.

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.

Nature of the genome of the saprophytic spirochete Spirochaeta aurantia and its ribosomal RNA operons

Using restriction endonucleases DraI, AseI, and I-CeuI in conjunction with pulsed-field gel electrophoresis, we have shown that Spirochaeta aurantia M1 possesses a circular 3.98-Mb genome. This is the second largest spirochete chromosome yet analyzed. The observation that the latter enzyme cuts in 3 places suggests the presence of 3 copies of the large subunit (23S) rRNA gene (rrl), which was confirmed by Southern hybridizations. The complete sequence of 2 of the ribosomal RNA operons was determined, revealing that their structure resembled that of the typical member of the bacterial superkingdom: rrs (16S; 1561 bp), tRNA, rrl (23S; 2972 bp), and rrf (5S; 110 bp). The S. aurantia rrs-rrl intergenic regions, as with Treponema denticola, contain genes specifying a 73-nt tRNA(Ala) (anticodon TGC) and a 77-nt tRNA(Ile) (anticodon GAT).

Single-nucleotide polymorphisms in the rRNA operon and variable numbers of tandem repeats in the lipoprotein gene among Mycoplasma genitalium strains from clinical specimens

We characterized the single-nucleotide polymorphisms in the rRNA operon and variable numbers of tandem repeats in the lipoprotein gene MG309 among Mycoplasma genitalium strains from clinical specimens by PCR and sequencing. Analysis of 31 M. genitalium-infected patient specimens and 7 American Type Culture Collection strains identified six types of rRNA sequences and 11 different numbers of MG309 repeats. Examination of sequential specimens from 10 patients showed that these genotypes were stable for at least 5 weeks. These data suggest the potential usefulness of the rRNA genotypes and the MG309 repeats for genotyping of M. genitalium.

Characterization of an anonymous molecular marker strongly linked to Escherichia coli strains causing neonatal meningitis

An anonymous 14.9-kb rrn-containing HindIII fragment is strongly linked to Escherichia coli strains causing neonatal meningitis. We show in this report that this fragment does not encode new virulence factors but lacks arpA, a gene common in avirulent E. coli strains, and we developed a PCR test to detect this fragment.

Comparative analysis of the four rRNA operons in Finegoldia magna ATCC29328

There are four rRNA operons rrnA, rrnB, rrnC and rrnD on the genome of Finegoldia magna (formerly Peptostreptococcus magnus) ATCC29328, which, in contrast to those of Clostridia, are dispersed around the chromosome. Using a BAC library we determined the nucleotide sequences and structures of all four operons, including their flanking regions, and performed comparative analyses. We identified putative boxA sequences in the operons, which should be required for rRNA transcription antitermination, as well as their respective tandem promoters, AT-rich UP elements in the upstream region and Rho-independent terminators in the downstream region. The mosaic features of the operons were revealed. Multiple tRNAs were identified in the downstream region of two operons, 18 in rrnC and 11 in rrnD. They were presumed to form transcription units together with rRNAs. rrnA and rrnB had repeat units with Rho-independent terminators instead of tRNAs in the downstream region. rrnB and rrnC were the most similar in rrn upstream promoter region. Focusing on the sequence variations of rRNA genes, rrnB alone was heterogeneous. In light of previous reports, we also assessed the correlation between intercistronic rRNA sequence differences and distances between the operons, but no positive correlation was seen in this strain.

Quantitative relationships for specific growth rates and macromolecular compositions of Mycobacterium tuberculosis, Streptomyces coelicolor A3(2) and Escherichia coli B/r: an integrative theoretical approach

Further understanding of the physiological states of Mycobacterium tuberculosis and other mycobacteria was sought through comparisons with the genomic properties and macromolecular compositions of Streptomyces coelicolor A3(2), grown at 30 °C, and Escherichia coli B/r, grown at 37 °C. A frame of reference was established based on quantitative relationships observed between specific growth rates (μ) of cells and their macromolecular compositions. The concept of a schematic cell based on transcription/translation coupling, average genes and average proteins was developed to provide an instantaneous view of macromolecular synthesis carried out by cells growing at their maximum rate. It was inferred that the ultra-fast growth of E. coli results from its ability to increase the average number of rRNA (rrn) operons per cell through polyploidy, thereby increasing its capacity for ribosome synthesis. The maximum growth rate of E. coli was deduced to be limited by the rate of uptake and consumption of nutrients providing energy. Three characteristic properties of S. coelicolor A3(2) growing optimally (μ=0·30 h−1) were identified. First, the rate of DNA replication was found to approach the rate reported for E. coli (μ=1·73 h−1); secondly, all rrn operons were calculated to be fully engaged in precursor-rRNA synthesis; thirdly, compared with E. coli, protein synthesis was found to depend on higher concentrations of ribosomes and lower concentrations of aminoacyl-tRNA and EF-Tu. An equation was derived for E. coli B/r relating μ to the number of rrn operons per genome. Values of μ=0·69 h−1 and μ=1·00 h−1 were obtained respectively for cells with one or two rrn operons per genome. Using the author's equation relating the number of rrn operons per genome to maximum growth rate, it is expected that M. tuberculosis with one rrn operon should be capable of growing much faster than it actually does. Therefore, it is suggested that the high number of insertion sequences in this species attenuates growth rate to still lower values.

Gene expression in Chromobacterium violaceum

The repertoire of 4,431 open reading frames (ORFs), eight rRNA operons and 98 tRNA genes of Chromobacterium violaceum must be expressed in a regulated manner for successful adaptation to a wide variety of environmental conditions. To accomplish this feat, the organism relies on protein machineries involved in transcription, RNA processing and translation. Analysis of the C. violaceum genome showed that transcription initiation, elongation and termination are performed by the five well-known RNA polymerase subunits, five categories of sigma 70 factors, one sigma 54 factor, as well as six auxiliary elongation and termination factors. RNA processing is performed by a variety of endonucleases and exonucleases, such as ribonuclease H, ribonuclease E, ribonuclease P, and ribonuclease III, in addition to poly(A) polymerase and specific methyltransferases and pseudouridine synthases. ORFs for all ribosomal proteins, except S22, were found. Only 19 aminoacyl-tRNA synthetases were found, in addition to three aminoacyl-tRNA synthetase-related proteins. Asparaginyl-tRNA (Asn) is probably obtained by enzymatic modification of a mischarged aminoacyl-tRNA. The translation factors IF-1, IF-2, IF-3, EF-Ts, EF-Tu, EF-G, RF-1, RF-2 and RF-3 are all present in the C. violaceum genome, although the absence of selB suggests that C. violaceum does not synthesize selenoproteins. The components of trans-translation, tmRNA and associated proteins, are present in the C. violaceum genome. Finally, a large number of ORFs related to regulation of gene expression were also found, which was expected, considering the apparent adaptability of this bacterium.

Evaluation of the rrn operon copy number in Bifidobacterium using real-time PCR

AIMS

A real-time PCR-based method was developed to evaluate the Bifidobacterium rRNA operon copy number. As a result of their repetitive nature, rRNA operons are very suitable targets for chromosomal integration of heterologous genes.

METHODS AND RESULTS

The rrn operon multiplicity per chromosome was determined by real-time PCR quantification of the 16S rRNA amplicons obtained from genomic DNA. The values obtained in several bifidobacterial strains of human origin ranged from 1 to 5. The reliability of the method developed was confirmed by Southern hybridization technique.

CONCLUSIONS

In the Bifidobacterium genus the rrn operon copies showed variability at species and strain level. The identification of Bifidobacterium strains with high rRNA multiplicity allowed the selection of potential hosts for chromosomal integration.

SIGNIFICANCE AND IMPACT OF THE STUDY

The methodology here proposed represents a rapid, reliable and sensitive new tool for the quantification of rrn operon copy number in bacteria.

Fis stabilizes the interaction between RNA polymerase and the ribosomal promoter rrnB P1, leading to transcriptional activation

It has been shown that Fis activates transcription of the ribosomal promoter rrnB P1; however, the mechanism by which Fis activates rrnB P1 transcription is not fully understood. Paradoxically, although Fis activates transcription of rrnB P1 in vitro, transcription from the promoter containing Fis sites (as measured from rrnB P1-lacZ fusions) is not reduced in a fis null mutant strain. In this study, we further investigated the mechanism by which Fis activates transcription of the rrnB P1 promoter and the role of Fis in rRNA synthesis and cell growth in Escherichia coli. Like all other stringent promoters investigated so far, open complex of rrnB P1 has been shown to be intrinsically unstable, making open complex stability a potential regulatory step in transcription of this class of promoters. Our results show that Fis acts at this regulatory step by stabilizing the interaction between RNA polymerase and rrnB P1 in the absence of NTPs. Mutational analysis of the Fis protein demonstrates that there is a complete correlation between Fis-mediated transcriptional activation of rrnB P1 and Fis-mediated stabilization of preinitiation complexes of the promoter. Thus, our study indicates that Fis-mediated stabilization of RNA polymerase-rrnB P1 preinitiation complexes, presumably at the open complex step, contributes prominently to transcriptional activation. Furthermore, our in vivo results show that rRNA synthesis from the P1 promoters of several rRNA operons are reduced 2-fold in a fis null mutant compared with the wild type strain, indicating that Fis plays an important role in the establishment of robust rRNA synthesis when E. coli cells are emerging from a growth-arrested phase to a rapid growth phase. Thus, our results resolve an apparent paradox of the role of Fis in vitro and in vivo in the field.

Atypical 16S rRNA gene copies in Ochrobactrum intermedium strains reveal a large genomic rearrangement by recombination between rrn copies

Ochrobactrum intermedium is an opportunistic human pathogen belonging to the alpha 2 subgroup of proteobacteria. The 16S rDNA sequences of nine O. intermedium isolates from a collection of clinical and environmental isolates exhibited a 46-bp insertion at position 187, which was present in only one sequence among the 82 complete or partial 16S rDNA sequences of Ochrobactrum spp. available in data banks. Reverse transcription-PCR experiments showed that the 46-bp insertion remained in the 16S rRNA. The inserted sequence folded into a stem-loop structure, which took place in and prolonged helix H184 of the 16S rRNA molecule. Helix H184 has been described as conserved in length among eubacteria, suggesting the idiosyncratic character of the 46-bp insertion. Pulsed-field gel electrophoresis experiments showed that seven of the clinical isolates carrying the 46-bp insertion belonged to the same clone. Insertion and rrn copy numbers were determined by hybridization and I-CeuI digestion. In the set of clonal isolates, the loss of two insertion copies revealed the deletion of a large genomic fragment of 150 kb, which included one rrn copy; deletion occurred during the in vivo evolution of the clone. Determination of the rrn skeleton suggested that the large genomic rearrangement occurred during events involving homologous recombination between rrn copies. The loss of insertion copies suggested a phenomenon of concerted evolution among heterogeneous rrn copies.

Understanding bias in microbial community analysis techniques due to rrn operon copy number heterogeneity

Molecular tools based on rRNA (rrn) genes are valuable techniques for the study of microbial communities. However, the presence of operon copy number heterogeneity represents a source of systematic error in community analysis. To understand the types and magnitude of such bias, four commonly used rrn-based techniques were used to perform an in silico analysis of a hypothetical community comprised organisms from the Comprehensive Microbial Resource database. Community profiles were generated, and diversity indices were calculated for length heterogeneity PCR, automated ribosomal integenic spacer analysis, denaturing gradient gel electrophoresis, and terminal RFLP (using RsaI, MspI, and HhaI). The results demonstrate that all techniques present a quantitative bias toward organisms with higher copy numbers. In addition, techniques may underestimate diversity by grouping similar ribotypes or overestimate diversity by allowing multiple signals for one organism. The results of this study suggest that caution should be used when interpreting rrn-based community analysis techniques.

Eyeing bacterial genomes

The density of information in a bacterial genome allows its history, organization and encoded functions to be distilled into a single graphical representation. These features have made it possible to discern the forces acting in and on bacterial genomes at levels not attainable in eukaryotes.

Analysis of conserved non-rRNA genes of Tropheryma whipplei

The causative agent of Whipple's disease, Tropheryma whipplei, is a slow-growing bacterium that remains poorly-understood. Genetic characterization of this organism has relied heavily upon rRNA sequence analysis. Pending completion of a complete genome sequencing effort, we have characterized several conserved non-rRNA genes from T. whipplei directly from infected tissue using broad-range PCR and a genome-walking strategy. Our goals were to evaluate its phylogenetic relationships, and to find ways to expand the strain typing scheme, based on rDNA sequence comparisons. The genes coding for the ATP synthase beta subunit (atpD), elongation factor Tu (tuf), heat shock protein GroEL (groEL), beta subunit of DNA-dependent RNA polymerase (rpoB), and RNase P RNA (rnpB) were analyzed, as well as the regions upstream and downstream of the rRNA operon. Phylogenetic analyses with all non-rRNA marker molecules consistently placed T. whipplei within the class, Actinobacteria. The arrangement of genes in the atpD and rpoB chromosomal regions was also consistent with other actinomycete genomes. Tandem sequence repeats were found upstream and downstream of the rRNA operon, and downstream of the groEL gene. These chromosomal sites and the 16S-23S rRNA intergenic spacer regions were examined in the specimens of 11 patients, and a unique combination of tandem repeat numbers and spacer polymorphisms was found in each patient. These data provide the basis for a more discriminatory typing method for T. whipplei.

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.

Unique architecture of the plastid ribosomal RNA operon promoter recognized by the multisubunit RNA polymerase in tobacco and other higher plants

Expression of the plastid rRNA operon (rrn) during development is highly regulated at the level of transcription. The plastid rrn operon in most higher plants is transcribed by the plastid-encoded RNA polymerase (PEP), the multisubunit plastid RNA polymerase from PrrnP1, a sigma(70)-type promoter with conserved -10 and -35 core promoter elements. To identify functionally important sequences, the tobacco PrrnP1 was dissected in vivo and in vitro. Based on in vivo deletion analysis, sequences upstream of nucleotide -83 do not significantly contribute to promoter function. The in vitro analyses identified an essential hexameric sequence upstream of the -35 element (GTGGGA; the rRNA operon upstream activator [RUA]) that is conserved in monocot and dicot species and suggested that the -10 element plays only a limited role in PrrnP1 recognition. Mutations in the initial transcribed sequence (+9 to +14) enhanced transcription, the characteristic of strong promoters in prokaryotes. We propose that sigma interaction with the -10 element in PrrnP1 is replaced in part by direct PEP-RUA (protein-DNA) interaction or by protein-protein interaction between the PEP and an RUA binding transcription factor.

Regulation of the Escherichia coli rrnB P2 promoter

The seven rRNA operons in Escherichia coli each contain two promoters, rrn P1 and rrn P2. Most previous studies have focused on the rrn P1 promoters. Here we report a systematic analysis of the activity and regulation of the rrnB P2 promoter in order to define the intrinsic properties of rrn P2 promoters and to understand better their contributions to rRNA synthesis when they are in their natural setting downstream of rrn P1 promoters. In contrast to the conclusions reached in some previous studies, we find that rrnB P2 is regulated: it displays clear responses to amino acid availability (stringent control), rRNA gene dose (feedback control), and changes in growth rate (growth rate-dependent control). Stringent control of rrnB P2 requires the alarmone ppGpp, but growth rate-dependent control of rrnB P2 does not require ppGpp. The rrnB P2 core promoter sequence (-37 to +7) is sufficient to serve as the target for growth rate-dependent regulation.

Chloroplast transformation with modified accD operon increases acetyl-CoA carboxylase and causes extension of leaf longevity and increase in seed yield in tobacco

Acetyl-CoA carboxylase (ACCase) in plastids is a key enzyme regulating the rate of de novo fatty acid biosynthesis in plants. Plastidic ACCase is composed of three nuclear-encoded subunits and one plastid-encoded accD subunit. To boost ACCase levels, we examined whether overexpression of accD elevates ACCase production. Using homologous recombination, we replaced the promoter of the accD operon in the tobacco plastid genome with a plastid rRNA-operon (rrn) promoter that directs enhanced expression in photosynthetic and non-photosynthetic organs, and successfully raised the total ACCase levels in plastids. This result suggests that the level of the accD subunit is a determinant of ACCase levels, and that enzyme levels are in part controlled post-transcriptionally at the level of subunit assembly. The resultant transformants grew normally and the fatty acid content was significantly increased in leaves, but not significantly in seeds. However, the transformants displayed extended leaf longevity and a twofold increase of seed yield over the control value, which eventually almost doubled the fatty acid production per plant of the transformants relative to control and wild-type plants. These findings offer a potential method for raising plant productivity and oil production.

Ribosomal and non-ribosomal resistance to oxazolidinones: species-specific idiosyncrasy of ribosomal alterations

A derivative of Mycobacterium smegmatis, which carries only one functional rRNA (rrn) operon, was used to isolate mutants resistant to the ribosome-targeted antibiotic linezolid. Isolation and characterization of linezolid-resistant clones revealed two classes of mutants. Ribosomes from class I mutants are resistant to oxazolidinones in an in vitro peptidyl transferase assay, indicating that resistance maps to the ribosome component. In contrast, ribosomes from class II mutants show wild-type susceptibility to a linezolid derivative in vitro, pointing to a non-ribosomal mechanism of resistance. Introduction of a wild-type ribosomal RNA operon into linezolid-resistant strains restored linezolid sensitivity in class I mutants, indicating that resistance (i) maps to the rRNA and (ii) is recessive. Sequencing of the entire rrn operon identified a single nucleotide alteration in 23S rRNA of class I mutant strains, 2447G --> T (Escherichia coli numbering). Introduction of mutant rrl2447T into M. smegmatis rrn- resulted in a linezolid-resistant phenotype, demonstrating a cause-effect relationship of the 2447G --> T alteration. The 2447G --> T mutation, which renders M. smegmatis linezolid resistant, confers lethality in E. coli. This finding is strong evidence of structural and pos-sibly functional differences between the ribosomes of Gram-positive and Gram-negative bacteria. In agreement with the results of the in vitro assay, class II mutants show a wild-type sequence of the complete rRNA operon. The lack of cross-resistance of the class II mutants to other antibiotics suggests a resistance mechanism other than activation of a broad-spectrum multidrug transporter.

The genome structure of Pseudomonas putida: high-resolution mapping and microarray analysis

As part of a collaborative project aimed at sequencing and functionally analysing the entire genome of Pseudomonas putida strain KT2440, a physical clone map was produced as an initial resource. To this end, a high-coverage cosmid library was arrayed and ordered by clone hybridizations. Restriction fragments generated by rare-cutting enzymes and plasmids containing the rrn operon and 23S rDNA of Pseudomonas aeruginosa were used as probes and, parts of the cosmids were end-sequenced. This provided the information necessary for merging and comparing the macro-restriction map, cosmid clone order and sequence information, thereby assuring co-linearity of the eventual sequence assembly with the actual genome. A tiling path of clones was selected, from the shotgun clones used for sequencing, for the production of DNA microarrays that represent the entire genome including its non-coding portions.

Interaction of avilamycin with ribosomes and resistance caused by mutations in 23S rRNA

The antibiotic growth promoter avilamycin inhibits protein synthesis by binding to bacterial ribosomes. Here the binding site is further characterized on Escherichia coli ribosomes. The drug interacts with domain V of 23S rRNA, giving a chemical footprint at nucleotides A2482 and A2534. Selection of avilamycin-resistant Halobacterium halobium cells revealed mutations in helix 89 of 23S rRNA. Furthermore, mutations in helices 89 and 91, which have previously been shown to confer resistance to evernimicin, give cross-resistance to avilamycin. These data place the binding site of avilamycin on 23S rRNA close to the elbow of A-site tRNA. It is inferred that avilamycin interacts with the ribosomes at the ribosomal A-site interfering with initiation factor IF2 and tRNA binding in a manner similar to evernimicin.

Sequence diversity of the internal transcribed spacer (ITS) region of the rRNA operons among different serogroups of Legionella pneumophila isolates

The genus Legionella is represented by 48 species and Legionella pneumophila includes 15 serogroups. In this work, we have studied the intergenic 16S-23S spacer region (ITS) in L. pneumophila to determine the feasability of using amplification polymorphisms in this region, to establish intraspecies differences and to discriminate Legionella species. The amplification of this region, using 16S14F and 23S0R primers, and the analysis of amplicons by the analysis of fragments technique showed that all the L. pneumophila serogroups studied presented the same electrophoretic pattern. Moreover, the analysis of different Legionella species showed that the amplification polymorphisms obtained were species-specific. In order to study the sequence variability of this region, the existence in L. pneumophila of three ribosomal operons was determined by pulsed field gel eletrophoresis (PFGE). Two of the 16S-23S rRNA ITS presented a tRNA Ala and the third one a tRNA Ile. Nevertheless, the variability expected in this region of the operon was not found and the rest of the ITS contained only punctual mutations.

[Transcription regulation in the ribosomal RNA operon of mycobacteria]

Tuberculosis is a extremely important infectious disease, caused by the bacilli Mycobacterium tuberculosis. One of the characteristic of this bacteria is its very slow rate of growth, that allows it to survive for long periods of time inside the host cells. Among the genetic elements involved in growth regulation the operon rrn is of extreme importance. This operon contains the genes that code the three rRNA molecules, essential components of the bacterial ribosome. The tuberculosis bacilli, differently from most of the microorganisms, has a single copy of the rrn operon per genome, meaning that it must be submitted to very strict control mechanisms. Another important conclusion is that the sequences of the rrn operon constitute ideal targets for anti-mycobacterial drugs. In this work we have studied some of the elements involved in transcription control in M. tuberculosis, particularly those present in the leader region of the operon. By using basic molecular biology techniques we have identified sequence elements in the leader region that seem to be involved in the control of transcription elongation, by a mechanism related to anti-termination.

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.