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

Spore-forming bacteria in gelatin: Characterization, identification by 16S rRNA and MALDI-TOF mass spectrometry (MS), and presence of heat resistance and virulence genes

Gelatin, a versatile protein derived from collagen, is widely used in the food, pharmaceutical and medical sectors. However, bacterial contamination by spore-forming bacteria during gelatin processing represents a significant concern for product safety and quality. In this study, an investigation was carried out to explore the heat and chemical resistance, as well as the identification and characterization of spore-forming bacteria isolated from gelatin processing. The methodologies involved chemical resistance tests with drastic pH in microplates and thermal resistance tests in capillary tubes of various isolates obtained at different processing stages. In addition, phenotypic and genotypic analyses were carried out to characterize the most resistant isolates of spore-forming bacteria. The findings of this study revealed the presence of several species, including Bacillus cereus, Bacillus licheniformis, Bacillus sonorensis, Bacillus subtilis, Geobacillus stearothermophilus, and Clostridium sporogenes, with some isolates exhibiting remarkable chemical and heat resistances. In addition, a significant proportion of the most resistant isolates showed gelatinase activity (n = 19/21; 90.5 %) and the presence of heat resistance (n = 5/21; 23.8 %), and virulence genes (n = 11/21; 52.4 %). The results of this study suggest that interventions should be done in quality control practices and that process parameter adjustments and effective contamination reduction strategies should be implemented through gelatin processing.

Intermolecular Gene Conversion for the Equalization of Genome Copies in the Polyploid Haloarchaeon Haloferax volcanii: Identification of Important Proteins

The model haloarchaeon Haloferax volcanii is polyploid with about 20 copies of its major chromosome. Recently it has been described that highly efficient intermolecular gene conversion operates in H. volcanii to equalize the chromosomal copies. In the current study, 24 genes were selected that encode proteins with orthologs involved in gene conversion or homologous recombination in archaea, bacteria, or eukaryotes. Single gene deletion strains of 22 genes and a control gene were constructed in two parent strains for a gene conversion assay; only radA and radB were shown to be essential. Protoplast fusions were used to generate strains that were heterozygous for the gene HVO_2528, encoding an enzyme for carotinoid biosynthesis. It was revealed that a lack of six of the proteins did not influence the efficiency of gene conversion, while sixteen mutants had severe gene conversion defects. Notably, lack of paralogous proteins of gene families had very different effects, e.g., mutant Δrad25b had no phenotype, while mutants Δrad25a, Δrad25c, and Δrad25d were highly compromised. Generation of a quadruple rad25 and a triple sph deletion strain also indicated that the paralogs have different functions, in contrast to sph2 and sph4, which cannot be deleted simultaneously. There was no correlation between the severity of the phenotypes and the respective transcript levels under non-stressed conditions, indicating that gene expression has to be induced at the onset of gene conversion. Phylogenetic trees of the protein families Rad3/25, MutL/S, and Sph/SMC/Rad50 were generated to unravel the history of the paralogous proteins of H. volcanii. Taken together, unselected intermolecular gene conversion in H. volcanii involves at least 16 different proteins, the molecular roles of which can be studied in detail in future projects.

High-throughput single-cell transcriptomics of bacteria using combinatorial barcoding

Microbial split-pool ligation transcriptomics (microSPLiT) is a high-throughput single-cell RNA sequencing method for bacteria. With four combinatorial barcoding rounds, microSPLiT can profile transcriptional states in hundreds of thousands of Gram-negative and Gram-positive bacteria in a single experiment without specialized equipment. As bacterial samples are fixed and permeabilized before barcoding, they can be collected and stored ahead of time. During the first barcoding round, the fixed and permeabilized bacteria are distributed into a 96-well plate, where their transcripts are reverse transcribed into cDNA and labeled with the first well-specific barcode inside the cells. The cells are mixed and redistributed two more times into new 96-well plates, where the second and third barcodes are appended to the cDNA via in-cell ligation reactions. Finally, the cells are mixed and divided into aliquot sub-libraries, which can be stored until future use or prepared for sequencing with the addition of a fourth barcode. It takes 4 days to generate sequencing-ready libraries, including 1 day for collection and overnight fixation of samples. The standard plate setup enables single-cell transcriptional profiling of up to 1 million bacterial cells and up to 96 samples in a single barcoding experiment, with the possibility of expansion by adding barcoding rounds. The protocol requires experience in basic molecular biology techniques, handling of bacterial samples and preparation of DNA libraries for next-generation sequencing. It can be performed by experienced undergraduate or graduate students. Data analysis requires access to computing resources, familiarity with Unix command line and basic experience with Python or R.

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

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

The overlooked evolutionary dynamics of 16S rRNA revises its role as the "gold standard" for bacterial species identification

The role of 16S rRNA has been and largely remains crucial for the identification of microbial organisms. Although 16S rRNA could certainly be described as one of the most studied sequences ever, the current view of it remains somewhat ambiguous. While some consider 16S rRNA to be a variable marker with resolution power down to the strain level, others consider them to be living fossils that carry information about the origin of domains of cellular life. We show that 16S rRNA is clearly an evolutionarily very rigid sequence, making it a largely unique and irreplaceable marker, but its applicability beyond the genus level is highly limited. Interestingly, it seems that the evolutionary rigidity is not driven by functional constraints of the sequence (RNA-protein interactions), but rather results from the characteristics of the host organism. Our results suggest that, at least in some lineages, Horizontal Gene Transfer (HGT) within genera plays an important role for the evolutionary non-dynamics (stasis) of 16S rRNA. Such genera exhibit an apparent lack of diversification at the 16S rRNA level in comparison to the rest of a genome. However, why it is limited specifically and solely to 16S rRNA remains enigmatic.

Genetic disruption of the bacterial raiA motif noncoding RNA causes defects in sporulation and aggregation

Several structured noncoding RNAs in bacteria are essential contributors to fundamental cellular processes. Thus, discoveries of additional ncRNA classes provide opportunities to uncover and explore biochemical mechanisms relevant to other major and potentially ancient processes. A candidate structured ncRNA named the "raiA motif" has been found via bioinformatic analyses in over 2,500 bacterial species. The gene coding for the RNA typically resides between the raiA and comFC genes of many species of Bacillota and Actinomycetota. Structural probing of the raiA motif RNA from the Gram-positive anaerobe Clostridium acetobutylicum confirms key features of its sophisticated secondary structure model. Expression analysis of raiA motif RNA reveals that the RNA is constitutively produced but reaches peak abundance during the transition from exponential growth to stationary phase. The raiA motif RNA becomes the fourth most abundant RNA in C. acetobutylicum, excluding ribosomal RNAs and transfer RNAs. Genetic disruption of the raiA motif RNA causes cells to exhibit substantially decreased spore formation and diminished ability to aggregate. Restoration of normal cellular function in this knock-out strain is achieved by expression of a raiA motif gene from a plasmid. These results demonstrate that raiA motif RNAs normally participate in major cell differentiation processes by operating as a trans-acting factor.

Microbial diversity and abundance vary along salinity, oxygen, and particle size gradients in the Chesapeake Bay

Marine snow and other particles are abundant in estuaries, where they drive biogeochemical transformations and elemental transport. Particles range in size, thereby providing a corresponding gradient of habitats for marine microorganisms. We used standard normalized amplicon sequencing, verified with microscopy, to characterize taxon-specific microbial abundances, (cells per litre of water and per milligrams of particles), across six particle size classes, ranging from 0.2 to 500 μm, along the main stem of the Chesapeake Bay estuary. Microbial communities varied in salinity, oxygen concentrations, and particle size. Many taxonomic groups were most densely packed on large particles (in cells/mg particles), yet were primarily associated with the smallest particle size class, because small particles made up a substantially larger portion of total particle mass. However, organisms potentially involved in methanotrophy, nitrite oxidation, and sulphate reduction were found primarily on intermediately sized (5-180 μm) particles, where species richness was also highest. All abundant ostensibly free-living organisms, including SAR11 and Synecococcus, appeared on particles, albeit at lower abundance than in the free-living fraction, suggesting that aggregation processes may incorporate them into particles. Our approach opens the door to a more quantitative understanding of the microscale and macroscale biogeography of marine microorganisms.

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

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

Harnessing Variabilities in Digital Melt Curves for Accurate Identification of Bacteria

Digital PCR combined with high resolution melt (HRM) is an emerging method for identifying pathogenic bacteria with single cell resolution via species-specific digital melt curves. Currently, the development of such digital PCR-HRM assays entails first identifying PCR primers to target hypervariable gene regions within the target bacteria panel, next performing bulk-based PCR-HRM to examine whether the resulting species-specific melt curves possess sufficient interspecies variability (i.e., variability between bacterial species), and then digitizing the bulk-based PCR-HRM assays with melt curves that have high interspecies variability via microfluidics. In this work, we first report our discovery that the current development workflow can be inadequate because a bulk-based PCR-HRM assay that produces melt curves with high interspecies variability can, in fact, lead to a digital PCR-HRM assay that produces digital melt curves with unwanted intraspecies variability (i.e., variability within the same bacterial species), consequently hampering bacteria identification accuracy. Our subsequent investigation reveals that such intraspecies variability in digital melt curves can arise from PCR primers that target nonidentical gene copies or amplify nonspecifically. We then show that computational in silico HRM opens a window to inspect both interspecies and intraspecies variabilities and thus provides the missing link between bulk-based PCR-HRM and digital PCR-HRM. Through this new development workflow, we report a new digital PCR-HRM assay with improved bacteria identification accuracy. More broadly, this work can serve as the foundation for enhancing the development of future digital PCR-HRM assays toward identifying causative pathogens and combating infectious diseases.

Ribosomal proteins can hold a more accurate record of bacterial thermal adaptation compared to rRNA

Ribosomal genes are widely used as 'molecular clocks' to infer evolutionary relationships between species. However, their utility as 'molecular thermometers' for estimating optimal growth temperature of microorganisms remains uncertain. Previously, some estimations were made using the nucleotide composition of ribosomal RNA (rRNA), but the universal application of this approach was hindered by numerous outliers. In this study, we aimed to address this problem by identifying additional indicators of thermal adaptation within the sequences of ribosomal proteins. By comparing sequences from 2021 bacteria with known optimal growth temperature, we identified novel indicators among the metal-binding residues of ribosomal proteins. We found that these residues serve as conserved adaptive features for bacteria thriving above 40°C, but not at lower temperatures. Furthermore, the presence of these metal-binding residues exhibited a stronger correlation with the optimal growth temperature of bacteria compared to the commonly used correlation with the 16S rRNA GC content. And an even more accurate correlation was observed between the optimal growth temperature and the YVIWREL amino acid content within ribosomal proteins. Overall, our work suggests that ribosomal proteins contain a more accurate record of bacterial thermal adaptation compared to rRNA. This finding may simplify the analysis of unculturable and extinct species.

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

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

Macroevolutionary constraints on global microbial diversity

Biologists have long sought to quantify the number of species on Earth. Often missing from these efforts is the contribution of microorganisms, the smallest but most abundant form of life on the planet. Despite recent large‐scale sampling efforts, estimates of global microbial diversity span many orders of magnitude. It is important to consider how speciation and extinction over the last 4 billion years constrain inventories of biodiversity. We parameterized macroevolutionary models based on birth–death processes that assume constant and universal speciation and extinction rates. The models reveal that richness beyond 1012 species is feasible and in agreement with empirical predictions. Additional simulations suggest that mass extinction events do not place hard limits on modern‐day microbial diversity. Together, our study provides independent support for a massive global‐scale microbiome while shedding light on the upper limits of life on Earth.

Ensifer canadensis sp. nov. strain T173T isolated from Melilotus albus (sweet clover) in Canada possesses recombinant plasmid pT173b harbouring symbiosis and type IV secretion system genes apparently acquired from Ensifer medicae

A bacterial strain, designated T173^T, was previously isolated from a root-nodule of a Melilotus albus plant growing in Canada and identified as a novel Ensifer lineage that shared a clade with the non-symbiotic species, Ensifer adhaerens. Strain T173^T was also previously found to harbour a symbiosis plasmid and to elicit root-nodules on Medicago and Melilotus species but not fix nitrogen. Here we present data for the genomic and taxonomic description of strain T173^T. Phylogenetic analyses including the analysis of whole genome sequences and multiple locus sequence analysis (MLSA) of 53 concatenated ribosome protein subunit (rps) gene sequences confirmed placement of strain T173^T in a highly supported lineage distinct from named Ensifer species with E. morelensis Lc04^T as the closest relative. The highest digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values of genome sequences of strain T173^T compared with closest relatives (35.7 and 87.9%, respectively) are well below the respective threshold values of 70% and 95-96% for bacterial species circumscription. The genome of strain T173^T has a size of 8,094,229 bp with a DNA G + C content of 61.0 mol%. Six replicons were detected: a chromosome (4,051,102 bp) and five plasmids harbouring plasmid replication and segregation (repABC) genes. These plasmids were also found to possess five apparent conjugation systems based on analysis of TraA (relaxase), TrbE/VirB4 (part of the Type IV secretion system (T4SS)) and TraG/VirD4 (coupling protein). Ribosomal RNA operons encoding 16S, 23S, and 5S rRNAs that are usually restricted to bacterial chromosomes were detected on plasmids pT173d and pT173e (946,878 and 1,913,930 bp, respectively) as well as on the chromosome of strain T173^T. Moreover, plasmid pT173b (204,278 bp) was found to harbour T4SS and symbiosis genes, including nodulation (nod, noe, nol) and nitrogen fixation (nif, fix) genes that were apparently acquired from E. medicae by horizontal transfer. Data for morphological, physiological and symbiotic characteristics complement the sequence-based characterization of strain T173^T. The data presented support the description of a new species for which the name Ensifer canadensis sp. nov. is proposed with strain T173^T (= LMG 32374^T = HAMBI 3766^T) as the species type strain.

Seasonal and spatial variations of antibiotic resistance genes and bacterial biodiversity in biofilms covering the equipment at successive stages of drinking water purification

The presence of ARGs (antibiotic resistance genes) in the aquatic environment is a serious threat to public health especially in environmental biofilms as natural reservoirs of ARGs in water treatment plants (WTP). It has been shown that the treatment technology and source of water have a significant impact on the abundance and type of genes determining antibiotic resistance. The following indicator genes were proposed that should absolutely be controlled in environmental biofilms: intl1, sul2, sul1, tetA, blaOXA, and blaTEM. In both studied WTPs, the highest number of copies was determined for the intI1 gene. Among the tested ARGs, the highest values were obtained for genes sul1 and tetA. The qPCR analysis also revealed that the amounts of determined ARGs decreased in the following order: sulphonamides>carbapenems >tetracyclines > β-lactams >macrolides. The dominant bacterial types in all analysed samples were Proteobacteria and Bacteroidetes. Both ARGs and bacterial biodiversity was determined rather by sampling site (spatial variation) than seasonality. The obtained results show that biofilms are reservoirs of ARGs. This may affect the microbiological quality of water entering the water system. It is therefore necessary to include their analysis in the classical studies of water quality.

Assembling the perfect bacterial genome using Oxford Nanopore and Illumina sequencing

A perfect bacterial genome assembly is one where the assembled sequence is an exact match for the organism's genome-each replicon sequence is complete and contains no errors. While this has been difficult to achieve in the past, improvements in long-read sequencing, assemblers, and polishers have brought perfect assemblies within reach. Here, we describe our recommended approach for assembling a bacterial genome to perfection using a combination of Oxford Nanopore Technologies long reads and Illumina short reads: Trycycler long-read assembly, Medaka long-read polishing, Polypolish short-read polishing, followed by other short-read polishing tools and manual curation. We also discuss potential pitfalls one might encounter when assembling challenging genomes, and we provide an online tutorial with sample data (github.com/rrwick/perfect-bacterial-genome-tutorial).

Comparative genomics and informational content analysis uncovered internal regions of the core genes rpoD, pepN and gltX for an MLSA with genome-level resolving power within the genus Pseudomonas

In the field of prokaryotic taxonomy, there has been a recent transition towards phylogenomics as the gold standard approach. However, genome-based phylogenetics is still restrictive for its cost when managing large amounts of isolates. Fast, cheap, and taxonomically competent alternatives, like multilocus sequence analysis (MLSA) are thus recommendable. Nevertheless, the criteria for selecting the conserved genes for MLSA have not been explicit for different bacterial taxa, including the broadly diverse Pseudomonas genus. Here, we have carried out an unbiased and rational workflow to select internal sequence regions of Pseudomonas core genes (CG) for a MLSA with the best phylogenetic power, and with a resolution comparable to the genome-based ANI approach. A computational workflow was established to inspect 126 complete genomes of representatives from over 60 Pseudomonas species and subspecies, in order to identify the most informative CG internal regions and determine which combinations in sets of three partial CG sequences have comparable phylogenetic resolution to that of the current ANI standard. We found that the rpoD^346-1196-pepN^1711-2571-gltX^86-909 concatenated sequences were the best performing in terms of phylogenetic robustness and resulted highly sensitive and specific when contrasted with ANI. The rpoD-pepN-gltX MLSA was validated in silico and in vitro. Altogether, the results presented here supports the proposal of the rpoD-pepN-gltX MLSA as a fast, affordable, and robust phylogenetic tool for members of the Pseudomonas genus.

Enterococcus faecium 129 BIO 3B is classified as Enterococcus lactis 129 BIO 3B

Enterococcus faecium 129 BIO 3B is a lactic acid bacterium that has been safely used as a probiotic product for over 100 years. Recently, concerns about its safety have arisen because some species of E. faecium belong to the vancomycin-resistant enterococci. The groups of E. faecium with less pathogenic potential have been split into a separate species (Enterococcus lactis). In this study, I investigated the phylogenetic classification and safety of E. faecium 129 BIO 3B as well as E. faecium 129 BIO 3B-R, which is naturally resistant to ampicillin. Mass spectrometry and basic local alignment search tool analysis using specific gene regions failed to differentiate 3B and 3B-R into E. faecium or E. lactis. However, multilocus sequence typing successfully identified 3B and 3B-R as the same sequence types as E. lactis. Overall genome relatedness indices showed that 3B and 3B-R have high degrees of homology with E. lactis. Gene amplification was confirmed for 3B and 3B-R with E. lactis species-specific primers. The minimum inhibitory concentration of ampicillin was confirmed to be 2 µg/mL for 3B, which is within the safety standard for E. faecium set by the European Food Safety Authority. Based on the above results, E. faecium 129 BIO 3B and E. faecium 129 BIO 3B-R were classified as E. lactis. The absence of pathogenic genes except for fms21 in this study demonstrates that these bacteria are safe for use as probiotics.

Macroevolutionary patterns in intragenomic rDNA variability among planktonic foraminifera

Ribosomal intragenomic variability in prokaryotes and eukaryotes is a genomic feature commonly studied for its inflationary impact on molecular diversity assessments. However, the evolutionary mechanisms and distribution of this phenomenon within a microbial group are rarely explored. Here, we investigate the intragenomic variability in 33 species of planktonic foraminifera, calcifying marine protists, by inspecting 2,403 partial SSU sequences obtained from single-cell clone libraries. Our analyses show that polymorphisms are common among planktonic foraminifera species, but the number of polymorphic sites significantly differs among clades. With our molecular simulations, we could assess that most of these mutations are located in paired regions that do not affect the secondary structure of the SSU fragment. Finally, by mapping the number of polymorphic sites on the phylogeny of the clades, we were able to discuss the evolution and potential sources of intragenomic variability in planktonic foraminifera, linking this trait to the distinctive nuclear and genomic dynamics of this microbial group.

Ribosomal RNA operons define a central functional compartment in the Streptomyces chromosome

Streptomyces are prolific producers of specialized metabolites with applications in medicine and agriculture. These bacteria possess a large linear chromosome genetically compartmentalized: core genes are grouped in the central part, while terminal regions are populated by poorly conserved genes. In exponentially growing cells, chromosome conformation capture unveiled sharp boundaries formed by ribosomal RNA (rrn) operons that segment the chromosome into multiple domains. Here we further explore the link between the genetic distribution of rrn operons and Streptomyces genetic compartmentalization. A large panel of genomes of species representative of the genus diversity revealed that rrn operons and core genes form a central skeleton, the former being identifiable from their core gene environment. We implemented a new nomenclature for Streptomyces genomes and trace their rrn-based evolutionary history. Remarkably, rrn operons are close to pericentric inversions. Moreover, the central compartment delimited by rrn operons has a very dense, nearly invariant core gene content. Finally, this compartment harbors genes with the highest expression levels, regardless of gene persistence and distance to the origin of replication. Our results highlight that rrn operons are structural boundaries of a central functional compartment prone to transcription in Streptomyces.

Spatio-temporal changes in endosymbiont diversity and composition in the African cassava whitefly, Bemisia tabaci SSA1

Sap-sucking insects, including whiteflies, are amongst the most devastating and widely distributed organisms on the planet. They are often highly invasive and endosymbiont communities within these insects help them adapt to new or changing environments. Bemisia tabaci (Gennadius; Hemiptera: Aleyrodidae) whitefly species are vectors of more than 500 known plant-viruses and harbour highly diverse endosymbionts communities. To date, however, whitefly–endosymbiont interactions, community structure and their spatio-temporal changes are still poorly understood. In this study, we investigated the spatio-temporal changes in the composition and diversity of bacterial endosymbionts in the agricultural crop pest whitefly species, Bemisia tabaci sub-Saharan Africa 1-subgroup 1 and 2 (SSA1-SG1 and SSA1-SG2). 16S rRNA amplicon sequencing analysis was carried out to characterise endosymbiont compositionsin field-collected SSA1 (SSA1-SG1 and SSA1-SG2) populations infesting cassava in Uganda in 1997 and 2017. We detected Portiera, Arsenophonus, Wolbachia, Hamiltonella and Hemipteriphilus, with Arsenophonus and Wolbachia infections being predominant. Hemipteriphilus and Hamiltonella frequencies were very low and were detected in seven and two samples, respectively. Bacterial diversity based on three independent parameters including Simpson index, number of haplotypes and Bray–Curtis dissimilarity matrix was significantly higher in 1997 than in 2017. This period also coincided with the advent of super-abundant cassava-whitefly populations on cassava crops in Uganda. We discuss how endosymbionts may influence the biology and behaviour of whiteflies leading to population explosions.

Effect of Intragenomic Sequence Heterogeneity among Multiple 16S rRNA Genes on Species Identification of Elizabethkingia

Accurate identification of Elizabethkingia species mostly requires the use of molecular techniques, and 16S rRNA gene sequencing is generally considered the method of choice. In this study, we evaluated the effect of intraspecific diversity among the multiple copies of the 16S rRNA gene on the accuracy of species identification in the genus Elizabethkingia. Sequences of 16S rRNA genes obtained from the 32 complete whole-genome sequences of Elizabethkingia deposited in GenBank and from 218 clinical isolates collected from 5 hospitals in Taiwan were analyzed. Four or five copies of 16S rRNA were identified in the Elizabethkingia species with complete genome sequences. The dissimilarity among the copies of the16S rRNA gene was <1% in all Elizabethkingia strains. E. meningoseptica demonstrated a significantly higher rate of nucleotide variations in the 16S rRNA than did E. anophelis (P = 0.011). Nucleotide alterations occurred more frequently in regions V2 and V6 than in other hypervariable regions (P < 0.001). E. meningoseptica, E. anophelis, and E. argenteiflava strains were clustered distinctly in the phylogenetic tree inferred from 16S rRNA genes, and the intragenomic variation of gene sequences had no profound effect on the classification of taxa. However, E. miricola, E. bruuniana, E. ursingii, and E. occulta were grouped closely in the phylogenetic analysis, and the variation among the multiple copies of the 16S rRNA in one E. ursingii strain affected species classification. Other marker genes may be required to supplement the species classification of closely related taxa in the genus Elizabethkingia. IMPORTANCE Incorrect identification of bacterial species would influence the epidemiology and clinical analysis of patients infected with Elizabethkingia. The results of the present study suggest that 16S rRNA gene sequencing should not be considered the gold standard for the accurate identification of Elizabethkingia species.

High-resolution phylogenetic and population genetic analysis of microbial communities with RoC-ITS

Microbial communities are inter-connected systems of incredible complexity and dynamism that play crucial roles in health, energy, and the environment. To better understand microbial communities and how they respond to change, it is important to know which microbes are present and their relative abundances at the greatest taxonomic resolution possible. Here, we describe a novel protocol (RoC-ITS) that uses the single-molecule Nanopore sequencing platform to assay the composition of microbial communities at the subspecies designation. Using rolling-circle amplification, this methodology produces long-read sequences from a circular construct containing the complete 16S ribosomal gene and the neighboring internally transcribed spacer (ITS). These long reads can be used to generate a high-fidelity circular consensus sequence. Generally, the ribosomal 16S gene provides phylogenetic information down to the species-level, while the much less conserved ITS region contains strain-level information. When linked together, this combination of markers allows for the identification of individual ribosomal units within a specific organism and the assessment of their relative stoichiometry, as well as the ability to monitor subtle shifts in microbial community composition with a single generic assay. We applied RoC-ITS to an artificial microbial community that was also sequenced using the Illumina platform, to assess its accuracy in quantifying the relative abundance and identity of each species.

Spectral imaging and nucleic acid mimics fluorescence in situ hybridization (SI-NAM-FISH) for multiplex detection of clinical pathogens

The application of nucleic acid mimics (NAMs), such as locked nucleic acid (LNA) and 2′-O-methyl-RNA (2’OMe), has improved the performance of fluorescence in situ hybridization (FISH) methods for the detection/location of clinical pathogens since they provide design versatility and thermodynamic control. However, an important limitation of FISH techniques is the low number of distinguishable targets. The use of filters in fluorescence image acquisition limits the number of fluorochromes that can be simultaneously differentiated. Recent advances in fluorescence spectral image acquisition have allowed the unambiguous identification of several microorganisms in a single sample. In this work, we aimed to combine NAM-FISH and spectral image analysis to develop and validate a new FISH variant, the spectral imaging-NAM-FISH (SI-NAM-FISH), that allows a multiplexed, robust and rapid detection of clinical pathogens. In the first stage, to implement/validate the method, we have selected seven fluorochromes with distinct spectral properties and seven bacterial species (Pseudomonas aeruginosa, Citrobacter freundii, Staphylococcus aureus, Enterococcus faecalis, Klebsiella pneumoniae, Escherichia coli, and Acinetobacter calcoaceticus). As a strong variation in fluorescence intensities is found between species and between fluorochromes, seven versions of a EUB LNA/2’OMe probe, each conjugated to one of seven fluorochromes, were used to rank species/fluorochromes by FISH and then optimize species/fluorochrome pairing. Then, final validation tests were performed using mixed populations to evaluate the potential of the technique for separating/quantifying the different targets. Overall, validation tests with different proportions of bacteria labeled with the respective fluorochrome have shown the ability of the method to correctly distinguish the species.

The Bacillaceae-1 RNA motif comprises two distinct classes

Non-coding RNAs are key regulatory players in bacteria. Many computationally predicted non-coding RNAs, however, lack functional associations. An example is the Bacillaceae-1 RNA motif, whose Rfam model consists of two hairpin loops. We find the motif conserved in nine of 13 non-pathogenic strains of the genus Bacillus but only in one pathogenic strain. To elucidate functional characteristics, we studied 118 hits of the Rfam model in 11 Bacillus spp. and found two distinct classes based on the ensemble diversity of their RNA secondary structure and the genomic context concerning the ribosomal RNA (rRNA) cluster. Forty hits are associated with the rRNA cluster, of which all 19 hits upstream flanking of 16S rRNA have a reverse complementary structure of low structural diversity. Fifty-two hits have large ensemble diversity, of which 38 are located between two coding genes. For eight hits in Bacillus subtilis, we investigated public expression data under various conditions and observed either the forward or the reverse complementary motif expressed. Five hits are associated with the rRNA cluster. Four of them are located upstream of the 16S rRNA and are not transcriptionally active, but instead, their reverse complements with low structural diversity are expressed together with the rRNA cluster. The three other hits are located between two coding genes in non-conserved genomic loci. Two of them are independently expressed from their surrounding genes and are structurally diverse. In summary, we found that Bacillaceae-1 RNA motifs upstream flanking of ribosomal RNA clusters tend to have one stable structure with the reverse complementary motif expressed in B. subtilis. In contrast, a subgroup of intergenic motifs has the thermodynamic potential for structural switches.

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.

Comamonas fluminis sp. nov., isolated from the Han River, Republic of Korea

A Gram-stain-negative, aerobic and motile bacterial strain, designated CJ34^T, was isolated from Han River water in the Republic of Korea. Strain CJ34^T grew optimally on tryptic soy agar at 30 °C and pH 7.0 in the absence of NaCl. Results of phylogenetic analysis based on 16S rRNA gene sequence showed that strain CJ34^T belonged to the genus Comamonas within the family Comamonadaceae and was most closely related to Comamonas testosteroni ATCC 11996^T and Comamonas thiooxydans DSM 17888^T (both 98.63 % similarity). The average nucleotide identity values between strain CJ34^T and two closely related type strains C. testosteroni ATCC 11996^T and C. thiooxydans DSM 17888^T were 82.77 and 82.73 %, respectively. The major isoprenoid quinone of strain CJ34^T was ubiquinone Q-8. The major cellular fatty acids of strain CJ34^T were C_16 : 0, C_16 : 1  ω6c and/or C_16 : 1  ω7c and C_18 : 1  ω6c and/or C_18 : 1  ω7c. The predominant polar lipids of strain CJ34^T were diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylglycerol and an unidentified aminophospholipid. Whole genome sequencing revealed that strain CJ34^T had a genome of 4.9 Mbp and the G+C content of the genomic DNA was 59.73 mol%. On the basis of the results of this polyphasic taxonomy study, strain CJ34^T represents a novel species in the genus Comamonas, for which the name Comamonas fluminis sp. nov. is proposed. The type strain is CJ34^T (=KACC 22237^T=JCM 34454^T).

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.

Curtobacterium spp. and Curtobacterium flaccumfaciens: Phylogeny, Genomics-Based Taxonomy, Pathogenicity, and Diagnostics

The genus of Curtobacterium, belonging to the Microbacteriaceae family of the Actinomycetales order, includes economically significant pathogenic bacteria of soybeans and other agricultural crops. Thorough phylogenetic and full-genome analysis using the latest genomic data has demonstrated a complex and contradictory taxonomic picture within the group of organisms classified as the Curtobacterium species. Based on these data, it is possible to delineate about 50 new species and to reclassify a substantial part of the Curtobacterium strains. It is suggested that 53 strains, including most of the Curtobacterium flaccumfaciens pathovars, can compose a monophyletic group classified as C. flaccumfaciens. A genomic analysis using the most recent inventory of bacterial chromosomal and plasmid genomes deposited to GenBank confirmed the possible role of Microbacteriaceae plasmids in pathogenicity and demonstrated the existence of a group of related plasmids carrying virulence factors and possessing a gene distantly related to DNA polymerase found in bacteriophages and archaeal and eukaryotic viruses. A PCR diagnostic assay specific to the genus Curtobacterium was developed and tested. The presented results assist in the understanding of the evolutionary relations within the genus and can lay the foundation for further taxonomic updates.

Recent advances in understanding the ecology of the filamentous bacteria responsible for activated sludge bulking

Activated sludge bulking caused by filamentous bacteria is still a problem in wastewater treatment plants around the world. Bulking is a microbiological problem, and so its solution on species-specific basis is likely to be reached only after their ecology, physiology and metabolism is better understood. Culture-independent molecular methods have provided much useful information about this group of organisms, and in this review, the methods employed and the information they provide are critically assessed. Their application to understanding bulking caused by the most frequently seen filament in Japan, 'Ca. Kouleothrix', is used here as an example of how these techniques might be used to develop control strategies. Whole genome sequences are now available for some of filamentous bacteria responsible for bulking, and so it is possible to understand why these filaments might thrive in activated sludge plants, and provide clues as to how eventually they might be controlled specifically.

A critical analysis of research methods and experimental models to study the root canal microbiome

Endodontic microbiology deals with the study of the microbial aetiology and pathogenesis of pulpal and periradicular inflammatory diseases. Research in endodontic microbiology started almost 130 years ago and since then has mostly focussed on establishing and confirming the infectious aetiology of apical periodontitis, identifying the microbial species associated with the different types of endodontic infections and determining the efficacy of treatment procedures in eradicating or controlling infection. Diverse analytical methods have been used over the years, each one with their own advantages and limitations. In this review, the main features and applications of the most used technologies are discussed, and advice is provided to improve study designs in order to properly address the scientific questions and avoid setbacks that can compromise the results. Finally, areas of future research are described.

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.

Unexpected Discovery of Hypermutator Phenotype Sounds the Alarm for Quality Control Strains

Microbial strains with high genomic stability are particularly sought after for testing the quality of commercial microbiological products, such as biological media and antibiotics. Yet, using mutation-accumulation experiments and de novo assembled complete genomes based on Nanopore long-read sequencing, we find that the widely used quality-control strain Shewanella putrefaciens ATCC-8071, also a facultative pathogen, is a hypermutator, with a base-pair substitution mutation rate of 2.42 × 10-8 per nucleotide site per cell division, ∼146-fold greater than that of the wild-type strain CGMCC-1.6515. Using complementation experiments, we confirm that mutL dysfunction, which was a recent evolutionary event, is the cause for the high mutation rate of ATCC-8071. Further analyses also give insight into possible relationships between mutation and genome evolution in this important bacterium. This discovery of a well-known strain being a hypermutator necessitates screening the mutation rate of bacterial strains before any quality control or experiments.

Exploring the Microdiversity Within Marine Bacterial Taxa: Toward an Integrated Biogeography in the Southern Ocean

Most of the microbial biogeographic patterns in the oceans have been depicted at the whole community level, leaving out finer taxonomic resolution (i.e., microdiversity) that is crucial to conduct intra-population phylogeographic study, as commonly done for macroorganisms. Here, we present a new approach to unravel the bacterial phylogeographic patterns combining community-wide survey by 16S rRNA gene metabarcoding and intra-species resolution through the oligotyping method, allowing robust estimations of genetic and phylogeographic indices, and migration parameters. As a proof-of-concept, we focused on the bacterial genus Spirochaeta across three distant biogeographic provinces of the Southern Ocean; maritime Antarctica, sub-Antarctic Islands, and Patagonia. Each targeted Spirochaeta operational taxonomic units were characterized by a substantial intrapopulation microdiversity, and significant genetic differentiation and phylogeographic structure among the three provinces. Gene flow estimations among Spirochaeta populations support the role of the Antarctic Polar Front as a biogeographic barrier to bacterial dispersal between Antarctic and sub-Antarctic provinces. Conversely, the Antarctic Circumpolar Current appears as the main driver of gene flow, connecting sub-Antarctic Islands with Patagonia and maritime Antarctica. Additionally, historical processes (drift and dispersal limitation) govern up to 86% of the spatial turnover among Spirochaeta populations. Overall, our approach bridges the gap between microbial and macrobial ecology by revealing strong congruency with macroorganisms distribution patterns at the populational level, shaped by the same oceanographic structures and ecological processes.

Characterization of Non-selected Intermolecular Gene Conversion in the Polyploid Haloarchaeon Haloferax volcanii

Gene conversion is defined as the non-reciprocal transfer of genetic information from one site to a homologous, but not identical site of the genome. In prokaryotes, gene conversion can increase the variance of sequences, like in antigenic variation, but can also lead to a homogenization of sequences, like in the concerted evolution of multigene families. In contrast to these intramolecular mechanisms, the intermolecular gene conversion in polyploid prokaryotes, which leads to the equalization of the multiple genome copies, has hardly been studied. We have previously shown the intermolecular gene conversion in halophilic and methanogenic archaea is so efficient that it can be studied without selecting for conversion events. Here, we have established an approach to characterize unselected intermolecular gene conversion in Haloferax volcanii making use of two genes that encode enzymes involved in carotenoid biosynthesis. Heterozygous strains were generated by protoplast fusion, and gene conversion was quantified by phenotype analysis or/and PCR. It was verified that unselected gene conversion is extremely efficient and it was shown that gene conversion tracts are much longer than in antigenic variation or concerted evolution in bacteria. Two sites were nearly always co-converted when they were 600 bp apart, and more than 30% co-conversion even occurred when two sites were 5 kbp apart. The gene conversion frequency was independent from the extent of genome differences, and even a one nucleotide difference triggered conversion.

Genomic analysis of high copy-number sequences for the targeted detection of Listeria species using a flow-through surveillance system

The bacterial foodborne pathogen Listeria monocytogenes has been implicated in fresh produce outbreaks with a significant economic impact. Given that L. monocytogenes is widespread in the environment, food production facilities constantly monitor for the presence of Listeria species. To develop a surveillance platform for food processing facilities, this study conducted a comparative genomic analysis for the identification of conserved high copy sequences in the ribosomal RNA of Listeria species. Simulated folding was performed to assess RNA accessibility in the identified genomic regions targeted for detection, and the developed singleplex assay accurately detected cell amounts lower than 5 cells, while no signals were detected for non-targeted bacteria. The singleplex assay was subsequently tested with a flow-through system, consisting of a DNA aptamer-capture step, followed by sample concentration and mechanical lysis for the detection of Listeria species. Validation experiments indicated the continuous flow-through system accurately detected Listeria species at low cell concentrations.

The maximum growth rate hypothesis is correct for eukaryotic photosynthetic organisms, but not cyanobacteria

The (maximum) growth rate (µ_max ) hypothesis predicts that cellular and tissue phosphorus (P) concentrations should increase with increasing growth rate, and RNA should also increase as most of the P is required to make ribosomes. Using published data, we show that though there is a strong positive relationship between the µ_max of all photosynthetic organisms and their P content (% dry weight), leading to a relatively constant P productivity, the relationship with RNA content is more complex. In eukaryotes there is a strong positive relationship between µ_max and RNA content expressed as % dry weight, and RNA constitutes a relatively constant 25% of total P. In prokaryotes the rRNA operon copy number is the important determinant of the amount of RNA present in the cell. The amount of phospholipid expressed as % dry weight increases with increasing µ_max in microalgae. The relative proportions of each of the five major P-containing constituents is remarkably constant, except that the proportion of RNA is greater and phospholipids smaller in prokaryotic than eukaryotic photosynthetic organisms. The effect of temperature differences between studies was minor. The evidence for and against P-containing constituents other than RNA being involved with ribosome synthesis and functioning is discussed.

Monoassociation with bacterial isolates reveals the role of colonization, community complexity and abundance on locomotor behavior in larval zebrafish

Background

Across taxa, animals with depleted intestinal microbiomes show disrupted behavioral phenotypes. Axenic (i.e., microbe-free) mice, zebrafish, and fruit flies exhibit increased locomotor behavior, or hyperactivity. The mechanism through which bacteria interact with host cells to trigger normal neurobehavioral development in larval zebrafish is not well understood. Here, we monoassociated zebrafish with either one of six different zebrafish-associated bacteria, mixtures of these host-associates, or with an environmental bacterial isolate.

Results

As predicted, the axenic cohort was hyperactive. Monoassociation with three different host-associated bacterial species, as well as with the mixtures, resulted in control-like locomotor behavior. Monoassociation with one host-associate and the environmental isolate resulted in the hyperactive phenotype characteristic of axenic larvae, while monoassociation with two other host-associated bacteria partially blocked this phenotype. Furthermore, we found an inverse relationship between the total concentration of bacteria per larvae and locomotor behavior. Lastly, in the axenic and associated cohorts, but not in the larvae with complex communities, we detected unexpected bacteria, some of which may be present as facultative predators.

Conclusions

These data support a growing body of evidence that individual species of bacteria can have different effects on host behavior, potentially related to their success at intestinal colonization. Specific to the zebrafish model, our results suggest that differences in the composition of microbes in fish facilities could affect the results of behavioral assays within pharmacological and toxicological studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-020-00069-x.

Plasmid localization of sole rrn operon in genomes of Oecophyllibacter saccharovorans (Acetobacteraceae)

The bacterium Oecophyllibacter saccharovorans of family Acetobacteraceae is a symbiont of weaver ant Oecophylla smaragdina. In our previous study, we published the finding of novel O. saccharovorans strains Ha5^T, Ta1 and Jb2 (Chua et al. 2020) but their plasmid sequences have not been reported before. Here, we demonstrate for the first time that the sole rrn operon of their genomes was detected on a 6.6 kb circular replicon. This replicon occurred in high copy number, much smaller size and lower G + C content than the main chromosome. Based on these features, the 6.6 kb circular replicon was regarded as rrn operon-containing plasmid. Further restriction analysis on the plasmids confirmed their circular conformation. A Southern hybridization analysis also corroborated the presence of 16S rRNA gene and thus the rrn operon on a single locus in the genome of the O. saccharovorans strains. However, similar genome architecture was not observed in other closely related bacterial strains. Additional survey also detected no plasmid-borne rrn operon in available genomes of validly described taxa of family Acetobacteraceae. To date, plasmid localization of rrn operon is rarely documented. This study reports the occurrence of rrn operon on the smallest bacterial plasmid in three O. saccharovorans strains and discusses its possible importance in enhancing their competitive fitness as bacterial symbiont of O. smaragdina.

Production of the antimicrobial compound tetrabromopyrrole and the Pseudomonas quinolone system precursor, 2-heptyl-4-quinolone, by a novel marine species Pseudoalteromonas galatheae sp. nov

Novel antimicrobials are urgently needed due to the rapid spread of antibiotic resistant bacteria. In a genome-wide analysis of Pseudoalteromonas strains, one strain (S4498) was noticed due to its potent antibiotic activity. It did not produce the yellow antimicrobial pigment bromoalterochromide, which was produced by several related type strains with which it shared less than 95% average nucleotide identity. Also, it produced a sweet-smelling volatile not observed from other strains. Mining the genome of strain S4498 using the secondary metabolite prediction tool antiSMASH led to eight biosynthetic gene clusters with no homology to known compounds, and synteny analyses revealed that the yellow pigment bromoalterochromide was likely lost during evolution. Metabolome profiling of strain S4498 using HPLC-HRMS analyses revealed marked differences to the type strains. In particular, a series of quinolones known as pseudanes were identified and verified by NMR. The characteristic odor of the strain was linked to the pseudanes. The highly halogenated compound tetrabromopyrrole was detected as the major antibacterial component by bioassay-guided fractionation. Taken together, the polyphasic analysis demonstrates that strain S4498 belongs to a novel species within the genus Pseudoalteromonas, and we propose the name Pseudoalteromonas galatheae sp. nov. (type strain S4498T = NCIMB 15250T = LMG 31599T).

Enhanced Symbiotic Characteristics in Bacterial Genomes with the Disruption of rRNA Operon

Simple Summary

Exploring the genomic changes that organisms have undergone to adapt to their specific environment is one of the most important processes in ecology and evolutionary biology. Here, we found that almost all rRNA operon-unlinked bacteria are symbiotic bacteria, which could be evidence of specific selective pressures in symbionts like genome reduction. This is meaningful and suggests that not only does the copy number variation of the rRNA operon sensitively respond to the bacterial lifestyle, but structural modification can also strongly reflect adaptation to the surrounding environmental conditions.

Abstract

Ribosomal RNA is an indispensable molecule in living organisms that plays an essential role in protein synthesis. Especially in bacteria, 16S, 23S, and 5S rRNAs are usually co-transcribed as operons. Despite the positive effects of rRNA co-transcription on growth and reproduction rate, a recent study revealed that bacteria with unlinked rRNA operons are more widespread than expected. However, it is still unclear why the rRNA operon is broken. Here, we explored rRNA operon linkage status in 15,898 bacterial genomes and investigated whether they have common features or lifestyles; 574 genomes were found to have unlinked rRNA operons and tended to be phylogenetically conserved. Most of them were symbionts and showed enhanced symbiotic genomic features such as reduced genome size and high adenine–thymine (AT) content. In an eggNOG-mapper analysis, they were also found to have significantly fewer genes than rRNA operon-linked bacteria in the “transcription” and “energy production and conversion in metabolism” categories. These genomes also tend to decrease RNases related to the synthesis of ribosomes and tRNA processing. Based on these results, the disruption of the rRNA operon seems to be one of the tendencies associated with the characteristics of bacteria requiring a low dynamic range.

Genomics Reveals the Metabolic Potential and Functions in the Redistribution of Dissolved Organic Matter in Marine Environments of the Genus Thalassotalea

Members of the bacterial genus Thalassotalea have been isolated recently from various marine environments, including marine invertebrates. A metagenomic study of the Deepwater Horizon oil plume has identified genes involved in aromatic hydrocarbon degradation in the Thalassotalea genome, shedding light on its potential role in the degradation of crude oils. However, the genomic traits of the genus are not well-characterized, despite the ability of the species to degrade complex natural compounds, such as agar, gelatin, chitin, or starch. Here, we obtained a complete genome of a new member of the genus, designated PS06, isolated from marine sediments containing dead marine benthic macroalgae. Unexpectedly, strain PS06 was unable to grow using most carbohydrates as sole carbon sources, which is consistent with the finding of few ABC transporters in the PS06 genome. A comparative analysis of 12 Thalassotalea genomes provided insights into their metabolic potential (e.g., microaerobic respiration and carbohydrate utilization) and evolutionary stability [including a low abundance of clustered regularly interspaced short palindromic repeats (CRISPR) loci and prophages]. The diversity and frequency of genes encoding extracellular enzymes for carbohydrate metabolism in the 12 genomes suggest that members of Thalassotalea contribute to nutrient cycling by the redistribution of dissolved organic matter in marine environments. Our study improves our understanding of the ecological and genomic properties of the genus Thalassotalea.

Phylogenetic comparison between Type IX Secretion System (T9SS) protein components suggests evidence of horizontal gene transfer

Porphyromonas gingivalis is one of the major bacteria that causes periodontitis. Chronic periodontitis is a severe form of periodontal disease that ultimately leads to tooth loss. Virulence factors that contribute to periodontitis are secreted by Type IX Secretion System (T9SS). There are aspects of T9SS protein components that have yet to be characterised. Thus, the aim of this study is to investigate the phylogenetic relationship between members of 20 T9SS component protein families. The Bayesian Inference (BI) trees for 19 T9SS protein components exhibit monophyletic clades for all major classes under Bacteroidetes with strong support for the monophyletic clades or its subclades that is consistent with phylogeny exhibited by the constructed BI tree of 16S rRNA. The BI tree of PorR is different from the 19 BI trees of T9SS protein components as it does not exhibit monophyletic clades for all major classes under Bacteroidetes. There is strong support for the phylogeny exhibited by the BI tree of PorR which deviates from the phylogeny based on 16S rRNA. Hence, it is possible that the porR gene is subjected to horizontal transfer as it is known that virulence factor genes could be horizontally transferred. Seven genes (porR included) that are involved in the biosynthesis of A-LPS are found to be flanked by insertion sequences (IS5 family transposons). Therefore, the intervening DNA segment that contains the porR gene might be transposed and subjected to conjugative transfer. Thus, the seven genes can be co-transferred via horizontal gene transfer. The BI tree of UgdA does not exhibit monophyletic clades for all major classes under Bacteroidetes which is similar to the BI tree of PorR (both are a part of the seven genes). Both BI trees also exhibit similar topology as the four identified clusters with strong support and have similar relative positions to each other in both BI trees. This reinforces the possibility that porR and the other six genes might be horizontally transferred. Other than the BI tree of PorR, the 19 other BI trees of T9SS protein components also exhibit evidence of horizontal gene transfer. However, their genes might undergo horizontal gene transfer less frequently compared to porR because the intervening DNA segment that contains porR is easily exchanged between bacteria under Bacteroidetes due to the presence of insertion sequences (IS5 family transposons) that flank it. In conclusion, this study can provide a better understanding about the phylogeny of T9SS protein components.

Effects of Habitat Partitioning on the Distribution of Bacterioplankton in Deep Lakes

In deep lakes, many investigations highlighted the existence of exclusive groups of bacteria adapted to deep oxygenated and hypoxic and anoxic hypolimnia. Nevertheless, the extent of bacterial strain diversity has been much less scrutinized. This aspect is essential for an unbiased estimation of genetic variation, biodiversity, and population structure, which are essential for studying important research questions such as biogeographical patterns, temporal and spatial variability and the environmental factors affecting this variability. This study investigated the bacterioplankton community in the epilimnetic layers and in the oxygenated and hypoxic/anoxic hypolimnia of five large and deep lakes located at the southern border of the Alps using high throughput sequencing (HTS) analyses (16S rDNA) and identification of amplicon sequence variants (ASVs) resolving reads differing by as little as one nucleotide. The study sites, which included two oligomictic (Garda and Como) and three meromictic lakes (Iseo, Lugano, and Idro) with maximum depths spanning from 124 to 410 m, were chosen among large lakes to represent an oxic-hypoxic gradient. The analyses showed the existence of several unique ASVs in the three layers of the five lakes. In the case of cyanobacteria, this confirmed previous analyses made at the level of strains or based on oligotyping methods. As expected, the communities in the hypoxic/anoxic monimolimnia showed a strong differentiation from the oxygenated layer, with the exclusive presence in single lakes of several unique ASVs. In the meromictic lakes, results supported the hypothesis that the formation of isolated monimolimnia sustained the development of highly diversified bacterial communities through ecological selection, leading to the establishment of distinctive biodiversity zones. The genera identified in these layers are well-known to activate a wide range of redox reactions at low O₂ conditions. As inferred from 16S rDNA data, the highly diversified and coupled processes sustained by the monimolimnetic microbiota are essential ecosystem services that enhance mineralization of organic matter and formation of reduced compounds, and also abatement of undesirable greenhouse gasses.

Conventional and organic soil management as divergent drivers of resident and active fractions of major soil food web constituents

Conventional agricultural production systems, typified by large inputs of mineral fertilizers and pesticides, reduce soil biodiversity and may negatively affect ecosystem services such as carbon fixation, nutrient cycling and disease suppressiveness. Organic soil management is thought to contribute to a more diverse and stable soil food web, but data detailing this effect are sparse and fragmented. We set out to map both the resident (rDNA) and the active (rRNA) fractions of bacterial, fungal, protozoan and metazoan communities under various soil management regimes in two distinct soil types with barley as the main crop. Contrasts between resident and active communities explained 22%, 14%, 21% and 25% of the variance within the bacterial, fungal, protozoan, and metazoan communities. As the active fractions of organismal groups define the actual ecological functioning of soils, our findings underline the relevance of characterizing both resident and active pools. All four major organismal groups were affected by soil management (p < 0.01), and most taxa showed both an increased presence and an enlarged activity under the organic regime. Hence, a prolonged organic soil management not only impacts the primary decomposers, bacteria and fungi, but also major representatives of the next trophic level, protists and metazoa.

The Phylogeny, Biodiversity, and Ecology of the Chloroflexi in Activated Sludge

It is now clear that several of the filamentous bacteria in activated sludge wastewater treatment plants globally, are members of the phylum Chloroflexi. They appear to be more commonly found in treatment plants designed to remove nitrogen (N) and phosphorus (P), most of which operate at long sludge ages and expose the biomass to anaerobic conditions. The Chloroflexi seem to play an important beneficial role in providing the filamentous scaffolding around which flocs are formed, to feed on the debris from lysed bacterial cells, to ferment carbohydrates and to degrade other complex polymeric organic compounds to low molecular weight substrates to support their growth and that of other bacterial populations. A few commonly extend beyond the floc surface, while others can align in bundles, which may facilitate interfloc bridging and hence generate a bulking sludge. Although several recent papers have examined the phylogeny and in situ physiology of Chloroflexi in activated sludge plants in Denmark, this review takes a wider look at what we now know about these filaments, especially their global distribution in activated sludge plants, and what their functional roles there might be. It also attempts to outline why such information might provide us with clues as to how their population levels may be manipulated, and the main research questions that need addressing to achieve these outcomes.

Targeting 16S rDNA for Stable Recombinant Gene Expression in Pseudomonas

Ribosomal RNA (rRNA) operons have recently been identified as promising sites for chromosomal integration of genetic elements in Pseudomonas putida, a bacterium that has gained considerable popularity as a microbial cell factory. We have developed a tool for targeted integration of recombinant genes into the rRNA operons of various Pseudomonas strains, where the native context of the rRNA clusters enables effective transcription of heterologous genes. However, a sufficient translation of foreign mRNA  transcriptionally fused to rRNA required optimization of RNA secondary structures, which was achieved utilizing synthetic ribozymes and a bicistronic design. The generated tool further enabled the characterization of the six rRNA promoter units of P. putida S12 under different growth conditions. The presence of multiple, almost identical rRNA operons in Pseudomonas also allowed the integration of multiple copies of heterologous genetic elements. The integration of two expression cassettes and the resulting disruption of rRNA units only moderately affects growth rates, and the constructs were highly stable over more than 160 generations.

Delineation of ecologically distinct units of marine Bacteroidetes in the Northwestern Mediterranean Sea

Bacteroidetes is one of the dominant phyla of ocean bacterioplankton, yet its diversity and population structure is poorly understood. To advance in the delineation of ecologically meaningful units within this group, we constructed near full-length 16S rRNA gene clone libraries from contrasting marine environments in the NW Mediterranean. Based on phylogeny and the associated ecological variables (depth and season), 24 different Bacteroidetes clades were delineated. By considering their relative abundance (from iTag amplicon sequencing studies), we described the distribution patterns of each of these clades, delimiting them as Ecologically Significant Taxonomic Units (ESTUs). Spatially, there was almost no overlap among ESTUs at different depths. In deep waters there was predominance of Owenweeksia, Leeuwenhoekiella, Muricauda-related genera, and some depth-associated ESTUs within the NS5 and NS2b marine clades. Seasonally, multi-annual dynamics of recurring ESTUs were present with dominance of some ESTUs within the NS4, NS5 and NS2b marine clades along most of the year, but with variable relative frequencies between months. A drastic change towards the predominance of Formosa-related ESTUs and one ESTU from the NS5 marine clade was typically present after the spring bloom. Even though there are no isolates available for these ESTUs to determine their physiology, correlation analyses identified the environmental preference of some of them. Overall, our results suggest that there is a high degree of niche specialisation within these closely related clades. This work constitutes a step forward in disentangling the ecology of marine Bacteroidetes, which are essential players in organic matter processing in the oceans.