Global RNA interactome of Salmonella discovers a 5' UTR sponge for the MicF small RNA that connects membrane permeability to transport capacity.

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.

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.

Global Mapping of Small RNA-Target Interactions in Bacteria.

Small RNAs (sRNAs) associated with the RNA chaperon protein Hfq are key posttranscriptional regulators of gene expression in bacteria. Deciphering the sRNA-target interactome is an essential step toward understanding the roles of sRNAs in the cellular networks. We developed a broadly applicable methodology termed RIL-seq (RNA interaction by ligation and sequencing), which integrates experimental and computational tools for in vivo transcriptome-wide identification of interactions involving Hfq-associated sRNAs. By applying this methodology to Escherichia coli we discovered an extensive network of interactions involving RNA pairs showing sequence complementarity. We expand the ensemble of targets for known sRNAs, uncover additional Hfq-bound sRNAs encoded in various genomic regions along with their trans encoded targets, and provide insights into binding and possible cycling of RNAs on Hfq. Comparison of the sRNA interactome under various conditions has revealed changes in the sRNA repertoire as well as substantial re-wiring of the network between conditions.

Common Adaptive Strategies Underlie Within-Host Evolution of Bacterial Pathogens.

Within-host adaptation is a hallmark of chronic bacterial infections, involving substantial genomic changes. Recent large-scale genomic data from prolonged infections allow the examination of adaptive strategies employed by different pathogens and open the door to investigate whether they converge toward similar strategies. Here, we compiled extensive data of whole-genome sequences of bacterial isolates belonging to miscellaneous species sampled at sequential time points during clinical infections. Analysis of these data revealed that different species share some common adaptive strategies, achieved by mutating various genes. Although the same genes were often mutated in several strains within a species, different genes related to the same pathway, structure, or function were changed in other species utilizing the same adaptive strategy (e.g., mutating flagellar genes). Strategies exploited by various bacterial species were often predicted to be driven by the host immune system, a powerful selective pressure that is not species specific. Remarkably, we find adaptive strategies identified previously within single species to be ubiquitous. Two striking examples are shifts from siderophore-based to heme-based iron scavenging (previously shown for Pseudomonas aeruginosa) and changes in glycerol-phosphate metabolism (previously shown to decrease sensitivity to antibiotics in Mycobacterium tuberculosis). Virulence factors were often adaptively affected in different species, indicating shifts from acute to chronic virulence and virulence attenuation during infection. Our study presents a global view on common within-host adaptive strategies employed by different bacterial species and provides a rich resource for further studying these processes.

Assessing the functional association of intronic miRNAs with their host genes.

In human, nearly half of the known microRNAs (miRNAs) are encoded within the introns of protein-coding genes. The embedment of these miRNA genes within the sequences of protein-coding genes alludes to a possible functional relationship between intronic miRNAs and their hosting genes. Several studies, using predicted targets, suggested that intronic miRNAs influence their hosts' function either antagonistically or synergistically. New experimental data of miRNA expression patterns and targets enable exploring this putative association by relying on actual data rather than on predictions. Here, our analysis based on currently available experimental data implies that the potential functional association between intronic miRNAs and their hosting genes is limited. For host-miRNA examples where functional associations were detected, it was manifested by either autoregulation, common targets of the miRNA and hosting gene, or through the targeting of transcripts participating in pathways in which the host gene is involved. This low prevalence of functional association is consistent with our observation that many intronic miRNAs have independent transcription start sites and are not coexpressed with the hosting gene. Yet, the intronic miRNAs that do show functional association with their hosts were found to be more evolutionarily conserved compared to other intronic miRNAs. This might suggest a selective pressure to maintain this architecture when it has a functional consequence.

In vivo cleavage rules and target repertoire of RNase III in Escherichia coli.

Bacterial RNase III plays important roles in the processing and degradation of RNA transcripts. A major goal is to identify the cleavage targets of this endoribonuclease at a transcriptome-wide scale and delineate its in vivo cleavage rules. Here we applied to Escherichia coli grown to either exponential or stationary phase a tailored RNA-seq-based technology, which allows transcriptome-wide mapping of RNase III cleavage sites at a nucleotide resolution. Our analysis of the large-scale in vivo cleavage data substantiated the established cleavage pattern of a double cleavage in an intra-molecular stem structure, leaving 2-nt-long 3' overhangs, and refined the base-pairing preferences in the cleavage site vicinity. Intriguingly, we observed that the two stem positions between the cleavage sites are highly base-paired, usually involving at least one G-C or C-G base pair. We present a clear distinction between intra-molecular stem structures that are RNase III substrates and intra-molecular stem structures randomly selected across the transcriptome, emphasizing the in vivo specificity of RNase III. Our study provides a comprehensive map of the cleavage sites in both intra-molecular and inter-molecular duplex substrates, providing novel insights into the involvement of RNase III in post-transcriptional regulation in the bacterial cell.

A defense-offense multi-layered regulatory switch in a pathogenic bacterium.

Cells adapt to environmental changes by efficiently adjusting gene expression programs. Staphylococcus aureus, an opportunistic pathogenic bacterium, switches between defensive and offensive modes in response to quorum sensing signal. We identified and studied the structural characteristics and dynamic properties of the core regulatory circuit governing this switch by deterministic and stochastic computational methods, as well as experimentally. This module, termed here Double Selector Switch (DSS), comprises the RNA regulator RNAIII and the transcription factor Rot, defining a double-layered switch involving both transcriptional and post-transcriptional regulations. It coordinates the inverse expression of two sets of target genes, immuno-modulators and exotoxins, expressed during the defensive and offensive modes, respectively. Our computational and experimental analyses show that the DSS guarantees fine-tuned coordination of the inverse expression of its two gene sets, tight regulation, and filtering of noisy signals. We also identified variants of this circuit in other bacterial systems, suggesting it is used as a molecular switch in various cellular contexts and offering its use as a template for an effective switching device in synthetic biology studies.

Evolutionary patterns of Escherichia coli small RNAs and their regulatory interactions.

Most bacterial small RNAs (sRNAs) are post-transcriptional regulators of gene expression, exerting their regulatory function by base-pairing with their target mRNAs. While it has become evident that sRNAs play central regulatory roles in the cell, little is known about their evolution and the evolution of their regulatory interactions. Here we used the prokaryotic phylogenetic tree to reconstruct the evolutionary history of Escherichia coli sRNAs and their binding sites on target mRNAs. We discovered that sRNAs currently present in E. coli mainly accumulated inside the Enterobacteriales order, succeeding the appearance of other types of noncoding RNAs and concurrently with the evolution of a variant of the Hfq protein exhibiting a longer C-terminal region. Our analysis of the evolutionary ages of sRNA-mRNA interactions revealed that while all sRNAs were evolutionarily older than most of their known binding sites on mRNA targets, for quite a few sRNAs there was at least one binding site that coappeared with or preceded them. It is conceivable that the establishment of these first interactions forced selective pressure on the sRNAs, after which additional targets were acquired by fitting a binding site to the active region of the sRNA. This conjecture is supported by the appearance of many binding sites on target mRNAs only after the sRNA gain, despite the prior presence of the target gene in ancestral genomes. Our results suggest a selective mechanism that maintained the sRNAs across the phylogenetic tree, and shed light on the evolution of E. coli post-transcriptional regulatory network.

Auxiliary tRNAs: large-scale analysis of tRNA genes reveals patterns of tRNA repertoire dynamics.

Decoding of all codons can be achieved by a subset of tRNAs. In bacteria, certain tRNA species are mandatory, while others are auxiliary and are variably used. It is currently unknown how this variability has evolved and whether it provides an adaptive advantage. Here we shed light on the subset of auxiliary tRNAs, using genomic data from 319 bacteria. By reconstructing the evolution of tRNAs we show that the auxiliary tRNAs are highly dynamic, being frequently gained and lost along the phylogenetic tree, with a clear dominance of loss events for most auxiliary tRNA species. We reveal distinct co-gain and co-loss patterns for subsets of the auxiliary tRNAs, suggesting that they are subjected to the same selection forces. Controlling for phylogenetic dependencies, we find that the usage of these tRNA species is positively correlated with GC content and may derive directly from nucleotide bias or from preference of Watson-Crick codon-anticodon interactions. Our results highlight the highly dynamic nature of these tRNAs and their complicated balance with codon usage.

Global regulation of transcription by a small RNA: a quantitative view.

Small RNAs are integral regulators of bacterial gene expression, the majority of which act posttranscriptionally by basepairing with target mRNAs, altering translation or mRNA stability. 6S RNA, however, is a small RNA that is a transcriptional regulator, acting by binding directly to σ(70)-RNA polymerase (σ(70)-RNAP) and preventing its binding to gene promoters. At the transition from exponential to stationary phase, 6S RNA accumulates and globally downregulates the transcription of hundreds of genes. At the transition from stationary to exponential phase (outgrowth), 6S RNA is released from σ(70)-RNAP, resulting in a fast increase in free σ(70)-RNAP and transcription of many genes. The transition from stationary to exponential phase is sharp, and is thus accessible for experimental study. However, the transition from exponential to stationary phase is gradual and complicated by changes in other factors, making it more difficult to isolate 6S RNA effects experimentally at this transition. Here, we use mathematical modeling and simulation to study the dynamics of 6S RNA-dependent regulation, focusing on transitions in growth mediated by altered nutrient availability. We first show that our model reproduces the sharp increase in σ(70)-RNAP at outgrowth, as well as the behavior of two experimentally tested mutants, thus justifying its use for characterizing the less accessible dynamics of the transition from exponential to stationary phase. We characterize the dynamics of the two transitions for Escherichia coli wild-type, as well as for mutants with various 6S RNA-RNAP affinities, demonstrating that the 6S RNA regulation mechanism is generally robust to a wide range of such mutations, although the level of regulation at single promoters and their resulting expression fold change will be altered with changes in affinity. Our results provide insight into the potential advantage of transcription regulation by 6S RNA, as it enables storage and efficient release of σ(70)-RNAP during transitions in nutrient availability, which is likely to give a competitive advantage to cells encountering diverse environmental conditions.

Interactions between distant ceRNAs in regulatory networks.

Competing endogenous RNAs (ceRNAs) were recently introduced as RNA transcripts that affect each other's expression level through competition for their microRNA (miRNA) coregulators. This stems from the bidirectional effects between miRNAs and their target RNAs, where a change in the expression level of one target affects the level of the miRNA regulator, which in turn affects the level of other targets. By the same logic, miRNAs that share targets compete over binding to their common targets and therefore also exhibit ceRNA-like behavior. Taken together, perturbation effects could propagate in the posttranscriptional regulatory network through a path of coregulated targets and miRNAs that share targets, suggesting the existence of distant ceRNAs. Here we study the prevalence of distant ceRNAs and their effect in cellular networks. Analyzing the network of miRNA-target interactions deciphered experimentally in HEK293 cells, we show that it is a dense, intertwined network, suggesting that many nodes can act as distant ceRNAs of one another. Indeed, using gene expression data from a perturbation experiment, we demonstrate small, yet statistically significant, changes in gene expression caused by distant ceRNAs in that network. We further characterize the magnitude of the propagated perturbation effect and the parameters affecting it by mathematical modeling and simulations. Our results show that the magnitude of the effect depends on the generation and degradation rates of involved miRNAs and targets, their interaction rates, the distance between the ceRNAs and the topology of the network. Although demonstrated for a miRNA-mRNA regulatory network, our results offer what to our knowledge is a new view on various posttranscriptional cellular networks, expanding the concept of ceRNAs and implying possible distant cross talk within the network, with consequences for the interpretation of indirect effects of gene perturbation.

A genome-wide view of the expression and processing patterns of Thermus thermophilus HB8 CRISPR RNAs.

The CRISPR-Cas system represents an RNA-based adaptive immune response system in prokaryotes and archaea. CRISPRs (clustered regularly interspaced short palindromic repeats) consist of arrays of short repeat-sequences interspaced by nonrepetitive short spacers, some of which show sequence similarity to foreign phage genetic elements. Their cistronic transcripts are processed to produce the mature CRISPR RNAs (crRNAs), the elements that confer immunity by base-pairing with exogenous nucleic acids. We characterized the expression and processing patterns of Thermus thermophilus HB8 CRISPRs by using differential deep-sequencing, which differentiates between 5' monophosphate and 5' non-monophosphate-containing RNAs and/or between 3' hydroxyl and 3' non-hydroxyl-containing RNAs. The genome of T. thermophilus HB8 encodes 11 CRISPRs, classified into three distinct repeat-sequence types, all of which were constitutively expressed without deliberately infecting the bacteria with phage. Analysis of the differential deep sequencing data suggested that crRNAs are generated by endonucleolytic cleavage, leaving fragments with 5' hydroxyl and 3' phosphate or 2',3'-cyclic phosphate termini. The 5' ends of all crRNAs are generated by site-specific cleavage 8 nucleotides upstream of the spacer first position; however, the 3' ends are generated by two alternative, repeat-sequence-type-dependent mechanisms. These observations are consistent with the operation of multiple crRNA processing systems within a bacterial strain.

In vivo mapping of RNA-RNA interactions in Staphylococcus aureus using the endoribonuclease III.

Ribonucleases play key roles in gene regulation and in the expression of virulence factors in Staphylococcus aureus. Among these enzymes, the double-strand specific endoribonuclease III (RNase III) is a key mediator of mRNA processing and degradation. Recently, we have defined, direct target sites for RNase III processing on a genome-wide scale in S. aureus. Our approach is based on deep sequencing of cDNA libraries obtained from RNAs isolated by in vivo co-immunoprecipitation with wild-type RNase III and two cleavage-defective mutants. The use of such catalytically inactivated enzymes, which still retain their RNA binding capacity, allows the identification of novel RNA substrates of RNase III. In this report, we will summarize the diversity of RNase III functions, discuss the advantages and the limitations of the approach, and how this strategy identifies novel mRNA targets of small non-coding RNAs in S. aureus.

Codon usage bias in prokaryotic pyrimidine-ending codons is associated with the degeneracy of the encoded amino acids.

Synonymous codons are unevenly distributed among genes, a phenomenon termed codon usage bias. Understanding the patterns of codon bias and the forces shaping them is a major step towards elucidating the adaptive advantage codon choice can confer at the level of individual genes and organisms. Here, we perform a large-scale analysis to assess codon usage bias pattern of pyrimidine-ending codons in highly expressed genes in prokaryotes. We find a bias pattern linked to the degeneracy of the encoded amino acid. Specifically, we show that codon-pairs that encode two- and three-fold degenerate amino acids are biased towards the C-ending codon while codons encoding four-fold degenerate amino acids are biased towards the U-ending codon. This codon usage pattern is widespread in prokaryotes, and its strength is correlated with translational selection both within and between organisms. We show that this bias is associated with an improved correspondence with the tRNA pool, avoidance of mis-incorporation errors during translation and moderate stability of codon-anticodon interaction, all consistent with more efficient translation.

A functional selection model explains evolutionary robustness despite plasticity in regulatory networks.

Evolutionary rewiring of regulatory networks is an important source of diversity among species. Previous evidence suggested substantial divergence of regulatory networks across species. However, systematically assessing the extent of this plasticity and its functional implications has been challenging due to limited experimental data and the noisy nature of computational predictions. Here, we introduce a novel approach to study cis-regulatory evolution, and use it to trace the regulatory history of 88 DNA motifs of transcription factors across 23 Ascomycota fungi. While motifs are conserved, we find a pervasive gain and loss in the regulation of their target genes. Despite this turnover, the biological processes associated with a motif are generally conserved. We explain these trends using a model with a strong selection to conserve the overall function of a transcription factor, and a much weaker selection over the specific genes it targets. The model also accounts for the turnover of bound targets measured experimentally across species in yeasts and mammals. Thus, selective pressures on regulatory networks mostly tolerate local rewiring, and may allow for subtle fine-tuning of gene regulation during evolution.

A dynamic view of domain-motif interactions.

Many protein-protein interactions are mediated by domain-motif interaction, where a domain in one protein binds a short linear motif in its interacting partner. Such interactions are often involved in key cellular processes, necessitating their tight regulation. A common strategy of the cell to control protein function and interaction is by post-translational modifications of specific residues, especially phosphorylation. Indeed, there are motifs, such as SH2-binding motifs, in which motif phosphorylation is required for the domain-motif interaction. On the contrary, there are other examples where motif phosphorylation prevents the domain-motif interaction. Here we present a large-scale integrative analysis of experimental human data of domain-motif interactions and phosphorylation events, demonstrating an intriguing coupling between the two. We report such coupling for SH3, PDZ, SH2 and WW domains, where residue phosphorylation within or next to the motif is implied to be associated with switching on or off domain binding. For domains that require motif phosphorylation for binding, such as SH2 domains, we found coupled phosphorylation events other than the ones required for domain binding. Furthermore, we show that phosphorylation might function as a double switch, concurrently enabling interaction of the motif with one domain and disabling interaction with another domain. Evolutionary analysis shows that co-evolution of the motif and the proximal residues capable of phosphorylation predominates over other evolutionary scenarios, in which the motif appeared before the potentially phosphorylated residue, or vice versa. Our findings provide strengthening evidence for coupled interaction-regulation units, defined by a domain-binding motif and a phosphorylated residue.

RepTar: a database of predicted cellular targets of host and viral miRNAs.

Computational identification of putative microRNA (miRNA) targets is an important step towards elucidating miRNA functions. Several miRNA target-prediction algorithms have been developed followed by publicly available databases of these predictions. Here we present a new database offering miRNA target predictions of several binding types, identified by our recently developed modular algorithm RepTar. RepTar is based on identification of repetitive elements in 3'-UTRs and is independent of both evolutionary conservation and conventional binding patterns (i.e. Watson-Crick pairing of 'seed' regions). The modularity of RepTar enables the prediction of targets with conventional seed sites as well as rarer targets with non-conventional sites, such as sites with seed wobbles (G-U pairing in the seed region), 3'-compensatory sites and the newly discovered centered sites. Furthermore, RepTar's independence of conservation enables the prediction of cellular targets of the less evolutionarily conserved viral miRNAs. Thus, the RepTar database contains genome-wide predictions of human and mouse miRNAs as well as predictions of cellular targets of human and mouse viral miRNAs. These predictions are presented in a user-friendly database, which allows browsing through the putative sites as well as conducting simple and advanced queries including data intersections of various types. The RepTar database is available at http://reptar.ekmd.huji.ac.il.

Legionella pneumophila 6S RNA optimizes intracellular multiplication.

Legionella pneumophila is a Gram-negative opportunistic human pathogen that infects and multiplies in a broad range of phagocytic protozoan and mammalian phagocytes. Based on the observation that small regulatory RNAs (sRNAs) play an important role in controlling virulence-related genes in several pathogenic bacteria, we attempted to identify sRNAs expressed by L. pneumophila. We used computational prediction followed by experimental verification to identify and characterize sRNAs encoded in the L. pneumophila genome. A 50-mer probe microarray was constructed to test the expression of predicted sRNAs in bacteria grown under a variety of conditions. This strategy successfully identified 22 expressed RNAs, out of which 6 were confirmed by northern blot and RACE. One of the identified sRNAs is highly expressed in postexponential phase, and computational prediction of its secondary structure reveals a striking similarity to the structure of 6S RNA, a widely distributed prokaryotic sRNA, known to regulate the activity of sigma(70)-containing RNA polymerase. A 70-mer probe microarray was used to identify genes affected by L. pneumophila 6S RNA in stationary phase. The 6S RNA positively regulates expression of genes encoding type IVB secretion system effectors, stress response genes such as groES and recA, as well as many genes involved in acquisition of nutrients and genes with unknown or hypothetical functions. Deletion of 6S RNA significantly reduced L. pneumophila intracellular multiplication in both protist and mammalian host cells, but had no detectable effect on growth in rich media.

Convergent Within-Host Adaptation of Pseudomonas aeruginosa through the Transcriptional Regulatory Network.

Bacteria adapt to their host by mutating specific genes and by reprogramming their gene expression. Different strains of a bacterial species often mutate the same genes during infection, demonstrating convergent genetic adaptation. However, there is limited evidence for convergent adaptation at the transcriptional level. To this end, we utilize genomic data of 114 Pseudomonas aeruginosa strains, derived from patients with chronic pulmonary infection, and the P. aeruginosa transcriptional regulatory network. Relying on loss-of-function mutations in genes encoding transcriptional regulators and predicting their effects through the network, we demonstrate predicted expression changes of the same genes in different strains through different paths in the network, implying convergent transcriptional adaptation. Furthermore, through the transcription lens we associate yet-unknown processes, such as ethanol oxidation and glycine betaine catabolism, with P. aeruginosa host adaptation. We also find that known adaptive phenotypes, including antibiotic resistance, which were identified before as achieved by specific mutations, are achieved also through transcriptional changes. Our study has revealed novel interplay between the genetic and transcriptional levels in host adaptation, demonstrating the versatility of the adaptive arsenal of bacterial pathogens and their ability to adapt to the host conditions in a myriad of ways. IMPORTANCE Pseudomonas aeruginosa causes significant morbidity and mortality. The pathogen's remarkable ability to establish chronic infections greatly depends on its adaptation to the host environment. Here, we use the transcriptional regulatory network to predict expression changes during adaptation. We expand the processes and functions known to be involved in host adaptation. We show that the pathogen modulates the activity of genes during adaptation, including genes implicated in antibiotic resistance, both directly via genomic mutations and indirectly via mutations in transcriptional regulators. Furthermore, we detect a subgroup of genes whose predicted changes in expression are associated with mucoid strains, a major adaptive phenotype in chronic infections. We propose that these genes constitute the transcriptional arm of the mucoid adaptive strategy. Identification of different adaptive strategies utilized by pathogens during chronic infection has major promise in the treatment of persistent infections and opens the door to personalized tailored antibiotic treatment in the future.