Publications

Quorum sensing (QS) orchestrates collective responses in bacteria. Yet, how individuality arises and functions within a system designed to unify behavior remains poorly understood. Here, we use imaging-transcriptomics to comprehensively profile the QS transcriptional response of the opportunistic pathogen Pseudomonas aeruginosa at single-cell resolution. We observe heterogeneity across all stages of QS, from signal production to public goods expression. While most cells cooperate, the extent of public goods contribution varies substantially, primarily due to transcriptional noise. Additionally, cells expressing high levels of a given exoproduct were often enriched for multiple other public goods. Expression of QS signal synthases, particularly in the Las and PQS systems, exhibits extreme cell-cell variability indicative of regulation, and exceeding levels observed in known differentiation into hypervirulent or motile subpopulations. These patterns are dynamic, peaking during early signaling, suggesting that a subset of highly committed cells may influence key decision points in the QS response. We find that cellular memory carried from prior growth cycles shapes QS dynamics via signaling-primed cells but does not influence de novo formation of hyper-signaling cells. Differentiation into hyper-signaling subpopulations is also robust to exogenous autoinducers across a wide concentration range, suggesting it is likely mediated by an internal mechanism. Despite variation in QS network outputs between strains, hyper-signaling subpopulations consistently emerge across diverse isolates, pointing to an evolutionarily conserved strategy. Our findings recast QS as a system that embeds a deliberately heterogeneous entry point within an otherwise synchronizing program, offering new insight into how bacterial populations balance cooperation and conflict.Competing Interest StatementThe authors have declared no competing interest.European Research Council, https://ror.org/0472cxd90, 101117863

Bacterial anti-phage systems are frequently clustered in microbial genomes, forming defense islands. This property enabled the recent discovery of multiple defense systems based on their genomic co-localization with known systems, but the full arsenal of anti-phage mechanisms remains unknown. We report the discovery of 21 defense systems that protect bacteria from phages, based on computational genomic analyses and phage-infection experiments. We identified multiple systems with domains involved in eukaryotic antiviral immunity, including those homologous to the ubiquitin-like ISG15 protein, dynamin-like domains, and SEFIR domains, and show their participation in bacterial defenses. Additional systems include domains predicted to manipulate DNA and RNA molecules, alongside toxin-antitoxin systems shown here to function in anti-phage defense. These systems are widely distributed in microbial genomes, and in some bacteria, they form a considerable fraction of the immune arsenal. Our data substantially expand the inventory of defense systems utilized by bacteria to counteract phage infection.

Capturing the heterogeneous phenotypes of microbial populations at relevant spatiotemporal scales is highly challenging. Here, we present par-seqFISH (parallel sequential fluorescence in situ hybridization), a transcriptome-imaging approach that records gene expression and spatial context within microscale assemblies at a single-cell and molecule resolution. We applied this approach to the opportunistic pathogen Pseudomonas aeruginosa, analyzing about 600,000 individuals across dozens of conditions in planktonic and biofilm cultures. We identified numerous metabolic- and virulence-related transcriptional states that emerged dynamically during planktonic growth, as well as highly spatially resolved metabolic heterogeneity in sessile populations. Our data reveal that distinct physiological states can coexist within the same biofilm just several micrometers away, underscoring the importance of the microenvironment. Our results illustrate the complex dynamics of microbial populations and present a new way of studying them at high resolution.

Phenazines are natural bacterial antibiotics that can protect crops from disease. However, for most crops it is unknown which producers and specific phenazines are ecologically relevant, and whether phenazine biodegradation can counter their effects. To better understand their ecology, we developed and environmentally-validated a quantitative metagenomic approach to mine for phenazine biosynthesis and biodegradation genes, applying it to >800 soil and plant-associated shotgun-metagenomes. We discover novel producer-crop associations and demonstrate that phenazine biosynthesis is prevalent across habitats and preferentially enriched in rhizospheres, whereas biodegrading bacteria are rare. We validate an association between maize and Dyella japonica, a putative producer abundant in crop microbiomes. D. japonica upregulates phenazine biosynthesis during phosphate limitation and robustly colonizes maize seedling roots. This work provides a global picture of phenazines in natural environments and highlights plant-microbe associations of agricultural potential. Our metagenomic approach may be extended to other metabolites and functional traits in diverse ecosystems.

N-Acyl homoserine lactone (AHL) quorum sensing (QS) controls expression of over 200 genes in Pseudomonas aeruginosa. There are two AHL regulatory systems: The LasR-LasI circuit and the RhlR-RhlI system. We mapped transcription termination sites affected by AHL QS in P. aeruginosa, and in doing so we identified AHL-regulated small RNAs (sRNAs). Of interest, we noted that one particular sRNA was located within the rhlI locus. We found that rhlI, which encodes the enzyme that produces the AHL N-butanoyl-homoserine lactone (C4-HSL), is controlled by a 5 = untranslated region (UTR)-derived sRNA we name RhlS. We also identified an antisense RNA encoded opposite the beginning of the rhlI open reading frame, which we name asRhlS. RhlS accumulates as wild-type cells enter stationary phase and is required for the production of normal levels of C4-HSL through activation of rhlI translation. RhlS also directly posttranscriptionally regulates at least one other unlinked gene, fpvA. The asRhlS appears to be expressed at maximal levels during logarithmic growth, and we suggest RhlS may act antagonistically to the asRhlS to regulate rhlI translation. The rhlI-encoded sRNAs represent a novel aspect of RNAmediated tuning of P. aeruginosa QS. IMPORTANCE The opportunistic human pathogen Pseudomonas aeruginosa possesses multiple quorum sensing systems that regulate and coordinate production of virulence factors and adaptation to different environments. Despite extensive research, the regulatory elements that play a role in this complex network are still not fully understood. By using several RNA sequencing techniques, we were able to identify a small regulatory RNA we named RhlS. RhlS increases translation of RhlI, a key enzyme in the quorum sensing pathway, and represses the fpvA mRNA encoding one of the siderophore pyoverdine receptors. Our results highlight a new regulatory layer of P. aeruginosa quorum sensing and contribute to the growing understanding of the role regulatory RNAs play in bacterial physiology.

A signal transducer and activator of transcription-1-deficient patient presented with prolonged fever, cachexia, anemia, hypoalbuminemia and finally relapsing debilitating mycobacterial osteomyelitis while receiving a previously effective antimycobacterial treatment. Progression despite rigorous workup and multiple antibiotics prompted shotgun metagenomics revealing adenovirus in liver samples. Brincidofovir led to a complete, sustained clinical recovery, including osteomyelitis, probably attributed to reversal of adenovirus-induced immune dysregulation.

To overcome the action of antibiotics, bacteria have evolved a variety of different strategies, such as drug modification, target mutation, and efflux pumps. Recently, we performed a genome-wide analysis of Listeria monocytogenes gene expression after growth in the presence of antibiotics, identifying genes that are up-regulated upon antibiotic treatment. One of them, lmo0762, is a homologof hflX, which encodes a heat shock protein that rescues stalled ribosomes by separating their two subunits. To our knowledge, ribosome splitting has never been described as an antibiotic resistance mechanism. We thus investigated the role of lmo0762 in antibiotic resistance. First, we demonstrated that lmo0762 is an antibiotic resistance gene that confers protection against lincomycin and erythromycin, and that we renamed hflXr (hflX resistance). We show that hflXr expression is regulated by a transcription attenuation mechanism relying on the presence of alternative RNA structures and a small ORF encoding a 14 amino acid peptide containing the RLR motif, characteristic of macrolide resistance genes. We also provide evidence that HflXr is involved in ribosome recycling in presence of antibiotics. Interestingly, L. monocytogenes possesses another copy of hflX, lmo1296, that is not involved in antibiotic resistance. Phylogenetic analysis shows several events of hflXr duplication in prokaryotes and widespread presence of hflXr in Firmicutes. Overall, this study reveals the Listeria hflXr as the founding member of a family of antibiotic resistance genes. The resistance conferred by this gene is probably of importance in the environment and within microbial communities.

Prokaryotic genomes encode a plethora of small noncoding RNAs (ncRNAs) that fine-tune the expression of specific genes. The vast majority of known bacterial ncRNAs are encoded from within intergenic regions, where their expression is controlled by promoter and terminator elements, similarly to protein-coding genes. In addition, recent studies have shown that functional ncRNAs can also be derived from gene 3' untranslated regions (3' UTRs) via an alternative biogenesis pathway, in which the ncRNA segment is separated from the mRNA via RNase cleavage. Here, we report the detection of a large set of decay-generated noncoding RNAs (decRNAs), many of which are completely embedded within protein-coding mRNA regions rather than in the UTRs. We show that these decRNAs are "carved out" of the mRNA through the action of RNase E and that they are predicted to fold into highly stable RNA structures, similar to those of known ncRNAs. A subset of these decRNAs is predicted to interact with Hfq or ProQ or both, which act as ncRNA chaperones, and some decRNAs display evolutionarily conserved sequences and conserved expression patterns between different species. These results suggest that mRNA protein-coding regions may harbor a large set of potentially functional small RNAs

Transcription termination in bacteria can occur either via Rho-dependent or independent (intrinsic) mechanisms. Intrinsic terminators are composed of a stem-loop RNA structure followed by a uridine stretch and are known to terminate in a precise manner. In contrast, Rho-dependent terminators have more loosely defined characteristics and are thought to terminate in a diffuse manner. While transcripts ending in an intrinsic terminator are protected from 3'-5' exonuclease digestion due to the stem-loop structure of the terminator, it remains unclear what protects Rho-dependent transcripts from being degraded. In this study, we mapped the exact steady-state RNA 3' ends of hundreds of Escherichia coli genes terminated either by Rho-dependent or independent mechanisms. We found that transcripts generated from Rho-dependent termination have precise 3'-ends at steady state. These termini were localized immediately downstream of energetically stable stem-loop structures, which were not followed by uridine rich sequences. We provide evidence that these structures protect Rho-dependent transcripts from 3'-5' exonucleases such as PNPase and RNase II, and present data localizing the Rho-utilization (rut) sites immediately downstream of these protective structures. This study represents the first extensive in-vivo map of exact RNA 3'-ends of Rho-dependent transcripts in E. coli.

Bacterial operons synchronize the expression of multiple genes by placing them under the control of a shared promoter. It was previously shown that polycistronic transcripts can undergo differential RNA decay, leaving some genes within the polycistron more stable than others, but the extent of regulation by differential mRNA decay or its evolutionary conservation remains unknown. Here, we find that a substantial fraction of E. coli genes display non-uniform mRNA stoichiometries despite being coded from the same operon. We further show that these altered operon stoichiometries are shaped post-transcriptionally by differential mRNA decay, which is regulated by RNA structures that protect specific regions in the transcript from degradation. These protective RNA structures are generally coded within the protein-coding regions of the regulated genes and are frequently evolutionarily conserved. Furthermore, we provide evidence that differences in ribosome densities across polycistronic transcript segments, together with the conserved structural RNA elements, play a major role in the differential decay process. Our results highlight a major role for differential mRNA decay in shaping bacterial transcriptomes.

Antibiotic resistance genes are commonly regulated by sophisticated mechanisms that activate gene expression in response to antibiotic exposure. Growing evidence suggest that cis-acting non-coding RNAs play a major role in regulating the expression of many resistance genes, specifically those which counteract the effects of translation-inhibiting antibiotics. These ncRNAs reside in the 5'UTR of the regulated gene, and sense the presence of the antibiotics by recruiting translating ribosomes onto short upstream open reading frames (uORFs) embedded in the ncRNA. In the presence of translation-inhibiting antibiotics ribosomes arrest over the uORF, altering the RNA structure of the regulator and switching the expression of the resistance gene to 'ON'. The specificity of these riboregulators is tuned to sense-specific classes of antibiotics based on the length and composition of the respective uORF. Here we review recent work describing new types of antibiotic-sensing RNA-based regulators and elucidating the molecular mechanisms by which they function to control antibiotic resistance in bacteria.

A common strategy for regulation of gene expression in bacteria is conditional transcription termination. This strategy is frequently employed by 5'UTR cis-acting RNA elements (riboregulators), including riboswitches and attenuators. Such riboregulators can assume two mutually exclusive RNA structures, one of which forms a transcriptional terminator and results in premature termination, and the other forms an antiterminator that allows read-through into the coding sequence to produce a full-length mRNA. We developed a machine-learning based approach, which, given a 5'UTR of a gene, predicts whether it can form the two alternative structures typical to riboregulators employing conditional termination. Using a large positive training set of riboregulators derived from 89 human microbiome bacteria, we show high specificity and sensitivity for our classifier. We further show that our approach allows the discovery of previously unidentified riboregulators, as exemplified by the detection of new LeuA leaders and T-boxes in Streptococci. Finally, we developed PASIFIC (www.weizmann.ac.il/molgen/Sorek/PASIFIC/), an online web-server that, given a user-provided 5'UTR sequence, predicts whether this sequence can adopt two alternative structures conforming with the conditional termination paradigm. This webserver is expected to assist in the identification of new riboswitches and attenuators in the bacterial pan-genome.

Transcription termination sets the 3' end boundaries of RNAs and plays key roles in gene regulation. Although termination has been well studied in bacteria, the signals that mediate termination in archaea remain poorly understood. Here, we applied term-seq to comprehensively map RNA 3' termini, with single-base precision, in two phylogenetically distant archaea: Methanosarcina mazei and Sulfolobus acidocaldarius. Comparison of RNA 3' ends across hundreds of genes revealed the sequence composition of transcriptional terminators in each organism, highlighting both common and divergent characteristics between the different archaeal phyla. We find that, in contrast to bacteria, a considerable portion of archaeal genes are controlled by multiple consecutive terminators, generating several alternative 3' untranslated region isoforms for >30% of the genes. These alternative isoforms often present marked length differences, implying that archaea can employ regulation via alternative 3' untranslated regions, similar to eukaryotes. Although most of the terminators are intergenic, we discover numerous cases in which termination of one gene occurs within the coding region of a downstream gene, implying that leaky termination may tune inter-transcript stoichiometry in multi-gene operons. These results provide the first high-throughput maps of transcriptional terminators in archaea and point to an evolutionary path linking bacterial and eukaryal non-coding regulatory strategies.

Riboswitches and attenuators are cis-regulatory RNA elements, most of which control bacterial gene expression via metabolite-mediated, premature transcription termination. We developed an unbiased experimental approach for genome-wide discovery of such ribo-regulators in bacteria. We also devised an experimental platform that quantitatively measures the in vivo activity of all such regulators in parallel and enables rapid screening for ribo-regulators that respond to metabolites of choice. Using this approach, we detected numerous antibioticresponsive ribo-regulators that control antibiotic resistance genes in pathogens and in the human microbiome. Studying one such regulator in Listeria monocytogenes revealed an attenuation mechanism mediated by antibiotic-stalled ribosomes. Our results expose broad roles for conditional termination in regulating antibiotic resistance and provide a tool for discovering riboswitches and attenuators that respond to previously unknown ligands.

Whole-transcriptome sequencing studies from recent years revealed an unexpected complexity in transcriptomes of bacteria and archaea, including abundant non-coding RNAs, cis-antisense transcription and regulatory untranslated regions (UTRs). Understanding the functional relevance of the plethora of non-coding RNAs in a given organism is challenging, especially since some of these RNAs were attributed to 'transcriptional noise'. To allow the search for conserved transcriptomic elements we produced comparative transcriptome maps for multiple species across the microbial tree of life. These transcriptome maps are detailed in annotations, comparable by gene families, and BLAST-searchable by user provided sequences. Our transcriptome collection includes 18 model organisms spanning 10 phyla/subphyla of bacteria and archaea that were sequenced using standardized RNA-seq methods. The utility of the comparative approach, as implemented in our web server, is demonstrated by highlighting genes with exceptionally long 5UTRs across species, which correspond to many known riboswitches and further suggest novel putative regulatory elements. Our study provides a standardized reference transcriptome to major clinically and environmentally important microbial phyla.

Riboswitches are ligand-binding elements contained within the 5 untranslated regions of bacterial transcripts, which generally regulate expression of downstream open reading frames. Here, we show that in Listeria monocytogenes, a riboswitch that binds vitamin B₁₂ controls expression of a noncoding regulatory RNA, Rli55. Rli55, in turn, controls expression of the eut genes, whose products enable ethanolamine utilization and require B₁₂ as a cofactor. Defects in ethanolamine utilization, or in its regulation by Rli55, significantly attenuate Listeria virulence in mice. Rli55 functions by sequestering the two-component response regulator EutV by means of a EutV-binding site contained within the RNA. Thus, Rli55 is a riboswitch-regulated member of the small group of regulatory RNAs that function by sequestering a protein and reveals a distinctive mechanism of signal integration in bacterial gene regulation.

The Gram-negative bacterium Pseudomonas syringae pv. tomato (Pst) is the causal agent of speck disease in tomato. Pst pathogenicity depends on a type III secretion system that delivers effector proteins into host cells, where they promote disease by manipulating processes to the advantage of the pathogen. Previous studies identified seven Pst effectors that inhibit growth when expressed in yeast under normal growth conditions, suggesting that they interfere with cellular processes conserved in yeast and plants. We hypothesized that effectors also target conserved cellular processes that are required for yeast growth only under stress conditions. We therefore examined phenotypes induced by expression of Pst effectors in yeast grown in the presence of various stressors. Out of 29 effectors tested, five (HopX1, HopG1, HopT1-1, HopH1 and AvrPtoB) displayed growth inhibition phenotypes only in combination with stress conditions. Viability assays revealed that the HopX1 effector caused loss of cell viability under prolonged osmotic stress. Using transcription reporters, we found that HopX1 attenuated the activation of the high osmolarity glycerol (HOG) mitogen-activated protein kinase (MAPK) pathway, which is responsible for yeast survival under osmotic stress, while other MAPK pathways were mildly affected by HopX1. Interestingly, HopX1-mediated phenotypes in yeast were dependent on the putative transglutaminase catalytic triad of the effector. This study enlarges the pool of phenotypes available for the functional analysis of Pst type III effectors in yeast, and exemplifies how analysis of phenotypes detected in yeast under stress conditions can lead to the identification of eukaryotic cellular processes affected by bacterial effectors.

The gram-negative bacterium Xanthomonas campestris pv. vesicatoria is the causal agent of spot disease in tomato and pepper. X. campestris pv. vesicatoria pathogenicity depends on a type III secretion system delivering effector proteins into the host cells. We hypothesized that some X. campestris pv. vesicatoria effectors target conserved eukaryotic cellular processes and examined phenotypes induced by their expression in yeast. Out of 21 effectors tested, 14 inhibited yeast growth in normal or stress conditions. Viability assay revealed that XopB and XopF2 attenuated cell proliferation, while AvrRxo1, XopX, and XopE1 were cytotoxic. Inspection of morphological features and DNA content of yeast cells indicated that cytotoxicity caused by XopX and AvrRxo1 was associated with cell-cycle arrest at G0/1. Interestingly, XopB, XopE1, XopF2, XopX, and AvrRxo1 that inhibited growth in yeast also caused phenotypes, such as chlorosis and cell death, when expressed in either host or nonhost plants. Finally, the ability of several effectors to cause phenotypes in yeast and plants was dependent on their putative catalytic residues or localization motifs. This study supports the use of yeast as a heterologous system for functional analysis of X. campestris pv. vesicatoria type III effectors, and sets the stage for identification of their eukaryotic molecular targets and modes of action.