Comparative genomics boosts target prediction for bacterial small RNAs

Exploring the targetome of IsrR, an iron-regulated sRNA controlling the synthesis of iron-containing proteins in Staphylococcus aureus

Staphylococcus aureus is a common colonizer of the skin and nares of healthy individuals, but also a major cause of severe human infections. During interaction with the host, pathogenic bacteria must adapt to a variety of adverse conditions including nutrient deprivation. In particular, they encounter severe iron limitation in the mammalian host through iron sequestration by haptoglobin and iron-binding proteins, a phenomenon called "nutritional immunity." In most bacteria, including S. aureus, the ferric uptake regulator (Fur) is the key regulator of iron homeostasis, which primarily acts as a transcriptional repressor of genes encoding iron acquisition systems. Moreover, Fur can control the expression of trans-acting small regulatory RNAs that play an important role in the cellular iron-sparing response involving major changes in cellular metabolism under iron-limiting conditions. In S. aureus, the sRNA IsrR is controlled by Fur, and most of its predicted targets are iron-containing proteins and other proteins related to iron metabolism and iron-dependent pathways. To characterize the IsrR targetome on a genome-wide scale, we combined proteomics-based identification of potential IsrR targets using S. aureus strains either lacking or constitutively expressing IsrR with an in silico target prediction approach, thereby suggesting 21 IsrR targets, of which 19 were negatively affected by IsrR based on the observed protein patterns. These included several Fe-S cluster- and heme-containing proteins, such as TCA cycle enzymes and catalase encoded by katA. IsrR affects multiple metabolic pathways connected to the TCA cycle as well as the oxidative stress response of S. aureus and links the iron limitation response to metabolic remodeling. In contrast to the majority of target mRNAs, the IsrR-katA mRNA interaction is predicted upstream of the ribosome binding site, and further experiments including mRNA half-life measurements demonstrated that IsrR, in addition to inhibiting translation initiation, can downregulate target protein levels by affecting mRNA stability.

Small RNA-modulated anaerobic respiration allows bacteria to survive under antibiotic stress conditions

Despite extensive knowledge of antibiotic-targeted bacterial cell death, deeper understanding of antibiotic tolerance mechanisms is necessary to combat multi-drug resistance in the global healthcare settings. Regulatory RNAs in bacteria control important cellular processes such as cell division, cellular respiration, metabolism, and virulence. Here, we investigated how exposing Escherichia coli to the moderately effective first-generation antibiotic cephalothin alters transcriptional and post-transcriptional dynamics. Bacteria switched from active aerobic respiration to anaerobic adaptation via an FnrS and Tp2 small RNA-mediated post-transcriptional regulatory circuit. From the early hours of antibiotic exposure, FnrS was involved in regulating reactive oxygen species levels, and delayed oxygen consumption in bacteria. We demonstrated that bacteria strive to maintain cellular homeostasis via sRNA-mediated sudden respiratory changes upon sublethal antibiotic exposure.

sRNA molecules participate in hyperosmotic stress response regulation in Sphingomonas melonis TY

Drought and salinity are ubiquitous environmental factors that pose hyperosmotic threats to microorganisms and impair their efficiency in performing environmental functions. However, bacteria have developed various responses and regulatory systems to cope with these abiotic challenges. Posttranscriptional regulation plays vital roles in regulating gene expression and cellular homeostasis, as hyperosmotic stress conditions can lead to the induction of specific small RNA molecules (sRNAs) that participate in stress response regulation. Here, we report a candidate functional sRNA landscape of Sphingomonas melonis TY under hyperosmotic stress, and 18 sRNAs were found with a clear response to hyperosmotic stress. These findings will help in the comprehensive analysis of sRNA regulation in Sphingomonas species. Weighted correlation network analysis revealed a 263 nucleotide sRNA, SNC251, which was transcribed from its own promoter and showed the most significant correlation with hyperosmotic response factors. Deletion of snc251 affected biofilm formation and multiple cellular processes, including ribosome-related pathways, aromatic compound degradation, and the nicotine degradation capacity of S. melonis TY, while overexpression of SNC251 facilitated biofilm formation by TY under hyperosmotic stress. Two genes involved in the TonB system were further verified to be activated by SNC251, which also indicated that SNC251 is a trans-acting sRNA. Briefly, this research reports a landscape of sRNAs participating in the hyperosmotic stress response in S. melonis and reveals a novel sRNA, SNC251, which contributes to the S. melonis TY biofilm formation and thus enhances its hyperosmotic stress response ability.IMPORTANCESphingomonas species play a vital role in plant defense and pollutant degradation and survive extensively under drought or salinity. Previous studies have focused on the transcriptional and translational responses of Sphingomonas under hyperosmotic stress, but the posttranscriptional regulation of small RNA molecules (sRNAs) is also crucial for quickly modulating cellular processes to adapt dynamically to osmotic environments. In addition, the current knowledge of sRNAs in Sphingomonas is extremely scarce. This research revealed a novel sRNA landscape of Sphingomonas melonis and will greatly enhance our understanding of sRNAs' acting mechanisms in the hyperosmotic stress response.

Balanced cell division is secured by two different regulatory sites in OxyS RNA

The hydrogen peroxide-induced small RNA OxyS has been proposed to originate from the 3' UTR of a peroxide mRNA. Unexpectedly, phylogenetic OxyS targetome predictions indicate that most OxyS targets belong to the category of "cell cycle," including cell division and cell elongation. Previously, we reported that Escherichia coli OxyS inhibits cell division by repressing expression of the essential transcription termination factor nusG, thereby leading to the expression of the KilR protein, which interferes with the function of the major cell division protein, FtsZ. By interfering with cell division, OxyS brings about cell-cycle arrest, thus allowing DNA damage repair. Cell division and cell elongation are opposing functions to the extent that inhibition of cell division requires a parallel inhibition of cell elongation for the cells to survive. In this study, we report that in addition to cell division, OxyS inhibits mepS, which encodes an essential peptidoglycan endopeptidase that is responsible for cell elongation. Our study indicates that cell-cycle arrest and balancing between cell division and cell elongation are important and conserved functions of the oxidative stress-induced sRNA OxyS.

How to do RNA-RNA Interaction Prediction? A Use-Case Driven Handbook Using IntaRNA

Computational prediction of RNA-RNA interactions (RRI) is a central methodology for the specific investigation of inter-molecular RNA interactions and regulatory effects of non-coding RNAs like eukaryotic microRNAs or prokaryotic small RNAs. Available methods can be classified according to their underlying prediction strategies, each implicating specific capabilities and restrictions often not transparent to the non-expert user. Within this work, we review seven classes of RRI prediction strategies and discuss the advantages and limitations of respective tools, since such knowledge is essential for selecting the right tool in the first place.Among the RRI prediction strategies, accessibility-based approaches have been shown to provide the most reliable predictions. Here, we describe how IntaRNA, as one of the state-of-the-art accessibility-based tools, can be applied in various use cases for the task of computational RRI prediction. Detailed hands-on examples for individual RRI predictions as well as large-scale target prediction scenarios are provided. We illustrate the flexibility and capabilities of IntaRNA through the examples. Each example is designed using real-life data from the literature and is accompanied by instructions on interpreting the respective results from IntaRNA output. Our use-case driven instructions enable non-expert users to comprehensively understand and utilize IntaRNA's features for effective RRI predictions.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

The Mar, Sox, and Rob Systems

Environments inhabited by Enterobacteriaceae are diverse and often stressful. This is particularly true for Escherichia coli and Salmonella during host association in the gastrointestinal systems of animals. There, E. coli and Salmonella must survive exposure to various antimicrobial compounds produced or ingested by their host. A myriad of changes to cellular physiology and metabolism are required to achieve this feat. A central regulatory network responsible for sensing and responding to intracellular chemical stressors like antibiotics are the Mar, Sox, and Rob systems found throughout the Enterobacteriaceae. Each of these distinct regulatory networks controls expression of an overlapping set of downstream genes whose collective effects result in increased resistance to a wide array of antimicrobial compounds. This collection of genes is known as the mar-sox-rob regulon. This review will provide an overview of the mar-sox-rob regulon and molecular architecture of the Mar, Sox, and Rob systems.

Recent insights into the world of dual-function bacterial sRNAs

Dual-function sRNAs refer to a small subgroup of small regulatory RNAs that merges base-pairing properties of antisense RNAs with peptide-encoding properties of mRNA. Both functions can be part of either same or in another metabolic pathway. Here, we want to update the knowledge of to the already known dual-function sRNAs and review the six new sRNAs found since 2017 regarding their structure, functional mechanisms, evolutionary conservation, and role in the regulation of distinct biological/physiological processes. The increasing identification of dual-function sRNAs through bioinformatics approaches, RNomics and RNA-sequencing and the associated increase in regulatory understanding will likely continue to increase at the same rate in the future. This may improve our understanding of the physiology, virulence and resistance of bacteria, as well as enable their use in technical applications. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

TargetRNA3: predicting prokaryotic RNA regulatory targets with machine learning

Small regulatory RNAs pervade prokaryotes, with the best-studied family of these non-coding genes corresponding to trans-acting regulators that bind via base pairing to their message targets. Given the increasing frequency with which these genes are being identified, it is important that methods for illuminating their regulatory targets keep pace. Using a machine learning approach, we investigate thousands of interactions between small RNAs and their targets, and we interrogate more than a hundred features indicative of these interactions. We present a new method, TargetRNA3, for predicting targets of small RNA regulators and show that it outperforms existing approaches. TargetRNA3 is available at https://cs.wellesley.edu/~btjaden/TargetRNA3 .

The heme-responsive PrrH sRNA regulates Pseudomonas aeruginosa pyochelin gene expression

Pseudomonas aeruginosa is an opportunistic pathogen that requires iron for growth and virulence, yet this nutrient is sequestered by the innate immune system during infection. When iron is limiting, P. aeruginosa expresses the PrrF1 and PrrF2 small RNAs (sRNAs), which post-transcriptionally repress expression of nonessential iron-containing proteins, thus sparing this nutrient for more critical processes. The genes for the PrrF1 and PrrF2 sRNAs are arranged in tandem on the chromosome, allowing for the transcription of a longer heme-responsive sRNA, termed PrrH. While the functions of PrrF1 and PrrF2 have been extensively studied, the role of PrrH in P. aeruginosa physiology and virulence is not well understood. In this study, we performed transcriptomic and proteomic studies to identify the PrrH regulon. In shaking cultures, the pyochelin synthesis proteins were increased in two distinct prrH mutants compared to the wild type, while the mRNAs for these proteins were not affected by the prrH mutation. We identified complementarity between the PrrH sRNA and the sequence upstream of the pchE mRNA, suggesting the potential for PrrH to directly regulate the expression of genes for pyochelin synthesis. We further showed that pchE mRNA levels were increased in the prrH mutants when grown in static but not shaking conditions. Moreover, we discovered that controlling for the presence of light was critical for examining the impact of PrrH on pchE expression. As such, our study reports on the first likely target of the PrrH sRNA and highlights key environmental variables that will allow for future characterization of PrrH function. IMPORTANCE In the human host, iron is predominantly in the form of heme, which Pseudomonas aeruginosa can acquire as an iron source during infection. We previously showed that the iron-responsive PrrF small RNAs (sRNAs) are critical for mediating iron homeostasis during P. aeruginosa infection; however, the function of the heme-responsive PrrH sRNA remains unclear. In this study, we identified genes for pyochelin siderophore biosynthesis, which mediates uptake of inorganic iron, as a novel target of PrrH regulation. This study therefore highlights a novel relationship between heme availability and siderophore biosynthesis in P. aeruginosa.

LinearCoFold and LinearCoPartition: linear-time algorithms for secondary structure prediction of interacting RNA molecules

Many RNAs function through RNA-RNA interactions. Fast and reliable RNA structure prediction with consideration of RNA-RNA interaction is useful, however, existing tools are either too simplistic or too slow. To address this issue, we present LinearCoFold, which approximates the complete minimum free energy structure of two strands in linear time, and LinearCoPartition, which approximates the cofolding partition function and base pairing probabilities in linear time. LinearCoFold and LinearCoPartition are orders of magnitude faster than RNAcofold. For example, on a sequence pair with combined length of 26,190 nt, LinearCoFold is 86.8× faster than RNAcofold MFE mode, and LinearCoPartition is 642.3× faster than RNAcofold partition function mode. Surprisingly, LinearCoFold and LinearCoPartition's predictions have higher PPV and sensitivity of intermolecular base pairs. Furthermore, we apply LinearCoFold to predict the RNA-RNA interaction between SARS-CoV-2 genomic RNA (gRNA) and human U4 small nuclear RNA (snRNA), which has been experimentally studied, and observe that LinearCoFold's prediction correlates better with the wet lab results than RNAcofold's.

Small RNAs, Large Networks: Posttranscriptional Regulons in Gram-Negative Bacteria

Small regulatory RNA (sRNAs) are key mediators of posttranscriptional gene control in bacteria. Assisted by RNA-binding proteins, a single sRNA often modulates the expression of dozens of genes, and thus sRNAs frequently adopt central roles in regulatory networks. Posttranscriptional regulation by sRNAs comes with several unique features that cannot be achieved by transcriptional regulators. However, for optimal network performance, transcriptional and posttranscriptional control mechanisms typically go hand-in-hand. This view is reflected by the ever-growing class of mixed network motifs involving sRNAs and transcription factors, which are ubiquitous in biology and whose regulatory properties we are beginning to understand. In addition, sRNA activity can be antagonized by base-pairing with sponge RNAs, adding yet another layer of complexity to these networks. In this article, we summarize the regulatory concepts underlying sRNA-mediated gene control in bacteria and discuss how sRNAs shape the output of a network, focusing on several key examples.

The Small RNA NcS25 Regulates Biological Amine-Transporting Outer Membrane Porin BCAL3473 in Burkholderia cenocepacia

Regulation of porin expression in bacteria is complex and often involves small-RNA regulators. Several small-RNA regulators have been described for Burkholderia cenocepacia, and this study aimed to characterize the biological role of the conserved small RNA NcS25 and its cognate target, outer membrane protein BCAL3473. The B. cenocepacia genome carries a large number of genes encoding porins with yet-uncharacterized functions. Expression of the porin BCAL3473 is strongly repressed by NcS25 and activated by other factors, such as a LysR-type regulator and nitrogen-depleted growth conditions. The porin is involved in transport of arginine, tyrosine, tyramine, and putrescine across the outer membrane. Porin BCAL3473, with NcS25 as a major regulator, plays an important role in the nitrogen metabolism of B. cenocepacia. IMPORTANCE Burkholderia cenocepacia is a Gram-negative bacterium which causes infections in immunocompromised individuals and in people with cystic fibrosis. A low outer membrane permeability is one of the factors giving it a high level of innate resistance to antibiotics. Porins provide selective permeability for nutrients, and antibiotics can also traverse the outer membrane by this means. Knowing the properties and specificities of porin channels is therefore important for understanding resistance mechanisms and for developing new antibiotics and could help in overcoming permeability issues in antibiotic treatment.

Web Services for RNA-RNA Interaction Prediction

Non-coding RNAs have various biological functions such as translational regulation, and RNA-RNA interactions play essential roles in the mechanisms of action of these RNAs. Therefore, RNA-RNA interaction prediction is an important problem in bioinformatics, and many tools have been developed for the computational prediction of RNA-RNA interactions. In addition to the development of novel algorithms with high accuracy, the development and maintenance of web services is essential for enhancing usability by experimental biologists. In this review, we survey web services for RNA-RNA interaction predictions and introduce how to use primary web services. We present various prediction tools, including general interaction prediction tools, prediction tools for specific RNA classes, and RNA-RNA interaction-based RNA design tools. Additionally, we discuss the future perspectives of the development of RNA-RNA interaction prediction tools and the sustainability of web services.

The heme-responsive PrrH sRNA regulates Pseudomonas aeruginosa pyochelin gene expression

Pseudomonas aeruginosa is an opportunistic pathogen that requires iron for growth and virulence, yet this nutrient is sequestered by the innate immune system during infection. When iron is limiting, P. aeruginosa expresses the PrrF1 and PrrF2 small regulatory RNAs (sRNAs), which post-transcriptionally repress expression of non-essential iron-containing proteins thus sparing this nutrient for more critical processes. The genes for the PrrF1 and PrrF2 sRNAs are arranged in tandem on the chromosome, allowing for the transcription of a longer heme-responsive sRNA, termed PrrH. While the functions of PrrF1 and PrrF2 have been studied extensively, the role of PrrH in P. aeruginosa physiology and virulence is not well understood. In this study, we performed transcriptomic and proteomic studies to identify the PrrH regulon. In shaking cultures, the pyochelin synthesis proteins were increased in two distinct prrH mutants compared to wild type, while the mRNAs for these proteins were not affected by prrH mutation. We identified complementarity between the PrrH sRNA and sequence upstream of the pchE mRNA, suggesting potential for PrrH to directly regulate expression of genes for pyochelin synthesis. We further showed that pchE mRNA levels were increased in the prrH mutants when grown in static but not shaking conditions. Moreover, we discovered controlling for the presence of light was critical for examining the impact of PrrH on pchE expression. As such, our study reports on the first likely target of the PrrH sRNA and highlights key environmental variables that will allow for future characterization of PrrH function.

RyhB Paralogs Downregulate the Expressions of Multiple Survival-Associated Genes and Attenuate the Survival of Salmonella Enteritidis in the Chicken Macrophage HD11

RyhB-1 and RyhB-2 are small non-coding RNAs in Salmonella that act as regulators of iron homeostasis by sensing the environmental iron concentration. Expressions of RyhB paralogs from Salmonella Typhimurium are increased within microphages. RyhB paralogs restrain the growth of S. Typhimurium in RAW264.7 macrophages by modulating the expression of Salmonella pathogenicity island 1 (SPI-1) genes sicA and rtsB. However, little is known about the regulatory role of RyhBs and their virulence-associated targets in Salmonella Enteritidis. We studied candidate targets of RyhB paralogs via RNA-Seq in conditions of iron limitation and hypoxia. RyhB paralogs were expressed when the S. Enteritidis strain CMCC(B)50336 (SE50336) interacted with the chicken macrophage line HD11. We analyzed gene expression associated with Salmonella survival and replication in macrophages in wild-type strain SE50336 and the RyhB deletion mutants after co-incubation with HD11 and screened out targets regulated by RyhBs. The expressions of both RyhB-1 and RyhB-2 were increased after co-incubation with HD11 for 8 h and several survival-associated genes within macrophages, such as ssaI, sseA, pagC, sodC, mgtC, yaeB, pocR, and hns, were upregulated in the ryhB-1 deletion mutant. Specifically, ssaI, the type-three secretion system 2 (T3SS-2) effector encoded by SPI-2, which promoted the survival of Salmonella in macrophages, was upregulated more than 3-fold in the ryhB-1 deletion mutant. We confirmed that both RyhB-1 and RyhB-2 downregulated the expression of ssaI to repress its mRNA translation by directly interacting with its coding sequence (CDS) region via an incomplete complementary base-pairing mechanism. The SPI-2 gene sseA was indirectly modulated by RyhB-1. The survival assays in macrophages showed that the ability of intracellular survival of ryhB-1 and/or ryhB-2 deletion mutants in HD11 was higher than that of the wild-type strain. These results indicate that RyhB paralogs downregulate survival-related virulence factors and attenuate the survival of S. Enteritidis inside chicken macrophage HD11.

Hfq protein and GcvB small RNA tailoring of oppA target mRNA to levels allowing translation activation by MicF small RNA in Escherichia coli

Traffic of molecules across the bacterial membrane mainly relies on porins and transporters, whose expression must adapt to environmental conditions. To ensure bacterial fitness, synthesis and assembly of functional porins and transporters are regulated through a plethora of mechanisms. Among them, small regulatory RNAs (sRNAs) are known to be powerful post-transcriptional regulators. In Escherichia coli, the MicF sRNA is known to regulate only four targets, a very narrow targetome for a sRNA responding to various stresses, such as membrane stress, osmotic shock, or thermal shock. Using an in vivo pull-down assay combined with high-throughput RNA sequencing, we sought to identify new targets of MicF to better understand its role in the maintenance of cellular homoeostasis. Here, we report the first positively regulated target of MicF, the oppA mRNA. The OppA protein is the periplasmic component of the Opp ATP-binding cassette (ABC) oligopeptide transporter and regulates the import of short peptides, some of them bactericides. Mechanistic studies suggest that oppA translation is activated by MicF through a mechanism of action involving facilitated access to a translation-enhancing region in oppA 5'UTR. Intriguingly, MicF activation of oppA translation depends on cross-regulation by negative trans-acting effectors, the GcvB sRNA and the RNA chaperone protein Hfq.

Small RNA Targets: Advances in Prediction Tools and High-Throughput Profiling

MicroRNAs (miRNAs) are an abundant class of small non-coding RNAs that regulate gene expression at the post-transcriptional level. They are suggested to be involved in most biological processes of the cell primarily by targeting messenger RNAs (mRNAs) for cleavage or translational repression. Their binding to their target sites is mediated by the Argonaute (AGO) family of proteins. Thus, miRNA target prediction is pivotal for research and clinical applications. Moreover, transfer-RNA-derived fragments (tRFs) and other types of small RNAs have been found to be potent regulators of Ago-mediated gene expression. Their role in mRNA regulation is still to be fully elucidated, and advancements in the computational prediction of their targets are in their infancy. To shed light on these complex RNA-RNA interactions, the availability of good quality high-throughput data and reliable computational methods is of utmost importance. Even though the arsenal of computational approaches in the field has been enriched in the last decade, there is still a degree of discrepancy between the results they yield. This review offers an overview of the relevant advancements in the field of bioinformatics and machine learning and summarizes the key strategies utilized for small RNA target prediction. Furthermore, we report the recent development of high-throughput sequencing technologies, and explore the role of non-miRNA AGO driver sequences.

Global profiling of the RNA and protein complexes of Escherichia coli by size exclusion chromatography followed by RNA sequencing and mass spectrometry (SEC-seq)

New methods for the global identification of RNA-protein interactions have led to greater recognition of the abundance and importance of RNA-binding proteins (RBPs) in bacteria. Here, we expand this tool kit by developing SEC-seq, a method based on a similar concept as the established Grad-seq approach. In Grad-seq, cellular RNA and protein complexes of a bacterium of interest are separated in a glycerol gradient, followed by high-throughput RNA-sequencing and mass spectrometry analyses of individual gradient fractions. New RNA-protein complexes are predicted based on the similarity of their elution profiles. In SEC-seq, we have replaced the glycerol gradient with separation by size exclusion chromatography, which shortens operation times and offers greater potential for automation. Applying SEC-seq to Escherichia coli, we find that the method provides a higher resolution than Grad-seq in the lower molecular weight range up to ~500 kDa. This is illustrated by the ability of SEC-seq to resolve two distinct, but similarly sized complexes of the global translational repressor CsrA with either of its antagonistic small RNAs, CsrB and CsrC. We also characterized changes in the SEC-seq profiles of the small RNA MicA upon deletion of its RNA chaperones Hfq and ProQ and investigated the redistribution of these two proteins upon RNase treatment. Overall, we demonstrate that SEC-seq is a tractable and reproducible method for the global profiling of bacterial RNA-protein complexes that offers the potential to discover yet-unrecognized associations between bacterial RNAs and proteins.

sRNA-controlled iron sparing response in Staphylococci

Staphylococcus aureus, a human opportunist pathogen, adjusts its metabolism to cope with iron deprivation within the host. We investigated the potential role of small non-coding RNAs (sRNAs) in dictating this process. A single sRNA, named here IsrR, emerged from a competition assay with tagged-mutant libraries as being required during iron starvation. IsrR is iron-repressed and predicted to target mRNAs expressing iron-containing enzymes. Among them, we demonstrated that IsrR down-regulates the translation of mRNAs of enzymes that catalyze anaerobic nitrate respiration. The IsrR sequence reveals three single-stranded C-rich regions (CRRs). Mutational and structural analysis indicated a differential contribution of these CRRs according to targets. We also report that IsrR is required for full lethality of S. aureus in a mouse septicemia model, underscoring its role as a major contributor to the iron-sparing response for bacterial survival during infection. IsrR is conserved among staphylococci, but it is not ortholog to the proteobacterial sRNA RyhB, nor to other characterized sRNAs down-regulating mRNAs of iron-containing enzymes. Remarkably, these distinct sRNAs regulate common targets, illustrating that RNA-based regulation provides optimal evolutionary solutions to improve bacterial fitness when iron is scarce.

Thirty Years of sRNA-Mediated Regulation in Staphylococcus aureus: From Initial Discoveries to In Vivo Biological Implications

Staphylococcus aureus is a widespread livestock and human pathogen that colonizes diverse microenvironments within its host. Its adaptation to the environmental conditions encountered within humans relies on coordinated gene expression. This requires a sophisticated regulatory network, among which regulatory RNAs (usually called sRNAs) have emerged as key players over the last 30 years. In S. aureus, sRNAs regulate target genes at the post-transcriptional level through base–pair interactions. The functional characterization of a subset revealed that they participate in all biological processes, including virulence, metabolic adaptation, and antibiotic resistance. In this review, we report 30 years of S. aureus sRNA studies, from their discovery to the in-depth characterizations of some of them. We also discuss their actual in vivo contribution, which is still lagging behind, and their place within the complex regulatory network. These shall be key aspects to consider in order to clearly uncover their in vivo biological functions.

Synthesis of the NarP response regulator of nitrate respiration in Escherichia coli is regulated at multiple levels by Hfq and small RNAs

Two-component systems (TCS) and small RNAs (sRNA) are widespread regulators that participate in the response and the adaptation of bacteria to their environments. TCSs and sRNAs mostly act at the transcriptional and post-transcriptional levels, respectively, and can be found integrated in regulatory circuits, where TCSs control sRNAs transcription and/or sRNAs post-transcriptionally regulate TCSs synthesis. In response to nitrate and nitrite, the paralogous NarQ-NarP and NarX-NarL TCSs regulate the expression of genes involved in anaerobic respiration of these alternative electron acceptors to oxygen. In addition to the previously reported repression of NarP synthesis by the SdsN137 sRNA, we show here that RprA, another Hfq-dependent sRNA, also negatively controls narP. Interestingly, the repression of narP by RprA actually relies on two independent mechanisms of control. The first is via the direct pairing of the central region of RprA to the narP translation initiation region and presumably occurs at the translation initiation level. In contrast, the second requires only the very 5' end of the narP mRNA, which is targeted, most likely indirectly, by the full-length or the shorter, processed, form of RprA. In addition, our results raise the possibility of a direct role of Hfq in narP control, further illustrating the diversity of post-transcriptional regulation mechanisms in the synthesis of TCSs.

A Sporulation-Specific sRNA Bvs196 Contributing to the Developing Spore in Bacillus velezensis

Many putative sRNAs have been characterized using bioinformatic analysis and high-throughput sequencing in Gram-positive Bacillus strains, but there are only a few functional studies on the sRNAs involved in the spore formation developmental process. In particular, there is no sRNA confirmed experimentally to regulate the late stages of sporulation. Bvs196 is an sRNA with a length of 294 nucleotides that is abundantly expressed in the stationary phase of several media and independently transcribed in Bacillus velezensis strain PEBA20, as validated by RNA-seq and Northern blot,. It is also confirmed, by qRT-PCR, that Bvs196 is transcribed abundantly throughout the intermediate and late stages of sporulation. Using the gfpmut3a gene transcriptional reporter demonstrates that Bvs196 is expressed specifically in the forespore during sporulation and controlled by σF and σG (mainly by σG). This was observed by fluorescence microscopy and multi-function microplate reader. Further evolutionary conservation analysis found that Bvs196 is widely present in Bacillus with a strongly conserved and stable secondary structure. Resistance phenotypic assays of spores formed from the Bvs196 deletion mutant, the overexpressed Bvs196 mutant, and the wild-type strain revealed that the absence of Bvs196 led to reduced heat and UV resistance and enhanced formaldehyde resistance. We determined, by MST analysis, that Bvs196 can directly interact with spo0A and sspN-tlp mRNAs in vitro, and that short incomplete complementary paired bases affect the binding affinity of Bvs196 to target mRNAs. Our results suggest that Bvs196 is a novel sporulation-specific sRNA of B. velezensis, 294 nt in length, independently transcribed under the control of σF and σG in the forespore during sporulation, and that it affects spore resistance, and is able to directly interact with spo0A and sspN-tlp mRNAs. The remarkable conservation and impressive expression level of Bvs196 imply that it acts as an important conservative regulator, presumably by interacting with many other unknown targets in the forespore, and therefore contributing to spore properties. This work provides new clues for further understanding of the spore formation regulatory network.

An overview of gene regulation in bacteria by small RNAs derived from mRNA 3' ends

Over the past two decades, small noncoding RNAs (sRNAs) that regulate mRNAs by short base pairing have gone from a curiosity to a major class of post-transcriptional regulators in bacteria. They are integral to many stress responses and regulatory circuits, affecting almost all aspects of bacterial life. Following pioneering sRNA searches in the early 2000s, the field quickly focused on conserved sRNA genes in the intergenic regions of bacterial chromosomes. Yet, it soon emerged that there might be another rich source of bacterial sRNAs—processed 3′ end fragments of mRNAs. Several such 3′ end-derived sRNAs have now been characterized, often revealing unexpected, conserved functions in diverse cellular processes. Here, we review our current knowledge of these 3′ end-derived sRNAs—their biogenesis through ribonucleases, their molecular mechanisms, their interactions with RNA-binding proteins such as Hfq or ProQ and their functional scope, which ranges from acting as specialized regulators of single metabolic genes to constituting entire noncoding arms in global stress responses. Recent global RNA interactome studies suggest that the importance of functional 3′ end-derived sRNAs has been vastly underestimated and that this type of cross-regulation between genes at the mRNA level is more pervasive in bacteria than currently appreciated.

Inferring the contribution of small RNAs to changes in gene expression in response to stress

A main strategy of bacteria adapting to environmental changes is the remodeling of their transcriptome. Changes in the transcript levels of specific genes are due to combined effects of various regulators, including small RNAs (sRNAs). sRNAs are post-transcriptional regulators of gene expression that mainly control translation, but also directly and indirectly affect the levels of their target transcripts. Yet, the relative contribution of an sRNA to the total change in the transcript level of a gene upon an environmental change has not been assessed. We present a design of differential gene expression analysis by RNA-seq that allows extracting the contribution of an sRNA to the total change in the transcript level of each gene in response to an environmental change by fitting a linear model to the data. We exemplify this for the sRNA RyhB in cells growing under iron limitation and show a variation among genes in the relative contribution of RyhB to the change in their transcript level upon iron limitation, from subtle to very substantial. Extracting the relative contribution of an sRNA to the total change in expression of genes is important for understanding the integration of regulation by sRNAs with other regulatory mechanisms in the cell.

A modular RNA interference system for multiplexed gene regulation

The rational design and realisation of simple-to-use genetic control elements that are modular, orthogonal and robust is essential to the construction of predictable and reliable biological systems of increasing complexity. To this effect, we introduce modular Artificial RNA interference (mARi), a rational, modular and extensible design framework that enables robust, portable and multiplexed post-transcriptional regulation of gene expression in Escherichia coli. The regulatory function of mARi was characterised in a range of relevant genetic contexts, demonstrating its independence from other genetic control elements and the gene of interest, and providing new insight into the design rules of RNA based regulation in E. coli, while a range of cellular contexts also demonstrated it to be independent of growth-phase and strain type. Importantly, the extensibility and orthogonality of mARi enables the simultaneous post-transcriptional regulation of multi-gene systems as both single-gene cassettes and poly-cistronic operons. To facilitate adoption, mARi was designed to be directly integrated into the modular BASIC DNA assembly framework. We anticipate that mARi-based genetic control within an extensible DNA assembly framework will facilitate metabolic engineering, layered genetic control, and advanced genetic circuit applications.

The sRNA NsiR4 fine-tunes arginine synthesis in the cyanobacterium Synechocystis sp. PCC 6803 by post-transcriptional regulation of PirA

As the only oxygenic phototrophs among prokaryotes, cyanobacteria employ intricate mechanisms to regulate common metabolic pathways. These mechanisms include small protein inhibitors exerting their function by protein-protein interaction with key metabolic enzymes and regulatory small RNAs (sRNAs). Here we show that the sRNA NsiR4, which is highly expressed under nitrogen limiting conditions, interacts with the mRNA of the recently described small protein PirA in the model strain Synechocystis sp. PCC 6803. In particular, NsiR4 targets the pirA 5'UTR close to the ribosome binding site. Heterologous reporter assays confirmed that this interaction interferes with pirA translation. PirA negatively impacts arginine synthesis under ammonium excess by competing with the central carbon/nitrogen regulator P_II that binds to and thereby activates the key enzyme of arginine synthesis, N-acetyl-L-glutamate-kinase (NAGK). Consistently, ectopic nsiR4 expression in Synechocystis resulted in lowered PirA accumulation in response to ammonium upshifts, which also affected intracellular arginine pools. As NsiR4 and PirA are inversely regulated by the global nitrogen transcriptional regulator NtcA, this regulatory axis enables fine tuning of arginine synthesis and conveys additional metabolic flexibility under highly fluctuating nitrogen regimes. Pairs of small protein inhibitors and of sRNAs that control the abundance of these enzyme effectors at the post-transcriptional level appear as fundamental building blocks in the regulation of primary metabolism in cyanobacteria.

RNA Structure Prediction, Analysis, and Design: An Introduction to Web-Based Tools

Understanding RNA structure has become critical in the study of RNA in their roles as mediators of biological processes. To aid in these studies, computational algorithms that utilize thermodynamics have been developed to predict RNA secondary structure. Due to the importance of intermolecular interactions, the algorithms have been expanded to determine and predict RNA-RNA hybridization. This chapter discusses popular webservers with the tools for RNA secondary structure prediction, RNA-RNA hybridization, and design. We address key features that distinguish common-functioning programs and their purposes for the interests of the user. Ultimately, we hope this review elucidates web-based tools researchers may take advantage of in their investigations of RNA structure and function.

Glutamine synthetase mRNA releases sRNA from its 3'UTR to regulate carbon/nitrogen metabolic balance in Enterobacteriaceae

Glutamine synthetase (GS) is the key enzyme of nitrogen assimilation induced under nitrogen limiting conditions. The carbon skeleton of glutamate and glutamine, 2-oxoglutarate, is supplied from the TCA cycle, but how this metabolic flow is controlled in response to nitrogen availability remains unknown. We show that the expression of the E1o component of 2-oxoglutarate dehydrogenase, SucA, is repressed under nitrogen limitation in Salmonella enterica and Escherichia coli. The repression is exerted at the post-transcriptional level by an Hfq-dependent sRNA GlnZ generated from the 3'UTR of the GS-encoding glnA mRNA. Enterobacterial GlnZ variants contain a conserved seed sequence and primarily regulate sucA through base-pairing far upstream of the translation initiation region. During growth on glutamine as the nitrogen source, the glnA 3'UTR deletion mutants expressed SucA at higher levels than the S. enterica and E. coli wild-type strains, respectively. In E. coli, the transcriptional regulator Nac also participates in the repression of sucA. Lastly, this study clarifies that the release of GlnZ from the glnA mRNA by RNase E is essential for the post-transcriptional regulation of sucA. Thus, the mRNA coordinates the two independent functions to balance the supply and demand of the fundamental metabolites.

sRNARFTarget: a fast machine-learning-based approach for transcriptome-wide sRNA target prediction

Bacterial small regulatory RNAs (sRNAs) are key regulators of gene expression in many processes related to adaptive responses. A multitude of sRNAs have been identified in many bacterial species; however, their function has yet to be elucidated. A key step to understand sRNAs function is to identify the mRNAs these sRNAs bind to. There are several computational methods for sRNA target prediction, and the most accurate one is CopraRNA which is based on comparative-genomics. However, species-specific sRNAs are quite common and CopraRNA cannot be used for these sRNAs. The most commonly used transcriptome-wide sRNA target prediction method and second-most-accurate method is IntaRNA. However, IntaRNA can take hours to run on a bacterial transcriptome. Here we present sRNARFTarget, a machine-learning-based method for transcriptome-wide sRNA target prediction applicable to any sRNA. We comparatively assessed the performance of sRNARFTarget, CopraRNA and IntaRNA in three bacterial species. Our results show that sRNARFTarget outperforms IntaRNA in terms of accuracy, ranking of true interacting pairs, and running time. However, CopraRNA substantially outperforms the other two programsin terms of accuracy. Thus, we suggest using CopraRNA when homolog sequences of the sRNA are available, and sRNARFTarget for transcriptome-wide prediction or for species-specific sRNAs. sRNARFTarget is available at https://github.com/BioinformaticsLabAtMUN/sRNARFTarget.

Post-Transcriptional Regulation of RseA by Small RNAs RyhB and FnrS in Escherichia coli

RseA is the critical central regulator of the σ^E-dependent stress response in E. coli and other related bacteria. The synthesis of RseA is controlled at the transcriptional level by several promoters and transcriptional regulators, including σ^E itself at two σ^E-dependent promoters: rpoE_P and rseA_P3. The presence of these two independent polycistrons encoding rseA is potentially redundant. We hypothesized that post-transcriptional control of the rseA_P3 transcript was necessary to overcome this redundancy. However, to date, nothing is known about the post-transcriptional control of the rseA_P3 transcript. We executed a targeted genetic screen to identify small RNA regulators of the rseA_P3 transcript and identified RyhB and FnrS as small RNA activators of the RseA P3 transcript. Through genetic analysis, we confirmed that a direct interaction occurs between RyhB and RseA. We also identified sequences within the 5′ untranslated region (UTR) of RseA that were inhibitory for RseA expression. Point mutations predicted to prevent an interaction between RyhB and RseA resulted in increased RseA expression. Taken together, this suggests that the 5’ UTR of the RseAP3 transcript prevents optimal expression of RseA, preventing redundancy due to RseA expression from the σ^E-dependent rpoE_P, and this is overcome by the stimulatory activity of RyhB and FnrS.

Riboregulation in bacteria: From general principles to novel mechanisms of the trp attenuator and its sRNA and peptide products

Gene expression strategies ensuring bacterial survival and competitiveness rely on cis- and trans-acting RNA-regulators (riboregulators). Among the cis-acting riboregulators are transcriptional and translational attenuators, and antisense RNAs (asRNAs). The trans-acting riboregulators are small RNAs (sRNAs) that bind proteins or base pairs with other RNAs. This classification is artificial since some regulatory RNAs act both in cis and in trans, or function in addition as small mRNAs. A prominent example is the archetypical, ribosome-dependent attenuator of tryptophan (Trp) biosynthesis genes. It responds by transcription attenuation to two signals, Trp availability and inhibition of translation, and gives rise to two trans-acting products, the attenuator sRNA rnTrpL and the leader peptide peTrpL. In Escherichia coli, rnTrpL links Trp availability to initiation of chromosome replication and in Sinorhizobium meliloti, it coordinates regulation of split tryptophan biosynthesis operons. Furthermore, in S. meliloti, peTrpL is involved in mRNA destabilization in response to antibiotic exposure. It forms two types of asRNA-containing, antibiotic-dependent ribonucleoprotein complexes (ARNPs), one of them changing the target specificity of rnTrpL. The posttranscriptional role of peTrpL indicates two emerging paradigms: (1) sRNA reprograming by small molecules and (2) direct involvement of antibiotics in regulatory RNPs. They broaden our view on RNA-based mechanisms and may inspire new approaches for studying, detecting, and using antibacterial compounds. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

A new role for SR1 from Bacillus subtilis: regulation of sporulation by inhibition of kinA translation

SR1 is a dual-function sRNA from Bacillus subtilis. It inhibits translation initiation of ahrC mRNA encoding the transcription activator of the arginine catabolic operons. Base-pairing is promoted by the RNA chaperone CsrA, which induces a slight structural change in the ahrC mRNA to facilitate SR1 binding. Additionally, SR1 encodes the small protein SR1P that interacts with glyceraldehyde-3P dehydrogenase A to promote binding to RNase J1 and enhancing J1 activity. Here, we describe a new target of SR1, kinA mRNA encoding the major histidine kinase of the sporulation phosphorelay. SR1 and kinA mRNA share 7 complementary regions. Base-pairing between SR1 and kinA mRNA decreases kinA translation without affecting kinA mRNA stability and represses transcription of the KinA/Spo0A downstream targets spoIIE, spoIIGA and cotA. The initial interaction between SR1 and kinA mRNA occurs 10 nt downstream of the kinA start codon and is decisive for inhibition. The sr1 encoded peptide SR1P is dispensable for kinA regulation. Deletion of sr1 accelerates sporulation resulting in low quality spores with reduced stress resistance and altered coat protein composition which can be compensated by sr1 overexpression. Neither CsrA nor Hfq influence sporulation or spore properties.

A nitrogen stress-inducible small RNA regulates CO2 fixation in Nostoc

In the absence of fixed nitrogen, some filamentous cyanobacteria differentiate heterocysts, specialized cells devoted to fixing atmospheric nitrogen (N2). This differentiation process is controlled by the global nitrogen regulator NtcA and involves extensive metabolic reprogramming, including shutdown of photosynthetic CO2 fixation in heterocysts, to provide a microaerobic environment suitable for N2 fixation. Small regulatory RNAs (sRNAs) are major post-transcriptional regulators of gene expression in bacteria. In cyanobacteria, responding to nitrogen deficiency involves transcribing several nitrogen-regulated sRNAs. Here, we describe the participation of nitrogen stress-inducible RNA 4 (NsiR4) in post-transcriptionally regulating the expression of two genes involved in CO2 fixation via the Calvin cycle: glpX, which encodes bifunctional sedoheptulose-1,7-bisphosphatase/fructose-1,6-bisphosphatase (SBPase), and pgk, which encodes phosphoglycerate kinase (PGK). Using a heterologous reporter assay in Escherichia coli, we show that NsiR4 interacts with the 5'-untranslated region (5'-UTR) of glpX and pgk mRNAs. Overexpressing NsiR4 in Nostoc sp. PCC 7120 resulted in a reduced amount of SBPase protein and reduced PGK activity, as well as reduced levels of both glpX and pgk mRNAs, further supporting that NsiR4 negatively regulates these two enzymes. In addition, using a gfp fusion to the nsiR4 promoter, we show stronger expression of NsiR4 in heterocysts than in vegetative cells, which could contribute to the heterocyst-specific shutdown of Calvin cycle flux. Post-transcriptional regulation of two Calvin cycle enzymes by NsiR4, a nitrogen-regulated sRNA, represents an additional link between nitrogen control and CO2 assimilation.

Involvement of non-coding RNAs during infection of rice by Acidovorax oryzae

The expression of non-coding RNAs (ncRNAs) has been observed in a variety of bacteria. However, the function of ncRNAs and their regulatory targets are largely unknown, and few ncRNAs are found to be associated with bacterial virulence. The bacterial brown stripe pathogen Acidovorax oryzae (Ao) RS-1 shows a high level of condition-dependent differential expression of ncRNA, which we identified in a genome wide screen. We experimentally validated 66 differentially expressed ncRNAs using an integrative analysis of conservative genome sequences and transcriptomic data during in vivo interaction of the bacterial pathogen with the rice plant. To test the relevance of the differentially expressed ncRNAs, we chose four with different positions within the genome, and with different secondary structures and promoter activities. The results show that the overexpression of the four ncRNAs caused a significant change in virulence-related phenotypes, resistance to various environmental stresses, expression of secretion systems and effector proteins, while changing the expression of ncRNA putative target genes. We conclude that these ncRNAs are examples for the inherent regulatory roles for many of the observed ncRNAs in response to changing conditions such as host interaction or environmental adaption.

A key antisense sRNA modulates the oxidative stress response and virulence in Xanthomonas oryzae pv. oryzicola

Pathogens integrate multiple environmental signals to navigate the host and control the expression of virulence genes. In this process, small regulatory noncoding RNAs (sRNAs) may function in gene expression as post-transcriptional regulators. In this study, the sRNA Xonc3711 functioned in the response of the rice pathogen, Xanthomonas oryzae pv. oryzicola (Xoc), to oxidative stress. Xonc3711 repressed production of the DNA-binding protein Xoc_3982 by binding to the xoc_3982 mRNA within the coding region. Mutational analysis showed that regulation required an antisense interaction between Xonc3711 and xoc_3982 mRNA, and RNase E was needed for degradation of the xoc_3982 transcript. Deletion of Xonc3711 resulted in a lower tolerance to oxidative stress due to the repression of flagella-associated genes and reduced biofilm formation. Furthermore, ChIP-seq and electrophoretic mobility shift assays showed that Xoc_3982 repressed the transcription of effector xopC2, which contributes to virulence in Xoc BLS256. This study describes how sRNA Xonc3711 modulates multiple traits in Xoc via signals perceived from the external environment.

Small, stable RNA species perform diverse functions in both prokaryotes and eukaryotes. In this study, the sRNA Xonc3711 decreased the production of DNA-binding protein Xoc_3982 in the bacterium Xanthomonas oryzae pv. oryzicola (Xoc) by base pairing with the xoc_3982 transcript. When Xonc3711 was mutated, Xoc was impaired in its ability to form flagella and produce biofilms, which reduced Xoc tolerance to oxidative stress. We also discovered that the DNA-binding protein Xoc_3982 represses the expression of xopC2, which encodes an effector protein, and reduces its expression. Our results show that Xonc3711 modulates and integrates multiple systems in Xoc to protect cells from oxidative damage.

Identification of an RNA sponge that controls the RoxS riboregulator of central metabolism in Bacillus subtilis

sRNAs are a taxonomically-restricted but transcriptomically-abundant class of post-transcriptional regulators. While of major importance for adaption to the environment, we currently lack global-scale methodology enabling target identification, especially in species without known RNA hub proteins (e.g. Hfq). Using psoralen RNA cross-linking and Illumina-sequencing we identify RNA-RNA interacting pairs in vivo in Bacillus subtilis, resolving previously well-described interactants. Although sRNA-sRNA pairings are rare (compared with sRNA-mRNA), we identify a robust example involving the conserved sRNA RoxS and an unstudied sRNA RosA (Regulator of sRNA A). We show RosA to be the first confirmed RNA sponge described in a Gram-positive bacterium. RosA interacts with at least two sRNAs, RoxS and FsrA. The RosA/RoxS interaction not only affects the levels of RoxS but also its processing and regulatory activity. We also found that the transcription of RosA is repressed by CcpA, the key regulator of carbon-metabolism in B. subtilis. Since RoxS is already known to be transcriptionally controlled by malate via the transcriptional repressor Rex, its post-transcriptional regulation by CcpA via RosA places RoxS in a key position to control central metabolism in response to varying carbon sources.

An Integrated Database of Small RNAs and Their Interplay With Transcriptional Gene Regulatory Networks in Corynebacteria

Small RNAs (sRNAs) are one of the key players in the post-transcriptional regulation of bacterial gene expression. These molecules, together with transcription factors, form regulatory networks and greatly influence the bacterial regulatory landscape. Little is known concerning sRNAs and their influence on the regulatory machinery in the genus Corynebacterium, despite its medical, veterinary and biotechnological importance. Here, we expand corynebacterial regulatory knowledge by integrating sRNAs and their regulatory interactions into the transcriptional regulatory networks of six corynebacterial species, covering four human and animal pathogens, and integrate this data into the CoryneRegNet database. To this end, we predicted sRNAs to regulate 754 genes, including 206 transcription factors, in corynebacterial gene regulatory networks. Amongst them, the sRNA Cd-NCTC13129-sRNA-2 is predicted to directly regulate ydfH, which indirectly regulates 66 genes, including the global regulator glxR in C. diphtheriae. All of the sRNA-enriched regulatory networks of the genus Corynebacterium have been made publicly available in the newest release of CoryneRegNet(www.exbio.wzw.tum.de/coryneregnet/) to aid in providing valuable insights and to guide future experiments.

Transcriptional and Post-transcriptional Control of the Nitrate Respiration in Bacteria

In oxygen (O2) limiting environments, numerous aerobic bacteria have the ability to shift from aerobic to anaerobic respiration to release energy. This process requires alternative electron acceptor to replace O2 such as nitrate (NO3 -), which has the next best reduction potential after O2. Depending on the organism, nitrate respiration involves different enzymes to convert NO3 - to ammonium (NH4 +) or dinitrogen (N2). The expression of these enzymes is tightly controlled by transcription factors (TFs). More recently, bacterial small regulatory RNAs (sRNAs), which are important regulators of the rapid adaptation of microorganisms to extremely diverse environments, have also been shown to control the expression of genes encoding enzymes or TFs related to nitrate respiration. In turn, these TFs control the synthesis of multiple sRNAs. These results suggest that sRNAs play a central role in the control of these metabolic pathways. Here we review the complex interplay between the transcriptional and the post-transcriptional regulators to efficiently control the respiration on nitrate.

Positive regulation of Type III secretion effectors and virulence by RyhB paralogs in Salmonella enterica serovar Enteritidis

Small non-coding RNA RyhB is a key regulator of iron homeostasis in bacteria by sensing iron availability in the environment. Although RyhB is known to influence bacterial virulence by interacting with iron metabolism related regulators, its interaction with virulence genes, especially the Type III secretion system (T3SS), has not been reported. Here, we demonstrate that two RyhB paralogs of Salmonella enterica serovar Enteritidis upregulate Type III secretion system (T3SS) effectors, and consequently affect Salmonella invasion into intestinal epithelial cells. Specifically, we found that RyhB-1 modulate Salmonella response to stress condition of iron deficiency and hypoxia, and stress in simulated intestinal environment (SIE). Under SIE culture conditions, both RyhB-1 and RyhB-2 are drastically induced and directly upregulate the expression of T3SS effector gene sipA by interacting with its 5′ untranslated region (5′ UTR) via an incomplete base-pairing mechanism. In addition, the RyhB paralogs upregulate the expression of T3SS effector gene sopE. By regulating the invasion-related genes, RyhBs in turn affect the ability of S. Enteritidis to adhere to and invade into intestinal epithelial cells. Our findings provide evidence that RyhBs function as critical virulence factors by directly regulating virulence-related gene expression. Thus, inhibition of RyhBs may be a potential strategy to attenuate Salmonella.

Bioinformatics as a Tool for the Structural and Evolutionary Analysis of Proteins

This chapter deals with the topic of bioinformatics, computational, mathematics, and statistics tools applied to biology, essential for the analysis and characterization of biological molecules, in particular proteins, which play an important role in all cellular and evolutionary processes of the organisms. In recent decades, with the next generation sequencing technologies and bioinformatics, it has facilitated the collection and analysis of a large amount of genomic, transcriptomic, proteomic, and metabolomic data from different organisms that have allowed predictions on the regulation of expression, transcription, translation, structure, and mechanisms of action of proteins as well as homology, mutations, and evolutionary processes that generate structural and functional changes over time. Although the information in the databases is greater every day, all bioinformatics tools continue to be constantly modified to improve performance that leads to more accurate predictions regarding protein functionality, which is why bioinformatics research remains a great challenge.

Bioinformatic prediction reveals posttranscriptional regulation of the chromosomal replication initiator gene dnaA by the attenuator sRNA rnTrpL in Escherichia coli

DnaA is the initiator protein of chromosome replication, but the regulation of its homoeostasis in enterobacteria is not well understood. The DnaA level remains stable at different growth rates, suggesting a link between metabolism and dnaA expression. In a bioinformatic prediction, which we made to unravel targets of the sRNA rnTrpL in Enterobacteriaceae, the dnaA mRNA was the most conserved target candidate. The sRNA rnTrpL is derived from the transcription attenuator of the tryptophan biosynthesis operon. In Escherichia coli, its level is higher in minimal than in rich medium due to derepressed transcription without external tryptophan supply. Overexpression and deletion of the rnTrpL gene decreased and increased, respectively, the levels of dnaA mRNA. The decrease of the dnaA mRNA level upon rnTrpL overproduction was dependent on hfq and rne. Base pairing between rnTrpL and dnaA mRNA in vivo was validated. In minimal medium, the oriC level was increased in the ΔtrpL mutant, in line with the expected DnaA overproduction and increased initiation of chromosome replication. In line with this, chromosomal rnTrpL mutation abolishing the interaction with dnaA increased both the dnaA mRNA and the oriC level. Moreover, upon addition of tryptophan to minimal medium cultures, the oriC level in the wild type was increased. Thus, rnTrpL is a base-pairing sRNA that posttranscriptionally regulates dnaA in E. coli. Furthermore, our data suggest that rnTrpL contributes to the DnaA homoeostasis in dependence on the nutrient availability, which is represented by the tryptophan level in the cell.

The mRNA derived MalH sRNA contributes to alternative carbon source utilization by tuning maltoporin expression in E. coli

Previous high-throughput studies in Gram-negative bacteria identified a large number of 3ʹUTR fragments that potentially function as sRNAs. Here we extensively characterize the MalH sRNA. We show that MalH is a stable degradation intermediate derived from the 3ʹ end of malG, which is part of the maltose uptake operon transcript malEFG. Unlike the majority of bacterial sRNAs, MalH is transiently expressed during the transition from the exponential to the stationary growth phase, suggesting that it contributes to adaptation to changes in nutrient availability. Over-expression of MalH reduces expression of general outer membrane porins and MicA, a repressor of the high-affinity maltose/maltodextrin transporter LamB. Disrupting MalH production and function significantly reduces lamB accumulation when maltose is the only available carbon source, presumably due to the accumulation of the MicA repressor. We propose that MalH is part of a regulatory network that, during the transition phase, directly or indirectly promotes accumulation of high-affinity maltose transporters in the outer membrane by dampening competing pathways.

NsiR3, a nitrogen stress-inducible small RNA, regulates proline oxidase expression in the cyanobacterium Nostoc sp. PCC 7120

NsiR3 (nitrogen stress-inducible RNA 3) is a small noncoding RNA strongly conserved in heterocyst-forming cyanobacteria. In Nostoc sp. PCC 7120, transcription of NsiR3 is induced by nitrogen starvation and depends on the global nitrogen regulator NtcA. A conserved NtcA-binding site is centered around position -42.5 with respect to the transcription start site of NsiR3 homologs, and NtcA binds in vitro to a DNA fragment containing this sequence. In the absence of combined nitrogen, NsiR3 expression is induced in all cells along the Nostoc filament but much more strongly in heterocysts, differentiated cells devoted to nitrogen fixation. Co-expression analysis of transcriptomic data obtained from microarrays hybridized with RNA obtained from Nostoc wild-type or mutant strains grown in the presence of ammonium or in the absence of combined nitrogen revealed that the expression profile of gene putA (proline oxidase) correlates negatively with that of NsiR3. Using a heterologous system in Escherichia coli, we show that NsiR3 binds to the 5'-UTR of putA mRNA, resulting in reduced expression of a reporter gene. Overexpression of NsiR3 in Nostoc resulted in strong reduction of putA mRNA accumulation, further supporting the negative regulation of putA by NsiR3. The higher expression of NsiR3 in heterocysts versus vegetative cells of the N₂ -fixing filament could contribute to the previously described absence of putA mRNA and of the catabolic pathway to produce glutamate from arginine via proline specifically in heterocysts. Post-transcriptional regulation by NsiR3 represents an indirect NtcA-operated regulatory mechanism of putA expression. DATABASE: Microarray data are available in GEO database under accession numbers GSE120377 and GSE150191.

A framework for the computational prediction and analysis of non-coding RNAs in microbial environmental populations and their experimental validation

Small regulatory RNAs and antisense RNAs play important roles in the regulation of gene expression in bacteria but are underexplored, especially in natural populations. While environmentally relevant microbes often are not amenable to genetic manipulation or cannot be cultivated in the laboratory, extensive metagenomic and metatranscriptomic datasets for these organisms might be available. Hence, dedicated workflows for specific analyses are needed to fully benefit from this information. Here, we identified abundant sRNAs from oceanic environmental populations of the ecologically important primary producer Prochlorococcus starting from a metatranscriptomic differential RNA-Seq (mdRNA-Seq) dataset. We tracked their homologs in laboratory isolates, and we provide a framework for their further detailed characterization. Several of the experimentally validated sRNAs responded to ecologically relevant changes in cultivation conditions. The expression of the here newly discovered sRNA Yfr28 was highly stimulated in low-nitrogen conditions. Its predicted top targets include mRNAs encoding cell division proteins, a sigma factor, and several enzymes and transporters, suggesting a pivotal role of Yfr28 in the coordination of primary metabolism and cell division. A cis-encoded antisense RNA was identified as a possible positive regulator of atpF encoding subunit b' of the ATP synthase complex. The presented workflow will also be useful for other environmentally relevant microorganisms for which experimental validation abilities are frequently limiting although there is wealth of sequence information available.

Low iron-induced small RNA BrrF regulates central metabolism and oxidative stress responses in Burkholderia cenocepacia

Regulatory small RNAs play an essential role in maintaining cell homeostasis in bacteria in response to environmental stresses such as iron starvation. Prokaryotes generally encode a large number of RNA regulators, yet their identification and characterisation is still in its infancy for most bacterial species. Burkholderia cenocepacia is an opportunistic pathogen with high innate antimicrobial resistance, which can cause the often fatal cepacia syndrome in individuals with cystic fibrosis. In this study we characterise a small RNA which is involved in the response to iron starvation, a condition that pathogenic bacteria are likely to encounter in the host. BrrF is a small RNA highly upregulated in Burkholderia cenocepacia under conditions of iron depletion and with a genome context consistent with Fur regulation. Its computationally predicted targets include iron-containing enzymes of the tricarboxylic acid (TCA) cycle such as aconitase and succinate dehydrogenase, as well as iron-containing enzymes responsible for the oxidative stress response, such as superoxide dismutase and catalase. Phenotypic and gene expression analysis of BrrF deletion and overexpression mutants show that the regulation of these genes is BrrF-dependent. Expression of acnA, fumA, sdhA and sdhC was downregulated during iron depletion in the wild type strain, but not in a BrrF deletion mutant. TCA cycle genes not predicted as target for BrrF were not affected in the same manner by iron depletion. Likewise, expression of sodB and katB was dowregulated during iron depletion in the wild type strain, but not in a BrrF deletion mutant. BrrF overexpression reduced aconitase and superoxide dismutase activities and increased sensitivity to hydrogen peroxide. All phenotypes and gene expression changes of the BrrF deletion mutant could be complemented by overexpressing BrrF in trans. Overall, BrrF acts as a regulator of central metabolism and oxidative stress response, possibly as an iron-sparing measure to maintain iron homeostasis under conditions of iron starvation.

Plugging Small RNAs into the Network

Small RNAs (sRNAs) have been discovered in every bacterium examined and have been shown to play important roles in the regulation of a diverse range of behaviors, from metabolism to infection. However, despite a wide range of available techniques for discovering and validating sRNA regulatory interactions, only a minority of these molecules have been well characterized. In part, this is due to the nature of posttranscriptional regulation: the activity of an sRNA depends on the state of the transcriptome as a whole, so characterization is best carried out under the conditions in which it is naturally active. In this issue of mSystems, Arrieta-Ortiz and colleagues (M. L. Arrieta-Ortiz, C. Hafemeister, B. Shuster, N. S. Baliga, et al., mSystems 5:e00057-20, 2020, https://doi.org/10.1128/mSystems.00057-20) present a network inference approach based on estimating sRNA activity across transcriptomic compendia. This shows promise not only for identifying new sRNA regulatory interactions but also for pinpointing the conditions in which these interactions occur, providing a new avenue toward functional characterization of sRNAs.

Inference of Bacterial Small RNA Regulatory Networks and Integration with Transcription Factor-Driven Regulatory Networks

Small noncoding RNAs (sRNAs) are key regulators of bacterial gene expression. Through complementary base pairing, sRNAs affect mRNA stability and translation efficiency. Here, we describe a network inference approach designed to identify sRNA-mediated regulation of transcript levels. We use existing transcriptional data sets and prior knowledge to infer sRNA regulons using our network inference tool, the Inferelator. This approach produces genome-wide gene regulatory networks that include contributions by both transcription factors and sRNAs. We show the benefits of estimating and incorporating sRNA activities into network inference pipelines using available experimental data. We also demonstrate how these estimated sRNA regulatory activities can be mined to identify the experimental conditions where sRNAs are most active. We uncover 45 novel experimentally supported sRNA-mRNA interactions in Escherichia coli, outperforming previous network-based efforts. Additionally, our pipeline complements sequence-based sRNA-mRNA interaction prediction methods by adding a data-driven filtering step. Finally, we show the general applicability of our approach by identifying 24 novel, experimentally supported, sRNA-mRNA interactions in Pseudomonas aeruginosa, Staphylococcus aureus, and Bacillus subtilis. Overall, our strategy generates novel insights into the functional context of sRNA regulation in multiple bacterial species.

IMPORTANCE Individual bacterial genomes can have dozens of small noncoding RNAs with largely unexplored regulatory functions. Although bacterial sRNAs influence a wide range of biological processes, including antibiotic resistance and pathogenicity, our current understanding of sRNA-mediated regulation is far from complete. Most of the available information is restricted to a few well-studied bacterial species; and even in those species, only partial sets of sRNA targets have been characterized in detail. To close this information gap, we developed a computational strategy that takes advantage of available transcriptional data and knowledge about validated and putative sRNA-mRNA interactions for inferring expanded sRNA regulons. Our approach facilitates the identification of experimentally supported novel interactions while filtering out false-positive results. Due to its data-driven nature, our method prioritizes biologically relevant interactions among lists of candidate sRNA-target pairs predicted in silico from sequence analysis or derived from sRNA-mRNA binding experiments.