Mining RNA-seq data reveals the massive regulon of GcvB small RNA and its physiological significance in maintaining amino acid homeostasis in Escherichia coli

Post-transcriptional regulation of aromatic amino acid metabolism by GcvB small RNA in Escherichia coli

Escherichia coli synthesizes aromatic amino acids (AAAs) through the common pathway to produce the precursor, chorismate, and the three terminal pathways to convert chorismate into Phe, Tyr, and Trp. E. coli also imports exogenous AAAs through five transporters. GcvB small RNA post-transcriptionally regulates more than 50 genes involved in amino acid uptake and biosynthesis in E. coli, but the full extent of GcvB regulon is still underestimated. This study examined all genes involved in AAA biosynthesis and transport using translation reporter assay and qRT-PCR analysis. In addition to previously verified targets, aroC, aroP, and trpE, we identified new target genes that were significantly repressed by GcvB primarily via the R1 seed region. Exceptionally, GcvB strongly inhibits the expression of aroG, which encodes the major isozyme of the first reaction in the common pathway, through direct base pairing between the aroG translation initiation region and the GcvB R3 seed sequence. RNase E mediates the degradation of target mRNAs except aroC and aroP via its C-terminal domain. GcvB overexpression prolongs the lag phase and reduces the growth rate in minimal media supplemented with AAAs and confers resistance to an antibiotic compound, azaserine, by repressing AAA transporters.IMPORTANCEE. coli strains have been genetically modified in relevant transcription factors and biosynthetic enzymes for industrial use in the fermentative production of aromatic amino acids (AAAs) and their derivative compounds. This study focuses on GcvB small RNA, a global regulator of amino acid metabolism in E. coli, and identifies new GcvB targets involved in AAA biosynthesis and uptake. GcvB represses the expression of the first and last enzymes of the common pathway and the first enzymes of Trp and Phe terminal pathways. GcvB also limits import of AAAs. This paper documents the impact of RNA-mediated regulation on AAA metabolism in E. coli.

Regulation of an RNA toxin-antitoxin system, SdsR-RyeA, by a small RNA GcvB

The toxin-antitoxin (TA) system regulates many physiological processes in free-living bacteria. One such TA system in Escherichia coli comprises an RNA toxin SdsR and an antitoxin RyeA. An overabundance of SdsR is toxic to the cells. RyeA normalizes SdsR abundance and helps the cells to adapt to altered conditions. The current study showed that a novel small RNA (sRNA) regulator GcvB directly interacts with RyeA to maintain its abundance in the cells under normal or low pH conditions. The deletion of the gcvB allele in the E. coli chromosome resulted in a ∼3-fold decrease in intrabacterial RyeA accumulation. An ectopic expression of GcvB in ΔgcvB strain reinstated RyeA abundance to its normal level. Induction of GcvB in the cells upon exposure to low pH resulted in a simultaneous increase in intracellular RyeA. While GcvB increases RyeA abundance in the cells, SdsR accumulation is divergently regulated by GcvB. The absence of the gcvB gene in E. coli leads to upregulation of SdsR and vice versa. The GcvB-mediated decrease of SdsR accumulation stems from the increased RyeA-driven normalization of SdsR. This study delineates a novel mechanism for the regulation of the expression of an RNA toxin SdsR by another sRNA regulator GcvB through a feed-forward control.

Revealing sRNA expression profiles of NDM-5-producing CRKP and explore the role of sRNA207 in NDM-5-producing CRKP resistance

New Delhi metallo-beta-lactamase-5 (NDM-5)-producing carbapenem-resistant Klebsiella pneumoniae (CRKP) is characterized by high virulence, high morbidity, and mortality, and the detection rate in children has increased in recent years. Therefore, it is urgent to develop new therapeutic targets and strategies. Non-coding small RNA (sRNA)-mediated RNA-based therapies offer a new direction for the treatment of bacterial infections, especially resistant bacteria. This study first analyzed the transcriptional expression profiles of NDM-5-producing CRKP and Carbapenem-susceptible Klebsiella pneumoniae (CSKP) isolates from the clinic by RNA-seq. A total of 4,623 genes were obtained, of which 307 genes were differentially expressed in NDM-5-producing CRKP, and these differentially expressed genes are mainly related to metabolism. Then, by analyzing the length and secondary structure of genes that could not match the reference gene and non-redundant protein database, we obtained 268 sRNAs, of which 13 sRNAs were differentially expressed in NDM-5-producing CRKP. After the expression level of differentially expressed sRNA was verified by RT-PCR to be consistent with that of RNA-seq, we chose sRNA207 as our research target. By knockdown of sRNA207 and smf-1 (the predicted target mRNA of sRNA207) in the strain, we found that increased expression of sRNA207 promoted biofilm formation by stabilizing expression of smf-1, which in turn affected the resistance of NDM-5-producing CRKP to carbapenems. This provides a new approach to treat CRKP infection.

IMPORTANCE

sRNAs form a regulatory network that regulates bacterial virulence, drug resistance, and other functions by targeting mRNAs. However, sRNA expression profile and function of NDM-5-producing carbapenem-resistant Klebsiella pneumoniae (CRKP) are still unknown. In this study, we analyzed the sRNA expression profiles of NDM-5-producing CRKP obtained from clinical by referring to the methods of previous articles. A total of 268 candidates sRNAs were obtained, of which 248 were newly discovered. More importantly, 13 sRNAs were differentially expressed in NDM-5-producing CRKP compared with CSKP. We knocked down sRNA207 in NDM-5-producing CRKP to validate its effect on smf-1, biofilm, and resistance of strains. We also confirmed the role of smf-1 in biofilm formation and drug resistance of NDM-5-producing CRKP by constructing smf-1-knockdown strain. The results suggest that smf-1 is the target gene of sRNA207. Increased expression of sRNA207 in NDM-5-producing CRKP stabilizes smf-1 expression, which in turn affects the resistance of the strains through biofilm formation.

Stochastic nature and physiological implications of 5'-NAD RNA cap in bacteria

RNA 5'-modification with NAD+/NADH (oxidized/reduced nicotinamide adenine dinucleotide) has been found in bacteria, eukaryotes and viruses. 5'-NAD is incorporated into RNA by RNA polymerases (RNAPs) during the initiation of synthesis. It is unknown (i) which factors and physiological conditions permit substantial NAD incorporation into RNA in vivo and (ii) how 5'-NAD impacts gene expression and the fate of RNA in bacteria. Here we show in Escherichia coli that RNA NADylation is stimulated by low cellular concentration of the competing substrate ATP, and by weakening ATP contacts with RNAP active site. Additionally, RNA NADylation may be influenced by DNA supercoiling. RNA NADylation does not interfere with posttranscriptional RNA processing by major ribonuclease RNase E. It does not impact the base-pairing between RNAI, the repressor of plasmid replication, and its antisense target, RNAII. Leaderless NADylated model mRNA cI-lacZ is recognized by the 70S ribosome and is translated with the same efficiency as triphosphorylated cI-lacZ mRNA. Translation exposes the 5'-NAD of this mRNA to de-capping by NudC enzyme. We suggest that NADylated mRNAs are rapidly degraded, consistent with their low abundance in published datasets. Furthermore, we observed that ppGpp inhibits NudC de-capping activity, contributing to the growth phase-dependency of NADylated RNA levels.

Unraveling the interplay between a small RNA and RNase E in bacteria

Small RNAs (sRNAs) are major regulators of gene expression in bacteria, exerting their regulation primarily via base pairing with their target transcripts and modulating translation. Accumulating evidence suggest that sRNAs can also affect the stability of their target transcripts by altering their accessibility to endoribonucleases. Yet, the effects of sRNAs on transcript stability and the mechanisms underlying them have not been studied in wide scale. Here we employ large-scale RNA-seq-based methodologies in the model bacterium Escherichia coli to quantitatively study the functional interaction between a sRNA and an endoribonuclease in regulating gene expression, using the well-established sRNA, GcvB, and the major endoribonuclease, RNase E. Studying single and double mutants of gcvB and rne and analysing their RNA-seq results by the Double Mutant Cycle approach, we infer distinct modes of the interplay between GcvB and RNase E. Transcriptome-wide mapping of RNase E cleavage sites provides further support to the results of the RNA-seq analysis, identifying cleavage sites in targets in which the functional interaction between GcvB and RNase E is evident. Together, our results indicate that the most dominant mode of GcvB-RNase E functional interaction is GcvB enhancement of RNase E cleavage, which varies in its magnitude between different targets.

Multiscale regulation of nutrient stress responses in Escherichia coli from chromatin structure to small regulatory RNAs

Recent research has indicated the presence of heterochromatin-like regions of extended protein occupancy and transcriptional silencing of bacterial genomes. We utilized an integrative approach to track chromatin structure and transcription in E. coli K-12 across a wide range of nutrient conditions. In the process, we identified multiple loci which act similarly to facultative heterochromatin in eukaryotes, normally silenced but permitting expression of genes under specific conditions. We also found a strong enrichment of small regulatory RNAs (sRNAs) among the set of differentially expressed transcripts during nutrient stress. Using a newly developed bioinformatic pipeline, the transcription factors regulating sRNA expression were bioinformatically predicted, with experimental follow-up revealing novel relationships for 36 sRNA-transcription factors candidates. Direct regulation of sRNA expression was confirmed by mutational analysis for five sRNAs of metabolic interest: IsrB, CsrB and CsrC, GcvB, and GadY. Our integrative analysis thus reveals additional layers of complexity in the nutrient stress response in E. coli and provides a framework for revealing similar poorly understood regulatory logic in other organisms.

Structural and functional insights into transcription activation of the essential LysR-type transcriptional regulators

The enormous LysR-type transcriptional regulators (LTTRs), which are diversely distributed amongst prokaryotes, play crucial roles in transcription regulation of genes involved in basic metabolic pathways, virulence and stress resistance. However, the precise transcription activation mechanism of these genes by LTTRs remains to be explored. Here, we determine the cryo-EM structure of a LTTR-dependent transcription activation complex comprising of Escherichia coli RNA polymerase (RNAP), an essential LTTR protein GcvA and its cognate promoter DNA. Structural analysis shows two N-terminal DNA binding domains of GcvA (GcvA_DBD) dimerize and engage the GcvA activation binding sites, presenting the -35 element for specific recognition with the conserved σ70R4. In particular, the versatile C-terminal domain of α subunit of RNAP directly interconnects with GcvA_DBD, σ70R4 and promoter DNA, providing more interfaces for stabilizing the complex. Moreover, molecular docking supports glycine as one potential inducer of GcvA, and single molecule photobleaching experiments kinetically visualize the occurrence of tetrameric GcvA-engaged transcription activation complex as suggested for the other LTTR homologs. Thus, a general model for tetrameric LTTR-dependent transcription activation is proposed. These findings will provide new structural and functional insights into transcription activation of the essential LTTRs.

Global Hfq-mediated RNA interactome of nitrogen starved Escherichia coli uncovers a conserved post-transcriptional regulatory axis required for optimal growth recovery

The RNA binding protein Hfq has a central role in the post-transcription control of gene expression in many bacteria. Numerous studies have mapped the transcriptome-wide Hfq-mediated RNA-RNA interactions in growing bacteria or bacteria that have entered short-term growth-arrest. To what extent post-transcriptional regulation underpins gene expression in growth-arrested bacteria remains unknown. Here, we used nitrogen (N) starvation as a model to study the Hfq-mediated RNA interactome as Escherichia coli enter, experience, and exit long-term growth arrest. We observe that the Hfq-mediated RNA interactome undergoes extensive changes during N starvation, with the conserved SdsR sRNA making the most interactions with different mRNA targets exclusively in long-term N-starved E. coli. Taking a proteomics approach, we reveal that in growth-arrested cells SdsR influences gene expression far beyond its direct mRNA targets. We demonstrate that the absence of SdsR significantly compromises the ability of the mutant bacteria to recover growth competitively from the long-term N-starved state and uncover a conserved post-transcriptional regulatory axis which underpins this process.

Multilayered regulation of amino acid metabolism in Escherichia coli

Amino acid metabolism in Escherichia coli has long been studied and has established the basis for regulatory mechanisms at the transcriptional, posttranscriptional, and posttranslational levels. In addition to the classical signal transduction cascade involving posttranslational modifications (PTMs), novel PTMs in the two primary nitrogen assimilation pathways have recently been uncovered. The regulon of the master transcriptional regulator NtrC is further expanded by a small RNA derived from the 3´UTR of glutamine synthetase mRNA, which coordinates central carbon and nitrogen metabolism. Furthermore, recent advances in sequencing technologies have revealed the global regulatory networks of transcriptional and posttranscriptional regulators, Lrp and GcvB. This review provides an update of the multilayered and interconnected regulatory networks governing amino acid metabolism in E. coli.

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

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

Central Role of Sibling Small RNAs NgncR_162 and NgncR_163 in Main Metabolic Pathways of Neisseria gonorrhoeae

Small bacterial regulatory RNAs (sRNAs) have been implicated in the regulation of numerous metabolic pathways. In most of these studies, sRNA-dependent regulation of mRNAs or proteins of enzymes in metabolic pathways has been predicted to affect the metabolism of these bacteria. However, only in a very few cases has the role in metabolism been demonstrated. Here, we performed a combined transcriptome and metabolome analysis to define the regulon of the sibling sRNAs NgncR_162 and NgncR_163 (NgncR_162/163) and their impact on the metabolism of Neisseria gonorrhoeae. These sRNAs have been reported to control genes of the citric acid and methylcitric acid cycles by posttranscriptional negative regulation. By transcriptome analysis, we now expand the NgncR_162/163 regulon by several new members and provide evidence that the sibling sRNAs act as both negative and positive regulators of target gene expression. Newly identified NgncR_162/163 targets are mostly involved in transport processes, especially in the uptake of glycine, phenylalanine, and branched-chain amino acids. NgncR_162/163 also play key roles in the control of serine-glycine metabolism and, hence, probably affect biosyntheses of nucleotides, vitamins, and other amino acids via the supply of one-carbon (C₁) units. Indeed, these roles were confirmed by metabolomics and metabolic flux analysis, which revealed a bipartite metabolic network with glucose degradation for the supply of anabolic pathways and the usage of amino acids via the citric acid cycle for energy metabolism. Thus, by combined deep RNA sequencing (RNA-seq) and metabolomics, we significantly extended the regulon of NgncR_162/163 and demonstrated the role of NgncR_162/163 in the regulation of central metabolic pathways of the gonococcus. IMPORTANCE Neisseria gonorrhoeae is a major human pathogen which infects more than 100 million people every year. An alarming development is the emergence of gonococcal strains that are resistant against virtually all antibiotics used for their treatment. Despite the medical importance and the vanishing treatment options of gonococcal infections, the bacterial metabolism and its regulation have been only weakly defined until today. Using RNA-seq, metabolomics, and ¹³C-guided metabolic flux analysis, we here investigated the gonococcal metabolism and its regulation by the previously studied sibling sRNAs NgncR_162/163. The results demonstrate the regulation of transport processes and metabolic pathways involved in the biosynthesis of nucleotides, vitamins, and amino acids by NgncR_162/163. In particular, the combination of transcriptome and metabolic flux analyses provides a heretofore unreached depth of understanding the core metabolic pathways and their regulation by the neisserial sibling sRNAs. This integrative approach may therefore also be suitable for the functional analysis of a growing number of other bacterial metabolic sRNA regulators.

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

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

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

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

GcvB Regulon Revealed by Transcriptomic and Proteomic Analysis in Vibrio alginolyticus

Vibrio alginolyticus is a widely distributed marine bacterium that is a threat to the aquaculture industry as well as human health. Evidence has revealed critical roles for small RNAs (sRNAs) in bacterial physiology and cellular processes by modulating gene expression post-transcriptionally. GcvB is one of the most conserved sRNAs that is regarded as the master regulator of amino acid uptake and metabolism in a wide range of Gram-negative bacteria. However, little information about GcvB-mediated regulation in V. alginolyticus is available. Here we first characterized GcvB in V. alginolyticus ZJ-T and determined its regulon by integrated transcriptome and quantitative proteome analysis. Transcriptome analysis revealed 40 genes differentially expressed (DEGs) between wild-type ZJ-T and gcvB mutant ZJ-T-ΔgcvB, while proteome analysis identified 50 differentially expressed proteins (DEPs) between them, but only 4 of them displayed transcriptional differences, indicating that most DEPs are the result of post-transcriptional regulation of gcvB. Among the differently expressed proteins, 21 are supposed to be involved in amino acid biosynthesis and transport, and 11 are associated with type three secretion system (T3SS), suggesting that GcvB may play a role in the virulence besides amino acid metabolism. RNA-EMSA showed that Hfq binds to GcvB, which promotes its stability.

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

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

Checkpoints That Regulate Balanced Biosynthesis of Lipopolysaccharide and Its Essentiality in Escherichia coli

The outer membrane (OM) of Gram-negative bacteria, such as Escherichia coli, is essential for their viability. Lipopolysaccharide (LPS) constitutes the major component of OM, providing the permeability barrier, and a tight balance exists between LPS and phospholipids amounts as both of these essential components use a common metabolic precursor. Hence, checkpoints are in place, right from the regulation of the first committed step in LPS biosynthesis mediated by LpxC through its turnover by FtsH and HslUV proteases in coordination with LPS assembly factors LapB and LapC. After the synthesis of LPS on the inner leaflet of the inner membrane (IM), LPS is flipped by the IM-located essential ATP-dependent transporter to the periplasmic face of IM, where it is picked up by the LPS transport complex spanning all three components of the cell envelope for its delivery to OM. MsbA exerts its intrinsic hydrocarbon ruler function as another checkpoint to transport hexa-acylated LPS as compared to underacylated LPS. Additional checkpoints in LPS assembly are: LapB-assisted coupling of LPS synthesis and translocation; cardiolipin presence when LPS is underacylated; the recruitment of RfaH transcriptional factor ensuring the transcription of LPS core biosynthetic genes; and the regulated incorporation of non-stoichiometric modifications, controlled by the stress-responsive RpoE sigma factor, small RNAs and two-component systems.