A major challenge for RNA-seq analysis of gene expression is to achieve sufficient coverage of informative nonribosomal transcripts. In eukaryotic samples, this is typically achieved by selective oligo(dT)-priming of messenger RNAs to exclude ribosomal RNA (rRNA) during cDNA synthesis. However, this strategy is not compatible with prokaryotes in which functional transcripts are generally not polyadenylated. To overcome this, we adopted DASH (depletion of abundant sequences by hybridization), initially developed for eukaryotic cells, to improve both the sensitivity and depth of bacterial RNA-seq. DASH uses the Cas9 nuclease to remove unwanted cDNA sequences prior to library amplification. We report the design, evaluation, and optimization of DASH experiments for standard bacterial short-read sequencing approaches, including software for automated guide RNA (gRNA) design for Cas9-mediated cleavage in bacterial rDNA sequences. Using these gRNA pools, we effectively removed rRNA reads (56%-86%) in RNA-seq libraries from two different model bacteria, the Gram-negative pathogen Salmonella enterica and the anaerobic gut commensal Bacteroides thetaiotaomicron DASH works robustly, even with subnanogram amounts of input RNA. Its efficiency, high sensitivity, ease of implementation, and low cost (∼$5 per sample) render DASH an attractive alternative to rRNA removal protocols, in particular for material-constrained studies where conventional ribodepletion techniques fail.

Capturing an individual cell's transcriptional history is a challenge exacerbated by the functional heterogeneity of cellular communities. Here, we leverage reprogrammed tracrRNAs (Rptrs) to record selected cellular transcripts as stored DNA edits in single living bacterial cells. Rptrs are designed to base pair with sensed transcripts, converting them into guide RNAs. The guide RNAs then direct a Cas9 base editor to target an introduced DNA target. The extent of base editing can then be read in the future by sequencing. We use this approach, called TIGER (transcribed RNAs inferred by genetically encoded records), to record heterologous and endogenous transcripts in individual bacterial cells. TIGER can quantify relative expression, distinguish single-nucleotide differences, record multiple transcripts simultaneously and read out single-cell phenomena. We further apply TIGER to record metabolic bet hedging and antibiotic resistance mobilization in Escherichia coli as well as host cell invasion by Salmonella. Through RNA recording, TIGER connects current cellular states with past transcriptional states to decipher complex cellular responses in single cells.

Bacterial small RNAs (sRNAs) are key elements of regulatory networks that modulate gene expression. The sRNA RydC of Salmonella sp. and Escherichia coli is an example of this class of riboregulators. Like many other sRNAs, RydC bears a 'seed' region that recognises specific transcripts through base-pairing, and its activities are facilitated by the RNA chaperone Hfq. The crystal structure of RydC in complex with E. coli Hfq at a 3.48 Å resolution illuminates how the protein interacts with and presents the sRNA for target recognition. Consolidating the protein-RNA complex is a host of distributed interactions mediated by the natively unstructured termini of Hfq. Based on the structure and other data, we propose a model for a dynamic effector complex comprising Hfq, small RNA, and the cognate mRNA target.

Many pathogenic bacteria utilize specialized secretion systems to deliver proteins called effectors into eukaryotic cells for manipulation of host pathways. The vast majority of known effector targets are host proteins, whereas a potential targeting of host nucleic acids remains little explored. There is only one family of effectors known to target DNA directly, and effectors binding host RNA are unknown. Here, we take a two-pronged approach to search for RNA-binding effectors, combining biocomputational prediction of RNA-binding domains (RBDs) in a newly assembled comprehensive dataset of bacterial secreted proteins, and experimental screening for RNA binding in mammalian cells. Only a small subset of effectors were predicted to carry an RBD, indicating that if RNA targeting was common, it would likely involve new types of RBDs. Our experimental evaluation of effectors with predicted RBDs further argues for a general paucity of RNA binding activities amongst bacterial effectors. We obtained evidence that PipB2 and Lpg2844, effector proteins of Salmonella and Legionella species, respectively, may harbor novel biochemical activities. Our study presenting the first systematic evaluation of the RNA-targeting potential of bacterial effectors offers a basis for discussion of whether or not host RNA is a prominent target of secreted bacterial proteins.

To understand the gene regulation of an organism of interest, a comprehensive genome annotation is essential. While some features, such as coding sequences, can be computationally predicted with high accuracy based purely on the genomic sequence, others, such as promoter elements or noncoding RNAs, are harder to detect. RNA sequencing (RNA-seq) has proven to be an efficient method to identify these genomic features and to improve genome annotations. However, processing and integrating RNA-seq data in order to generate high-resolution annotations is challenging, time consuming, and requires numerous steps. We have constructed a powerful and modular tool called ANNOgesic that provides the required analyses and simplifies RNA-seq-based bacterial and archaeal genome annotation. It can integrate data from conventional RNA-seq and differential RNA-seq and predicts and annotates numerous features, including small noncoding RNAs, with high precision. The software is available under an open source license (ISCL) at https://pypi.org/project/ANNOgesic/.

Bulk RNA sequencing technologies have provided invaluable insights into host and bacterial gene expression and associated regulatory networks. Nevertheless, the majority of these approaches report average expression across cell populations, hiding the true underlying expression patterns that are often heterogeneous in nature. Due to technical advances, single-cell transcriptomics in bacteria has recently become reality, allowing exploration of these heterogeneous populations, which are often the result of environmental changes and stressors. In this work, we have improved our previously published bacterial single-cell RNA sequencing (scRNA-seq) protocol that is based on multiple annealing and deoxycytidine (dC) tailing-based quantitative scRNA-seq (MATQ-seq), achieving a higher throughput through the integration of automation. We also selected a more efficient reverse transcriptase, which led to reduced cell loss and higher workflow robustness. Moreover, we successfully implemented a Cas9-based rRNA depletion protocol into the MATQ-seq workflow. Applying our improved protocol on a large set of single Salmonella cells sampled over different growth conditions revealed improved gene coverage and a higher gene detection limit compared to our original protocol and allowed us to detect the expression of small regulatory RNAs, such as GcvB or CsrB at a single-cell level. In addition, we confirmed previously described phenotypic heterogeneity in Salmonella in regard to expression of pathogenicity-associated genes. Overall, the low percentage of cell loss and high gene detection limit makes the improved MATQ-seq protocol particularly well suited for studies with limited input material, such as analysis of small bacterial populations in host niches or intracellular bacteria. IMPORTANCE Gene expression heterogeneity among isogenic bacteria is linked to clinically relevant scenarios, like biofilm formation and antibiotic tolerance. The recent development of bacterial single-cell RNA sequencing (scRNA-seq) enables the study of cell-to-cell variability in bacterial populations and the mechanisms underlying these phenomena. Here, we report a scRNA-seq workflow based on MATQ-seq with increased robustness, reduced cell loss, and improved transcript capture rate and gene coverage. Use of a more efficient reverse transcriptase and the integration of an rRNA depletion step, which can be adapted to other bacterial single-cell workflows, was instrumental for these improvements. Applying the protocol to the foodborne pathogen Salmonella, we confirmed transcriptional heterogeneity across and within different growth phases and demonstrated that our workflow captures small regulatory RNAs at a single-cell level. Due to low cell loss and high transcript capture rates, this protocol is uniquely suited for experimental settings in which the starting material is limited, such as infected tissues.

While an increasing number of conserved small regulatory RNAs (sRNAs) are known to function in general bacterial physiology, the roles and modes of action of sRNAs from horizontally acquired genomic regions remain little understood. The IsrK sRNA of Gifsy-1 prophage of Salmonella belongs to the latter class. This regulatory RNA exists in two isoforms. The first forms, when a portion of transcripts originating from isrK promoter reads-through the IsrK transcription-terminator producing a translationally inactive mRNA target. Acting in trans, the second isoform, short IsrK RNA, binds the inactive transcript rendering it translationally active. By switching on translation of the first isoform, short IsrK indirectly activates the production of AntQ, an antiterminator protein located upstream of isrK. Expression of antQ globally interferes with transcription termination resulting in bacterial growth arrest and ultimately cell death. Escherichia coli and Salmonella cells expressing AntQ display condensed chromatin morphology and localization of UvrD to the nucleoid. The toxic phenotype of AntQ can be rescued by co-expression of the transcription termination factor, Rho, or RNase H, which protects genomic DNA from breaks by resolving R-loops. We propose that AntQ causes conflicts between transcription and replication machineries and thus promotes DNA damage. The isrK locus represents a unique example of an island-encoded sRNA that exerts a highly complex regulatory mechanism to tune the expression of a toxic protein.

RNA-binding proteins (RBPs) play important roles in bacterial gene expression and physiology but their true number and functional scope remain little understood even in model microbes. To advance global RBP discovery in bacteria, we here establish glycerol gradient sedimentation with RNase treatment and mass spectrometry (GradR). Applied to Salmonella enterica, GradR confirms many known RBPs such as CsrA, Hfq, and ProQ by their RNase-sensitive sedimentation profiles, and discovers the FopA protein as a new member of the emerging family of FinO/ProQ-like RBPs. FopA, encoded on resistance plasmid pCol1B9, primarily targets a small RNA associated with plasmid replication. The target suite of FopA dramatically differs from the related global RBP ProQ, revealing context-dependent selective RNA recognition by FinO-domain RBPs. Numerous other unexpected RNase-induced changes in gradient profiles suggest that cellular RNA helps to organize macromolecular complexes in bacteria. By enabling poly(A)-independent generic RBP discovery, GradR provides an important element in the quest to build a comprehensive catalog of microbial RBPs.

CRISPR/Cas systems constitute a widespread class of immunity systems that protect bacteria and archaea against phages and plasmids, and commonly use repeat/spacer-derived short crRNAs to silence foreign nucleic acids in a sequence-specific manner. Although the maturation of crRNAs represents a key event in CRISPR activation, the responsible endoribonucleases (CasE, Cas6, Csy4) are missing in many CRISPR/Cas subtypes. Here, differential RNA sequencing of the human pathogen Streptococcus pyogenes uncovered tracrRNA, a trans-encoded small RNA with 24-nucleotide complementarity to the repeat regions of crRNA precursor transcripts. We show that tracrRNA directs the maturation of crRNAs by the activities of the widely conserved endogenous RNase III and the CRISPR-associated Csn1 protein; all these components are essential to protect S. pyogenes against prophage-derived DNA. Our study reveals a novel pathway of small guide RNA maturation and the first example of a host factor (RNase III) required for bacterial RNA-mediated immunity against invaders.

Successful elimination of bacteria in phagocytes occurs in the phago-lysosomal system, but also depends on mitochondrial pathways. Yet, how these two organelle systems communicate is largely unknown. Here we identify the lysosomal biogenesis factor transcription factor EB (TFEB) as regulator for phago-lysosome-mitochondria crosstalk in macrophages. By combining cellular imaging and metabolic profiling, we find that TFEB activation, in response to bacterial stimuli, promotes the transcription of aconitate decarboxylase (Acod1, Irg1) and synthesis of its product itaconate, a mitochondrial metabolite with antimicrobial activity. Activation of the TFEB-Irg1-itaconate signalling axis reduces the survival of the intravacuolar pathogen Salmonella enterica serovar Typhimurium. TFEB-driven itaconate is subsequently transferred via the Irg1-Rab32-BLOC3 system into the Salmonella-containing vacuole, thereby exposing the pathogen to elevated itaconate levels. By activating itaconate production, TFEB selectively restricts proliferating Salmonella, a bacterial subpopulation that normally escapes macrophage control, which contrasts TFEB's role in autophagy-mediated pathogen degradation. Together, our data define a TFEB-driven metabolic pathway between phago-lysosomes and mitochondria that restrains Salmonella Typhimurium burden in macrophages in vitro and in vivo.

While mucosal inflammation is a major source of stress during enteropathogen infection, it remains to be fully elucidated how the host benefits from this environment to clear the pathogen. Here, we show that host stress induced by different stimuli mimicking inflammatory conditions strongly reduces the binding of Shigella flexneri to epithelial cells. Mechanistically, stress activates acid sphingomyelinase leading to host membrane remodeling. Consequently, knockdown or pharmacological inhibition of the acid sphingomyelinase blunts the stress-dependent inhibition of Shigella binding to host cells. Interestingly, stress caused by intracellular Shigella replication also results in remodeling of the host cell membrane, in vitro and in vivo, which precludes re-infection by this and other non-motile pathogens. In contrast, Salmonella Typhimurium overcomes the shortage of permissive entry sites by gathering effectively at the remaining platforms through its flagellar motility. Overall, our findings reveal host membrane remodeling as a novel stress-responsive cell-autonomous defense mechanism that protects epithelial cells from infection by non-motile bacterial pathogens.

The transcriptome is a powerful proxy for the physiological state of a cell, healthy or diseased. As a result, transcriptome analysis has become a key tool in understanding the molecular changes that accompany bacterial infections of eukaryotic cells. Until recently, such transcriptomic studies have been technically limited to analyzing mRNA expression changes in either the bacterial pathogen or the infected eukaryotic host cell. However, the increasing sensitivity of high-throughput RNA sequencing now enables "dual RNA-seq" studies, simultaneously capturing all classes of coding and noncoding transcripts in both the pathogen and the host. In the five years since the concept of dual RNA-seq was introduced, the technique has been applied to a range of infection models. This has not only led to a better understanding of the physiological changes in pathogen and host during the course of an infection but has also revealed hidden molecular phenotypes of virulence-associated small noncoding RNAs that were not visible in standard infection assays. Here, we use the knowledge gained from these recent studies to suggest experimental and computational guidelines for the design of future dual RNA-seq studies. We conclude this review by discussing prospective applications of the technique.

The envelope of Gram-negative bacteria is a vital barrier that must balance protection and nutrient uptake. Small RNAs are crucial regulators of the envelope composition and function. Here, using RIL-seq to capture the Hfq-mediated RNA-RNA interactome in Salmonella enterica, we discover envelope-related riboregulators, including OppX. We show that OppX acts as an RNA sponge of MicF sRNA, a prototypical porin repressor. OppX originates from the 5' UTR of oppABCDF, encoding the major inner-membrane oligopeptide transporter, and sequesters MicF's seed region to derepress the synthesis of the porin OmpF. Intriguingly, OppX operates as a true sponge, storing MicF in an inactive complex without affecting its levels or stability. Conservation of the opp-OppX-MicF-ompF axis in related bacteria suggests that it serves an important mechanism, adjusting envelope porosity to specific transport capacity. These data also highlight the resource value of this Salmonella RNA interactome, which will aid in unraveling RNA-centric regulation in enteric pathogens.

New methods for the global identification of RNA-protein interactions have led to greater recognition of the abundance and importance of RNA-binding proteins (RBPs) in bacteria. Here, we expand this tool kit by developing SEC-seq, a method based on a similar concept as the established Grad-seq approach. In Grad-seq, cellular RNA and protein complexes of a bacterium of interest are separated in a glycerol gradient, followed by high-throughput RNA-sequencing and mass spectrometry analyses of individual gradient fractions. New RNA-protein complexes are predicted based on the similarity of their elution profiles. In SEC-seq, we have replaced the glycerol gradient with separation by size exclusion chromatography, which shortens operation times and offers greater potential for automation. Applying SEC-seq to Escherichia coli, we find that the method provides a higher resolution than Grad-seq in the lower molecular weight range up to ~500 kDa. This is illustrated by the ability of SEC-seq to resolve two distinct, but similarly sized complexes of the global translational repressor CsrA with either of its antagonistic small RNAs, CsrB and CsrC. We also characterized changes in the SEC-seq profiles of the small RNA MicA upon deletion of its RNA chaperones Hfq and ProQ and investigated the redistribution of these two proteins upon RNase treatment. Overall, we demonstrate that SEC-seq is a tractable and reproducible method for the global profiling of bacterial RNA-protein complexes that offers the potential to discover yet-unrecognized associations between bacterial RNAs and proteins.

The protein ProQ has recently been identified as a global small noncoding RNA-binding protein in Salmonella, and a similar role is anticipated for its numerous homologs in divergent bacterial species. We report the solution structure of Escherichia coli ProQ, revealing an N-terminal FinO-like domain, a C-terminal domain that unexpectedly has a Tudor domain fold commonly found in eukaryotes, and an elongated bridging intradomain linker that is flexible but nonetheless incompressible. Structure-based sequence analysis suggests that the Tudor domain was acquired through horizontal gene transfer and gene fusion to the ancestral FinO-like domain. Through a combination of biochemical and biophysical approaches, we have mapped putative RNA-binding surfaces on all three domains of ProQ and modeled the protein's conformation in the apo and RNA-bound forms. Taken together, these data suggest how the FinO, Tudor, and linker domains of ProQ cooperate to recognize complex RNA structures and serve to promote RNA-mediated regulation.

Antisense oligomers (ASOs), such as peptide nucleic acids (PNAs), designed to inhibit the translation of essential bacterial genes, have emerged as attractive sequence- and species-specific programmable RNA antibiotics. Yet, potential drawbacks include unwanted side effects caused by their binding to transcripts other than the intended target. To facilitate the design of PNAs with minimal off-target effects, we developed MASON (make antisense oligomers now), a web server for the design of PNAs that target bacterial mRNAs. MASON generates PNA sequences complementary to the translational start site of a bacterial gene of interest and reports critical sequence attributes and potential off-target sites. We based MASON's off-target predictions on experiments in which we treated Salmonella enterica serovar Typhimurium with a series of 10-mer PNAs derived from a PNA targeting the essential gene acpP but carrying two serial mismatches. Growth inhibition and RNA-sequencing (RNA-seq) data revealed that PNAs with terminal mismatches are still able to target acpP, suggesting wider off-target effects than anticipated. Comparison of these results to an RNA-seq data set from uropathogenic Escherichia coli (UPEC) treated with eleven different PNAs confirmed that our findings are not unique to Salmonella We believe that MASON's off-target assessment will improve the design of specific PNAs and other ASOs.

Small non-coding RNAs (sRNAs) are an emerging class of post-transcriptional regulators of bacterial gene expression. To study sRNAs and their potential protein interaction partners, it is desirable to purify sRNAs from cells in their native form. Here, we used RNA-based affinity chromatography to purify sRNAs following their expression as aptamer-tagged variants in vivo. To this end, we developed a family of plasmids to express sRNAs with any of three widely used aptamer sequences (MS2, boxB, eIF4A), and systematically tested how the aptamer tagging impacted on intracellular accumulation and target regulation of the Salmonella GcvB, InvR or RybB sRNAs. In addition, we successfully tagged the chromosomal rybB gene with MS2 to observe that RybB-MS2 is fully functional as an envelope stress-induced repressor of ompN mRNA following induction of sigmaE. We further demonstrate that the common sRNA-binding protein, Hfq, co-purifies with MS2-tagged sRNAs of Salmonella. The presented affinity purification strategy may facilitate the isolation of in vivo assembled sRNA-protein complexes in a wide range of bacteria.

The RNA chaperone, Hfq, plays a diverse role in bacterial physiology beyond its original role as a host factor required for replication of Qbeta RNA bacteriophage. In this study, we show that Hfq is involved in the expression and secretion of virulence factors in the facultative intracellular pathogen, Salmonella typhimurium. A Salmonella hfq deletion strain is highly attenuated in mice after both oral and intraperitoneal infection, and shows a severe defect in invasion of epithelial cells and a growth defect in both epithelial cells and macrophages in vitro. Surprisingly, we find that these phenotypes are largely independent of the previously reported requirement of Hfq for expression of the stationary phase sigma factor, RpoS. Our results implicate Hfq as a key regulator of multiple aspects of virulence including regulation of motility and outer membrane protein (OmpD) expression in addition to invasion and intracellular growth. These pleiotropic effects are suggested to involve a network of regulatory small non-coding RNAs, placing Hfq at the centre of post-transcriptional regulation of virulence gene expression in Salmonella. In addition, the hfq mutation appears to cause a chronic activation of the RpoE-mediated envelope stress response which is likely due to a misregulation of membrane protein expression.

Small proteins are an emerging class of gene products with diverse roles in bacterial physiology. However, a full understanding of their importance has been hampered by insufficient genome annotations and a lack of comprehensive characterization in microbes other than Escherichia coli. We have taken an integrative approach to accelerate the discovery of small proteins and their putative virulence-associated functions in Salmonella Typhimurium. We merged the annotated small proteome of Salmonella with new small proteins predicted with in silico and experimental approaches. We then exploited existing and newly generated global datasets that provide information on small open reading frame expression during infection of epithelial cells (dual RNA-seq), contribution to bacterial fitness inside macrophages (Transposon-directed insertion sequencing), and potential engagement in molecular interactions (Grad-seq). This integrative approach suggested a new role for the small protein MgrB beyond its known function in regulating PhoQ. We demonstrate a virulence and motility defect of a Salmonella ΔmgrB mutant and reveal an effect of MgrB in regulating the Salmonella transcriptome and proteome under infection-relevant conditions. Our study highlights the power of interpreting available 'omics' datasets with a focus on small proteins, and may serve as a blueprint for a data integration-based survey of small proteins in diverse bacteria.

FinO domain proteins such as ProQ of the model pathogen Salmonella enterica have emerged as a new class of major RNA-binding proteins in bacteria. ProQ has been shown to target hundreds of transcripts, including mRNAs from many virulence regions, but its role, if any, in bacterial pathogenesis has not been studied. Here, using a Dual RNA-seq approach to profile ProQ-dependent gene expression changes as Salmonella infects human cells, we reveal dysregulation of bacterial motility, chemotaxis, and virulence genes which is accompanied by altered MAPK (mitogen-activated protein kinase) signaling in the host. Comparison with the other major RNA chaperone in Salmonella, Hfq, reinforces the notion that these two global RNA-binding proteins work in parallel to ensure full virulence. Of newly discovered infection-associated ProQ-bound small noncoding RNAs (sRNAs), we show that the 3'UTR-derived sRNA STnc540 is capable of repressing an infection-induced magnesium transporter mRNA in a ProQ-dependent manner. Together, this comprehensive study uncovers the relevance of ProQ for Salmonella pathogenesis and highlights the importance of RNA-binding proteins in regulating bacterial virulence programs.IMPORTANCE The protein ProQ has recently been discovered as the centerpiece of a previously overlooked "third domain" of small RNA-mediated control of gene expression in bacteria. As in vitro work continues to reveal molecular mechanisms, it is also important to understand how ProQ affects the life cycle of bacterial pathogens as these pathogens infect eukaryotic cells. Here, we have determined how ProQ shapes Salmonella virulence and how the activities of this RNA-binding protein compare with those of Hfq, another central protein in RNA-based gene regulation in this and other bacteria. To this end, we apply global transcriptomics of pathogen and host cells during infection. In doing so, we reveal ProQ-dependent transcript changes in key virulence and host immune pathways. Moreover, we differentiate the roles of ProQ from those of Hfq during infection, for both coding and noncoding transcripts, and provide an important resource for those interested in ProQ-dependent small RNAs in enteric bacteria.