Rates and consequences of recombination between rRNA operons

Examination of Large Chromosomal Inversions in the Genome of Erwinia amylovora Strains Reveals Worldwide Distribution and North America-Specific Types

Erwinia amylovora is a relatively homogeneous species with low genetic diversity at the nucleotide level. However, phenotypic differences and genomic structural variations among E. amylovora strains have been documented. In this study, we identified 10 large chromosomal inversion (LCI) types in the Spiraeoideae-infecting (SI) E. amylovora strains by combining whole genome sequencing and PCR-based molecular markers. It was found that LCIs were mainly caused by homologous recombination events among seven rRNA operons (rrns) in SI E. amylovora strains. Although ribotyping results identified inter- and intra-variations in the internal transcribed spacer (ITS1 and ITS2) regions among rrns, LCIs tend to occur between rrns transcribed in the opposite directions and with the same tRNA content (tRNA-Glu or tRNA-Ile/Ala) in ITS1. Based on the LCI types, physical/estimated replichore imbalance (PRI/ERI) was examined and calculated. Among the 117 SI strains evaluated, the LCI types of Ea1189, CFBP1430, and Ea273 were the most common, with ERI values at 1.31, 7.87, and 4.47°, respectively. These three LCI types had worldwide distribution, whereas the remaining seven LCI types were restricted to North America (or certain regions of the United States). Our results indicated ongoing chromosomal recombination events in the SI E. amylovora population and showed that LCI events are mostly symmetrical, keeping the ERI less than 15°. These findings provide initial evidence about the prevalence of certain LCI types in E. amylovora strains, how LCI occurs, and its potential evolutionary advantage and history, which might help track the movement of the pathogen.

Markerless bacterial artificial chromosome manipulation method by red proteins of phage λ mediated homologous recombination utilizing fluorescent proteins for both positive and counter selection

Manipulating viral genomes is an essential technique in reverse genetics and recombinant vaccine development. A strategy for manipulating large viral genomes involves introducing their entire genome into bacterial artificial chromosomes and employing Escherichia coli genetic tools. For sequence manipulation on bacterial artificial chromosomes (bacterial artificial chromosomes recombineering), a well-established method that relies on the Escherichia coli strain GS1783, and the template plasmid, pEPKan-S, is often used. This method, known as markerless DNA manipulation, allows for the generation of a recombinant bacterial artificial chromosome that does not retain the selection markers used during recombination. Although this method is highly innovative, there remains room for improvement as the plasmid is currently only available for positive selection. Additionally, differentiating true recombinants from false negatives often proves time-consuming. Consequently, an improved method for bacterial artificial chromosomes recombineering, which utilizes fluorescent proteins, has been developed. This method's core comprises three plasmids containing the I-SceI recognition site, antibiotic resistance genes (ampicillin, kanamycin, and zeocin), and fluorescent genes (YPet, mOrange, and mScarlet). The success or failure of Red recombination can be confirmed via fluorescent signals. To validate this method, the Lassa virus genes were introduced into the bacterial artificial chromosomes, containing the entire genome of the vaccinia virus strain LC16m8. Consequently, the expression of fluorescent protein genes contributed to positive selection, such as blue-white screening and counter-selection during the first and second Red recombination.

Improved site-specific mutagenesis in Rhodococcus opacus using a novel conditional suicide plasmid

Rhodococcus opacus PD630 is a biotechnologically important bacterium with metabolic capability for bioremediation, metal recovery, and storage of triacylglycerols. Genome editing by homologous recombination in R. opacus is hampered by a very low combined frequency of DNA transfer and recombination. To improve recombination in the species, a conjugative, conditional suicide plasmid based on the replicon derived from the Corynebacterium glutamicum plasmid pGA1 was constructed and evaluated in R. opacus. The replication of this plasmid is controlled by a dual inducible and repressible promoter system originally developed for Mycobacterium spp. Next, we demonstrated that a derivative of this plasmid containing sacB as a counterselection marker and homologous regions of R. opacus could be used for homologous recombination, and that the problem of obtaining recombinants had been solved. Like for other Corynebacteriales, the cell wall of Rhodococcus spp. contains mycolic acids which form a hydrophobic and impermeable outer layer. Mycolic acids are essential for Mycobacterium smegmatis, but not for Corynebacterium glutamicum, and the new vector was used to study if mycolic acid is essential for R. opacus. We found that accD3 that is necessary for mycolic acid synthesis could only be deleted from the chromosome in strains containing a plasmid-encoded copy of accD3. This indicates that mycolic acid is important for R. opacus viability. The conditional suicide vector should be useful for homologous recombination or for delivering gene products like recombinases or Cas proteins and gRNA to Rhodococcus and related genera, while the approach should be applicable for any plasmid needing a plasmid-encoded protein for replication. KEY POINTS: • Improved vector for homologous recombination in R. opacus. • Mycolic acid is important for survival of R. opacus like it is for Mycobacterium. • Similar conditional suicide plasmids may be constructed for other bacteria.

Phylogenies of the 16S rRNA gene and its hypervariable regions lack concordance with core genome phylogenies

BACKGROUND

The 16S rRNA gene is used extensively in bacterial phylogenetics, in species delineation, and now widely in microbiome studies. However, the gene suffers from intragenomic heterogeneity, and reports of recombination and an unreliable phylogenetic signal are accumulating. Here, we compare core gene phylogenies to phylogenies constructed using core gene concatenations to estimate the strength of signal for the 16S rRNA gene, its hypervariable regions, and all core genes at the intra- and inter-genus levels. Specifically, we perform four intra-genus analyses (Clostridium, n = 65; Legionella, n = 47; Staphylococcus, n = 36; and Campylobacter, n = 17) and one inter-genus analysis [41 core genera of the human gut microbiome (31 families, 17 orders, and 12 classes), n = 82].

RESULTS

At both taxonomic levels, the 16S rRNA gene was recombinant and subject to horizontal gene transfer. At the intra-genus level, the gene showed one of the lowest levels of concordance with the core genome phylogeny (50.7% average). Concordance for hypervariable regions was lower still, with entropy masking providing little to no benefit. A major factor influencing concordance was SNP count, which showed a positive logarithmic association. Using this relationship, we determined that 690 ± 110 SNPs were required for 80% concordance (average 16S rRNA gene SNP count was 254). We also found a wide range in 16S-23S-5S rRNA operon copy number among genomes (1-27). At the inter-genus level, concordance for the whole 16S rRNA gene was markedly higher (73.8% - 10th out of 49 loci); however, the most concordant hypervariable regions (V4, V3-V4, and V1-V2) ranked in the third quartile (62.5 to 60.0%).

CONCLUSIONS

Ramifications of a poor phylogenetic performance for the 16S rRNA gene are far reaching. For example, in addition to incorrect species/strain delineation and phylogenetic inference, it has the potential to confound community diversity metrics if phylogenetic information is incorporated - for example, with popular approaches such as Faith's phylogenetic diversity and UniFrac. Our results highlight the problematic nature of these approaches and their use (along with entropy masking) is discouraged. Lastly, the wide range in 16S rRNA gene copy number among genomes also has a strong potential to confound diversity metrics. Video Abstract.

Natural Chromosome-Chromid Fusion across rRNA Operons in a Burkholderiaceae Bacterium

Chromids (secondary chromosomes) in bacterial genomes that are present in addition to the main chromosome appear to be evolutionarily conserved in some specific bacterial groups. In rare cases among these groups, a small number of strains from Rhizobiales and Vibrionales were shown to possess a naturally fused single chromosome that was reported to have been generated through intragenomic homologous recombination between repeated sequences on the chromosome and chromid. Similar examples have never been reported in the family Burkholderiaceae, a well-documented group that conserves chromids. Here, an in-depth genomic characterization was performed on a Burkholderiaceae bacterium that was isolated from a soil bacterial consortium maintained on diesel fuel and mutagenic benzo[a]pyrene. This organism, Cupriavidus necator strain KK10, was revealed to carry a single chromosome with unexpectedly large size (>6.6 Mb), and results of comparative genomics with the genome of C. necator N-1T indicated that the single chromosome of KK10 was generated through fusion of the prototypical chromosome and chromid at the rRNA operons. This fusion hypothetically occurred through homologous recombination with a crossover between repeated rRNA operons on the chromosome and chromid. Some metabolic functions that were likely expressed from genes on the prototypical chromid region were indicated to be retained. If this phenomenon-the bacterial chromosome-chromid fusion across the rRNA operons through homologous recombination-occurs universally in prokaryotes, the multiple rRNA operons in bacterial genomes may not only contribute to the robustness of ribosome function, but also provide more opportunities for genomic rearrangements through frequent recombination. IMPORTANCE A bacterial chromosome that was naturally fused with the secondary chromosome, or "chromid," and presented as an unexpectedly large single replicon was discovered in the genome of Cupriavidus necator strain KK10, a biotechnologically useful member of the family Burkholderiaceae. Although Burkholderiaceae is a well-documented group that conserves chromids in their genomes, this chromosomal fusion event has not been previously reported for this family. This fusion has hypothetically occurred through intragenomic homologous recombination between repeated rRNA operons and, if so, provides novel insight into the potential of multiple rRNA operons in bacterial genomes to lead to chromosome-chromid fusion. The harsh conditions under which strain KK10 was maintained-a genotoxic hydrocarbon-enriched milieu-may have provided this genotype with a niche in which to survive.

Exploiting Polyploidy for Markerless and Plasmid-Free Genome Engineering in Cyanobacteria

Here we describe a universal approach for plasmid-free genome engineering in cyanobacteria that exploits the polyploidy of their chromosomes as a natural counterselection system. Rather than being delivered via replicating plasmids, genes encoding for DNA modifying enzymes are instead integrated into essential genes on the chromosome by allelic exchange, as facilitated by antibiotic selection, a process that occurs readily and with only minor fitness defects. By virtue of the essentiality of these integration sites, full segregation is never achieved, with the strain instead remaining as a merodiploid so long as antibiotic selection is maintained. As a result, once the desired genome modification is complete, removal of antibiotic selection results in the gene encoding for the DNA modifying enzyme to then be promptly eliminated from the population. Proof of concept of this new and generalizable strategy is provided using two different site-specific recombination systems, CRE-lox and DRE-rox, in the fast-growing cyanobacterium Synechococcus sp. PCC 7002, as well as CRE-lox in the model cyanobacterium Synechocystis sp. PCC 6803. Reusability of the method, meanwhile, is demonstrated by constructing a high-CO₂ requiring and markerless Δndh3 Δndh4 ΔbicA ΔsbtA mutant of Synechococcus sp. PCC 7002. Overall, this method enables the simple and efficient construction of stable and unmarked mutants in cyanobacteria without the need to develop additional shuttle vectors nor counterselection systems.

High-resolution microbiome analysis enabled by linking of 16S rRNA gene sequences with adjacent genomic contexts

Sequence-based characterization of bacterial communities has long been a hostage of limitations of both 16S rRNA gene and whole metagenome sequencing. Neither approach is universally applicable, and the main efforts to resolve constraints have been devoted to improvement of computational prediction tools. Here, we present semi-targeted 16S rRNA sequencing (st16S-seq), a method designed for sequencing V1-V2 regions of the 16S rRNA gene along with the genomic locus upstream of the gene. By in silico analysis of 13 570 bacterial genome assemblies, we show that genome-linked 16S rRNA sequencing is superior to individual hypervariable regions or full-length gene sequences in terms of classification accuracy and identification of gene copy numbers. Using mock communities and soil samples we experimentally validate st16S-seq and benchmark it against the established microbial classification techniques. We show that st16S-seq delivers accurate estimation of 16S rRNA gene copy numbers, enables taxonomic resolution at the species level and closely approximates community structures obtainable by whole metagenome sequencing.

Bioinformatic pipelines combining denoising and clustering tools allow for more comprehensive prokaryotic and eukaryotic metabarcoding

Environmental DNA metabarcoding is a powerful tool for studying biodiversity. However, bioinformatic approaches need to adjust to the diversity of taxonomic compartments targeted as well as to each barcode gene specificities. We built and tested a pipeline based on read correction with DADA2 allowing analysing metabarcoding data from prokaryotic (16S) and eukaryotic (18S, COI) life compartments. We implemented the option to cluster amplicon sequence variants (ASVs) into operational taxonomic units (OTUs) with swarm, a network-based clustering algorithm, and the option to curate ASVs/OTUs using LULU. Finally, taxonomic assignment was implemented via the Ribosomal Database Project Bayesian classifier (RDP) and BLAST. We validated this pipeline with ribosomal and mitochondrial markers using metazoan mock communities and 42 deep-sea sediment samples. The results show that ASVs and OTUs describe different levels of biotic diversity, the choice of which depends on the research questions. They underline the advantages and complementarity of clustering and LULU-curation for producing metazoan biodiversity inventories at a level approaching the one obtained using morphological criteria. While clustering removes intraspecific variation, LULU effectively removes spurious clusters, originating from errors or intragenomic variability. Swarm clustering affected alpha and beta diversity differently depending on genetic marker. Specifically, d-values > 1 appeared to be less appropriate with 18S for metazoans. Similarly, increasing LULU's minimum ratio level proved essential to avoid losing species in sample-poor data sets. Comparing BLAST and RDP underlined that accurate assignments of deep-sea species can be obtained with RDP, but highlighted the need for a concerted effort to build comprehensive, ecosystem-specific databases.

Rates of gene conversions between Escherichia coli ribosomal operons

Due to their universal presence and high sequence conservation, ribosomal RNA (rRNA) sequences are used widely in phylogenetics for inferring evolutionary relationships between microbes and in metagenomics for analyzing the composition of microbial communities. Most microbial genomes encode multiple copies of rRNA genes to supply cells with sufficient capacity for protein synthesis. These copies typically undergo concerted evolution that keeps their sequences identical, or nearly so, due to gene conversion, a type of intragenomic recombination that changes one copy of a homologous sequence to exactly match another. Widely varying rates of rRNA gene conversion have previously been estimated by comparative genomics methods and using genetic reporter assays. To more directly measure rates of rRNA intragenomic recombination, we sequenced the seven Escherichia coli rRNA operons in 15 lineages that were evolved for ∼13,750 generations with frequent single-cell bottlenecks that reduce the effects of selection. We identified 38 gene conversion events and estimated an overall rate of intragenomic recombination within the 16S and 23S genes between rRNA copies of 3.6 × 10⁻⁴ per genome per generation or 8.6 × 10⁻⁶ per rRNA operon per homologous donor operon per generation. This rate varied only slightly from random expectations at different sites within the rRNA genes and between rRNA operons located at different positions in the genome. Our accurate estimate of the rate of rRNA gene conversions fills a gap in our quantitative understanding of how ribosomal sequences and other multicopy elements diversify and homogenize during microbial genome evolution.

CRISPR with a Happy Ending: Non-Templated DNA Repair for Prokaryotic Genome Engineering

The exploration of microbial metabolism is expected to support the development of a sustainable economy and tackle several problems related to the burdens of human consumption. Microorganisms have the potential to catalyze processes that are currently unavailable, unsustainable and/or inefficient. Their metabolism can be optimized and further expanded using tools like the clustered regularly interspaced short palindromic repeats and their associated proteins (CRISPR-Cas) systems. These tools have revolutionized the field of biotechnology, as they greatly streamline the genetic engineering of organisms from all domains of life. CRISPR-Cas and other nucleases mediate double-strand DNA breaks, which must be repaired to prevent cell death. In prokaryotes, these breaks can be repaired through either homologous recombination, when a DNA repair template is available, or through template-independent end joining, of which two major pathways are known. These end joining pathways depend on different sets of proteins and mediate DNA repair with different outcomes. Understanding these DNA repair pathways can be advantageous to steer the results of genome engineering experiments. In this review, we discuss different strategies for the genetic engineering of prokaryotes through either non-homologous end joining (NHEJ) or alternative end joining (AEJ), both of which are independent of exogenous DNA repair templates.

Endogenous rRNA Sequence Variation Can Regulate Stress Response Gene Expression and Phenotype

Prevailing dogma holds that ribosomes are uniform in composition and function. Here, we show that nutrient limitation-induced stress in E. coli changes the relative expression of rDNA operons to alter the rRNA composition within the actively translating ribosome pool. The most upregulated operon encodes the unique 16S rRNA, rrsH, distinguished by conserved sequence variation within the small ribosomal subunit. rrsH-bearing ribosomes affect the expression of functionally coherent gene sets and alter the levels of the RpoS sigma factor, the master regulator of the general stress response. These impacts are associated with phenotypic changes in antibiotic sensitivity, biofilm formation, and cell motility and are regulated by stress response proteins, RelA and RelE, as well as the metabolic enzyme and virulence-associated protein, AdhE. These findings establish that endogenously encoded, naturally occurring rRNA sequence variation can modulate ribosome function, central aspects of gene expression regulation, and cellular physiology.

Most organisms encode multiple, distinct copies of rRNA genes, rendering the composition of the ribosome pool intrinsically heterogeneous. Here, Kurylo et al. show that nutrient limitation in E. coli upregulates the expression of ribosomes bearing conserved sequence variation in 16S rRNA that can regulate gene expression and phenotype.

IGS sequences in Cestrum present AT- and GC-rich conserved domains, with strong regulatory potential for 5S rDNA

The 35S and 5S ribosomal DNA (rDNA) organized in thousands of copies in genomes, have been widely used in numerous comparative cytogenetic studies. Nevertheless, several questions related to the diversity and organization of regulatory motifs in 5S rDNA remain to be addressed. The 5S rDNA unit is composed of a conserved 120 bp length coding region and an intergenic spacer (IGS) containing potential regulatory motifs (Poly-T, AT-rich and GC-rich) differing in number, redundancy and position along the IGS. The Cestrum species (Solanaceae) have large genomes (about 10 pg/1C) and conserved 2n = 16 karyotypes. Strikingly, these genomes show high diversity of heterochromatin distribution, variability in 35S rDNA loci and the occurrence of B chromosomes. However, the 5S rDNA loci are highly conserved in the proximal region of chromosome 8. Comparison of seventy-one IGS sequences in plants revealed several conserved motifs with potential regulatory function. The AT- and GC-rich domains appeared highly conserved in Cestrum chromosomes. The 5S genic and the GC-rich IGS probe produced FISH signals in both A (pair 8) and B chromosomes. The GC-rich domain presented a strong potential for regulation because it may be associated with CpG islands organization, as well as to hairpin and loop organization. Another interesting aspect was the ability of AT- and GC-rich motifs to produce non-heterochromatic CMA/DAPI signals. While the length of the 5S rDNA IGS region varied in size between the Cestrum species, the individual sequence motifs seem to be conserved suggesting their regulatory function. The most striking feature was the conserved GC-rich domain in Cestrum, which is recognized as a signature trait of the proximal region of chromosome pair 8.

An update of the suicide plasmid-mediated genome editing system in Corynebacterium glutamicum

Corynebacterium glutamicum is an important industrial microorganism, but the availability of tools for its genetic modification has lagged compared to other model microorganisms such as Escherichia coli. Despite great progress in CRISPR-based technologies, the most feasible genome editing method in C. glutamicum is suicide plasmid-mediated, the editing efficiency of which is low due to high false-positive rates of sacB counter selection, and the requirement for tedious two-round selection and verification of rare double-cross-over events. In this study, an rpsL mutant conferring streptomycin resistance was harnessed for counter selection, significantly increasing the positive selection rate. More importantly, with the aid of high selection efficiencies through the use of antibiotics, namely kanamycin and streptomycin, the two-step verification strategy can be simplified to just one-step verification of the final edited strain. As proof of concept, a 2.5-kb DNA fragment comprising aroGfbr pheAfbr expressing cassettes was integrated into the genome of C. glutamicum, with an efficiency of 20% out of the theoretical 50%. The resulting strain produced 110 mg l-1  l-tyrosine in shake-flask fermentation. This updated suicide plasmid-mediated genome editing system will greatly facilitate genetic manipulations including single nucleotide mutation, gene deletion and gene insertion in C. glutamicum and can be easily applied to other microbes.

Molecular identification and phylogenetic relationships of clinical Nocardia isolates

Species identification of Nocardia is difficult because of a complex and rapidly evolving taxonomy. In this study, gene sequencing (16S rRNA, gyrB, secA1, hsp65, rpoB), single 16S rRNA gene sequence phylogenetic analysis, and MALDI-TOF analysis were used to accurately identify 46 clinical Nocardia isolates to the species level. This provided a basis for establishing a routine method of multi-locus sequence analysis (MLSA) for molecular identification of Nocardia species. Genetic polymorphism analysis showed that MLSA was a powerfully discriminating method compared with the 16S rRNA gene to identify clinical Nocardia isolates. However, five-locus (gyrB-16S rRNA-secA1-hsp65-rpoB) MLSA led to misidentifications of all of the five Nocardia abscessus, which were confirmed by digital DNA-DNA hybridization (DDH) analysis. Interestingly, four strains identified as Nocardia beijingensis by a 16S rRNA gene phylogenetic tree may be novel species as suggested by DDH studies. For the purpose of achieving both accuracy and discrimination, the data of MLSA were reanalyzed. A three-locus MLSA with concatenated gyrB-16S rRNA-secA1 sequences was used to construct the phylogenetic tree with high accuracy and powerful discrimination. Therefore, a routine method of MLSA was developed to identify clinical Nocardia species.

Multiple Ribosomal RNA Operons in Bacteria; Their Concerted Evolution and Potential Consequences on the Rate of Evolution of Their 16S rRNA

Bacterial species differ greatly in the number and location of the rRNA operons which may be present in the bacterial chromosomes and plasmids. Most bacterial species contain more than one ribosomal RNA operon copy in their genomes, with some species containing up to 15 such copies. We review the number and location of the rRNA operons and discuss evolution of 16S rRNA (rrs) genes -which are considered as ultimate chronometers for phylogenetic classification- in bacteria with multiple copies of these genes. In these bacterial species, the rrs genes must evolve in concert and sequence changes generated by mutation or horizontal gene transfer must be either erased or spread to every gene copy to avoid divergence, as it occurs when they are present in different species. Analysis of polymorphic sites in intra-genomic rrs copies identifies putative conversion events and demonstrates that sequence conversion is patchy and occurs in small conversion tracts. Sequence conversion probably arises by a non-reciprocal transfer between two or more copies where one copy contributes only a small contiguous segment of DNA, whereas the other copy contributes the rest of the genome in a fairly well understood molecular process. Because concerted evolution implies that a mutation in any of the rrs copies is either eliminated or transferred to every rrs gene in the genome, this process should slow their evolution rate relative to that of single copy genes. However, available data on the rrs genes in bacterial genomes do not show a clear relationship between their evolution rates and the number of their copies in the genome.

Genetic tools for reliable gene expression and recombineering in Pseudomonas putida

Pseudomonas putida is a promising bacterial host for producing natural products, such as polyketides and nonribosomal peptides. In these types of projects, researchers need a genetic toolbox consisting of plasmids, characterized promoters, and techniques for rapidly editing the genome. Past reports described constitutive promoter libraries, a suite of broad host range plasmids that replicate in P. putida, and genome-editing methods. To augment those tools, we have characterized a set of inducible promoters and discovered that IPTG-inducible promoter systems have poor dynamic range due to overexpression of the LacI repressor. By replacing the promoter driving lacI expression with weaker promoters, we increased the fold induction of an IPTG-inducible promoter in P. putida KT2440 to 80-fold. Upon discovering that gene expression from a plasmid was unpredictable when using a high-copy mutant of the BBR1 origin, we determined the copy numbers of several broad host range origins and found that plasmid copy numbers are significantly higher in P. putida KT2440 than in the synthetic biology workhorse, Escherichia coli. Lastly, we developed a λRed/Cas9 recombineering method in P. putida KT2440 using the genetic tools that we characterized. This method enabled the creation of scarless mutations without the need for performing classic two-step integration and marker removal protocols that depend on selection and counterselection genes. With the method, we generated four scarless deletions, three of which we were unable to create using a previously established genome-editing technique.

Phylogenetic Network Analysis Revealed the Occurrence of Horizontal Gene Transfer of 16S rRNA in the Genus Enterobacter

Horizontal gene transfer (HGT) is a ubiquitous genetic event in bacterial evolution, but it seldom occurs for genes involved in highly complex supramolecules (or biosystems), which consist of many gene products. The ribosome is one such supramolecule, but several bacteria harbor dissimilar and/or chimeric 16S rRNAs in their genomes, suggesting the occurrence of HGT of this gene. However, we know little about whether the genes actually experience HGT and, if so, the frequency of such a transfer. This is primarily because the methods currently employed for phylogenetic analysis (e.g., neighbor-joining, maximum likelihood, and maximum parsimony) of 16S rRNA genes assume point mutation-driven tree-shape evolution as an evolutionary model, which is intrinsically inappropriate to decipher the evolutionary history for genes driven by recombination. To address this issue, we applied a phylogenetic network analysis, which has been used previously for detection of genetic recombination in homologous alleles, to the 16S rRNA gene. We focused on the genus Enterobacter, whose phylogenetic relationships inferred by multi-locus sequence alignment analysis and 16S rRNA sequences are incompatible. All 10 complete genomic sequences were retrieved from the NCBI database, in which 71 16S rRNA genes were included. Neighbor-joining analysis demonstrated that the genes residing in the same genomes clustered, indicating the occurrence of intragenomic recombination. However, as suggested by the low bootstrap values, evolutionary relationships between the clusters were uncertain. We then applied phylogenetic network analysis to representative sequences from each cluster. We found three ancestral 16S rRNA groups; the others were likely created through recursive recombination between the ancestors and chimeric descendants. Despite the large sequence changes caused by the recombination events, the RNA secondary structures were conserved. Successive intergenomic and intragenomic recombination thus shaped the evolution of 16S rRNA genes in the genus Enterobacter.

Comparative Pan-Genome Analysis of Piscirickettsia salmonis Reveals Genomic Divergences within Genogroups

Piscirickettsia salmonis is the etiological agent of salmonid rickettsial septicemia, a disease that seriously affects the salmonid industry. Despite efforts to genomically characterize P. salmonis, functional information on the life cycle, pathogenesis mechanisms, diagnosis, treatment, and control of this fish pathogen remain lacking. To address this knowledge gap, the present study conducted an in silico pan-genome analysis of 19 P. salmonis strains from distinct geographic locations and genogroups. Results revealed an expected open pan-genome of 3,463 genes and a core-genome of 1,732 genes. Two marked genogroups were identified, as confirmed by phylogenetic and phylogenomic relationships to the LF-89 and EM-90 reference strains, as well as by assessments of genomic structures. Different structural configurations were found for the six identified copies of the ribosomal operon in the P. salmonis genome, indicating translocation throughout the genetic material. Chromosomal divergences in genomic localization and quantity of genetic cassettes were also found for the Dot/Icm type IVB secretion system. To determine divergences between core-genomes, additional pan-genome descriptions were compiled for the so-termed LF and EM genogroups. Open pan-genomes composed of 2,924 and 2,778 genes and core-genomes composed of 2,170 and 2,228 genes were respectively found for the LF and EM genogroups. The core-genomes were functionally annotated using the Gene Ontology, KEGG, and Virulence Factor databases, revealing the presence of several shared groups of genes related to basic function of intracellular survival and bacterial pathogenesis. Additionally, the specific pan-genomes for the LF and EM genogroups were defined, resulting in the identification of 148 and 273 exclusive proteins, respectively. Notably, specific virulence factors linked to adherence, colonization, invasion factors, and endotoxins were established. The obtained data suggest that these genes could be directly associated with inter-genogroup differences in pathogenesis and host-pathogen interactions, information that could be useful in designing novel strategies for diagnosing and controlling P. salmonis infection.

Highly divergent 16S rRNA sequences in ribosomal operons of Scytonema hyalinum (Cyanobacteria)

A highly divergent 16S rRNA gene was found in one of the five ribosomal operons present in a species complex currently circumscribed as Scytonema hyalinum (Nostocales, Cyanobacteria) using clone libraries. If 16S rRNA sequence macroheterogeneity among ribosomal operons due to insertions, deletions or truncation is excluded, the sequence heterogeneity observed in S. hyalinum was the highest observed in any prokaryotic species thus far (7.3-9.0%). The secondary structure of the 16S rRNA molecules encoded by the two divergent operons was nearly identical, indicating possible functionality. The 23S rRNA gene was examined for a few strains in this complex, and it was also found to be highly divergent from the gene in Type 2 operons (8.7%), and likewise had nearly identical secondary structure between the Type 1 and Type 2 operons. Furthermore, the 16S-23S ITS showed marked differences consistent between operons among numerous strains. Both operons have promoter sequences that satisfy consensus requirements for functional prokaryotic transcription initiation. Horizontal gene transfer from another unknown heterocytous cyanobacterium is considered the most likely explanation for the origin of this molecule, but does not explain the ultimate origin of this sequence, which is very divergent from all 16S rRNA sequences found thus far in cyanobacteria. The divergent sequence is highly conserved among numerous strains of S. hyalinum, suggesting adaptive advantage and selective constraint of the divergent sequence.

A Robust Framework for Microbial Archaeology

Microbial archaeology is flourishing in the era of high-throughput sequencing, revealing the agents behind devastating historical plagues, identifying the cryptic movements of pathogens in prehistory, and reconstructing the ancestral microbiota of humans. Here, we introduce the fundamental concepts and theoretical framework of the discipline, then discuss applied methodologies for pathogen identification and microbiome characterization from archaeological samples. We give special attention to the process of identifying, validating, and authenticating ancient microbes using high-throughput DNA sequencing data. Finally, we outline standards and precautions to guide future research in the field.

Colocalization of distant chromosomal loci in space in E. coli: a bacterial nucleolus

Gaal et al. examined the relative positions of the ribosomal RNA operons in space. The results show that E. coli bacterial chromosome folding in three dimensions is not dictated entirely by genetic position but rather includes functionally related, genetically distant loci that come into close proximity, with rRNA operons forming a structure reminiscent of the eukaryotic nucleolus.

The spatial organization of DNA within the bacterial nucleoid remains unclear. To investigate chromosome organization in Escherichia coli, we examined the relative positions of the ribosomal RNA (rRNA) operons in space. The seven rRNA operons are nearly identical and separated from each other by as much as 180° on the circular genetic map, a distance of ≥2 million base pairs. By inserting binding sites for fluorescent proteins adjacent to the rRNA operons and then examining their positions pairwise in live cells by epifluorescence microscopy, we found that all but rrnC are in close proximity. Colocalization of the rRNA operons required the rrn P1 promoter region but not the rrn P2 promoter or the rRNA structural genes and occurred with and without active transcription. Non-rRNA operon pairs did not colocalize, and the magnitude of their physical separation generally correlated with that of their genetic separation. Our results show that E. coli bacterial chromosome folding in three dimensions is not dictated entirely by genetic position but rather includes functionally related, genetically distant loci that come into close proximity, with rRNA operons forming a structure reminiscent of the eukaryotic nucleolus.

Comparative analyses of phenotypic methods and 16S rRNA, khe, rpoB genes sequencing for identification of clinical isolates of Klebsiella pneumoniae

The present work aimed to evaluate 16S rRNA, khe and rpoB gene sequencing for the identification of Klebsiella pneumoniae in comparison with phenotypic methods. Fifteen clinical isolates were examined, which were initially identified as K. pneumoniae subsp. pneumoniae using the automated VITEK 32 system in two hospitals in Enshi City, China. Their identity was further supported by conventional phenotypic methods on the basis of morphological and biochemical characteristics. Using Bayesian phylogenetic analyses and haplotypes network reconstruction, 13 isolates were identified as K. pneumoniae, whereas the other two isolates (K19, K24) were classified as Shigella sp. and Enterobacter sp., respectively. Of the three genes, 16S rRNA and khe gene could discriminate the clinical isolates at the genus level, whereas rpoB could discriminate Klebsiella at the species and even subspecies level. Overall, the gene tree based on rpoB is more compatible with the currently accepted classification of Klebsiella than those based on 16S rRNA and khe genes, showing that rpoB can be a powerful tool for identification of K. pneumoniae isolates. Above all, our study challenges the utility of khe as a species-specific marker for identification of K. pneumoniae.

Misguided phylogenetic comparisons using DGGE excised bands may contaminate public sequence databases

Controversy surrounding bacterial phylogenies has become one of the most important challenges for microbial ecology. Comparative analyses with nucleotide databases and phylogenetic reconstruction of the amplified 16S rRNA genes from DGGE (Denaturing Gradient Gel Electrophoresis) excised bands have been used by several researchers for the identification of organisms in complex samples. Here, we individually analyzed DGGE-excised 16S rRNA gene bands from 10 certified bacterial strains of different species, and demonstrated that this kind of approach can deliver erroneous outcomes to researchers, besides causing/emphasizing errors in public databases.

Diversity and Evolution in the Genome of Clostridium difficile

Clostridium difficile infection (CDI) is the leading cause of antimicrobial and health care-associated diarrhea in humans, presenting a significant burden to global health care systems. In the last 2 decades, PCR- and sequence-based techniques, particularly whole-genome sequencing (WGS), have significantly furthered our knowledge of the genetic diversity, evolution, epidemiology, and pathogenicity of this once enigmatic pathogen. C. difficile is taxonomically distinct from many other well-known clostridia, with a diverse population structure comprising hundreds of strain types spread across at least 6 phylogenetic clades. The C. difficile species is defined by a large diverse pangenome with extreme levels of evolutionary plasticity that has been shaped over long time periods by gene flux and recombination, often between divergent lineages. These evolutionary events are in response to environmental and anthropogenic activities and have led to the rapid emergence and worldwide dissemination of virulent clonal lineages. Moreover, genome analysis of large clinically relevant data sets has improved our understanding of CDI outbreaks, transmission, and recurrence. The epidemiology of CDI has changed dramatically over the last 15 years, and CDI may have a foodborne or zoonotic etiology. The WGS era promises to continue to redefine our view of this significant pathogen.

Alterations in rRNA-mRNA interaction during plastid evolution

Translation initiation depends on the recognition of mRNA by a ribosome. For this to occur, prokaryotes primarily use the Shine-Dalgarno (SD) interaction, where the 3'-tail of small subunit rRNA (core motif: 3'CCUCC) forms base pairs with a complementary signal sequence in the 5'-untranslated region of mRNA. Here, we examined what happened to SD interactions during the evolution of a cyanobacterial endosymbiont into modern plastids (including chloroplasts). Our analysis of available complete plastid genome sequences revealed that the majority of plastids retained SD interactions but with varying levels of usage. Parallel losses of SD interactions took place in plastids of Chlorophyta, Euglenophyta, and Chromerida/Apicomplexa lineages, presumably related to their extensive reductive evolution. Interestingly, we discovered that the classical SD interaction (3'CCUCC/5'GGAGG [rRNA/mRNA]) was replaced by an altered SD interaction (3'CCCU/5'GGGA or 3'CUUCC/5'GAAGG) through coordinated changes in the sequences of the core rRNA motif and its paired mRNA signal. These changes in plastids of Chlorophyta and Euglenophyta proceeded through intermediate stages that allowed both the classical and altered SD interactions. This coevolution between the rRNA motif and the mRNA signal demonstrates unexpected plasticity in the translation initiation machinery.

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.

The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses

16S ribosomal RNA currently represents the most important target of study in bacterial ecology. Its use for the description of bacterial diversity is, however, limited by the presence of variable copy numbers in bacterial genomes and sequence variation within closely related taxa or within a genome. Here we use the information from sequenced bacterial genomes to explore the variability of 16S rRNA sequences and copy numbers at various taxonomic levels and apply it to estimate bacterial genome and DNA abundances. In total, 7,081 16S rRNA sequences were in silico extracted from 1,690 available bacterial genomes (1-15 per genome). While there are several phyla containing low 16S rRNA copy numbers, in certain taxa, e.g., the Firmicutes and Gammaproteobacteria, the variation is large. Genome sizes are more conserved at all tested taxonomic levels than 16S rRNA copy numbers. Only a minority of bacterial genomes harbors identical 16S rRNA gene copies, and sequence diversity increases with increasing copy numbers. While certain taxa harbor dissimilar 16S rRNA genes, others contain sequences common to multiple species. Sequence identity clusters (often termed operational taxonomic units) thus provide an imperfect representation of bacterial taxa of a certain phylogenetic rank. We have demonstrated that the information on 16S rRNA copy numbers and genome sizes of genome-sequenced bacteria may be used as an estimate for the closest related taxon in an environmental dataset to calculate alternative estimates of the relative abundance of individual bacterial taxa in environmental samples. Using an example from forest soil, this procedure would increase the abundance estimates of Acidobacteria and decrease these of Firmicutes. Using the currently available information, alternative estimates of bacterial community composition may be obtained in this way if the variation of 16S rRNA copy numbers among bacteria is considered.

Structure of rrn operons in pathogenic non-cultivable treponemes: sequence but not genomic position of intergenic spacers correlates with classification of Treponema pallidum and Treponema paraluiscuniculi strains

This study examined the sequences of the two rRNA (rrn) operons of pathogenic non-cultivable treponemes, comprising 11 strains of T. pallidum ssp. pallidum (TPA), five strains of T. pallidum ssp. pertenue (TPE), two strains of T. pallidum ssp. endemicum (TEN), a simian Fribourg-Blanc strain and a rabbit T. paraluiscuniculi (TPc) strain. PCR was used to determine the type of 16S–23S ribosomal intergenic spacers in the rrn operons from 30 clinical samples belonging to five different genotypes. When compared with the TPA strains, TPc Cuniculi A strain had a 17 bp deletion, and the TPE, TEN and Fribourg-Blanc isolates had a deletion of 33 bp. Other than these deletions, only 17 heterogeneous sites were found within the entire region (excluding the 16S–23S intergenic spacer region encoding tRNA-Ile or tRNA-Ala). The pattern of nucleotide changes in the rrn operons corresponded to the classification of treponemal strains, whilst two different rrn spacer patterns (Ile/Ala and Ala/Ile) appeared to be distributed randomly across species/subspecies classification, time and geographical source of the treponemal strains. It is suggested that the random distribution of tRNA genes is caused by reciprocal translocation between repetitive sequences mediated by a recBCD-like system.

Diversity of 5S rRNA genes within individual prokaryotic genomes

We examined intragenomic variation of paralogous 5S rRNA genes to evaluate the concept of ribosomal constraints. In a dataset containing 1161 genomes from 779 unique species, 96 species exhibited > 3% diversity. Twenty-seven species with > 10% diversity contained a total of 421 mismatches between all pairs of the most dissimilar copies of 5S rRNA genes. The large majority (401 of 421) of the diversified positions were conserved at the secondary structure level. The high diversity was associated with partial rRNA operon, split operon, or spacer length-related divergence. In total, these findings indicated that there are tight ribosomal constraints on paralogous 5S rRNA genes in a genome despite of the high degree of diversity at the primary structure level.

A novel method for simultaneous Enterococcus species identification/typing and van genotyping by high resolution melt analysis

In order to develop a typing and identification method for van gene containing Enterococcus faecium, two multiplex PCR reactions were developed for use in HRM-PCR (High Resolution Melt-PCR): (i) vanA, vanB, vanC, vanC23 to detect van genes from different Enterococcus species; (ii) ISR (intergenic spacer region between the 16S and 23S rRNA genes) to detect all Enterococcus species and obtain species and isolate specific HRM curves. To test and validate the method three groups of isolates were tested: (i) 1672 Enterococcus species isolates from January 2009 to December 2009; (ii) 71 isolates previously identified and typed by PFGE (pulsed-field gel electrophoresis) and MLST (multi-locus sequence typing); and (iii) 18 of the isolates from (i) for which ISR sequencing was done. As well as successfully identifying 2 common genotypes by HRM from the Austin Hospital clinical isolates, this study analysed the sequences of all the vanB genes deposited in GenBank and developed a numerical classification scheme for the standardised naming of these vanB genotypes. The identification of Enterococcus faecalis from E. faecium was reliable and stable using ISR PCR. The typing of E. faecium by ISR PCR: (i) detected two variable peaks corresponding to different copy numbers of insertion sequences I and II corresponding to peak I and II respectively; (ii) produced 7 melt profiles for E. faecium with variable copy numbers of sequences I and II; (iii) demonstrated stability and instability of peak heights with equal frequency within the patient sample (36.4±4.5 days and 38.6±5.8 days respectively for 192 patients); (iv) detected ISR-HRM types with as much discrimination as PFGE and more than MLST; and (v) detected ISR-HRM types that differentiated some isolates that were identical by PFGE and MLST. In conjunction with the rapid and accurate van genotyping method described here, this ISR-HRM typing and identification method can be used as a stable identification and typing method with predictable instability based on recombination and concerted evolution of the rrn operon that will complement existing typing methods.

Variable copy number, intra-genomic heterogeneities and lateral transfers of the 16S rRNA gene in Pseudomonas

Even though the 16S rRNA gene is the most commonly used taxonomic marker in microbial ecology, its poor resolution is still not fully understood at the intra-genus level. In this work, the number of rRNA gene operons, intra-genomic heterogeneities and lateral transfers were investigated at a fine-scale resolution, throughout the Pseudomonas genus. In addition to nineteen sequenced Pseudomonas strains, we determined the 16S rRNA copy number in four other Pseudomonas strains by Southern hybridization and Pulsed-Field Gel Electrophoresis, and studied the intra-genomic heterogeneities by Denaturing Gradient Gel Electrophoresis and sequencing. Although the variable copy number (from four to seven) seems to be correlated with the evolutionary distance, some close strains in the P. fluorescens lineage showed a different number of 16S rRNA genes, whereas all the strains in the P. aeruginosa lineage displayed the same number of genes (four copies). Further study of the intra-genomic heterogeneities revealed that most of the Pseudomonas strains (15 out of 19 strains) had at least two different 16S rRNA alleles. A great difference (5 or 19 nucleotides, essentially grouped near the V1 hypervariable region) was observed only in two sequenced strains. In one of our strains studied (MFY30 strain), we found a difference of 12 nucleotides (grouped in the V3 hypervariable region) between copies of the 16S rRNA gene. Finally, occurrence of partial lateral transfers of the 16S rRNA gene was further investigated in 1803 full-length sequences of Pseudomonas available in the databases. Remarkably, we found that the two most variable regions (the V1 and V3 hypervariable regions) had probably been laterally transferred from another evolutionary distant Pseudomonas strain for at least 48.3 and 41.6% of the 16S rRNA sequences, respectively. In conclusion, we strongly recommend removing these regions of the 16S rRNA gene during the intra-genus diversity studies.

Large variations in bacterial ribosomal RNA genes

Ribosomal RNA (rRNA) genes, essential to all forms of life, have been viewed as highly conserved and evolutionarily stable, partly because very little is known about their natural variations. Here, we explored large-scale variations of rRNA genes through bioinformatic analyses of available complete bacterial genomic sequences with an emphasis on formation mechanisms and biological significance. Interestingly, we found bacterial genomes in which no 16S rRNA genes harbor the conserved core of the anti–Shine-Dalgarno sequence (5′-CCTCC-3′). This loss was accompanied by elimination of Shine-Dalgarno–like sequences upstream of their protein-coding genes. Those genomes belong to 1 or 2 of the following categories: primary symbionts, hemotropic Mycoplasma, and Flavobacteria. We also found many rearranged rRNA genes and reconstructed their history. Conjecturing the underlying mechanisms, such as inversion, partial duplication, transposon insertion, deletion, and substitution, we were able to infer their biological significance, such as co-orientation of rRNA transcription and chromosomal replication, lateral transfer of rRNA gene segments, and spread of rRNA genes with an apparent structural defect through gene conversion. These results open the way to understanding dynamic evolutionary changes of rRNA genes and the translational machinery.

Intragenomic heterogeneity of the 16S rRNA gene in strain UFO1 caused by a 100-bp insertion in helix 6

Two different versions of the 16S rRNA gene, one of which contained an unusual 100-bp insertion in helix 6, were detected in isolate UFO1 acquired from the Oak Ridge Integrated Field-Research Challenge (ORIFRC) site in Tennessee. rRNA was extracted from UFO1 and analyzed by reverse transcriptase-quantitative PCR with insert- and non-insert-specific primers; only the noninsert 16S rRNA gene sequence was detected. Similarly, PCR-based screening of a cDNA library (190 clones) constructed from reverse-transcribed rRNA from UFO1 did not detect any clones containing the 100-bp insert. Examination of cDNA with primers specific to the insert-bearing 16S rRNA gene, but downstream of the insert, suggests that the insert was excised from rRNA. Inspection of other 16S rRNA genes in the GenBank database revealed that a homologous insert sequence, also found in helix 6, has been reported in other environmental clones, including those acquired from ORIFRC enrichments. These findings demonstrate the existence of widely divergent copies of the 16S rRNA gene within the same organism, which may confound 16S rRNA gene-based methods of estimating microbial diversity in environmental samples.

Diversity of 16S rRNA genes within individual prokaryotic genomes

Analysis of intragenomic variation of 16S rRNA genes is a unique approach to examining the concept of ribosomal constraints on rRNA genes; the degree of variation is an important parameter to consider for estimation of the diversity of a complex microbiome in the recently initiated Human Microbiome Project (http://nihroadmap.nih.gov/hmp). The current GenBank database has a collection of 883 prokaryotic genomes representing 568 unique species, of which 425 species contained 2 to 15 copies of 16S rRNA genes per genome (2.22 +/- 0.81). Sequence diversity among the 16S rRNA genes in a genome was found in 235 species (from 0.06% to 20.38%; 0.55% +/- 1.46%). Compared with the 16S rRNA-based threshold for operational definition of species (1 to 1.3% diversity), the diversity was borderline (between 1% and 1.3%) in 10 species and >1.3% in 14 species. The diversified 16S rRNA genes in Haloarcula marismortui (diversity, 5.63%) and Thermoanaerobacter tengcongensis (6.70%) were highly conserved at the 2 degrees structure level, while the diversified gene in B. afzelii (20.38%) appears to be a pseudogene. The diversified genes in the remaining 21 species were also conserved, except for a truncated 16S rRNA gene in "Candidatus Protochlamydia amoebophila." Thus, this survey of intragenomic diversity of 16S rRNA genes provides strong evidence supporting the theory of ribosomal constraint. Taxonomic classification using the 16S rRNA-based operational threshold could misclassify a number of species into more than one species, leading to an overestimation of the diversity of a complex microbiome. This phenomenon is especially seen in 7 bacterial species associated with the human microbiome or diseases.

Homozygous triplicate mutations in three 16S rRNA genes responsible for high-level aminoglycoside resistance in Nocardia farcinica clinical isolates from a Canada-wide bovine mastitis epizootic

Nocardia farcinica strains showing high-level resistance to amikacin were isolated from clinical cases in a Canada-wide bovine mastitis epizootic. Shotgun cloning of the resistance genes in the amikacin-resistant mastitis isolate N. farcinica IFM 10580 (W6220 [Centers for Disease Control and Prevention]) using a multicopy vector system revealed that the 16S rRNA gene with an A-to-G single-point mutation at position 1408 (in Escherichia coli numbering) conferred "moderate" cross-resistance to amikacin and other aminoglycosides to an originally susceptible N. farcinica strain IFM 10152. Subsequent DNA sequence analyses revealed that, in contrast to the susceptible strain, all three chromosomal 16S rRNA genes of IFM 10580, the epizootic clinical strain, contained the same A1408G point mutations. Mutant colonies showing high-level aminoglycoside resistance were obtained when the susceptible strain N. farcinica IFM 10152 was transformed with a multicopy plasmid carrying the A1408G mutant 16S rRNA gene and was cultured in the presence of aminoglycosides for 3 to 5 days. Of these transformants, at least two of the three chromosomal 16S rRNA genes contained A1408G mutations. A triple mutant was easily obtained from a strain carrying the two chromosomal A1408G mutant genes and one wild-type gene, even in the absence of the plasmid. The triple mutant showed the highest level of resistance to aminoglycosides, even in the absence of the plasmid carrying the mutant 16S rRNA gene. These results suggest that the homozygous mutations in the three 16S rRNA genes are responsible for the high-level aminoglycoside resistance found in N. farcinica isolates of the bovine mastitis epizootic.

Reverse line blot hybridization and DNA sequencing studies of the 16S-23S rRNA gene intergenic spacer regions of five emerging pathogenic Nocardia species

The objective of this study was to examine DNA sequence polymorphisms in the 16S-23S rRNA gene intergenic spacer (ITS) regions of five emerging pathogenic Nocardia species: Nocardia beijingensis, Nocardia blacklockiae, Nocardia thailandica, Nocardia elegans and Nocardia vinacea. A set of six isolates belonging to the species of interest and 135 isolates belonging to other Nocardia species was studied. A PCR-based reverse line blot (RLB) hybridization assay incorporating species- or intraspecies ITS rRNA gene operon-specific probes was then developed for species identification. Substantial intraspecies sequence variation among different ITS operons was identified. Four sequence types of N. thailandica, eight sequence types of N. beijingensis (four types for each of two strains) and five sequence types of N. blacklockiae, N. elegans and N. vinacea were found. The results represent the first evidence of ITS sequence heterogeneity in emerging species of Nocardia. By incorporating species/operon-specific probes into a RLB assay, unique RLB patterns were identified for each of the species and every sequence type. The PCR/RLB assay demonstrated high specificity and showed promise in both the identification and genotyping of Nocardia species. More detailed studies of the polymorphism within the ITS locus may further advance our capacity to reliably identify and subtype medically important Nocardia species.

Utilization of an unstable plasmid and the I-SceI endonuclease to generate routine markerless deletion mutants in Francisella tularensis

We engineered an efficient system to make Francisella tularensis deletion mutations using an unstable, poorly maintained plasmid to enhance the likelihood of homologous recombination. For counterselection, we adapted a strategy using I-SceI, which causes a double-stranded break in the integrated suicide vector, forcing a second recombination to mediate allelic replacement.

Understanding the differences between genome sequences of Escherichia coli B strains REL606 and BL21(DE3) and comparison of the E. coli B and K-12 genomes

Each difference between the genome sequences of Escherichia coli B strains REL606 and BL21(DE3) can be interpreted in light of known laboratory manipulations plus a gene conversion between ribosomal RNA operons. Two treatments with 1-methyl-3-nitro-1-nitrosoguanidine in the REL606 lineage produced at least 93 single-base-pair mutations ( approximately 90% GC-to-AT transitions) and 3 single-base-pair GC deletions. Two UV treatments in the BL21(DE3) lineage produced only 4 single-base-pair mutations but 16 large deletions. P1 transductions from K-12 into the two B lineages produced 317 single-base-pair differences and 9 insertions or deletions, reflecting differences between B DNA in BL21(DE3) and integrated restriction fragments of K-12 DNA inherited by REL606. Two sites showed selective enrichment of spontaneous mutations. No unselected spontaneous single-base-pair mutations were evident. The genome sequences revealed that a progenitor of REL606 had been misidentified, explaining initially perplexing differences. Limited sequencing of other B strains defined characteristic properties of B and allowed assembly of the inferred genome of the ancestral B of Delbrück and Luria. Comparison of the B and K-12 genomes shows that more than half of the 3793 proteins of their basic genomes are predicted to be identical, although approximately 310 appear to be functional in either B or K-12 but not in both. The ancestral basic genome appears to have had approximately 4039 coding sequences occupying approximately 4.0 Mbp. Repeated horizontal transfer from diverged Escherichia coli genomes and homologous recombination may explain the observed variable distribution of single-base-pair differences. Fifteen sites are occupied by phage-related elements, but only six by comparable elements at the same site. More than 50 sites are occupied by IS elements in both B and K, 16 in common, and likely founding IS elements are identified. A signature of widespread cryptic phage P4-type mobile elements was identified. Complex deletions (dense clusters of small deletions and substitutions) apparently removed nonessential genes from approximately 30 sites in the basic genomes.

Divergence among genes encoding the elongation factor Tu of Yersinia Species

Elongation factor Tu (EF-Tu), encoded by tuf genes, carries aminoacyl-tRNA to the ribosome during protein synthesis. Duplicated tuf genes (tufA and tufB), which are commonly found in enterobacterial species, usually coevolve via gene conversion and are very similar to one another. However, sequence analysis of tuf genes in our laboratory has revealed highly divergent copies in 72 strains spanning the genus Yersinia (representing 12 Yersinia species). The levels of intragenomic divergence between tufA and tufB sequences ranged from 8.3 to 16.2% for the genus Yersinia, which is significantly greater than the 0.0 to 3.6% divergence observed for other enterobacterial genera. We further explored tuf gene evolution in Yersinia and other Enterobacteriaceae by performing directed sequencing and phylogenetic analyses. Phylogenetic trees constructed using concatenated tufA and tufB sequences revealed a monophyletic genus Yersinia in the family Enterobacteriaceae. Moreover, Yersinia strains form clades within the genus that mostly correlate with their phenotypic and genetic classifications. These genetic analyses revealed an unusual divergence between Yersinia tufA and tufB sequences, a feature unique among sequenced Enterobacteriaceae and indicative of a genus-wide loss of gene conversion. Furthermore, they provided valuable phylogenetic information for possible reclassification and identification of Yersinia species.

rrnDB: documenting the number of rRNA and tRNA genes in bacteria and archaea

A dramatic exception to the general pattern of single-copy genes in bacterial and archaeal genomes is the presence of 1–15 copies of each ribosomal RNA encoding gene. The original version of the Ribosomal RNA Database (rrnDB) cataloged estimates of the number of 16S rRNA-encoding genes; the database now includes the number of genes encoding each of the rRNAs (5S, 16S and 23S), an internally transcribed spacer region, and the number of tRNA genes. The rrnDB has been used largely by microbiologists to predict the relative rate at which microbial populations respond to favorable growth conditions, and to interpret 16S rRNA-based surveys of microbial communities. To expand the functionality of the rrnDB (http://ribosome.mmg.msu.edu/rrndb/index.php), the search engine has been redesigned to allow database searches based on 16S rRNA gene copy number, specific organisms or taxonomic subsets of organisms. The revamped database also computes average gene copy numbers for any collection of entries selected. Curation tools now permit rapid updates, resulting in an expansion of the database to include data for 785 bacterial and 69 archaeal strains. The rrnDB continues to serve as the authoritative, curated source that documents the phylogenetic distribution of rRNA and tRNA genes in microbial genomes.

Ectopic gene conversions in bacterial genomes

We characterized the gene conversions found between the duplicated genes of 75 bacterial genomes from five species groups (archaea, nonpathogenic and pathogenic firmicutes, and nonpathogenic and pathogenic proteobacteria). The number of gene conversions is positively correlated with the size of multigene families and the size of multigene families is not significantly different between pathogenic and nonpathogenic taxa. However, gene conversions occur twice as frequently in pathogenic species as in nonpathogenic species. Comparisons between closely related species also indicate a trend towards increased gene conversion in pathogenic species. Whereas the length of the conversions is positively correlated with flanking sequence similarity in all five groups, these correlations are smaller for pathogenic firmicutes and proteobacteria than for nonpathogenic firmicutes and proteobacteria. These results are consistent with our previous work on E. coli genomes and suggest that pathogenic bacteria allow recombination between more divergent gene sequences. This higher permissiveness is likely adaptive because it allows them to generate more genetic variability.

Evolution of Streptococcus pneumoniae and its close commensal relatives

Streptococcus pneumoniae is a member of the Mitis group of streptococci which, according to 16S rRNA-sequence based phylogenetic reconstruction, includes 12 species. While other species of this group are considered prototypes of commensal bacteria, S. pneumoniae is among the most frequent microbial killers worldwide. Population genetic analysis of 118 strains, supported by demonstration of a distinct cell wall carbohydrate structure and competence pheromone sequence signature, shows that S. pneumoniae is one of several hundred evolutionary lineages forming a cluster separate from Streptococcus oralis and Streptococcus infantis. The remaining lineages of this distinct cluster are commensals previously collectively referred to as Streptococcus mitis and each represent separate species by traditional taxonomic standard. Virulence genes including the operon for capsule polysaccharide synthesis and genes encoding IgA1 protease, pneumolysin, and autolysin were randomly distributed among S. mitis lineages. Estimates of the evolutionary age of the lineages, the identical location of remnants of virulence genes in the genomes of commensal strains, the pattern of genome reductions, and the proportion of unique genes and their origin support the model that the entire cluster of S. pneumoniae, S. pseudopneumoniae, and S. mitis lineages evolved from pneumococcus-like bacteria presumably pathogenic to the common immediate ancestor of hominoids. During their adaptation to a commensal life style, most of the lineages gradually lost the majority of genes determining virulence and became genetically distinct due to sexual isolation in their respective hosts.

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.

Mesorhizobium albiziae sp. nov., a novel bacterium that nodulates Albizia kalkora in a subtropical region of China

A novel Mesorhizobium group associated with Albizia kalkora [Wang et al. (2006), Syst Appl Microbiol 29, 502-517] was further characterized. The seven strains in this group showed similar protein patterns and were different from defined Mesorhizobium species in SDS-PAGE of whole-cell proteins. The representative strain CCBAU 61158(T) formed a novel Mesorhizobium lineage in phylogenetic analyses of 16S rRNA, atpD, glnII and nifH genes. However, its nodC gene sequence was more similar to that of Rhizobium gallicum R602sp(T) than to those of Mesorhizobium species. DNA-DNA relatedness between CCBAU 61158(T) and reference strains of defined Mesorhizobium species was lower than 34.1 %. These results indicated that this Mesorhizobium group was a unique genomic species. The subtropical distribution, host origin, PCR-RFLP patterns of 16S rRNA genes, fatty acid profile and a series of phenotypic characteristics could be used as distinctive features of this group. This group is therefore proposed as a novel species, Mesorhizobium albiziae sp. nov., with CCBAU 61158(T) (=LMG 23507(T)=USDA 4964(T)) as the type strain. Strain CCBAU 61158(T) could form effective nodules on Albizia julibrissin, Glycine max, Leucaena leucocephala and Phaseolus vulgaris.