Regulatory small RNAs from the 3' regions of bacterial mRNAs

Interplay of two small RNAs fine-tunes hierarchical flagella gene expression in Campylobacter jejuni

Like for many bacteria, flagella are crucial for Campylobacter jejuni motility and virulence. Biogenesis of the flagellar machinery requires hierarchical transcription of early, middle (RpoN-dependent), and late (FliA-dependent) genes. However, little is known about post-transcriptional regulation of flagellar biogenesis by small RNAs (sRNAs). Here, we characterized two sRNAs with opposing effects on C. jejuni filament assembly and motility. We demonstrate that CJnc230 sRNA (FlmE), encoded downstream of the flagellar hook protein, is processed from the RpoN-dependent flgE mRNA by RNase III, RNase Y, and PNPase. We identify mRNAs encoding a flagella-interaction regulator and the anti-sigma factor FlgM as direct targets of CJnc230 repression. CJnc230 overexpression upregulates late genes, including the flagellin flaA, culminating in longer flagella and increased motility. In contrast, overexpression of the FliA-dependent sRNA CJnc170 (FlmR) reduces flagellar length and motility. Overall, our study demonstrates how the interplay of two sRNAs post-transcriptionally fine-tunes flagellar biogenesis through balancing of the hierarchically-expressed components.

A novel ionizing radiation-induced small RNA, DrsS, promotes the detoxification of reactive oxygen species in Deinococcus radiodurans

A plethora of gene regulatory mechanisms with eccentric attributes in Deinoccocus radiodurans confer it to possess a distinctive ability to survive under ionizing radiation. Among the many regulatory processes, small RNA (sRNA)-mediated regulation of gene expression is prevalent in bacteria but barely investigated in D. radiodurans. In the current study, we identified a novel sRNA, DrsS, through RNA-seq analysis in D. radiodurans cells while exposed to ionizing radiation. Initial sequence analysis for promoter identification revealed that drsS is potentially co-transcribed with sodA and dr_1280 from a single operon. Elimination of the drsS allele in D. radiodurans chromosome resulted in an impaired growth phenotype under γ-radiation. DrsS has also been found to be upregulated under oxidative and genotoxic stresses. Deletion of the drsS gene resulted in the depletion of intracellular concentration of both Mn2+ and Fe2+ by ~70% and 40%, respectively, with a concomitant increase in carbonylation of intracellular protein. Complementation of drsS gene in ΔdrsS cells helped revert its intracellular Mn2+ and Fe2+ concentration and alleviated carbonylation of intracellular proteins. Cells with deleted drsS gene exhibited higher sensitivity to oxidative stress than wild-type cells. Extrachromosomally expressed drsS in ΔdrsS cells retrieved its oxidative stress resistance properties by catalase-mediated detoxification of reactive oxygen species (ROS). In vitro binding assays indicated that DsrS directly interacts with the coding region of the katA transcript, thus possibly protecting it from cellular endonucleases in vivo. This study identified a novel small RNA DrsS and investigated its function under oxidative stress in D. radiodurans.

IMPORTANCE

Deinococcus radiodurans possesses an idiosyncratic quality to survive under extreme ionizing radiation and, thus, has evolved with diverse mechanisms which promote the mending of intracellular damages caused by ionizing radiation. As sRNAs play a pivotal role in modulating gene expression to adapt to altered conditions and have been delineated to participate in almost all physiological processes, understanding the regulatory mechanism of sRNAs will unearth many pathways that lead to radioresistance in D. radiodurans. In that direction, DrsS has been identified to be a γ-radiation-induced sRNA, which is also induced by oxidative and genotoxic stresses. DrsS appeared to activate catalase under oxidative stress and detoxify intracellular ROS. This sRNA has also been shown to balance intracellular Mn(II) and Fe concentrations protecting intracellular proteins from carbonylation. This novel mechanism of DrsS identified in D. radiodurans adds substantially to our knowledge of how this bacterium exploits sRNA for its survival under stresses.

The global RNA-RNA interactome of Klebsiella pneumoniae unveils a small RNA regulator of cell division

The ubiquitous RNA chaperone Hfq is involved in the regulation of key biological processes in many species across the bacterial kingdom. In the opportunistic human pathogen Klebsiella pneumoniae, deletion of the hfq gene affects the global transcriptome, virulence, and stress resistance; however, the ligands of the major RNA-binding protein in this species have remained elusive. In this study, we have combined transcriptomic, co-immunoprecipitation, and global RNA interactome analyses to compile an inventory of conserved and species-specific RNAs bound by Hfq and to monitor Hfq-mediated RNA-RNA interactions. In addition to dozens of RNA-RNA pairs, our study revealed an Hfq-dependent small regulatory RNA (sRNA), DinR, which is processed from the 3' terminal portion of dinI mRNA. Transcription of dinI is controlled by the master regulator of the SOS response, LexA. As DinR accumulates in K. pneumoniae in response to DNA damage, the sRNA represses translation of the ftsZ transcript by occupation of the ribosome binding site. Ectopic overexpression of DinR causes depletion of ftsZ mRNA and inhibition of cell division, while deletion of dinR antagonizes cell elongation in the presence of DNA damage. Collectively, our work highlights the important role of RNA-based gene regulation in K. pneumoniae and uncovers the central role of DinR in LexA-controlled division inhibition during the SOS response.

In vivo RNA interactome profiling reveals 3'UTR-processed small RNA targeting a central regulatory hub

Small noncoding RNAs (sRNAs) are crucial regulators of gene expression in bacteria. Acting in concert with major RNA chaperones such as Hfq or ProQ, sRNAs base-pair with multiple target mRNAs and form large RNA-RNA interaction networks. To systematically investigate the RNA-RNA interactome in living cells, we have developed a streamlined in vivo approach iRIL-seq (intracellular RIL-seq). This generic approach is highly robust, illustrating the dynamic sRNA interactomes in Salmonella enterica across multiple stages of growth. We have identified the OmpD porin mRNA as a central regulatory hub that is targeted by a dozen sRNAs, including FadZ cleaved from the conserved 3'UTR of fadBA mRNA. Both ompD and FadZ are activated by CRP, constituting a type I incoherent feed-forward loop in the fatty acid metabolism pathway. Altogether, we have established an approach to profile RNA-RNA interactomes in live cells, highlighting the complexity of RNA regulatory hubs and RNA networks.

Ribonuclease D Processes a Small RNA Regulator of Multicellular Development in Myxobacteria

By targeting mRNA transcripts, non-coding small RNAs (sRNAs) regulate the expression of genes governing a wide range of bacterial functions. In the social myxobacterium Myxococcus xanthus, the sRNA Pxr serves as a gatekeeper of the regulatory pathway controlling the life-cycle transition from vegetative growth to multicellular fruiting body development. When nutrients are abundant, Pxr prevents the initiation of the developmental program, but Pxr-mediated inhibition is alleviated when cells starve. To identify genes essential for Pxr function, a developmentally defective strain in which Pxr-mediated blockage of development is constitutively active (strain "OC") was transposon-mutagenized to identify suppressor mutations that inactivate or bypass Pxr inhibition and thereby restore development. One of the four loci in which a transposon insertion restored development is rnd, encoding the Ribonuclease D protein (RNase D). RNase D is an exonuclease important for tRNA maturation. Here, we show that disruption of rnd abolishes the accumulation of Pxr-S, the product of Pxr processing from a longer precursor form (Pxr-L) and the active inhibitor of development. Additionally, the decrease in Pxr-S caused by rnd disruption was associated with increased accumulation primarily of a longer novel Pxr-specific transcript (Pxr-XL) rather than of Pxr-L. The introduction of a plasmid expressing rnd reverted cells back to OC-like phenotypes in development and Pxr accumulation, indicating that a lack of RNase D alone suppresses the developmental defect of OC. Moreover, an in vitro Pxr-processing assay demonstrated that RNase D processes Pxr-XL into Pxr-L; this implies that overall, Pxr sRNA maturation requires a sequential two-step processing. Collectively, our results indicate that a housekeeping ribonuclease plays a central role in a model form of microbial aggregative development. To our knowledge, this is the first evidence implicating RNase D in sRNA processing.

Transcription start site mapping and small RNA profiling of Leptospira biflexa serovar Patoc

Leptospirosis is an emerging zoonotic disease caused by bacterial species of the genus Leptospira. However, the regulatory mechanisms and pathways underlying the adaptation of pathogenic and non-pathogenic Leptospira spp. in different environmental conditions remain elusive. Leptospira biflexa is a non-pathogenic species of Leptospira that lives exclusively in a natural environment. It is an ideal model not only for exploring molecular mechanisms underlying the environmental survival of Leptospira species but also for identifying virulence factors unique to Leptospira's pathogenic species. In this study, we aim to establish the transcription start site (TSS) landscape and the small RNA (sRNA) profile of L. biflexa serovar Patoc grown to exponential and stationary phases via differential RNA-seq (dRNA-seq) and small RNA-seq (sRNA-seq) analyses, respectively. Our dRNA-seq analysis uncovered a total of 2726 TSSs, which are also used to identify other elements, e.g., promoter and untranslated regions (UTRs). Besides, our sRNA-seq analysis revealed a total of 603 sRNA candidates, comprising 16 promoter-associated sRNAs, 184 5'UTR-derived sRNAs, 230 true intergenic sRNAs, 136 5'UTR-antisense sRNAs, and 130 open reading frame (ORF)-antisense sRNAs. In summary, these findings reflect the transcriptional complexity of L. biflexa serovar Patoc under different growth conditions and help to facilitate our understanding of regulatory networks in L. biflexa. To the best of our knowledge, this is the first study reporting the TSS landscape of L. biflexa. The TSS and sRNA landscapes of L. biflexa can also be compared with its pathogenic counterparts, e.g., L. borgpetersenii and L. interrogans, to identify features contributing to their environmental survival and virulence.

RNA Preparation and RNA-Seq Bioinformatics for Comparative Transcriptomics

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

Sense and antisense RNA products of the uxuR gene can affect motility and chemotaxis acting independent of the UxuR protein

Small non-coding and antisense RNAs are widespread in all kingdoms of life, however, the diversity of their functions in bacteria is largely unknown. Here, we study RNAs synthesised from divergent promoters located in the 3'-end of the uxuR gene, encoding transcription factor regulating hexuronate metabolism in Escherichia coli. These overlapping promoters were predicted in silico with rather high scores, effectively bound RNA polymerase in vitro and in vivo and were capable of initiating transcription in sense and antisense directions. The genome-wide correlation between in silico promoter scores and RNA polymerase binding in vitro and in vivo was higher for promoters located on the antisense strands of the genes, however, sense promoters within the uxuR gene were more active. Both regulatory RNAs synthesised from the divergent promoters inhibited expression of genes associated with the E. coli motility and chemotaxis independent of a carbon source on which bacteria had been grown. Direct effects of these RNAs were confirmed for the fliA gene encoding σ²⁸ subunit of RNA polymerase. In addition to intracellular sRNAs, promoters located within the uxuR gene could initiate synthesis of transcripts found in the fraction of RNAs secreted in the extracellular medium. Their profile was also carbon-independent suggesting that intragenic uxuR transcripts have a specific regulatory role not directly related to the function of the protein in which gene they are encoded.

The potential of cold-shock promoters for the expression of recombinant proteins in microbes and mammalian cells

BACKGROUND

Low-temperature expression of recombinant proteins may be advantageous to support their proper folding and preserve bioactivity. The generation of expression vectors regulated under cold conditions can improve the expression of some target proteins that are difficult to express in different expression systems. The cspA encodes the major cold-shock protein from Escherichia coli (CspA). The promoter of cspA has been widely used to develop cold shock-inducible expression platforms in E. coli. Moreover, it is often necessary to employ expression systems other than bacteria, particularly when recombinant proteins require complex post-translational modifications. Currently, there are no commercial platforms available for expressing target genes by cold shock in eukaryotic cells. Consequently, genetic elements that respond to cold shock offer the possibility of developing novel cold-inducible expression platforms, particularly suitable for yeasts, and mammalian cells.

CONCLUSIONS

This review covers the importance of the cellular response to low temperatures and the prospective use of cold-sensitive promoters to direct the expression of recombinant proteins. This concept may contribute to renewing interest in applying white technologies to produce recombinant proteins that are difficult to express.

Identification of Attenuators of Transcriptional Termination: Implications for RNA Regulation in Escherichia coli

The regulatory function of many bacterial small RNAs (sRNAs) requires the binding of the RNA chaperone Hfq to the 3' portion of the sRNA intrinsic terminator, and therefore sRNA signaling might be regulated by modulating its terminator. Here, using a multicopy screen developed with the terminator of sRNA SgrS, we identified an sRNA gene (cyaR) and three protein-coding genes (cspD, ygjH, and rof) that attenuate SgrS termination in Escherichia coli. Analyses of CyaR and YgjH, a putative tRNA binding protein, suggested that the CyaR activity was indirect and the effect of YgjH was moderate. Overproduction of the protein attenuators CspD and Rof resulted in more frequent readthrough at terminators of SgrS and two other sRNAs, and regulation by SgrS of target mRNAs was reduced. The effect of Rof, a known inhibitor of Rho, was mimicked by bicyclomycin or by a rho mutant, suggesting an unexpected role for Rho in sRNA termination. CspD, a member of the cold shock protein family, bound both terminated and readthrough transcripts, stabilizing them and attenuating termination. By RNA sequencing analysis of the CspD overexpression strain, we found global effects of CspD on gene expression across some termination sites. We further demonstrated effects of endogenous CspD under slow growth conditions where cspD is highly expressed. These findings provided evidence of changes in the efficiency of intrinsic termination, confirming this as an additional layer of the regulation of sRNA signaling. IMPORTANCE Growing evidence suggests that the modulation of intrinsic termination and readthrough of transcription is more widespread than previously appreciated. For small RNAs, proper termination plays a critical role in their regulatory function. Here, we present a multicopy screen approach to identify factors that attenuate small RNA termination and therefore abrogate signaling dependent on the small RNA. This study highlights a new aspect of regulation of small RNA signaling as well as the modulation of intrinsic termination.

High-Resolution Small RNAs Landscape Provides Insights into Alkane Adaptation in the Marine Alkane-Degrader Alcanivorax dieselolei B-5

Alkanes are widespread in the ocean, and Alcanivorax is one of the most ubiquitous alkane-degrading bacteria in the marine ecosystem. Small RNAs (sRNAs) are usually at the heart of regulatory pathways, but sRNA-mediated alkane metabolic adaptability still remains largely unknown due to the difficulties of identification. Here, differential RNA sequencing (dRNA-seq) modified with a size selection (~50-nt to 500-nt) strategy was used to generate high-resolution sRNAs profiling in the model species Alcanivorax dieselolei B-5 under alkane (n-hexadecane) and non-alkane (acetate) conditions. As a result, we identified 549 sRNA candidates at single-nucleotide resolution of 5'-ends, 63.4% of which are with transcription start sites (TSSs), and 36.6% of which are with processing sites (PSSs) at the 5'-ends. These sRNAs originate from almost any location in the genome, regardless of intragenic (65.8%), antisense (20.6%) and intergenic (6.2%) regions, and RNase E may function in the maturation of sRNAs. Most sRNAs locally distribute across the 15 reference genomes of Alcanivorax, and only 7.5% of sRNAs are broadly conserved in this genus. Expression responses to the alkane of several core conserved sRNAs, including 6S RNA, M1 RNA and tmRNA, indicate that they may participate in alkane metabolisms and result in more actively global transcription, RNA processing and stresses mitigation. Two novel CsrA-related sRNAs are identified, which may be involved in the translational activation of alkane metabolism-related genes by sequestering the global repressor CsrA. The relationships of sRNAs with the characterized genes of alkane sensing (ompS), chemotaxis (mcp, cheR, cheW2), transporting (ompT1, ompT2, ompT3) and hydroxylation (alkB1, alkB2, almA) were created based on the genome-wide predicted sRNA-mRNA interactions. Overall, the sRNA landscape lays the ground for uncovering cryptic regulations in critical marine bacterium, among which both the core and species-specific sRNAs are implicated in the alkane adaptive metabolisms.

The stationary phase-specific sRNA FimR2 is a multifunctional regulator of bacterial motility, biofilm formation and virulence

Bacterial pathogens employ a plethora of virulence factors for host invasion, and their use is tightly regulated to maximize infection efficiency and manage resources in a nutrient-limited environment. Here we show that during Escherichia coli stationary phase the 3' UTR-derived small non-coding RNA FimR2 regulates fimbrial and flagellar biosynthesis at the post-transcriptional level, leading to biofilm formation as the dominant mode of survival under conditions of nutrient depletion. FimR2 interacts with the translational regulator CsrA, antagonizing its functions and firmly tightening control over motility and biofilm formation. Generated through RNase E cleavage, FimR2 regulates stationary phase biology by fine-tuning target mRNA levels independently of the chaperones Hfq and ProQ. The Salmonella enterica orthologue of FimR2 induces effector protein secretion by the type III secretion system and stimulates infection, thus linking the sRNA to virulence. This work reveals the importance of bacterial sRNAs in modulating various aspects of bacterial physiology including stationary phase and virulence.

A 3' UTR-derived small RNA connecting nitrogen and carbon metabolism in enteric bacteria

Increasing numbers of small, regulatory RNAs (sRNAs) corresponding to 3' untranslated regions (UTR) are being discovered in bacteria. One such sRNA, denoted GlnZ, corresponds to the 3' UTR of the Escherichia coli glnA mRNA encoding glutamine synthetase. Several forms of GlnZ, processed from the glnA mRNA, are detected in cells growing with limiting ammonium. GlnZ levels are regulated transcriptionally by the NtrC transcription factor and post-transcriptionally by RNase III. Consistent with the expression, E. coli cells lacking glnZ show delayed outgrowth from nitrogen starvation compared to wild type cells. Transcriptome-wide RNA-RNA interactome datasets indicated that GlnZ binds to multiple target RNAs. Immunoblots and assays of fusions confirmed GlnZ-mediated repression of glnP and sucA, encoding proteins that contribute to glutamine transport and the citric acid cycle, respectively. Although the overall sequences of GlnZ from E. coli K-12, Enterohemorrhagic E. coli and Salmonella enterica have significant differences due to various sequence insertions, all forms of the sRNA were able to regulate the two targets characterized. Together our data show that GlnZ impacts growth of E. coli under low nitrogen conditions by modulating genes that affect carbon and nitrogen flux.

Genome-Wide Identification of Novel sRNAs in Streptococcus mutans

Streptococcus mutans is a major pathobiont involved in the development of dental caries. Its ability to utilize numerous sugars and to effectively respond to environmental stress promotes S. mutans proliferation in oral biofilms. Because of their quick action and low energetic cost, noncoding small RNAs (sRNAs) represent an ideal mode of gene regulation in stress response networks, yet their roles in oral pathogens have remained largely unexplored. We identified 15 novel sRNAs in S. mutans and show that they respond to four stress-inducing conditions commonly encountered by the pathogen in human mouth: sugar-phosphate stress, hydrogen peroxide exposure, high temperature, and low pH. To better understand the role of sRNAs in S. mutans, we further explored the function of the novel sRNA SmsR4. Our data demonstrate that SmsR4 regulates the enzyme IIA (EIIA) component of the sorbitol phosphotransferase system, which transports and phosphorylates the sugar alcohol sorbitol. The fine-tuning of EIIA availability by SmsR4 likely promotes S. mutans growth while using sorbitol as the main carbon source. Our work lays a foundation for understanding the role of sRNAs in regulating gene expression in stress response networks in S. mutans and highlights the importance of the underexplored phenomenon of posttranscriptional gene regulation in oral bacteria. IMPORTANCE Small RNAs (sRNAs) are important gene regulators in bacteria, but the identities and functions of sRNAs in Streptococcus mutans, the principal bacterium involved in the formation of dental caries, are unknown. In this study, we identified 15 putative sRNAs in S. mutans and show that they respond to four common stress-inducing conditions present in human mouth: sugar-phosphate stress, hydrogen peroxide exposure, high temperature, and low pH. We further show that the novel sRNA SmsR4 likely modulates sorbitol transport into the cell by regulating SMU_313 mRNA, which encodes the EIIA subunit of the sorbitol phosphotransferase system. Gaining a better understanding of sRNA-based gene regulation may provide new opportunities to develop specific inhibitors of S. mutans growth, thereby improving oral health.

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.

A Peroxide-Responding sRNA Evolved from a Peroxidase mRNA

Small RNAs (sRNAs) are important gene regulators in bacteria, but it is unclear how new sRNAs originate and become part of regulatory networks that coordinate bacterial response to environmental stimuli. Using a covariance modeling-based approach, we analyzed the presence of hundreds of sRNAs in more than a thousand genomes across Enterobacterales, a bacterial order with a confluence of factors that allows robust genome-scale sRNA analyses: several well-studied organisms with fairly conserved genome structures, an established phylogeny, and substantial nucleotide diversity within a narrow evolutionary space. We discovered that a majority of sRNAs arose recently, and uncovered protein-coding genes as a potential source from which new sRNAs arise. A detailed investigation of the emergence of OxyS, a peroxide-responding sRNA, revealed that it evolved from a fragment of a peroxidase messenger RNA. Importantly, although it replaced the ancestral peroxidase, OxyS continues to be part of the ancestral peroxide-response regulon, indicating that an sRNA that arises from a protein-coding gene would inherently be part of the parental protein's regulatory network. This new insight provides a fresh framework for understanding sRNA origin and regulatory integration in bacteria.

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.

RNase III-mediated processing of a trans-acting bacterial sRNA and its cis-encoded antagonist

Bacterial small RNAs (sRNAs) are important post-transcriptional regulators in stress responses and virulence. They can be derived from an expanding list of genomic contexts, such as processing from parental transcripts by RNase E. The role of RNase III in sRNA biogenesis is less well understood despite its well-known roles in rRNA processing, RNA decay, and cleavage of sRNA-mRNA duplexes. Here, we show that RNase III processes a pair of cis-encoded sRNAs (CJnc190 and CJnc180) of the food-borne pathogen Campylobacter jejuni. While CJnc180 processing by RNase III requires CJnc190, RNase III processes CJnc190 independent of CJnc180 via cleavage of an intramolecular duplex. We also show that CJnc190 directly represses translation of the colonization factor PtmG by targeting a G-rich ribosome-binding site, and uncover that CJnc180 is a cis-acting antagonist of CJnc190, indirectly affecting ptmG regulation. Our study highlights a role for RNase III in sRNA biogenesis and adds cis-encoded RNAs to the expanding diversity of transcripts that can antagonize bacterial sRNAs.

Small RNAs that target G-rich sequences are generated by diverse biogenesis pathways in Epsilonproteobacteria

Bacterial small RNAs (sRNAs) are widespread post-transcriptional regulators controlling bacterial stress responses and virulence. Nevertheless, little is known about how they arise and evolve. Homologues can be difficult to identify beyond the strain level using sequence-based approaches, and similar functionalities can arise by convergent evolution. Here, we found that the virulence-associated CJnc190 sRNA of the foodborne pathogen Campylobacter jejuni resembles the RepG sRNA from the gastric pathogen Helicobacter pylori. However, while both sRNAs bind G-rich sites in their target mRNAs using a C/U-rich loop, they largely differ in their biogenesis. RepG is transcribed from a stand-alone gene and does not require processing, whereas CJnc190 is transcribed from two promoters as precursors that are processed by RNase III and also has a cis-encoded antagonist, CJnc180. By comparing CJnc190 homologues in diverse Campylobacter species, we show that RNase III-dependent processing of CJnc190 appears to be a conserved feature even outside of C. jejuni. We also demonstrate the CJnc180 antisense partner is expressed in C. coli, yet here might be derived from the 3'UTR of the upstream flagella-related gene. Our analysis of G-tract targeting sRNAs in Epsilonproteobacteria demonstrates that similar sRNAs can have markedly different biogenesis pathways.

The 3'UTR-derived sRNA RsaG coordinates redox homeostasis and metabolism adaptation in response to glucose-6-phosphate uptake in Staphylococcus aureus

Staphylococcus aureus RsaG is a 3'-untranslated region (3'UTR) derived sRNA from the conserved uhpT gene encoding a glucose-6-phosphate (G6P) transporter expressed in response to extracellular G6P. The transcript uhpT-RsaG undergoes degradation from 5'- to 3'-end by the action of the exoribonucleases J1/J2, which are blocked by a stable hairpin structure at the 5'-end of RsaG, leading to its accumulation. RsaG together with uhpT is induced when bacteria are internalized into host cells or in the presence of mucus-secreting cells. Using MS2-affinity purification coupled with RNA sequencing, several RNAs were identified as targets including mRNAs encoding the transcriptional factors Rex, CcpA, SarA, and the sRNA RsaI. Our data suggested that RsaG contributes to the control of redox homeostasis and adjusts metabolism to changing environmental conditions. RsaG uses different molecular mechanisms to stabilize, degrade, or repress the translation of its mRNA targets. Although RsaG is conserved only in closely related species, the uhpT 3'UTR of the ape pathogen S. simiae harbors an sRNA, whose sequence is highly different, and which does not respond to G6P levels. Our results hypothesized that the 3'UTRs from UhpT transporter encoding mRNAs could have rapidly evolved to enable adaptation to host niches.

Diversity and Versatility in Small RNA-Mediated Regulation in Bacterial Pathogens

Bacterial gene expression is under the control of a large set of molecules acting at multiple levels. In addition to the transcription factors (TFs) already known to be involved in global regulation of gene expression, small regulatory RNAs (sRNAs) are emerging as major players in gene regulatory networks, where they allow environmental adaptation and fitness. Developments in high-throughput screening have enabled their detection in the entire bacterial kingdom. These sRNAs influence a plethora of biological processes, including but not limited to outer membrane synthesis, metabolism, TF regulation, transcription termination, virulence, and antibiotic resistance and persistence. Almost always noncoding, they regulate target genes at the post-transcriptional level, usually through base-pair interactions with mRNAs, alone or with the help of dedicated chaperones. There is growing evidence that sRNA-mediated mechanisms of actions are far more diverse than initially thought, and that they go beyond the so-called cis- and trans-encoded classifications. These molecules can be derived and processed from 5' untranslated regions (UTRs), coding or non-coding sequences, and even from 3' UTRs. They usually act within the bacterial cytoplasm, but recent studies showed sRNAs in extracellular vesicles, where they influence host cell interactions. In this review, we highlight the various functions of sRNAs in bacterial pathogens, and focus on the increasing examples of widely diverse regulatory mechanisms that might compel us to reconsider what constitute the sRNA.

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.

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.

Prediction of Novel Bacterial Small RNAs From RIL-Seq RNA-RNA Interaction Data

The genomic revolution and subsequent advances in large-scale genomic and transcriptomic technologies highlighted hidden genomic treasures. Among them stand out non-coding small RNAs (sRNAs), shown to play important roles in post-transcriptional regulation of gene expression in both pro- and eukaryotes. Bacterial sRNA-encoding genes were initially identified in intergenic regions, but recent evidence suggest that they can be encoded within other, well-defined, genomic elements. This notion was strongly supported by data generated by RIL-seq, a RNA-seq-based methodology we recently developed for deciphering chaperon-dependent sRNA-target networks in bacteria. Applying RIL-seq to Hfq-bound RNAs in Escherichia coli, we found that ∼64% of the detected RNA pairs involved known sRNAs, suggesting that yet unknown sRNAs may be included in the ∼36% remaining pairs. To determine the latter, we first tested and refined a set of quantitative features derived from RIL-seq data, which distinguish between Hfq-dependent sRNAs and “other RNAs”. We then incorporated these features in a machine learning-based algorithm that predicts novel sRNAs from RIL-seq data, and identified high-scoring candidates encoded in various genomic regions, mostly intergenic regions and 3′ untranslated regions, but also 5′ untranslated regions and coding sequences. Several candidates were further tested and verified by northern blot analysis as Hfq-dependent sRNAs. Our study reinforces the emerging concept that sRNAs are encoded within various genomic elements, and provides a computational framework for the detection of additional sRNAs in Hfq RIL-seq data of E. coli grown under different conditions and of other bacteria manifesting Hfq-mediated sRNA-target interactions.

Three Ribosomal Operons of Escherichia coli Contain Genes Encoding Small RNAs That Interact With Hfq and CsrA in vitro

Three out of the seven ribosomal RNA operons in Escherichia coli end in dual terminator structures. Between the two terminators of each operon is a short sequence that we report here to be an sRNA gene, transcribed as part of the ribosomal RNA primary transcript by read-through of the first terminator. The sRNA genes (rrA, rrB and rrF) from the three operons (rrnA, rrnB and rrnD) are more than 98% identical, and pull-down experiments show that their transcripts interact with Hfq and CsrA. Deletion of rrA, B, F, as well as overexpression of rrB, only modestly affect known CsrA-regulated phenotypes like biofilm formation, pgaA translation and glgC translation, and the role of the sRNAs in vivo may not yet be fully understood. Since RrA, B, F are short-lived and transcribed along with the ribosomal RNA components, their concentration reflect growth-rate regulation at the ribosomal RNA promoters and they could function to fine-tune other growth-phase-dependent processes in the cell. The primary and secondary structure of these small RNAs are conserved among species belonging to different genera of Enterobacteriales.

Inverse regulation of light harvesting and photoprotection is mediated by a 3'-end-derived sRNA in cyanobacteria

Phycobilisomes (PBSs), the principal cyanobacterial antenna, are among the most efficient macromolecular structures in nature, and are used for both light harvesting and directed energy transfer to the photosynthetic reaction center. However, under unfavorable conditions, excess excitation energy needs to be rapidly dissipated to avoid photodamage. The orange carotenoid protein (OCP) senses light intensity and induces thermal energy dissipation under stress conditions. Hence, its expression must be tightly controlled; however, the molecular mechanism of this regulation remains to be elucidated. Here, we describe the discovery of a posttranscriptional regulatory mechanism in Synechocystis sp. PCC 6803 in which the expression of the operon encoding the allophycocyanin subunits of the PBS is directly and in an inverse fashion linked to the expression of OCP. This regulation is mediated by ApcZ, a small regulatory RNA that is derived from the 3'-end of the tetracistronic apcABC-apcZ operon. ApcZ inhibits ocp translation under stress-free conditions. Under most stress conditions, apc operon transcription decreases and ocp translation increases. Thus, a key operon involved in the collection of light energy is functionally connected to the expression of a protein involved in energy dissipation. Our findings support the view that regulatory RNA networks in bacteria evolve through the functionalization of mRNA 3'-UTRs.

An mRNA-mRNA Interaction Couples Expression of a Virulence Factor and Its Chaperone in Listeria monocytogenes

Bacterial pathogens often employ RNA regulatory elements located in the 5' untranslated regions (UTRs) to control gene expression. Using a comparative structural analysis, we examine the structure of 5' UTRs at a global scale in the pathogenic bacterium Listeria monocytogenes under different conditions. In addition to discovering an RNA thermoswitch and detecting simultaneous interaction of ribosomes and small RNAs with mRNA, we identify structural changes in the 5' UTR of an mRNA encoding the post-translocation chaperone PrsA2 during infection conditions. We demonstrate that the 5' UTR of the prsA2 mRNA base pairs with the 3' UTR of the full-length hly mRNA encoding listeriolysin O, thus preventing RNase J1-mediated degradation of the prsA2 transcript. Mutants lacking the hly-prsA2 interaction exhibit reduced virulence properties. This work highlights an additional level of RNA regulation, where the mRNA encoding a chaperone is stabilized by the mRNA encoding its substrate.

Critical Role of 3'-Downstream Region of pmrB in Polymyxin Resistance in Escherichia coli BL21(DE3)

Polymyxins, such as colistin and polymyxin B, are the drugs used as a last resort to treat multidrug-resistant Gram-negative bacterial infections in humans. Increasing colistin resistance has posed a serious threat to human health, warranting in-depth mechanistic research. In this study, using a functional cloning approach, we examined the molecular basis of colistin resistance in Escherichia coli BL21(DE3). Five transformants with inserts ranging from 3.8 to 10.7 kb displayed significantly increased colistin resistance, three of which containing pmrB locus and two containing pmrD locus. Stepwise subcloning indicated that both the pmrB with a single G361A mutation and at least a 103 bp downstream region of pmrB are essential for conferring colistin resistance. Analysis of the mRNA level and stability showed that the length of the downstream region drastically affected the pmrB mRNA level but not its half-life. Lipid A analysis, by mass spectrometry, revealed that the constructs containing pmrB with a longer downstream region (103 or 126 bp) have charge-altering l-4-aminoarabinose (Ara4N) and phosphoethanolamine (pEtN) modifications in lipid A, which were not observed in both vector control and the construct containing pmrB with an 86 bp downstream region. Together, the findings from this study indicate that the 3′-downstream region of pmrB is critical for the PmrB-mediated lipid A modifications and colistin resistance in E. coli BL21(DE3), suggesting a novel regulatory mechanism of PmrB-mediated colistin resistance in E. coli.

Toward a Comprehensive Analysis of Posttranscriptional Regulatory Networks: a New Tool for the Identification of Small RNA Regulators of Specific mRNAs

A number of computational or experimental tools have been developed to identify targets of small RNA (sRNA) regulation. Here, we modified one of these methods, based on in vivo proximity ligation of sRNAs bound to their targets, referred to as rGRIL-seq, that can be used to capture sRNA regulators of a gene of interest. Intracellular expression of bacteriophage T4 RNA ligase leads to a covalent linking of sRNAs base-paired with mRNAs, and the chimeras are captured using oligonucleotides complementary to the mRNA, followed by sequencing. This allows the identification of known as well as novel sRNAs. We applied rGRIL-seq toward finding sRNA regulators of expression of the stress response sigma factor RpoS in Escherichia coli, Pseudomonas aeruginosa, and Vibrio cholerae. In E. coli, we confirmed the regulatory role of known sRNAs and discovered a new negative regulator, asYbiE. When applied to P. aeruginosa and V. cholerae, we identified two novel sRNAs (s03661 and s0223) in P. aeruginosa and two known sRNAs (TfoR and Vcr043) in V. cholerae as direct regulators of rpoS. The use of rGRIL-seq for defining multiple posttranscriptional regulatory inputs into individual mRNAs represents a step toward a more comprehensive understanding of the workings of bacterial regulatory networks.

A Grad-seq View of RNA and Protein Complexes in Pseudomonas aeruginosa under Standard and Bacteriophage Predation Conditions

The Gram-negative rod-shaped bacterium Pseudomonas aeruginosa is not only a major cause of nosocomial infections but also serves as a model species of bacterial RNA biology. While its transcriptome architecture and posttranscriptional regulation through the RNA-binding proteins Hfq, RsmA, and RsmN have been studied in detail, global information about stable RNA-protein complexes in this human pathogen is currently lacking. Here, we implement gradient profiling by sequencing (Grad-seq) in exponentially growing P. aeruginosa cells to comprehensively predict RNA and protein complexes, based on glycerol gradient sedimentation profiles of >73% of all transcripts and ∼40% of all proteins. As to benchmarking, our global profiles readily reported complexes of stable RNAs of P. aeruginosa, including 6S RNA with RNA polymerase and associated product RNAs (pRNAs). We observe specific clusters of noncoding RNAs, which correlate with Hfq and RsmA/N, and provide a first hint that P. aeruginosa expresses a ProQ-like FinO domain-containing RNA-binding protein. To understand how biological stress may perturb cellular RNA/protein complexes, we performed Grad-seq after infection by the bacteriophage ΦKZ. This model phage, which has a well-defined transcription profile during host takeover, displayed efficient translational utilization of phage mRNAs and tRNAs, as evident from their increased cosedimentation with ribosomal subunits. Additionally, Grad-seq experimentally determines previously overlooked phage-encoded noncoding RNAs. Taken together, the Pseudomonas protein and RNA complex data provided here will pave the way to a better understanding of RNA-protein interactions during viral predation of the bacterial cell.

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.

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.

An RNA-centric view on gut Bacteroidetes

Bacteria employ noncoding RNAs to maintain cellular physiology, adapt global gene expression to fluctuating environments, sense nutrients, coordinate their interaction with companion microbes and host cells, and protect themselves against bacteriophages. While bacterial RNA research has made fundamental contributions to biomedicine and biotechnology, the bulk of our knowledge of RNA biology stems from the study of a handful of aerobic model species. In comparison, RNA research is lagging in many medically relevant obligate anaerobic species, in particular the numerous commensal bacteria comprising our gut microbiota. This review presents a guide to RNA-based regulatory mechanisms in the phylum Bacteroidetes, focusing on the most abundant bacterial genus in the human gut, Bacteroides spp. This includes recent case reports on riboswitches, an mRNA leader, cis- and trans-encoded small RNAs (sRNAs) in Bacteroides spp., and a survey of CRISPR-Cas systems across Bacteroidetes. Recent work from our laboratory now suggests the existence of hundreds of noncoding RNA candidates in Bacteroides thetaiotaomicron, the emerging model organism for functional microbiota research. Based on these collective observations, we predict mechanistic and functional commonalities and differences between Bacteroides sRNAs and those of other model bacteria, and outline open questions and tools needed to boost Bacteroidetes RNA research.

Adaptor protein RapZ activates endoribonuclease RNase E by protein-protein interaction to cleave a small regulatory RNA

In Escherichia coli, endoribonuclease RNase E initiates degradation of many RNAs and represents a hub for post-transcriptional regulation. The tetrameric adaptor protein RapZ targets the small regulatory RNA GlmZ to degradation by RNase E. RapZ binds GlmZ through a domain located at the carboxyl terminus and interacts with RNase E, promoting GlmZ cleavage in the base-pairing region. When necessary, cleavage of GlmZ is counteracted by the homologous small RNA GlmY, which sequesters RapZ through molecular mimicry. In the current study, we addressed the molecular mechanism employed by RapZ. We show that RapZ mutants impaired in RNA-binding but proficient in binding RNase E are able to stimulate GlmZ cleavage in vivo and in vitro when provided at increased concentrations. In contrast, a truncated RapZ variant retaining RNA-binding activity but incapable of contacting RNase E lacks this activity. In agreement, we find that tetrameric RapZ binds the likewise tetrameric RNase E through direct interaction with its large globular domain within the catalytic amino terminus, independent of RNA. Although RapZ stimulates cleavage of at least one non-cognate RNA by RNase E in vitro, its activity is restricted to GlmZ in vivo as revealed by RNA sequencing, suggesting that certain features within the RNA substrate are also required for cleavage. In conclusion, RapZ boosts RNase E activity through interaction with its catalytic domain, which represents a novel mechanism of RNase E activation. In contrast, RNA-binding has a recruiting role, increasing the likelihood that productive RapZ/GlmZ/RNase E complexes form.

Gene autoregulation by 3' UTR-derived bacterial small RNAs

Negative feedback regulation, that is the ability of a gene to repress its own synthesis, is the most abundant regulatory motif known to biology. Frequently reported for transcriptional regulators, negative feedback control relies on binding of a transcription factor to its own promoter. Here, we report a novel mechanism for gene autoregulation in bacteria relying on small regulatory RNA (sRNA) and the major endoribonuclease, RNase E. TIER-seq analysis (transiently-inactivating-an-endoribonuclease-followed-by-RNA-seq) revealed ~25,000 RNase E-dependent cleavage sites in Vibrio cholerae, several of which resulted in the accumulation of stable sRNAs. Focusing on two examples, OppZ and CarZ, we discovered that these sRNAs are processed from the 3' untranslated region (3' UTR) of the oppABCDF and carAB operons, respectively, and base-pair with their own transcripts to inhibit translation. For OppZ, this process also triggers Rho-dependent transcription termination. Our data show that sRNAs from 3' UTRs serve as autoregulatory elements allowing negative feedback control at the post-transcriptional level.

A high-resolution transcriptome map identifies small RNA regulation of metabolism in the gut microbe Bacteroides thetaiotaomicron

Bacteria of the genus Bacteroides are common members of the human intestinal microbiota and important degraders of polysaccharides in the gut. Among them, the species Bacteroides thetaiotaomicron has emerged as the model organism for functional microbiota research. Here, we use differential RNA sequencing (dRNA-seq) to generate a single-nucleotide resolution transcriptome map of B. thetaiotaomicron grown under defined laboratory conditions. An online browser, called ‘Theta-Base’ (www.helmholtz-hiri.de/en/datasets/bacteroides), is launched to interrogate the obtained gene expression data and annotations of ~4500 transcription start sites, untranslated regions, operon structures, and 269 noncoding RNA elements. Among the latter is GibS, a conserved, 145 nt-long small RNA that is highly expressed in the presence of N-acetyl-D-glucosamine as sole carbon source. We use computational predictions and experimental data to determine the secondary structure of GibS and identify its target genes. Our results indicate that sensing of N-acetyl-D-glucosamine induces GibS expression, which in turn modifies the transcript levels of metabolic enzymes.

Bacteroides thetaiotaomicron is a human gut microbe and an emergent model organism. Here, Ryan et al. generate single-nucleotide resolution RNA-seq data for this bacterium and map transcription start sites and noncoding RNAs, one of which modulates expression of metabolic enzymes.

Differential evolution in 3'UTRs leads to specific gene expression in Staphylococcus

The evolution of gene expression regulation has contributed to species differentiation. The 3′ untranslated regions (3′UTRs) of mRNAs include regulatory elements that modulate gene expression; however, our knowledge of their implications in the divergence of bacterial species is currently limited. In this study, we performed genome-wide comparative analyses of mRNAs encoding orthologous proteins from the genus Staphylococcus and found that mRNA conservation was lost mostly downstream of the coding sequence (CDS), indicating the presence of high sequence diversity in the 3′UTRs of orthologous genes. Transcriptomic mapping of different staphylococcal species confirmed that 3′UTRs were also variable in length. We constructed chimeric mRNAs carrying the 3′UTR of orthologous genes and demonstrated that 3′UTR sequence variations affect protein production. This suggested that species-specific functional 3′UTRs might be specifically selected during evolution. 3′UTR variations may occur through different processes, including gene rearrangements, local nucleotide changes, and the transposition of insertion sequences. By extending the conservation analyses to specific 3′UTRs, as well as the entire set of Escherichia coli and Bacillus subtilis mRNAs, we showed that 3′UTR variability is widespread in bacteria. In summary, our work unveils an evolutionary bias within 3′UTRs that results in species-specific non-coding sequences that may contribute to bacterial diversity.

Hfq CLASH uncovers sRNA-target interaction networks linked to nutrient availability adaptation

By shaping gene expression profiles, small RNAs (sRNAs) enable bacteria to efficiently adapt to changes in their environment. To better understand how Escherichia coli acclimatizes to nutrient availability, we performed UV cross-linking, ligation and sequencing of hybrids (CLASH) to uncover Hfq-associated RNA-RNA interactions at specific growth stages. We demonstrate that Hfq CLASH robustly captures bona fide RNA-RNA interactions. We identified hundreds of novel sRNA base-pairing interactions, including many sRNA-sRNA interactions and involving 3’UTR-derived sRNAs. We rediscovered known and identified novel sRNA seed sequences. The sRNA-mRNA interactions identified by CLASH have strong base-pairing potential and are highly enriched for complementary sequence motifs, even those supported by only a few reads. Yet, steady state levels of most mRNA targets were not significantly affected upon over-expression of the sRNA regulator. Our results reinforce the idea that the reproducibility of the interaction, not base-pairing potential, is a stronger predictor for a regulatory outcome.

New sequencing methodologies reveal interplay between multiple RNA-binding proteins and their RNAs

It is now established that base-pairing regulatory RNAs are key players in post-transcriptional regulatory networks where they affect the translation and/or stability of their target RNAs. In many cases, the base-pairing between two RNAs is facilitated by an RNA-binding protein (RBP) that serves as an RNA chaperone. Recent advances in sequencing methods have revealed the RNA populations bound by the RBPs, yielding insights valuable into regulatory networks. Further analyses of these networks can improve our understanding of the roles played by RBPs in the regulation of gene expression by regulatory RNAs, especially when multiple RBPs are involved. For example, using an RNA sequencing-based methodology that captures RNA-RNA pairs on RBP, an interplay between two RBPs in bacteria that compete on the same RNA-RNA pair was revealed. In this case, one protein promotes negative regulation of the target RNA, while the second protein can block this regulation. In this mini-review, I outline the exciting future directions that can be taken to deepen our understanding of the roles played by RBPs in post-transcriptional regulation, and discuss how the different sequencing methods can assist in deciphering the relationships among RBPs, and between the RBPs and the RNAs they bind. Having a more detailed picture of the RBPs-RNAs network will elucidate how bacteria can have nuanced control of gene expression, critical for survival in the varied environments in which bacteria live.

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.

The conserved 3' UTR-derived small RNA NarS mediates mRNA crossregulation during nitrate respiration

Small noncoding RNAs (sRNAs) from mRNA 3′ UTRs seem to present a previously unrecognized layer of bacterial post-transcriptional control whereby mRNAs influence each other's expression, independently of transcriptional control. Studies in Escherichia coli and Salmonella enterica showed that such sRNAs are natural products of RNase E-mediated mRNA decay and associate with major RNA-binding proteins (RBPs) such as Hfq and ProQ. If so, there must be additional sRNAs from mRNAs that accumulate only under specific physiological conditions. We test this prediction by characterizing candidate NarS that represents the 3′ UTR of nitrate transporter NarK whose gene is silent during standard aerobic growth. We find that NarS acts by Hfq-dependent base pairing to repress the synthesis of the nitrite transporter, NirC, resulting in mRNA cross-regulation of nitrate and nitrite transporter genes. Interestingly, the NarS-mediated repression selectively targets the nirC cistron of the long nirBDC-cysG operon, an observation that we rationalize as a mechanism to protect the bacterial cytoplasm from excessive nitrite toxicity during anaerobic respiration with abundant nitrate. Our successful functional assignment of a 3′ UTR sRNA from a non-standard growth condition supports the notion that mRNA crossregulation is more pervasive than currently appreciated.

Trans-Acting Small RNAs and Their Effects on Gene Expression in Escherichia coli and Salmonella enterica

The last few decades have led to an explosion in our understanding of the major roles that small regulatory RNAs (sRNAs) play in regulatory circuits and the responses to stress in many bacterial species. Much of the foundational work was carried out with Escherichia coli and Salmonella enterica serovar Typhimurium. The studies of these organisms provided an overview of how the sRNAs function and their impact on bacterial physiology, serving as a blueprint for sRNA biology in many other prokaryotes. They also led to the development of new technologies. In this chapter, we first summarize how these sRNAs were identified, defining them in the process. We discuss how they are regulated and how they act and provide selected examples of their roles in regulatory circuits and the consequences of this regulation. Throughout, we summarize the methodologies that were developed to identify and study the regulatory RNAs, most of which are applicable to other bacteria. Newly updated databases of the known sRNAs in E. coli K-12 and S. enterica Typhimurium SL1344 serve as a reference point for much of the discussion and, hopefully, as a resource for readers and for future experiments to address open questions raised in this review.

Differential quantitative proteomics reveals the functional difference of two yigP locus products, UbiJ and EsrE

The yigP (ubiJ) locus has been shown to be associated with many phenotypic changes in Escherichia coli, while the individual function of its two products, EsrE small RNA and UbiJ protein, is still elusive. In this study, we constructed two single-element mutants, EsrE mutant strain Mut and UbiJ mutant strain Ter, on the basis of the base substitution programs. The variable antibiotics resistance and ubiquinone (UQ, coenzyme Q) yield and the similar cell growth between mutants revealed the division of labor and collaboration of EsrE and UbiJ in JM83. Furthermore, we detected the concentration of intracellular proteins of Mut and Ter by stable isotope-labeled quantitative proteomics. The results demonstrate that both EsrE and UbiJ are involved in the aerobic growth of E. coli, while EsrE preferentially contributes to the amino acid-related pathway, and UbiJ is an indispensable factor in the biosynthesis of UQ. Moreover, we uncovered a potential regulatory circuit of d-cycloserine (DCS) that composed of EsrE, GcvA, and GcvB by proteomic analysis.

The Transcription Unit Architecture of Streptomyces lividans TK24

Streptomyces lividans is an attractive host for production of heterologous proteins and secondary metabolites of other Streptomyces species. To fully harness the industrial potential of S. lividans, understanding its metabolism and genetic regulatory elements is essential. This study aimed to determine its transcription unit (TU) architecture and elucidate its diverse regulatory elements, including promoters, ribosome binding sites, 5′-untranslated regions, and transcription terminators. Total 1,978 transcription start sites and 1,640 transcript 3′-end positions were identified, which were integrated to determine 1,300 TUs, consistent with transcriptomic profiles. The conserved promoter sequences were found as 5′-TANNNT and 5′-TGAC, representing the −10 and −35 elements, respectively. Analysis of transcript 3′-end positions revealed the presence of distinctive terminator sequences and the RNA stem structure responsible for the determination of the 3′-boundary of a transcript. Functionally related genes are likely to be regulated simultaneously by using similar promoters and being transcribed as a poly-cistronic TU. Poly-cistronic TUs were further processed or alternatively transcribed into multiple TUs to fine-regulate individual genes in response to environmental conditions. The TU information and regulatory elements identified will serve as invaluable resources for understanding the complex regulatory mechanisms of S. lividans and to elevate its industrial potential.

Broadening the Definition of Bacterial Small RNAs: Characteristics and Mechanisms of Action

The first report of trans-acting RNA-based regulation in bacterial cells dates back to 1984. Subsequent studies in diverse bacteria unraveled shared properties of trans-acting small regulatory RNAs, forming a clear definition of these molecules. These shared characteristics have been used extensively to identify new small RNAs (sRNAs) and their interactomes. Recently however, emerging technologies able to resolve RNA-RNA interactions have identified new types of regulatory RNAs. In this review, we present a broader definition of trans-acting sRNA regulators and discuss their newly discovered intrinsic characteristics.

Functional expansion of a TCA cycle operon mRNA by a 3' end-derived small RNA

Global RNA profiling studies in bacteria have predicted the existence of many of small noncoding RNAs (sRNAs) that are processed off mRNA 3′ ends to regulate other mRNAs via the RNA chaperones Hfq and ProQ. Here, we present targets of SdhX (RybD), an Hfq-dependent sRNA that is generated by RNase E mediated 3′ processing of the ∼10 000-nt mRNA of the TCA cycle operon sdhCDAB-sucABCD in enteric bacteria. An in silico search predicted ackA mRNA, which encodes acetate kinase, as a conserved primary target of SdhX. Through base pairing, SdhX represses AckA synthesis during growth of Salmonella on acetate. Repression can be achieved by a naturally occurring 38-nucleotide SdhX variant, revealing the shortest functional Hfq-associated sRNA yet. Salmonella SdhX also targets the mRNAs of fumB (anaerobic fumarase) and yfbV, a gene of unknown function adjacent to ackA. Instead, through a slightly different seed sequence, SdhX can repress other targets in Escherichia coli, namely katG (catalase) and fdoG (aerobic formate dehydrogenase). This study illustrates how a key operon from central metabolism is functionally connected to other metabolic pathways through a 3′ appended sRNA, and supports the notion that mRNA 3′UTRs are a playground for the evolution of regulatory RNA networks in bacteria.

Hfq Globally Binds and Destabilizes sRNAs and mRNAs in Yersinia pestis

Hfq is a ubiquitous Sm-like RNA-binding protein in bacteria involved in physiological fitness and pathogenesis, while its in vivo binding nature remains elusive. Here we reported genome-wide Hfq-bound RNAs in Yersinia pestis, a causative agent of plague, by using cross-linking immunoprecipitation coupled with deep sequencing (CLIP-seq) approach. We show that the Hfq binding density is enriched in more than 80% mRNAs of Y. pestis and that Hfq also globally binds noncoding small RNAs (sRNAs) encoded by the intergenic, antisense, and 3' regions of mRNAs. An Hfq U-rich stretch is highly enriched in sRNAs, while motifs partially complementary to AGAAUAA and GGGGAUUA are enriched in both mRNAs and sRNAs. Hfq-binding motifs are enriched at both terminal sites and in the gene body of mRNAs. Surprisingly, a large fraction of the sRNA and mRNA regions bound by Hfq and those downstream are destabilized, likely via a 5'P-activated RNase E degradation pathway, which is consistent with a model in which Hfq facilitates sRNA-mRNA base pairing and the coupled degradation in Y. pestis These results together have presented a high-quality Hfq-RNA interaction map in Y. pestis, which should be important for further deciphering the regulatory role of Hfq-sRNAs in Y. pestis IMPORTANCE Discovered in 1968 as an Escherichia coli host factor that was essential for replication of the bacteriophage Qβ, the Hfq protein is a ubiquitous and highly abundant RNA-binding protein in many bacteria. With the assistance of Hfq, small RNAs in bacteria play important roles in regulating the stability and translation of mRNAs by base pairing. In this study, we want to elucidate the Hfq-assisted sRNA-mRNA regulation in Yersinia pestis A global map of Hfq interaction sites in Y. pestis was obtained by sequencing cDNAs converted from the Hfq-bound RNA fragments using UV cross-linking coupled immunoprecipitation technology. We demonstrate that Hfq could bind to hundreds of sRNAs and the majority of mRNAs in Y. pestis The enriched binding motifs in sRNAs and mRNAs are complementary to each other, suggesting a general base-pairing mechanism for sRNA-mRNA interaction. The Hfq-bound sRNA and mRNA regions were both destabilized. The results suggest that Hfq binding facilitates sRNA-mRNA base pairing and coordinates their degradation, which might enable Hfq to surveil the homeostasis of most mRNAs in bacteria.

Coordinate regulation of the expression of SdsR toxin and its downstream pphA gene by RyeA antitoxin in Escherichia coli

In Escherichia coli, SdsR and RyeA, a unique pair of mutually cis-encoded small RNAs (sRNAs), act as toxin and antitoxin, respectively. SdsR and RyeA expression are reciprocally regulated; however, how each regulates the synthesis of the other remains unclear. Here, we characterized the biosynthesis of the two sRNAs during growth and investigated their coordinate regulation using sdsR and ryeA promoter mutant strains. We found that RyeA transcription occurred even upon entry of cells into the stationary phase, but its apparent expression was restricted to exponentially growing cells because of its degradation by SdsR. Likewise, the appearance of SdsR was delayed owing to its RyeA-mediated degradation. We also found that the sdsR promoter was primarily responsible for transcription of the downstream pphA gene encoding a phosphatase and that pphA mRNA was synthesized by transcriptional read-through over the sdsR terminator. Transcription from the σ⁷⁰-dependent ryeA promoter inhibited transcription from the σ^S-dependent sdsR promoter through transcriptional interference. This transcriptional inhibition also downregulated pphA expression, but RyeA itself did not downregulate pphA expression.