Assembly of a complete genome sequence for <Emphasis Type="Italic">Gemmata obscuriglobus</Emphasis> reveals a novel prokaryotic rRNA operon gene architecture

Gemmata obscuriglobus is a Gram-negative bacterium with several intriguing biological features. Here, we present a complete, de novo whole genome assembly for G. obscuriglobus which consists of a single, circular 9 Mb chromosome, with no plasmids detected. The genome was annotated using the NCBI Prokaryotic Genome Annotation pipeline to generate common gene annotations. Analysis of the rRNA genes revealed three interesting features for a bacterium. First, linked G. obscuriglobus rrn operons have a unique gene order, 23S–5S–16S, compared to typical prokaryotic rrn operons (16S–23S–5S). Second, G. obscuriglobus rrn operons can either be linked or unlinked (a 16S gene is in a separate genomic location from a 23S and 5S gene pair). Third, all of the 23S genes (5 in total) have unique polymorphisms. Genome analysis of a different Gemmata species (SH-PL17), revealed a similar 23S–5S–16S gene order in all of its linked rrn operons and the presence of an unlinked operon. Together, our findings show that unique and rare features in Gemmata rrn operons among prokaryotes provide a means to better define the evolutionary relatedness of Gemmata species and the divergence time for different Gemmata species. Additionally, these rrn operon differences provide important insights into the rrn operon architecture of common ancestors of the planctomycetes.     

Intra- and Intergenomic Variation of Ribosomal RNA Operons in Concurrent Alteromonas macleodii Strains

Biodiversity estimates based on ribosomal operon sequence diversity rely on the premise that a sequence is characteristic of a single specific taxon or operational taxonomic unit (OTU). Here, we have studied the sequence diversity of 14 ribosomal RNA operons (rrn) contained in the genomes of two isolates (five operons in each genome) and four metagenomic fosmids, all from the same seawater sample. Complete sequencing of the isolate genomes and the fosmids establish that they represent strains of the same species, Alteromonas macleodii, with average nucleotide identity (ANI) values >97 %. Nonetheless, we observed high levels of intragenomic heterogeneity (i.e., variability between operons of a single genome) affecting multiple regions of the 16S and 23S rRNA genes as well as the internally transcribed spacer 1 (ITS-1) region. Furthermore, the ribosomal operons exhibited intergenomic heterogeneity (i.e., variability between operons located in separate genomes) in each of these regions, compounding the variability. Our data reveal the extensive heterogeneity observed in natural populations of A. macleodii at a single point in time and support the idea that distinct lineages of A. macleodii exist in the deep Mediterranean. These findings highlight the potential of rRNA fingerprinting methods to misrepresent species diversity while simultaneously failing to recognize the ecological significance of individual strains.     

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. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Reviewing Editor: Dr. Margaret Riley     

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.     

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.     

Sequence homologies between Escherichia coli and chloroplast ribosomal DNA as seen by heteroduplex analysis

Circular Vicia faba (broad bean) chloroplast DNA was hybridized to the restriction fragment BamHI B from the DNA of the transducing phage lambda rif^d18, which carries the Escherichia coli ribosomal RNA operon rrnB. Cytochrome spreadings of the heteroduplexes show homologies in the 16 S and 23 S rRNA regions, but none in the spacer. The same lambda rif^d18 fragment was hybridized to the Vicia cpDNA ²SalI fragment 3, which contains the Vicia rBNA operon, resulting in an analogous heteroduplex configuration. Cytochrome spreadings of this heteroduplex in increasing concentrations of formamide reveal regions of incomplete homologies. Heteroduplexes between the E. coli rrnD operon, obtained from the recombinant plasmid pBK8, and circular Vicia cpDNA revealed homologies in the spacer region as well as in the 16 S and 23 S rRNA region. Hybrids between all three types of rDNA and their homologous rRNAs were prepared using the mica adsorption technique. They show that the 23 S, 16 S, and 5 S rRNAs are transcribed from the same strand of Vicia cpDNA. The positions of the rRNAs were measured and compared to the heteroduplex structure. It was observed that the E. coli rrnD operon in the plasmid pBK8 contains two 5 S rRNA sequences near the distal end.     

Sequence variation in the putative replication origins of the five genetically distinct basicEuoenothera plastid chromosomes (plastomes)

Approximately 4200 nucleotides of the 16S/23S rDNA spacer and the 5′ region flanking therrn operon of the plastid chromosomes representing the five basic, phylogenetically relatedEuoenothera plastomes were sequenced and compared. The sequences that harbor the putative replication origins are almost identical except for a 785-bp intercistronic segment between the genes for the 16S rRNA andtrnI. Differences are mainly caused by insertions/deletions and duplications; the predicted potential for formation of quite extensive secondary structure differs among the plastomes. Unexpected intraplastome variation has also been noted. Furthermore, the sequence-based and published genetically deduced plastome pedigrees differ significantly.     

Sequence and characterisation of a ribosomal RNA operon fromAgrobacterium vitis

One of the four ribosomal RNA operons (rrnA) from theAgrobacterium vitis vitopine strain S4 was sequenced.rrnA is most closely related to therrn operons ofBradyrhizobium japonicum andRhodobacter sphaeroides and carries an fMet-tRNA gene downstream of its 5S gene, as in the case ofR. sphaeroides. The 16S rRNA sequence of S4 differs from theA. vitis K309 type strain sequence by only one nucleotide, in spite of the fact that S4 and K309 have very different Ti plasmids. The predicted secondary structure of the S4 23S rRNA shows several features that are specific for the alpha proteobacteria, and an unusual branched structure in the universal B8 stem. The 3′ ends of the three otherrrn copies of S4 were also cloned and sequenced. Sequence comparison delimits the 3′ ends of the four repeats and defines two groups:rrnA/rrnB andrrnC/rrnD.     

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.     

Structure and rearrangements of rRNA genes in chloroplast DNA in two strains of Euglena gracilis

Summary The organisation of the rRNA genes in the chloroplast genomes of two strains of Euglena gracilis were analyzed and compared. It was previously shown that the bacillaris strain contains three complete rrn (rRNA) operons (7) and that the Z-S strain contains one operon (21). Using heteroduplex analysis it was found that the bacillaris strain contains, apart from the three complete rrn operons, an extra 16S rRNA gene, an extra partial 23S rRNA gene sequence and an inverted duplication of a stretch within the 5S–16S spacer. In addition a short (<100 bp) inverted repeat sequence (13) which forms a stem/loop structure in single-stranded cpDNA was located between the 3-end of the extra 16S rRNA gene and the partial 23 S rRNA sequence.The Z-S strain differs from the bacillaris strain by a deletion of two units of the complete rrn operons. The region upstream of the single complete rrn operon, including the inverted repeats, the partial 23S and the extra 16S rRNA sequences is identical with the bacillaris strain.The only non-homology found in heteroduplexes between the SalI fragments of B of the two strains is the deletion-insertion loop which represents the two rrn operons. A small deletion loop was found occasionally in hetero-and in homoduplexes of both strands in the region of variable size. Apart from the deletion/insertion of two rrn operons the two genomes appear to be colinear as can be seen from partial denaturation mapping. The organisation of the rRNA genes of the two strains is compared with those of the Z strain and the bacillaris-ATCC strain.     

Isolation of Spectinomycin Resistance Mutations in the 16S rRNA of Salmonella enterica serovar Typhimurium and Expression in Escherichia coli and Salmonella

Two single-base mutations in 16S rRNA conferring high-level resistance to spectinomycin were isolated on a plasmid-borne copy of the rrnD operon from Salmonella enterica serovar Typhimurium. Neither of the mutations (C1066U and C1192U) had appreciable effects on cell growth, but each had differential effects on resistance to spectinomycin and fusidic acid. Both mutations also conferred resistance to spectinomycin in Escherichia coli strains containing deletions of all seven chromosomal rrn operons and expressing plasmid-encoded Salmonella rRNA exclusively. In contrast, when expressed in E. coli strains containing intact chromosomal rrn operons, the strains were sensitive to spectinomycin. However, chromosomal mutations arose that allowed expression of the rRNA-dependent spectinomycin resistance phenotype. It is proposed that in heterogeneous rRNA populations, the native E. coli rRNA out-competes the heterologous Salmonella rRNA for binding to ribosomal proteins, translation factors, or ribosome assembly, thus limiting entry of the antibiotic-resistant 30S subunits into the functioning ribosome pool. Received: 28 September 2001 / Accepted: 26 March 2002     

Ribosomal RNAs are tolerant toward genetic insertions: evolutionary origin of the expansion segments

Ribosomal RNAs (rRNAs), assisted by ribosomal proteins, form the basic structure of the ribosome, and play critical roles in protein synthesis. Compared to prokaryotic ribosomes, eukaryotic ribosomes contain elongated rRNAs with several expansion segments and larger numbers of ribosomal proteins. To investigate architectural evolution and functional capability of rRNAs, we employed a Tn5 transposon system to develop a systematic genetic insertion of an RNA segment 31 nt in length into Escherichia coli rRNAs. From the plasmid library harboring a single rRNA operon containing random insertions, we isolated surviving clones bearing rRNAs with functional insertions that enabled rescue of the E. coli strain (Δ7rrn) in which all chromosomal rRNA operons were depleted. We identified 51 sites with functional insertions, 16 sites in 16S rRNA and 35 sites in 23S rRNA, revealing the architecture of E. coli rRNAs to be substantially flexible. Most of the insertion sites show clear tendency to coincide with the regions of the expansion segments found in eukaryotic rRNAs, implying that eukaryotic rRNAs evolved from prokaryotic rRNAs suffering genetic insertions and selections.     

Phylogenetic analysis of nitrite, nitric oxide, and nitrous oxide respiratory enzymes reveal a complex evolutionary history for denitrification

Denitrification is a facultative respiratory pathway in which nitrite (NO2(-)), nitric oxide (NO), and nitrous oxide (N2O) are successively reduced to nitrogen gas (N(2)), effectively closing the nitrogen cycle. The ability to denitrify is widely dispersed among prokaryotes, and this polyphyletic distribution has raised the possibility of horizontal gene transfer (HGT) having a substantial role in the evolution of denitrification. Comparisons of 16S rRNA and denitrification gene phylogenies in recent studies support this possibility; however, these results remain speculative as they are based on visual comparisons of phylogenies from partial sequences. We reanalyzed publicly available nirS, nirK, norB, and nosZ partial sequences using Bayesian and maximum likelihood phylogenetic inference. Concomitant analysis of denitrification genes with 16S rRNA sequences from the same organisms showed substantial differences between the trees, which were supported by examining the posterior probability of monophyletic constraints at different taxonomic levels. Although these differences suggest HGT of denitrification genes, the presence of structural variants for nirK, norB, and nosZ makes it difficult to determine HGT from other evolutionary events. Additional analysis using phylogenetic networks and likelihood ratio tests of phylogenies based on full-length sequences retrieved from genomes also revealed significant differences in tree topologies among denitrification and 16S rRNA gene phylogenies, with the exception of the nosZ gene phylogeny within the data set of the nirK-harboring genomes. However, inspection of codon usage and G + C content plots from complete genomes gave no evidence for recent HGT. Instead, the close proximity of denitrification gene copies in the genomes of several denitrifying bacteria suggests duplication. Although HGT cannot be ruled out as a factor in the evolution of denitrification genes, our analysis suggests that other phenomena, such gene duplication/divergence and lineage sorting, may have differently influenced the evolution of each denitrification gene.     

Identification of Cold Shock Gene Loci in Sinorhizobium meliloti by Using a luxAB Reporter Transposon

Using a luxAB reporter transposon, seven mutants of Sinorhizobium meliloti were identified as containing insertions in four cold shock loci. LuxAB activity was strongly induced (25- to 160-fold) after a temperature shift from 30 to 15°C. The transposon and flanking host DNA from each mutant was cloned, and the nucleic acid sequence of the insertion site was determined. Unexpectedly, five of the seven luxAB mutants contained transposon insertions in the 16S and 23S rRNA genes of two of the three rrn operons of S. meliloti. Directed insertion of luxAB genes into each of the three rrn operons revealed that all three operons were similarly affected by cold shock. Two other insertions were found to be located downstream of a homolog of the major Escherichia coli cold shock gene, cspA. Although the cold shock loci were highly induced in response to a shift to low temperature, none of the insertions resulted in a statistically significant decrease in growth rate at 15°C.     

Polyphasic analysis reveals correlation between phenotypic and genotypic analysis in soybean bradyrhizobia (Bradyrhizobium spp.)

Soybean bradyrhizobia (Bradyrhizobium spp.) are bacteria that fix atmospheric nitrogen within the root nodules of soybean, a crop critical for meeting global nutritional protein demand. Members of this group differ in symbiotic effectiveness, and historically both phenotypic and genotypic approaches have been used to assess bradyrhizobial diversity. However, agreement between various approaches of assessment is poorly known. A collection (n = 382) of soybean bradyrhizobia (Bradyrhizobium japonicum, B. diazoefficiens, and B. elkanii) were characterized by Internal Transcribed Spacer – Restriction Fragment Length Polymorphism (ITS-RFLP), cellular fatty acid composition (fatty acid methyl esters, FAME), and serological reactions to assess agreement between phenotypic and genotypic methods. Overall, 76% of the accessions demonstrated identical clustering with each of these techniques. FAME was able to identify all 382 accessions, whereas 14% were non-reactive serologically. One ITS-RFLP group, containing 36 Delaware isolates, produced multiple ITS amplicons indicating they possess multiple ribosomal RNA (rrn) operons. Cloning and sequencing revealed that these strains contained as many as three heterogenous rrn operons, a trait previously unknown in bradyrhizobia. A representative subset of 96 isolates was further characterized using 16S rRNA and Internal Transcribed Spacer (ITS) amplicon sequencing. ITS sequences showed better inter- and intra-species discrimination (65–99% identity) than 16S sequences (96–99% identity). This study shows that phenotypic and genotypic approaches are strongly correlated at the species level but should be approached with caution. We also suggest using combined 16S and ITS genotyping data to obtain better inter- and intra-species resolution in bradyrhizobia classification.