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.     

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

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

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.     

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     

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     

Intragenomic heterogeneity and intergenomic recombination among haloarchaeal rRNA genes

More than one copy of rRNA operons, which code for both the small-subunit (SSU) and large-subunit (LSU) rRNA, are often found in prokaryotes. It is generally assumed that all rRNA operons within a single cell are almost identical. A notable exception is the extremely halophilic archaeal genus Haloarcula, most species of which are known to harbor highly divergent rRNA operons that differ at approximately 5% of the nucleotide positions in the SSU gene and at 1 to 2% of the nucleotide positions in the LSU gene. We report that such intragenomic heterogeneity is not unique to Haloarcula, as high levels of intragenomic sequence variation have been observed for the SSU genes of two other genera of extreme halophiles, Halosimplex and Natrinema. To investigate this in detail, the two rRNA operons of Halosimplex carlsbadense and the four operons of Natrinema sp. strain XA3-1 were cloned and completely sequenced. The SSU and LSU genes of H. carlsbadense show the highest levels of intragenomic heterogeneity observed so far in archaea (6.7 and 2.6%). The operons of Natrinema sp. strain XA3-1 have additional unusual characteristics, such as identical internal transcribed spacers, while one of four SSU genes is 5% divergent and all LSU genes differ from each other by 0.9 to 1.9%. The heterogeneity among the Natrinema sp. strain XA3-1 LSU genes is localized in hot spots, and one of these regions is shown to be the result of a recombination event with a distantly related halophile. This is the first example of interspecies recombination between rRNA genes in archaea, and the recombination occurred over one of the largest phylogenetic distances ever reported for such an event. We suggest that intragenomic heterogeneity of rRNA operons is an ancient and stable trait in several lineages of the Halobacteriales. The impact of this phenomenon on the taxonomy of extremely halophilic archaea is discussed.     

Insights into the Origin of Clostridiumbotulinum Strains: Evolution of Distinct Restriction Endonuclease Sites in rrs (16S rRNA gene)

Diversity analysis of Clostridium botulinum strains is complicated by high microheterogeneity caused by the presence of 9–22 copies of rrs (16S rRNA gene). The need is to mine genetic markers to identify very closely related strains. Multiple alignments of the nucleotide sequences of the 212 rrs of 13 C. botulinum strains revealed intra- and inter-genomic heterogeneity. Low intragenomic heterogeneity in rrs was evident in strains 230613, Alaska E43, Okra, Eklund 17B, Langeland, 657, Kyoto, BKT015925, and Loch Maree. The most heterogenous rrs sequences were those of C. botulinum strains ATCC 19397, Hall, H04402065, and ATCC 3502. In silico restriction mapping of these rrs sequences was observable with 137 type II Restriction endonucleases (REs). Nucleotide changes (NC) at these RE sites resulted in appearance of distinct and additional sites, and loss in certain others. De novo appearances of RE sites due to NC were recorded at different positions in rrs gene. A nucleotide transition A>G in rrs of C. botulinum Loch Maree and 657 resulted in the generation of 4 and 10 distinct RE sites, respectively. Transitions A>G, G>A, and T>C led to the loss of RE sites. A perusal of the entire NC and in silico RE mapping of rrs of C. botulinum strains provided insights into their evolution. Segregation of strains on the basis of RE digestion patterns of rrs was validated by the cladistic analysis involving six house keeping genes: dnaN, gyrB, metG, prfA, pyrG, and Rho.

Electronic supplementary material

The online version of this article (doi:10.1007/s12088-015-0514-z) contains supplementary material, which is available to authorized users.     

Intragenomic 16S rDNA Divergence in <Emphasis Type="Italic">Haloarcula marismortui</Emphasis> Is an Adaptation to Different Temperatures

The halophilic archaeon Haloarcula marismortui contains three ribosomal RNA operons, designated rrnA, rrnB, and rrnC. Operons A and C are virtually identical, whereas operon B presents a high divergence in nucleotide sequence, having up to 135 nucleotide polymorphisms among the three 16S, 23S, and 5S ribosomal RNA genes. Quantitative PCR and structural analyses have been performed to elucidate whether the presence of this intragenomic heterogeneity could be an adaptation to the variable environmental conditions in the natural habitat of H. marismortui. Variation in salt concentration did not affect expression but variation in incubation temperature did produce significant changes, with operon B displaying an expression level four times higher than the other two together at 50°C and three times lower at 15°C. We show that the putative promoter region of operon B is also different. In addition, the predicted secondary structure of these genes indicated that they have distinct stabilities at different temperatures and a mutant strain lacking operon B grew slower at high temperatures. This study supports the idea that divergent rRNA genes can be adaptive, with different variants being functional under different environmental conditions (e.g., temperature). The same phenomenon could take place in other halophiles or thermophiles with intragenomic rDNA heterogeneity, where the use of 16S rDNA as a phylogenetic marker and indicator of biodiversity should be used with caution.     

Use of 16S rRNA and rpoB Genes as Molecular Markers for Microbial Ecology Studies

Several characteristics of the 16S rRNA gene, such as its essential function, ubiquity, and evolutionary properties, have allowed it to become the most commonly used molecular marker in microbial ecology. However, one fact that has been overlooked is that multiple copies of this gene are often present in a given bacterium. These intragenomic copies can differ in sequence, leading to identification of multiple ribotypes for a single organism. To evaluate the impact of such intragenomic heterogeneity on the performance of the 16S rRNA gene as a molecular marker, we compared its phylogenetic and evolutionary characteristics to those of the single-copy gene rpoB. Full-length gene sequences and gene fragments commonly used for denaturing gradient gel electrophoresis were compared at various taxonomic levels. Heterogeneity found between intragenomic 16S rRNA gene copies was concentrated in specific regions of rRNA secondary structure. Such “heterogeneity hot spots” occurred within all gene fragments commonly used in molecular microbial ecology. This intragenomic heterogeneity influenced 16S rRNA gene tree topology, phylogenetic resolution, and operational taxonomic unit estimates at the species level or below. rpoB provided comparable phylogenetic resolution to that of the 16S rRNA gene at all taxonomic levels, except between closely related organisms (species and subspecies levels), for which it provided better resolution. This is particularly relevant in the context of a growing number of studies focusing on subspecies diversity, in which single-copy protein-encoding genes such as rpoB could complement the information provided by the 16S rRNA gene.     

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.     

Intragenomic Heterogeneity of 16S rRNA Genes Causes Overestimation of Prokaryotic Diversity

Ever since Carl Woese introduced the use of 16S rRNA genes for determining the phylogenetic relationships of prokaryotes, this method has been regarded as the “gold standard” in both microbial phylogeny and ecology studies. However, intragenomic heterogeneity within 16S rRNA genes has been reported in many investigations and is believed to bias the estimation of prokaryotic diversity. In the current study, 2,013 completely sequenced genomes of bacteria and archaea were analyzed and intragenomic heterogeneity was found in 952 genomes (585 species), with 87.5% of the divergence detected being below the 1% level. In particular, some extremophiles (thermophiles and halophiles) were found to harbor highly divergent 16S rRNA genes. Overestimation caused by 16S rRNA gene intragenomic heterogeneity was evaluated at different levels using the full-length and partial 16S rRNA genes usually chosen as targets for pyrosequencing. The result indicates that, at the unique level, full-length 16S rRNA genes can produce an overestimation of as much as 123.7%, while at the 3% level, an overestimation of 12.9% for the V6 region may be introduced. Further analysis showed that intragenomic heterogeneity tends to concentrate in specific positions, with the V1 and V6 regions suffering the most intragenomic heterogeneity and the V4 and V5 regions suffering the least intragenomic heterogeneity in bacteria. This is the most up-to-date overview of the diversity of 16S rRNA genes within prokaryotic genomes. It not only provides general guidance on how much overestimation can be introduced when applying 16S rRNA gene-based methods, due to its intragenomic heterogeneity, but also recommends that, for bacteria, this overestimation be minimized using primers targeting the V4 and V5 regions.     

Evaluation of four molecular typing methodologies as tools for determining taxonomy relations and for identifying species among Yersinia isolates

In the last few decades, molecular typing has become an important tool in taxonomic, phylogenetic and identification studies of numerous groups of bacteria, including the yersiniae. In this study, Enterobacterial Repetitive Intergenic Consensus PCR (ERIC-PCR), Pulsed-Field Gel Electrophoresis (PFGE), 16S rRNA gene sequencing and Multilocus Sequence Analysis (MLSA) were performed to determine the ability of these techniques to be used in taxonomy and identification of Yersinia strains. A total of 60 Yersinia strains were genotyped by ERIC-PCR and PFGE. Moreover, an in silico analysis was carried out for 16S rRNA gene sequencing and MLSA, using 68 and 49 Yersinia strains, respectively. A phylogenetic tree constructed from the ERIC-PCR, 16S rRNA gene sequencing and MLSA data grouped most of the Yersinia species into distinct species-specific clusters. In the PFGE assay these clusters were not observed. On this basis, ERIC-PCR, 16S rRNA gene sequencing and MLSA seem to be valuable techniques for use in taxonomic and identification studies of the genus Yersinia, whereas PFGE does not. Furthermore, ERIC-PCR has the advantage of being a cheaper, easier and faster assay than 16S rRNA gene sequencing or MLSA, and for these reasons can be considerate an alternative tool in taxonomic studies of yersiniae.     

Detection of simple mutations and polymorphisms in large genomic regions

We have developed a novel technology that makes it possible to detect simple nucleotide polymorphisms directly within a sample of total genomic DNA. It allows, in a single Southern blot experiment, the determination of sequence identity of genomic regions with a combined length of hundreds of kilobases. This technology does not require PCR amplification of the target DNA regions, but exploits preparative size-fractionation of restriction-digested genomic DNA and a newly discovered property of the mismatch-specific endonuclease CEL I to cleave heteroduplex DNA with a very high specificity and sensitivity. We have used this technique to detect various simple mutations directly in the genomic DNA of isogenic pairs of recombinant Pseudomonas aeruginosa, Escherichia coli and Salmonella isolates. Also, by using a cosmid DNA library and genomic fractions as hybridization probes, we have compared total genomic DNA of two clinical P.aeruginosa clones isolated from the same patient, but exhibiting divergent phenotypes. The mutation scan correctly detected a GA insertion in the quorum-sensing regulator gene rhlR and, in addition, identified a novel intragenomic polymorphism in rrn operons, indicating very high stability of the bacterial genomes under natural non-mutator conditions.     

Comparative Analysis of Physical Maps of Four Bacillus subtilis (natto) Genomes

The complete SfiI and I-CeuI physical maps of four Bacillus subtilis (natto) strains, which were previously isolated as natto (fermented soybean) starters, were constructed to elucidate the genome structure. Not only the similarity in genome size and organization but also the microheterogeneity of the gene context was revealed. No large-scale genome rearrangements among the four strains were indicated by mapping of the genes, including 10 rRNA operons (rrn) and relevant genes required for natto production, to the loci corresponding to those of the B. subtilis strain Marburg 168. However, restriction fragment length polymorphism and the presence or absence of strain-specific DNA sequences, such as the prophages SPβ, skin element, and PBSX, as well as the insertion element IS4Bsu1, could be used to identify one of these strains as a Marburg type and the other three strains as natto types. The genome structure and gene heterogeneity were also consistent with the type of indigenous plasmids harbored by the strains.