The impact of Hfq-mediated sRNA-mRNA interactome on the virulence of enteropathogenic Escherichia coli

A trans-acting sRNA SaaS targeting hilD, cheA and csgA to inhibit biofilm formation of S. Enteritidis

INTRODUCTION

Salmonella Enteritidis has brought great harm to public health, animal production and food safety worldwide. The biofilm formed by Salmonella Enteritidis plays a critical role in microbial cross-contamination. Small non-coding RNAs (sRNAs) have been demonstrated to be responsible for regulating the formation of biofilm. The sRNA SaaS has been identified previously, that promotes pathogenicity by regulating invasion and virulence factors. However, whether the SaaS is implicated in regulating biofilm formation in abiotic surfaces remains unclear.

OBJECTIVES

This study aimed to clarify the effect of SaaS in Salmonella Enteritidis and explore the modulatory mechanism on the biofilm formation.

METHODS

Motility characteristics and total biomass of biofilm of test strains were investigated by the phenotypes in three soft agar plates and crystal violet staining in polystyrene microplates. Studies of microscopic structure and extracellular polymeric substances (EPS) of biofilm on solid surfaces were carried out using confocal laser scanning microscope (CLSM) and Raman spectra. Transcriptomics and proteomics were applied to analyze the changes of gene expression and EPS component. The RNA-protein pull-down and promoter-reporter β-galactosidase activity assays were employed to analyze RNA binding proteins and identify target mRNAs, respectively.

RESULTS

SaaS inhibits biofilm formation by repressing the adhesion potential and the secretion of EPS components. Integration of transcriptomics and proteomics analysis revealed that SaaS strengthened the expression of the flagellar synthesis system and downregulated the expression of curli amyloid fibers. Furthermore, RNA-protein pull-down interactome datasets indicated that SaaS binds to Hfq (an RNA molecular chaperone protein, known as a host factor for phage Qbeta RNA replication) uniquely among 193 candidate proteins, and promoter-reporter β-galactosidase activity assay confirmed target mRNAs including hilD, cheA, and csgA.

CONCLUSION

SaaS inhibits the properties of bacterial mobility, perturbs the secretion of EPS, and contributes to the inhibition of biofilm formation by interacting with target mRNA (hilD, cheA, and csgA) through the Hfq-mediated pathway.

ChimericFragments: computation, analysis and visualization of global RNA networks

RNA-RNA interactions are a key feature of post-transcriptional gene regulation in all domains of life. While ever more experimental protocols are being developed to study RNA duplex formation on a genome-wide scale, computational methods for the analysis and interpretation of the underlying data are lagging behind. Here, we present ChimericFragments, an analysis framework for RNA-seq experiments that produce chimeric RNA molecules. ChimericFragments implements a novel statistical method based on the complementarity of the base-pairing RNAs around their ligation site and provides an interactive graph-based visualization for data exploration and interpretation. ChimericFragments detects true RNA-RNA interactions with high precision and is compatible with several widely used experimental procedures such as RIL-seq, LIGR-seq or CLASH. We further demonstrate that ChimericFragments enables the systematic detection of novel RNA regulators and RNA-target pairs with crucial roles in microbial physiology and virulence. ChimericFragments is written in Julia and available at: https://github.com/maltesie/ChimericFragments.

RIL-seq reveals extensive involvement of small RNAs in virulence and capsule regulation in hypervirulent Klebsiella pneumoniae

Hypervirulent Klebsiella pneumoniae (hvKp) can infect healthy individuals, in contrast to classical strains that commonly cause nosocomial infections. The recent convergence of hypervirulence with carbapenem-resistance in K. pneumoniae can potentially create 'superbugs' that are challenging to treat. Understanding virulence regulation of hvKp is thus critical. Accumulating evidence suggest that posttranscriptional regulation by small RNAs (sRNAs) plays a role in bacterial virulence, but it has hardly been studied in K. pneumoniae. We applied RIL-seq to a prototypical clinical isolate of hvKp to unravel the Hfq-dependent RNA-RNA interaction (RRI) network. The RRI network is dominated by sRNAs, including predicted novel sRNAs, three of which we validated experimentally. We constructed a stringent subnetwork composed of RRIs that involve at least one hvKp virulence-associated gene and identified the capsule gene loci as a hub target where multiple sRNAs interact. We found that the sRNA OmrB suppressed both capsule production and hypermucoviscosity when overexpressed. Furthermore, OmrB base-pairs within kvrA coding region and partially suppresses translation of the capsule regulator KvrA. This agrees with current understanding of capsule as a major virulence and fitness factor. It emphasizes the intricate regulatory control of bacterial phenotypes by sRNAs, particularly of genes critical to bacterial physiology and virulence.

Hfq-Antisense RNA I Binding Regulates RNase E-Dependent RNA Stability and ColE1 Plasmid Copy Number

The mechanisms and consequences of gene regulation by Hfq on trans-encoded small RNAs (sRNAs) have been well studied and documented. Recent employment of Genomic SELEX to search for Hfq-binding motifs has indicated that Hfq might frequently regulate gene expression controlled by cis-antisense RNAs. Here, we use the classic ColE1 plasmid antisense RNA-based regulation model (i.e., RNA I) to study the role of Hfq in controlling antisense regulatory functions. We show that Hfq exhibits a high binding affinity for RNA I and that binding limits RNase E cleavage, thereby stabilizing RNA I and reducing the plasmid copy number. Full-length RNA I displays a binding affinity for Hfq in the sub-micromolar range. In vivo overexpression of Hfq prolongs RNA I stability and reduces the ColE1 plasmid copy number, whereas deletion of hfq reduces RNA I stability and increases the plasmid copy number. RNA I predominantly binds to the proximal face of Hfq and exhibits competitive ability against a chromosome-borne proximal face-bound sRNA (DsrA) for Hfq binding. Through its strong promoter and high gene dosage features, plasmid-encoded antisense RNA I results in high RNA I expression, so it may antagonize the effects of trans-encoded RNAs in controlling target gene expression.

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.

Agnodice: indexing experimentally supported bacterial sRNA-RNA interactions

In the last decade, the immense growth in the field of bacterial small RNAs (sRNAs), along with the biotechnological breakthroughs in Deep Sequencing permitted the deeper understanding of sRNA-RNA interactions. However, microbiology is currently lacking a thoroughly curated collection of this rapidly expanding universe. We present Agnodice (https://dianalab.e-ce.uth.gr/agnodice), our effort to systematically catalog and annotate experimentally supported bacterial sRNA-RNA interactions. Agnodice, for the first time, incorporates thousands of bacterial sRNA-RNA interactions derived from a diverse set of experimental methodologies including state-of-the-art Deep Sequencing interactome identification techniques. It comprises 39,600 entries which are annotated at strain-level resolution and pertain to 399 sRNAs and 12,137 target RNAs identified in 71 bacterial strains. The database content is exclusively experimentally supported, incorporating interactions derived via low yield as well as state-of-the-art high-throughput methods. The entire content of the database is freely accessible and can be directly downloaded for further analysis. Agnodice will serve as a valuable source, enabling microbiologists to form novel hypotheses, design/identify novel sRNA-based drug targets, and explore the therapeutic potential of microbiomes from the perspective of small regulatory RNAs.IMPORTANCEAgnodice (https://dianalab.e-ce.uth.gr/agnodice) is an effort to systematically catalog and annotate experimentally supported bacterial small RNA (sRNA)-RNA interactions. Agnodice, for the first time, incorporates thousands of bacterial sRNA-RNA interactions derived from a diverse set of experimental methodologies including state-of-the-art Next Generation Sequencing interactome identification techniques.

Defining Bacterial RNA-RNA Interactomes Using CLASH

Methicillin-resistant Staphylococcus aureus (MRSA) is a bacterial pathogen accounting for high mortality rates among infected patients. Transcriptomic regulation by small RNAs (sRNAs) has been shown to regulate networks promoting antibiotic resistance and virulence in S. aureus. Yet, the biological role of most sRNAs during MRSA host infection remains unknown. To fill this gap, in collaboration with the lab of Jai Tree, we performed comprehensive RNA-RNA interactome analyses in MRSA using CLASH under conditions that mimic the host environment. Here we present a detailed version of this optimized CLASH (cross-linking, ligation, and sequencing of hybrids) protocol we recently developed, which has been tailored to explore the RNA interactome in S. aureus as well as other Gram-positive bacteria. Alongside, we introduce a compilation of helpful Python functions for analyzing folding energies of putative RNA-RNA interactions and streamlining sRNA and mRNA seed discovery in CLASH data. In the accompanying computational demonstration, we aim to establish a standardized strategy to evaluate the likelihood that observed chimeras arise from true RNA-RNA interactions.

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.

TRS: a method for determining transcript termini from RNAtag-seq sequencing data

In bacteria, determination of the 3' termini of transcripts plays an essential role in regulation of gene expression, affecting the functionality and stability of the transcript. Several experimental approaches were developed to identify the 3' termini of transcripts, however, these were applied only to a limited number of bacteria and growth conditions. Here we present a straightforward approach to identify 3' termini from widely available RNA-seq data without the need for additional experiments. Our approach relies on the observation that the RNAtag-seq sequencing protocol results in overabundance of reads mapped to transcript 3' termini. We present TRS (Termini by Read Starts), a computational pipeline exploiting this property to identify 3' termini in RNAtag-seq data, and show that the identified 3' termini are highly reliable. Since RNAtag-seq data are widely available for many bacteria and growth conditions, our approach paves the way for studying bacterial transcription termination in an unprecedented scope.

σ28-dependent small RNA regulation of flagella biosynthesis

Flagella are important for bacterial motility as well as for pathogenesis. Synthesis of these structures is energy intensive and, while extensive transcriptional regulation has been described, little is known about the posttranscriptional regulation. Small RNAs (sRNAs) are widespread posttranscriptional regulators, most base pairing with mRNAs to affect their stability and/or translation. Here, we describe four UTR-derived sRNAs (UhpU, MotR, FliX and FlgO) whose expression is controlled by the flagella sigma factor σ28 (fliA) in Escherichia coli. Interestingly, the four sRNAs have varied effects on flagellin protein levels, flagella number and cell motility. UhpU, corresponding to the 3´ UTR of a metabolic gene, likely has hundreds of targets including a transcriptional regulator at the top flagella regulatory cascade connecting metabolism and flagella synthesis. Unlike most sRNAs, MotR and FliX base pair within the coding sequences of target mRNAs and act on ribosomal protein mRNAs connecting ribosome production and flagella synthesis. The study shows how sRNA-mediated regulation can overlay a complex network enabling nuanced control of flagella synthesis.

Advantages and limitations of UV cross-linking analysis of protein-RNA interactomes in microbes

RNA-binding proteins (RBPs) govern the lifespan of nearly all transcripts and play key roles in adaptive responses in microbes. A robust approach to examine protein-RNA interactions involves irradiating cells with UV light to form covalent adducts between RBPs and their cognate RNAs. Combined with RNA or protein purification, these procedures can provide global RBP censuses or transcriptomic maps for all target sequences of a single protein in living cells. The recent development of novel methods has quickly populated the RBP landscape in microorganisms. Here, we provide an overview of prominent UV cross-linking techniques which have been applied to investigate RNA interactomes in microbes. By assessing their advantages and caveats, this technical evaluation intends to guide the selection of appropriate methods and experimental design as well as to encourage the use of complementary UV-dependent techniques to inspect RNA-binding activity.

Hfq-licensed RNA-RNA interactome in Pseudomonas aeruginosa reveals a keystone sRNA

The RNA chaperone Hfq plays important regulatory roles in many bacteria by facilitating the base pairing between small RNAs (sRNAs) and their cognate mRNA targets. In the gram-negative opportunistic pathogen Pseudomonas aeruginosa, over a hundred putative sRNAs have been identified but for most, their regulatory targets remained unknown. Using RIL-seq with Hfq in P. aeruginosa, we identified the mRNA targets for dozens of previously known and unknown sRNAs. Strikingly, hundreds of the RNA-RNA interactions we discovered involved PhrS. This sRNA was thought to mediate its effects by pairing with a single target mRNA and regulating the abundance of the transcription regulator MvfR required for the synthesis of the quorum sensing signal PQS. We present evidence that PhrS controls many transcripts by pairing with them directly and employs a two-tiered mechanism for governing PQS synthesis that involves control of an additional transcription regulator called AntR. Our findings in P. aeruginosa expand the repertoire of targets for previously known sRNAs, reveal potential regulatory targets for previously unknown sRNAs, and suggest that PhrS may be a keystone sRNA with the ability to pair with an unusually large number of transcripts in this organism.

Involvement of RNA chaperone Hfq in the regulation of antibiotic resistance and virulence in Shigella sonnei

The host factor for RNA phage Qβ replicase (Hfq) is a crucial post-transcriptional regulator in many bacterial pathogens, facilitating the interaction between small non-coding RNAs (sRNAs) and their target mRNAs. Studies have suggested that Hfq plays a role in antibiotic resistance and virulence in bacteria, although its functions in Shigella are not fully understood. In this study, we investigated the functional roles of Hfq in Shigella sonnei (S. sonnei) by constructing an hfq deletion mutant. Our phenotypic assays showed that the hfq deletion mutant was more sensitivity to antibiotics and had impaired virulence. Transcriptome analyses supported the results concerning the phenotype of the hfq mutant and showed that differentially expressed genes were mainly enriched in the KEGG pathways two-component system, ABC transporters, ribosome, and Escherichia coli biofilm formation. Additionally, we predicted eleven novel Hfq-dependent sRNAs, which were potentially involved in the regulation of antibiotic resistance and/or virulence in S. sonnei. Our findings suggest that Hfq plays a post-transcriptional role in regulating antibiotic resistance and virulence in S. sonnei, and could provide a basis for future studies on Hfq-sRNA-mRNA regulatory networks in this important pathogen.

The phosphate-induced small RNA EsrL promotes E. coli virulence, biofilm formation, and intestinal colonization

Escherichia coli are part of the normal intestinal microbiome, but some enterohemorrhagic E. coli (EHEC) and enteropathogenic E. coli (EPEC) strains can cause potentially life-threatening gastroenteritis. Virulence factors underlying the ability of EHEC and EPEC to cause disease include those encoded in the locus of the enterocyte effacement (LEE) pathogenicity island. Here, we demonstrated that EsrL, a small RNA present in many E. coli strains, promoted pathogenicity, adhesion, and biofilm formation in EHEC and EPEC. PhoB, the response regulator of the two-component system that controls cellular responses to phosphate, directly repressed esrL expression under low-phosphate conditions. A phosphate-rich environment, similar to that of the human intestine, relieved PhoB-mediated repression of esrL. EsrL interacted with and stabilized the LEE-encoded regulator (ler) transcript, which encodes a transcription factor for LEE genes, leading to increased bacterial adhesion to cultured cells and colonization of the rabbit colon. EsrL also bound to and stabilized the fimC transcript, which encodes a chaperone that is required for the assembly of type 1 pili, resulting in enhanced cell adhesion in pathogenic E. coli and enhanced biofilm formation in pathogenic and nonpathogenic E. coli. Our findings demonstrate that EsrL stimulates the expression of virulence genes in both EHEC and EPEC under phosphate-rich conditions, thus promoting the pathogenicity of EHEC and EPEC in the nutrient-rich gut environment.

An RNA sponge controls quorum sensing dynamics and biofilm formation in Vibrio cholerae

Small regulatory RNAs (sRNAs) acting in concert with the RNA chaperone Hfq are prevalent in many bacteria and typically act by base-pairing with multiple target transcripts. In the human pathogen Vibrio cholerae, sRNAs play roles in various processes including antibiotic tolerance, competence, and quorum sensing (QS). Here, we use RIL-seq (RNA-interaction-by-ligation-and-sequencing) to identify Hfq-interacting sRNAs and their targets in V. cholerae. We find hundreds of sRNA-mRNA interactions, as well as RNA duplexes formed between two sRNA regulators. Further analysis of these duplexes identifies an RNA sponge, termed QrrX, that base-pairs with and inactivates the Qrr1-4 sRNAs, which are known to modulate the QS pathway. Transcription of qrrX is activated by QrrT, a previously uncharacterized LysR-type transcriptional regulator. Our results indicate that QrrX and QrrT are required for rapid conversion from individual to community behaviours in V. cholerae.

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

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

Hfq Is a Critical Modulator of Pathogenicity of Dickeya oryzae in Rice Seeds and Potato Tubers

The frequent outbreaks of soft-rot diseases caused by Dickeya oryzae have emerged as severe problems in plant production in recent years and urgently require the elucidation of the virulence mechanisms of D. oryzae. Here, we report that Hfq, a conserved RNA chaperone protein in bacteria, is involved in modulating a series of virulence-related traits and bacterial virulence in D. oryzae EC1. The findings show that the null mutation of the hfqEC1 gene totally abolished the production of zeamine phytotoxins and protease, significantly attenuated the production of two other types of cell wall degrading enzymes, i.e., pectate lyase and cellulase, as well as attenuating swarming motility, biofilm formation, the development of hypersensitive response to Nicotiana benthamiana, and bacterial infections in rice seeds and potato tubers. QRT-PCR analysis and promoter reporter assay further indicated that HfqEC1 regulates zeamine production via modulating the expression of the key zeamine biosynthesis (zms) cluster genes. Taken together, these findings highlight that the Hfq of D. oryzae is one of the key regulators in modulating the production of virulence determinants and bacterial virulence in rice seeds and potato tubers.

The RNA Chaperone Protein Hfq Regulates the Characteristic Sporulation and Insecticidal Activity of Bacillus thuringiensis

Bacillus thuringiensis (Bt) is one of the most widely used bio-insecticides at present. It can produce many virulence factors and insecticidal crystal proteins during growth and sporulation. Hfq, on the other hand, is a bacterial RNA chaperone that can regulate the function of different kinds of RNAs, thereby affecting various bacterial phenotypes. To further explore the physiological functions of Hfq in Bt, we took BMB171 as the starting strain, knocked out one, two, or three hfq genes in its genome in different combinations, and compared the phenotypic differences between the deletion mutant strains and the starting strain. We did observe significant changes in several phenotypes, including motility, biofilm formation, sporulation, and insecticidal activity against cotton bollworm, among others. Afterward, we found through transcriptome studies that when all hfq genes were deleted, 32.5% of the genes in Bt were differentially transcribed, with particular changes in the sporulation-related and virulence-related genes. The above data demonstrated that Hfq plays a pivotal role in Bt and can regulate its various physiological functions. Our study on the regulatory mechanism of Hfq in Bt, especially the mining of the regulatory network of its sporulation and insecticidal activity, could lay a theoretical foundation for the better utilization of Bt as an effective insecticide.

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.