Detection of 5'- and 3'-UTR-derived small RNAs and cis-encoded antisense RNAs in Escherichia coli

Small RNA Esr41 inversely regulates expression of LEE and flagellar genes in enterohaemorrhagic Escherichia coli

Enterohaemorrhagic Escherichia coli (EHEC) is a life-threatening human pathogen worldwide. The locus of enterocyte effacement (LEE) in EHEC encodes a type three secretion system and effector proteins, all of which are essential for bacterial adherence to host cells. When LEE expression is activated, flagellar gene expression is down-regulated because bacterial flagella induce the immune responses of host cells at the infection stage. Therefore, this inverse regulation is also important for EHEC infection. We report here that a small regulatory RNA (sRNA), Esr41, mediates LEE repression and flagellar gene activation. Multiple copies of esr41 abolished LEE expression by down-regulating the expression of ler and pch, which encode positive regulators of LEE. This regulation led to reduced EHEC adhesion to host cells. Translational gene-reporter fusion experiments revealed that Esr41 regulates ler expression at a post-transcriptional level, and pch transcription, probably via an unknown target of Esr41. Esr41-mediated ler and pch repression was not observed in cells lacking hfq, which encodes an RNA-binding protein essential for most sRNA functions, indicating that Esr41 acts in an Hfq-dependent manner. We previously reported an increase in cell motility induced by Esr41. This motility enhancement was also observed in EHEC lacking ler, showing that Esr41-mediated enhancement of cell motility is in a ler-independent manner. In addition, Esr41 activated the expression of flagellar Class 3 genes by indirectly inducing the transcription of fliA, which encodes the sigma factor for flagellar synthesis. These results suggest that Esr41 plays important roles in the inverse regulation of LEE and flagellar gene expression.

Pervasive Regulatory Functions of mRNA Structure Revealed by High-Resolution SHAPE Probing

mRNAs can fold into complex structures that regulate gene expression. Resolving such structures de novo has remained challenging and has limited our understanding of the prevalence and functions of mRNA structure. We use SHAPE-MaP experiments in living E. coli cells to derive quantitative, nucleotide-resolution structure models for 194 endogenous transcripts encompassing approximately 400 genes. Individual mRNAs have exceptionally diverse architectures, and most contain well-defined structures. Active translation destabilizes mRNA structure in cells. Nevertheless, mRNA structure remains similar between in-cell and cell-free environments, indicating broad potential for structure-mediated gene regulation. We find that the translation efficiency of endogenous genes is regulated by unfolding kinetics of structures overlapping the ribosome binding site. We discover conserved structured elements in 35% of UTRs, several of which we validate as novel protein binding motifs. RNA structure regulates every gene studied here in a meaningful way, implying that most functional structures remain to be discovered.

Translation Stress Positively Regulates MscL-Dependent Excretion of Cytoplasmic Proteins

The apparent mislocalization or excretion of cytoplasmic proteins is a commonly observed phenomenon in both bacteria and eukaryotes. However, reports on the mechanistic basis and the cellular function of this so-called “nonclassical protein secretion” are limited. Here we report that protein overexpression in recombinant cells and antibiotic-induced translation stress in wild-type Escherichia coli cells both lead to excretion of cytoplasmic protein (ECP). Condition-specific metabolomic and proteomic analyses, combined with genetic knockouts, indicate a role for both the large mechanosensitive channel (MscL) and the alternative ribosome rescue factor A (ArfA) in ECP. Collectively, the findings indicate that MscL-dependent protein excretion is positively regulated in response to both osmotic stress and arfA-mediated translational stress.

Protein translocation is an essential feature of cellular organisms. Bacteria, like all single-cell organisms, interact with their environment by translocation of proteins across their cell membranes via dedicated secretion pathways. Proteins destined for secretion are directed toward the secretion pathways by the presence of specific signal peptides. This study demonstrates that under conditions of both osmotic stress and translation stress, E. coli cells undergo an excretion phenomenon whereby signal peptide-less proteins are translocated across both the inner and outer cell membranes into the extracellular environment. Confirming the presence of alternative translocation/excretion pathways and understanding their function and regulation are thus important for fundamental microbiology and biotechnology.

The malS-5'UTR weakens the ability of Salmonella enterica serovar Typhi to survive in macrophages by increasing intracellular ATP levels

Bacterial non-coding RNAs (ncRNAs), as important regulatory factors, are involved in many cellular processes, including virulence and protection against environmental stress. The 5' untranslated region (UTR) of malS (named malS-5'UTR), a regulatory ncRNA, increases the invasive capacity and influences histidine biosynthesis in Salmonella enterica serovar Typhi (S. Typhi). In this study, we found that overexpression of the malS-5'UTR decreased S. Typhi survival within macrophages. A microarray analysis of a strain overexpressing the malS-5'UTR revealed a significant increase in the mRNA levels of the atp operon. The intracellular ATP levels were elevated in the malS-5'UTR overexpression strain. Quantitative real-time polymerase chain reaction results showed that the malS-5'UTR downregulated the mRNA levels of phoP, phoQ, and mgtC. MgtC, its expression is regulated by PhoP/PhoQ two-component regulatory system, inhibits the F1F0 ATP synthase, thereby preventing the accumulation of ATP to non-physiological levels and the acidification of the cytoplasm within macrophages. Thus, we propose that the malS-5'UTR weakens the ability of S. Typhi to survive in macrophages, probably because of the accumulation of ATP within macrophages, by regulating the mRNA levels of mgtC and the atp operon in a phoP-dependent manner.

Identification of sRNA mediated responses to nutrient depletion in Burkholderia pseudomallei

The Burkholderia genus includes many species that are known to survive in diverse environmental conditions including low nutrient environments. One species, Burkholderia pseudomallei is a versatile pathogen that can survive in a wide range of hosts and environmental conditions. In this study, we investigated how a nutrient depleted growth environment evokes sRNA mediated responses by B. pseudomallei. Computationally predicted B. pseudomallei D286 sRNAs were mapped to RNA-sequencing data for cultures grown under two conditions: (1) BHIB as a nutrient rich media reference environment and (2) M9 media as a nutrient depleted stress environment. The sRNAs were further selected to identify potentially cis-encoded systems by investigating their possible interactions with their flanking genes. The mappings of predicted sRNA genes and interactions analysis to their flanking genes identified 12 sRNA candidates that may possibly have cis-acting regulatory roles that are associated to a nutrient depleted growth environment. Our approach can be used for identifying novel sRNA genes and their possible role as cis-mediated regulatory systems.

Identification of small RNAs abundant in Burkholderia cenocepacia biofilms reveal putative regulators with a potential role in carbon and iron metabolism

Small RNAs play a regulatory role in many central metabolic processes of bacteria, as well as in developmental processes such as biofilm formation. Small RNAs of Burkholderia cenocepacia, an opportunistic pathogenic beta-proteobacterium, are to date not well characterised. To address that, we performed genome-wide transcriptome structure analysis of biofilm grown B. cenocepacia J2315. 41 unannotated short transcripts were identified in intergenic regions of the B. cenocepacia genome. 15 of these short transcripts, highly abundant in biofilms, widely conserved in Burkholderia sp. and without known function, were selected for in-depth analysis. Expression profiling showed that most of these sRNAs are more abundant in biofilms than in planktonic cultures. Many are also highly abundant in cells grown in minimal media, suggesting they are involved in adaptation to nutrient limitation and growth arrest. Their computationally predicted targets include a high proportion of genes involved in carbon metabolism. Expression and target genes of one sRNA suggest a potential role in regulating iron homoeostasis. The strategy used for this study to detect sRNAs expressed in B. cenocepacia biofilms has successfully identified sRNAs with a regulatory function.

A trans-acting leader RNA from a Salmonella virulence gene

Bacteria use flagella to move toward nutrients, find its host, or retract from toxic substances. Because bacterial flagellum is one of the ligands that activate the host innate immune system, its synthesis should be tightly regulated during host infection, which is largely unknown. Here, we report that a bacterial leader mRNA from the mgtCBR virulence operon in the intracellular pathogen Salmonella enterica serovar Typhimurium binds to the fljB coding region of mRNAs in the fljBA operon encoding the FljB phase 2 flagellin, a main component of bacterial flagella and the FljA repressor for the FliC phase 1 flagellin, and degrades fljBA mRNAs in an RNase E-dependent fashion during infection. A nucleotide substitution of the fljB flagellin gene that prevents the mgtC leader RNA-mediated down-regulation increases the fljB-encoded flagellin synthesis, leading to a hypermotile phenotype inside macrophages. Moreover, the fljB nucleotide substitution renders Salmonella hypervirulent, indicating that FljB-based motility must be compromised in the phagosomal compartment where Salmonella resides. This suggests that this pathogen promotes pathogenicity by producing a virulence protein and limits locomotion by a trans-acting leader RNA from the same virulence gene during infection.

Role of the terminator hairpin in the biogenesis of functional Hfq-binding sRNAs

Rho-independent transcription terminators of the genes encoding bacterial Hfq-binding sRNAs possess a set of seven or more T residues at the 3' end, as noted in previous studies. Here, we have studied the role of the terminator hairpin in the biogenesis of sRNAs focusing on SgrS and RyhB in Escherichia coli. We constructed variant sRNA genes in which the GC-rich inverted repeat sequences are extended to stabilize the terminator hairpins. We demonstrate that the extension of the hairpin stem leads to generation of heterogeneous transcripts in which the poly(U) tail is shortened. The transcripts with shortened poly(U) tails no longer bind to Hfq and lose the ability to repress the target mRNAs. The shortened transcripts are generated in an in vitro transcription system with purified RNA polymerase, indicating that the generation of shortened transcripts is caused by premature transcription termination. We conclude that the terminator structure of sRNA genes is optimized to generate functional sRNAs. Thus, the Rho-independent terminators of sRNA genes possess two common features: a long T residue stretch that is a prerequisite for generation of functional sRNAs and a moderate strength of hairpin structure that ensures the termination at the seventh or longer position within the consecutive T stretch. The modulation of the termination position at the Rho-independent terminators is critical for biosynthesis of functional sRNAs.

The 5΄ UTR of the type I toxin ZorO can both inhibit and enhance translation

Many bacterial type I toxin mRNAs possess a long 5΄ untranslated region (UTR) that serves as the target site of the corresponding antitoxin sRNA. This is the case for the zorO-orzO type I system where the OrzO antitoxin base pairs to the 174-nucleotide zorO 5΄ UTR. Here, we demonstrate that the full-length 5΄ UTR of the zorO type I toxin hinders its own translation independent of the sRNA whereas a processed 5΄ UTR (zorO Δ28) promotes translation. The full-length zorO 5΄ UTR folds into an extensive secondary structure sequestering the ribosome binding site (RBS). Processing of the 5΄ UTR does not alter the RBS structure, but opens a large region (EAP region) located upstream of the RBS. Truncation of this EAP region impairs zorO translation, but this defect can be rescued upon exposing the RBS. Additionally, the region spanning +35 to +50 of the zorO mRNA is needed for optimal translation of zorO. Importantly, the positive and negative effects on translation imparted by the 5΄ UTR can be transferred onto a reporter gene, indicative that the 5΄ UTR can solely drive regulation. Moreover, we show that the OrzO sRNA can inhibit zorO translation via base pairing to the of the EAP region.

RNA-based regulation in type I toxin-antitoxin systems and its implication for bacterial persistence

Bacterial dormancy is a valuable survival strategy upon challenging environmental conditions. Dormant cells tolerate the consequences of high stress levels and may re-populate the environment upon return to favorable conditions. Antibiotic-tolerant bacteria—termed persisters—regularly cause relapsing infections, increase the likelihood of antibiotic resistance, and, therefore, earn increasing attention. Their generation often depends on toxins from chromosomal toxin–antitoxin systems. Here, we review recent insights concerning RNA-based control of toxin synthesis, and discuss possible implications for persister generation.

Why so narrow: Distribution of anti-sense regulated, type I toxin-antitoxin systems compared with type II and type III systems

Toxin-antitoxin (TA) systems are gene modules that appear to be horizontally mobile across a wide range of prokaryotes. It has been proposed that type I TA systems, with an antisense RNA-antitoxin, are less mobile than other TAs that rely on direct toxin-antitoxin binding but no direct comparisons have been made. We searched for type I, II and III toxin families using iterative searches with profile hidden Markov models across phyla and replicons. The distribution of type I toxin families were comparatively narrow, but these patterns weakened with recently discovered families. We discuss how the function and phenotypes of TA systems as well as biases in our search methods may account for differences in their distribution.

Growth Phase-dependent Variation of RNase BN/Z Affects Small RNAs: REGULATION OF 6S RNA

RNase BN, the RNase Z family member in E. coli, can participate in the processing of tRNA precursors. However, this function only becomes apparent when other processing enzymes are absent, raising the question of its primary physiological role. Here, we show that RNase BN itself is subject to growth phase-dependent regulation, because both rbn mRNA and RNase BN protein are at their highest levels in early exponential phase, but then decrease dramatically and are essentially absent in stationary phase. As a consequence of this variation, certain small RNAs, such as 6S RNA, remain low in exponential phase cells, and increase greatly in stationary phase. RNase BN affects 6S RNA abundance by decreasing its stability in exponential phase. RNase BN levels increase rapidly as cells exit stationary phase and are primarily responsible for the decrease in 6S RNA that accompanies this process. Purified RNase BN directly cleaves 6S RNA as shown by in vitro assays, and the 6S RNA:pRNA duplex is an even more favorable substrate of RNase BN. The exoribonuclease activity of RNase BN is unnecessary because all its action on 6S RNA is due to endonucleolytic cleavages. These data indicate that RNase BN plays an important role in determining levels of the global transcription regulator, 6S RNA, throughout the growth cycle.

Mutually exclusive RNA secondary structures regulate translation initiation of DinQ in Escherichia coli

Protein translation can be affected by changes in the secondary structure of mRNA. The dinQ gene in Escherichia coli encodes a primary transcript (+1) that is inert to translation. Ribonucleolytic removal of the 44 first nucleotides converts the +1 transcript into a translationally active form, but the mechanism behind this structural change is unknown. Here we present experimental evidence for a mechanism where alternative RNA secondary structures in the two dinQ mRNA variants affect translation initiation by mediating opening or closing of the ribosome binding sequence. This structural switch is determined by alternative interactions of four sequence elements within the dinQ mRNA and also by the agrB antisense RNA. Additionally, the structural conformation of +1 dinQ suggests a locking mechanism comprised of an RNA stem that both stabilizes and prevents translation initiation from the full-length dinQ transcript. BLAST search and multiple sequence alignments define a new family of dinQ-like genes widespread in Enterobacteriaceae with close RNA sequence similarities in their 5' untranslated regions. Thus, it appears that a whole new family of genes is regulated by the same mechanism of alternative secondary RNA structures.

The Conserved Dcw Gene Cluster of R. sphaeroides Is Preceded by an Uncommonly Extended 5' Leader Featuring the sRNA UpsM

Cell division and cell wall synthesis mechanisms are similarly conserved among bacteria. Consequently some bacterial species have comparable sets of genes organized in the dcw (division andcellwall) gene cluster. Dcw genes, their regulation and their relative order within the cluster are outstandingly conserved among rod shaped and gram negative bacteria to ensure an efficient coordination of growth and division. A well studied representative is the dcw gene cluster of E. coli. The first promoter of the gene cluster (mraZ1p) gives rise to polycistronic transcripts containing a 38 nt long 5’ UTR followed by the first gene mraZ. Despite reported conservation we present evidence for a much longer 5’ UTR in the gram negative and rod shaped bacterium Rhodobacter sphaeroides and in the family of Rhodobacteraceae. This extended 268 nt long 5’ UTR comprises a Rho independent terminator, which in case of termination gives rise to a non-coding RNA (UpsM). This sRNA is conditionally cleaved by RNase E under stress conditions in an Hfq- and very likely target mRNA-dependent manner, implying its function in trans. These results raise the question for the regulatory function of this extended 5’ UTR. It might represent the rarely described case of a trans acting sRNA derived from a riboswitch with exclusive presence in the family of Rhodobacteraceae.

Structure, Biology, and Therapeutic Application of Toxin-Antitoxin Systems in Pathogenic Bacteria

Bacterial toxin–antitoxin (TA) systems have received increasing attention for their diverse identities, structures, and functional implications in cell cycle arrest and survival against environmental stresses such as nutrient deficiency, antibiotic treatments, and immune system attacks. In this review, we describe the biological functions and the auto-regulatory mechanisms of six different types of TA systems, among which the type II TA system has been most extensively studied. The functions of type II toxins include mRNA/tRNA cleavage, gyrase/ribosome poison, and protein phosphorylation, which can be neutralized by their cognate antitoxins. We mainly explore the similar but divergent structures of type II TA proteins from 12 important pathogenic bacteria, including various aspects of protein–protein interactions. Accumulating knowledge about the structure–function correlation of TA systems from pathogenic bacteria has facilitated a novel strategy to develop antibiotic drugs that target specific pathogens. These molecules could increase the intrinsic activity of the toxin by artificially interfering with the intermolecular network of the TA systems.

Global RNA recognition patterns of post-transcriptional regulators Hfq and CsrA revealed by UV crosslinking in vivo

The molecular roles of many RNA‐binding proteins in bacterial post‐transcriptional gene regulation are not well understood. Approaches combining in vivo UV crosslinking with RNA deep sequencing (CLIP‐seq) have begun to revolutionize the transcriptome‐wide mapping of eukaryotic RNA‐binding protein target sites. We have applied CLIP‐seq to chart the target landscape of two major bacterial post‐transcriptional regulators, Hfq and CsrA, in the model pathogen Salmonella Typhimurium. By detecting binding sites at single‐nucleotide resolution, we identify RNA preferences and structural constraints of Hfq and CsrA during their interactions with hundreds of cellular transcripts. This reveals 3′‐located Rho‐independent terminators as a universal motif involved in Hfq–RNA interactions. Additionally, Hfq preferentially binds 5′ to sRNA‐target sites in mRNAs, and 3′ to seed sequences in sRNAs, reflecting a simple logic in how Hfq facilitates sRNA–mRNA interactions. Importantly, global knowledge of Hfq sites significantly improves sRNA‐target predictions. CsrA binds AUGGA sequences in apical loops and targets many Salmonella virulence mRNAs. Overall, our generic CLIP‐seq approach will bring new insights into post‐transcriptional gene regulation by RNA‐binding proteins in diverse bacterial species.

Linking bacterial type I toxins with their actions

Bacterial type I toxin-antitoxin systems consist of stable toxin-encoding mRNAs whose expression is counteracted by unstable RNA antitoxins. Accumulating evidence suggests that these players belong to broad regulatory networks influencing overall bacterial physiology. The majority of known transmembrane type I toxic peptides have conserved structural characteristics. However, recent studies demonstrated that their mechanisms of toxicity are diverse and complex. To better assess the current state of the art, type I toxins can be grouped into two classes according to their location and mechanisms of action: membrane-associated toxins acting by pore formation and/or by nucleoid condensation; and cytosolic toxins inducing nucleic acid cleavage. This classification will evolve as a result of future investigations.

The hmsT 3' untranslated region mediates c-di-GMP metabolism and biofilm formation in Yersinia pestis

Yersinia pestis, the cause of plague, forms a biofilm in the proventriculus of its flea vector to enhance transmission. Biofilm formation in Y. pestis is regulated by the intracellular levels of cyclic diguanylate (c-di-GMP). In this study, we investigated the role of the 3' untranslated region (3'UTR) in hmsT mRNA, a transcript that encodes a diguanylate cyclase that stimulates biofilm formation in Y. pestis by synthesizing the second messenger c-di-GMP. Deletion of the 3'UTR increased the half-life of hmsT mRNA, thereby upregulating c-di-GMP levels and biofilm formation. Our findings indicate that multiple regulatory sequences might be present in the hmsT 3'UTR that function together to mediate mRNA turnover. We also found that polynucleotide phosphorylase is partially responsible for hmsT 3'UTR-mediated mRNA decay. In addition, the hmsT 3'UTR strongly repressed gene expression at 37°C and 26°C, but affected gene expression only slightly at 21°C. Our findings suggest that the 3'UTR might be involved in precise and rapid regulation of hmsT expression, allowing Y. pestis to fine-tune c-di-GMP synthesis and consequently regulate biofilm production to adapt to the changing host environment.

Differential expression of small RNAs under chemical stress and fed-batch fermentation in E. coli

Background

Bacterial small RNAs (sRNAs) are recognized as posttranscriptional regulators involved in the control of bacterial lifestyle and adaptation to stressful conditions. Although chemical stress due to the toxicity of precursor and product compounds is frequently encountered in microbial bioprocessing applications, the involvement of sRNAs in this process is not well understood. We have used RNA sequencing to map sRNA expression in E. coli under chemical stress and high cell density fermentation conditions with the aim of identifying sRNAs involved in the transcriptional response and those with potential roles in stress tolerance.

Results

RNA sequencing libraries were prepared from RNA isolated from E. coli K-12 MG1655 cells grown under high cell density fermentation conditions or subjected to chemical stress with twelve compounds including four organic solvent-like compounds, four organic acids, two amino acids, geraniol and decanoic acid. We have discovered 253 novel intergenic transcripts with this approach, adding to the roughly 200 intergenic sRNAs previously reported in E. coli. There are eighty-four differentially expressed sRNAs during fermentation, of which the majority are novel, supporting possible regulatory roles for these transcripts in adaptation during different fermentation stages. There are a total of 139 differentially expressed sRNAs under chemical stress conditions, where twenty-nine exhibit significant expression changes in multiple tested conditions, suggesting that they may be involved in a more general chemical stress response. Among those with known functions are sRNAs involved in regulation of outer membrane proteins, iron availability, maintaining envelope homeostasis, as well as sRNAs incorporated into complex networks controlling motility and biofilm formation.

Conclusions

This study has used deep sequencing to reveal a wealth of hitherto undescribed sRNAs in E. coli and provides an atlas of sRNA expression during seventeen different growth and stress conditions. Although the number of novel sRNAs with regulatory functions is unknown, several exhibit specific expression patterns during high cell density fermentation and are differentially expressed in the presence of multiple chemicals, suggesting they may play regulatory roles during these stress conditions. These novel sRNAs, together with specific known sRNAs, are candidates for improving stress tolerance and our understanding of the E. coli regulatory network during fed-batch fermentation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2231-8) contains supplementary material, which is available to authorized users.

A novel small regulatory RNA enhances cell motility in enterohemorrhagic Escherichia coli

Small regulatory RNAs (sRNAs) are conserved among a wide range of bacteria. They modulate the translational efficiency of target mRNAs through base-pairing with the help of RNA chaperone Hfq. The present study identified a novel sRNA, Esr41 (enterohemorrhagic Escherichia coli O157 small RNA #41), from an intergenic region of an enterohemorrhagic E. coli (EHEC) O157:H7 Sakai-specific sequence that is not present in the nonpathogenic E. coli K-12. Esr41 was detected as an RNA molecule approximately 70 nucleotides long with a 3' GC-rich palindrome sequence followed by a long poly(U), which is a characteristic of rho-independent terminators and is also a structural feature required for the action of Hfq. EHEC O157 harboring a multicopy plasmid carrying the esr41 gene increased cell motility and the expression of fliC, a gene encoding a major flagellar component. These results indicate that Esr41 stimulates fliC expression in EHEC O157. Furthermore, the increase in cell motility induced by Esr41 was also observed in the E. coli K-12, suggesting that target genes controlled by Esr41 are present in both EHEC O157 and K-12.

sRNAs in bacterial type I and type III toxin-antitoxin systems

Toxin-antitoxin (TA) loci consist of two genes: a stable toxin whose overexpression kills the cell or causes growth stasis and an unstable antitoxin that neutralizes the toxin action. Currently, five TA systems are known. Here, we review type I and type III systems in which the antitoxins are regulatory RNAs. Type I antitoxins act by a base-pairing mechanism on toxin mRNAs. By contrast, type III antitoxins are RNA pseudoknots that bind their cognate toxins directly in an RNA-protein interaction. Whereas for a number of plasmid-encoded systems detailed information on structural requirements, kinetics of interaction with their targets and regulatory mechanisms employed by the antitoxin RNAs is available, the investigation of chromosomal systems is still in its infancy. Here, we summarize our current knowledge on that topic. Furthermore, we compare factors and conditions that induce antitoxins or toxins and different mechanisms of toxin action. Finally, we discuss biological roles for chromosome-encoded TA systems.

A cluster of four homologous small RNAs modulates C1 metabolism and the pyruvate dehydrogenase complex in Rhodobacter sphaeroides under various stress conditions

In bacteria, regulatory RNAs play an important role in the regulation and balancing of many cellular processes and stress responses. Among these regulatory RNAs, trans-encoded small RNAs (sRNAs) are of particular interest since one sRNA can lead to the regulation of multiple target mRNAs. In the purple bacterium Rhodobacter sphaeroides, several sRNAs are induced by oxidative stress. In this study, we focused on the functional characterization of four homologous sRNAs that are cotranscribed with the gene for the conserved hypothetical protein RSP_6037, a genetic arrangement described for only a few sRNAs until now. Each of the four sRNAs is characterized by two stem-loops that carry CCUCCUCCC motifs in their loops. They are induced under oxidative stress, as well as by various other stress conditions, and were therefore renamed here sRNAs CcsR1 to CcsR4 (CcsR1-4) for conserved CCUCCUCCC motif stress-induced RNAs 1 to 4. Increased CcsR1-4 expression decreases the expression of genes involved in C1 metabolism or encoding components of the pyruvate dehydrogenase complex either directly by binding to their target mRNAs or indirectly. One of the CcsR1-4 target mRNAs encodes the transcriptional regulator FlhR, an activator of glutathione-dependent methanol/formaldehyde metabolism. Downregulation of this glutathione-dependent pathway increases the pool of glutathione, which helps to counteract oxidative stress. The FlhR-dependent downregulation of the pyruvate dehydrogenase complex reduces a primary target of reactive oxygen species and reduces aerobic electron transport, a main source of reactive oxygen species. Our findings reveal a previously unknown strategy used by bacteria to counteract oxidative stress.

IMPORTANCE

Phototrophic organisms have to cope with photo-oxidative stress due to the function of chlorophylls as photosensitizers for the formation of singlet oxygen. Our study assigns an important role in photo-oxidative stress resistance to a cluster of four homologous sRNAs in the anoxygenic phototrophic bacterium Rhodobacter sphaeroides. We reveal a function of these regulatory RNAs in the fine-tuning of C1 metabolism. A model that relates oxidative stress defense to C1 metabolism is presented.

Regulatory small RNAs from the 3' regions of bacterial mRNAs

Most studies of small regulatory RNAs in bacteria have focussed on conserved transcripts in intergenic regions. However, several recent developments including single-nucleotide resolution transcriptome profiling by RNA-seq and increased knowledge of the cellular targets of the RNA chaperone Hfq suggest that the bacterial world of functional small RNAs is more diverse. One emerging class are small RNAs that are identical to the 3' regions of known mRNAs, but are produced either by transcription from internal promoters or by mRNA processing. Using several recently discovered examples of such sRNAs, we discuss their biogenesis and modes of action, and illustrate how they can facilitate mRNA crosstalk in various physiological processes.

Cut to the chase--Regulating translation through RNA cleavage

Activation of toxin-antitoxin (TA) systems provides an important mechanism for bacteria to adapt to challenging and ever changing environmental conditions. Known TA systems are classified into five families based on the mechanisms of antitoxin inhibition and toxin activity. For type II TA systems, the toxin is inactivated in exponentially growing cells by tightly binding its antitoxin partner protein, which also serves to regulate cellular levels of the complex through transcriptional auto-repression. During cellular stress, however, the antitoxin is degraded thus freeing the toxin, which is then able to regulate central cellular processes, primarily protein translation to adjust cell growth to the new conditions. In this review, we focus on the type II TA pairs that regulate protein translation through cleavage of ribosomal, transfer, or messenger RNA.

The extracellular RNA complement of Escherichia coli

The secretion of biomolecules into the extracellular milieu is a common and well-conserved phenomenon in biology. In bacteria, secreted biomolecules are not only involved in intra-species communication but they also play roles in inter-kingdom exchanges and pathogenicity. To date, released products, such as small molecules, DNA, peptides, and proteins, have been well studied in bacteria. However, the bacterial extracellular RNA complement has so far not been comprehensively characterized. Here, we have analyzed, using a combination of physical characterization and high-throughput sequencing, the extracellular RNA complement of both outer membrane vesicle (OMV)-associated and OMV-free RNA of the enteric Gram-negative model bacterium Escherichia coli K-12 substrain MG1655 and have compared it to its intracellular RNA complement. Our results demonstrate that a large part of the extracellular RNA complement is in the size range between 15 and 40 nucleotides and is derived from specific intracellular RNAs. Furthermore, RNA is associated with OMVs and the relative abundances of RNA biotypes in the intracellular, OMV and OMV-free fractions are distinct. Apart from rRNA fragments, a significant portion of the extracellular RNA complement is composed of specific cleavage products of functionally important structural noncoding RNAs, including tRNAs, 4.5S RNA, 6S RNA, and tmRNA. In addition, the extracellular RNA pool includes RNA biotypes from cryptic prophages, intergenic, and coding regions, of which some are so far uncharacterised, for example, transcripts mapping to the fimA-fimL and ves-spy intergenic regions. Our study provides the first detailed characterization of the extracellular RNA complement of the enteric model bacterium E. coli. Analogous to findings in eukaryotes, our results suggest the selective export of specific RNA biotypes by E. coli, which in turn indicates a potential role for extracellular bacterial RNAs in intercellular communication.

Genome-wide analyses in bacteria show small-RNA enrichment for long and conserved intergenic regions

Interest in finding small RNAs (sRNAs) in bacteria has significantly increased in recent years due to their regulatory functions. Development of high-throughput methods and more sophisticated computational algorithms has allowed rapid identification of sRNA candidates in different species. However, given their various sizes (50 to 500 nucleotides [nt]) and their potential genomic locations in the 5' and 3' untranslated regions as well as in intergenic regions, identification and validation of true sRNAs have been challenging. In addition, the evolution of bacterial sRNAs across different species continues to be puzzling, given that they can exert similar functions with various sequences and structures. In this study, we analyzed the enrichment patterns of sRNAs in 13 well-annotated bacterial species using existing transcriptome and experimental data. All intergenic regions were analyzed by WU-BLAST to examine conservation levels relative to species within or outside their genus. In total, more than 900 validated bacterial sRNAs and 23,000 intergenic regions were analyzed. The results indicate that sRNAs are enriched in intergenic regions, which are longer and more conserved than the average intergenic regions in the corresponding bacterial genome. We also found that sRNA-coding regions have different conservation levels relative to their flanking regions. This work provides a way to analyze how noncoding RNAs are distributed in bacterial genomes and also shows conserved features of intergenic regions that encode sRNAs. These results also provide insight into the functions of regions surrounding sRNAs and into optimization of RNA search algorithms.

NAD captureSeq indicates NAD as a bacterial cap for a subset of regulatory RNAs

A distinctive feature of prokaryotic gene expression is the absence of 5'-capped RNA. In eukaryotes, 5',5'-triphosphate-linked 7-methylguanosine protects messenger RNA from degradation and modulates maturation, localization and translation. Recently, the cofactor nicotinamide adenine dinucleotide (NAD) was reported as a covalent modification of bacterial RNA. Given the central role of NAD in redox biochemistry, posttranslational protein modification and signalling, its attachment to RNA indicates that there are unknown functions of RNA in these processes and undiscovered pathways in RNA metabolism and regulation. The unknown identity of NAD-modified RNAs has so far precluded functional analyses. Here we identify NAD-linked RNAs from bacteria by chemo-enzymatic capture and next-generation sequencing (NAD captureSeq). Among those identified, specific regulatory small RNAs (sRNAs) and sRNA-like 5'-terminal fragments of certain mRNAs are particularly abundant. Analogous to a eukaryotic cap, 5'-NAD modification is shown in vitro to stabilize RNA against 5'-processing by the RNA-pyrophosphohydrolase RppH and against endonucleolytic cleavage by ribonuclease (RNase) E. The nudix phosphohydrolase NudC decaps NAD-RNA and thereby triggers RNase-E-mediated RNA decay, while being inactive against triphosphate-RNA. In vivo, ∼13% of the abundant sRNA RNAI is NAD-capped in the presence, and ∼26% in the absence, of functional NudC. To our knowledge, this is the first description of a cap-like structure and a decapping machinery in bacteria.

RNase J is required for processing of a small number of RNAs in Rhodobacter sphaeroides

All bacteria contain multiple exoribonucleases to ensure a fast breakdown of different RNA molecules, either for maturation or for complete degradation to the level of mononucleotides. This efficient RNA degradation plays pivotal roles in the post-transcriptional gene regulation, in RNA processing and maturation as well as in RNA quality control mechanisms and global adaption to stress conditions. Besides different 3'-to-5' exoribonucleases mostly with overlapping functions in vivo many bacteria additionally possess the 5'-to-3' exoribonuclease, RNase J, to date the only known bacterial ribonuclease with this activity. An RNA-seq approach was applied to identify specific targets of RNase J in the α-proteobacterium Rhodobacter sphaeroides. Only few transcripts were strongly affected by the lack of RNase J implying that its function is mostly required for specific processing/degradation steps in this bacterium. The accumulation of diverse RNA fragments in the RNase J deletion mutant points to RNA features that apparently cannot be targeted by the conventional 3'-exoribonucleases in Gram-negative bacteria.

Identification of non-coding RNAs with a new composite feature in the Hybrid Random Forest Ensemble algorithm

To identify non-coding RNA (ncRNA) signals within genomic regions, a classification tool was developed based on a hybrid random forest (RF) with a logistic regression model to efficiently discriminate short ncRNA sequences as well as long complex ncRNA sequences. This RF-based classifier was trained on a well-balanced dataset with a discriminative set of features and achieved an accuracy, sensitivity and specificity of 92.11%, 90.7% and 93.5%, respectively. The selected feature set includes a new proposed feature, SCORE. This feature is generated based on a logistic regression function that combines five significant features-structure, sequence, modularity, structural robustness and coding potential-to enable improved characterization of long ncRNA (lncRNA) elements. The use of SCORE improved the performance of the RF-based classifier in the identification of Rfam lncRNA families. A genome-wide ncRNA classification framework was applied to a wide variety of organisms, with an emphasis on those of economic, social, public health, environmental and agricultural significance, such as various bacteria genomes, the Arthrospira (Spirulina) genome, and rice and human genomic regions. Our framework was able to identify known ncRNAs with sensitivities of greater than 90% and 77.7% for prokaryotic and eukaryotic sequences, respectively. Our classifier is available at http://ncrna-pred.com/HLRF.htm.

Discovery of ethanol-responsive small RNAs in Zymomonas mobilis

Zymomonas mobilis is a bacterium that can produce ethanol by fermentation. Due to its unique metabolism and efficient ethanol production, Z. mobilis has attracted special interest for biofuel energy applications; an important area of study is the regulation of those specific metabolic pathways. Small RNAs (sRNAs) have been studied as molecules that function as transcriptional regulators in response to cellular stresses. While sRNAs have been discovered in various organisms by computational prediction and experimental approaches, their discovery in Z. mobilis has not yet been reported. In this study, we have applied transcriptome analysis and computational predictions to facilitate identification and validation of 15 novel sRNAs in Z. mobilis. We furthermore characterized their expression in the context of high and low levels of intracellular ethanol. Here, we report that 3 of the sRNAs (Zms2, Zms4, and Zms6) are differentially expressed under aerobic and anaerobic conditions, when low and high ethanol productions are observed, respectively. Importantly, when we tested the effect of ethanol stress on the expression of sRNAs in Z. mobilis, Zms2, Zms6, and Zms18 showed differential expression under 5% ethanol stress conditions. These data suggest that in this organism regulatory RNAs can be associated with metabolic functions involved in ethanol stress responses.

A eukaryotic-like 3' untranslated region in Salmonella enterica hilD mRNA

Long 3' untranslated regions (3'UTRs) are common in eukaryotic mRNAs. In contrast, long 3'UTRs are rare in bacteria, and have not been characterized in detail. We describe a 3'UTR of 310 nucleotides in hilD mRNA, a transcript that encodes a transcriptional activator of Salmonella enterica pathogenicity island 1 (SPI-1). Deletion of the hilD 3'UTR increases the hilD mRNA level, suggesting that the hilD 3'UTR may play a role in hilD mRNA turnover. Cloning of the hilD 3'UTR downstream of the green fluorescent protein (gfp) gene decreases green fluorescent protein (GFP) activity in both Escherichia coli and S. enterica, indicating that the hilD 3'UTR can act as an independent module. S. enterica mutants lacking either ribonuclease E or polynucleotide phosphorylase contain similar amounts of hilD and hilD Δ3'UTR mRNAs, suggesting that the hilD 3'UTR is a target for hilD mRNA degradation by the degradosome. The hilD 3'UTR is also necessary for modulation of hilD and SPI-1 expression by the RNA chaperone Hfq. Overexpression of SPI-1 in the absence of the hilD 3'UTR retards Salmonella growth and causes uncontrolled invasion of epithelial cells. Based on these observations, we propose that the S. enterica hilD 3'UTR is a cis-acting element that contributes to cellular homeostasis by promoting hilD mRNA turnover.

Genome-wide profiling of Hfq-binding RNAs uncovers extensive post-transcriptional rewiring of major stress response and symbiotic regulons in Sinorhizobium meliloti

The RNA chaperone Hfq is a global post-transcriptional regulator in bacteria. Here, we used RNAseq to analyze RNA populations from the legume symbiont Sinorhizobium meliloti that were co-immunoprecipitated (CoIP-RNA) with a FLAG-tagged Hfq in five growth/stress conditions. Hfq-bound transcripts (1315) were largely identified in stressed bacteria and derived from small RNAs (sRNAs), both trans-encoded (6.4%) and antisense (asRNAs; 6.3%), and mRNAs (86%). Pull-down with Hfq recovered a small proportion of annotated S. meliloti sRNAs (14% of trans-sRNAs and 2% of asRNAs) suggesting a discrete impact of this protein in sRNA pathways. Nonetheless, Hfq selectively stabilized CoIP-enriched sRNAs, anticipating that these interactions are functionally significant. Transcription of 26 Hfq-bound sRNAs was predicted to occur from promoters recognized by the major stress σ factors σ(E2) or σ(H1/2). Recovery rates of sRNAs in each of the CoIP-RNA libraries suggest a large impact of Hfq-assisted riboregulation in S. meliloti osmoadaptation. Hfq directly targeted 18% of the predicted S. meliloti mRNAs, which encode functionally diverse proteins involved in transport and metabolism, σ(E2)-dependent stress responses, quorum sensing, flagella biosynthesis, ribosome, and membrane assembly or symbiotic nitrogen fixation. Canonical targeting of the 5' regions of two of the ABC transporter mRNAs by the homologous Hfq-binding AbcR1 and AbcR2 sRNAs leading to inhibition of protein synthesis was confirmed in vivo. We therefore provide a comprehensive resource for the systems-level deciphering of hitherto unexplored S. meliloti stress and symbiotic post-transcriptional regulons and the identification of Hfq-dependent sRNA-mRNA regulatory pairs.

Base pairing interaction between 5'- and 3'-UTRs controls icaR mRNA translation in Staphylococcus aureus

The presence of regulatory sequences in the 3′ untranslated region (3′-UTR) of eukaryotic mRNAs controlling RNA stability and translation efficiency is widely recognized. In contrast, the relevance of 3′-UTRs in bacterial mRNA functionality has been disregarded. Here, we report evidences showing that around one-third of the mapped mRNAs of the major human pathogen Staphylococcus aureus carry 3′-UTRs longer than 100-nt and thus, potential regulatory functions. We selected the long 3′-UTR of icaR, which codes for the repressor of the main exopolysaccharidic compound of the S. aureus biofilm matrix, to evaluate the role that 3′-UTRs may play in controlling mRNA expression. We showed that base pairing between the 3′-UTR and the Shine-Dalgarno (SD) region of icaR mRNA interferes with the translation initiation complex and generates a double-stranded substrate for RNase III. Deletion or substitution of the motif (UCCCCUG) within icaR 3′-UTR was sufficient to abolish this interaction and resulted in the accumulation of IcaR repressor and inhibition of biofilm development. Our findings provide a singular example of a new potential post-transcriptional regulatory mechanism to modulate bacterial gene expression through the interaction of a 3′-UTR with the 5′-UTR of the same mRNA.

At both sides of the protein-coding region, the mRNA molecule contains sequences that are not translated to protein. In eukaryotes, the untranslated 3′ region (3′-UTR), which comprises from the last codon used in translation to the 3′ end of the mRNA, controls mRNA stability, location and translation efficiency. In contrast, knowledge about the functions of 3′-UTRs in bacterial physiology is scarce. Here, we demonstrate that bacterial 3′-UTRs might play regulatory functions that might resemble those already described in eukaryotes. Transcriptome analysis of the human pathogen Staphylococcus aureus revealed that at least 30% of mRNAs contain long 3′-UTRs. Using the 3′-UTR of the mRNA encoding the main biofilm repressor IcaR as a model, we show that the 3′-UTR interferes with the translation initiation complex and promotes mRNA decay through base pairing with the ribosome binding site. This event contributes to adjusting IcaR level and modulating exopolysaccharide production and biofilm development in S. aureus. Our data illustrate that bacterial 3′-UTRs can provide strategies for fine-tuning control of gene expression.

Constitutive expression of RyhB regulates the heme biosynthesis pathway and increases the 5-aminolevulinic acid accumulation in Escherichia coli

In the current study, the small RNA ryhB, which regulates the metabolism of iron in Escherichia coli, was constitutively expressed in engineered E. coli DALA. The resulting strain E. coli DALRA produced 16% more 5-aminolevulinic acid (ALA) than the parent strain E. coli DALA in batch fermentation. Meanwhile, we found that addition of iron in the medium increased heme formation and reduced ALA yield, whereas the presence of iron chelator in the medium decreased heme concentration and increased the ALA production efficiency (ALA yield per OD600). The qRT-PCR analysis showed that the mRNA levels of hemB and hemH were also decreased as well as the known RyhB target genes of acnAB, sdhAB, fumA, and cydAB in E. coli DALRA. These results indicated that small RNA can be used as a tool for regulating ALA accumulation in E. coli.

Identification of novel sRNAs in mycobacterial species

Bacterial small RNAs (sRNAs) are short transcripts that typically do not encode proteins and often act as regulators of gene expression through a variety of mechanisms. Regulatory sRNAs have been identified in many species, including Mycobacterium tuberculosis, the causative agent of tuberculosis. Here, we use a computational algorithm to predict sRNA candidates in the mycobacterial species M. smegmatis and M. bovis BCG and confirmed the expression of many sRNAs using Northern blotting. Thus, we have identified 17 and 23 novel sRNAs in M. smegmatis and M. bovis BCG, respectively. We have also applied a high-throughput technique (Deep-RACE) to map the 5' and 3' ends of many of these sRNAs and identified potential regulators of sRNAs by analysis of existing ChIP-seq datasets. The sRNAs identified in this work likely contribute to the unique biology of mycobacteria.

Prediction and characterization of small non-coding RNAs related to secondary metabolites in Saccharopolyspora erythraea

Saccharopolyspora erythraea produces a large number of secondary metabolites with biological activities, including erythromycin. Elucidation of the mechanisms through which the production of these secondary metabolites is regulated may help to identify new strategies for improved biosynthesis of erythromycin. In this paper, we describe the systematic prediction and analysis of small non-coding RNAs (sRNAs) in S. erythraea, with the aim to elucidate sRNA-mediated regulation of secondary metabolite biosynthesis. In silico and deep-sequencing technologies were applied to predict sRNAs in S. erythraea. Six hundred and forty-seven potential sRNA loci were identified, of which 382 cis-encoded antisense RNA are complementary to protein-coding regions and 265 predicted transcripts are located in intergenic regions. Six candidate sRNAs (sernc292, sernc293, sernc350, sernc351, sernc361, and sernc389) belong to four gene clusters (tpc3, pke, pks6, and nrps5) that are involved in secondary metabolite biosynthesis. Deep-sequencing data showed that the expression of all sRNAs in the strain HL3168 E3 (E3) was higher than that in NRRL23338 (M), except for sernc292 and sernc361 expression. The relative expression of six sRNAs in strain M and E3 were validated by qRT-PCR at three different time points (24, 48, and 72 h). The results showed that, at each time point, the transcription levels of sernc293, sernc350, sernc351, and sernc389 were higher in E3 than in M, with the largest difference observed at 72 h, whereas no signals for sernc292 and sernc361 were detected. sernc293, sernc350, sernc351, and sernc389 probably regulate iron transport, terpene metabolism, geosmin synthesis, and polyketide biosynthesis, respectively. The major significance of this study is the successful prediction and identification of sRNAs in genomic regions close to the secondary metabolism-related genes in S. erythraea. A better understanding of the sRNA-target interaction would help to elucidate the complete range of functions of sRNAs in S. erythraea, including sRNA-mediated regulation of erythromycin biosynthesis.

Bacterial Mg2+ homeostasis, transport, and virulence

Organisms must maintain physiological levels of Mg(2+) because this divalent cation is critical for the stabilization of membranes and ribosomes, for the neutralization of nucleic acids, and as a cofactor in a variety of enzymatic reactions. In this review, we describe the mechanisms that bacteria utilize to sense the levels of Mg(2+) both outside and inside the cytoplasm. We examine how bacteria achieve Mg(2+) homeostasis by adjusting the expression and activity of Mg(2+) transporters and by changing the composition of their cell envelope. We discuss the connections that exist between Mg(2+) sensing, Mg(2+) transport, and bacterial virulence. Additionally, we explore the logic behind the fact that bacterial genomes encode multiple Mg(2+) transporters and distinct sensing systems for cytoplasmic and extracytoplasmic Mg(2+). These analyses may be applicable to the homeostatic control of other cations.

Comparative analysis of non-coding RNAs in the antibiotic-producing Streptomyces bacteria

Background

Non-coding RNAs (ncRNAs) are key regulatory elements that control a wide range of cellular processes in all bacteria in which they have been studied. Taking advantage of recent technological innovations, we set out to fully explore the ncRNA potential of the multicellular, antibiotic-producing Streptomyces bacteria.

Results

Using a comparative RNA sequencing analysis of three divergent model streptomycetes (S. coelicolor, S. avermitilis and S. venezuelae), we discovered hundreds of novel cis-antisense RNAs and intergenic small RNAs (sRNAs). We identified a ubiquitous antisense RNA species that arose from the overlapping transcription of convergently-oriented genes; we termed these RNA species ‘cutoRNAs’, for convergent untranslated overlapping RNAs. Conservation between different classes of ncRNAs varied greatly, with sRNAs being more conserved than antisense RNAs. Many species-specific ncRNAs, including many distinct cutoRNA pairs, were located within antibiotic biosynthetic clusters, including the actinorhodin, undecylprodigiosin, and coelimycin clusters of S. coelicolor, the chloramphenicol cluster of S. venezuelae, and the avermectin cluster of S. avermitilis.

Conclusions

These findings indicate that ncRNAs, including a novel class of antisense RNA, may exert a previously unrecognized level of regulatory control over antibiotic production in these bacteria. Collectively, this work has dramatically expanded the ncRNA repertoire of three Streptomyces species and has established a critical foundation from which to investigate ncRNA function in this medically and industrially important bacterial genus.

Prediction of candidate small non-coding RNAs in Agrobacterium by computational analysis

Small non-coding RNAs with important regulatory roles are not confined to eukaryotes. Recent work has uncovered a growing number of bacterial small RNAs (sRNAs), some of which have been shown to regulate critical cellular processes. Computational approaches, in combination with molecular experiments, have played an important role in the identification of these sRNAs. At present, there is no information on the presence of small non-coding RNAs and their genes in the Agrobacterium tumefaciens genome. To identify potential sRNAs in this important bacterium, deep sequencing of the short RNA populations isolated from Agrobacterium tumefaciens C58 was carried out. From a data set of more than 10,000 short sequences, 16 candidate sRNAs have been tentatively identified based on computational analysis. All of these candidates can form stem-loop structures by RNA folding predictions and the majority of the secondary structures are rich in GC base pairs. Some are followed by a short stretch of U residues, indicative of a rho-independent transcription terminator, whereas some of the short RNAs are found in the stem region of the hairpin, indicative of eukaryotic-like sRNAs. Experimental strategies will need to be used to verify these candidates. The study of an expanded list of candidate sRNAs in Agrobacterium will allow a more complete understanding of the range of roles played by regulatory RNAs in prokaryotes.

Whole-Genome Sequencing and Comparative Analysis of Yersinia pestis, the Causative Agent of a Plague Outbreak in Northern Peru

The plague is a zoonotic disease caused by the bacterium Yersinia pestis. Here, we report the complete genome sequence of the Y. pestis strain INS, which was isolated from swollen lymph gland aspirate (bubo aspirate) of an infected patient from a pneumonic outbreak in 2010 in northern Peru.

Structural overview of toxin-antitoxin systems in infectious bacteria: a target for developing antimicrobial agents

The bacterial toxin-antitoxin (TA) system is a module that may play a role in cell survival under stress conditions. Generally, toxin molecules act as negative regulators in cell survival and antitoxin molecules as positive regulators. Thus, the expression levels and interactions between toxins and antitoxins should be systematically harmonized so that bacteria can escape such harmful conditions. Since TA systems are able to control the fate of bacteria, they are considered potent targets for the development of new antimicrobial agents. TA systems are widely prevalent with a variety of systems existing in bacteria: there are three types of bacterial TA systems depending on the property of the antitoxin which binds either the protein toxin or mRNA coding the toxin protein. Moreover, the multiplicity of TA genes has been observed even in species of bacteria. Therefore, knowledge on TA systems such as the individual characteristics of TA systems, integrative working mechanisms of various TA systems in bacteria, interactions between toxin molecules and cellular targets, and so on is currently limited due to their complexity. In this regard, it would be helpful to know the structural characteristics of TA modules for understanding TA systems in bacteria. Until now, 85 out of the total structures deposited in PDB have been bacterial TA system proteins including TA complexes or isolated toxins/antitoxins. Here, we summarized the structural information of TA systems and analyzed the structural characteristics of known TA modules from several bacteria, especially focusing on the TA modules of several infectious bacteria.

New type I toxin-antitoxin families from "wild" and laboratory strains of E. coli: Ibs-Sib, ShoB-OhsC and Zor-Orz

Type I toxin-antitoxin loci consist of two genes: one encodes a small, toxic protein and the second encodes a small RNA antitoxin that represses toxin gene expression. These pairs were first described on plasmids where they regulate plasmid maintenance. However, recent discoveries have found novel type I loci, with no homology to plasmid sequences, in the chromosome of Escherichia coli and closely related species. The Ibs-Sib, ShoB-OhsC and Zor-Orz loci are examples of these new loci. For these toxic proteins, much more is known about how their expression is regulated than their biological function. Although all are found in E. coli and closely related bacteria, there is great variation among species as to which loci they possess. Herein, I discuss how these sRNA antitoxins prevent toxin production and how the distribution of these loci across species may be providing insights into their true function.

Divergently overlapping cis-encoded antisense RNA regulating toxin-antitoxin systems from E. coli: hok/sok, ldr/rdl, symE/symR

Toxin-antitoxin (TA) systems are categorized into three classes based on the type of antitoxin. In type I TA systems, the antitoxin is a small antisense RNA that inhibits translation of small toxic proteins by binding to the corresponding mRNAs. Those type I TA systems were originally identified as plasmid stabilization modules rendering a post-segregational killing (PSK) effect on the host cells. The type I TA loci also exist on the Escherichia coli chromosome but their biological functions are less clear. Genetic organization and regulatory elements of hok/sok and ldr/rdl families are very similar and the toxins are predicted to contain a transmembrane domain, but otherwise share no detectable sequence similarity. This review will give an overview of the type I TA modules of E. coli K-12, especially hok/sok, ldr/rdl and SOS-inducible symE/symR systems, which are regulated by divergently overlapping cis-encoded antisense RNAs.

An atlas of Hfq-bound transcripts reveals 3' UTRs as a genomic reservoir of regulatory small RNAs

The small RNAs associated with the protein Hfq constitute one of the largest classes of post-transcriptional regulators known to date. Most previously investigated members of this class are encoded by conserved free-standing genes. Here, deep sequencing of Hfq-bound transcripts from multiple stages of growth of Salmonella typhimurium revealed a plethora of new small RNA species from within mRNA loci, including DapZ, which overlaps with the 3' region of the biosynthetic gene, dapB. Synthesis of the DapZ small RNA is independent of DapB protein synthesis, and is controlled by HilD, the master regulator of Salmonella invasion genes. DapZ carries a short G/U-rich domain similar to that of the globally acting GcvB small RNA, and uses GcvB-like seed pairing to repress translation of the major ABC transporters, DppA and OppA. This exemplifies double functional output from an mRNA locus by the production of both a protein and an Hfq-dependent trans-acting RNA. Our atlas of Hfq targets suggests that the 3' regions of mRNA genes constitute a rich reservoir that provides the Hfq network with new regulatory small RNAs.

Deep sequencing of RNA from immune cell-derived vesicles uncovers the selective incorporation of small non-coding RNA biotypes with potential regulatory functions

Cells release RNA-carrying vesicles and membrane-free RNA/protein complexes into the extracellular milieu. Horizontal vesicle-mediated transfer of such shuttle RNA between cells allows dissemination of genetically encoded messages, which may modify the function of target cells. Other studies used array analysis to establish the presence of microRNAs and mRNA in cell-derived vesicles from many sources. Here, we used an unbiased approach by deep sequencing of small RNA released by immune cells. We found a large variety of small non-coding RNA species representing pervasive transcripts or RNA cleavage products overlapping with protein coding regions, repeat sequences or structural RNAs. Many of these RNAs were enriched relative to cellular RNA, indicating that cells destine specific RNAs for extracellular release. Among the most abundant small RNAs in shuttle RNA were sequences derived from vault RNA, Y-RNA and specific tRNAs. Many of the highly abundant small non-coding transcripts in shuttle RNA are evolutionary well-conserved and have previously been associated to gene regulatory functions. These findings allude to a wider range of biological effects that could be mediated by shuttle RNA than previously expected. Moreover, the data present leads for unraveling how cells modify the function of other cells via transfer of specific non-coding RNA species.

Comparative profiling of Pseudomonas aeruginosa strains reveals differential expression of novel unique and conserved small RNAs

Pseudomonas aeruginosa is a highly adaptable bacterium that thrives in a broad range of ecological niches and can infect multiple hosts as diverse as plants, nematodes and mammals. In humans, it is an important opportunistic pathogen. This wide adaptability correlates with its broad genetic diversity. In this study, we used a deep-sequencing approach to explore the complement of small RNAs (sRNAs) in P. aeruginosa as the number of such regulatory molecules previously identified in this organism is relatively low, considering its genome size, phenotypic diversity and adaptability. We have performed a comparative analysis of PAO1 and PA14 strains which share the same host range but differ in virulence, PA14 being considerably more virulent in several model organisms. Altogether, we have identified more than 150 novel candidate sRNAs and validated a third of them by Northern blotting. Interestingly, a number of these novel sRNAs are strain-specific or showed strain-specific expression, strongly suggesting that they could be involved in determining specific phenotypic traits.

A genome-wide identification analysis of small regulatory RNAs in Mycobacterium tuberculosis by RNA-Seq and conservation analysis

We propose a new method for smallRNAs (sRNAs) identification. First we build an effective target genome (ETG) by means of a strand-specific procedure. Then we propose a new bioinformatic pipeline based mainly on the combination of two types of information: the first provides an expression map based on RNA-seq data (Reads Map) and the second applies principles of comparative genomics leading to a Conservation Map. By superimposing these two maps, a robust method for the search of sRNAs is obtained. We apply this methodology to investigate sRNAs in Mycobacterium tuberculosis H37Rv. This bioinformatic procedure leads to a total list of 1948 candidate sRNAs. The size of the candidate list is strictly related to the aim of the study and to the technology used during the verification process. We provide performance measures of the algorithm in identifying annotated sRNAs reported in three recent published studies.

Dynamics of Salmonella small RNA expression in non-growing bacteria located inside eukaryotic cells

Small non-coding regulatory RNAs (sRNAs) have been studied in many bacterial pathogens during infection. However, few studies have focused on how intracellular pathogens modulate sRNA expression inside eukaryotic cells. Here, we monitored expression of all known sRNAs of Salmonella enterica serovar Typhimurium (S. Typhimurium) in bacteria located inside fibroblasts, a host cell type in which this pathogen restrains growth. sRNA sequences known in S. Typhimurium and Escherichia coli were searched in the genome of S. Typhimurium virulent strain SL1344, the subject of this study. Expression of 84 distinct sRNAs was compared in extra- and intracellular bacteria. Non-proliferating intracellular bacteria upregulated six sRNAs, including IsrA, IsrG, IstR-2, RyhB-1, RyhB-2 and RseX while repressed the expression of the sRNAs DsrA, GlmZ, IsrH-1, IsrI, SraL, SroC, SsrS(6S) and RydC. Interestingly, IsrH-1 was previously reported as an sRNA induced by S. Typhimurium inside macrophages. Kinetic analyses unraveled changing expression patterns for some sRNAs along the infection. InvR and T44 expression dropped after an initial induction phase while IstR-2 was induced exclusively at late infection times (> 6 h). Studies focused on the Salmonella-specific sRNA RyhB-2 revealed that intracellular bacteria use this sRNA to regulate negatively YeaQ, a cis-encoded protein of unknown function. RyhB-2, together with RyhB-1, contributes to attenuate intracellular bacterial growth. To our knowledge, these data represent the first comprehensive study of S. Typhimurium sRNA expression in intracellular bacteria and provide the first insights into sRNAs that may direct pathogen adaptation to a non-proliferative state inside the host cell.