Unprecedented high-resolution view of bacterial operon architecture revealed by RNA sequencing

Pathological R-loops in bacteria from engineered expression of endogenous antisense RNAs whose synthesis is ordinarily terminated by Rho

In many bacteria, the essential factors Rho and NusG mediate termination of synthesis of nascent transcripts (including antisense RNAs) that are not being simultaneously translated. It has been proposed that in Rho's absence toxic RNA-DNA hybrids (R-loops) may be generated from nascent untranslated transcripts, and genome-wide mapping studies in Escherichia coli have identified putative loci of R-loop formation from more than 100 endogenous antisense transcripts that are synthesized only in a Rho-deficient strain. Here we provide evidence that engineered expression in wild-type E. coli of several such individual antisense regions on a plasmid or the chromosome generates R-loops that, in an RNase H-modulated manner, serve to disrupt genome integrity. Rho inhibition was associated with increased prevalence of antisense R-loops also in Xanthomonas oryzae pv. oryzae and Caulobacter crescentus. Our results confirm the essential role of Rho in several bacterial genera for prevention of toxic R-loops from pervasive yet cryptic endogenous antisense transcripts. Engineered antisense R-looped regions may be useful for studies on both site-specific impediments to bacterial chromosomal replication and the mechanisms of their resolution.

Deciphering transcript architectural complexity in bacteria and archaea

RNA transcripts are potential therapeutic targets, yet bacterial transcripts have uncharacterized biodiversity. We developed an algorithm for transcript prediction called tp.py using it to predict transcripts (mRNA and other RNAs) in Escherichia coli K12 and E2348/69 strains (Bacteria:gamma-Proteobacteria), Listeria monocytogenes strains Scott A and RO15 (Bacteria:Firmicute), Pseudomonas aeruginosa strains SG17M and NN2 strains (Bacteria:gamma-Proteobacteria), and Haloferax volcanii (Archaea:Halobacteria). From >5 million E. coli K12 and >3 million E. coli E2348/69 newly generated Oxford Nanopore Technologies direct RNA sequencing reads, 2,487 K12 mRNAs and 1,844 E2348/69 mRNAs were predicted, with the K12 mRNAs containing more than half of the predicted E. coli K12 proteins. While the number of predicted transcripts varied by strain based on the amount of sequence data used, across all strains examined, the predicted average size of the mRNAs was 1.6-1.7 kbp, while the median size of the 5'- and 3'-untranslated regions (UTRs) were 30-90 bp. Given the lack of bacterial and archaeal transcript annotation, most predictions were of novel transcripts, but we also predicted many previously characterized mRNAs and ncRNAs, including post-transcriptionally generated transcripts and small RNAs associated with pathogenesis in the E. coli E2348/69 LEE pathogenicity islands. We predicted small transcripts in the 100-200 bp range as well as >10 kbp transcripts for all strains, with the longest transcript for two of the seven strains being the nuo operon transcript, and for another two strains it was a phage/prophage transcript. This quick, easy, and reproducible method will facilitate the presentation of transcripts, and UTR predictions alongside coding sequences and protein predictions in bacterial genome annotation as important resources for the research community.IMPORTANCEOur understanding of bacterial and archaeal genes and genomes is largely focused on proteins since there have only been limited efforts to describe bacterial/archaeal RNA diversity. This contrasts with studies on the human genome, where transcripts were sequenced prior to the release of the human genome over two decades ago. We developed software for the quick, easy, and reproducible prediction of bacterial and archaeal transcripts from Oxford Nanopore Technologies direct RNA sequencing data. These predictions are urgently needed for more accurate studies examining bacterial/archaeal gene regulation, including regulation of virulence factors, and for the development of novel RNA-based therapeutics and diagnostics to combat bacterial pathogens, like those with extreme antimicrobial resistance.

Analysis of bacterial transcriptome and epitranscriptome using nanopore direct RNA sequencing

Bacterial gene expression is a complex process involving extensive regulatory mechanisms. Along with growing interests in this field, Nanopore Direct RNA Sequencing (DRS) provides a promising platform for rapid and comprehensive characterization of bacterial RNA biology. However, the DRS of bacterial RNA is currently deficient in the yield of mRNA-mapping reads and has yet to be exploited for transcriptome-wide RNA modification mapping. Here, we showed that pre-processing of bacterial total RNA (size selection followed by ribosomal RNA depletion and polyadenylation) guaranteed high throughputs of sequencing data and considerably increased the amount of mRNA reads. This way, complex transcriptome architectures were reconstructed for Escherichia coli and Staphylococcus aureus and extended the boundaries of 225 known E. coli operons and 89 defined S. aureus operons. Utilizing unmodified in vitro-transcribed (IVT) RNA libraries as a negative control, several Nanopore-based computational tools globally detected putative modification sites in the E. coli and S. aureus transcriptomes. Combined with Next-Generation Sequencing-based N6-methyladenosine (m6A) detection methods, 75 high-confidence m6A candidates were identified in the E. coli protein-coding transcripts, while none were detected in S. aureus. Altogether, we demonstrated the potential of Nanopore DRS in systematic and convenient transcriptome and epitranscriptome analysis.

PIPETS: A statistically informed, gene-annotation agnostic analysis method to study bacterial termination using 3'-end sequencing

Background

Over the last decade the drop in short-read sequencing costs have allowed experimental techniques utilizing sequencing to address specific biological questions to proliferate, oftentimes outpacing standardized or effective analysis approaches for the data generated. There are growing amounts of bacterial 3'-end sequencing data, yet there is currently no commonly accepted analysis methodology for this datatype. Most data analysis approaches are somewhat ad hoc and, despite the presence of substantial signal within annotated genes, focus on genomic regions outside the annotated genes (e.g. 3' or 5' UTRs). Furthermore, lack of a systematic approach to analyzing such data makes it impossible to compare conclusions generated by different labs, using different organisms.

Results

We present PIPETS (Poisson Identification of PEaks from Term-Seq data), an R package available on Bioconductor that provides a novel analysis method for 3'-end sequencing data. PIPETS is a statistically informed, gene-annotation agnostic methodology. Across two different datasets PIPETS identified significant 3'-end termination signal across a wider range of annotated genomic contexts than existing analysis approaches, suggesting that existing approaches may miss biologically relevant signal. Furthermore, assessment of the previously called 3'-end positions not captured by PIPETS showed that they were uniformly very low coverage.

Conclusions

PIPETS provides a broadly applicable platform to compare 3'-end sequencing data sets across different organisms. It requires only the 3'-end sequencing data, and is broadly accessible to non-expert users.

Genetic evidence for functional diversification of gram-negative intermembrane phospholipid transporters

The outer membrane of gram-negative bacteria is a barrier to chemical and physical stress. Phospholipid transport between the inner and outer membranes has been an area of intense investigation and, in E. coli K-12, it has recently been shown to be mediated by YhdP, TamB, and YdbH, which are suggested to provide hydrophobic channels for phospholipid diffusion, with YhdP and TamB playing the major roles. However, YhdP and TamB have different phenotypes suggesting distinct functions. It remains unclear whether these functions are related to phospholipid metabolism. We investigated a synthetic cold sensitivity caused by deletion of fadR, a transcriptional regulator controlling fatty acid degradation and unsaturated fatty acid production, and yhdP, but not by ΔtamB ΔfadR or ΔydbH ΔfadR. Deletion of tamB recuses the ΔyhdP ΔfadR cold sensitivity further demonstrating the phenotype is related to functional diversification between these genes. The ΔyhdP ΔfadR strain shows a greater increase in cardiolipin upon transfer to the non-permissive temperature and genetically lowering cardiolipin levels can suppress cold sensitivity. These data also reveal a qualitative difference between cardiolipin synthases in E. coli, as deletion of clsA and clsC suppresses cold sensitivity but deletion of clsB does not. Moreover, increased fatty acid saturation is necessary for cold sensitivity and lowering this level genetically or through supplementation of oleic acid suppresses the cold sensitivity of the ΔyhdP ΔfadR strain. Together, our data clearly demonstrate that the diversification of function between YhdP and TamB is related to phospholipid metabolism. Although indirect regulatory effects are possible, we favor the parsimonious hypothesis that YhdP and TamB have differential phospholipid-substrate transport preferences. Thus, our data provide a potential mechanism for independent control of the phospholipid composition of the inner and outer membranes in response to changing conditions based on regulation of abundance or activity of YhdP and TamB.

Cross-evaluation of E. coli's operon structures via a whole-cell model suggests alternative cellular benefits for low- versus high-expressing operons

Many bacteria use operons to coregulate genes, but it remains unclear how operons benefit bacteria. We integrated E. coli's 788 polycistronic operons and 1,231 transcription units into an existing whole-cell model and found inconsistencies between the proposed operon structures and the RNA-seq read counts that the model was parameterized from. We resolved these inconsistencies through iterative, model-guided corrections to both datasets, including the correction of RNA-seq counts of short genes that were misreported as zero by existing alignment algorithms. The resulting model suggested two main modes by which operons benefit bacteria. For 86% of low-expression operons, adding operons increased the co-expression probabilities of their constituent proteins, whereas for 92% of high-expression operons, adding operons resulted in more stable expression ratios between the proteins. These simulations underscored the need for further experimental work on how operons reduce noise and synchronize both the expression timing and the quantity of constituent genes. A record of this paper's transparent peer review process is included in the supplemental information.

Polyphosphate kinase regulates LPS structure and polymyxin resistance during starvation in E. coli

Polyphosphates (polyP) are chains of inorganic phosphates that can reach over 1,000 residues in length. In Escherichia coli, polyP is produced by the polyP kinase (PPK) and is thought to play a protective role during the response to cellular stress. However, the molecular pathways impacted by PPK activity and polyP accumulation remain poorly characterized. In this work, we used label-free mass spectrometry to study the response of bacteria that cannot produce polyP (Δppk) during starvation to identify novel pathways regulated by PPK. In response to starvation, we found 92 proteins significantly differentially expressed between wild-type and Δppk mutant cells. Wild-type cells were enriched for proteins related to amino acid biosynthesis and transport, while Δppk mutants were enriched for proteins related to translation and ribosome biogenesis, suggesting that without PPK, cells remain inappropriately primed for growth even in the absence of the required building blocks. From our data set, we were particularly interested in Arn and EptA proteins, which were down-regulated in Δppk mutants compared to wild-type controls, because they play a role in lipid A modifications linked to polymyxin resistance. Using western blotting, we confirm differential expression of these and related proteins in K-12 strains and a uropathogenic isolate, and provide evidence that this mis-regulation in Δppk cells stems from a failure to induce the BasRS two-component system during starvation. We also show that Δppk mutants unable to up-regulate Arn and EptA expression lack the respective L-Ara4N and pEtN modifications on lipid A. In line with this observation, loss of ppk restores polymyxin sensitivity in resistant strains carrying a constitutively active basR allele. Overall, we show a new role for PPK in lipid A modification during starvation and provide a rationale for targeting PPK to sensitize bacteria towards polymyxin treatment. We further anticipate that our proteomics work will provide an important resource for researchers interested in the diverse pathways impacted by PPK.

Improved prediction of bacterial CRISPRi guide efficiency from depletion screens through mixed-effect machine learning and data integration

CRISPR interference (CRISPRi) is the leading technique to silence gene expression in bacteria; however, design rules remain poorly defined. We develop a best-in-class prediction algorithm for guide silencing efficiency by systematically investigating factors influencing guide depletion in genome-wide essentiality screens, with the surprising discovery that gene-specific features substantially impact prediction. We develop a mixed-effect random forest regression model that provides better estimates of guide efficiency. We further apply methods from explainable AI to extract interpretable design rules from the model. This study provides a blueprint for predictive models for CRISPR technologies where only indirect measurements of guide activity are available.

Genomic SELEX Screening of Regulatory Targets of Transcription Factors

The genome of Escherichia coli K-12 is transcribed by a single species of RNA polymerase. The selectivity of transcriptional targets is determined via interaction with one of seven species of the sigma subunit and a total of approximately 300 species of transcription factor (TFs). For comprehensive identification of the regulatory targets of these two groups of regulatory proteins on the genome, we developed an in vitro approach, "Genomic SELEX" (gSELEX) screening. Here we describe a detailed protocol of the gSELEX screening system, which uses purified regulatory proteins and fragments of genomic DNA from E. coli. Moreover, we describe methods and examples of results using cell-free synthetic proteins.

The timing of transcription of RpoS-dependent genes varies across multiple stresses in Escherichia coli K-12

IMPORTANCE

Bacteria adapt to changing environments by altering the transcription of their genes. Specific proteins can regulate these changes. This study explored how a single protein called RpoS controls how many genes change expression during adaptation to three stresses. We found that: (i) RpoS is responsible for activating different genes in different stresses; (ii) that during a stress, the timing of gene activation depends on the what stress it is; and (iii) that how much RpoS a gene needs in order to be activated can predict when that gene will be activated during the stress of stationary phase.

Evidence that OLE RNA is a component of a major stress-responsive ribonucleoprotein particle in extremophilic bacteria

OLE RNA is a ~600-nucleotide noncoding RNA present in many Gram-positive bacteria that thrive mostly in extreme environments, including elevated temperature, salt, and pH conditions. The precise biochemical functions of this highly conserved RNA remain unknown, but it forms a ribonucleoprotein (RNP) complex that localizes to cell membranes. Genetic disruption of the RNA or its essential protein partners causes reduced cell growth under various stress conditions. These phenotypes include sensitivity to short-chain alcohols, cold intolerance, reduced growth on sub-optimal carbon sources, and intolerance of even modest concentrations of Mg2+ . Thus, many bacterial species appear to employ OLE RNA as a component of an intricate RNP apparatus to monitor fundamental cellular processes and make physiological and metabolic adaptations. Herein we hypothesize that the OLE RNP complex is functionally equivalent to the eukaryotic TOR complexes, which integrate signals from various diverse pathways to coordinate processes central to cell growth, replication, and survival.

RNA polymerase sliding on DNA can couple the transcription of nearby bacterial operons

DNA transcription initiates after an RNA polymerase (RNAP) molecule binds to the promoter of a gene. In bacteria, the canonical picture is that RNAP comes from the cytoplasmic pool of freely diffusing RNAP molecules. Recent experiments suggest the possible existence of a separate pool of polymerases, competent for initiation, which freely slide on the DNA after having terminated one round of transcription. Promoter-dependent transcription reinitiation from this pool of posttermination RNAP may lead to coupled initiation at nearby operons, but it is unclear whether this can occur over the distance and timescales needed for it to function widely on a bacterial genome in vivo. Here, we mathematically model the hypothesized reinitiation mechanism as a diffusion-to-capture process and compute the distances over which significant interoperon coupling can occur and the time required. These quantities depend on molecular association and dissociation rate constants between DNA, RNAP, and the transcription initiation factor σ70; we measure these rate constants using single-molecule experiments in vitro. Our combined theory/experimental results demonstrate that efficient coupling can occur at physiologically relevant σ70 concentrations and on timescales appropriate for transcript synthesis. Coupling is efficient over terminator-promoter distances up to ∼1,000 bp, which includes the majority of terminator-promoter nearest neighbor pairs in the Escherichia coli genome. The results suggest a generalized mechanism that couples the transcription of nearby operons and breaks the paradigm that each binding of RNAP to DNA can produce at most one messenger RNA.

Improved RNA stability estimation indicates that transcriptional interference is frequent in diverse bacteria

We used stochastic simulations and experimental data from E. coli, K. aerogenes, Synechococcus PCC 7002 and Synechocystis PCC 6803 to provide evidence that transcriptional interference via the collision mechanism is likely a prevalent mechanism for bacterial gene regulation. Rifampicin time-series data can be used to globally monitor and quantify collision between sense and antisense transcription-complexes. Our findings also highlight that transcriptional events, such as differential RNA decay, partial termination, and internal transcriptional start sites often deviate from gene annotations. Consequently, within a single gene annotation, there exist transcript segments with varying half-lives and transcriptional properties. To address these complexities, we introduce 'rifi', an R-package that analyzes transcriptomic data from rifampicin time series. 'rifi' employs a dynamic programming-based segmentation approach to identify individual transcripts, enabling accurate assessment of RNA stability and detection of diverse transcriptional events.

OmpC-Dependent Bile Tolerance Contributes to E. coli Colonization of the Mammalian Intestine

Escherichia coli persistently colonizes the mammalian intestine by mechanisms that are not fully understood. Previously, we found when streptomycin-treated mice were fed E. coli MG1655, the intestine selected for envZ missense mutants that outcompeted the wild type. The better-colonizing envZ mutants had a higher level of OmpC and reduced OmpF. This suggested the EnvZ/OmpR two-component system and outer membrane proteins play a role in colonization. In this study, we show that wild-type E. coli MG1655 outcompetes an envZ-ompR knockout mutant. Moreover, ompA and ompC knockout mutants are outcompeted by the wild type, while an ompF knockout mutant colonizes better than the wild type. Outer membrane protein gels show the ompF mutant overproduces OmpC. An ompC mutant is more sensitive to bile salts than the wild type and ompF mutant. The ompC mutant initiates colonization slowly because it is sensitive to physiological concentrations of bile salts in the intestine. Overexpression of ompC under the control of a constitutive promoter confers a colonization advantage only when ompF is deleted. These results indicate that fine-tuning of OmpC and OmpF levels is needed to maximize competitive fitness in the intestine. RNA sequencing reveals the EnvZ/OmpR two-component system is active in the intestine: ompC is upregulated and ompF is downregulated. While other factors could also contribute to the advantage provided by OmpC, we provide evidence that OmpC is important for E. coli to colonize the intestine because its smaller pore size excludes bile salts or other unknown toxic substances, while OmpF is deleterious because its larger pore size allows bile salts or other unknown toxic substances to enter the periplasm. IMPORTANCE Every mammalian intestine is colonized with Escherichia coli. Although E. coli is one of the most studied model organisms, how it colonizes the intestine is not fully understood. Here, we investigated the role of the EnvZ/OmpR two-component system and outer membrane proteins in colonization of the mouse intestine by E. coli. We report that an ompC mutant is a poor colonizer, while an ompF mutant, which overproduces OmpC, outcompetes the wild type. OmpF has a larger pore size that allows toxic bile salts or other toxic compounds into the cell and is deleterious for colonization of the intestine. OmpC has a smaller pore size and excludes bile salts. Our findings provide insights into why E. coli fine-tunes the levels of OmpC and OmpF during colonization via the EnvZ/OmpR two-component system.

Glacier-Fed Stream Biofilms Harbor Diverse Resistomes and Biosynthetic Gene Clusters

Antimicrobial resistance (AMR) is a universal phenomenon the origins of which lay in natural ecological interactions such as competition within niches, within and between micro- to higher-order organisms. To study these phenomena, it is crucial to examine the origins of AMR in pristine environments, i.e., limited anthropogenic influences. In this context, epilithic biofilms residing in glacier-fed streams (GFSs) are an excellent model system to study diverse, intra- and inter-domain, ecological crosstalk. We assessed the resistomes of epilithic biofilms from GFSs across the Southern Alps (New Zealand) and the Caucasus (Russia) and observed that both bacteria and eukaryotes encoded twenty-nine distinct AMR categories. Of these, beta-lactam, aminoglycoside, and multidrug resistance were both abundant and taxonomically distributed in most of the bacterial and eukaryotic phyla. AMR-encoding phyla included Bacteroidota and Proteobacteria among the bacteria, alongside Ochrophyta (algae) among the eukaryotes. Additionally, biosynthetic gene clusters (BGCs) involved in the production of antibacterial compounds were identified across all phyla in the epilithic biofilms. Furthermore, we found that several bacterial genera (Flavobacterium, Polaromonas, Superphylum Patescibacteria) encode both atimicrobial resistance genes (ARGs) and BGCs within close proximity of each other, demonstrating their capacity to simultaneously influence and compete within the microbial community. Our findings help unravel how naturally occurring BGCs and AMR contribute to the epilithic biofilms mode of life in GFSs. Additionally, we report that eukaryotes may serve as AMR reservoirs owing to their potential for encoding ARGs. Importantly, these observations may be generalizable and potentially extended to other environments that may be more or less impacted by human activity. IMPORTANCE Antimicrobial resistance is an omnipresent phenomenon in the anthropogenically influenced ecosystems. However, its role in shaping microbial community dynamics in pristine environments is relatively unknown. Using metagenomics, we report the presence of antimicrobial resistance genes and their associated pathways in epilithic biofilms within glacier-fed streams. Importantly, we observe biosynthetic gene clusters associated with antimicrobial resistance in both pro- and eukaryotes in these biofilms. Understanding the role of resistance in the context of this pristine environment and complex biodiversity may shed light on previously uncharacterized mechanisms of cross-domain interactions.

Evolved bacterial resistance to the chemotherapy gemcitabine modulates its efficacy in co-cultured cancer cells

Drug metabolism by the microbiome can influence anticancer treatment success. We previously suggested that chemotherapies with antimicrobial activity can select for adaptations in bacterial drug metabolism that can inadvertently influence the host's chemoresistance. We demonstrated that evolved resistance against fluoropyrimidine chemotherapy lowered its efficacy in worms feeding on drug-evolved bacteria (Rosener et al., 2020). Here, we examine a model system that captures local interactions that can occur in the tumor microenvironment. Gammaproteobacteria-colonizing pancreatic tumors can degrade the nucleoside-analog chemotherapy gemcitabine and, in doing so, can increase the tumor's chemoresistance. Using a genetic screen in Escherichia coli, we mapped all loss-of-function mutations conferring gemcitabine resistance. Surprisingly, we infer that one third of top resistance mutations increase or decrease bacterial drug breakdown and therefore can either lower or raise the gemcitabine load in the local environment. Experiments in three E. coli strains revealed that evolved adaptation converged to inactivation of the nucleoside permease NupC, an adaptation that increased the drug burden on co-cultured cancer cells. The two studies provide complementary insights on the potential impact of microbiome adaptation to chemotherapy by showing that bacteria-drug interactions can have local and systemic influence on drug activity.

RNA Preparation and RNA-Seq Bioinformatics for Comparative Transcriptomics

The principal transcriptome analysis is the determination of differentially expressed genes across experimental conditions. For this, the next-generation sequencing of RNA (RNA-seq) has several advantages over other techniques, such as the capability of detecting all the transcripts in one assay over RT-qPCR, such as its higher accuracy and broader dynamic range over microarrays or the ability to detect novel transcripts, including non-coding RNA molecules, at nucleotide-level resolution over both techniques. Despite these advantages, many microbiology laboratories have not yet applied RNA-seq analyses to their investigations. The high cost of the equipment for next-generation sequencing is no longer an issue since this intermediate part of the analysis can be provided by commercial or central services. Here, we detail a protocol for the first part of the analysis, the RNA extraction and an introductory protocol to the bioinformatics analysis of the sequencing data to generate the differential expression results.

INRI-seq enables global cell-free analysis of translation initiation and off-target effects of antisense inhibitors

Ribosome profiling (Ribo-seq) is a powerful method for the transcriptome-wide assessment of protein synthesis rates and the study of translational control mechanisms. Yet, Ribo-seq also has limitations. These include difficulties with the analysis of translation-modulating molecules such as antibiotics, which are often toxic or challenging to deliver into living cells. Here, we have developed in vitro Ribo-seq (INRI-seq), a cell-free method to analyze the translational landscape of a fully customizable synthetic transcriptome. Using Escherichia coli as an example, we show how INRI-seq can be used to analyze the translation initiation sites of a transcriptome of interest. We also study the global impact of direct translation inhibition by antisense peptide nucleic acid (PNA) to analyze PNA off-target effects. Overall, INRI-seq presents a scalable, sensitive method to study translation initiation in a transcriptome-wide manner without the potentially confounding effects of extracting ribosomes from living cells.

Regulatory Interplay between RNase III and Antisense RNAs in E. coli: the Case of AsflhD and FlhD, Component of the Master Regulator of Motility

In order to respond to ever-changing environmental cues, bacteria display resilient regulatory mechanisms controlling gene expression. At the post-transcriptional level, this is achieved by a combination of RNA-binding proteins, such as ribonucleases (RNases), and regulatory RNAs, including antisense RNAs (asRNAs). Bound to their complementary mRNA, asRNAs are primary targets for the double-strand-specific endoribonuclease, RNase III. Taking advantage of our own and previously published transcriptomic data sets obtained in strains inactivated for RNase III, we selected several candidate asRNAs and confirmed the existence of RNase III-sensitive asRNAs for crp, ompR, phoP, and flhD genes, all encoding global regulators of gene expression in Escherichia coli. Using FlhD, a component of the master regulator of motility (FlhD₄C₂), as our model, we demonstrate that the asRNA AsflhD, transcribed from the coding sequence of flhD, is involved in the fine-tuning of flhD expression and thus participates in the control of motility.

Mapping the Complex Transcriptional Landscape of the Phytopathogenic Bacterium Dickeya dadantii

Dickeya dadantii is a phytopathogenic bacterium that causes soft rot in a wide range of plant hosts worldwide and a model organism for studying virulence gene regulation. The present study provides a comprehensive and annotated transcriptomic map of D. dadantii obtained by a computational method combining five independent transcriptomic data sets: (i) paired-end RNA sequencing (RNA-seq) data for a precise reconstruction of the RNA landscape; (ii) DNA microarray data providing transcriptional responses to a broad variety of environmental conditions; (iii) long-read Nanopore native RNA-seq data for isoform-level transcriptome validation and determination of transcription termination sites; (iv) differential RNA sequencing (dRNA-seq) data for the precise mapping of transcription start sites; (v) in planta DNA microarray data for a comparison of gene expression profiles between in vitro experiments and the early stages of plant infection. Our results show that transcription units sometimes coincide with predicted operons but are generally longer, most of them comprising internal promoters and terminators that generate alternative transcripts of variable gene composition. We characterize the occurrence of transcriptional read-through at terminators, which might play a basal regulation role and explain the extent of transcription beyond the scale of operons. We finally highlight the presence of noncontiguous operons and excludons in the D. dadantii genome, novel genomic arrangements that might contribute to the basal coordination of transcription. The highlighted transcriptional organization may allow D. dadantii to finely adjust its gene expression program for a rapid adaptation to fast-changing environments. IMPORTANCE This is the first transcriptomic map of a Dickeya species. It may therefore significantly contribute to further progress in the field of phytopathogenicity. It is also one of the first reported applications of long-read Nanopore native RNA-seq in prokaryotes. Our findings yield insights into basal rules of coordination of transcription that might be valid for other bacteria and may raise interest in the field of microbiology in general. In particular, we demonstrate that gene expression is coordinated at the scale of transcription units rather than operons, which are larger functional genomic units capable of generating transcripts with variable gene composition for a fine-tuning of gene expression in response to environmental changes. In line with recent studies, our findings indicate that the canonical operon model is insufficient to explain the complexity of bacterial transcriptomes.

RegulonDB 11.0: Comprehensive high-throughput datasets on transcriptional regulation in Escherichia coli K-12

Genomics has set the basis for a variety of methodologies that produce high-throughput datasets identifying the different players that define gene regulation, particularly regulation of transcription initiation and operon organization. These datasets are available in public repositories, such as the Gene Expression Omnibus, or ArrayExpress. However, accessing and navigating such a wealth of data is not straightforward. No resource currently exists that offers all available high and low-throughput data on transcriptional regulation in Escherichia coli K-12 to easily use both as whole datasets, or as individual interactions and regulatory elements. RegulonDB (https://regulondb.ccg.unam.mx) began gathering high-throughput dataset collections in 2009, starting with transcription start sites, then adding ChIP-seq and gSELEX in 2012, with up to 99 different experimental high-throughput datasets available in 2019. In this paper we present a radical upgrade to more than 2000 high-throughput datasets, processed to facilitate their comparison, introducing up-to-date collections of transcription termination sites, transcription units, as well as transcription factor binding interactions derived from ChIP-seq, ChIP-exo, gSELEX and DAP-seq experiments, besides expression profiles derived from RNA-seq experiments. For ChIP-seq experiments we offer both the data as presented by the authors, as well as data uniformly processed in-house, enhancing their comparability, as well as the traceability of the methods and reproducibility of the results. Furthermore, we have expanded the tools available for browsing and visualization across and within datasets. We include comparisons against previously existing knowledge in RegulonDB from classic experiments, a nucleotide-resolution genome viewer, and an interface that enables users to browse datasets by querying their metadata. A particular effort was made to automatically extract detailed experimental growth conditions by implementing an assisted curation strategy applying Natural language processing and machine learning. We provide summaries with the total number of interactions found in each experiment, as well as tools to identify common results among different experiments. This is a long-awaited resource to make use of such wealth of knowledge and advance our understanding of the biology of the model bacterium E. coli K-12.

Genome-wide promoter assembly in E. coli measured at single-base resolution

When detected at single-base-pair resolution, the genome-wide location, occupancy level, and structural organization of DNA-binding proteins provide mechanistic insights into genome regulation. Here we use ChIP-exo to provide a near-base-pair resolution view of the epigenomic organization of the Escherichia coli transcription machinery and nucleoid structural proteins at the time when cells are growing exponentially and upon rapid reprogramming (acute heat shock). We examined the site specificity of three sigma factors (RpoD/σ70, RpoH/σ32, and RpoN/σ54), RNA polymerase (RNAP or RpoA, -B, -C), and two nucleoid proteins (Fis and IHF). We suggest that DNA shape at the flanks of cognate motifs helps drive site specificity. We find that although RNAP and sigma factors occupy active cognate promoters, RpoH and RpoN can occupy quiescent promoters without the presence of RNAP. Thus, promoter-bound sigma factors can be triggered to recruit RNAP by a mechanism that is distinct from an obligatory cycle of free sigma binding RNAP followed by promoter binding. These findings add new dimensions to how sigma factors achieve promoter specificity through DNA sequence and shape, and further define mechanistic steps in regulated genome-wide assembly of RNAP at promoters in E. coli.

Relationship between the Chromosome Structural Dynamics and Gene Expression—A Chicken and Egg Dilemma?

Prokaryotic transcription was extensively studied over the last half-century. A great deal of data has been accumulated regarding the control of gene expression by transcription factors regulating their target genes by binding at specific DNA sites. However, there is a significant gap between the mechanistic description of transcriptional control obtained from in vitro biochemical studies and the complexity of transcriptional regulation in the context of the living cell. Indeed, recent studies provide ample evidence for additional levels of complexity pertaining to the regulation of transcription in vivo, such as, for example, the role of the subcellular localization and spatial organization of different molecular components involved in the transcriptional control and, especially, the role of chromosome configurational dynamics. The question as to how the chromosome is dynamically reorganized under the changing environmental conditions and how this reorganization is related to gene expression is still far from being clear. In this article, we focus on the relationships between the chromosome structural dynamics and modulation of gene expression during bacterial adaptation. We argue that spatial organization of the bacterial chromosome is of central importance in the adaptation of gene expression to changing environmental conditions and vice versa, that gene expression affects chromosome dynamics.

Direct measurements of mRNA translation kinetics in living cells

Ribosome mediated mRNA translation is central to life. The cycle of translation, however, has been characterized mostly using reconstituted systems, with only few techniques applicable for studies in the living cell. Here we describe a live-cell ribosome-labeling method, which allows us to characterize the whole processes of finding and translating an mRNA, using single-molecule tracking techniques. We find that more than 90% of both bacterial ribosomal subunits are engaged in translation at any particular time, and that the 30S and 50S ribosomal subunits spend the same average time bound to an mRNA, revealing that 30S re-initiation on poly-cistronic mRNAs is not prevalent in E. coli. Instead, our results are best explained by substantial 70S re-initiation of translation of poly-cistronic mRNAs, which is further corroborated by experiments with translation initiation inhibitors. Finally, we find that a variety of previously described orthogonal ribosomes, with altered anti-Shine-Dalgarno sequences, show significant binding to endogenous mRNAs.

Transcriptional organization and regulation of the Pseudomonas putida flagellar system

A single region of the Pseudomonas putida genome, designated the flagellar cluster, includes 59 genes potentially involved in the biogenesis and function of the flagellar system. Here, we combine bioinformatics and in vivo gene expression analyses to clarify the transcriptional organization and regulation of the flagellar genes in the cluster. We have identified 11 flagellar operons and characterized 22 primary and internal promoter regions. Our results indicate that synthesis of the flagellar apparatus and core chemotaxis machinery is regulated by a three-tier cascade in which fleQ is a Class I gene, standing at the top of the transcriptional hierarchy. FleQ- and σ⁵⁴ -dependent Class II genes encode most components of the flagellar structure, part of the chemotaxis machinery and multiple regulatory elements, including the flagellar σ factor FliA. FliA activation of Class III genes enables synthesis of the filament, one stator complex and completion of the chemotaxis apparatus. Accessory regulatory proteins and an intricate operon architecture add complexity to the regulation by providing feedback and feed-forward loops to the main circuit. Because of the high conservation of the gene arrangement and promoter motifs, we believe that the regulatory circuit presented here may also apply to other environmental pseudomonads.

Assembling stable syntrophic Escherichia coli communities by comprehensively identifying beneficiaries of secreted goods

Metabolic cross-feeding frequently underlies mutualistic relationships in natural microbial communities and is often exploited to assemble synthetic microbial consortia. We systematically identified all single-gene knockouts suitable for imposing cross-feeding in Escherichia coli and used this information to assemble syntrophic communities. Most strains benefiting from shared goods were dysfunctional in biosynthesis of amino acids, nucleotides, and vitamins or mutants in central carbon metabolism. We tested cross-feeding potency in 1,444 strain pairs and mapped the interaction network between all functional groups of mutants. This network revealed that auxotrophs for vitamins are optimal cooperators. Lastly, we monitored how assemblies composed of dozens of auxotrophs change over time and observed that they rapidly and repeatedly coalesced to seven strain consortia composed primarily from vitamin auxotrophs. The composition of emerging consortia suggests that they were stabilized by multiple cross-feeding interactions. We conclude that vitamins are ideal shared goods since they optimize consortium growth while still imposing member co-dependence.

The Complex Transcriptional Landscape of Magnetosome Gene Clusters in Magnetospirillum gryphiswaldense

Magnetosomes are complex membrane organelles synthesized by magnetotactic bacteria (MTB) for navigation in the Earth’s magnetic field. In the alphaproteobacterium Magnetospirillum gryphiswaldense, all steps of magnetosome formation are tightly controlled by >30 specific genes arranged in several gene clusters. However, the transcriptional organization of the magnetosome gene clusters has remained poorly understood. Here, by applying Cappable-seq and whole-transcriptome shotgun RNA sequencing, we show that mamGFDCop and feoAB1op are transcribed as single transcriptional units, whereas multiple transcription start sites (TSS) are present in mms6op, mamXYop, and the long (>16 kb) mamABop. Using a bioluminescence reporter assay and promoter knockouts, we demonstrate that most of the identified TSS originate from biologically meaningful promoters which mediate production of multiple transcripts and are functionally relevant for proper magnetosome biosynthesis. In addition, we identified a strong promoter in a large intergenic region within mamXYop, which likely drives transcription of a noncoding RNA important for gene expression in this operon. In summary, our data suggest a more complex transcriptional architecture of the magnetosome operons than previously recognized, which is largely conserved in other magnetotactic Magnetospirillum species and, thus, is likely fundamental for magnetosome biosynthesis in these organisms.

IMPORTANCE Magnetosomes have emerged as a model system to study prokaryotic organelles and a source of biocompatible magnetic nanoparticles for various biomedical applications. However, the lack of knowledge about the transcriptional organization of magnetosome gene clusters has severely impeded the engineering, manipulation, and transfer of this highly complex biosynthetic pathway into other organisms. Here, we provide a high-resolution image of the previously unappreciated transcriptional landscape of the magnetosome operons. Our findings are important for further unraveling the complex genetic framework of magnetosome biosynthesis. In addition, they will facilitate the rational reengineering of magnetic bacteria for improved bioproduction of tunable magnetic nanoparticles, as well as transplantation of magnetosome biosynthesis into foreign hosts by synthetic biology approaches. Overall, our study exemplifies how a genetically complex pathway is orchestrated at the transcriptional level to ensure the balanced expression of the numerous constituents required for the proper assembly of one of the most intricate prokaryotic organelles.

Cis-Regulatory Logic Produces Gene-Expression Noise Describing Phenotypic Heterogeneity in Bacteria

Gene transcriptional process is random. It occurs in bursts and follows single-molecular kinetics. Intermittent bursts are measured based on their frequency and size. They influence temporal fluctuations in the abundance of total mRNA and proteins by generating distinct transcriptional variations referred to as “noise”. Noisy expression induces uncertainty because the association between transcriptional variation and the extent of gene expression fluctuation is ambiguous. The promoter architecture and remote interference of different cis-regulatory elements are the crucial determinants of noise, which is reflected in phenotypic heterogeneity. An alternative perspective considers that cellular parameters dictating genome-wide transcriptional kinetics follow a universal pattern. Research on noise and systematic perturbations of promoter sequences reinforces that both gene-specific and genome-wide regulation occur across species ranging from bacteria and yeast to animal cells. Thus, deciphering gene-expression noise is essential across different genomics applications. Amidst the mounting conflict, it is imperative to reconsider the scope, progression, and rational construction of diversified viewpoints underlying the origin of the noise. Here, we have established an indication connecting noise, gene expression variations, and bacterial phenotypic variability. This review will enhance the understanding of gene-expression noise in various scientific contexts and applications.

Enhancement of lysine biosynthesis confers high-temperature stress tolerance to Escherichia coli cells

Abstract

Lysine, a nutritionally important amino acid, is involved in adaptation and tolerance to environmental stresses in various organisms. Previous studies reported that lysine accumulation occurs in response to stress and that lysine supplementation enhances stress tolerance; however, the effect of lysine biosynthesis enhancement on stress tolerance has yet to be elucidated. In this study, we confirmed that lysine supplementation to the culture medium increased intracellular lysine content and improved cell growth of Escherichia coli at high temperature (42.5 °C). Lysine-overproducing strains were then isolated from the lysine analogue S-adenosylmethionine-resistant mutants by conventional mutagenesis and exhibited higher tolerance to high-temperature stress than the wild-type strain. We identified novel amino acid substitutions Gly474Asp and Cys554Tyr on ThrA, a bifunctional aspartate kinase/homoserine dehydrogenase (AK/HSDH), in the lysine-overproducing mutants. Interestingly, the Gly474Asp and Cys554Tyr variants of ThrA induced lysine accumulation and conferred high-temperature stress tolerance to E. coli cells. Enzymatic analysis revealed that the Gly474Asp substitution in ThrA reduced HSDH activity, suggesting that the intracellular level of aspartate semialdehyde, which is a substrate for HSDH and an intermediate for lysine biosynthesis, is elevated by the loss of HSDH activity and converted to lysine in E. coli. The present study demonstrated that both lysine supplementation and lysine biosynthesis enhancement improved the high-temperature stress tolerance of E. coli cells. Our findings suggest that lysine-overproducing strains have the potential as stress-tolerant microorganisms and can be applied to robust host cells for microbial production of useful compounds.

Key points

• Lysine supplementation improved the growth of E. coli cells at high temperature.

• The G474D and C554Y variant ThrA increased lysine productivity in E. coli cells.

• The G474D substitution in ThrA reduced homoserine dehydrogenase activity.

• E. coli cells that overproduce lysine exhibited high-temperature stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-021-11519-0.

Clusters of hairpins induce intrinsic transcription termination in bacteria

Intrinsic transcription termination (ITT) sites are currently identified by locating single and double-adjacent RNA hairpins downstream of the stop codon. ITTs for a limited number of genes/operons in only a few bacterial genomes are currently known. This lack of coverage is a lacuna in the existing ITT inference methods. We have studied the inter-operon regions of 13 genomes covering all major phyla in bacteria, for which good quality public RNA-seq data exist. We identify ITT sites in 87% of cases by predicting hairpin(s) and validate against 81% of cases for which the RNA-seq derived sites could be calculated. We identify 72% of these sites correctly, with 98% of them located ≤ 80 bases downstream of the stop codon. The predicted hairpins form a cluster (when present < 15 bases) in two-thirds of the cases, the remaining being single hairpins. The largest number of clusters is formed by two hairpins, and the occurrence decreases exponentially with an increasing number of hairpins in the cluster. Our study reveals that hairpins form an effective ITT unit when they act in concert in a cluster. Their pervasiveness along with single hairpin terminators corroborates a wider utilization of ITT mechanisms for transcription control across bacteria.

Quantitative Control for Stoichiometric Protein Synthesis

Bacterial protein synthesis rates have evolved to maintain preferred stoichiometries at striking precision, from the components of protein complexes to constituents of entire pathways. Setting relative protein production rates to be well within a factor of two requires concerted tuning of transcription, RNA turnover, and translation, allowing many potential regulatory strategies to achieve the preferred output. The last decade has seen a greatly expanded capacity for precise interrogation of each step of the central dogma genome-wide. Here, we summarize how these technologies have shaped the current understanding of diverse bacterial regulatory architectures underpinning stoichiometric protein synthesis. We focus on the emerging expanded view of bacterial operons, which encode diverse primary and secondary mRNA structures for tuning protein stoichiometry. Emphasis is placed on how quantitative tuning is achieved. We discuss the challenges and open questions in the application of quantitative, genome-wide methodologies to the problem of precise protein production. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Predicting synchronized gene coexpression patterns from fibration symmetries in gene regulatory networks in bacteria

BACKGROUND

Gene regulatory networks coordinate the expression of genes across physiological states and ensure a synchronized expression of genes in cellular subsystems, critical for the coherent functioning of cells. Here we address the question whether it is possible to predict gene synchronization from network structure alone. We have recently shown that synchronized gene expression can be predicted from symmetries in the gene regulatory networks described by the concept of symmetry fibrations. We showed that symmetry fibrations partition the genes into groups called fibers based on the symmetries of their 'input trees', the set of paths in the network through which signals can reach a gene. In idealized dynamic gene expression models, all genes in a fiber are perfectly synchronized, while less idealized models-with gene input functions differencing between genes-predict symmetry breaking and desynchronization.

RESULTS

To study the functional role of gene fibers and to test whether some of the fiber-induced coexpression remains in reality, we analyze gene fibrations for the gene regulatory networks of E. coli and B. subtilis and confront them with expression data. We find approximate gene coexpression patterns consistent with symmetry fibrations with idealized gene expression dynamics. This shows that network structure alone provides useful information about gene synchronization, and suggest that gene input functions within fibers may be further streamlined by evolutionary pressures to realize a coexpression of genes.

CONCLUSIONS

Thus, gene fibrations provide a sound conceptual tool to describe tunable coexpression induced by network topology and shaped by mechanistic details of gene expression.

Hierarchy of transcription factor network in Escherichia coli K-12: H-NS-mediated silencing and Anti-silencing by global regulators

Transcriptional regulation for genome expression determines growth and adaptation of single-cell bacteria that are directly exposed to environment. The transcriptional apparatus in Escherichia coli K-12 is composed of RNA polymerase core enzyme and two groups of its regulatory proteins, seven species of promoter-recognition subunit sigma and about 300 species of transcription factors. The identification of regulatory targets for all these regulatory proteins is critical toward understanding the genome regulation as a whole. For this purpose, we performed a systematic search in vitro of the whole set of binding sites for each factor by gSELEX system. This review summarizes the accumulated knowledge of regulatory targets for more than 150 TFs from E. coli K-12. Overall TFs could be classified into four families: nucleoid-associated bifunctional TFs; global regulators; local regulators; and single-target regulators, in which the regulatory functions remain uncharacterized for the nucleoid-associated TFs. Here we overview the regulatory targets of two nucleoid-associated TFs, H-NS and its paralog StpA, both together playing the silencing role of a set of non-essential genes. Participation of LeuO and other global regulators have been indicated for the anti-silencing. Finally, we propose the hierarchy of TF network as a key framework of the bacterial genome regulation.

Triple RNA-Seq characterizes aphid gene expression in response to infection with unequally virulent strains of the endosymbiont Hamiltonella defensa

Background

Secondary endosymbionts of aphids provide benefits to their hosts, but also impose costs such as reduced lifespan and reproductive output. The aphid Aphis fabae is host to different strains of the secondary endosymbiont Hamiltonella defensa, which encode different putative toxins. These strains have very different phenotypes: They reach different densities in the host, and the costs and benefits (protection against parasitoid wasps) they confer to the host vary strongly.

Results

We used RNA-Seq to generate hypotheses on why four of these strains inflict such different costs to A. fabae. We found different H. defensa strains to cause strain-specific changes in aphid gene expression, but little effect of H. defensa on gene expression of the primary endosymbiont, Buchnera aphidicola. The highly costly and over-replicating H. defensa strain H85 was associated with strongly reduced aphid expression of hemocytin, a marker of hemocytes in Drosophila. The closely related strain H15 was associated with downregulation of ubiquitin-related modifier 1, which is related to nutrient-sensing and oxidative stress in other organisms. Strain H402 was associated with strong differential regulation of a set of hypothetical proteins, the majority of which were only differentially regulated in presence of H402.

Conclusions

Overall, our results suggest that costs of different strains of H. defensa are likely caused by different mechanisms, and that these costs are imposed by interacting with the host rather than the host’s obligatory endosymbiont B. aphidicola.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07742-8.

Different csrA Expression Levels in C versus K-12 E. coli Strains Affect Biofilm Formation and Impact the Regulatory Mechanism Presided by the CsrB and CsrC Small RNAs

Escherichia coli C is a strong biofilm producer in comparison to E. coli K-12 laboratory strains due to higher expression of the pgaABCD operon encoding the enzymes for the biosynthesis of the extracellular polysaccharide poly-β-1,6-N-acetylglucosamine (PNAG). The pgaABCD operon is negatively regulated at the post-transcriptional level by two factors, namely CsrA, a conserved RNA-binding protein controlling multiple pathways, and the RNA exonuclease polynucleotide phosphorylase (PNPase). In this work, we investigated the molecular bases of different PNAG production in C-1a and MG1655 strains taken as representative of E. coli C and K-12 strains, respectively. We found that pgaABCD operon expression is significantly lower in MG1655 than in C-1a; consistently, CsrA protein levels were much higher in MG1655. In contrast, we show that the negative effect exerted by PNPase on pgaABCD expression is much stronger in C-1a than in MG1655. The amount of CsrA and of the small RNAs CsrB, CsrC, and McaS sRNAs regulating CsrA activity is dramatically different in the two strains, whereas PNPase level is similar. Finally, the compensatory regulation acting between CsrB and CsrC in MG1655 does not occur in E. coli C. Our results suggest that PNPase preserves CsrA-dependent regulation by indirectly modulating csrA expression.

Calcium-Regulated Protein CarP Responds to Multiple Host Signals and Mediates Regulation of Pseudomonas aeruginosa Virulence by Calcium

Pseudomonas aeruginosa is an opportunistic pathogen causing life-threatening infections. Previously, we showed that elevated calcium (Ca2+) levels increase the production of virulence factors in P. aeruginosa In an effort to characterize the Ca2+ regulatory network, we identified a Ca2+-regulated β-propeller protein, CarP, and showed that expression of the encoding gene is controlled by the Ca2+-regulated two-component system CarSR. Here, by using a Galleria melonella model, we showed that CarP plays a role in regulating P. aeruginosa virulence. By using transcriptome sequencing (RNA-Seq), reverse transcription (RT)-PCR, quantitative RT-PCR (RT-qPCR), and promoter fusions, we determined that carP is transcribed into at least two transcripts and regulated by several bacterial and host factors. The transcription of carP is elevated in response to Ca2+ in P. aeruginosa cystic fibrosis isolates and PAO1 laboratory strain. Elevated Fe2+ also induces carP The simultaneous addition of Ca2+ and Fe2+ increased the carP promoter activity synergistically, which requires the presence of CarR. In silico analysis of the intergenic sequence upstream of carP predicted recognition sites of RhlR/LasR, OxyR, and LexA, suggesting regulation by quorum sensing (QS) and oxidative stress. In agreement, the carP promoter was activated in response to stationary-phase PAO1 supernatant and required the presence of elevated Ca2+ and CarR but remained silent in the triple mutant lacking rhlI, lasI, and pqsA synthases. We also showed that carP transcription is regulated by oxidative stress and that CarP contributes to P. aeruginosa Ca2+-dependent H2O2 tolerance. The multifactorial regulation of carP suggests that CarP plays an important role in P. aeruginosa adaptations to host environments.IMPORTANCEP. aeruginosa is a human pathogen causing life-threatening infections. It is particularly notorious for its ability to adapt to diverse environments within the host. Understanding the signals and the signaling pathways enabling P. aeruginosa adaptation is imperative for developing effective therapies to treat infections caused by this organism. One host signal of particular importance is calcium. Previously, we identified a component of the P. aeruginosa calcium-signaling network, CarP, whose expression is induced by elevated levels of calcium. Here, we show that carP plays an important role in P. aeruginosa virulence and is upregulated in P. aeruginosa strains isolated from sputa of patients with cystic fibrosis. We also identified several bacterial and host factors that regulate the transcription of carP Such multifactorial regulation highlights the interconnectedness between regulatory circuits and, together with the pleotropic effect of CarP on virulence, suggests the importance of this protein in P. aeruginosa adaptations to the host.

The role of two-component regulatory systems in environmental sensing and virulence in Salmonella

Adaptation to environments with constant fluctuations imposes challenges that are only overcome with sophisticated strategies that allow bacteria to perceive environmental conditions and develop an appropriate response. The gastrointestinal environment is a complex ecosystem that is home to trillions of microorganisms. Termed microbiota, this microbial ensemble plays important roles in host health and provides colonization resistance against pathogens, although pathogens have evolved strategies to circumvent this barrier. Among the strategies used by bacteria to monitor their environment, one of the most important are the sensing and signalling machineries of two-component systems (TCSs), which play relevant roles in the behaviour of all bacteria. Salmonella enterica is no exception, and here we present our current understanding of how this important human pathogen uses TCSs as an integral part of its lifestyle. We describe important aspects of these systems, such as the stimuli and responses involved, the processes regulated, and their roles in virulence. We also dissect the genomic organization of histidine kinases and response regulators, as well as the input and output domains for each TCS. Lastly, we explore how these systems may be promising targets for the development of antivirulence therapeutics to combat antibiotic-resistant infections.

Rapid transcriptional and metabolic adaptation of intestinal microbes to host immune activation

The gut microbiota produces metabolites that regulate host immunity, thereby impacting disease resistance and susceptibility. The extent to which commensal bacteria reciprocally respond to immune activation, however, remains largely unexplored. Herein, we colonized mice with four anaerobic symbionts and show that acute immune responses result in dramatic transcriptional reprogramming of these commensals with minimal changes in their relative abundance. Transcriptomic changes include induction of stress-response mediators and downregulation of carbohydrate-degrading factors such as polysaccharide utilization loci (PULs). Flagellin and anti-CD3 antibody, two distinct immune stimuli, induced similar transcriptional profiles, suggesting that commensal bacteria detect common effectors or activate shared pathways when facing different host responses. Immune activation altered the intestinal metabolome within 6 hours, decreasing luminal short-chain fatty acid and increasing aromatic metabolite concentrations. Thus, intestinal bacteria, prior to detectable shifts in community composition, respond to acute host immune activation by rapidly changing gene transcription and immunomodulatory metabolite production.

Local genetic context shapes the function of a gene regulatory network

Gene expression levels are influenced by multiple coexisting molecular mechanisms. Some of these interactions such as those of transcription factors and promoters have been studied extensively. However, predicting phenotypes of gene regulatory networks (GRNs) remains a major challenge. Here, we use a well-defined synthetic GRN to study in Escherichia coli how network phenotypes depend on local genetic context, i.e. the genetic neighborhood of a transcription factor and its relative position. We show that one GRN with fixed topology can display not only quantitatively but also qualitatively different phenotypes, depending solely on the local genetic context of its components. Transcriptional read-through is the main molecular mechanism that places one transcriptional unit (TU) within two separate regulons without the need for complex regulatory sequences. We propose that relative order of individual TUs, with its potential for combinatorial complexity, plays an important role in shaping phenotypes of GRNs.

BSGatlas: a unified Bacillus subtilis genome and transcriptome annotation atlas with enhanced information access

A large part of our current understanding of gene regulation in Gram-positive bacteria is based on Bacillus subtilis, as it is one of the most well studied bacterial model systems. The rapid growth in data concerning its molecular and genomic biology is distributed across multiple annotation resources. Consequently, the interpretation of data from further B. subtilis experiments becomes increasingly challenging in both low- and large-scale analyses. Additionally, B. subtilis annotation of structured RNA and non-coding RNA (ncRNA), as well as the operon structure, is still lagging behind the annotation of the coding sequences. To address these challenges, we created the B. subtilis genome atlas, BSGatlas, which integrates and unifies multiple existing annotation resources. Compared to any of the individual resources, the BSGatlas contains twice as many ncRNAs, while improving the positional annotation for 70 % of the ncRNAs. Furthermore, we combined known transcription start and termination sites with lists of known co-transcribed gene sets to create a comprehensive transcript map. The combination with transcription start/termination site annotations resulted in 717 new sets of co-transcribed genes and 5335 untranslated regions (UTRs). In comparison to existing resources, the number of 5' and 3' UTRs increased nearly fivefold, and the number of internal UTRs doubled. The transcript map is organized in 2266 operons, which provides transcriptional annotation for 92 % of all genes in the genome compared to the at most 82 % by previous resources. We predicted an off-target-aware genome-wide library of CRISPR-Cas9 guide RNAs, which we also linked to polycistronic operons. We provide the BSGatlas in multiple forms: as a website (https://rth.dk/resources/bsgatlas/), an annotation hub for display in the UCSC genome browser, supplementary tables and standardized GFF3 format, which can be used in large scale -omics studies. By complementing existing resources, the BSGatlas supports analyses of the B. subtilis genome and its molecular biology with respect to not only non-coding genes but also genome-wide transcriptional relationships of all genes.

Prediction and analysis of metagenomic operons via MetaRon: a pipeline for prediction of Metagenome and whole-genome opeRons

Background

Efficient regulation of bacterial genes in response to the environmental stimulus results in unique gene clusters known as operons. Lack of complete operonic reference and functional information makes the prediction of metagenomic operons a challenging task; thus, opening new perspectives on the interpretation of the host-microbe interactions.

Results

In this work, we identified whole-genome and metagenomic operons via MetaRon (Metagenome and whole-genome opeRon prediction pipeline). MetaRon identifies operons without any experimental or functional information. MetaRon was implemented on datasets with different levels of complexity and information. Starting from its application on whole-genome to simulated mixture of three whole-genomes (E. coli MG1655, Mycobacterium tuberculosis H37Rv and Bacillus subtilis str. 16), E. coli c20 draft genome extracted from chicken gut and finally on 145 whole-metagenome data samples from human gut. MetaRon consistently achieved high operon prediction sensitivity, specificity and accuracy across E. coli whole-genome (97.8, 94.1 and 92.4%), simulated genome (93.7, 75.5 and 88.1%) and E. coli c20 (87, 91 and 88%,), respectively. Finally, we identified 1,232,407 unique operons from 145 paired-end human gut metagenome samples. We also report strong association of type 2 diabetes with Maltose phosphorylase (K00691), 3-deoxy-D-glycero-D-galacto-nononate 9-phosphate synthase (K21279) and an uncharacterized protein (K07101).

Conclusion

With MetaRon, we were able to remove two notable limitations of existing whole-genome operon prediction methods: (1) generalizability (ability to predict operons in unrelated bacterial genomes), and (2) whole-genome and metagenomic data management. We also demonstrate the use of operons as a subset to represent the trends of secondary metabolites in whole-metagenome data and the role of secondary metabolites in the occurrence of disease condition. Using operonic data from metagenome to study secondary metabolic trends will significantly reduce the data volume to more precise data. Furthermore, the identification of metabolic pathways associated with the occurrence of type 2 diabetes (T2D) also presents another dimension of analyzing the human gut metagenome. Presumably, this study is the first organized effort to predict metagenomic operons and perform a detailed analysis in association with a disease, in this case type 2 diabetes. The application of MetaRon to metagenomic data at diverse scale will be beneficial to understand the gene regulation and therapeutic metagenomics.

Characterizing Genetic Parts and Devices Using RNA Sequencing

Synthetic genetic circuits are composed of many parts that must interact and function together to produce a desired pattern of gene expression. A challenge when assembling circuits is that genetic parts often behave differently within a circuit, potentially impacting the desired functionality. Existing debugging methods based on fluorescent reporter proteins allow for only a few internal states to be monitored simultaneously, making diagnosis of the root cause impossible for large systems. Here, we present a tool called the Genetic Analyzer which uses RNA sequencing data to simultaneously characterize all transcriptional parts (e.g., promoters and terminators) and devices (e.g., sensors and logic gates) in complex genetic circuits. This provides a complete picture of the inner workings of a genetic circuit enabling faults to be easily identified and fixed. We construct a complete workflow to coordinate the execution of the various data processing and analysis steps and explain the options available when adapting these for the characterization of new systems.

Protein dosage of the lldPRD operon is correlated with RNase E-dependent mRNA processing

The ability of Escherichia coli to grow on L-lactate as a sole carbon source depends on the expression of the lldPRD operon. A striking feature of this operon is that the transcriptional regulator (LldR) encoding gene is located between the permease (LldP) and the dehydrogenase (LldD) encoding genes. In this study we report that dosage of the LldP, LldR, and LldD proteins is not modulated on the transcriptional level. Instead, modulation of protein dosage is primarily correlated with RNase E-dependent mRNA processing events that take place within the lldR mRNA, leading to the immediate inactivation of lldR, to differential segmental stabilities of the resulting cleavage products, and to differences in the translation efficiencies of the three cistrons. A model for the processing events controlling the molar quantities of the proteins in the lldPRD operon is presented and discussed.ImportanceAdjustment of gene expression is critical for proper cell function. For the case of polycistronic transcripts, posttranscriptional regulatory mechanisms can be used to fine-tune the expression of individual cistrons. Here, we elucidate how protein dosage of the Escherichia coli lldPRD operon, which presents the paradox of having the gene encoding a regulator protein located between genes that code for a permease and an enzyme, is regulated. Our results demonstrate that the key event in this regulatory mechanism involves the RNase E-dependent cleavage of the primary lldPRD transcript at internal site(s) located within the lldR cistron, resulting in a drastic decrease of intact lldR mRNA, to differential segmental stabilities of the resulting cleavage products, and to differences in the translation efficiencies of the three cistrons.

Evolved bacterial resistance against fluoropyrimidines can lower chemotherapy impact in the Caenorhabditis elegans host

Metabolism of host-targeted drugs by the microbiome can substantially impact host treatment success. However, since many host-targeted drugs inadvertently hamper microbiome growth, repeated drug administration can lead to microbiome evolutionary adaptation. We tested if evolved bacterial resistance against host-targeted drugs alters their drug metabolism and impacts host treatment success. We used a model system of Caenorhabditis elegans, its bacterial diet, and two fluoropyrimidine chemotherapies. Genetic screens revealed that most of loss-of-function resistance mutations in Escherichia coli also reduced drug toxicity in the host. We found that resistance rapidly emerged in E. coli under natural selection and converged to a handful of resistance mechanisms. Surprisingly, we discovered that nutrient availability during bacterial evolution dictated the dietary effect on the host - only bacteria evolving in nutrient-poor media reduced host drug toxicity. Our work suggests that bacteria can rapidly adapt to host-targeted drugs and by doing so may also impact the host.

Advances in bacterial transcriptome understanding: From overlapping transcription to the excludon concept

In the last decade, the implementation of high‐throughput methods for RNA profiling has uncovered that a large part of the bacterial genome is transcribed well beyond the boundaries of known genes. Therefore, the transcriptional space of a gene very often invades the space of a neighbouring gene, creating large regions of overlapping transcription. The biological significance of these findings was initially regarded with scepticism. However, mounting evidence suggests that overlapping transcription between neighbouring genes conforms to regulatory purposes and provides new strategies for coordinating bacterial gene expression. In this MicroReview, considering the discoveries made in a pioneering transcriptome analysis performed on Listeria monocytogenes as a starting point, we discuss the progress in understanding the biological meaning of overlapping transcription that has given rise to the excludon concept. We also discuss new conditional transcriptional termination events that create antisense RNAs depending on the metabolite concentrations and new genomic arrangements, known as noncontiguous operons, which contain an interspersed gene that is transcribed in the opposite direction to the rest of the operon.

Redefining fundamental concepts of transcription initiation in bacteria

Despite enormous progress in understanding the fundamentals of bacterial gene regulation, our knowledge remains limited when compared with the number of bacterial genomes and regulatory systems to be discovered. Derived from a small number of initial studies, classic definitions for concepts of gene regulation have evolved as the number of characterized promoters has increased. Together with discoveries made using new technologies, this knowledge has led to revised generalizations and principles. In this Expert Recommendation, we suggest precise, updated definitions that support a logical, consistent conceptual framework of bacterial gene regulation, focusing on transcription initiation. The resulting concepts can be formalized by ontologies for computational modelling, laying the foundation for improved bioinformatics tools, knowledge-based resources and scientific communication. Thus, this work will help researchers construct better predictive models, with different formalisms, that will be useful in engineering, synthetic biology, microbiology and genetics.

Ribosome recycling is not critical for translational coupling in Escherichia coli

We used ribosome profiling to characterize the biological role of ribosome recycling factor (RRF) in Escherichia coli. As expected, RRF depletion leads to enrichment of post-termination 70S complexes in 3′-UTRs. We also observe that elongating ribosomes are unable to complete translation because they are blocked by non-recycled ribosomes at stop codons. Previous studies have suggested a role for recycling in translational coupling within operons; if a ribosome remains bound to an mRNA after termination, it may re-initiate downstream. We found, however, that RRF depletion did not significantly affect coupling efficiency in reporter assays or in ribosome density genome-wide. These findings argue that re-initiation is not a major mechanism of translational coupling in E. coli. Finally, RRF depletion has dramatic effects on the activity of ribosome rescue factors tmRNA and ArfA. Our results provide a global view of the effects of the loss of ribosome recycling on protein synthesis in E. coli.

Deciphering the regulatory genome of Escherichia coli, one hundred promoters at a time

Advances in DNA sequencing have revolutionized our ability to read genomes. However, even in the most well-studied of organisms, the bacterium Escherichia coli, for ≈65% of promoters we remain ignorant of their regulation. Until we crack this regulatory Rosetta Stone, efforts to read and write genomes will remain haphazard. We introduce a new method, Reg-Seq, that links massively parallel reporter assays with mass spectrometry to produce a base pair resolution dissection of more than a E. coli promoters in 12 growth conditions. We demonstrate that the method recapitulates known regulatory information. Then, we examine regulatory architectures for more than 80 promoters which previously had no known regulatory information. In many cases, we also identify which transcription factors mediate their regulation. This method clears a path for highly multiplexed investigations of the regulatory genome of model organisms, with the potential of moving to an array of microbes of ecological and medical relevance.

Whole-genome Identification of Transcriptional Start Sites by Differential RNA-seq in Bacteria

Gene transcription in bacteria often starts some nucleotides upstream of the start codon. Identifying the specific Transcriptional Start Site (TSS) is essential for genetic manipulation, as in many cases upstream of the start codon there are sequence elements that are involved in gene expression regulation. Taken into account the classical gene structure, we are able to identify two kinds of transcriptional start site: primary and secondary. A primary transcriptional start site is located some nucleotides upstream of the translational start site, while a secondary transcriptional start site is located within the gene encoding sequence. Here, we present a step by step protocol for genome-wide transcriptional start sites determination by differential RNA-sequencing (dRNA-seq) using the enteric pathogen Shigella flexneri serotype 5a strain M90T as model. However, this method can be employed in any other bacterial species of choice. In the first steps, total RNA is purified from bacterial cultures using the hot phenol method. Ribosomal RNA (rRNA) is specifically depleted via hybridization probes using a commercial kit. A 5'-monophosphate-dependent exonuclease (TEX)-treated RNA library enriched in primary transcripts is then prepared for comparison with a library that has not undergone TEX-treatment, followed by ligation of an RNA linker adaptor of known sequence allowing the determination of TSS with single nucleotide precision. Finally, the RNA is processed for Illumina sequencing library preparation and sequenced as purchased service. TSS are identified by in-house bioinformatic analysis. Our protocol is cost-effective as it minimizes the use of commercial kits and employs freely available software.

Transcriptional Organization of the Stability Module of Broad-Host-Range Plasmid RA3, from the IncU Group

The broad-host-range (BHR) conjugative plasmids have developed diverse adaptive mechanisms defining the range of their promiscuity. The BHR conjugative RA3 plasmid, the archetype of the IncU group, can transfer between, replicate in, and be maintained in representatives of Alpha-, Beta-, and Gammaproteobacteria Its stability module encompasses ten open reading frames (ORFs) apparently organized into five operons, all transcribed in the same direction from several strong promoters that are tightly regulated either by autorepressors or by global plasmid-encoded regulators. In this paper, we demonstrate that owing to an efficient RNA polymerase (RNAP) read-through, the transcription from the first promoter, orf02p, may continue through the whole module. Moreover, an analysis of mRNA produced from the wild-type (WT) stability module and its deletion variants deprived of particular internal transcription initiation sites reveals that in fact each operon may be transcribed from any upstream promoter, giving rise to multicistronic transcripts of variable length and creating an additional level of gene expression control by transcript dosage adjustment. The gene expression patterns differ among various hosts, indicating that promoter recognition, regulation, and the RNAP read-through mechanisms are modulated in a species-specific manner.IMPORTANCE The efficiently disseminating conjugative or mobilizable BHR plasmids play key roles in the horizontal spread of genetic information between closely related and phylogenetically distant species, which can be harmful from the medical, veterinary, or industrial point of view. Understanding the mechanisms determining the plasmid's ability to function in diverse hosts is essential to help limit the spread of undesirable plasmid-encoded traits, e.g., antibiotic resistance. The range of a plasmid's promiscuity depends on the adaptations of its transfer, replication, and stability functions to the various hosts. IncU plasmids, with the archetype plasmid RA3, are considered to constitute a reservoir of antibiotic resistance genes in aquatic environments; however, the molecular mechanisms determining their adaptability to a broad range of hosts are rather poorly characterized. Here, we present the transcriptional organization of the stability module and show that the gene transcript dosage effect is an important determinant of the stable maintenance of RA3 in different hosts.