Evolutionary dynamics of small RNAs in 27 Escherichia coli and Shigella genomes

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.

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.

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.

Natural selection shapes maintenance of orthologous sRNAs in divergent host-restricted bacterial genomes

Historically it has been difficult to study the evolution of bacterial small RNAs (sRNAs) across distantly related species. For example, identifying homologs of sRNAs is often difficult in genomes that have undergone multiple structural rearrangements. Also, some types of regulatory sRNAs evolve at rapid rates. The high degree of genomic synteny among divergent host-restricted bacterial lineages, including intracellular symbionts, is conducive to sRNA maintenance and homolog identification. In turn, symbiont genomes can provide us with novel insights into sRNA evolution. Here, we examine the sRNA expression profile of the obligate symbiont of psyllids, Carsonella ruddii, which has one of the smallest cellular genomes described. Using RNA-seq, we identified 36 and 32 antisense sRNAs (asRNAs) expressed by Carsonella from the psyllids Bactericera cockerelli (Carsonella-BC) and Diaphorina citri (Carsonella-DC), respectively. The majority of these asRNAs were associated with genes that are involved in essential amino acid biosynthetic pathways. Eleven of the asRNAs were conserved in both Carsonella lineages and the majority were maintained by selection. Notably, five of the corresponding coding sequences are also the targets of conserved asRNAs in a distantly related insect symbiont, Buchnera. We detected differential expression of two asRNAs for genes involved in arginine and leucine biosynthesis occurring between two distinct Carsonella-BC life stages. Using asRNAs identified in Carsonella, Buchnera, and Profftella which are all endosymbionts, and Escherichia coli, we determined that regions upstream of these asRNAs encode unique conserved patterns of AT/GC richness, GC skew, and sequence motifs which may be involved in asRNA regulation.

Metabolic remodeling in Escherichia coli MG1655. A prophage e14-encoded small RNA, co293, post-transcriptionally regulates transcription factors HcaR and FadR

In previous studies, we have shown the existence of metabolic remodeling in glucose-grown Escherichia coli MG1655 cells expressing the esterase Orf306 from the opd island of Sphingobium fuliginis. We now show that Orf306-dependent metabolic remodeling is due to regulation of a novel small RNA (sRNA). Endogenous propionate, produced due to the esterase/lipase activity of Orf306, repressed expression of a novel E. coli sRNA, co293. This sRNA post-transcriptionally regulates expression of the transcription factors HcaR and FadR either by inhibiting translation or by destabilizing their transcripts. Hence, repression of co293 expression elevates the levels of HcaR and FadR with consequent activation of alternative carbon catabolic pathways. HcaR activates the hca and MHP operons leading to upregulation of the phenyl propionate and hydroxy phenyl propionate (HPP) degradation pathways. Similarly, FadR stimulates the expression of the transcription factor IclR which negatively regulates the glyoxylate bypass pathway genes, aceBAK.

RNA-Sequencing Analyses of Small Bacterial RNAs and their Emergence as Virulence Factors in Host-Pathogen Interactions

Proteins have long been considered to be the most prominent factors regulating so-called invasive genes involved in host-pathogen interactions. The possible role of small non-coding RNAs (sRNAs), either intracellular, secreted or packaged in outer membrane vesicles (OMVs), remained unclear until recently. The advent of high-throughput RNA-sequencing (RNA-seq) techniques has accelerated sRNA discovery. RNA-seq radically changed the paradigm on bacterial virulence and pathogenicity to the point that sRNAs are emerging as an important, distinct class of virulence factors in both gram-positive and gram-negative bacteria. The potential of OMVs, as protectors and carriers of these functional, gene regulatory sRNAs between cells, has also provided an additional layer of complexity to the dynamic host-pathogen relationship. Using a non-exhaustive approach and through examples, this review aims to discuss the involvement of sRNAs, either free or loaded in OMVs, in the mechanisms of virulence and pathogenicity during bacterial infection. We provide a brief overview of sRNA origin and importance, and describe the classical and more recent methods of identification that have enabled their discovery, with an emphasis on the theoretical lower limit of RNA sizes considered for RNA sequencing and bioinformatics analyses.

NfiS, a species-specific regulatory noncoding RNA of Pseudomonas stutzeri, enhances oxidative stress tolerance in Escherichia coli

Noncoding RNAs (ncRNAs) can finely control the expression of target genes at the posttranscriptional level in prokaryotes. Regulatory small RNAs (sRNAs) designed to control target gene expression for applications in metabolic engineering and synthetic biology have been successfully developed and used. However, the effect on the heterologous expression of species- or strain-specific ncRNAs in other bacterial strains remains poorly understood. In this work, a Pseudomonas stutzeri species-specific regulatory ncRNA, NfiS, which has been shown to play an important role in the response to oxidative stress as well as osmotic stress in P. stutzeri A1501, was cloned and transferred to the Escherichia coli strain Trans10. Recombinant NfiS-expressing E. coli, namely, Trans10-nfiS, exhibited significant enhancement of tolerance to oxidative stress. To map the possible gene regulatory networks mediated by NfiS in E. coli under oxidative stress, a microarray assay was performed to delineate the transcriptomic differences between Trans10-nfiS and wild-type E. coli under H₂O₂ shock treatment conditions. In all, 1184 genes were found to be significantly altered, and these genes were divided into mainly five functional categories: stress response, regulation, metabolism related, transport or membrane protein and unknown function. Our results suggest that the P. stutzeri species-specific ncRNA NfiS acts as a regulator that integrates adaptation to H₂O₂ with other cellular stress responses and helps protect E. coli cells against oxidative damage.

Origin, Evolution, and Loss of Bacterial Small RNAs

Despite the central role of bacterial noncoding small RNAs (sRNAs) in posttranscriptional regulation, little is understood about their evolution. Here we compile what has been studied to date and trace a life cycle of sRNAs-from their mechanisms of emergence, through processes of change and frequent neofunctionalization, to their loss from bacterial lineages. Because they possess relatively unrestrictive structural requirements, we find that sRNA origins are varied, and include de novo emergence as well as formation from preexisting genetic elements via duplication events and horizontal gene transfer. The need for only partial complementarity to their mRNA targets facilitates apparent rapid change, which also contributes to significant challenges in tracing sRNAs across broad evolutionary distances. We document that recently emerged sRNAs in particular evolve quickly, mirroring dynamics observed in microRNAs, their functional analogs in eukaryotes. Mutations in mRNA-binding regions, transcriptional regulator or sigma factor binding sites, and protein-binding regions are all likely sources of shifting regulatory roles of sRNAs. Finally, using examples from the few evolutionary studies available, we examine cases of sRNA loss and describe how these may be the result of adaptive in addition to neutral processes. We highlight the need for more-comprehensive analyses of sRNA evolutionary patterns as a means to improve novel sRNA detection, enhance genome annotation, and deepen our understanding of regulatory networks in bacteria.

Functional annotation and distribution overview of RNA families in 27 Streptococcus agalactiae genomes

Background

Streptococcus agalactiae, also known as Group B Streptococcus (GBS), is a Gram-positive bacterium that colonizes the gastrointestinal and genitourinary tract of humans. This bacterium has also been isolated from various animals, such as fish and cattle. Non-coding RNAs (ncRNAs) can act as regulators of gene expression in bacteria, such as Streptococcus pneumoniae and Streptococcus pyogenes. However, little is known about the genomic distribution of ncRNAs and RNA families in S. agalactiae.

Results

Comparative genome analysis of 27 S. agalactiae strains showed more than 5 thousand genomic regions identified and classified as Core, Exclusive, and Shared genome sequences. We identified 27 to 89 RNA families per genome distributed over these regions, from these, 25 were in Core regions while Shared and Exclusive regions showed variations amongst strains. We propose that the amount and type of ncRNA present in each genome can provide a pattern to contribute in the identification of the clonal types.

Conclusions

The identification of RNA families provides an insight over ncRNAs, sRNAs and ribozymes function, that can be further explored as targets for antibiotic development or studied in gene regulation of cellular processes. RNA families could be considered as markers to determine infection capabilities of different strains. Lastly, pan-genome analysis of GBS including the full range of functional transcripts provides a broader approach in the understanding of this pathogen.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4951-z) contains supplementary material, which is available to authorized users.

Unraveling the evolution and coevolution of small regulatory RNAs and coding genes in Listeria

BACKGROUND

Small regulatory RNAs (sRNAs) are widely found in bacteria and play key roles in many important physiological and adaptation processes. Studying their evolution and screening for events of coevolution with other genomic features is a powerful way to better understand their origin and assess a common functional or adaptive relationship between them. However, evolution and coevolution of sRNAs with coding genes have been sparsely investigated in bacterial pathogens.

RESULTS

We designed a robust and generic phylogenomics approach that detects correlated evolution between sRNAs and protein-coding genes using their observed and inferred patterns of presence-absence in a set of annotated genomes. We applied this approach on 79 complete genomes of the Listeria genus and identified fifty-two accessory sRNAs, of which most were present in the Listeria common ancestor and lost during Listeria evolution. We detected significant coevolution between 23 sRNA and 52 coding genes and inferred the Listeria sRNA-coding genes coevolution network. We characterized a main hub of 12 sRNAs that coevolved with genes encoding cell wall proteins and virulence factors. Among them, an sRNA specific to L. monocytogenes species, rli133, coevolved with genes involved either in pathogenicity or in interaction with host cells, possibly acting as a direct negative post-transcriptional regulation.

CONCLUSIONS

Our approach allowed the identification of candidate sRNAs potentially involved in pathogenicity and host interaction, consistent with recent findings on known pathogenicity actors. We highlight four sRNAs coevolving with seven internalin genes, some of which being important virulence factors in Listeria.

Divergence of functional effects among bacterial sRNA paralogs

BACKGROUND

Non-coding small RNAs (sRNAs) regulate a variety of important biological processes across all life domains, including bacteria. However, little is known about the functional evolution of sRNAs in bacteria, which might occur via changes in sRNA structure and/or stability or changes in interactions between sRNAs and their associated regulatory networks, including target mRNAs. The sRNA Pxr functions as a developmental gatekeeper in the model cooperative bacterium Myxococcus xanthus. Specifically, Pxr prevents the initiation of fruiting body development when nutrients are abundant. Previous work has shown that Pxr appears to have a recent origin within a sub-clade of the myxobacteria, which allowed us to infer the most recent common ancestor of pxr and examine the divergence of Pxr since its origin.

RESULTS

To test for inter-specific divergence in functional effects, extant pxr homologs from several species and their inferred ancestor were introduced into an M. xanthus deletion mutant lacking pxr. Both the inferred ancestral pxr and all extant alleles from species containing only one copy of pxr were found to control development in M. xanthus in a qualitatively similar manner to the native M. xanthus allele. However, multiple paralogs present in Cystobacter species exhibited divergent effects, with two paralogs controlling M. xanthus development but two others failing to do so. These differences may have occurred through changes in gene expression caused by apparent structural differences in the sRNA variants encoded by these paralogs.

CONCLUSIONS

Taken together, our results suggest that Pxr plays a common fundamental role in developmental gene regulation across diverse species of myxobacteria but also that the functional effects of some Pxr variants may be evolving in some lineages.

Emergence of New sRNAs in Enteric Bacteria is Associated with Low Expression and Rapid Evolution

Non-coding small RNAs (sRNAs) are critical to post-transcriptional gene regulation in bacteria. However, unlike for protein-coding genes, the evolutionary forces that shape sRNAs are not understood. We investigated sRNAs in enteric bacteria and discovered that recently emerged sRNAs evolve at significantly faster rates than older sRNAs. Concomitantly, younger sRNAs are expressed at significantly lower levels than older sRNAs. This process could potentially facilitate the integration of newly emerged sRNAs into bacterial regulatory networks. Furthermore, it has previously been difficult to trace the evolutionary histories of sRNAs because rapid evolution obscures their original sources. We overcame this challenge by identifying a recently evolved sRNA in Escherichia coli, which allowed us to determine that novel sRNAs could emerge from vestigial bacteriophage genes, the first known source for sRNA origination.

Systematic characterization of artificial small RNA-mediated inhibition of Escherichia coli growth

A new screening system for artificial small RNAs (sRNAs) that inhibit the growth of Escherichia coli was constructed. In this system, we used a plasmid library to express RNAs of ∼120 nucleotides, each with a random 30-nucleotide sequence that can recognize its target mRNA(s). After approximately 60,000 independent colonies were screened, several plasmids that inhibited bacterial growth were isolated. To understand the inhibitory mechanism, we focused on one sRNA, S-20, that exerted a strong inhibitory effect. A time-course analysis of the proteome of S-20-expressing E. coli and a bioinformatic analysis were used to identify potential S-20 target mRNAs, and suggested that S-20 binds the translation initiation sites of several mRNAs encoding enzymes such as peroxiredoxin (osmC), glycyl-tRNA synthetase α subunit (glyQ), uncharacterized protein ygiM, and tryptophan synthase β chain (trpB). An in vitro translation analysis of chimeric luciferase-encoding mRNAs, each containing a potential S-20 target sequence, indicated that the translation of these mRNAs was inhibited in the presence of S-20. A gel shift analysis combined with the analysis of a series of S-20 mutants suggested that S-20 targets multiple mRNAs that are responsible for inhibiting E. coli growth. These data also suggest that S-20 acts like an endogenous sRNA and that E. coli can utilize artificial sRNAs.

Small Regulatory RNAs of Rickettsia conorii

Small regulatory RNAs comprise critically important modulators of gene expression in bacteria, yet very little is known about their prevalence and functions in Rickettsia species. R. conorii, the causative agent of Mediterranean spotted fever, is a tick-borne pathogen that primarily infects microvascular endothelium in humans. We have determined the transcriptional landscape of R. conorii during infection of Human Microvascular Endothelial Cells (HMECs) by strand-specific RNA sequencing to identify 4 riboswitches, 13 trans-acting (intergenic), and 22 cis-acting (antisense) small RNAs (termed ‘Rc_sR’s). Independent expression of four novel trans-acting sRNAs (Rc_sR31, Rc_sR33, Rc_sR35, and Rc_sR42) and known bacterial sRNAs (6S, RNaseP_bact_a, ffs, and α-tmRNA) was next confirmed by Northern hybridization. Comparative analysis during infection of HMECs vis-à-vis tick AAE2 cells revealed significantly higher expression of Rc_sR35 and Rc_sR42 in HMECs, whereas Rc_sR31 and Rc_sR33 were expressed at similar levels in both cell types. We further predicted a total of 502 genes involved in all important biological processes as potential targets of Rc_sRs and validated the interaction of Rc_sR42 with cydA (cytochrome d ubiquinol oxidase subunit I). Our findings constitute the first evidence of the existence of post-transcriptional riboregulatory mechanisms in R. conorii and interactions between a novel Rc_sR and its target mRNA.

Exploring lateral genetic transfer among microbial genomes using TF-IDF

Many microbes can acquire genetic material from their environment and incorporate it into their genome, a process known as lateral genetic transfer (LGT). Computational approaches have been developed to detect genomic regions of lateral origin, but typically lack sensitivity, ability to distinguish donor from recipient, and scalability to very large datasets. To address these issues we have introduced an alignment-free method based on ideas from document analysis, term frequency-inverse document frequency (TF-IDF). Here we examine the performance of TF-IDF on three empirical datasets: 27 genomes of Escherichia coli and Shigella, 110 genomes of enteric bacteria, and 143 genomes across 12 bacterial and three archaeal phyla. We investigate the effect of k-mer size, gap size and delineation of groups on the inference of genomic regions of lateral origin, finding an interplay among these parameters and sequence divergence. Because TF-IDF identifies donor groups and delineates regions of lateral origin within recipient genomes, aggregating these regions by gene enables us to explore, for the first time, the mosaic nature of lateral genes including the multiplicity of biological sources, ancestry of transfer and over-writing by subsequent transfers. We carry out Gene Ontology enrichment tests to investigate which biological processes are potentially affected by LGT.

The novel regulatory ncRNA, NfiS, optimizes nitrogen fixation via base pairing with the nitrogenase gene nifK mRNA in Pseudomonas stutzeri A1501

Unlike most Pseudomonas, the root-associated bacterium Pseudomonas stutzeri A1501 fixes nitrogen after the horizontal acquisition of a nitrogen-fixing (nif) island. A genome-wide search for small noncoding RNAs (ncRNAs) in P. stutzeri A1501 identified the novel P. stutzeri-specific ncRNA NfiS in the core genome, whose synthesis was significantly induced under nitrogen fixation or sorbitol stress conditions. The expression of NfiS was RNA chaperone Hfq-dependent and activated by the sigma factor RpoN/global nitrogen activator NtrC/nif-specific activator NifA regulatory cascade. The nfiS-deficient mutant displayed reduced nitrogenase activity, as well as increased sensitivity to multiple stresses, such as osmotic and oxidative stresses. Secondary structure prediction and complementation studies confirmed that a stem-loop structure was essential for NfiS to regulate the nitrogenase gene nifK mRNA synthesis and thus nitrogenase activity. Microscale thermophoresis and physiological analysis showed that NfiS directly pairs with nifK mRNA and ultimately enhances nitrogenase activity by increasing the translation efficiency and the half-life of nifK mRNA. Our data also suggest structural and functional divergence of NfiS evolution in diazotrophic and nondiazotrophic backgrounds. It is proposed that NfiS was recruited by nifK mRNA as a novel regulator to integrate the horizontally acquired nif island into host global networks.

The role of mRNA structure in bacterial translational regulation

The characteristics of bacterial messenger RNAs (mRNAs) that influence translation efficiency provide many convenient handles for regulation of gene expression, especially when coupled with the processes of transcription termination and mRNA degradation. An mRNA's structure, especially near the site of initiation, has profound consequences for how readily it is translated. This property allows bacterial gene expression to be altered by changes to mRNA structure induced by temperature, or interactions with a wide variety of cellular components including small molecules, other RNAs (such as sRNAs and tRNAs), and RNA-binding proteins. This review discusses the links between mRNA structure and translation efficiency, and how mRNA structure is manipulated by conditions and signals within the cell to regulate gene expression. The range of RNA regulators discussed follows a continuum from very complex tertiary structures such as riboswitch aptamers and ribosomal protein-binding sites to thermosensors and mRNA:sRNA interactions that involve only base-pairing interactions. Furthermore, the high degrees of diversity observed for both mRNA structures and the mechanisms by which inhibition of translation occur have significant consequences for understanding the evolution of bacterial translational regulation. WIREs RNA 2017, 8:e1370. doi: 10.1002/wrna.1370 For further resources related to this article, please visit the WIREs website.

A Prophage-Encoded Small RNA Controls Metabolism and Cell Division in Escherichia coli

Hundreds of small RNAs (sRNAs) have been identified in diverse bacterial species, and while the functions of most remain unknown, some regulate key processes, particularly stress responses. The sRNA DicF was identified over 25 years ago as an inhibitor of cell division but since then has remained uncharacterized. DicF consists of 53 nucleotides and is encoded by a gene carried on a prophage (Qin) in the genomes of many Escherichia coli strains. We demonstrated that DicF inhibits cell division via direct base pairing with ftsZ mRNA to repress translation and prevent new synthesis of the bacterial tubulin homolog FtsZ. Systems analysis using computational and experimental methods identified additional mRNA targets of DicF: xylR and pykA mRNAs, encoding the xylose uptake and catabolism regulator and pyruvate kinase, respectively. Genetic analyses showed that DicF directly base pairs with and represses translation of these targets. Phenotypes of cells expressing DicF variants demonstrated that DicF-associated growth inhibition is not solely due to repression of ftsZ, indicating that the physiological consequences of DicF-mediated regulation extend beyond effects on cell division caused by reduced FtsZ synthesis. IMPORTANCE sRNAs are ubiquitous and versatile regulators of bacterial gene expression. A number of well-characterized examples in E. coli are highly conserved and present in the E. coli core genome. In contrast, the sRNA DicF (identified over 20 years ago but remaining poorly characterized) is encoded by a gene carried on a defective prophage element in many E. coli genomes. Here, we characterize DicF in order to better understand how horizontally acquired sRNA regulators impact bacterial gene expression and physiology. Our data confirm the long-hypothesized DicF-mediated regulation of ftsZ, encoding the bacterial tubulin homolog required for cell division. We further uncover DicF-mediated posttranscriptional control of metabolic gene expression. Ectopic production of DicF is highly toxic to E. coli cells, but the toxicity is not attributable to DicF regulation of ftsZ. Further work is needed to reveal the biological roles of and benefits for the host conferred by DicF and other products encoded by defective prophages.

Riboregulators: Fine-Tuning Virulence in Shigella

Within the past several years, RNA-mediated regulation (ribo-regulation) has become increasingly recognized for its importance in controlling critical bacterial processes. Regulatory RNA molecules, or riboregulators, are perpetually responsive to changes within the micro-environment of a bacterium. Notably, several characterized riboregulators control virulence in pathogenic bacteria, as is the case for each riboregulator characterized to date in Shigella. The timing of virulence gene expression and the ability of the pathogen to adapt to rapidly changing environmental conditions is critical to the establishment and progression of infection by Shigella species; ribo-regulators mediate each of these important processes. This mini review will present the current state of knowledge regarding RNA-mediated regulation in Shigella by detailing the characterization and function of each identified riboregulator in these pathogens.

Diverse mechanisms of post-transcriptional repression by the small RNA regulator of glucose-phosphate stress

The Escherichia coli small RNA SgrS controls a metabolic stress response that occurs upon accumulation of certain glycolytic intermediates. SgrS base pairs with and represses translation of ptsG and manXYZ mRNAs, which encode sugar transporters, and activates translation of yigL mRNA, encoding a sugar phosphatase. This study defines four new genes as direct targets of E. coli SgrS. These new targets, asd, adiY, folE and purR, encode transcription factors or enzymes of diverse metabolic pathways, including aspartate semialdehyde dehydrogenase, arginine decarboxylase gene activator, GTP cyclohydrolase I and a repressor of purine biosynthesis, respectively. SgrS represses translation of each of the four target mRNAs via distinct mechanisms. SgrS binding sites overlapping the Shine-Dalgarno sequences of adiY and folE mRNAs suggest that SgrS pairing with these targets directly occludes ribosome binding and prevents translation initiation. SgrS binding within the purR coding sequence recruits the RNA chaperone Hfq to directly repress purR translation. Two separate SgrS binding sites were found on asd mRNA, and both are required for full translational repression. Ectopic overexpression of asd, adiY and folE is specifically detrimental to cells experiencing glucose-phosphate stress, suggesting that SgrS-dependent repression of the metabolic functions encoded by these targets promotes recovery from glucose-phosphate stress.

Advantages of mixing bioinformatics and visualization approaches for analyzing sRNA-mediated regulatory bacterial networks

The revolution in high-throughput sequencing technologies has enabled the acquisition of gigabytes of RNA sequences in many different conditions and has highlighted an unexpected number of small RNAs (sRNAs) in bacteria. Ongoing exploitation of these data enables numerous applications for investigating bacterial transacting sRNA-mediated regulation networks. Focusing on sRNAs that regulate mRNA translation in trans, recent works have noted several sRNA-based regulatory pathways that are essential for key cellular processes. Although the number of known bacterial sRNAs is increasing, the experimental validation of their interactions with mRNA targets remains challenging and involves expensive and time-consuming experimental strategies. Hence, bioinformatics is crucial for selecting and prioritizing candidates before designing any experimental work. However, current software for target prediction produces a prohibitive number of candidates because of the lack of biological knowledge regarding the rules governing sRNA-mRNA interactions. Therefore, there is a real need to develop new approaches to help biologists focus on the most promising predicted sRNA-mRNA interactions. In this perspective, this review aims at presenting the advantages of mixing bioinformatics and visualization approaches for analyzing predicted sRNA-mediated regulatory bacterial networks.

How do base-pairing small RNAs evolve?

The increasing numbers of characterized base-pairing small RNAs (sRNAs) and the identification of these regulators in a broad range of bacteria are allowing comparisons between species and explorations of sRNA evolution. In this review, we describe some examples of trans-encoded base-pairing sRNAs that are species-specific and others that are more broadly distributed. We also describe examples of sRNA orthologs where different features are conserved. These examples provide the background for a discussion of mechanisms of sRNA evolution and selective pressures on the sRNAs and their mRNA target(s).

Discovery of putative small non-coding RNAs from the obligate intracellular bacterium Wolbachia pipientis

Wolbachia pipientis is an endosymbiotic bacterium that induces a wide range of effects in its insect hosts, including manipulation of reproduction and protection against pathogens. Little is known of the molecular mechanisms underlying the insect-Wolbachia interaction, though it is likely to be mediated via the secretion of proteins or other factors. There is an increasing amount of evidence that bacteria regulate many cellular processes, including secretion of virulence factors, using small non-coding RNAs (sRNAs), but sRNAs have not previously been described from Wolbachia. We have used two independent approaches, one based on comparative genomics and the other using RNA-Seq data generated for gene expression studies, to identify candidate sRNAs in Wolbachia. We experimentally characterized the expression of one of these candidates in four Wolbachia strains, and showed that it is differentially regulated in different host tissues and sexes. Given the roles played by sRNAs in other host-associated bacteria, the conservation of the candidate sRNAs between different Wolbachia strains, and the sex- and tissue-specific differential regulation we have identified, we hypothesise that sRNAs may play a significant role in the biology of Wolbachia, and in particular in its interactions with its host.

Systematic analysis of the role of bacterial Hfq-interacting sRNAs in the response to antibiotics

OBJECTIVES

To systematically analyse the interplay between the expression of Hfq-associated small non-coding RNAs (sRNAs) and antibiotic susceptibility in Gram-negative bacteria.

METHODS

To identify the roles of sRNAs in the antibiotic susceptibility of Escherichia coli and Salmonella species, susceptibility tests, growth analyses and viability assays were performed using E. coli Hfq-associated sRNAs from overexpression libraries. Prediction, susceptibility testing of gene knockouts and expression analysis of target genes under conditions of sRNA overexpression or knockout were performed to identify candidate targets for modulating antibiotic susceptibility.

RESULTS

The susceptibilities of E. coli strains overexpressing each of the 26 known Hfq-dependent sRNAs to major classes of antibiotics were determined. Induced expression of 17 sRNAs modulated the susceptibility of E. coli to antibiotics. Among them, four sRNA knockout strains partially or completely reversed susceptibility phenotypes of sRNA overexpression. The phenotype of OxyS, RseX or MicF was not entirely dependent on the presence of Hfq protein, in contrast to the dependency of previously characterized roles. The function of eight of nine sRNAs was found to be conserved in the response to antibiotics in Salmonella. Some MicF- or RyeB-mediated cellular target genes and pathways that may be important for the regulation of antibiotic susceptibility were identified. Finally, the overexpression of RyeB potentiated the efficacy of levofloxacin against MDR strains.

CONCLUSIONS

Our data indicate that Hfq-associated sRNAs potentially enable bacteria to adapt to antibiotic challenges via multifaceted approaches. Therefore, sRNA-based applications will form a new antibiotic arsenal for combating the rise in antibiotic resistance.

Genome rearrangements can make and break small RNA genes

Small RNAs (sRNAs) are short, transcribed regulatory elements that are typically encoded in the intergenic regions (IGRs) of bacterial genomes. Several sRNAs, first recognized in Escherichia coli, are conserved among enteric bacteria, but because of the regulatory roles of sRNAs, differences in sRNA repertoires might be responsible for features that differentiate closely related species. We scanned the E. coli MG1655 and Salmonella enterica Typhimurium genomes for nonsyntenic IGRs as a potential source of uncharacterized, species-specific sRNAs and found that genome rearrangements have reconfigured several IGRs causing the disruption and formation of sRNAs. Within an IGR that is present in E. coli but was disrupted in Salmonella by a translocation event is an sRNA that is associated with the FNR/CRP global regulators and influences E. coli biofilm formation. A Salmonella-specific sRNA evolved de novo through point mutations that generated a σ⁷⁰ promoter sequence in an IGR that arose through genome rearrangement events. The differences in the sRNA pools among bacterial species have previously been ascribed to duplication, deletion, or horizontal acquisition. Here, we show that genomic rearrangements also contribute to this process by either disrupting sRNA-containing IGRs or creating IGRs in which novel sRNAs may evolve.

Sibling rivalry: related bacterial small RNAs and their redundant and non-redundant roles

Small RNA molecules (sRNAs) are now recognized as key regulators controlling bacterial gene expression, as sRNAs provide a quick and efficient means of positively or negatively altering the expression of specific genes. To date, numerous sRNAs have been identified and characterized in a myriad of bacterial species, but more recently, a theme in bacterial sRNAs has emerged: the presence of more than one highly related sRNAs produced by a given bacterium, here termed sibling sRNAs. Sibling sRNAs are those that are highly similar at the nucleotide level, and while it might be expected that sibling sRNAs exert identical regulatory functions on the expression of target genes based on their high degree of relatedness, emerging evidence is demonstrating that this is not always the case. Indeed, there are several examples of bacterial sibling sRNAs with non-redundant regulatory functions, but there are also instances of apparent regulatory redundancy between sibling sRNAs. This review provides a comprehensive overview of the current knowledge of bacterial sibling sRNAs, and also discusses important questions about the significance and evolutionary implications of this emerging class of regulators.

Large-scale mapping of sequence-function relations in small regulatory RNAs reveals plasticity and modularity

Two decades into the genomics era the question of mapping sequence to function has evolved from identifying functional elements to characterizing their quantitative properties including, in particular, their specificity and efficiency. Here, we use a large-scale approach to establish a quantitative map between the sequence of a bacterial regulatory RNA and its efficiency in modulating the expression of its targets. Our approach generalizes the sort-seq method, introduced recently to analyze promoter sequences, in order to accurately quantify the efficiency of a large library of sequence variants. We focus on two small RNAs (sRNAs) in E. coli, DsrA and RyhB, and their regulation of both repressed and activated targets. In addition to precisely identifying functional elements in the sRNAs, our data establish quantitative relationships between structural and energetic features of the sRNAs and their regulatory activity, and characterize a large set of direct and indirect interactions between nucleotides. A core of these interactions supports a model where specificity can be enhanced by a rigid molecular structure. Both sRNAs exhibit a modular design with limited cross-interactions, dividing the requirements for structural stability and target binding among modules.

Evolutionary patterns of Escherichia coli small RNAs and their regulatory interactions

Most bacterial small RNAs (sRNAs) are post-transcriptional regulators of gene expression, exerting their regulatory function by base-pairing with their target mRNAs. While it has become evident that sRNAs play central regulatory roles in the cell, little is known about their evolution and the evolution of their regulatory interactions. Here we used the prokaryotic phylogenetic tree to reconstruct the evolutionary history of Escherichia coli sRNAs and their binding sites on target mRNAs. We discovered that sRNAs currently present in E. coli mainly accumulated inside the Enterobacteriales order, succeeding the appearance of other types of noncoding RNAs and concurrently with the evolution of a variant of the Hfq protein exhibiting a longer C-terminal region. Our analysis of the evolutionary ages of sRNA-mRNA interactions revealed that while all sRNAs were evolutionarily older than most of their known binding sites on mRNA targets, for quite a few sRNAs there was at least one binding site that coappeared with or preceded them. It is conceivable that the establishment of these first interactions forced selective pressure on the sRNAs, after which additional targets were acquired by fitting a binding site to the active region of the sRNA. This conjecture is supported by the appearance of many binding sites on target mRNAs only after the sRNA gain, despite the prior presence of the target gene in ancestral genomes. Our results suggest a selective mechanism that maintained the sRNAs across the phylogenetic tree, and shed light on the evolution of E. coli post-transcriptional regulatory network.

RNA-seq analysis of the influence of anaerobiosis and FNR on Shigella flexneri

Background

Shigella flexneri is an important human pathogen that has to adapt to the anaerobic environment in the gastrointestinal tract to cause dysentery. To define the influence of anaerobiosis on the virulence of Shigella, we performed deep RNA sequencing to identify transcriptomic differences that are induced by anaerobiosis and modulated by the anaerobic Fumarate and Nitrate Reduction regulator, FNR.

Results

We found that 528 chromosomal genes were differentially expressed in response to anaerobic conditions; of these, 228 genes were also influenced by FNR. Genes that were up-regulated in anaerobic conditions are involved in carbon transport and metabolism (e.g. ptsG, manX, murQ, cysP, cra), DNA topology and regulation (e.g. ygiP, stpA, hns), host interactions (e.g. yciD, nmpC, slyB, gapA, shf, msbB) and survival within the gastrointestinal tract (e.g. shiA, ospI, adiY, cysP). Interestingly, there was a marked effect of available oxygen on genes involved in Type III secretion system (T3SS), which is required for host cell invasion and pathogenesis. These genes, located on the large Shigella virulence plasmid, were down regulated in anaerobiosis in an FNR-dependent manner. We also confirmed anaerobic induction of csrB and csrC small RNAs in an FNR-independent manner.

Conclusions

Anaerobiosis promotes survival and adaption strategies of Shigella, while modulating virulence plasmid genes involved in T3SS-mediated host cell invasion. The influence of FNR on this process is more extensive than previously appreciated, although aside from the virulence plasmid, this transcriptional regulator does not govern expression of genes on other horizontally acquired sequences on the chromosome such as pathogenicity islands.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-438) contains supplementary material, which is available to authorized users.

Small RNAs in the Vibrionaceae: an ocean still to be explored

In bacteria, the discovery of noncoding small RNAs (sRNAs) as modulators of gene expression in response to environmental signals has brought new insights into bacterial gene regulation, including control of pathogenicity. The Vibrionaceae constitute a family of marine bacteria of which many are responsible for infections affecting not only humans, such as Vibrio cholerae but also fish and marine invertebrates, representing the major cause of mortality in farmed marine species. They are able to colonize many habitats, existing as planktonic forms, in biofilms or associated with various hosts. This high adaptability is linked to their capacity to generate genetic diversity, in part through lateral gene transfer, but also by varying gene expression control. In the recent years, several major studies have illustrated the importance of small regulatory sRNAs in the Vibrionaceae for the control of pathogenicity and adaptation to environment and nutrient sources such as chitin, especially in V. cholerae and Vibrio harveyi. The existence of a complex regulatory network controlled by quorum sensing has been demonstrated in which sRNAs play central roles. This review covers major advances made in the discovery and elucidation of functions of Vibrionaceae sRNAs within the last 10 years.

Towards a processual microbial ontology

Standard microbial evolutionary ontology is organized according to a nested hierarchy of entities at various levels of biological organization. It typically detects and defines these entities in relation to the most stable aspects of evolutionary processes, by identifying lineages evolving by a process of vertical inheritance from an ancestral entity. However, recent advances in microbiology indicate that such an ontology has important limitations. The various dynamics detected within microbiological systems reveal that a focus on the most stable entities (or features of entities) over time inevitably underestimates the extent and nature of microbial diversity. These dynamics are not the outcome of the process of vertical descent alone. Other processes, often involving causal interactions between entities from distinct levels of biological organisation, or operating at different time scales, are responsible not only for the destabilisation of pre-existing entities, but also for the emergence and stabilisation of novel entities in the microbial world. In this article we consider microbial entities as more or less stabilised functional wholes, and sketch a network-based ontology that can represent a diverse set of processes including, for example, as well as phylogenetic relations, interactions that stabilise or destabilise the interacting entities, spatial relations, ecological connections, and genetic exchanges. We use this pluralistic framework for evaluating (i) the existing ontological assumptions in evolution (e.g. whether currently recognized entities are adequate for understanding the causes of change and stabilisation in the microbial world), and (ii) for identifying hidden ontological kinds, essentially invisible from within a more limited perspective. We propose to recognize additional classes of entities that provide new insights into the structure of the microbial world, namely "processually equivalent" entities, "processually versatile" entities, and "stabilized" entities.

Transcriptomic profiling of the oyster pathogen Vibrio splendidus opens a window on the evolutionary dynamics of the small RNA repertoire in the Vibrio genus

Work in recent years has led to the recognition of the importance of small regulatory RNAs (sRNAs) in bacterial regulation networks. New high-throughput sequencing technologies are paving the way to the exploration of an expanding sRNA world in nonmodel bacteria. In the Vibrio genus, compared to the enterobacteriaceae, still a limited number of sRNAs have been characterized, mostly in Vibrio cholerae, where they have been shown to be important for virulence, as well as in Vibrio harveyi. In addition, genome-wide approaches in V. cholerae have led to the discovery of hundreds of potential new sRNAs. Vibrio splendidus is an oyster pathogen that has been recently associated with massive mortality episodes in the French oyster growing industry. Here, we report the first RNA-seq study in a Vibrio outside of the V. cholerae species. We have uncovered hundreds of candidate regulatory RNAs, be it cis-regulatory elements, antisense RNAs, and trans-encoded sRNAs. Conservation studies showed the majority of them to be specific to V. splendidus. However, several novel sRNAs, previously unidentified, are also present in V. cholerae. Finally, we identified 28 trans sRNAs that are conserved in all the Vibrio genus species for which a complete genome sequence is available, possibly forming a Vibrio "sRNA core."

Comparative analysis of regulatory elements between Escherichia coli and Klebsiella pneumoniae by genome-wide transcription start site profiling

Genome-wide transcription start site (TSS) profiles of the enterobacteria Escherichia coli and Klebsiella pneumoniae were experimentally determined through modified 5′ RACE followed by deep sequencing of intact primary mRNA. This identified 3,746 and 3,143 TSSs for E. coli and K. pneumoniae, respectively. Experimentally determined TSSs were then used to define promoter regions and 5′ UTRs upstream of coding genes. Comparative analysis of these regulatory elements revealed the use of multiple TSSs, identical sequence motifs of promoter and Shine-Dalgarno sequence, reflecting conserved gene expression apparatuses between the two species. In both species, over 70% of primary transcripts were expressed from operons having orthologous genes during exponential growth. However, expressed orthologous genes in E. coli and K. pneumoniae showed a strikingly different organization of upstream regulatory regions with only 20% identical promoters with TSSs in both species. Over 40% of promoters had TSSs identified in only one species, despite conserved promoter sequences existing in the other species. 662 conserved promoters having TSSs in both species resulted in the same number of comparable 5′ UTR pairs, and that regulatory element was found to be the most variant region in sequence among promoter, 5′ UTR, and ORF. In K. pneumoniae, 48 sRNAs were predicted and 36 of them were expressed during exponential growth. Among them, 34 orthologous sRNAs between two species were analyzed in depth, and the analysis showed that many sRNAs of K. pneumoniae, including pleiotropic sRNAs such as rprA, arcZ, and sgrS, may work in the same way as in E. coli. These results reveal a new dimension of comparative genomics such that a comparison of two genomes needs to be comprehensive over all levels of genome organization.

In order to investigate similarities and differences of closely related species, most of the comparative genomics studies focus on comparing the gene contents either shared or specific for each genome. However, it is also important to investigate the differences in non-coding regulatory elements because they influence the transcriptional and post-transcriptional processes. Thus, we performed a genome-wide profiling of transcription start sites (TSSs) in two species, E. coli K-12 MG1655 and K. pneumoniae MGH78578. Experimental identification of TSSs is important for precise definition of promoter regions and 5′ untranslated regions upstream of coding genes. Comparative analysis of these regulatory elements revealed the use of multiple TSSs, identical sequence motifs of promoter and Shine-Dalgarno sequence. However, we observed that the upstream regulatory regions of the majority of operons having orthologous genes were organized with different usage of promoters and TSSs, resulting in diverse and complex gene regulation. We also found that the 5′ UTR is the least conserved regulatory element in sequence between the two species. Moreover, 34 orthologous sRNAs between E. coli and K. pneumoniae were analyzed in depth. The analysis suggested many of K. pneumoniae sRNAs might regulate the target genes as in E. coli.