Extensive reshaping of bacterial operons by programmed mRNA decay

Resolving spatiotemporal dynamics in bacterial multicellular populations: approaches and challenges

SUMMARYThe development of multicellularity represents a key evolutionary transition that is crucial for the emergence of complex life forms. Although multicellularity has traditionally been studied in eukaryotes, it originates in prokaryotes. Coordinated aggregation of individual cells within the confines of a colony results in emerging, higher-level functions that benefit the population as a whole. During colony differentiation, an almost infinite number of ecological and physiological population-forming forces are at work, creating complex, intricate colony structures with divergent functions. Understanding the assembly and dynamics of such populations requires resolving individual cells or cell groups within such macroscopic structures. Addressing how each cell contributes to the collective action requires pushing the resolution boundaries of key technologies that will be presented in this review. In particular, single-cell techniques provide powerful tools for studying bacterial multicellularity with unprecedented spatial and temporal resolution. These advancements include novel microscopic techniques, mass spectrometry imaging, flow cytometry, spatial transcriptomics, single-bacteria RNA sequencing, and the integration of spatiotemporal transcriptomics with microscopy, alongside advanced microfluidic cultivation systems. This review encourages exploring the synergistic potential of the new technologies in the study of bacterial multicellularity, with a particular focus on individuals in differentiated bacterial biofilms (colonies). It highlights how resolving population structures at the single-cell level and understanding their respective functions can elucidate the overarching functions of bacterial multicellular populations.

Predicting synthetic mRNA stability using massively parallel kinetic measurements, biophysical modeling, and machine learning

mRNA degradation is a central process that affects all gene expression levels, though it remains challenging to predict the stability of a mRNA from its sequence, due to the many coupled interactions that control degradation rate. Here, we carried out massively parallel kinetic decay measurements on over 50,000 bacterial mRNAs, using a learn-by-design approach to develop and validate a predictive sequence-to-function model of mRNA stability. mRNAs were designed to systematically vary translation rates, secondary structures, sequence compositions, G-quadruplexes, i-motifs, and RppH activity, resulting in mRNA half-lives from about 20 seconds to 20 minutes. We combined biophysical models and machine learning to develop steady-state and kinetic decay models of mRNA stability with high accuracy and generalizability, utilizing transcription rate models to identify mRNA isoforms and translation rate models to calculate ribosome protection. Overall, the developed model quantifies the key interactions that collectively control mRNA stability in bacterial operons and predicts how changing mRNA sequence alters mRNA stability, which is important when studying and engineering bacterial genetic systems.

The expression of integron arrays is shaped by the translation rate of cassettes

Integrons are key elements in the rise and spread of multidrug resistance in Gram-negative bacteria. These genetic platforms capture cassettes containing promoterless genes and stockpile them in arrays of variable length. In the current integron model, expression of cassettes is granted by the Pc promoter in the platform and is assumed to decrease as a function of its distance. Here we explored this model using a large collection of 136 antibiotic resistance cassettes and show the effect of distance is in fact negligible. Instead, cassettes have a strong impact in the expression of downstream genes because their translation rate affects the stability of the whole polycistronic mRNA molecule. Hence, cassettes with reduced translation rates decrease the expression and resistance phenotype of cassettes downstream. Our data puts forward an integron model in which expression is contingent on the translation of cassettes upstream, rather than on the distance to the Pc.

Metabolic engineering for single-cell protein production from renewable feedstocks and its applications

Proteins are indispensable for maintaining a healthy diet and performing crucial functions in a multitude of physiological processes. The growth of the global population and the emergence of environmental concerns have significantly increased the demand for protein-rich foods such as meat and dairy products, exerting considerable pressure on global food supplies. Single-cell proteins (SCP) have emerged as a promising alternative source, characterized by their high protein content and essential amino acids, lipids, carbohydrates, nucleic acids, inorganic salts, vitamins, and trace elements. SCP offers several advantages over the traditional animal and plant proteins. These include shorter production cycles, the use of diverse raw material sources, high energy efficiency, and minimal environmental impact. This review is primarily concerned with the microbial species employed in SCP production, utilization of non-food renewable materials as a source of feedstock, and application of rational and non-rational metabolic engineering strategies to increase SCP biomass and protein content. Moreover, the current applications, production shortages, and safety concerns associated with SCP are discussed.

Exploring the transcription start sites and other genomic features facilitates the accurate identification and annotation of small RNAs across multiple stress conditions in Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB) is a pathogen that is known for its ability to persist in harsh environments and cause chronic infections. Understanding the regulatory networks of MTB is crucial for developing effective treatments. Small regulatory RNAs (sRNAs) play important roles in gene expression regulation in all kingdoms of life, and their classification based solely on genomic location can be imprecise due to the computational-based prediction of protein-coding genes in bacteria, which often neglects segments of mRNA such as 5'UTRs, 3'UTRs, and intercistronic regions of operons. To address this issue, our study simultaneously discovered genomic features such as TSSs, UTRs, and operons together with sRNAs in the M. tuberculosis H37Rv strain (ATCC 27294) across multiple stress conditions. Our analysis identified 1,376 sRNA candidates and 8,173 TSSs in MTB, providing valuable insights into its complex regulatory landscape. TSS mapping enabled us to classify these sRNAs into more specific categories, including promoter-associated sRNAs, 5'UTR-derived sRNAs, 3'UTR-derived sRNAs, true intergenic sRNAs, and antisense sRNAs. Three of these sRNA candidates were experimentally validated using 3'-RACE-PCR: predictedRNA_0240, predictedRNA_0325, and predictedRNA_0578. Future characterization and validation are necessary to fully elucidate the functions and roles of these sRNAs in MTB. Our study is the first to simultaneously unravel TSSs and sRNAs in MTB and demonstrate that the identification of other genomic features, such as TSSs, UTRs, and operons, allows for more accurate and specific classification of sRNAs.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

Heterogeneous lineage-specific arginine deiminase expression within dental microbiome species

Arginine catabolism by the bacterial arginine deiminase system (ADS) has anticariogenic properties through the production of ammonia, which modulates the pH of the oral environment. Given the potential protective capacity of the ADS pathway, the exploitation of ADS-competent oral microbes through pre- or probiotic applications is a promising therapeutic target to prevent tooth decay. To date, most investigations of the ADS in the oral cavity and its relation to caries have focused on indirect measures of activity or on specific bacterial groups, yet the pervasiveness and rate of expression of the ADS operon in diverse mixed microbial communities in oral health and disease remain an open question. Here, we use a multivariate approach, combining ultra-deep metatranscriptomic sequencing with paired metataxonomic and in vitro citrulline quantification to characterize the microbial community and ADS operon expression in healthy and late-stage cavitated teeth. While ADS activity is higher in healthy teeth, we identify multiple bacterial lineages with upregulated ADS activity on cavitated teeth that are distinct from those found on healthy teeth using both reference-based mapping and de novo assembly methods. Our dual metataxonomic and metatranscriptomic approach demonstrates the importance of species abundance for gene expression data interpretation and that patterns of differential expression can be skewed by low-abundance groups. Finally, we identify several potential candidate probiotic bacterial lineages within species that may be useful therapeutic targets for the prevention of tooth decay and propose that the development of a strain-specific, mixed-microbial probiotic may be a beneficial approach given the heterogeneity of taxa identified here across health groups.

IMPORTANCE

Tooth decay is the most common preventable chronic disease, affecting more than two billion people globally. The development of caries on teeth is primarily a consequence of acid production by cariogenic bacteria that inhabit the plaque microbiome. Other bacterial strains in the oral cavity may suppress or prevent tooth decay by producing ammonia as a byproduct of the arginine deiminase metabolic pathway, increasing the pH of the plaque biofilm. While the benefits of arginine metabolism on oral health have been extensively documented in specific bacterial groups, the prevalence and consistency of arginine deiminase system (ADS) activity among oral bacteria in a community context remain an open question. In the current study, we use a multi-omics approach to document the pervasiveness of the expression of the ADS operon in both health and disease to better understand the conditions in which ADS activity may prevent tooth decay.

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.

Selective RNA Processing and Stabilization are Multi-Layer and Stoichiometric Regulators of Gene Expression in Escherichia coli

Selective RNA processing and stabilization (SRPS) facilitates the differential expression of multiple genes in polycistronic operons. However, how the coordinated actions of SRPS-related enzymes affect stoichiometric regulation remains unclear. In the present study, the first genome-wide targetome analysis is reported of these enzymes in Escherichia coli, at a single-nucleotide resolution. A strictly linear relationship is observed between the RNA pyrophosphohydrolase processing ratio and scores assigned to the first three nucleotides of the primary transcript. Stem-loops associated with PNPase targetomes exhibit a folding free energy that is negatively correlated with the termination ratio of PNPase at the 3' end. More than one-tenth of the RNase E processing sites in the 5'-untranslated regions（UTR） form different stem-loops that affect ribosome-binding and translation efficiency. The effectiveness of the SRPS elements is validated using a dual-fluorescence reporter system. The findings highlight a multi-layer and quantitative regulatory method for optimizing the stoichiometric expression of genes in bacteria and promoting the application of SRPS in synthetic biology.

Methodologies for bacterial ribonuclease characterization using RNA-seq

Bacteria adjust gene expression at the post-transcriptional level through an intricate network of small regulatory RNAs and RNA-binding proteins, including ribonucleases (RNases). RNases play an essential role in RNA metabolism, regulating RNA stability, decay, and activation. These enzymes exhibit species-specific effects on gene expression, bacterial physiology, and different strategies of target recognition. Recent advances in high-throughput RNA sequencing (RNA-seq) approaches have provided a better understanding of the roles and modes of action of bacterial RNases. Global studies aiming to identify direct targets of RNases have highlighted the diversity of RNase activity and RNA-based mechanisms of gene expression regulation. Here, we review recent RNA-seq approaches used to study bacterial RNases, with a focus on the methods for identifying direct RNase targets.

Complex Regulation of Gamma-Hemolysin Expression Impacts Staphylococcus aureus Virulence

Staphylococcus aureus gamma-hemolysin CB (HlgCB) is a core-genome-encoded pore-forming toxin that targets the C5a receptor, similar to the phage-encoded Panton-Valentine leucocidin (PVL). Absolute quantification by mass spectrometry of HlgCB in 39 community-acquired pneumonia (CAP) isolates showed considerable variations in the HlgC and HlgB yields between isolates. Moreover, although HlgC and HlgB are encoded on a single operon, their levels were dissociated in 10% of the clinical strains studied. To decipher the molecular basis for the variation in hlgCB expression and protein production among strains, different regulation levels were analyzed in representative clinical isolates and reference strains. Both the HlgCB level and the HlgC/HlgB ratio were found to depend on hlgC promoter activity and mRNA processing and translation. Strikingly, only one single nucleotide polymorphism (SNP) in the 5' untranslated region (UTR) of hlgCB mRNA strongly impaired hlgC translation in the USA300 strain, leading to a strong decrease in the level of HlgC but not in HlgB. Finally, we found that high levels of HlgCB synthesis led to mortality in a rabbit model of pneumonia, correlated with the implication of the role of HlgCB in severe S. aureus CAP. Taken together, this work illustrates the complexity of virulence factor expression in clinical strains and demonstrates a butterfly effect where subtle genomic variations have a major impact on phenotype and virulence. IMPORTANCE S. aureus virulence in pneumonia results in its ability to produce several virulence factors, including the leucocidin PVL. Here, we demonstrate that HlgCB, another leucocidin, which targets the same receptors as PVL, highly contributes to S. aureus virulence in pvl-negative strains. In addition, considerable variations in HlgCB quantities are observed among clinical isolates from patients with CAP. Biomolecular analyses have revealed that a few SNPs in the promoter sequences and only one SNP in the 5' UTR of hlgCB mRNA induce the differential expression of hlgCB, drastically impacting hlgC mRNA translation. This work illustrates the subtlety of regulatory mechanisms in bacteria, especially the sometimes major effects on phenotypes of single nucleotide variation in noncoding regions.

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.

Envelope Stress Activates Expression of the Twin Arginine Translocation (Tat) System in Salmonella

The twin arginine translocation system (Tat) is a protein export system that is conserved in bacteria, archaea, and plants. In Gram-negative bacteria, it is required for the export of folded proteins from the cytoplasm to the periplasm. In Salmonella, there are 30 proteins that are predicted substrates of Tat, and among these are enzymes required for anaerobic respiration and peptidoglycan remodeling. We have demonstrated that some conditions that induce bacterial envelope stress activate expression of a ΔtatABC-lacZ fusion in Salmonella enterica serovar Typhimurium. Particularly, the addition of bile salts to the growth medium causes a 3-fold induction of a ΔtatABC-lacZ reporter fusion. Our data demonstrate that this induction is mediated via the phage shock protein (Psp) stress response system protein PspA. Further, we show that deletion of tatABC increases the induction of tatABC expression in bile salts. Indeed, the data suggest significant interaction between PspA and the Tat system in the regulatory response to bile salts. Although we have not identified the precise mechanism of Psp regulation of tatABC, our work shows that PspA is involved in the activation of tatABC expression by bile salts and adds another layer of complexity to the Salmonella response to envelope stress. IMPORTANCE Salmonella species cause an array of diseases in a variety of hosts. This research is significant in showing induction of the Tat system as a defense against periplasmic stress. Understanding the underlying mechanism of this regulation broadens our understanding of the Salmonella stress response, which is critical to the ability of the organism to cause infection.

Staphylococcus aureus ftnA 3'-Untranslated Region Modulates Ferritin Production Facilitating Growth Under Iron Starvation Conditions

Iron acquisition and modulation of its intracellular concentration are critical for the development of all living organisms. So far, several proteins have been described to be involved in iron homeostasis. Among them, ferritins act as the major iron storage proteins, sequestering internalized iron and modulating its concentration inside bacterial cells. We previously described that the deletion of the 3’-untranslated region (3’UTR) of the ftnA gene, which codes for ferritin in Staphylococcus aureus, increased the ftnA mRNA and ferritin levels. Here, we show that the ferritin levels are affected by RNase III and PNPase, which target the ftnA 3’UTR. Rifampicin mRNA stability experiments revealed that the half-life of the ftnA mRNA is affected by both RNase III and the ftnA 3’UTR. A transcriptional fusion of the ftnA 3’UTR to the gfp reporter gene decreased green fluorescent protein (GFP) expression, indicating that the ftnA 3’UTR could work as an independent module. Additionally, a chromosomal deletion of the ftnA 3’UTR impaired S. aureus growth under conditions of iron starvation. Overall, this work highlights the biological relevance of the ftnA 3’UTR for iron homeostasis in S. aureus.

Ubiquitous mRNA decay fragments in E. coli redefine the functional transcriptome

Bacterial mRNAs have short life cycles, in which transcription is rapidly followed by translation and degradation within seconds to minutes. The resulting diversity of mRNA molecules across different life-cycle stages impacts their functionality but has remained unresolved. Here we quantitatively map the 3’ status of cellular RNAs in Escherichia coli during steady-state growth and report a large fraction of molecules (median>60%) that are fragments of canonical full-length mRNAs. The majority of RNA fragments are decay intermediates, whereas nascent RNAs contribute to a smaller fraction. Despite the prevalence of decay intermediates in total cellular RNA, these intermediates are underrepresented in the pool of ribosome-associated transcripts and can thus distort quantifications and differential expression analyses for the abundance of full-length, functional mRNAs. The large heterogeneity within mRNA molecules in vivo highlights the importance in discerning functional transcripts and provides a lens for studying the dynamic life cycle of mRNAs.

Reassembling a cannon in the DNA defense arsenal: Genetics of StySA, a BREX phage exclusion system in Salmonella lab strains

Understanding mechanisms that shape horizontal exchange in prokaryotes is a key problem in biology. A major limit on DNA entry is imposed by restriction-modification (RM) processes that depend on the pattern of DNA modification at host-specified sites. In classical RM, endonucleolytic DNA cleavage follows detection of unprotected sites on entering DNA. Recent investigation has uncovered BREX (BacteRiophage EXclusion) systems. These RM-like activities employ host protection by DNA modification, but immediate replication arrest occurs without evident of nuclease action on unmodified phage DNA. Here we show that the historical stySA RM locus of Salmonella enterica sv Typhimurium is a variant BREX system. A laboratory strain disabled for both the restriction and methylation activity of StySA nevertheless has wild type sequence in pglX, the modification gene homolog. Instead, flanking genes pglZ and brxC each carry multiple mutations (μ) in their C-terminal domains. We further investigate this system in situ, replacing the mutated pglZμ and brxCμ genes with the WT counterpart. PglZ-WT supports methylation in the presence of either BrxCμ or BrxC-WT but not in the presence of a deletion/insertion allele, ΔbrxC::cat. Restriction requires both BrxC-WT and PglZ-WT, implicating the BrxC C-terminus specifically in restriction activity. These results suggests that while BrxC, PglZ and PglX are principal components of the BREX modification activity, BrxL is required for restriction only. Furthermore, we show that a partial disruption of brxL disrupts transcription globally.

Horizontal gene transfer is a major driver of evolution and adaptation in bacteria. Genes from outside may be beneficial or dangerous to the receiving cell. Benefits include new food sources such as sugars, or new homes by adhesion, or new resistances, as to antibiotics. Dangers are posed by bacteriophages--viruses that take over the cell machinery, multiply, and release progeny to kill sister cells. Host-dependent restriction-modification systems enable defense that distinguishes relatives from strangers: using a modification pattern (M) carried by DNA bases added by the host cell to prevent restriction (R). Sisters and cousin cells will have the same protective pattern on DNA, while DNA of foreign origin will have the wrong M pattern and be restricted (R, rejected). Typically, restriction involves nuclease digestion. Here we address the enigmatic StySA RM system, one of the earliest to be genetically characterized. It is a variant of the newly recognized defense mechanism, BREX. BREX systems also track DNA history via modification pattern, but restrict by a novel, uncharacterized mechanism. Like other BREX family systems, StySA-BREX modification requires multiple components. When StySA-BREX transcription is unbalanced, we find global disruption of gene transcription. The disruption pattern does not suggest SOS-inducing damage to DNA.

Synthetic 3'-UTR valves for optimal metabolic flux control in Escherichia coli

As the design of genetic circuitry for synthetic biology becomes more sophisticated, diverse regulatory bioparts are required. Despite their importance, well-characterized 3′-untranslated region (3′-UTR) bioparts are limited. Thus, transcript 3′-ends require further investigation to understand the underlying regulatory role and applications of the 3′-UTR. Here, we revisited the use of Term-Seq in the Escherichia coli strain K-12 MG1655 to enhance our understanding of 3′-UTR regulatory functions and to provide a diverse collection of tunable 3′-UTR bioparts with a wide termination strength range. Comprehensive analysis of 1,629 transcript 3′-end positions revealed multiple 3′-termini classes generated through transcription termination and RNA processing. The examination of individual Rho-independent terminators revealed a reduction in downstream gene expression over a wide range, which led to the design of novel synthetic metabolic valves that control metabolic fluxes in branched pathways. These synthetic metabolic valves determine the optimal balance of heterologous pathways for maximum target biochemical productivity. The regulatory strategy using 3′-UTR bioparts is advantageous over promoter- or 5′-UTR-based transcriptional control as it modulates gene expression at transcription levels without trans-acting element requirements (e.g. transcription factors). Our results provide a foundational platform for 3′-UTR engineering in synthetic biology applications.

OperonSEQer: A set of machine-learning algorithms with threshold voting for detection of operon pairs using short-read RNA-sequencing data

Operon prediction in prokaryotes is critical not only for understanding the regulation of endogenous gene expression, but also for exogenous targeting of genes using newly developed tools such as CRISPR-based gene modulation. A number of methods have used transcriptomics data to predict operons, based on the premise that contiguous genes in an operon will be expressed at similar levels. While promising results have been observed using these methods, most of them do not address uncertainty caused by technical variability between experiments, which is especially relevant when the amount of data available is small. In addition, many existing methods do not provide the flexibility to determine the stringency with which genes should be evaluated for being in an operon pair. We present OperonSEQer, a set of machine learning algorithms that uses the statistic and p-value from a non-parametric analysis of variance test (Kruskal-Wallis) to determine the likelihood that two adjacent genes are expressed from the same RNA molecule. We implement a voting system to allow users to choose the stringency of operon calls depending on whether your priority is high recall or high specificity. In addition, we provide the code so that users can retrain the algorithm and re-establish hyperparameters based on any data they choose, allowing for this method to be expanded as additional data is generated. We show that our approach detects operon pairs that are missed by current methods by comparing our predictions to publicly available long-read sequencing data. OperonSEQer therefore improves on existing methods in terms of accuracy, flexibility, and adaptability.

Bacteria and archaea, single-cell organisms collectively known as prokaryotes, live in all imaginable environments and comprise the majority of living organisms on this planet. Prokaryotes play a critical role in the homeostasis of multicellular organisms (such as animals and plants) and ecosystems. In addition, bacteria can be pathogenic and cause a variety of diseases in these same hosts and ecosystems. In short, understanding the biology and molecular functions of bacteria and archaea and devising mechanisms to engineer and optimize their properties are critical scientific endeavors with significant implications in healthcare, agriculture, manufacturing, and climate science among others. One major molecular difference between unicellular and multicellular organisms is the way they express genes–multicellular organisms make individual RNA molecules for each gene while, prokaryotes express operons (i.e., a group of genes coding functionally related proteins) in contiguous polycistronic RNA molecules. Understanding which genes exist within operons is critical for elucidating basic biology and for engineering organisms. In this work, we use a combination of statistical and machine learning-based methods to use next-generation sequencing data to predict operon structure across a range of prokaryotes. Our method provides an easily implemented, robust, accurate, and flexible way to determine operon structure in an organism-agnostic manner using readily available data.

Burning the Candle at Both Ends: Have Exoribonucleases Driven Divergence of Regulatory RNA Mechanisms in Bacteria?

Regulatory RNAs have emerged as ubiquitous gene regulators in all bacterial species studied to date. The combination of sequence-specific RNA interactions and malleable RNA structure has allowed regulatory RNA to adopt different mechanisms of gene regulation in a diversity of genetic backgrounds. In the model GammaproteobacteriaEscherichia coli and Salmonella, the regulatory RNA chaperone Hfq appears to play a global role in gene regulation, directly controlling ∼20 to 25% of the entire transcriptome. While the model FirmicutesBacillus subtilis and Staphylococcus aureus encode a Hfq homologue, its role has been significantly depreciated. These bacteria also have marked differences in RNA turnover. E. coli and Salmonella degrade RNA through internal endonucleolytic and 3′→5′ exonucleolytic cleavage that appears to allow transient accumulation of mRNA 3′ UTR cleavage fragments that contain stabilizing 3′ structures. In contrast, B. subtilis and S. aureus are able to exonucleolytically attack internally cleaved RNA from both the 5′ and 3′ ends, efficiently degrading mRNA 3′ UTR fragments. Here, we propose that the lack of 5′→3′ exoribonuclease activity in Gammaproteobacteria has allowed the accumulation of mRNA 3′ UTR ends as the “default” setting. This in turn may have provided a larger pool of unconstrained RNA sequences that has fueled the expansion of Hfq function and small RNA (sRNA) regulation in E. coli and Salmonella. Conversely, the exoribonuclease RNase J may be a significant barrier to the evolution of 3′ UTR sRNAs in B. subtilis and S. aureus that has limited the pool of RNA ligands available to Hfq and other sRNA chaperones, depreciating their function in these model Firmicutes.

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.

The essential role of mRNA degradation in understanding and engineering E. coli metabolism

Metabolic engineering strategies are crucial for the development of bacterial cell factories with improved performance. Until now, optimal metabolic networks have been designed based on systems biology approaches integrating large-scale data on the steady-state concentrations of mRNA, protein and metabolites, sometimes with dynamic data on fluxes, but rarely with any information on mRNA degradation. In this review, we compile growing evidence that mRNA degradation is a key regulatory level in E. coli that metabolic engineering strategies should take into account. We first discuss how mRNA degradation interacts with transcription and translation, two other gene expression processes, to balance transcription regulation and remove poorly translated mRNAs. The many reciprocal interactions between mRNA degradation and metabolism are also highlighted: metabolic activity can be controlled by changes in mRNA degradation and in return, the activity of the mRNA degradation machinery is controlled by metabolic factors. The mathematical models of the crosstalk between mRNA degradation dynamics and other cellular processes are presented and discussed with a view towards novel mRNA degradation-based metabolic engineering strategies. We show finally that mRNA degradation-based strategies have already successfully been applied to improve heterologous protein synthesis. Overall, this review underlines how important mRNA degradation is in regulating E. coli metabolism and identifies mRNA degradation as a key target for innovative metabolic engineering strategies in biotechnology.

Halobacterium salinarum and Haloferax volcanii Comparative Transcriptomics Reveals Conserved Transcriptional Processing Sites

Post-transcriptional processing of messenger RNA is an important regulatory strategy that allows relatively fast responses to changes in environmental conditions. In halophile systems biology, the protein perspective of this problem (i.e., ribonucleases which implement the cleavages) is generally more studied than the RNA perspective (i.e., processing sites). In the present in silico work, we mapped genome-wide transcriptional processing sites (TPS) in two halophilic model organisms, Halobacterium salinarum NRC-1 and Haloferax volcanii DS2. TPS were established by reanalysis of publicly available differential RNA-seq (dRNA-seq) data, searching for non-primary (monophosphorylated RNAs) enrichment. We found 2093 TPS in 43% of H. salinarum genes and 3515 TPS in 49% of H. volcanii chromosomal genes. Of the 244 conserved TPS sites found, the majority were located around start and stop codons of orthologous genes. Specific genes are highlighted when discussing antisense, ribosome and insertion sequence associated TPS. Examples include the cell division gene ftsZ2, whose differential processing signal along growth was detected and correlated with post-transcriptional regulation, and biogenesis of sense overlapping transcripts associated with IS200/IS605. We hereby present the comparative, transcriptomics-based processing site maps with a companion browsing interface.

Global RNA profiles show target selectivity and physiological effects of peptide-delivered antisense antibiotics

Antisense peptide nucleic acids (PNAs) inhibiting mRNAs of essential genes provide a straight-forward way to repurpose our knowledge of bacterial regulatory RNAs for development of programmable species-specific antibiotics. While there is ample proof of PNA efficacy, their target selectivity and impact on bacterial physiology are poorly understood. Moreover, while antibacterial PNAs are typically designed to block mRNA translation, effects on target mRNA levels are not well-investigated. Here, we pioneer the use of global RNA-seq analysis to decipher PNA activity in a transcriptome-wide manner. We find that PNA-based antisense oligomer conjugates robustly decrease mRNA levels of the widely-used target gene, acpP, in Salmonella enterica, with limited off-target effects. Systematic analysis of several different PNA-carrier peptides attached not only shows different bactericidal efficiency, but also activation of stress pathways. In particular, KFF-, RXR- and Tat-PNA conjugates especially induce the PhoP/Q response, whereas the latter two additionally trigger several distinct pathways. We show that constitutive activation of the PhoP/Q response can lead to Tat-PNA resistance, illustrating the utility of RNA-seq for understanding PNA antibacterial activity. In sum, our study establishes an experimental framework for the design and assessment of PNA antimicrobials in the long-term quest to use these for precision editing of microbiota.

Ribosomal RNA degradation induced by the bacterial RNA polymerase inhibitor rifampicin

Rifampicin, a broad-spectrum antibiotic, inhibits bacterial RNA polymerase. Here we show that rifampicin treatment of Escherichia coli results in a 50% decrease in cell size due to a terminal cell division. This decrease is a consequence of inhibition of transcription as evidenced by an isogenic rifampicin-resistant strain. There is also a 50% decrease in total RNA due mostly to a 90% decrease in 23S and 16S rRNA levels. Control experiments showed this decrease is not an artifact of our RNA purification protocol and therefore due to degradation in vivo. Since chromosome replication continues after rifampicin treatment, ribonucleotides from rRNA degradation could be recycled for DNA synthesis. Rifampicin-induced rRNA degradation occurs under different growth conditions and in different strain backgrounds. However, rRNA degradation is never complete thus permitting the re-initiation of growth after removal of rifampicin. The orderly shutdown of growth under conditions where the induction of stress genes is blocked by rifampicin is noteworthy. Inhibition of protein synthesis by chloramphenicol resulted in a partial decrease in 23S and 16S rRNA levels whereas kasugamycin treatment had no effect. Analysis of temperature-sensitive mutant strains implicate RNase E, PNPase and RNase R in rifampicin-induced rRNA degradation. We cannot distinguish between a direct role for RNase E in rRNA degradation versus an indirect role involving a slowdown of mRNA degradation. Since mRNA and rRNA appear to be degraded by the same ribonucleases, competition by rRNA is likely to result in slower mRNA degradation rates in the presence of rifampicin than under normal growth conditions.

Elucidating the Regulatory Elements for Transcription Termination and Posttranscriptional Processing in the Streptomyces clavuligerus Genome

Identification of transcriptional regulatory elements in the GC-rich Streptomyces genome is essential for the production of novel biochemicals from secondary metabolite biosynthetic gene clusters (smBGCs). Despite many efforts to understand the regulation of transcription initiation in smBGCs, information on the regulation of transcription termination and posttranscriptional processing remains scarce. In this study, we identified the transcriptional regulatory elements in β-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27064 by determining a total of 1,427 transcript 3'-end positions (TEPs) using the term-seq method. Termination of transcription was governed by three classes of TEPs, of which each displayed unique sequence features. The data integration with transcription start sites and transcriptome data generated 1,648 transcription units (TUs) and 610 transcription unit clusters (TUCs). TU architecture showed that the transcript abundance in TU isoforms of a TUC was potentially affected by the sequence context of their TEPs, suggesting that the regulatory elements of TEPs could control the transcription level in additional layers. We also identified TU features of a xenobiotic response element (XRE) family regulator and DUF397 domain-containing protein, particularly showing the abundance of bidirectional TEPs. Finally, we found that 189 noncoding TUs contained potential cis- and trans-regulatory elements that played a major role in regulating the 5' and 3' UTR. These findings highlight the role of transcriptional regulatory elements in transcription termination and posttranscriptional processing in Streptomyces sp.IMPORTANCE Streptomyces sp. is a great source of bioactive secondary metabolites, including antibiotics, antifungal agents, antiparasitic agents, immunosuppressant compounds, and other drugs. Secondary metabolites are synthesized via multistep conversions of the precursor molecules from primary metabolism, governed by multicomplex enzymes from secondary metabolite biosynthetic gene clusters. As their production is closely related with the growth phase and dynamic cellular status in response to various intra- and extracellular signals, complex regulatory systems tightly control the gene expressions related to secondary metabolism. In this study, we determined genome-wide transcript 3'-end positions and transcription units in the β-lactam antibiotic producer Streptomyces clavuligerus ATCC 27064 to elucidate the transcriptional regulatory elements in transcription termination and posttranscriptional processing by integration of multiomics data. These unique features, such as transcript 3'-end sequence, potential riboregulators, and potential 3'-untranslated region (UTR) cis-regulatory elements, can be potentially used to design engineering tools that can regulate the transcript abundance of genes for enhancing secondary metabolite production.

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.

A cold shock protein promotes high-temperature microbial growth through binding to diverse RNA species

Endowing mesophilic microorganisms with high-temperature resistance is highly desirable for industrial microbial fermentation. Here, we report a cold-shock protein (CspL) that is an RNA chaperone protein from a lactate producing thermophile strain (Bacillus coagulans 2–6), which is able to recombinantly confer strong high-temperature resistance to other microorganisms. Transgenic cspL expression massively enhanced high-temperature growth of Escherichia coli (a 2.4-fold biomass increase at 45 °C) and eukaryote Saccharomyces cerevisiae (a 2.6-fold biomass increase at 36 °C). Importantly, we also found that CspL promotes growth rates at normal temperatures. Mechanistically, bio-layer interferometry characterized CspL’s nucleotide-binding functions in vitro, while in vivo we used RNA-Seq and RIP-Seq to reveal CspL’s global effects on mRNA accumulation and CspL’s direct RNA binding targets, respectively. Thus, beyond establishing how a cold-shock protein chaperone provides high-temperature resistance, our study introduces a strategy that may facilitate industrial thermal fermentation.

Another layer of complexity in Staphylococcus aureus methionine biosynthesis control: unusual RNase III-driven T-box riboswitch cleavage determines met operon mRNA stability and decay

In Staphylococcus aureus, de novo methionine biosynthesis is regulated by a unique hierarchical pathway involving stringent-response controlled CodY repression in combination with a T-box riboswitch and RNA decay. The T-box riboswitch residing in the 5′ untranslated region (met leader RNA) of the S. aureus metICFE-mdh operon controls downstream gene transcription upon interaction with uncharged methionyl-tRNA. met leader and metICFE-mdh (m)RNAs undergo RNase-mediated degradation in a process whose molecular details are poorly understood. Here we determined the secondary structure of the met leader RNA and found the element to harbor, beyond other conserved T-box riboswitch structural features, a terminator helix which is target for RNase III endoribonucleolytic cleavage. As the terminator is a thermodynamically highly stable structure, it also forms posttranscriptionally in met leader/ metICFE-mdh read-through transcripts. Cleavage by RNase III releases the met leader from metICFE-mdh mRNA and initiates RNase J-mediated degradation of the mRNA from the 5′-end. Of note, metICFE-mdh mRNA stability varies over the length of the transcript with a longer lifespan towards the 3′-end. The obtained data suggest that coordinated RNA decay represents another checkpoint in a complex regulatory network that adjusts costly methionine biosynthesis to current metabolic requirements.

Characterization of Bicistronic Transcription in Budding Yeast

Bicistronic transcripts (operon-like transcripts) have occasionally been reported in eukaryotes, including unicellular yeasts, plants, and humans, despite the fact that they lack trans-splice mechanisms. However, the characteristics of eukaryotic bicistronic transcripts are poorly understood, except for those in nematodes. Here, we describe the genomic, transcriptomic, and ribosome profiling features of bicistronic transcripts in unicellular yeasts. By comparing the expression level of bicistronic transcripts with their monocistronic equivalents, we identify two main categories of bicistronic transcripts: highly and lowly expressed. These two categories exhibit quite different features. First, highly expressed bicistronic transcripts have higher conservation within and between strains and shorter intergenic spacers with higher GC content and less stable secondary structure. Second, genes in highly expressed bicistronic transcripts have lower translation efficiency, with the second gene showing statistically significant lower translation efficiency than the first. Finally, the genes found in these highly expressed bicistronic transcripts tend to be younger, with more recent origins. Together, these results suggest that bicistronic transcripts in yeast are heterogeneous. We further propose that at least some highly expressed bicistronic transcripts appear to play a role in modulating monocistronic translation.IMPORTANCE Operons, where a single mRNA transcript encodes multiple adjacent proteins, are a widespread feature of bacteria and archaea. In contrast, the genes of eukaryotes are generally considered monocistronic. However, a number of studies have revealed the presence of bicistronic transcripts in eukaryotes, including humans. The basic features of these transcripts are largely unknown in eukaryotes, especially in organisms lacking trans-splice mechanisms. Our analyses characterize bicistronic transcripts in one such eukaryotic group, yeasts. We show that highly expressed bicistronic transcripts have unusual features compared to lowly expressed bicistronic transcripts, with several features influencing translational modulation.

Systematic Quantification of Sequence and Structural Determinants Controlling mRNA stability in Bacterial Operons

mRNA degradation is a central process that affects all gene expression levels, and yet, the determinants that control mRNA decay rates remain poorly characterized. Here, we applied a synthetic biology, learn-by-design approach to elucidate the sequence and structural determinants that control mRNA stability in bacterial operons. We designed, constructed, and characterized 82 operons in Escherichia coli, systematically varying RNase binding site characteristics, translation initiation rates, and transcriptional terminator efficiencies in the 5' untranslated region (UTR), intergenic, and 3' UTR regions, followed by measuring their mRNA levels using reverse transcription quantitative polymerase chain reaction (RT-qPCR) assays during exponential growth. We show that introducing long single-stranded RNA into 5' UTRs reduced mRNA levels by up to 9.4-fold and that lowering translation rates reduced mRNA levels by up to 11.8-fold. We also found that RNase binding sites in intergenic regions had much lower effects on mRNA levels. Surprisingly, changing the transcriptional termination efficiency or introducing long single-stranded RNA into 3' UTRs had no effect on upstream mRNA levels. From these measurements, we developed and validated biophysical models of ribosome protection and RNase activity with excellent quantitative agreement. We also formulated design rules to rationally control a mRNA's stability, facilitating the automated design of engineered genetic systems with desired functionalities.

Bacterial 3'UTRs: A Useful Resource in Post-transcriptional Regulation

Bacterial messenger RNAs (mRNAs) are composed of 5′ and 3′ untranslated regions (UTRs) that flank the coding sequences (CDSs). In eukaryotes, 3′UTRs play key roles in post-transcriptional regulatory mechanisms. Shortening or deregulation of these regions is associated with diseases such as cancer and metabolic disorders. Comparatively, little is known about the functions of 3′UTRs in bacteria. Over the past few years, 3′UTRs have emerged as important players in the regulation of relevant bacterial processes such as virulence, iron metabolism, and biofilm formation. This MiniReview is an update for the different 3′UTR-mediated mechanisms that regulate gene expression in bacteria. Some of these include 3′UTRs that interact with the 5′UTR of the same transcript to modulate translation, 3′UTRs that are targeted by specific ribonucleases, RNA-binding proteins and small RNAs (sRNAs), and 3′UTRs that act as reservoirs of trans-acting sRNAs, among others. In addition, recent findings regarding a differential evolution of bacterial 3′UTRs and its impact in the species-specific expression of orthologous genes are also discussed.

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.

CpeT is the phycoerythrobilin lyase for Cys-165 on β-phycoerythrin from Fremyella diplosiphon and the chaperone-like protein CpeZ greatly improves its activity

Bilin lyases are enzymes which ligate linear tetrapyrrole chromophores to specific cysteine residues on light harvesting proteins present in cyanobacteria and red algae. The lyases responsible for chromophorylating the light harvesting protein phycoerythrin (PE) have not been fully characterized. In this study, we explore the role of CpeT, a putative bilin lyase, in the biosynthesis of PE in the cyanobacterium Fremyella diplosiphon. Recombinant protein studies show that CpeT alone can bind phycoerythrobilin (PEB), but CpeZ, a chaperone-like protein, is needed in order to correctly and efficiently attach PEB to the β-subunit of PE. MS analyses of the recombinant β-subunit of PE coexpressed with CpeT and CpeZ show that PEB is attached at Cys-165. Purified phycobilisomes from a cpeT knockout mutant and wild type (WT) samples from F. diplosiphon were analyzed and compared. The cpeT mutant contained much less PE and more phycocyanin than WT cells grown under green light, conditions which should maximize the production of PE. In addition, Northern blot analyses showed that the cpeCDESTR operon mRNAs were upregulated while the cpeBcpeA mRNAs were downregulated in the cpeT mutant strain when compared with WT, suggesting that CpeT may also play a direct or indirect regulatory role in transcription of these operons or their mRNA stability, in addition to its role as a PEB lyase for Cys-165 on β-PE.

Competitive effects in bacterial mRNA decay

In living organisms, the same enzyme catalyses the degradation of thousands of different mRNAs, but the possible influence of competing substrates has been largely ignored so far. We develop a simple mechanistic model of the coupled degradation of all cell mRNAs using the total quasi-steady-state approximation of the Michaelis-Menten framework. Numerical simulations of the model using carefully chosen parameters and analyses of rate sensitivity coefficients show how substrate competition alters mRNA decay. The model predictions reproduce and explain a number of experimental observations on mRNA decay following transcription arrest, such as delays before the onset of degradation, the occurrence of variable degradation profiles with increased non linearities and the negative correlation between mRNA half-life and concentration. The competition acts at different levels, through the initial concentration of cell mRNAs and by modifying the enzyme affinity for its targets. The consequence is a global slow down of mRNA decay due to enzyme titration and the amplification of its apparent affinity. Competition happens to stabilize weakly affine mRNAs and to destabilize the most affine ones. We believe that this mechanistic model is an interesting alternative to the exponential models commonly used for the determination of mRNA half-lives. It allows analysing regulatory mechanisms of mRNA degradation and its predictions are directly comparable to experimental data.

Translation Initiation Control of RNase E-Mediated Decay of Polycistronic gal mRNA

In bacteria, mRNA decay is a major mechanism for regulating gene expression. In Escherichia coli, mRNA decay initiates with endonucleolytic cleavage by RNase E. Translating ribosomes impede RNase E cleavage, thus providing stability to mRNA. In transcripts containing multiple cistrons, the translation of each cistron initiates separately. The effect of internal translation initiations on the decay of polycistronic transcripts remains unknown, which we have investigated here using the four-cistron galETKM transcript. We find that RNase E cleaves a few nucleotides (14-36) upstream of the translation initiation site of each cistron, generating decay intermediates galTKM, galKM, and galM mRNA with fewer but full cistrons. Blocking translation initiation reduced stability, particularly of the mutated cistrons and when they were the 5'-most cistrons. This indicates that, together with translation failure, the location of the cistron is important for its elimination. The instability of the 5'-most cistron did not propagate to the downstream cistrons, possibly due to translation initiation there. Cistron elimination from the 5' end was not always sequential, indicating that RNase E can also directly access a ribosome-free internal cistron. The finding in gal operon of mRNA decay by cistron elimination appears common in E. coli and Salmonella.

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.

The Leader Peptide peTrpL Forms Antibiotic-Containing Ribonucleoprotein Complexes for Posttranscriptional Regulation of Multiresistance Genes

Bacterial ribosome-dependent attenuators are widespread posttranscriptional regulators. They harbor small upstream open reading frames (uORFs) encoding leader peptides, for which no functions in trans are known yet. In the plant symbiont Sinorhizobium meliloti, the tryptophan biosynthesis gene trpE(G) is preceded by the uORF trpL and is regulated by transcription attenuation according to tryptophan availability. However, trpLE(G) transcription is initiated independently of the tryptophan level in S. meliloti, thereby ensuring a largely tryptophan-independent production of the leader peptide peTrpL. Here, we provide evidence for a tryptophan-independent role of peTrpL in trans We found that peTrpL increases the resistance toward tetracycline, erythromycin, chloramphenicol, and the flavonoid genistein, which are substrates of the major multidrug efflux pump SmeAB. Coimmunoprecipitation with a FLAG-peTrpL suggested smeR mRNA, which encodes the transcription repressor of smeABR, as a peptide target. Indeed, upon antibiotic exposure, smeR mRNA was destabilized and smeA stabilized in a peTrpL-dependent manner, showing that peTrpL acts in the differential regulation of smeABR Furthermore, smeR mRNA was coimmunoprecipitated with peTrpL in antibiotic-dependent ribonucleoprotein (ARNP) complexes, which, in addition, contained an antibiotic-induced antisense RNA complementary to smeRIn vitro ARNP reconstitution revealed that the above-mentioned antibiotics and genistein directly support complex formation. A specific region of the antisense RNA was identified as a seed region for ARNP assembly in vitro Altogether, our data show that peTrpL is involved in a mechanism for direct utilization of antimicrobial compounds in posttranscriptional regulation of multiresistance genes. Importantly, this role of peTrpL in resistance is conserved in other AlphaproteobacteriaIMPORTANCE Leader peptides encoded by transcription attenuators are widespread small proteins that are considered nonfunctional in trans We found that the leader peptide peTrpL of the soil-dwelling plant symbiont Sinorhizobium meliloti is required for differential, posttranscriptional regulation of a multidrug resistance operon upon antibiotic exposure. Multiresistance achieved by efflux of different antimicrobial compounds ensures survival and competitiveness in nature and is important from both evolutionary and medical points of view. We show that the leader peptide forms antibiotic- and flavonoid-dependent ribonucleoprotein complexes (ARNPs) for destabilization of smeR mRNA encoding the transcription repressor of the major multidrug resistance operon. The seed region for ARNP assembly was localized in an antisense RNA, whose transcription is induced by antimicrobial compounds. The discovery of ARNP complexes as new players in multiresistance regulation opens new perspectives in understanding bacterial physiology and evolution and potentially provides new targets for antibacterial control.

CRISPR Tools To Control Gene Expression in Bacteria

CRISPR-Cas systems have been engineered as powerful tools to control gene expression in bacteria. The most common strategy relies on the use of Cas effectors modified to bind target DNA without introducing DNA breaks. These effectors can either block the RNA polymerase or recruit it through activation domains. Here, we discuss the mechanistic details of how Cas effectors can modulate gene expression by blocking transcription initiation or acting as transcription roadblocks. CRISPR-Cas tools can be further engineered to obtain fine-tuned control of gene expression or target multiple genes simultaneously. Several caveats in using these tools have also been revealed, including off-target effects and toxicity, making it important to understand the design rules of engineered CRISPR-Cas effectors in bacteria. Alternatively, some types of CRISPR-Cas systems target RNA and could be used to block gene expression at the posttranscriptional level. Finally, we review applications of these tools in high-throughput screens and the progress and challenges in introducing CRISPR knockdown to other species, including nonmodel bacteria with industrial or clinical relevance. A deep understanding of how CRISPR-Cas systems can be harnessed to control gene expression in bacteria and build powerful tools will certainly open novel research directions.

Dynamic insights on transcription initiation and RNA processing during bacterial adaptation

Transcription initiation and RNA processing govern gene expression and enable bacterial adaptation by reshaping the RNA landscape. The aim of this study was to simultaneously observe these two fundamental processes in a transcriptome responding to an environmental signal. A controlled σ^E system in E. coli was coupled to our previously described tagRNA-seq method to yield process kinetics information. Changes in transcription initiation frequencies (TIF) and RNA processing frequencies (PF) were followed using 5' RNA tags. Changes in TIF showed a binary increased/decreased pattern that alternated between transcriptionally activated and repressed promoters, providing the bacterial population with transcriptional oscillation. PF variation fell into three categories of cleavage activity: (i) constant and independent of RNA levels, (ii) increased once RNA has accumulated, and (iii) positively correlated to changes in TIF. This work provides a comprehensive and dynamic view of major events leading to transcriptomic reshaping during bacterial adaptation. It unveils an interplay between transcription initiation and the activity of specific RNA cleavage sites. This study utilized a well-known genetic system to analyze fundamental processes and can serve as a blueprint for comprehensive studies that exploit the RNA metabolism to decipher and understand bacterial gene expression control.

RNA helicase-regulated processing of the Synechocystis rimO-crhR operon results in differential cistron expression and accumulation of two sRNAs

The arrangement of functionally-related genes in operons is a fundamental element of how genetic information is organized in prokaryotes. This organization ensures coordinated gene expression by co-transcription. Often, however, alternative genetic responses to specific stress conditions demand the discoordination of operon expression. During cold temperature stress, accumulation of the gene encoding the sole Asp-Glu-Ala-Asp (DEAD)-box RNA helicase in Synechocystis sp. PCC 6803, crhR (slr0083), increases 15-fold. Here, we show that crhR is expressed from a dicistronic operon with the methylthiotransferase rimO/miaB (slr0082) gene, followed by rapid processing of the operon transcript into two monocistronic mRNAs. This cleavage event is required for and results in destabilization of the rimO transcript. Results from secondary structure modeling and analysis of RNase E cleavage of the rimO-crhR transcript in vitro suggested that CrhR plays a role in enhancing the rate of the processing in an auto-regulatory manner. Moreover, two putative small RNAs are generated from additional processing, degradation, or both of the rimO transcript. These results suggest a role for the bacterial RNA helicase CrhR in RNase E-dependent mRNA processing in Synechocystis and expand the known range of organisms possessing small RNAs derived from processing of mRNA transcripts.

Switching fatty acid metabolism by an RNA-controlled feed forward loop

Hfq (host factor for phage Q beta) is key for posttranscriptional gene regulation in many bacteria. Hfq's function is to stabilize sRNAs and to facilitate base-pairing with trans-encoded target mRNAs. Loss of Hfq typically results in pleiotropic phenotypes, and, in the major human pathogen Vibrio cholerae, Hfq inactivation has been linked to reduced virulence, failure to produce biofilms, and impaired intercellular communication. However, the RNA ligands of Hfq in V. cholerae are currently unknown. Here, we used RIP-seq (RNA immunoprecipitation followed by high-throughput sequencing) analysis to identify Hfq-bound RNAs in V. cholerae Our work revealed 603 coding and 85 noncoding transcripts associated with Hfq, including 44 sRNAs originating from the 3' end of mRNAs. Detailed investigation of one of these latter transcripts, named FarS (fatty acid regulated sRNA), showed that this sRNA is produced by RNase E-mediated maturation of the fabB 3'UTR, and, together with Hfq, inhibits the expression of two paralogous fadE mRNAs. The fabB and fadE genes are antagonistically regulated by the major fatty acid transcription factor, FadR, and we show that, together, FadR, FarS, and FadE constitute a mixed feed-forward loop regulating the transition between fatty acid biosynthesis and degradation in V. cholerae Our results provide the molecular basis for studies on Hfq in V. cholerae and highlight the importance of a previously unrecognized sRNA for fatty acid metabolism in this major human pathogen.

Phosphate Transporter PstSCAB of Campylobacter jejuni Is a Critical Determinant of Lactate-Dependent Growth and Colonization in Chickens

Campylobacter jejuni causes acute gastroenteritis worldwide and is transmitted primarily through poultry, in which it is often a commensal member of the intestinal microbiota. Previous transcriptome sequencing (RNA-Seq) experiment showed that transcripts from an operon encoding a high-affinity phosphate transporter (PstSCAB) of C. jejuni were among the most abundant when the bacterium was grown in chickens. Elevated levels of the pstSCAB mRNA were also identified in an RNA-Seq experiment from human infection studies. In this study, we explore the role of PstSCAB in the biology and colonization potential of C. jejuni Our results demonstrate that cells lacking PstSCAB survive poorly in stationary phase, in nutrient-limiting media, and under osmotic conditions reflective of those in the chicken. Polyphosphate levels in the mutant cells were elevated at stationary phase, consistent with alterations in expression of polyphosphate metabolism genes. The mutant strain was highly attenuated for colonization of newly hatched chicks, with levels of bacteria at several orders of magnitude below wild-type levels. Mutant and wild type grew similarly in complex media, but the pstS::kan mutant exhibited a significant growth defect in minimal medium supplemented with l-lactate, postulated as a carbon source in vivo Poor growth in lactate correlated with diminished expression of acetogenesis pathway genes previously demonstrated as important for colonizing chickens. The phosphate transport system is thus essential for diverse aspects of C. jejuni physiology and in vivo fitness and survival.IMPORTANCE Campylobacter jejuni causes millions of human gastrointestinal infections annually, with poultry a major source of infection. Due to the emergence of multidrug resistance in C. jejuni, there is need to identify alternative ways to control this pathogen. Genes encoding the high-affinity phosphate transporter PstSCAB are highly expressed by C. jejuni in chickens and humans. In this study, we address the role of PstSCAB on chicken colonization and other C. jejuni phenotypes. PstSCAB is required for colonization in chicken, metabolism and survival under different stress responses, and during growth on lactate, a potential growth substrate in chickens. Our study highlights that PstSCAB may be an effective target to develop mechanisms for controlling bacterial burden in both chicken and human.

Trends in Symbiont-Induced Host Cellular Differentiation

Bacteria participate in a wide diversity of symbiotic associations with eukaryotic hosts that require precise interactions for bacterial recognition and persistence. Most commonly, host-associated bacteria interfere with host gene expression to modulate the immune response to the infection. However, many of these bacteria also interfere with host cellular differentiation pathways to create a hospitable niche, resulting in the formation of novel cell types, tissues, and organs. In both of these situations, bacterial symbionts must interact with eukaryotic regulatory pathways. Here, we detail what is known about how bacterial symbionts, from pathogens to mutualists, control host cellular differentiation across the central dogma, from epigenetic chromatin modifications, to transcription and mRNA processing, to translation and protein modifications. We identify four main trends from this survey. First, mechanisms for controlling host gene expression appear to evolve from symbionts co-opting cross-talk between host signaling pathways. Second, symbiont regulatory capacity is constrained by the processes that drive reductive genome evolution in host-associated bacteria. Third, the regulatory mechanisms symbionts exhibit correlate with the cost/benefit nature of the association. And, fourth, symbiont mechanisms for interacting with host genetic regulatory elements are not bound by native bacterial capabilities. Using this knowledge, we explore how the ubiquitous intracellular Wolbachia symbiont of arthropods and nematodes may modulate host cellular differentiation to manipulate host reproduction. Our survey of the literature on how infection alters gene expression in Wolbachia and its hosts revealed that, despite their intermediate-sized genomes, different strains appear capable of a wide diversity of regulatory manipulations. Given this and Wolbachia's diversity of phenotypes and eukaryotic-like proteins, we expect that many symbiont-induced host differentiation mechanisms will be discovered in this system.

Transcript degradation and codon usage regulate gene expression in a lytic phage

Many viral genomes are small, containing only single- or double-digit numbers of genes and relatively few regulatory elements. Yet viruses successfully execute complex regulatory programs as they take over their host cells. Here, we propose that some viruses regulate gene expression via a carefully balanced interplay between transcription, translation, and transcript degradation. As our model system, we employ bacteriophage T7, whose genome of approximately sixty genes is well annotated and for which there is a long history of computational models of gene regulation. We expand upon prior modeling work by implementing a stochastic gene expression simulator that tracks individual transcripts, polymerases, ribosomes, and ribonucleases participating in the transcription, translation, and transcript-degradation processes occurring during a T7 infection. By combining this detailed mechanistic modeling of a phage infection with high-throughput gene expression measurements of several strains of bacteriophage T7, evolved and engineered, we can show that both the dynamic interplay between transcription and transcript degradation, and between these two processes and translation, appear to be critical components of T7 gene regulation. Our results point to targeted degradation as a generic gene regulation strategy that may have evolved in many other viruses. Further, our results suggest that detailed mechanistic modeling may uncover the biological mechanisms at work in both evolved and engineered virus variants.

Control of Type III Secretion System Effector/Chaperone Ratio Fosters Pathogen Adaptation to Host-Adherent Lifestyle

Host colonization by extracellular pathogens often entails the transition from a planktonic lifestyle to a host-attached state. Enteropathogenic E. coli (EPEC), a Gram-negative pathogen, attaches to the intestinal epithelium of the host and employs a type III secretion system (T3SS) to inject effector proteins into the cytoplasm of infected cells. The most abundant effector protein injected is Tir, whose translocation is dependent on the Tir-bound chaperon CesT. Upon Tir injection, the liberated CesT binds to and inhibits the posttranscriptional regulator CsrA, resulting in reprogramming of gene expression in the host-attached bacteria. Thus, adaptation to the host-attached state involves dynamic remodeling of EPEC gene expression, which is mediated by the relative levels of Tir and CesT. Fluctuating from the optimal Tir/CesT ratio results in a decrease in EPEC virulence. Here we elucidate a posttranscriptional circuit that prevents sharp variations from this ratio, thus improving host colonization.

The transition from a planktonic lifestyle to a host-attached state is often critical for bacterial virulence. Upon attachment to host cells, enteropathogenic Escherichia coli (EPEC) employs a type III secretion system (T3SS) to inject into the host cells ∼20 effector proteins, including Tir. CesT, which is encoded from the same operon downstream of tir, is a Tir-bound chaperone that facilitates Tir translocation. Upon Tir translocation, the liberated CesT remains in the bacterial cytoplasm and antagonizes the posttranscriptional regulator CsrA, thus eliciting global regulation in the infecting pathogen. Importantly, tight control of the Tir/CesT ratio is vital, since an uncontrolled surge in free CesT levels may repress CsrA in an untimely manner, thus abrogating colonization. We investigated how fluctuations in Tir translation affect the regulation of this ratio. By creating mutations that cause the premature termination of Tir translation, we revealed that the untranslated tir mRNA becomes highly unstable, resulting in a rapid drop in cesT mRNA levels and, thus, CesT levels. This mechanism couples Tir and CesT levels to ensure a stable Tir/CesT ratio. Our results expose an additional level of regulation that enhances the efficacy of the initial interaction of EPEC with the host cell, providing a better understanding of the bacterial switch from the planktonic to the cell-adherent lifestyle.

Full-length RNA profiling reveals pervasive bidirectional transcription terminators in bacteria

The ability to determine full-length nucleotide composition of individual RNA molecules is essential for understanding the architecture and function of a transcriptome. However, experimental approaches capable of capturing the sequences of both 5' and 3' termini of the same transcript remain scarce. In the present study, simultaneous 5' and 3' end sequencing (SEnd-seq)-a high-throughput and unbiased method that simultaneously maps transcription start and termination sites with single-nucleotide resolution-is presented. Using this method, a comprehensive view of the Escherichia coli transcriptome was obtained, which displays an unexpected level of complexity. SEnd-seq notably expands the catalogue of transcription start sites and termination sites, defines unique transcription units and detects prevalent antisense RNA. Strikingly, the results of the present study unveil widespread overlapping bidirectional terminators located between opposing gene pairs. Furthermore, it has been shown that convergent transcription is a major contributor to highly efficient bidirectional termination both in vitro and in vivo. This finding highlights an underappreciated role of RNA polymerase conflicts in shaping transcript boundaries and suggests an evolutionary strategy for modulating transcriptional output by arranging gene orientation.

Substrate Recognition and Autoinhibition in the Central Ribonuclease RNase E

The endoribonuclease RNase E is a principal factor in RNA turnover and processing that helps to exercise fine control of gene expression in bacteria. While its catalytic activity can be strongly influenced by the chemical identity of the 5′ end of RNA substrates, the enzyme can also cleave numerous substrates irrespective of the chemistry of their 5′ ends through a mechanism that has remained largely unexplained. We report structural and functional data illuminating details of both operational modes. Our crystal structure of RNase E in complex with the sRNA RprA reveals a duplex recognition site that saddles an inter-protomer surface to help present substrates for cleavage. Our data also reveal an autoinhibitory pocket that modulates the overall activity of the ribonuclease. Taking these findings together, we propose how RNase E uses versatile modes of RNA recognition to achieve optimal activity and specificity.

•

RNase E recognizes RNA secondary structure

•

Signature on the substrate 5′ end recognizes and activates RNase E

•

RNase E intrinsic activity is repressed by a conserved autoinhibition motif

Bacterial Noncoding RNAs Excised from within Protein-Coding Transcripts

Bacteria and archaea utilize regulatory small noncoding RNAs (ncRNAs) to control the expression of specific genetic programs. These ncRNAs are almost exclusively encoded within intergenic regions and are independently transcribed. Here, we report on a large set ncRNAs that are “carved out” from within the protein-coding regions of Escherichia coli mRNAs by cellular RNases. These protected mRNA fragments fold into energetically stable RNA structures, reminiscent of those of intergenic regulatory ncRNAs. In addition, a subset of these ncRNAs coprecipitate with the major ncRNA chaperones Hfq and ProQ and display evolutionarily conserved sequences and conserved expression patterns between different bacterial species. Our data suggest that protein-coding genes can potentially act as a reservoir of regulatory ncRNAs.

Prokaryotic genomes encode a plethora of small noncoding RNAs (ncRNAs) that fine-tune the expression of specific genes. The vast majority of known bacterial ncRNAs are encoded from within intergenic regions, where their expression is controlled by promoter and terminator elements, similarly to protein-coding genes. In addition, recent studies have shown that functional ncRNAs can also be derived from gene 3′ untranslated regions (3′UTRs) via an alternative biogenesis pathway, in which the ncRNA segment is separated from the mRNA via RNase cleavage. Here, we report the detection of a large set of decay-generated noncoding RNAs (decRNAs), many of which are completely embedded within protein-coding mRNA regions rather than in the UTRs. We show that these decRNAs are “carved out” of the mRNA through the action of RNase E and that they are predicted to fold into highly stable RNA structures, similar to those of known ncRNAs. A subset of these decRNAs is predicted to interact with Hfq or ProQ or both, which act as ncRNA chaperones, and some decRNAs display evolutionarily conserved sequences and conserved expression patterns between different species. These results suggest that mRNA protein-coding regions may harbor a large set of potentially functional small RNAs.

Neisseria gonorrhoeae Exposed to Sublethal Levels of Hydrogen Peroxide Mounts a Complex Transcriptional Response

Neisseria gonorrhoeae mounts a substantial transcriptional program in response to hydrogen peroxide (HP), a prominent reactive oxygen species (ROS) encountered during infection. We tested which strain FA1090 genes show differential transcript abundance in response to sublethal amounts of HP to differentiate HP-responsive signaling from widespread cellular death and dysregulation. RNA sequencing (RNA-Seq) revealed that 150 genes were significantly upregulated and 143 genes downregulated following HP exposure. We annotated HP-responsive operons and all transcriptional start sites (TSSs) and identified which TSSs responded to HP treatment. We compared the HP responses and other previously reported genes and found only partial overlapping of other regulatory networks, indicating that the response to HP involves multiple biological functions. Using a representative subset of responsive genes, we validated the RNA-Seq results and found that the HP transcriptome was similar to that of sublethal organic peroxide. None of the genes in the representative subset, however, responded to sublethal levels of HOCl or O2 -. These results support the idea that N. gonorrhoeae may use variations in HP levels as a signal for different stages of infection. IMPORTANCE The strict human pathogen Neisseria gonorrhoeae is the only causative agent of the sexually transmitted disease gonorrhea. This bacterium encounters hydrogen peroxide produced from host cells during infection, but the organism survives in the presence of this antimicrobial agent. This work shows that the bacterium responds to hydrogen peroxide by regulating the expression of many genes involved in multiple processes.