Regulation of bacterial metabolism by small RNAs using diverse mechanisms

Investigating the impact of insertion sequences and transposons in the genomes of the most significant phytopathogenic bacteria

Genetic variability in phytopathogens is one of the main problems encountered for effective plant disease control. This fact may be related to the presence of transposable elements (TEs), but little is known about their role in host genomes. Here, we performed the most comprehensive analysis of insertion sequences (ISs) and transposons (Tns) in the genomes of the most important bacterial plant pathogens. A total of 35 692 ISs and 71 transposons were identified in 270 complete genomes. The level of pathogen-host specialization was found to be a significant determinant of the element distribution among the species. Some Tns were identified as carrying virulence factors, such as genes encoding effector proteins of the type III secretion system and resistance genes for the antimicrobial streptomycin. Evidence for IS-mediated ectopic recombination was identified in Xanthomonas genomes. Moreover, we found that IS elements tend to be inserted in regions near virulence and fitness genes, such ISs disrupting avirulence genes in X. oryzae genomes. In addition, transcriptome analysis under different stress conditions revealed differences in the expression of genes encoding transposases in the Ralstonia solanacearum, X. oryzae, and P. syringae species. Lastly, we also investigated the role of Tns in regulation via small noncoding regulatory RNAs and found these elements may target plant-cell transcriptional activators. Taken together, the results indicate that TEs may have a fundamental role in variability and virulence in plant pathogenic bacteria.

An Integrated Affinity Chromatography-Based Approach to Unravel the sRNA Interactome in Nitrogen-Fixing Rhizobia

The activity mechanism and function of bacterial base-pairing small non-coding RNA regulators (sRNAs) are largely shaped by their main interacting cellular partners, i.e., proteins and mRNAs. We describe here an MS2 affinity chromatography-based procedure adapted to unravel the sRNA interactome in nitrogen-fixing legume endosymbiotic bacteria. The method consists of tagging of the bait sRNA at its 5'-end with the MS2 aptamer followed by pulse overexpression and immobilization of the chimeric transcript from cell lysates by an MS2-MBP fusion protein conjugated to an amylose resin. The sRNA-binding proteins and target mRNAs are further profiled by mass spectrometry and RNAseq, respectively.

Understanding the Effect of Different Glucose Concentrations in the Oligotrophic Bacterium Bacillus subtilis BS-G1 through Transcriptomics Analysis

Glucose is an important carbon source for microbial growth, and its content in infertile soils is essential for the growth of bacteria. Since the mechanism of oligotrophic bacterium adaptation in barren soils is unclear, this research employed RNA-seq technology to examine the impact of glucose concentration on the oligotrophic bacterium B. subtilis BS-G1 in soil affected by desertification. A global transcriptome analysis (RNA-Seq) revealed that the significantly differentially expressed genes (DEGs) histidine metabolism, glutamate synthesis, the HIF-1 signaling pathway, sporulation, and the TCA cycle pathway of B. subtilis BS-G1 were significantly enriched with a 0.015 g/L glucose concentration (L group), compared to a 10 g/L glucose concentration (H group). The DEGs amino acid system, two-component system, metal ion transport, and nitrogen metabolism system of B. subtilis BS-G1 were significantly enriched in the 5 g/L glucose concentration (M group), compared with the H group. In addition, the present study identified the regulation pattern and key genes under a low-glucose environment (7 mRNAs and 16 sRNAs). This study primarily investigates the variances in the regulatory pathways of the oligotrophic B. subtilis BS-G1, which holds substantial importance in comprehending the mechanism underlying the limited sugar tolerance of oligotrophic bacteria.

Direct interaction of small non-coding RNAs CjNC140 and CjNC110 optimizes expression of key pathogenic phenotypes of Campylobacter jejuni

Small non-coding RNAs (sRNAs) are important players in modulating gene expression in bacterial pathogens, but their functions are largely undetermined in Campylobacter jejuni, an important cause of foodborne gastroenteritis in humans. In this study, we elucidated the functions of sRNA CjNC140 and its interaction with CjNC110, a previously characterized sRNA involved in the regulation of several virulence phenotypes of C. jejuni. Inactivation of CjNC140 increased motility, autoagglutination, L-methionine concentration, autoinducer-2 production, hydrogen peroxide resistance, and early chicken colonization, indicating a primarily inhibitory role of CjNC140 for these phenotypes. Apart from motility, all these effects directly contrasted the previously demonstrated positive regulation by CjNC110, suggesting that CjNC110 and CjNC140 operate in an opposite manner to modulate physiologic processes in C. jejuni. RNAseq and northern blotting further demonstrated that expression of CjNC140 increased in the absence of CjNC110, while expression of CjNC110 decreased in the absence of CjNC140, suggesting a possibility of their direct interaction. Indeed, electrophoretic mobility shift assay demonstrated a direct binding between the two sRNAs via GA- (CjNC110) and CU- (CjNC140) rich stem-loops. Additionally, RNAseq and follow-up experiments identified that CjNC140 positively regulates p19, which encodes a key iron uptake transporter in Campylobacter. Furthermore, computational analysis revealed both CjNC140 and CjNC110 are highly conserved in C. jejuni, and the predicted secondary structures support CjNC140 as a functional homolog of the iron regulatory sRNA, RyhB. These findings establish CjNC140 and CjNC110 as a key checks-and- balances mechanism in maintaining homeostasis of gene expression and optimizing phenotypes critical for C. jejuni pathobiology. IMPORTANCE Gene regulation is critical to all aspects of pathogenesis of bacterial disease, and small non-coding RNAs (sRNAs) represent a new frontier in gene regulation of bacteria. In Campylobacter jejuni, the role of sRNAs remains largely unexplored. Here, we investigate the role of two highly conserved sRNAs, CjNC110 and CjNC140, and demonstrate that CjNC140 displays a primarily inhibitory role in contrast to a primarily activating role for CjNC110 for several key virulence-associated phenotypes. Our results also revealed that the sRNA regulatory pathway is intertwined with the iron uptake system, another virulence mechanism critical for in vivo colonization. These findings open a new direction for understanding C. jejuni pathobiology and identify potential targets for intervention for this major foodborne pathogen.

Small regulatory RNAs in Vibrio cholerae

Vibrio cholerae is a major human pathogen causing the diarrheal disease, cholera. Regulation of virulence in V. cholerae is a multifaceted process involving gene expression changes at the transcriptional and post-transcriptional level. Whereas various transcription factors have been reported to modulate virulence in V. cholerae, small regulatory RNAs (sRNAs) have now been established to also participate in virulence control and the regulation of virulence-associated processes, such as biofilm formation, quorum sensing, stress response, and metabolism. In most cases, these sRNAs act by base-pairing with multiple target transcripts and this process typically requires the aid of an RNA-binding protein, such as the widely conserved Hfq protein. This review article summarizes the functional roles of sRNAs in V. cholerae, their underlying mechanisms of gene expression control, and how sRNAs partner with transcription factors to modulate complex regulatory programs. In addition, we will discuss regulatory principles discovered in V. cholerae that not only apply to other Vibrio species, but further extend into the large field of RNA-mediated gene expression control in bacteria.

An integrated transcriptomics-functional genomics approach reveals a small RNA that modulates Bacteroides thetaiotaomicron sensitivity to tetracyclines

Gene expression plasticity allows bacteria to adapt to diverse environments, tie their metabolism to available nutrients, and cope with stress. This is particularly relevant in a niche as dynamic and hostile as the human intestinal tract, yet transcriptional networks remain largely unknown in gut Bacteroides spp. Here, we map transcriptional units and profile their expression levels in Bacteroides thetaiotaomicron over a suite of 15 defined experimental conditions that are relevant in vivo , such as variation of temperature, pH, and oxygen tension, exposure to antibiotic stress, and growth on simple carbohydrates or on host mucin-derived glycans. Thereby, we infer stress- and carbon source-specific transcriptional regulons, including conditional expression of capsular polysaccharides and polysaccharide utilization loci, and expand the annotation of small regulatory RNAs (sRNAs) in this organism. Integrating this comprehensive expression atlas with transposon mutant fitness data, we identify conditionally important sRNAs. One example is MasB, whose inactivation led to increased bacterial tolerance of tetracyclines. Using MS2 affinity purification coupled with RNA sequencing, we predict targets of this sRNA and discuss their potential role in the context of the MasB-associated phenotype. Together, this transcriptomic compendium in combination with functional sRNA genomics-publicly available through a new iteration of the 'Theta-Base' web browser (www.helmholtz-hiri.de/en/datasets/bacteroides-v2)-constitutes a valuable resource for the microbiome and sRNA research communities alike.

An Estuarine Cyanophage S-CREM1 Encodes Three Distinct Antitoxin Genes and a Large Number of Non-Coding RNA Genes

Cyanophages play important roles in regulating the population dynamics, community structure, metabolism, and evolution of cyanobacteria in aquatic ecosystems. Here, we report the genomic analysis of an estuarine cyanophage, S-CREM1, which represents a new genus of T4-like cyanomyovirus and exhibits new genetic characteristics. S-CREM1 is a lytic phage which infects estuarine Synechococcus sp. CB0101. In contrast to many cyanomyoviruses that usually have a broad host range, S-CREM1 only infected the original host strain. In addition to cyanophage-featured auxiliary metabolic genes (AMGs), S-CREM1 also contains unique AMGs, including three antitoxin genes, a MoxR family ATPase gene, and a pyrimidine dimer DNA glycosylase gene. The finding of three antitoxin genes in S-CREM1 implies a possible phage control of host cells during infection. One small RNA (sRNA) gene and three cis-regulatory RNA genes in the S-CREM1 genome suggest potential molecular regulations of host metabolism by the phage. In addition, S-CREM1 contains a large number of tRNA genes which may reflect a genomic adaption to the nutrient-rich environment. Our study suggests that we are still far from understanding the viral diversity in nature, and the complicated virus-host interactions remain to be discovered. The isolation and characterization of S-CREM1 further our understanding of the gene diversity of cyanophages and phage-host interactions in the estuarine environment.

Recognition and Binding of RsmE to an AGGAC Motif of RsmZ: Insights from Molecular Dynamics Simulations

CsrA/RsmE is a post-transcriptional regulator protein widely distributed in bacteria. It impedes the expression of target mRNAs by attaching their 5' untranslated region. The translation is restored by small, noncoding RNAs that sequester CsrA/RsmE acting as sponges. In both cases, the protein recognizes and attaches to specific AGGAX and AXGGAX motifs, where X refers to any nucleotide. RsmZ of Pseudomonas protegens is one of these small RNAs. The structures of some of its complexes with RsmE were disclosed a few years ago. We have used umbrella sampling simulations to force the unbinding of RsmE from the AGGAC motif located in the single-stranded region sited between stem loops 2 and 3 of RsmZ. The calculations unveiled the identity of the main residues and nucleotides involved in the process. They also showed that the region adopts a hairpin-like conformation during the initial stages of the binding. The ability to acquire this conformation requires that the region has a length of at least nine nucleotides. Besides, we performed standard molecular dynamics simulations of the isolated fragments, analyzed their typical conformations, and characterized their movements. This analysis revealed that the free molecules oscillate along specific collective coordinates that facilitate the initial stages of the binding. The results strongly suggest that the flexibility of the single-stranded region of RsmZ crucially affects the ability of its binding motif to catch RsmE.

Prediction of Bacterial sRNAs Using Sequence-Derived Features and Machine Learning

Small ribonucleic acid (sRNA) sequences are 50–500 nucleotide long, noncoding RNA (ncRNA) sequences that play an important role in regulating transcription and translation within a bacterial cell. As such, identifying sRNA sequences within an organism’s genome is essential to understand the impact of the RNA molecules on cellular processes. Recently, numerous machine learning models have been applied to predict sRNAs within bacterial genomes. In this study, we considered the sRNA prediction as an imbalanced binary classification problem to distinguish minor positive sRNAs from major negative ones within imbalanced data and then performed a comparative study with six learning algorithms and seven assessment metrics. First, we collected numerical feature groups extracted from known sRNAs previously identified in Salmonella typhimurium LT2 (SLT2) and Escherichia coli K12 (E. coli K12) genomes. Second, as a preliminary study, we characterized the sRNA-size distribution with the conformity test for Benford’s law. Third, we applied six traditional classification algorithms to sRNA features and assessed classification performance with seven metrics, varying positive-to-negative instance ratios, and utilizing stratified 10-fold cross-validation. We revisited important individual features and feature groups and found that classification with combined features perform better than with either an individual feature or a single feature group in terms of Area Under Precision-Recall curve (AUPR). We reconfirmed that AUPR properly measures classification performance on imbalanced data with varying imbalance ratios, which is consistent with previous studies on classification metrics for imbalanced data. Overall, eXtreme Gradient Boosting (XGBoost), even without exploiting optimal hyperparameter values, performed better than the other five algorithms with specific optimal parameter settings. As a future work, we plan to extend XGBoost further to a large amount of published sRNAs in bacterial genomes and compare its classification performance with recent machine learning models’ performance.

Regulation of Carbapenemase Gene Conjugation in Escherichia coli Clinical Isolates

Background: The purpose of this study is to raise awareness of the hazards of carbapenemase epidemics and provide theoretical support for preventing the spread of carbapenemase-producing organisms. Methods: A total of 893 non-duplicate E. coil strains were recruited from three major local hospitals. The carbapenemase genotype of each imipenem-resistant strain was analyzed. Molecular typing and homology analysis of the main carbapenemase-producing strains reveal the transmission mode of resistance genes. Through the conjugation experiment, the potential spreading risk of carbapenemase genes was analyzed. Extended-spectrum beta-lactamase genes and replicon detection of the conjugant carrying plasmid were performed. The unannotated Escherichia coli bacterial small non-coding RNAs (sRNAs) interacting with sdiA were predicted through a bioinformatics tool. The sRNAs overexpression and knockout strains were constructed, and the effect of sRNA on conjugation was analyzed. Results: A total of 8 carbapenemase-producing strains were detected (0.90%, 8/893). The main carbapenemase genotype was bla_{KPC -2} (7 strains). Multilocus sequence typing indicated that 7 E. coli isolates belonged to ST-10, ST-101, ST-131, ST-405, ST-410, and ST-1193, ST-2562, respectively. Homologous cluster analysis revealed that the sequence types among the 7 E. coli were high diversity. The bla_{KPC -2} genes were successfully transferred from these isolates to EC600 by conjugation. All transconjugant cells exhibited significantly reduced susceptibility to the imipenem. IncFII was the most common conjugative plasmid type (85.7%, 6/7). Bioinformatics predicted the interaction between RydB and sdiA. Further experiments found that the interaction between RydB and sdiA improved the bacterial conjugation rate between MG1655 and EC600. The regulation effect of RydB on E. coli conjugation was not affected by the replicon type and/or harboring resistance coding genotype in conjugative plasmids. Conclusion: Our findings emphasized the epidemiological characteristics of carbapenemase-resistant E. coli. A functional phenotype of the new sRNA RydB was identified, and the regulation effect of RydB on E. coli conjugation was improved.

Multi-layer CRISPRa/i circuits for dynamic genetic programs in cell-free and bacterial systems

CRISPR-Cas transcriptional circuits hold great promise as platforms for engineering metabolic networks and information processing circuits. Historically, prokaryotic CRISPR control systems have been limited to CRISPRi. Creating approaches to integrate CRISPRa for transcriptional activation with existing CRISPRi-based systems would greatly expand CRISPR circuit design space. Here, we develop design principles for engineering prokaryotic CRISPRa/i genetic circuits with network topologies specified by guide RNAs. We demonstrate that multi-layer CRISPRa/i cascades and feedforward loops can operate through the regulated expression of guide RNAs in cell-free expression systems and E. coli. We show that CRISPRa/i circuits can program complex functions by designing type 1 incoherent feedforward loops acting as fold-change detectors and tunable pulse-generators. By investigating how component characteristics relate to network properties such as depth, width, and speed, this work establishes a framework for building scalable CRISPRa/i circuits as regulatory programs in cell-free expression systems and bacterial hosts. A record of this paper's transparent peer review process is included in the supplemental information.

sRNA23, a novel small RNA, regulates to the pathogenesis of Streptococcus suis serotype 2

ABBREVIATION

sRNA: small noncoding RNA; FBA: fructose diphosphate aldolase; rplB: 50S ribosomal protein L2; RACE: rapid amplification of cDNA ends; EMSA: electrophoretic mobility shift assay; THB: Todd-Hewitt broth; FBS: fetal bovine serum; BIP: 2,2'-Bipyridine.

The "guiding" principles of noncoding RNA function

The human genome is pervasively transcribed and yet only a small fraction of these RNAs (less than 2%) are known to code for proteins. The vast majority of the RNAs are classified as noncoding RNAs (ncRNAs) and are further subgrouped as small (shorter than 200 bases) and long noncoding RNAs. The ncRNAs have been identified in all three domains of life and regulate diverse cellular processes through transcriptional and posttranscriptional gene regulation. Most of these RNAs work in conjunction with proteins forming a wide array of base pairing interactions. The determinants of these base pairing interactions are now becoming more evident and show striking similarities among the diverse group of ncRNAs. Here we present a mechanistic overview of pairing between RNA-RNA or RNA-DNA that dictates the function of ncRNAs; we provide examples to illustrate that ncRNAs work through shared evolutionary mechanisms that encompasses a guide-target interaction, involving not only classical Watson-Crick but also noncanonical Wobble and Hoogsteen base pairing. We also highlight the similarities in target selection, proofreading, and the ruler mechanism of ncRNA-protein complexes that confers target specificity and target site selection. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA-Based Catalysis > RNA-Mediated Cleavage RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

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

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

Prokaryotic Genome Annotation

In the last decade, the high-throughput and relatively low cost of short-read sequencing technologies have revolutionized prokaryotic genomics. This has led to an exponential increase in the number of bacterial and archaeal genome sequences available, as well as corresponding increase of genome assembly and annotation tools developed. Together, these hardware and software technologies have given scientists unprecedented options to study their chosen microbial systems without the need for large teams of bioinformaticists or supercomputing facilities. While these analysis tools largely fall into only a few categories, each may have different requirements, caveats and file formats, and some may be rarely updated or even abandoned. And so, despite the apparent ease in sequencing and analyzing a prokaryotic genome, it is no wonder that the budding genomicist may quickly find oneself overwhelmed. Here, we aim to provide the reader with an overview of genome annotation and its most important considerations, as well as an easy-to-follow protocol to get started with annotating a prokaryotic genome.

A dual-function RNA balances carbon uptake and central metabolism in Vibrio cholerae

Bacterial small RNAs (sRNAs) are well known to modulate gene expression by base pairing with trans‐encoded transcripts and are typically non‐coding. However, several sRNAs have been reported to also contain an open reading frame and thus are considered dual‐function RNAs. In this study, we discovered a dual‐function RNA from Vibrio cholerae, called VcdRP, harboring a 29 amino acid small protein (VcdP), as well as a base‐pairing sequence. Using a forward genetic screen, we identified VcdRP as a repressor of cholera toxin production and link this phenotype to the inhibition of carbon transport by the base‐pairing segment of the regulator. By contrast, we demonstrate that the VcdP small protein acts downstream of carbon transport by binding to citrate synthase (GltA), the first enzyme of the citric acid cycle. Interaction of VcdP with GltA results in increased enzyme activity and together VcdR and VcdP reroute carbon metabolism. We further show that transcription of vcdRP is repressed by CRP allowing us to provide a model in which VcdRP employs two different molecular mechanisms to synchronize central metabolism in V. cholerae.

Identification and initial characterization of a new pair of sibling sRNAs of Neisseria gonorrhoeae involved in type IV pilus biogenesis

Non-coding regulatory RNAs mediate post-transcriptional gene expression control by a variety of mechanisms relying mostly on base-pairing interactions with a target mRNA. Though a plethora of putative non-coding regulatory RNAs have been identified by global transcriptome analysis, knowledge about riboregulation in the pathogenic Neisseriae is still limited. Here we report the initial characterization of a pair of sRNAs of N. gonorrhoeae, TfpR1 and TfpR2, which exhibit a similar secondary structure and identical single-stranded seed regions, and therefore might be considered as sibling sRNAs. By combination of in silico target prediction and sRNA pulse expression followed by differential RNA sequencing we identified target genes of TfpR1 which are involved in type IV pilus biogenesis and DNA damage repair. We provide evidence that members of the TfpR1 regulon can also be targeted by the sibling TfpR2.

Molecular Dynamics Simulations Unveil the Basis of the Sequential Binding of RsmE to the Noncoding RNA RsmZ

CsrA/RsmE are dimeric proteins that bind to targeted mRNAs repressing translation. This mechanism modulates several metabolic pathways and allows bacteria to efficiently adjust their responses to environmental changes. In turn, small RNAs (sRNA) such as CsrB or RsmZ, restore translation by sequestering CsrA/RsmE dimers. Thus, these molecules act in tandem as a gene-expression regulatory system. Recently, a combined NMR-EPR approach solved the structure of part of RsmZ of Pseudomonas fluorescens, attached to three RsmE dimers. The study demonstrated that RsmE assembles onto RsmZ following a specific sequential order. The reasons underlying this peculiar behavior are still unclear. Here, we present a molecular dynamics analysis that explores the conformational diversity of RsmZ and RsmZ-RsmE complexes. The results reveal a clear pattern regarding the exposure of the alternative GGA binding motifs of RsmZ. This pattern is tuned by the attachment of RsmE dimers. Altogether, the observations provide a simple and convincing explanation for the order observed in the sequestration of RsmE dimers. Typical structures for RsmZ and RsmZ-RsmE complexes have been identified. Their characteristics concerning the exposure of the GGA sequences are presented and their most significant interactions are described.

Pervasive RNA Regulation of Metabolism Enhances the Root Colonization Ability of Nitrogen-Fixing Symbiotic α-Rhizobia

The rhizosphere and rhizoplane are nutrient-rich but selective environments for the root microbiome. Here, we deciphered a posttranscriptional network regulated by the homologous trans-small RNAs (sRNAs) AbcR1 and AbcR2, which rewire the metabolism of the nitrogen-fixing α-rhizobium Sinorhizobium meliloti during preinfection stages of symbiosis with its legume host alfalfa. The LysR-type regulator LsrB, which transduces the cell redox state, is indispensable for AbcR1 expression in actively dividing bacteria, whereas the stress-induced transcription of AbcR2 depends on the alternative σ factor RpoH1. MS2 affinity purification coupled with RNA sequencing unveiled exceptionally large and overlapping AbcR1/2 mRNA interactomes, jointly representing ⁓6% of the S. meliloti protein-coding genes. Most mRNAs encode transport/metabolic proteins whose translation is silenced by base pairing to two distinct anti-Shine Dalgarno motifs that function independently in both sRNAs. A metabolic model-aided analysis of the targetomes predicted changes in AbcR1/2 expression driven by shifts in carbon/nitrogen sources, which were confirmed experimentally. Low AbcR1/2 levels in some defined media anticipated overexpression growth phenotypes linked to the silencing of specific mRNAs. As a proof of principle, we confirmed AbcR1/2-mediated downregulation of the l-amino acid AapQ permease. AbcR1/2 interactomes are well represented in rhizosphere-related S. meliloti transcriptomic signatures. Remarkably, a lack of AbcR1 specifically compromised the ability of S. meliloti to colonize the root rhizoplane. The AbcR1 regulon likely ranks the utilization of available substrates to optimize metabolism, thus conferring on S. meliloti an advantage for efficient rhizosphere/rhizoplane colonization. AbcR1 regulation is predicted to be conserved in related α-rhizobia, which opens unprecedented possibilities for engineering highly competitive biofertilizers. IMPORTANCE Nitrogen-fixing root nodule symbioses between rhizobia and legume plants provide more than half of the combined nitrogen incorporated annually into terrestrial ecosystems, rendering plant growth independent of environmentally unfriendly chemical fertilizers. The success of symbiosis depends primarily on the capacity of rhizobia to establish competitive populations in soil and rhizosphere environments. Here, we provide insights into the regulation and architecture of an extensive RNA posttranscriptional network that fine-tunes the metabolism of the alfalfa symbiont S. meliloti, thereby enhancing the ability of this beneficial bacterium to colonize nutrient-rich but extremely selective niches, such as the rhizosphere of its host plant. This pervasive RNA regulation of metabolism is a major adaptive mechanism, predicted to operate in diverse rhizobial species. Because RNA regulation relies on modifiable base-pairing interactions, our findings open unexplored avenues for engineering the legumes rhizobiome within sustainable agricultural practices.

Genome-wide screen and functional analysis in Xanthomonas reveal a large number of mRNA-derived sRNAs, including the novel RsmA-sequester RsmU

Although bacterial small noncoding RNAs (sRNAs) are known to play a critical role in various cellular processes, including pathogenesis, the identity and action of such sRNAs are still poorly understood in many organisms. Here we have performed a genome-wide screen and functional analysis of the sRNAs in Xanthomonas campestris pv. campestris (Xcc), an important phytopathogen. The 50-500-nt RNA fragments isolated from the wild-type strain grown in a virulence gene-inducing condition were sequenced and a total of 612 sRNA candidates (SRCs) were identified. The majority (82%) of the SRCs were derived from mRNA, rather than specific sRNA genes. A representative panel of 121 SRCs were analysed by northern blotting; 117 SRCs were detected, supporting the contention that the overwhelming majority of the 612 SRCs identified are indeed sRNAs. Phenotypic analysis of strains overexpressing different candidates showed that a particular sRNA, RsmU, acts as a negative regulator of virulence, the hypersensitive response, and cell motility in Xcc. In vitro electrophoretic mobility shift assay and in vivo coimmunoprecipitation analyses indicated that RsmU interacted with the global posttranscriptional regulator RsmA, although sequence analysis displayed that RsmU is not a member of the sRNAs families known to antagonize RsmA. Northern blotting analyses demonstrated that RsmU has two isoforms that are processed from the 3'-untranslated region of the mRNA of XC1332 predicted to encode ComEA, a periplasmic protein required for DNA uptake in bacteria. This work uncovers an unexpected major sRNA biogenesis strategy in bacteria and a hidden layer of sRNA-mediated virulence regulation in Xcc.

A Novel sRNA in Shigella flexneri That Regulates Tolerance and Virulence Under Hyperosmotic Pressure

Regulation of the environmental stress response and virulence of Shigella flexneri may involve multiple signaling pathways; however, these mechanisms are not well-defined. In bacteria, small regulatory RNAs (sRNAs) regulate bacterial growth, metabolism, virulence, and environmental stress response. Therefore, identifying novel functional sRNAs in S. flexneri could help elucidate pathogenic adaptations to host micro-environmental stresses and associated virulence. The aim of this study was to confirm the presence of an sRNA, Ssr54, in S. flexneri and to determine its functions and possible mechanism of action. Ssr54 was found to regulate tolerance and virulence under hyperosmotic pressure. Its expression was verified by qRT-PCR and Northern blotting, and its genomic position was confirmed by 5'-rapid amplification of cDNA ends. Ssr54 expression was significantly decreased (~ 80%) under hyperosmotic conditions (680 mM NaCl), and the survival rate of the Ssr54 deletion strain increased by 20% under these conditions. This suggested that Ssr54 has been selected to promote host survival under hyperosmotic conditions. Additionally, virulence assessment, including guinea pig Sereny test and competitive invasion assays in mouse lungs, revealed that Ssr54 deletion significantly decreased S. flexneri virulence. Two-dimensional gel analyses suggest that Ssr54 may modulate the expression of tolC, ompA, and treF genes, which may affect the virulence and survival of S. flexneri under osmotic pressures. Furthermore, treF expression has been shown to improve the survival of S. flexneri under osmotic pressures. These results suggest that Ssr54 has a broad range of action in S. flexneri response to hyperosmotic environmental stresses and in controlling its virulence to adapt to environmental stresses encountered during host infection.

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

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

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

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.

In silico prediction and expression profile analysis of small non-coding RNAs in Herbaspirillum seropedicae SmR1

BACKGROUND

Herbaspirillum seropedicae is a diazotrophic bacterium from the β-proteobacteria class that colonizes endophytically important gramineous species, promotes their growth through phytohormone-dependent stimulation and can express nif genes and fix nitrogen inside plant tissues. Due to these properties this bacterium has great potential as a commercial inoculant for agriculture. The H. seropedicae SmR1 genome is completely sequenced and annotated but despite the availability of diverse structural and functional analysis of this genome, studies involving small non-coding RNAs (sRNAs) has not yet been done. We have conducted computational prediction and RNA-seq analysis to select and confirm the expression of sRNA genes in the H. seropedicae SmR1 genome, in the presence of two nitrogen independent sources and in presence of naringenin, a flavonoid secreted by some plants.

RESULTS

This approach resulted in a set of 117 sRNAs distributed in riboswitch, cis-encoded and trans-encoded categories and among them 20 have Rfam homologs. The housekeeping sRNAs tmRNA, ssrS and 4.5S were found and we observed that a large number of sRNAs are more expressed in the nitrate condition rather than the control condition and in the presence of naringenin. Some sRNAs expression were confirmed in vitro and this work contributes to better understand the post transcriptional regulation in this bacterium.

CONCLUSIONS

H. seropedicae SmR1 express sRNAs in the presence of two nitrogen sources and/or in the presence of naringenin. The functions of most of these sRNAs remains unknown but their existence in this bacterium confirms the evidence that sRNAs are involved in many different cellular activities to adapt to nutritional and environmental changes.

Trends in Symbiont-Induced Host Cellular Differentiation

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

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

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

Determinants of target prioritization and regulatory hierarchy for the bacterial small RNA SgrS

Small RNA (sRNA) regulators promote efficient responses to stress, but the mechanisms for prioritizing target mRNA regulation remain poorly understood. This study examines mechanisms underlying hierarchical regulation by the sRNA SgrS, found in enteric bacteria and produced under conditions of metabolic stress. SgrS posttranscriptionally coordinates a nine-gene regulon to restore growth and homeostasis. An in vivo reporter system quantified SgrS-dependent regulation of target genes and established that SgrS exhibits a clear target preference. Regulation of some targets is efficient even at low SgrS levels, whereas higher SgrS concentrations are required to regulate other targets. In vivo and in vitro analyses revealed that RNA structure and the number and position of base pairing sites relative to the start of translation impact the efficiency of regulation of SgrS targets. The RNA chaperone Hfq uses distinct modes of binding to different SgrS mRNA targets, which differentially influences positive and negative regulation. The RNA degradosome plays a larger role in regulation of some SgrS targets compared to others. Collectively, our results suggest that sRNA selection of target mRNAs and regulatory hierarchy are influenced by several molecular features and that the combination of these features precisely tunes the efficiency of regulation of multi-target sRNA regulons.

Mechanisms for Hfq-Independent Activation of rpoS by DsrA, a Small RNA, in Escherichia coli

Many small RNAs (sRNAs) regulate gene expression by base pairing to their target messenger RNAs (mRNAs) with the help of Hfq in Escherichia coli. The sRNA DsrA activates translation of the rpoS mRNA in an Hfq-dependent manner, but this activation ability was found to partially bypass Hfq when DsrA is overproduced. The precise mechanism by which DsrA bypasses Hfq is unknown. In this study, we constructed strains lacking all three rpoS-activating sRNAs (i.e., ArcZ, DsrA, and RprA) in hfq+ and Hfq- backgrounds, and then artificially regulated the cellular DsrA concentration in these strains by controlling its ectopic expression. We then examined how the expression level of rpoS was altered by a change in the concentration of DsrA. We found that the translation and stability of the rpoS mRNA are both enhanced by physiological concentrations of DsrA regardless of Hfq, but that depletion of Hfq causes a rapid degradation of DsrA and thereby decreases rpoS mRNA stability. These results suggest that the observed Hfq dependency of DsrA-mediated rpoS activation mainly results from the destabilization of DsrA in the absence of Hfq, and that DsrA itself contributes to the translational activation and stability of the rpoS mRNA in an Hfq-independent manner.

On a stake-out: Mycobacterial small RNA identification and regulation

Persistence of mycobacteria in the hostile environment of human macrophage is pivotal for its successful pathogenesis. Rapid adaptation to diverse stresses is the key aspect for their survival in the host cells. A range of heterogeneous mechanisms operate in bacteria to retaliate stress conditions. Small RNAs (sRNA) have been implicated in many of those mechanisms in either a single or multiple regulatory networks to post-transcriptionally modulate bacterial gene expression. Although small RNA profiling in mycobacteria by advanced technologies like deep sequencing, tilling microarray etc. have identified hundreds of sRNA, however, a handful of those small RNAs have been unearthed with precise regulatory mechanism. Extensive investigations on sRNA-mediated gene regulations in eubacteria like Escherichia coli revealed the existence of a plethora of distinctive sRNA mechanisms e.g. base pairing, protein sequestration, RNA decoy etc. Increasing studies on mycobacterial sRNA also discovered several eccentric mechanisms where sRNAs act at the posttranscriptional stage to either activate or repress target gene expression that lead to promote mycobacterial survival in stresses. Several intrinsic features like high GC content, absence of any homologue of abundant RNA chaperones, Hfq and ProQ, isolate sRNA mechanisms of mycobacteria from that of other bacteria. An insightful approach has been taken in this review to describe sRNA identification and its regulations in mycobacterial species especially in Mycobacterium tuberculosis.

mRNA dynamics and alternative conformations adopted under low and high arginine concentrations control polyamine biosynthesis in Salmonella

Putrescine belongs to the large group of polyamines, an essential class of metabolites that exists throughout all kingdoms of life. The Salmonella speF gene encodes an inducible ornithine decarboxylase that produces putrescine from ornithine. Putrescine can be also synthesized from arginine in a parallel metabolic pathway. Here, we show that speF expression is controlled at multiple levels through regulatory elements contained in a long leader sequence. At the heart of this regulation is a short open reading frame, orf34, which is required for speF production. Translation of orf34 interferes with Rho-dependent transcription termination and helps to unfold an inhibitory RNA structure sequestering speF ribosome-binding site. Two consecutive arginine codons in the conserved domain of orf34 provide a third level of speF regulation. Uninterrupted translation of orf34 under conditions of high arginine allows the formation of a speF mRNA structure that is degraded by RNase G, whereas ribosome pausing at the consecutive arginine codons in the absence of arginine enables the formation of an alternative structure that is resistant to RNase G. Thus, the rate of ribosome progression during translation of the upstream ORF influences the dynamics of speF mRNA folding and putrescine production. The identification of orf34 and its regulatory functions provides evidence for the evolutionary conservation of ornithine decarboxylase regulatory elements and putrescine production.

Polyamines are widely distributed in nature, they bind nucleic acids and proteins and although their exact mechanism of action is not clear, their effect on fundamental cellular functions is well documented. The canonical biosynthesis pathway of polyamines is conserved and begins with speF encoding ornithine decarboxylase, an inducible enzyme that produces putrescine from ornithine. Putrescine can also be produced from arginine in an alternative metabolic pathway. Here, we show that the rate of ribosome progression during translation of a short ORF (ORF34) upstream of speF influences the dynamics of speF mRNA folding and thus putrescine production. Uninterrupted translation of orf34 carrying two consecutive arginine codons, under conditions of high arginine, results in the formation of a speF mRNA structure that is degraded by RNase G, whereas ribosomes slow-down at the consecutive arginine codons in the absence of arginine enables the formation of an alternative structure that is unsusceptible to RNase G and thus results in putrescine production. The study of Salmonella speF regulation provides evidence that, despite variations in the mechanistic details, RNA-based regulation of putrescine biosynthesis and ornithine decarboxylase is conserved from bacteria to mammals.

Quorum sensing in Pseudomonas aeruginosa mediated by RhlR is regulated by a small RNA PhrD

Pseudomonas aeruginosa is a highly invasive human pathogen in spite of the absence of classical host specific virulence factors. Virulence factors regulated by quorum sensing (QS) in P. aeruginosa cause acute infections to shift to chronic diseases. Several small regulatory RNAs (sRNAs) mediate fine-tuning of bacterial responses to environmental signals and regulate quorum sensing. In this study, we show that the quorum sensing regulator RhlR is positively influenced upon over expression of the Hfq dependent small RNA PhrD in Pseudomonas. RhlR transcripts starting from two of the four different promoters have same sequence predicted to base pair with PhrD. Over expression of PhrD increased RhlR transcript levels and production of the biosurfactant rhamnolipid and the redox active pyocyanin pigment. A rhlR::lacZ translational fusion from one of the four promoters showed 2.5-fold higher expression and, a 9-fold increase in overall rhlR transcription was seen in the wild type when compared to the isogenic phrD disruption mutant. Expression, in an E. coli host background, of a rhlR::lacZ fusion in comparison to a construct that harboured a scrambled interaction region resulted in a 10-fold increase under phrD over expression. The interaction of RhlR-5′UTR with PhrD in E. coli indicated that this regulation could function without the involvement of any Pseudomonas specific proteins. Overall, this study demonstrates that PhrD has a positive effect on RhlR and its associated physiology in P. aeruginosa.

Cross-Regulation between Bacteria and Phages at a Posttranscriptional Level

The study of bacteriophages (phages) and prophages has provided key insights into almost every cellular process as well as led to the discovery of unexpected new mechanisms and the development of valuable tools. This is exemplified for RNA-based regulation. For instance, the characterization and exploitation of the antiphage CRISPR (clustered regularly interspaced short palindromic repeat) systems is revolutionizing molecular biology. Phage-encoded proteins such as the RNA-binding MS2 protein, which is broadly used to isolate tagged RNAs, also have been developed as valuable tools. Hfq, the RNA chaperone protein central to the function of many base-pairing small RNAs (sRNAs), was first characterized as a bacterial host factor required for Qβ phage replication. The ongoing studies of RNAs are continuing to reveal regulatory connections between infecting phages, prophages, and bacteria and to provide novel insights. There are bacterial and prophage sRNAs that regulate prophage genes, which impact bacterial virulence as well as bacterial cell killing. Conversely, phage- and prophage-encoded sRNAs modulate the expression of bacterial genes modifying metabolism. An interesting subcategory of the prophage-encoded sRNAs are sponge RNAs that inhibit the activities of bacterial-encoded sRNAs. Phages also affect posttranscriptional regulation in bacteria through proteins that inhibit or alter the activities of key bacterial proteins involved in posttranscriptional regulation. However, what is most exciting about phage and prophage research, given the millions of phage-encoded genes that have not yet been characterized, is the vast potential for discovering new RNA regulators and novel mechanisms and for gaining insight into the evolution of regulatory RNAs.

iTRAQ-based quantitative proteomics reveals insights into metabolic and molecular responses of glucose-grown cells of Rubrivivax benzoatilyticus JA2

Anoxygenic photosynthetic bacteria thrive under diverse habitats utilising an extended range of inorganic/organic compounds under different growth modes. Although they display incredible metabolic flexibility, their responses and adaptations to changing carbon regimes is largely unexplored. In the present study, we employed iTRAQ-based global proteomic profiling and physiological studies to uncover the adaptive strategies of a phototrophic bacterium, Rubrivivax benzoatilyticus JA2 to glucose. Strain JA2 displayed altered growth rates, reduced cell size and progressive loss of pigmentation when grown on glucose compared to malate under photoheterotrophic condition. A ten-fold increase in the saturated to unsaturated fatty acid ratio of glucose-grown cells indicates a possible membrane adaptation. Proteomic profiling revealed extensive metabolic remodelling in the glucose-grown cells wherein signal-transduction, selective-transcription, DNA-repair, transport and protein quality control processes were up-regulated to cope with the changing milieu. Proteins involved in DNA replication, translation, electron-transport, photosynthetic machinery were down-regulated possibly to conserve the energy. Glycolysis/gluconeogenesis, TCA cycle and pigment biosynthesis were also down-regulated. The cell has activated alternative energy metabolic pathways viz., fatty acid β-oxidation, glyoxylate, acetate-switch and Entner-Doudoroff pathways. Overall, the present study deciphered the molecular/metabolic events associated with glucose-grown cells of strain JA2 and also unraveled how a carbon source modulates the metabolic phenotypes. SIGNIFICANCE: Anoxygenic photosynthetic bacteria (APB) exhibit incredible metabolic flexibility leading to diverse phenotypes. They thrive under diverse habitat using an array of inorganic/organic compounds as carbon sources, yet their metabolic adaptation to varying carbon regime is mostly unexplored. Present study uncovered the proteomic insights of the cellular responses of strain JA2 to changing carbon sources viz. malate and glucose under photoheterotrophic conditions. Our study suggests that carbon source can also determine the metabolic fate of the cells and reshape the energy dynamics of APB. Here, for the first time study highlighted the plausible carbon source (glucose) mediated regulation of photosynthesis in APB. The study sheds light on the plausible cellular events and adaptive metabolic strategies employed by strain JA2 in presence of non-preferred carbon source. It also revealed new insights into the metabolic plasticity of APB to the changing milieu.

Regulation of glycolytic flux and overflow metabolism depending on the source of energy generation for energy demand

Overflow metabolism is a common phenomenon observed at higher glycolytic flux in many bacteria, yeast (known as Crabtree effect), and mammalian cells including cancer cells (known as Warburg effect). This phenomenon has recently been characterized as the trade-offs between protein costs and enzyme efficiencies based on coarse-graining approaches. Moreover, it has been recognized that the glycolytic flux increases as the source of energy generation changes from energetically efficient respiration to inefficient respiro-fermentative or fermentative metabolism causing overflow metabolism. It is highly desired to clarify the metabolic regulation mechanisms behind such phenomena. Metabolic fluxes are located on top of the hierarchical regulation systems, and represent the outcome of the integrated response of all levels of cellular regulation systems. In the present article, we discuss about the different levels of regulation systems for the modulation of fluxes depending on the growth rate, growth condition such as oxygen limitation that alters the metabolism towards fermentation, and genetic perturbation affecting the source of energy generation from respiration to respiro-fermentative metabolism in relation to overflow metabolism. The intracellular metabolite of the upper glycolysis such as fructose 1,6-bisphosphate (FBP) plays an important role not only for flux sensing, but also for the regulation of the respiratory activity either directly or indirectly (via transcription factors) at higher growth rate. The glycolytic flux regulation is backed up (enhanced) by unphosphorylated EIIA and HPr of the phosphotransferase system (PTS) components, together with the sugar-phosphate stress regulation, where the transcriptional regulation is further modulated by post-transcriptional regulation via the degradation of mRNA (stability of mRNA) in Escherichia coli. Moreover, the channeling may also play some role in modulating the glycolytic cascade reactions.

The novel EHEC gene asa overlaps the TEGT transporter gene in antisense and is regulated by NaCl and growth phase

Only a few overlapping gene pairs are known in the best-analyzed bacterial model organism Escherichia coli. Automatic annotation programs usually annotate only one out of six reading frames at a locus, allowing only small overlaps between protein-coding sequences. However, both RNAseq and RIBOseq show signals corresponding to non-trivially overlapping reading frames in antisense to annotated genes, which may constitute protein-coding genes. The transcription and translation of the novel 264 nt gene asa, which overlaps in antisense to a putative TEGT (Testis-Enhanced Gene Transfer) transporter gene is detected in pathogenic E. coli, but not in two apathogenic E. coli strains. The gene in E. coli O157:H7 (EHEC) was further analyzed. An overexpression phenotype was identified in two stress conditions, i.e. excess in salt or arginine. For this, EHEC overexpressing asa was grown competitively against EHEC with a translationally arrested asa mutant gene. RT-qPCR revealed conditional expression dependent on growth phase, sodium chloride, and arginine. Two potential promoters were computationally identified and experimentally verified by reporter gene expression and determination of the transcription start site. The protein Asa was verified by Western blot. Close homologues of asa have not been found in protein databases, but bioinformatic analyses showed that it may be membrane associated, having a largely disordered structure.

Poly(A) polymerase is required for RyhB sRNA stability and function in Escherichia coli

Small regulatory RNAs (sRNAs) are an important class of bacterial post-transcriptional regulators that control numerous physiological processes, including stress responses. In Gram-negative bacteria including Escherichia coli, the RNA chaperone Hfq binds many sRNAs and facilitates pairing to target transcripts, resulting in changes in mRNA transcription, translation, or stability. Here, we report that poly(A) polymerase (PAP I), which promotes RNA degradation by exoribonucleases through the addition of poly(A) tails, has a crucial role in the regulation of gene expression by Hfq-dependent sRNAs. Specifically, we show that deletion of pcnB, encoding PAP I, paradoxically resulted in an increased turnover of certain Hfq-dependent sRNAs, including RyhB. RyhB instability in the pcnB deletion strain was suppressed by mutations in hfq or ryhB that disrupt pairing of RyhB with target RNAs, by mutations in the 3' external transcribed spacer of the glyW-cysT-leuZ transcript (3'ETS^LeuZ) involved in pairing with RyhB, or an internal deletion in rne, which encodes the endoribonuclease RNase E. Finally, the reduced stability of RyhB in the pcnB deletion strain resulted in impaired regulation of some of its target mRNAs, specifically sodB and sdhCDAB. Altogether our data support a model where PAP I plays a critical role in ensuring the efficient decay of the 3'ETS^LeuZ In the absence of PAP I, the 3'ETS^LeuZ transcripts accumulate, bind Hfq, and pair with RyhB, resulting in its depletion via RNase E-mediated decay. This ultimately leads to a defect in RyhB function in a PAP I deficient strain.

The Novel Anaerobiosis-Responsive Overlapping Gene ano Is Overlapping Antisense to the Annotated Gene ECs2385 of Escherichia coli O157:H7 Sakai

Current notion presumes that only one protein is encoded at a given bacterial genetic locus. However, transcription and translation of an overlapping open reading frame (ORF) of 186 bp length were discovered by RNAseq and RIBOseq experiments. This ORF is almost completely embedded in the annotated L,D-transpeptidase gene ECs2385 of Escherichia coli O157:H7 Sakai in the antisense reading frame -3. The ORF is transcribed as part of a bicistronic mRNA, which includes the annotated upstream gene ECs2384, encoding a murein lipoprotein. The transcriptional start site of the operon resides 38 bp upstream of the ECs2384 start codon and is driven by a predicted σ⁷⁰ promoter, which is constitutively active under different growth conditions. The bicistronic operon contains a ρ-independent terminator just upstream of the novel gene, significantly decreasing its transcription. The novel gene can be stably expressed as an EGFP-fusion protein and a translationally arrested mutant of ano, unable to produce the protein, shows a growth advantage in competitive growth experiments compared to the wild type under anaerobiosis. Therefore, the novel antisense overlapping gene is named ano (anaerobiosis responsive overlapping gene). A phylostratigraphic analysis indicates that ano originated very recently de novo by overprinting after the Escherichia/Shigella clade separated from other enterobacteria. Therefore, ano is one of the very rare cases of overlapping genes known in the genus Escherichia.

A novel short L-arginine responsive protein-coding gene (laoB) antiparallel overlapping to a CadC-like transcriptional regulator in Escherichia coli O157:H7 Sakai originated by overprinting

Background

Due to the DNA triplet code, it is possible that the sequences of two or more protein-coding genes overlap to a large degree. However, such non-trivial overlaps are usually excluded by genome annotation pipelines and, thus, only a few overlapping gene pairs have been described in bacteria. In contrast, transcriptome and translatome sequencing reveals many signals originated from the antisense strand of annotated genes, of which we analyzed an example gene pair in more detail.

Results

A small open reading frame of Escherichia coli O157:H7 strain Sakai (EHEC), designated laoB (L-arginine responsive overlapping gene), is embedded in reading frame −2 in the antisense strand of ECs5115, encoding a CadC-like transcriptional regulator. This overlapping gene shows evidence of transcription and translation in Luria-Bertani (LB) and brain-heart infusion (BHI) medium based on RNA sequencing (RNAseq) and ribosomal-footprint sequencing (RIBOseq). The transcriptional start site is 289 base pairs (bp) upstream of the start codon and transcription termination is 155 bp downstream of the stop codon. Overexpression of LaoB fused to an enhanced green fluorescent protein (EGFP) reporter was possible. The sequence upstream of the transcriptional start site displayed strong promoter activity under different conditions, whereas promoter activity was significantly decreased in the presence of L-arginine. A strand-specific translationally arrested mutant of laoB provided a significant growth advantage in competitive growth experiments in the presence of L-arginine compared to the wild type, which returned to wild type level after complementation of laoB in trans. A phylostratigraphic analysis indicated that the novel gene is restricted to the Escherichia/Shigella clade and might have originated recently by overprinting leading to the expression of part of the antisense strand of ECs5115.

Conclusions

Here, we present evidence of a novel small protein-coding gene laoB encoded in the antisense frame −2 of the annotated gene ECs5115. Clearly, laoB is evolutionarily young and it originated in the Escherichia/Shigella clade by overprinting, a process which may cause the de novo evolution of bacterial genes like laoB.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1134-0) contains supplementary material, which is available to authorized users.

Identification of streptococcal small RNAs that are putative targets of RNase III through bioinformatics analysis of RNA sequencing data

Background

Small noncoding regulatory RNAs (sRNAs) are post-transcriptional regulators, regulating mRNAs, proteins, and DNA in bacteria. One class of sRNAs, trans-acting sRNAs, are the most abundant sRNAs transcribed from the intergenic regions (IGRs) of the bacterial genome. In Streptococcus pyogenes, a common and potentially deadly pathogen, many sRNAs have been identified, but only a few have been studied. The goal of this study is to identify trans-acting sRNAs that can be substrates of RNase III. The endoribonuclease RNase III cleaves double stranded RNAs, which can be formed during the interaction between an sRNA and target mRNAs.

Results

For this study, we created an RNase III null mutant of Streptococcus pyogenes and its RNA sequencing (RNA-Seq) data were analyzed and compared to that of the wild-type. First, we developed a custom script that can detect intergenic regions of the S. pyogenes genome. A differential expression analysis with Cufflinks and Stringtie was then performed to identify the intergenic regions whose expression was influenced by the RNase III gene deletion.

Conclusion

This analysis yielded 12 differentially expressed regions with >|2| fold change and p ≤ 0.05. Using Artemis and Bamview genome viewers, these regions were visually verified leaving 6 putative sRNAs. This study not only expanded our knowledge on novel sRNAs but would also give us new insight into sRNA degradation.

EsrE-A yigP Locus-Encoded Transcript-Is a 3' UTR sRNA Involved in the Respiratory Chain of E. coli

The yigP locus is widely conserved among γ-proteobacteria. Mutation of the yigP locus impacts aerobic growth of Gram-negative bacteria. However, the underlying mechanism of how the yigP locus influences aerobic growth remains largely unknown. Here, we demonstrated that the yigP locus in Escherichia coli encodes two transcripts; the mRNA of ubiquinone biosynthesis protein, UbiJ, and the 3′ untranslated region small regulatory RNA (sRNA), EsrE. EsrE is an independent transcript that is transcribed using an internal promoter of the yigP locus. Surprisingly, we found that both the EsrE sRNA and UbiJ protein were required for Q₈ biosynthesis, and were sufficient to rescue the growth defect ascribed to deletion of the yigP locus. Moreover, our data showed that EsrE targeted multiple mRNAs involved in several cellular processes including murein biosynthesis and the tricarboxylic acid cycle. Among these targets, sdhD mRNA that encodes one subunit of succinate dehydrogenase (SDH), was significantly activated. Our findings provided an insight into the important function of EsrE in bacterial adaptation to various environments, as well as coordinating different aspects of bacterial physiology.

Intermolecular base stacking mediates RNA-RNA interaction in a crystal structure of the RNA chaperone Hfq

The RNA-chaperone Hfq catalyses the annealing of bacterial small RNAs (sRNAs) with target mRNAs to regulate gene expression in response to environmental stimuli. Hfq acts on a diverse set of sRNA-mRNA pairs using a variety of different molecular mechanisms. Here, we present an unusual crystal structure showing two Hfq-RNA complexes interacting via their bound RNA molecules. The structure contains two Hfq₆:A₁₈ RNA assemblies positioned face-to-face, with the RNA molecules turned towards each other and connected via interdigitating base stacking interactions at the center. Biochemical data further confirm the observed interaction, and indicate that RNA-mediated contacts occur between Hfq-RNA complexes with various (ARN)_X motif containing RNA sequences in vitro, including the stress response regulator OxyS and its target, fhlA. A systematic computational survey also shows that phylogenetically conserved (ARN)_X motifs are present in a subset of sRNAs, some of which share similar modular architectures. We hypothesise that Hfq can co-opt RNA-RNA base stacking, an unanticipated structural trick, to promote the interaction of (ARN)_X motif containing sRNAs with target mRNAs on a “speed-dating” fashion, thereby supporting their regulatory function.

Crosstalk between the tricarboxylic acid cycle and peptidoglycan synthesis in Caulobacter crescentus through the homeostatic control of α-ketoglutarate

To achieve robust replication, bacteria must integrate cellular metabolism and cell wall growth. While these two processes have been well characterized, the nature and extent of cross-regulation between them is not well understood. Here, using classical genetics, CRISPRi, metabolomics, transcriptomics and chemical complementation approaches, we show that a loss of the master regulator Hfq in Caulobacter crescentus alters central metabolism and results in cell shape defects in a nutrient-dependent manner. We demonstrate that the cell morphology phenotype in the hfq deletion mutant is attributable to a disruption of α-ketoglutarate (KG) homeostasis. In addition to serving as a key intermediate of the tricarboxylic acid (TCA) cycle, KG is a by-product of an enzymatic reaction required for the synthesis of peptidoglycan, a major component of the bacterial cell wall. Accumulation of KG in the hfq deletion mutant interferes with peptidoglycan synthesis, resulting in cell morphology defects and increased susceptibility to peptidoglycan-targeting antibiotics. This work thus reveals a direct crosstalk between the TCA cycle and cell wall morphogenesis. This crosstalk highlights the importance of metabolic homeostasis in not only ensuring adequate availability of biosynthetic precursors, but also in preventing interference with cellular processes in which these intermediates arise as by-products.

Identification and characterisation of isoprene-degrading bacteria in an estuarine environment

Approximately one‐third of volatile organic compounds (VOCs) emitted to the atmosphere consists of isoprene, originating from the terrestrial and marine biosphere, with a profound effect on atmospheric chemistry. However, isoprene provides an abundant and largely unexplored source of carbon and energy for microbes. The potential for isoprene degradation in marine and estuarine samples from the Colne Estuary, UK, was investigated using DNA‐Stable Isotope Probing (DNA‐SIP). Analysis at two timepoints showed the development of communities dominated by Actinobacteria including members of the genera Mycobacterium, Rhodococcus, Microbacterium and Gordonia. Representative isolates, capable of growth on isoprene as sole carbon and energy source, were obtained from marine and estuarine locations, and isoprene‐degrading strains of Gordonia and Mycobacterium were characterised physiologically and their genomes were sequenced. Genes predicted to be required for isoprene metabolism, including four‐component isoprene monooxygenases (IsoMO), were identified and compared with previously characterised examples. Transcriptional and activity assays of strains growing on isoprene or alternative carbon sources showed that growth on isoprene is an inducible trait requiring a specific IsoMO. This study is the first to identify active isoprene degraders in estuarine and marine environments using DNA‐SIP and to characterise marine isoprene‐degrading bacteria at the physiological and molecular level.

Quantitative Proteomics Reveals Potential Crosstalk between a Small RNA CoaR and a Two-Component Regulator Slr1037 in Synechocystis sp. PCC6803

Bacterial small RNAs (sRNAs) and two-component systems (TCSs) were two vital regulatory mechanisms employed by microorganisms to respond to environmental changes and stresses. As a promising "autotrophic cell factory", photosynthetic cyanobacteria have attracted a lot of attention these years. Although most studies focused on studying the roles of sRNAs or TCS regulators in stress response in photosynthetic cyanobacteria, limited work has elucidated their potential crosstalk. Our previous work has identified a negative sRNA regulator CoaR and a positive response regulator Slr1037 both related to 1-butanol stress regulation in Synechocystis sp. PCC6803. In this work, the potential crosstalk between CoaR and Slr1307 (i.e., the coregulated genes mediated by CoaR and Slr1037) was identified and validated through quantitative proteomics and quantitative real-time PCR (qRT-PCR), respectively. The results showed that the sensitive phenotype to 1-butanol of Δslr1037 could be rescued by suppressing coaR in Δslr1037, probably due to the fact that some target genes of Slr1037 could be reactivated by repression of CoaR. Twenty-eight coregulated proteins mediated by CoaR and Slr1037 were found through quantitative proteomics, and 10 of the annotated proteins were validated via qRT-PCR. This study proved the existence of crosstalk between sRNAs and response regulators and provided new insights into the coregulation of biofuel resistance in cyanobacteria.

The Hfq regulon of Neisseria meningitidis

The conserved RNA‐binding protein, Hfq, has multiple regulatory roles within the prokaryotic cell, including promoting stable duplex formation between small RNAs and mRNAs, and thus hfq deletion mutants have pleiotropic phenotypes. Previous proteome and transcriptome studies of Neisseria meningitidis have generated limited insight into differential gene expression due to Hfq loss. In this study, reversed‐phase liquid chromatography combined with data‐independent alternate scanning mass spectrometry (LC‐MS^E) was utilized for rapid high‐resolution quantitative proteomic analysis to further elucidate the differentially expressed proteome of a meningococcal hfq deletion mutant. Whole‐cell lysates of N. meningitidis serogroup B H44/76 wild‐type (wt) and H44/76Δhfq (Δhfq) grown in liquid growth medium were subjected to tryptic digestion. The resulting peptide mixtures were separated by liquid chromatography (LC) prior to analysis by mass spectrometry (MS^E). Differential expression was analyzed by Student's t‐test with control for false discovery rate (FDR). Reliable quantitation of relative expression comparing wt and Δhfq was achieved with 506 proteins (20%). Upon FDR control at q ≤ 0.05, 48 up‐ and 59 downregulated proteins were identified. From these, 81 were identified as novel Hfq‐regulated candidates, while 15 proteins were previously found by SDS/PAGE/MS and 24 with microarray analyses. Thus, using LC‐MS^E we have expanded the repertoire of Hfq‐regulated proteins. In conjunction with previous studies, a comprehensive network of Hfq‐regulated proteins was constructed and differentially expressed proteins were found to be involved in a large variety of cellular processes. The results and comparisons with other gram‐negative model systems, suggest still unidentified sRNA analogs in N. meningitidis.

The Small Protein SgrT Controls Transport Activity of the Glucose-Specific Phosphotransferase System

The bacterial small RNA (sRNA) SgrS has been a fruitful model for discovery of novel RNA-based regulatory mechanisms and new facets of bacterial physiology and metabolism. SgrS is one of only a few characterized dual-function sRNAs. SgrS can control gene expression posttranscriptionally via sRNA-mRNA base-pairing interactions. Its second function is coding for the small protein SgrT. Previous work demonstrated that both functions contribute to relief of growth inhibition caused by glucose-phosphate stress, a condition characterized by disrupted glycolytic flux and accumulation of sugar phosphates. The base-pairing activity of SgrS has been the subject of numerous studies, but the activity of SgrT is less well characterized. Here, we provide evidence that SgrT acts to specifically inhibit the transport activity of the major glucose permease PtsG. Superresolution microscopy demonstrated that SgrT localizes to the cell membrane in a PtsG-dependent manner. Mutational analysis determined that residues in the N-terminal domain of PtsG are important for conferring sensitivity to SgrT-mediated inhibition of transport activity. Growth assays support a model in which SgrT-mediated inhibition of PtsG transport activity reduces accumulation of nonmetabolizable sugar phosphates and promotes utilization of alternative carbon sources by modulating carbon catabolite repression. The results of this study expand our understanding of a basic and well-studied biological problem, namely, how cells coordinate carbohydrate transport and metabolism. Further, this work highlights the complex activities that can be carried out by sRNAs and small proteins in bacteria.IMPORTANCE Sequencing, annotation and investigation of hundreds of bacterial genomes have identified vast numbers of small RNAs and small proteins, the majority of which have no known function. In this study, we explore the function of a small protein that acts in tandem with a well-characterized small RNA during metabolic stress to help bacterial cells maintain balanced metabolism and continue growing. Our results indicate that this protein acts on the glucose transport system, inhibiting its activity under stress conditions in order to allow cells to utilize alternative carbon sources. This work sheds new light on a key biological problem: how cells coordinate carbohydrate transport and metabolism. The study also expands our understanding of the functional capacities of small proteins.