rRNA operon multiplicity in Escherichia coli and the physiological implications of rrn inactivation

Adaptive traits of Nitrosocosmicus clade ammonia-oxidizing archaea

Nitrification is a core process in the global nitrogen (N) cycle mediated by ammonia-oxidizing microorganisms, including ammonia-oxidizing archaea (AOA) as a key player. Although much is known about AOA abundance and diversity across environments, the genetic drivers of the ecophysiological adaptations of the AOA are often less clearly defined. This is especially true for AOA within the genus Nitrosocosmicus, which have several unique physiological traits (e.g., high substrate tolerance, low substrate affinity, and large cell size). To better understand what separates the physiology of Nitrosocosmicus AOA, we performed comparative genomics with genomes from 39 cultured AOA, including five Nitrosocosmicus AOA. The absence of a canonical high-affinity type ammonium transporter and typical S-layer structural genes was found to be conserved across all Nitrosocosmicus AOA. In agreement, cryo-electron tomography confirmed the absence of a visible outermost S-layer structure, which has been observed in other AOA. In contrast to other AOA, the cryo-electron tomography highlighted the possibility that Nitrosocosmicus AOA may possess a glycoprotein or glycolipid-based glycocalyx cell covering outer layer. Together, the genomic, physiological, and metabolic properties revealed in this study provide insight into niche adaptation mechanisms and the overall ecophysiology of members of the Nitrosocosmicus clade in various terrestrial ecosystems.

IMPORTANCE

Nitrification is a vital process within the global biogeochemical nitrogen cycle but plays a significant role in the eutrophication of aquatic ecosystems and the production of the greenhouse gas nitrous oxide (N2O) from industrial agriculture ecosystems. While various types of ammonia-oxidizing microorganisms play a critical role in the N cycle, ammonia-oxidizing archaea (AOA) are often the most abundant nitrifiers in natural environments. Members of the genus Nitrosocosmicus are one of the prevalent AOA groups detected in undisturbed terrestrial ecosystems and have previously been reported to possess a range of physiological characteristics that set their physiology apart from other AOA species. This study provides significant progress in understanding these unique physiological traits and their genetic drivers. Our results highlight how physiological studies based on comparative genomics-driven hypotheses can contribute to understanding the unique niche of Nitrosocosmicus AOA.

RFGR: Repeat Finder for Complete and Assembled Whole Genomes and NGS Reads

Repetitive DNA sequences cause genomic instability and are important genetic markers. Identification of repeats is a critical step in genome annotation and analysis. On the other hand, repeats also pose a technical challenge for genome assembly and alignment programs using NGS data. RFGR is a comprehensive tool that can find exact repetitive sequences in complete genomes and assembled genomes, as well as NGS reads of prokaryotes. For complete genomes, RFGR uses a suffix trees to find seed repeats of repetitive sequences of fixed length with indels. For assembled genomes, RFGR uses a modified Bowtie aligner to find seed repeats of exact repetitive sequences in the contigs/ scaffolds, which are then extended to maximal repeats. The repeats are classified and for repeats near a gene, RFGR reports the gene as well. For the control dataset of E. coli UTI89 and E. coli K12, RFGR reports 35,141 and 49,352 repeats, respectively. For NGS reads, RFGR uses the frequency of the repetitive k-mers to determine FASTQ reads containing repetitive sequences and removes them from the dataset. An E. coli K12 NGS dataset pre-processed using RFGR, on comparison with the original dataset, gives an improved assembly. The N50 value improves by 22.86% with a decrease in size of the assembly graph by nearly 50%. Thus, with RFGR, we achieve a better assembly with reduced computation. RFGR can be improved in terms of the length of the minimum repeat found, extending to find approximate repeats and to be applicable to Eukaryotes as well.

Ribosome Pool Engineering Increases Protein Biosynthesis Yields

The biosynthetic capability of the bacterial ribosome motivates efforts to understand and harness sequence-optimized versions for synthetic biology. However, functional differences between natively occurring ribosomal RNA (rRNA) operon sequences remain poorly characterized. Here, we use an in vitro ribosome synthesis and translation platform to measure protein production capabilities of ribosomes derived from all unique combinations of 16S and 23S rRNAs from seven distinct Escherichia coli rRNA operon sequences. We observe that polymorphisms that distinguish native E. coli rRNA operons lead to significant functional changes in the resulting ribosomes, ranging from negligible or low gene expression to matching the protein production activity of the standard rRNA operon B sequence. We go on to generate strains expressing single rRNA operons and show that not only do some purified in vivo expressed homogeneous ribosome pools outperform the wild-type, heterogeneous ribosome pool but also that a crude cell lysate made from the strain expressing only operon A ribosomes shows significant yield increases for a panel of medically and industrially relevant proteins. We anticipate that ribosome pool engineering can be applied as a tool to increase yields across many protein biomanufacturing systems, as well as improve basic understanding of ribosome heterogeneity and evolution.

Structural basis of ribosomal 30S subunit degradation by RNase R

Protein synthesis is a major energy-consuming process of the cell that requires the controlled production1-3 and turnover4,5 of ribosomes. Although the past few years have seen major advances in our understanding of ribosome biogenesis, structural insight into the degradation of ribosomes has been lacking. Here we present native structures of two distinct small ribosomal 30S subunit degradation intermediates associated with the 3' to 5' exonuclease ribonuclease R (RNase R). The structures reveal that RNase R binds at first to the 30S platform to facilitate the degradation of the functionally important anti-Shine-Dalgarno sequence and the decoding-site helix 44. RNase R then encounters a roadblock when it reaches the neck region of the 30S subunit, and this is overcome by a major structural rearrangement of the 30S head, involving the loss of ribosomal proteins. RNase R parallels this movement and relocates to the decoding site by using its N-terminal helix-turn-helix domain as an anchor. In vitro degradation assays suggest that head rearrangement poses a major kinetic barrier for RNase R, but also indicate that the enzyme alone is sufficient for complete degradation of 30S subunits. Collectively, our results provide a mechanistic basis for the degradation of 30S mediated by RNase R, and reveal that RNase R targets orphaned 30S subunits using a dynamic mechanism involving an anchored switching of binding sites.

Genomic Insights into Paucibacter aquatile DH15, a Cyanobactericidal Bacterium, and Comparative Genomics of the Genus Paucibacter

Microcystis blooms threaten ecosystem function and cause substantial economic losses. Microorganism-based methods, mainly using cyanobactericidal bacteria, are considered one of the most ecologically sound methods to control Microcystis blooms. This study focused on gaining genomic insights into Paucibacter aquatile DH15 that exhibited excellent cyanobactericidal effects against Microcystis. Additionally, a pan-genome analysis of the genus Paucibacter was conducted to enhance our understanding of the ecophysiological significance of this genus. Based on phylogenomic analyses, strain DH15 was classified as a member of the species Paucibacter aquatile. The genome analysis supported that strain DH15 can effectively destroy Microcystis, possibly due to the specific genes involved in the flagellar synthesis, cell wall degradation, and the production of cyanobactericidal compounds. The pan-genome analysis revealed the diversity and adaptability of the genus Paucibacter, highlighting its potential to absorb external genetic elements. Paucibacter species were anticipated to play a vital role in the ecosystem by potentially providing essential nutrients, such as vitamins B7, B12, and heme, to auxotrophic microbial groups. Overall, our findings contribute to understanding the molecular mechanisms underlying the action of cyanobactericidal bacteria against Microcystis and shed light on the ecological significance of the genus Paucibacter.

Land use alters bacterial growth dynamics in soil

Microbial growth and mortality are major determinants of soil carbon cycling. We measured in situ growth dynamics of individual bacterial taxa in cropped and successional soils in response to a resource pulse. We hypothesized that land use imposes selection pressures on growth characteristics. We estimated growth and death for 453 and 73 taxa, respectively. The average generation time was 5.04 ± 6.28 (SD; range 0.7-63.5) days. Lag times were shorter in cultivated than successional soils and resource amendment decreased lag times. Taxa exhibiting the greatest growth response also exhibited the greatest mortality, indicative of boom-and-bust dynamics. We observed a bimodal growth rate distribution, representing fast- and slow-growing clusters. Both clusters grew more rapidly in successional soils, which had more organic matter, than cultivated soils. Resource amendment increased the growth rate of the slower growing but not the faster-growing cluster via a mixture of increased growth rates and species turnover, indicating that competitive dynamics constrain growth rates in situ. These two clusters show that copiotrophic bacteria in soils may be subdivided into different life history groups and that these subgroups respond independently to land use and resource availability.

Ribosome Abundance Control in Prokaryotes

Cell growth is an essential phenotype of any unicellular organism and it crucially depends on precise control of protein synthesis. We construct a model of the feedback mechanisms that regulate abundance of ribosomes in E. coli, a prototypical prokaryotic organism. Since ribosomes are needed to produce more ribosomes, the model includes a positive feedback loop central to the control of cell growth. Our analysis of the model shows that there can be only two coexisting equilibrium states across all 23 parameters. This precludes the existence of hysteresis, suggesting that the ribosome abundance changes continuously with parameters. These states are related by a transcritical bifurcation, and we provide an analytic formula for parameters that admit either state.

Ribosomal RNA gene operon copy number, a functional trait indicating the hydrocarbon degradation level of bacterial communities

The lack of universal indicators for predicting microbial biodegradation potential and assessing remediation effects limits the generalization of bioremediation. The community-level ribosomal RNA gene operon (rrn) copy number, an important functional trait, has the potential to serve as a key indicator of the bioremediation of organic pollutants. A meta-analysis based on 1275 samples from 26 hydrocarbon-related studies revealed a positive relationship between the microbial hydrocarbon biodegradation level and the community-level rrn copy number in soil, seawater and culture. Subsequently, a microcosm experiment was performed to decipher the community-level rrn copy number response mechanism during total petroleum hydrocarbon (TPH) biodegradation. The treatment combining straw with resuscitation-promoting factor (Rpf) exhibited the highest community-level rrn copy number and the most effective biodegradation compared with other treatments, and the initial TPH content (20,000 mg kg-1) was reduced by 67.67% after 77 days of incubation. TPH biodegradation rate was positively correlated with the average community-level rrn copy number (p = 0.001, R2 = 0.5781). Both meta and community analyses showed that rrn copy number may reflect the potential of hydrocarbon degradation and microbial dormancy. Our findings provide insight into the applicability of the community-level rrn copy number to assess bacterial biodegradation for pollution remediation.

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.

The layered costs and benefits of translational redundancy

The rate and accuracy of translation hinges upon multiple components - including transfer RNA (tRNA) pools, tRNA modifying enzymes, and rRNA molecules - many of which are redundant in terms of gene copy number or function. It has been hypothesized that the redundancy evolves under selection, driven by its impacts on growth rate. However, we lack empirical measurements of the fitness costs and benefits of redundancy, and we have poor a understanding of how this redundancy is organized across components. We manipulated redundancy in multiple translation components of Escherichia coli by deleting 28 tRNA genes, 3 tRNA modifying systems, and 4 rRNA operons in various combinations. We find that redundancy in tRNA pools is beneficial when nutrients are plentiful and costly under nutrient limitation. This nutrient-dependent cost of redundant tRNA genes stems from upper limits to translation capacity and growth rate, and therefore varies as a function of the maximum growth rate attainable in a given nutrient niche. The loss of redundancy in rRNA genes and tRNA modifying enzymes had similar nutrient-dependent fitness consequences. Importantly, these effects are also contingent upon interactions across translation components, indicating a layered hierarchy from copy number of tRNA and rRNA genes to their expression and downstream processing. Overall, our results indicate both positive and negative selection on redundancy in translation components, depending on a species' evolutionary history with feasts and famines.

The role of replication-induced chromosomal copy numbers in spatio-temporal gene regulation and evolutionary chromosome plasticity

For a coherent response to environmental changes, bacterial evolution has formed a complex transcriptional regulatory system comprising classical DNA binding proteins sigma factors and modulation of DNA topology. In this study, we investigate replication-induced gene copy numbers - a regulatory concept that is unlike the others not based on modulation of promoter activity but on replication dynamics. We show that a large fraction of genes are predominantly affected by transient copy numbers and identify cellular functions and central pathways governed by this mechanism in Escherichia coli. Furthermore, we show quantitatively that the previously observed spatio-temporal expression pattern between different growth phases mainly emerges from transient chromosomal copy numbers. We extend the analysis to the plant pathogen Dickeya dadantii and the biotechnologically relevant organism Vibrio natriegens. The analysis reveals a connection between growth phase dependent gene expression and evolutionary gene migration in these species. A further extension to the bacterial kingdom indicates that chromosome evolution is governed by growth rate related transient copy numbers.

A tRNA-derived fragment present in E. coli OMVs regulates host cell gene expression and proliferation

RNA-sequencing has led to a spectacular increase in the repertoire of bacterial sRNAs and improved our understanding of their biological functions. Bacterial sRNAs have also been found in outer membrane vesicles (OMVs), raising questions about their potential involvement in bacteria-host relationship, but few studies have documented this issue. Recent RNA-Sequencing analyses of bacterial RNA unveiled the existence of abundant very small RNAs (vsRNAs) shorter than 16 nt. These especially include tRNA fragments (tRFs) that are selectively loaded in OMVs and are predicted to target host mRNAs. Here, in Escherichia coli (E. coli), we report the existence of an abundant vsRNA, Ile-tRF-5X, which is selectively modulated by environmental stress, while remaining unaffected by inhibition of transcription or translation. Ile-tRF-5X is released through OMVs and can be transferred to human HCT116 cells, where it promoted MAP3K4 expression. Our findings provide a novel perspective and paradigm on the existing symbiosis between bacteria and human cells.

We previously outlined by RNA-Sequencing (RNA-seq) the existence of abundant very small (<16 nt) bacterial and eukaryote RNA (vsRNA) population with potential regulatory functions. However, it is not exceptional to see vsRNA species removed from the RNA-seq libraries or datasets because being considered as random degradation products. As a proof of concept, we present in this study a 13 nt in length isoleucine tRNA-derived fragment (Ile-tRF-5X) which is selectively modulated by nutritional and thermal stress while remaining unaffected by transcription and translation inhibitions. We also showed that OMVs and their Ile-tRF-5X vsRNAs are delivered into human HCT116 cells and both can promote host cell gene expression and proliferation. Ile-tRF-5X appears to regulate gene silencing properties of miRNAs by competition. Our findings provide a novel perspective and paradigm on the existing symbiosis between hosts and bacteria but also brings a new insight of host-pathogen interactions mediated by tRFs which remain so far poorly characterized in bacteria.

Gut microbiome reflect adaptation of earthworms to cave and surface environments

Background

Caves are special natural laboratories for most biota and the cave communities are unique. Establishing population in cave is accompanied with modifications in adaptability for most animals. To date, little is known about the survival mechanisms of soil animals in cave environments, albeit they play vital roles in most terrestrial ecosystems. Here, we investigated whether and how gut microbes would contribute to the adaptation of earthworms by comparing the gut microbiome of two earthworm species from the surface and caves.

Results

Two dominant earthworm species inhabited caves, i.e., Allolobophora chlorotica and Aporrectodea rosea. Compared with the counterparts on the surface, A. rosea significantly decreased population in the cave, while A. chlorotica didn’t change. Microbial taxonomic and phylogenetic diversities between the earthworm gut and soil environment were asynchronic with functional diversity, with functional gene diversity been always higher in earthworm gut than in soil, but species richness and phylogenetic diversity lower. In addition, earthworm gut microbiome were characterized by higher rrn operon numbers and lower network complexity than soil microbiota.

Conclusions

Different fitness of the two earthworm species in cave is likely to coincide with gut microbiota, suggesting interactions between host and gut microbiome are essential for soil animals in adapting to new environments. The functional gene diversity provided by gut microbiome is more important than taxonomic or phylogenetic diversity in regulating host adaptability. A stable and high-efficient gut microbiome, including microbiota and metabolism genes, encoded potential functions required by the animal hosts during the processes of adapting to and establishing in the cave environments. Our study also demonstrates how the applications of microbial functional traits analysis may advance our understanding of animal-microbe interactions that may aid animals to survive in extreme ecosystems.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-022-00200-0.

rRNA operon multiplicity as a bacterial genome stability insurance policy

Quick growth restart after upon encountering favourable environmental conditions is a major fitness contributor in natural environment. It is widely assumed that the time required to restart growth after nutritional upshift is determined by how long it takes for cells to synthesize enough ribosomes to produce the proteins required to reinitiate growth. Here we show that a reduction in the capacity to synthesize ribosomes by reducing number of ribosomal RNA (rRNA) operons (rrn) causes a longer transition from stationary phase to growth of Escherichia coli primarily due to high mortality rates. Cell death results from DNA replication blockage and massive DNA breakage at the sites of the remaining rrn operons that become overloaded with RNA polymerases (RNAPs). Mortality rates and growth restart duration can be reduced by preventing R-loop formation and improving DNA repair capacity. The same molecular mechanisms determine the duration of the recovery phase after ribosome-damaging stresses, such as antibiotics, exposure to bile salts or high temperature. Our study therefore suggests that a major function of rrn operon multiplicity is to ensure that individual rrn operons are not saturated by RNAPs, which can result in catastrophic chromosome replication failure and cell death during adaptation to environmental fluctuations.

Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales

The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel's lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.

Regulatory perturbations of ribosome allocation in bacteria reshape the growth proteome with a trade-off in adaptation capacity

Bacteria regulate their cellular resource allocation to enable fast growth-adaptation to a variety of environmental niches. We studied the ribosomal allocation, growth, and expression profiles of two sets of fast-growing mutants of Escherichia coli K-12 MG1655. Mutants with only three of the seven copies of ribosomal RNA operons grew faster than the wild-type strain in minimal media and show similar phenotype to previously studied fast-growing rpoB mutants. Comparing these two different regulatory perturbations (rRNA promoters or rpoB mutations), we show how they reshape the proteome for growth with a concomitant fitness cost. The fast-growing mutants shared downregulation of hedging functions and upregulated growth functions. They showed longer diauxic shifts and reduced activity of gluconeogenic promoters during glucose-acetate shifts, suggesting reduced availability of the RNA polymerase for expressing hedging proteome. These results show that the regulation of ribosomal allocation underlies the growth/hedging phenotypes obtained from laboratory evolution experiments.

•

Mutants with only three ribosomal operons grow faster than wild-type in minimal medium

•

Faster growth of mutants is achieved by increased ribosome content

•

Fast-growing mutants display reduced hedging expression and adaptation trade-offs

Sleeping ribosomes: Bacterial signaling triggers RaiA mediated persistence to aminoglycosides

Indole is a molecule proposed to be involved in bacterial signaling. We find that indole secretion is induced by sublethal tobramycin concentrations and increases persistence to aminoglycosides in V. cholerae. Indole transcriptomics showed increased expression of raiA, a ribosome associated factor. Deletion of raiA abolishes the appearance of indole dependent persisters to aminoglycosides, although its overexpression leads to 100-fold increase of persisters, and a reduction in lag phase, evocative of increased active 70S ribosome content, confirmed by sucrose gradient analysis. We propose that, under stress conditions, RaiA-bound inactive 70S ribosomes are stored as “sleeping ribosomes”, and are rapidly reactivated upon stress relief. Our results point to an active process of persister formation through ribosome protection during translational stress (e.g., aminoglycoside treatment) and reactivation upon antibiotic removal. Translation is a universal process, and these results could help elucidate a mechanism of persistence formation in a controlled, thus inducible way.

•

Indole is produced under sub-MIC tobramycin stress in V. cholerae and upregulates raiA

•

RaiA is involved in indole-dependent formation of aminoglycoside specific persisters

•

RaiA overexpression allows faster growth restart and increases 70S ribosome content

•

RaiA-bound inactive 70S ribosomes form intact and reactivable sleeping ribosome pools

Considerations on the Identity and Diversity of Organisms Affiliated with Sphingobacterium multivorum-Proposal for a New Species, Sphingobacterium paramultivorum

Plant biomass offers great potential as a sustainable resource, and microbial consortia are primordial in its bioconversion. The wheat-straw-biodegradative bacterial strain w15 has drawn much attention as a result of its biodegradative potential and superior degradation performance in bacterial-fungal consortia. Strain w15 was originally assigned to the species Sphingobacterium multivorum based on its 16S ribosomal RNA (rRNA) gene sequence. A closer examination of this taxonomic placement revealed that the sequence used has 98.9% identity with the ‘divergent’ 16S rRNA gene sequence of S. multivorum NCTC 11343^T, yet lower relatedness with the canonical 16S rRNA sequence. A specific region of the gene, located between positions 186 and 210, was found to be highly variable and determinative for the divergence. To solve the identity of strain w15, genome metrics and analyses of ecophysiological niches were undertaken on a selection of strains assigned to S. multivorum and related species. These analyses separated all strains into three clusters, with strain w15, together with strain BIGb0170, constituting a separate radiation, next to S. multivorum and S. siyangense. Moreover, the strains denoted FDAARGOS 1141 and 1142 were placed inside S. siyangense. We propose the renaming of strains w15 and BIGb0170 as members of the novel species, coined Sphingobacterium paramultivorum.

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.

Comparative genomic analysis of the genus Novosphingobium and the description of two novel species Novosphingobium aerophilum sp. nov. and Novosphingobium jiangmenense sp. nov

Members of the genus Novosphingobium are well known for their metabolically versatile and great application potential in pollution elimination. The three novel bacterial strains, designated 4Y4^T, 4Y9, and 1Y9A^T, were isolated from aquaculture water and characterized by using a polyphasic taxonomic approach. The 16S rRNA gene sequences phylogenetic analysis revealed that the three strains belonged to the genus Novosphingobium. The phylogenomic analysis indicated that the three strains formed two independent and robust branches distinct from all reference strains. The analyses of dDDH values and ANIs between the three strains and their relatives further demonstrated that the three strains represented two different novel genospecies. Comparative genomic analysis of the three isolates and 32 type strains of the genus Novosphingobium showed that the most important central metabolic pathways of these strains appeared to be similar, while specific and specialized metabolic pathways were flexible and variable among these strains. Chemotaxonomic characterization exhibited that the predominant cellular fatty acids were summed feature 8, summed feature 3, and C_14:0 2OH; the major polar lipids were diphosphatidylglycerol, phosphatidylethanolamine, phosphatidyldimethylethanolamine, phosphatidylglycerol, and sphingoglycolipid; the major respiratory quinone and polyamine were Q-10 and spermidine. The DNA G + C contents were 67.6 and 64.7 %. Based on the genotypic and phenotypic characteristics, strains 4Y4^T and 1Y9A^T are concluded to represent two novel species of the genus Novosphingobium, for which the names Novosphingobium aerophilum sp. nov. (type strain 4Y4^T = GDMCC 1.1828 ^T = KACC 21946 ^T) and Novosphingobium jiangmenense sp. nov. (type strain 1Y9A^T = GDMCC 1.1936 ^T = KACC 22085 ^T) are proposed.

Beyond horizontal gene transfer: the role of plasmids in bacterial evolution

Plasmids have a key role in bacterial ecology and evolution because they mobilize accessory genes by horizontal gene transfer. However, recent studies have revealed that the evolutionary impact of plasmids goes above and beyond their being mere gene delivery platforms. Plasmids are usually kept at multiple copies per cell, producing islands of polyploidy in the bacterial genome. As a consequence, the evolution of plasmid-encoded genes is governed by a set of rules different from those affecting chromosomal genes, and these rules are shaped by unusual concepts in bacterial genetics, such as genetic dominance, heteroplasmy or segregational drift. In this Review, we discuss recent advances that underscore the importance of plasmids in bacterial ecology and evolution beyond horizontal gene transfer. We focus on new evidence that suggests that plasmids might accelerate bacterial evolution, mainly by promoting the evolution of plasmid-encoded genes, but also by enhancing the adaptation of their host chromosome. Finally, we integrate the most relevant theoretical and empirical studies providing a global understanding of the forces that govern plasmid-mediated evolution in bacteria.

Translational regulation of environmental adaptation in bacteria

Bacteria must rapidly respond to both intracellular and environmental changes to survive. One critical mechanism to rapidly detect and adapt to changes in environmental conditions is control of gene expression at the level of protein synthesis. At each of the three major steps of translation-initiation, elongation, and termination-cells use stimuli to tune translation rate and cellular protein concentrations. For example, changes in nutrient concentrations in the cell can lead to translational responses involving mechanisms such as dynamic folding of riboswitches during translation initiation or the synthesis of alarmones, which drastically alter cell physiology. Moreover, the cell can fine-tune the levels of specific protein products using programmed ribosome pausing or inducing frameshifting. Recent studies have improved understanding and revealed greater complexity regarding long-standing paradigms describing key regulatory steps of translation such as start-site selection and the coupling of transcription and translation. In this review, we describe how bacteria regulate their gene expression at the three translational steps and discuss how translation is used to detect and respond to changes in the cellular environment. Finally, we appraise the costs and benefits of regulation at the translational level in bacteria.

In vitro thermal adaptation of mesophilic Acetobacter pasteurianus NBRC 3283 generates thermotolerant strains with evolutionary trade-offs

Thermotolerant strains are critical for low-cost high temperature fermentation. In this study, we carried out the thermal adaptation of A. pasteurianus IFO 3283–32 under acetic acid fermentation conditions using an experimental evolution approach from 37ºC to 40ºC. The adapted strain exhibited an increased growth and acetic acid fermentation ability at high temperatures, however, with the trade-off response of the opposite phenotype at low temperatures. Genome analysis followed by PCR sequencing showed that the most adapted strain had 11 mutations, a single 64-kb large deletion, and a single plasmid loss. Comparative phenotypic analysis showed that at least the large deletion (containing many ribosomal RNAs and tRNAs genes) and a mutation of DNA polymerase (one of the 11 mutations) critically contributed to this thermotolerance. The relationship between the phenotypic changes and the gene mutations are discussed, comparing with another thermally adapted A. pasteurianus strains obtained previously.

Bacterial growth physiology and RNA metabolism

Bacteria are sophisticated systems with high capacity and flexibility to adapt to various environmental conditions. Each prokaryote however possesses a defined metabolic network, which sets its overall metabolic capacity, and therefore the maximal growth rate that can be reached. To achieve optimal growth, bacteria adopt various molecular strategies to optimally adjust gene expression and optimize resource allocation according to the nutrient availability. The resulting physiological changes are often accompanied by changes in the growth rate, and by global regulation of gene expression. The growth-rate-dependent variation of the abundances in the cellular machineries, together with condition-specific regulatory mechanisms, affect RNA metabolism and fate and pose a challenge for rational gene expression reengineering of synthetic circuits. This article is part of a Special Issue entitled: RNA and gene control in bacteria, edited by Dr. M. Guillier and F. Repoila.

Genomic adaptations in information processing underpin trophic strategy in a whole-ecosystem nutrient enrichment experiment

Several universal genomic traits affect trade-offs in the capacity, cost, and efficiency of the biochemical information processing that underpins metabolism and reproduction. We analyzed the role of these traits in mediating the responses of a planktonic microbial community to nutrient enrichment in an oligotrophic, phosphorus-deficient pond in Cuatro Ciénegas, Mexico. This is one of the first whole-ecosystem experiments to involve replicated metagenomic assessment. Mean bacterial genome size, GC content, total number of tRNA genes, total number of rRNA genes, and codon usage bias in ribosomal protein sequences were all higher in the fertilized treatment, as predicted on the basis of the assumption that oligotrophy favors lower information-processing costs whereas copiotrophy favors higher processing rates. Contrasting changes in trait variances also suggested differences between traits in mediating assembly under copiotrophic versus oligotrophic conditions. Trade-offs in information-processing traits are apparently sufficiently pronounced to play a role in community assembly because the major components of metabolism-information, energy, and nutrient requirements-are fine-tuned to an organism's growth and trophic strategy.

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.

Time Series Genome-Centric Analysis Unveils Bacterial Response to Operational Disturbance in Activated Sludge

Understanding ecosystem response to disturbances and identifying the most critical traits for the maintenance of ecosystem functioning are important goals for microbial community ecology. In this study, we used 16S rRNA amplicon sequencing and metagenomics to investigate the assembly of bacterial populations in a full-scale municipal activated sludge wastewater treatment plant over a period of 3 years, including a 9-month period of disturbance characterized by short-term plant shutdowns. Following the reconstruction of 173 metagenome-assembled genomes, we assessed the functional potential, the number of rRNA gene operons, and the in situ growth rate of microorganisms present throughout the time series. Operational disturbances caused a significant decrease in bacteria with a single copy of the rRNA (rrn) operon. Despite moderate differences in resource availability, replication rates were distributed uniformly throughout time, with no differences between disturbed and stable periods. We suggest that the length of the growth lag phase, rather than the growth rate, is the primary driver of selection under disturbed conditions. Thus, the system could maintain its function in the face of disturbance by recruiting bacteria with the capacity to rapidly resume growth under unsteady operating conditions.IMPORTANCE Disturbance is a key determinant of community assembly and dynamics in natural and engineered ecosystems. Microbiome response to disturbance is thought to be influenced by bacterial growth traits and life history strategies. In this time series observational study, the response to disturbance of microbial communities in a full-scale activated sludge wastewater treatment plant was assessed by computing specific cellular traits of genomes retrieved from metagenomes. It was found that the genomes observed in disturbed periods have more copies of the rRNA operon than genomes observed in stable periods, whereas the in situ mean relative growth rates of bacteria present during stable and disturbed periods were indistinguishable. From these intriguing observations, we infer that the length of the lag phase might be a growth trait that affects the microbial response to disturbance. Further exploration of this hypothesis could contribute to better understanding of the adaptive response of microbiomes to unsteady environmental conditions.

Visualization and prediction of CRISPR incidence in microbial trait-space to identify drivers of antiviral immune strategy

Bacteria and archaea are locked in a near-constant battle with their viral pathogens. Despite previous mechanistic characterization of numerous prokaryotic defense strategies, the underlying ecological drivers of different strategies remain largely unknown and predicting which species will take which strategies remains a challenge. Here, we focus on the CRISPR immune strategy and develop a phylogenetically-corrected machine learning approach to build a predictive model of CRISPR incidence using data on over 100 traits across over 2600 species. We discover a strong but hitherto-unknown negative interaction between CRISPR and aerobicity, which we hypothesize may result from interference between CRISPR-associated proteins and non-homologous end-joining DNA repair due to oxidative stress. Our predictive model also quantitatively confirms previous observations of an association between CRISPR and temperature. Finally, we contrast the environmental associations of different CRISPR system types (I, II, III) and restriction modification systems, all of which act as intracellular immune systems.

Wastewater microbial community structure and functional traits change over short timescales

Wastewater contains microorganisms coming from various sources, e.g. feces discharges, soil infiltrations and sewer biofilms and sediments. The primary objective of this work was to determine if end-of-pipe wastewater microbial community structures exhibits short-timescale variation, and assess possible microbial origins. To this end, we measured hourly physicochemical characteristics of wastewater influent for 2 days and analyzed the microbial community at 4-h intervals using 16S rRNA gene amplicon sequencing. Results showed large variations in the microbial community composition at phylum and genus levels, i.e. Proteobacteria ranged from 44 to 63% of the total relative abundance and Arcobacter ranged from 11 to 22%. Diurnal patterns were observed in the alpha-diversity, beta-diversity and the prevalence of several taxa. Wastewater physicochemical characteristics explained 61% of the total microbial community variance by Canonical Correspondence Analysis (CCA), with flow rate being the main explanatory variable exhibiting a clear diurnal profile. Comparison with public databases using closed reference OTUs revealed that only 7.3% of the sequences were shared with human gut microbiota and 21.7% with soil microbiota, the majority being from the sewer biofilms and sediments. The functional trait, weighted average ribosomal RNA operon (rrn) copy number per genome, was found to be relatively high in the wastewater microbiota (average 3.6, soil 2.1, and human gut 2.6) and significantly correlated with flow, inferring active microbial enrichments in the sewer. The prevalence of Methylophilaceae, methanol oxidation genes and denitrification genes were related to high influent methanol and NO₃^- concentration in the influent wastewater. These functional organisms and genes indicate important carbon and nutrient removal related functions in the sewer. Together, the observed temporal patterns of the microbial community and functional traits suggest that high wastewater flow causes greater transport of active sewer microorganisms which are functionally important.

Elongation factor P is required to maintain proteome homeostasis at high growth rate

Elongation factor P (EF-P) is a universally conserved translation factor that alleviates ribosome pausing at polyproline (PPX) motifs by facilitating peptide bond formation. In the absence of EF-P, PPX peptide bond formation can limit translation rate, leading to pleotropic phenotypes including slowed growth, increased antibiotic sensitivity, and loss of virulence. In this study, we observe that many of these phenotypes are dependent on growth rate. Limiting growth rate suppresses a variety of detrimental phenotypes associated with ribosome pausing at PPX motifs in the absence of EF-P. Polysome levels are also similar to wild-type under slow growth conditions, consistent with global changes in ribosome queuing in cells without EF-P when growth rate is decreased. Inversely, under high protein synthesis demands, we observe that Escherichia coli lacking EF-P have reduced fitness. Our data demonstrate that EF-P-mediated relief of ribosome queuing is required to maintain proteome homeostasis under conditions of high translational demands.

Multiple cis-Acting rDNAs Contribute to Nucleoid Separation and Recruit the Bacterial Condensin Smc-ScpAB

Condensins load onto DNA to organize chromosomes. Smc-ScpAB clearly loads onto the parS sites bound by Spo0J, but other loading site(s) must operate independently of parS. In this study, we asked where and how Smc-ScpAB normally selects its loading site. Our results suggest that rDNA is also a loading site. A pull-down assay revealed that Smc-ScpAB preferentially loads onto rDNA in the wild-type cell and even in a Δspo0J mutant but not in a Δsmc mutant. Moreover, we showed that deletion mutants of rDNAs cause a defect in nucleoid separation, and at least two rDNAs near oriC are essential for separation. Full-length rDNA, including promoters, is required for loading and nucleoid separation. A synthetic defect by deletions of both rDNA and spo0J resulted in more aberrant nucleoid separation. We propose that a single-stranded segment of DNA that is exposed at highly transcribed rRNA operons would become a target for Smc-ScpAB loading.

Extrachromosomal Nucleolus-Like Compartmentalization by a Plasmid-Borne Ribosomal RNA Operon and Its Role in Nucleoid Compaction

In the fast-growing Escherichia coli cells, RNA polymerase (RNAP) molecules are concentrated and form foci at clusters of ribosomal RNA (rRNA) operons resembling eukaryotic nucleolus. The bacterial nucleolus-like organization, spatially compartmentalized at the surface of the compact bacterial chromosome (nucleoid), serves as transcription factories for rRNA synthesis and ribosome biogenesis, which influences the organization of the nucleoid. Unlike wild type that has seven rRNA operons in the genome in a mutant that has six (Δ6rrn) rRNA operons deleted in the genome, there are no apparent transcription foci and the nucleoid becomes uncompacted, indicating that formation of RNAP foci requires multiple copies of rRNA operons clustered in space and is critical for nucleoid compaction. It has not been determined, however, whether a multicopy plasmid-borne rRNA operon (prrnB) could substitute the multiple chromosomal rRNA operons for the organization of the bacterial nucleolus-like structure in the mutants of Δ6rrn and Δ7rrn that has all seven rRNA operons deleted in the genome. We hypothesized that extrachromosomal nucleolus-like structures are similarly organized and functional in trans from prrnB in these mutants. In this report, using multicolor images of three-dimensional superresolution Structured Illumination Microscopy (3D-SIM), we determined the distributions of both RNAP and NusB that are a transcription factor involved in rRNA synthesis and ribosome biogenesis, prrnB clustering, and nucleoid structure in these two mutants in response to environmental cues. Our results found that the extrachromosomal nucleolus-like organization tends to be spatially located at the poles of the mutant cells. In addition, formation of RNAP foci at the extrachromosomal nucleolus-like structure condenses the nucleoid, supporting the idea that active transcription at the nucleolus-like organization is a driving force in nucleoid compaction.

Bacterial strategies along nutrient and time gradients, revealed by metagenomic analysis of laboratory microcosms

There is considerable interest in the functional basis of ecological strategies amongst bacteria. We used laboratory microcosms based on culturing of elutant from soil to study the effects of varying initial nutrient concentration, and time succession, on the community metagenome. We found a distinct set of nutrient-related or time-related changes in the functional metagenome. For example, a high nutrient (copiotrophic) strategy was associated with greater abundance of genes related to cell division and cell cycle, while a low nutrient (oligotrophic) strategy had greater abundance of genes related to carbohydrate metabolism and virulence, disease and defense. We also found time-related changes in the functional metagenome, revealing a distinct 'r'-related strategy with greater abundance of genes related to regulation and cell signaling, and a 'K' strategy rich in motility and chemotaxis-related genes. These different gene-based strategies may help to explain how so many bacterial OTUs coexist in nature, and the functional principles dominating natural communities. In terms of diversity, both the OTU richness and the richness of species assignment of functional genes showed linear correlations with functional gene richness, supporting the hypothesis that greater taxonomic diversity is associated with greater functional diversity, with possible implications for ecosystem stability.

On the intrinsic constraint of bacterial growth rate: M. tuberculosis's view of the protein translation capacity

In nature, the maximal growth rates vary widely among different bacteria species. Fast-growing bacteria species such as Escherichia coli can have a shortest generation time of 20 min. Slow-growing bacteria species are perhaps best known for Mycobacterium tuberculosis, a human pathogen with a generation time being no less than 16 h. Despite of the significant progress made on understanding the pathogenesis of M. tuberculosis, we know little on the origin of its intriguingly slow growth. From a global view, the intrinsic constraint of the maximal growth rate of bacteria remains to be a fundamental question in microbiology. In this review, we analyze and discuss this issue from the angle of protein translation capacity, which is the major demand for cell growth. Based on quantitative analysis, we propose four parameters: rRNA chain elongation rate, abundance of RNA polymerase engaged in rRNA synthesis, polypeptide chain elongation rate, and active ribosome fraction, which potentially limit the maximal growth rate of bacteria. We further discuss the relation of these parameters with the growth rate for M. tuberculosis as well as other bacterial species. We highlight future comprehensive investigation of these parameters for different bacteria species to understand how bacteria set their own specific growth rates.

The insect pathogenic bacterium Xenorhabdus innexi has attenuated virulence in multiple insect model hosts yet encodes a potent mosquitocidal toxin

BACKGROUND

Xenorhabdus innexi is a bacterial symbiont of Steinernema scapterisci nematodes, which is a cricket-specialist parasite and together the nematode and bacteria infect and kill crickets. Curiously, X. innexi expresses a potent extracellular mosquitocidal toxin activity in culture supernatants. We sequenced a draft genome of X. innexi and compared it to the genomes of related pathogens to elucidate the nature of specialization.

RESULTS

Using green fluorescent protein-expressing X. innexi we confirm previous reports using culture-dependent techniques that X. innexi colonizes its nematode host at low levels (~3-8 cells per nematode), relative to other Xenorhabdus-Steinernema associations. We found that compared to the well-characterized entomopathogenic nematode symbiont X. nematophila, X. innexi fails to suppress the insect phenoloxidase immune pathway and is attenuated for virulence and reproduction in the Lepidoptera Galleria mellonella and Manduca sexta, as well as the dipteran Drosophila melanogaster. To assess if, compared to other Xenorhabdus spp., X. innexi has a reduced capacity to synthesize virulence determinants, we obtained and analyzed a draft genome sequence. We found no evidence for several hallmarks of Xenorhabdus spp. toxicity, including Tc and Mcf toxins. Similar to other Xenorhabdus genomes, we found numerous loci predicted to encode non-ribosomal peptide/polyketide synthetases. Anti-SMASH predictions of these loci revealed one, related to the fcl locus that encodes fabclavines and zmn locus that encodes zeamines, as a likely candidate to encode the X. innexi mosquitocidal toxin biosynthetic machinery, which we designated Xlt. In support of this hypothesis, two mutants each with an insertion in an Xlt biosynthesis gene cluster lacked the mosquitocidal compound based on HPLC/MS analysis and neither produced toxin to the levels of the wild type parent.

CONCLUSIONS

The X. innexi genome will be a valuable resource in identifying loci encoding new metabolites of interest, but also in future comparative studies of nematode-bacterial symbiosis and niche partitioning among bacterial pathogens.

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.

Ecological strategies and metabolic trade-offs of complex environmental biofilms

Microorganisms aggregated into matrix-enclosed biofilms dominate microbial life in most natural, engineered, and medical systems. Despite this, the ecological adaptations and metabolic trade-offs of the formation of complex biofilms are currently poorly understood. Here, exploring the dynamics of bacterial ribosomal RNA operon (rrn) copy numbers, we unravel the genomic underpinning of the formation and success of stream biofilms that contain hundreds of bacterial taxa. Experimenting with stream biofilms, we found that nascent biofilms in eutrophic systems had reduced lag phases and higher growth rates, and more taxa with higher rrn copy number than biofilms from oligotrophic systems. Based on these growth-related traits, our findings suggest that biofilm succession was dominated by slow-but-efficient bacteria likely with leaky functions, such as the production of extracellular polymeric substances at the cost of rapid growth. Expanding our experimental findings to biofilms from 140 streams, we found that rrn copy number distribution reflects functional trait allocation and ecological strategies of biofilms to be able to thrive in fluctuating environments. These findings suggest that alternative trade-offs dominating over rate-yield trade-offs contribute to the evolutionary success of stream biofilms.

Analyzing natural biofilms containing many types of bacteria yields insights into microbial strategies for success in complex biofilms. The ecological adaptations and metabolic trade-offs involved in the formation of multi-bacterial biofilms in the environment are not well understood. Researchers in Switzerland and Austria, led by Tom Battin and Hannes Peter at the École Polytechnique Fédérale de Lausanne, performed genetic analysis of biofilms sampled from 140 streams. The biofilms contained hundreds of types of bacteria, unlike the mono-bacterial biofilms examined in many laboratory studies. Genetic analysis techniques revealed a diversity of metabolic strategies that allow bacteria to survive within the rich ecology of natural biofilms. Slow-growing but metabolically efficient bacteria that release more extracellular biofilm components thrive better than those adapted for quick growth alone. The findings significantly improve understanding of biofilm ecology in the natural environment.

Expression and Function of Different Guanine-Plus-Cytosine Content 16S rRNA Genes in Haloarcula hispanica at Different Temperatures

The halophilic archaeon Haloarcula hispanica harbors three ribosomal RNA (rRNA) operons (rrnA, rrnB, and rrnC) that contain the 16S rRNA genes rrsA, rrsB, and rrsC, respectively. Although rrsB and rrsC (rrsBC) have almost identical sequences, the rrsA and rrsBC sequences differ by 5.4%, and they differ by 2.5% with respect to guanine-plus-cytosine content (P_GC). The strong correlation between the typical growth temperatures of archaea and P_GC of their 16S rRNA genes suggests that H. hispanica may harbor different 16S rRNA genes having different P_GC to maintain rapid growth in a wide range of temperatures. We therefore performed reverse transcription-coupled quantitative PCR to assess expression levels of rrsA (P_GC, 58.9%) and rrsBC (P_GC, 56.4-56.5%) at various temperatures. The expression ratio of rrsA to rrsBC increased with culture temperature. Mutants with complete deletions of one or two of the three rRNA operons were constructed and their growth rates at different temperatures compared to that of the wild-type. The growth characteristics of the rRNA operon single-mutant strains were indistinguishable from the wild-type. The rRNA operon double-mutant strains maintained the same temperature range as wild-type but displayed reduced growth rates. In particular, the double-mutant strains grew much slower than wild-type at low temperature related to minimum growth temperature of the wild-type. On the other hand, at physiologically high temperatures the wild-type and the double-mutant strain which harbors only rrnA with high-P_GCrrsA grew significantly faster than the double-mutant strain which harbors only rrnC with low-P_GCrrsC. These findings suggest the importance of 16S rRNAs transcribed from rrsA with high-P_GC in maintaining rapid growth of this halophilic archaeon at raised growth temperatures.

Improvement of identification methods for honeybee specific Lactic Acid Bacteria; future approaches

Honeybees face many parasites and pathogens and consequently rely on a diverse set of individual and group-level defenses to prevent disease. The crop microbiota of Apis mellifera, composed of 13 Lactic Acid Bacterial (LAB) species within the genera Lactobacillus and Bifidobacterium, form a beneficial symbiotic relationship with each other and the honeybee to protect their niche and their host. Possibly playing a vital role in honeybee health, it is important that these honeybee specific Lactic Acid Bacterial (hbs-LAB) symbionts can be correctly identified, isolated and cultured, to further investigate their health promoting properties. We have previously reported successful identification to the strain level by culture-dependent methods and we recently sequenced and annotated the genomes of the 13 hbs-LAB. However, the hitherto applied techniques are unfortunately very time consuming, expensive and not ideal when analyzing a vast quantity of samples. In addition, other researchers have constantly failed to identify the 13 hbs-LAB from honeybee samples by using inadequate media and/or molecular techniques based on 16S rRNA gene sequencing with insufficient discriminatory power. The aim of this study was to develop better and more suitable methods for the identification and cultivation of hbs-LAB. We compared currently used bacterial cultivation media and could for the first time demonstrate a significant variation in the hbs-LAB basic requirements for optimal growth. We also present a new bacterial identification approach based on amplicon sequencing of a region of the 16S rRNA gene using the Illumina platform and an error correction software that can be used to successfully differentiate and rapidly identify the 13 hbs-LAB to the strain level.

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

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

A hydrophobic ammonia-oxidizing archaeon of the Nitrosocosmicus clade isolated from coal tar-contaminated sediment

A wide diversity of ammonia-oxidizing archaea (AOA) within the phylum Thaumarchaeota exists and plays a key role in the N cycle in a variety of habitats. In this study, we isolated and characterized an ammonia-oxidizing archaeon, strain MY3, from a coal tar-contaminated sediment. Phylogenetically, strain MY3 falls in clade 'Nitrosocosmicus' of the thaumarchaeotal group I.1b. The cells of strain MY3 are large 'walnut-like' cocci, divide by binary fission along a central cingulum, and form aggregates. Strain MY3 is mesophilic and neutrophilic. An assay of ¹³ C-bicarbonate incorporation into archaeal membrane lipids indicated that strain MY3 is capable of autotrophy. In contrast to some other AOA, TCA cycle intermediates, i.e. pruvate, oxaloacetate and α-ketoglutarate, did not affect the growth rates and yields of strain MY3. The attachment of cells of strain MY3 to XAD-7 hydrophobic beads and to the adsorbent vermiculite demonstrated the potential of strain MY3 to form biofilms. The cell surface was confirmed to be hydrophobic by the extraction of strain MY3 from an aqueous medium with p-xylene. Our finding of a strong potential for surface attachment by strain MY3 may reflect an adaptation to the selective pressures in hydrophobic terrestrial environments.

rRNA Operon Copy Number Can Explain the Distinct Epidemiology of Hospital-Associated Methicillin-Resistant Staphylococcus aureus

The distinct epidemiology of original hospital-associated methicillin-resistant Staphylococcus aureus (HA-MRSA) and early community-associated MRSA (CA-MRSA) is largely unexplained. S. aureus carries either five or six rRNA operon copies. Evidence is provided for a scenario in which MRSA has adapted to the hospital environment by rRNA operon loss (six to five copies) due to antibiotic pressure. Early CA-MRSA, in contrast, results from wild-type methicillin-susceptible S. aureus (MSSA) that acquired mecA without loss of an rRNA operon. Of the HA-MRSA isolates (n = 77), 67.5% had five rRNA operon copies, compared to 23.2% of the CA-MRSA isolates (n = 69) and 7.7% of MSSA isolates (n = 195) (P < 0.001). In addition, 105 MSSA isolates from cystic fibrosis patients were tested, because these patients are repeatedly treated with antibiotics; 32.4% of these isolates had five rRNA operon copies. For all subsets, a correlation between resistance profile and rRNA copy number was found. Furthermore, we showed that in vitro antibiotic pressure may result in rRNA operon copy loss. We also showed that without antibiotic pressure, S. aureus isolates containing six rRNA copies are more fit than isolates with five copies. We conclude that HA-MRSA and cystic fibrosis isolates most likely have adapted to an environment with high antibiotic pressure by the loss of an rRNA operon copy. This loss has facilitated resistance development, which promoted survival in these niches. However, strain fitness decreased, which explains their lack of success in the community. In contrast, CA-MRSA isolates retained six rRNA operon copies, rendering them fitter and thereby able to survive and spread in the community.

Exploiting rRNA operon copy number to investigate bacterial reproductive strategies

The potential for rapid reproduction is a hallmark of microbial life, but microbes in nature must also survive and compete when growth is constrained by resource availability. Successful reproduction requires different strategies when resources are scarce compared to when they are abundant^1,2, but a systematic framework for predicting these reproductive strategies in bacteria has not been available. Here we show that the number of ribosomal RNA operons (rrn) in bacterial genomes predicts two important components of reproduction – growth rate and growth efficiency – which are favored under contrasting regimes of resource availability^3,4. We find that the maximum reproductive rate of bacteria doubles with a doubling of rrn copy number, while the efficiency of carbon use is inversely related to maximal growth rate and rrn copy number. We also identify a feasible explanation for these patterns: the rate and yield of protein synthesis mirror the overall pattern in maximum growth rate and growth efficiency. Furthermore, comparative analysis of genomes from 1,167 bacterial species reveals that rrn copy number predicts traits associated with resource availability, including chemotaxis and genome streamlining. Genome-wide patterns of orthologous gene content covary with rrn copy number, suggesting convergent evolution in response to resource availability. Our findings indicate that basic cellular processes adapt in contrasting ways to long-term differences in resource availability. They also establish a basis for predicting changes in bacterial community composition in response to resource perturbations using rrn copy number measurements⁵ or inferences^6,7.

Contrasting Patterns of rDNA Homogenization within the Zygosaccharomyces rouxii Species Complex

Arrays of repetitive ribosomal DNA (rDNA) sequences are generally expected to evolve as a coherent family, where repeats within such a family are more similar to each other than to orthologs in related species. The continuous homogenization of repeats within individual genomes is a recombination process termed concerted evolution. Here, we investigated the extent and the direction of concerted evolution in 43 yeast strains of the Zygosaccharomyces rouxii species complex (Z. rouxii, Z. sapae, Z. mellis), by analyzing two portions of the 35S rDNA cistron, namely the D1/D2 domains at the 5’ end of the 26S rRNA gene and the segment including the internal transcribed spacers (ITS) 1 and 2 (ITS regions). We demonstrate that intra-genomic rDNA sequence variation is unusually frequent in this clade and that rDNA arrays in single genomes consist of an intermixing of Z. rouxii, Z. sapae and Z. mellis-like sequences, putatively evolved by reticulate evolutionary events that involved repeated hybridization between lineages. The levels and distribution of sequence polymorphisms vary across rDNA repeats in different individuals, reflecting four patterns of rDNA evolution: I) rDNA repeats that are homogeneous within a genome but are chimeras derived from two parental lineages via recombination: Z. rouxii in the ITS region and Z. sapae in the D1/D2 region; II) intra-genomic rDNA repeats that retain polymorphisms only in ITS regions; III) rDNA repeats that vary only in their D1/D2 domains; IV) heterogeneous rDNA arrays that have both polymorphic ITS and D1/D2 regions. We argue that an ongoing process of homogenization following allodiplodization or incomplete lineage sorting gave rise to divergent evolutionary trajectories in different strains, depending upon temporal, structural and functional constraints. We discuss the consequences of these findings for Zygosaccharomyces species delineation and, more in general, for yeast barcoding.

A Mechanistic Model for Cooperative Behavior of Co-transcribing RNA Polymerases

In fast-transcribing prokaryotic genes, such as an rrn gene in Escherichia coli, many RNA polymerases (RNAPs) transcribe the DNA simultaneously. Active elongation of RNAPs is often interrupted by pauses, which has been observed to cause RNAP traffic jams; yet some studies indicate that elongation seems to be faster in the presence of multiple RNAPs than elongation by a single RNAP. We propose that an interaction between RNAPs via the torque produced by RNAP motion on helically twisted DNA can explain this apparent paradox. We have incorporated the torque mechanism into a stochastic model and simulated transcription both with and without torque. Simulation results illustrate that the torque causes shorter pause durations and fewer collisions between polymerases. Our results suggest that the torsional interaction of RNAPs is an important mechanism in maintaining fast transcription times, and that transcription should be viewed as a cooperative group effort by multiple polymerases.