Regulatory roles of Escherichia coli 5' UTR and ORF-internal RNAs detected by 3' end mapping

The LysR family transcriptional regulator ORF-L16 regulates spinosad biosynthesis in Saccharopolyspora spinosa

Spinosad, a potent broad-spectrum bioinsecticide produced by Saccharopolyspora spinosa, has significant market potential. Despite its effectiveness, the regulatory mechanisms of spinosad biosynthesis remain unclear. Our investigation identified the crucial role of the LysR family transcriptional regulator ORF-L16, located upstream of spinosad biosynthetic genes, in spinosad biosynthesis. Through reverse transcription PCR (RT-PCR) and 5'-rapid amplification of cDNA ends (5'-Race), we unveiled that the spinosad biosynthetic gene cluster (BGC) contains six transcription units and seven promoters. Electrophoretic mobility shift assays (EMSAs) demonstrated that ORF-L16 bound to seven promoters within the spinosad BGC, indicating its involvement in regulating spinosad biosynthesis. Notably, deletion of ORF-L16 led to a drastic reduction in spinosad production from 1818.73 mg/L to 1.69 mg/L, accompanied by decreased transcription levels of spinosad biosynthetic genes, confirming its positive regulatory function. Additionally, isothermal titration calorimetry (ITC) and EMSA confirmed that spinosyn A, the main product of the spinosad BGC, served as an effector of ORF-L16. Specifically, it decreased the binding affinity between ORF-L16 and spinosad BGC promoters, thus exerting negative feedback regulation on spinosad biosynthesis. This research enhances our comprehension of spinosad biosynthesis regulation and lays the groundwork for future investigations on transcriptional regulators in S. spinosa.

A simple BLASTn-based approach generates novel insights into the regulation and biological function of type I toxin-antitoxins

Bacterial chromosomal type I toxin-antitoxin systems consist of a small protein, typically under 60 amino acids, and a small RNA (sRNA) that represses toxin translation. These gene pairs have gained attention over the last decade for their contribution to antibiotic persistence and phage tolerance in bacteria. However, biological functions for many remain elusive as gene deletions often fail to produce an observable phenotype. For many pairs, it is still unknown when the toxin and/or antitoxin gene are natively expressed within the bacterium. We examined sequence conservation of three type I toxin-antitoxin systems, tisB/istR-1, shoB/ohsC, and zor/orz, in over 2,000 Escherichia coli strains, including pathogenic and commensal isolates. Using our custom database, we found that these gene pairs are widespread across E. coli and have expression potential via BLASTn. We identified an alternative, dominant sequence variant of TisB and confirmed that it is toxic upon overproduction. Additionally, analyses revealed a highly conserved sequence in the zorO mRNA untranslated region that is required for full toxicity. We further noted that over 30% of E. coli genomes contain an orz antitoxin gene only and confirmed its expression in a representative strain: the first confirmed report of a type I antitoxin without its cognate toxin. Our results add to our understanding of these systems, and our methodology is applicable for other type I loci to identify critical regulatory and functional features.IMPORTANCEChromosomal type I toxin-antitoxins are a class of genes that have gained increasing attention over the last decade for their roles in antibiotic persistence which may contribute to therapeutic failures. However, the control of many of these genes and when they function have remained elusive. We demonstrate that a simple genetic conservation-based approach utilizing free, publicly available data yields known and novel insights into the regulation and function of three chromosomal type I toxin-antitoxins in Escherichia coli. This study also provides a framework for how this approach could be applied to other genes of interest.

Genetic evidence for functional diversification of gram-negative intermembrane phospholipid transporters

The outer membrane of gram-negative bacteria is a barrier to chemical and physical stress. Phospholipid transport between the inner and outer membranes has been an area of intense investigation and, in E. coli K-12, it has recently been shown to be mediated by YhdP, TamB, and YdbH, which are suggested to provide hydrophobic channels for phospholipid diffusion, with YhdP and TamB playing the major roles. However, YhdP and TamB have different phenotypes suggesting distinct functions. It remains unclear whether these functions are related to phospholipid metabolism. We investigated a synthetic cold sensitivity caused by deletion of fadR, a transcriptional regulator controlling fatty acid degradation and unsaturated fatty acid production, and yhdP, but not by ΔtamB ΔfadR or ΔydbH ΔfadR. Deletion of tamB recuses the ΔyhdP ΔfadR cold sensitivity further demonstrating the phenotype is related to functional diversification between these genes. The ΔyhdP ΔfadR strain shows a greater increase in cardiolipin upon transfer to the non-permissive temperature and genetically lowering cardiolipin levels can suppress cold sensitivity. These data also reveal a qualitative difference between cardiolipin synthases in E. coli, as deletion of clsA and clsC suppresses cold sensitivity but deletion of clsB does not. Moreover, increased fatty acid saturation is necessary for cold sensitivity and lowering this level genetically or through supplementation of oleic acid suppresses the cold sensitivity of the ΔyhdP ΔfadR strain. Together, our data clearly demonstrate that the diversification of function between YhdP and TamB is related to phospholipid metabolism. Although indirect regulatory effects are possible, we favor the parsimonious hypothesis that YhdP and TamB have differential phospholipid-substrate transport preferences. Thus, our data provide a potential mechanism for independent control of the phospholipid composition of the inner and outer membranes in response to changing conditions based on regulation of abundance or activity of YhdP and TamB.

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.

Genetic evidence for functional diversification of gram-negative intermembrane phospholipid transporters

The outer membrane of Gram-negative bacteria is a barrier to chemical and physical stress. Phospholipid transport between the inner and outer membranes has been an area of intense investigation and, in E. coli K-12, it has recently been shown to be mediated by YhdP, TamB, and YdbH, which are suggested to provide hydrophobic channels for phospholipid diffusion, with YhdP and TamB playing the major roles. However, YhdP and TamB have different phenotypes suggesting distinct functions. We investigated these functions using synthetic cold sensitivity (at 30 °C) caused by deletion of yhdP and fadR, a transcriptional regulator controlling fatty acid degradation and unsaturated fatty acid production, but not by ΔtamB ΔfadR or ΔydbH ΔfadR,. Deletion of tamB suppresses the ΔyhdP ΔfadR cold sensitivity suggesting this phenotype is related to phospholipid transport. The ΔyhdP ΔfadR strain shows a greater increase in cardiolipin upon transfer to the non-permissive temperature and genetically lowering cardiolipin levels can suppress cold sensitivity. These data also reveal a qualitative difference between cardiolipin synthases in E. coli, as deletion of clsA and clsC suppresses cold sensitivity but deletion of clsB does not despite lower cardiolipin levels. In addition to increased cardiolipin, increased fatty acid saturation is necessary for cold sensitivity and lowering this level genetically or through supplementation of oleic acid suppresses the cold sensitivity of the ΔyhdP ΔfadR strain. Although indirect effects are possible, we favor the parsimonious hypothesis that YhdP and TamB have differential substrate transport preferences, most likely with YhdP preferentially transporting more saturated phospholipids and TamB preferentially transporting more unsaturated phospholipids. We envision cardiolipin contributing to this transport preference by sterically clogging TamB-mediated transport of saturated phospholipids. Thus, our data provide a potential mechanism for independent control of the phospholipid composition of the inner and outer membranes in response to changing conditions.

Re-defining how mRNA degradation is coordinated with transcription and translation in bacteria

In eukaryotic cells, transcription, translation, and mRNA degradation occur in distinct subcellular regions. How these mRNA processes are organized in bacteria, without employing membrane-bound compartments, remains unclear. Here, we present generalizable principles underlying coordination between these processes in bacteria. In Escherichia coli, we found that co-transcriptional degradation is rare for mRNAs except for those encoding inner membrane proteins, due to membrane localization of the main ribonuclease, RNase E. We further found, by varying ribosome binding sequences, that translation affects mRNA stability not because ribosomes protect mRNA from degradation, but because low translation leads to premature transcription termination in the absence of transcription-translation coupling. Extending our analyses to Bacillus subtilis and Caulobacter crescentus, we established subcellular localization of RNase E (or its homolog) and premature transcription termination in the absence of transcription-translation coupling as key determinants that explain differences in transcriptional and translational coupling to mRNA degradation across genes and species.

Genome-wide annotation of transcript boundaries using bacterial Rend-seq datasets

Accurate annotation to single-nucleotide resolution of the transcribed regions in genomes is key to optimally analyse RNA-seq data, understand regulatory events and for the design of experiments. However, currently most genome annotations provided by GenBank generally lack information about untranslated regions. Additionally, information regarding genomic locations of non-coding RNAs, such as sRNAs, or anti-sense RNAs is frequently missing. To provide such information, diverse RNA-seq technologies, such as Rend-seq, have been developed and applied to many bacterial species. However, incorporating this vast amount of information into annotation files has been limited and is bioinformatically challenging, resulting in UTRs and other non-coding elements being overlooked or misrepresented. To overcome this problem, we present pyRAP (python Rend-seq Annotation Pipeline), a software package that analyses Rend-seq datasets to accurately resolve transcript boundaries genome-wide. We report the use of pyRAP to find novel transcripts, transcript isoforms, and RNase-dependent sRNA processing events. In Bacillus subtilis we uncovered 63 novel transcripts and provide genomic coordinates with single-nucleotide resolution for 2218 5'UTRs, 1864 3'UTRs and 161 non-coding RNAs. In Escherichia coli, we report 117 novel transcripts, 2429 5'UTRs, 1619 3'UTRs and 91 non-coding RNAs, and in Staphylococcus aureus, 16 novel transcripts, 664 5'UTRs, 696 3'UTRs, and 81 non-coding RNAs. Finally, we use pyRAP to produce updated annotation files for B. subtilis 168, E. coli K-12 MG1655, and S. aureus 8325 for use in the wider microbial genomics research community.

Direct and indirect control of Rho-dependent transcription termination by the Escherichia coli lysC riboswitch

Bacterial riboswitches are molecular structures that play a crucial role in controlling gene expression to maintain cellular balance. The Escherichia coli lysC riboswitch has been previously shown to regulate gene expression through translation initiation and mRNA decay. Recent research suggests that lysC gene expression is also influenced by Rho-dependent transcription termination. Through a series of in silico, in vitro, and in vivo experiments, we provide experimental evidence that the lysC riboswitch directly and indirectly modulates Rho transcription termination. Our study demonstrates that Rho-dependent transcription termination plays a significant role in the cotranscriptional regulation of lysC expression. Together with previous studies, our work suggests that lysC expression is governed by a lysine-sensing riboswitch that regulates translation initiation, transcription termination, and mRNA degradation. Notably, both Rho and RNase E target the same region of the RNA molecule, implying that RNase E may degrade Rho-terminated transcripts, providing a means to selectively eliminate these incomplete messenger RNAs. Overall, this study sheds light on the complex regulatory mechanisms used by bacterial riboswitches, emphasizing the role of transcription termination in the control of gene expression and mRNA stability.

Refining the transcriptional landscapes for distinct clades of virulent phages infecting Pseudomonas aeruginosa

The introduction of high-throughput sequencing has resulted in a surge of available bacteriophage genomes, unveiling their tremendous genomic diversity. However, our current understanding of the complex transcriptional mechanisms that dictate their gene expression during infection is limited to a handful of model phages. Here, we applied ONT-cappable-seq to reveal the transcriptional architecture of six different clades of virulent phages infecting Pseudomonas aeruginosa. This long-read microbial transcriptomics approach is tailored to globally map transcription start and termination sites, transcription units, and putative RNA-based regulators on dense phage genomes. Specifically, the full-length transcriptomes of LUZ19, LUZ24, 14-1, YuA, PAK_P3, and giant phage phiKZ during early, middle, and late infection were collectively charted. Beyond pinpointing traditional promoter and terminator elements and transcription units, these transcriptional profiles provide insights in transcriptional attenuation and splicing events and allow straightforward validation of Group I intron activity. In addition, ONT-cappable-seq data can guide genome-wide discovery of novel regulatory element candidates, including noncoding RNAs and riboswitches. This work substantially expands the number of annotated phage-encoded transcriptional elements identified to date, shedding light on the intricate and diverse gene expression regulation mechanisms in Pseudomonas phages, which can ultimately be sourced as tools for biotechnological applications in phage and bacterial engineering.

Exploring the transcriptional landscape of phage-host interactions using novel high-throughput approaches

In the last decade, powerful high-throughput sequencing approaches have emerged to analyse microbial transcriptomes at a global scale. However, to date, applications of these approaches to microbial viruses such as phages remain scarce. Tailoring these techniques to virus-infected bacteria promises to obtain a detailed picture of the underexplored RNA biology and molecular processes during infection. In addition, transcriptome study of stress and perturbations induced by phages in their infected bacterial hosts is likely to reveal new fundamental mechanisms of bacterial metabolism and gene regulation. Here, we provide references and blueprints to implement emerging transcriptomic approaches towards addressing transcriptome architecture, RNA-RNA and RNA-protein interactions, RNA modifications, structures and heterogeneity of transcription profiles in infected cells that will provide guides for future directions in phage-centric therapeutic applications and microbial synthetic biology.

Directed Screening for sRNA Targets in E. coli Using a Plasmid Library

A large number of bacterial small regulatory RNAs (sRNAs) modulate gene expression by base pairing to a target mRNA, affecting its translation or stability. This posttranscriptional regulation has been shown to be essential and critical for bacterial physiology. One of the challenges of studying sRNA signaling is identifying the sRNA regulators of specific genes. Here, we describe a protocol for making an sRNA expression library and using this library to screen for sRNA regulators of genes of interest in E. coli. This library can be easily expanded and adapted to use in other bacteria.

Ribosome profiling reveals the fine-tuned response of Escherichia coli to mild and severe acid stress

IMPORTANCE

Bacteria react very differently to survive in acidic environments, such as the human gastrointestinal tract. Escherichia coli is one of the extremely acid-resistant bacteria and has a variety of acid-defense mechanisms. Here, we provide the first genome-wide overview of the adaptations of E. coli K-12 to mild and severe acid stress at both the transcriptional and translational levels. Using ribosome profiling and RNA sequencing, we uncover novel adaptations to different degrees of acidity, including previously hidden stress-induced small proteins and novel key transcription factors for acid defense, and report mRNAs with pH-dependent differential translation efficiency. In addition, we distinguish between acid-specific adaptations and general stress response mechanisms using denoising autoencoders. This workflow represents a powerful approach that takes advantage of next-generation sequencing techniques and machine learning to systematically analyze bacterial stress responses.

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.

The environmentally-regulated interplay between local three-dimensional chromatin organisation and transcription of proVWX in E. coli

Nucleoid associated proteins (NAPs) maintain the architecture of bacterial chromosomes and regulate gene expression. Thus, their role as transcription factors may involve three-dimensional chromosome re-organisation. While this model is supported by in vitro studies, direct in vivo evidence is lacking. Here, we use RT-qPCR and 3C-qPCR to study the transcriptional and architectural profiles of the H-NS (histone-like nucleoid structuring protein)-regulated, osmoresponsive proVWX operon of Escherichia coli at different osmolarities and provide in vivo evidence for transcription regulation by NAP-mediated chromosome re-modelling in bacteria. By consolidating our in vivo investigations with earlier in vitro and in silico studies that provide mechanistic details of how H-NS re-models DNA in response to osmolarity, we report that activation of proVWX in response to a hyperosmotic shock involves the destabilization of H-NS-mediated bridges anchored between the proVWX downstream and upstream regulatory elements (DRE and URE), and between the DRE and ygaY that lies immediately downstream of proVWX. The re-establishment of these bridges upon adaptation to hyperosmolarity represses the operon. Our results also reveal additional structural features associated with changes in proVWX transcript levels such as the decompaction of local chromatin upstream of the operon, highlighting that further complexity underlies the regulation of this model operon. H-NS and H-NS-like proteins are wide-spread amongst bacteria, suggesting that chromosome re-modelling may be a typical feature of transcriptional control in bacteria.

Uncovering the temporal dynamics and regulatory networks of thermal stress response in a hyperthermophile using transcriptomics and proteomics

Facing rapid fluctuations in their natural environment, extremophiles, like the hyperthermophilic archaeon Pyrococcus furiosus, exhibit remarkable adaptability to extreme conditions. However, our understanding of their dynamic cellular responses remains limited. This study integrates RNA-sequencing and mass spectrometry data, thereby elucidating transcriptomic and proteomic responses to heat and cold shock stress in P. furiosus. Our results reveal rapid and dynamic changes in gene and protein expression following these stress responses. Heat shock triggers extensive transcriptome reprogramming, orchestrated by the transcriptional regulator Phr, targeting a broader gene repertoire than previously demonstrated. For heat shock signature genes, RNA levels swiftly return to baseline upon recovery, while protein levels remain persistently upregulated, reflecting a rapid but sustained response. Intriguingly, cold shock at 4°C elicits distinct short- and long-term responses at both RNA and protein levels. Cluster analysis identified gene sets with either congruent or contrasting trends in RNA and protein changes, representing well-separated arCOG groups tailored to their individual cellular responses. Particularly, upregulation of ribosomal proteins and significant enrichment of 5'-leadered sequences in cold-shock responsive genes suggest that translation regulation is important during cold shock adaption. Further investigating transcriptomic features, we reveal that thermal stress genes are equipped with basal sequence elements, such as strong promoter and poly(U)-terminators, facilitating a regulated response of the respective transcription units. Our study provides a comprehensive overview of the cellular response to temperature stress, advancing our understanding of stress response mechanisms in hyperthermophilic archaea and providing valuable insights into the molecular adaptations that facilitate life in extreme environments.IMPORTANCEExtreme environments provide unique challenges for life, and the study of extremophiles can shed light on the mechanisms of adaptation to such conditions. Pyrococcus furiosus, a hyperthermophilic archaeon, is a model organism for studying thermal stress response mechanisms. In this study, we used an integrated analysis of RNA-sequencing and mass spectrometry data to investigate the transcriptomic and proteomic responses of P. furiosus to heat and cold shock stress and recovery. Our results reveal the rapid and dynamic changes in gene and protein expression patterns associated with these stress responses, as well as the coordinated regulation of different gene sets in response to different stressors. These findings provide valuable insights into the molecular adaptations that facilitate life in extreme environments and advance our understanding of stress response mechanisms in hyperthermophilic archaea.

Methodologies for bacterial ribonuclease characterization using RNA-seq

Bacteria adjust gene expression at the post-transcriptional level through an intricate network of small regulatory RNAs and RNA-binding proteins, including ribonucleases (RNases). RNases play an essential role in RNA metabolism, regulating RNA stability, decay, and activation. These enzymes exhibit species-specific effects on gene expression, bacterial physiology, and different strategies of target recognition. Recent advances in high-throughput RNA sequencing (RNA-seq) approaches have provided a better understanding of the roles and modes of action of bacterial RNases. Global studies aiming to identify direct targets of RNases have highlighted the diversity of RNase activity and RNA-based mechanisms of gene expression regulation. Here, we review recent RNA-seq approaches used to study bacterial RNases, with a focus on the methods for identifying direct RNase targets.

Internal transcription termination widely regulates differential expression of operon-organized genes including ribosomal protein and RNA polymerase genes in an archaeon

Genes organized within operons in prokaryotes benefit from coordinated expression. However, within many operons, genes are expressed at different levels, and the mechanisms for this remain obscure. By integrating PacBio-seq, dRNA-seq, Term-seq and Illumina-seq data of a representative archaeon Methanococcus maripaludis, internal transcription termination sites (ioTTSs) were identified within 38% of operons. Higher transcript and protein abundances were found for genes upstream than downstream of ioTTSs. For representative operons, these differences were confirmed by northern blotting, qRT-PCR and western blotting, demonstrating that these ioTTS terminations were functional. Of special interest, mutation of ioTTSs in ribosomal protein (RP)-RNA polymerase (RNAP) operons not only elevated expression of the downstream RNAP genes but also decreased production of the assembled RNAP complex, slowed whole cell transcription and translation, and inhibited growth. Overexpression of the RNAP subunits with a shuttle vector generated the similar physiological effects. Therefore, ioTTS termination is a general and physiologically significant regulatory mechanism of the operon gene expression. Because the RP-RNAP operons are found to be widely distributed in archaeal species, this regulatory mechanism could be commonly employed in archaea.

Improved RNA stability estimation indicates that transcriptional interference is frequent in diverse bacteria

We used stochastic simulations and experimental data from E. coli, K. aerogenes, Synechococcus PCC 7002 and Synechocystis PCC 6803 to provide evidence that transcriptional interference via the collision mechanism is likely a prevalent mechanism for bacterial gene regulation. Rifampicin time-series data can be used to globally monitor and quantify collision between sense and antisense transcription-complexes. Our findings also highlight that transcriptional events, such as differential RNA decay, partial termination, and internal transcriptional start sites often deviate from gene annotations. Consequently, within a single gene annotation, there exist transcript segments with varying half-lives and transcriptional properties. To address these complexities, we introduce 'rifi', an R-package that analyzes transcriptomic data from rifampicin time series. 'rifi' employs a dynamic programming-based segmentation approach to identify individual transcripts, enabling accurate assessment of RNA stability and detection of diverse transcriptional events.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen Borrelia burgdorferi

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. Here, we use several RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. We identify complex gene arrangements and operons, untranslated regions and small RNAs. We predict intrinsic terminators and experimentally test examples of Rho-dependent transcription termination. Remarkably, 63% of RNA 3' ends map upstream of or internal to open reading frames (ORFs), including genes involved in the unique infectious cycle of B. burgdorferi. We suggest these RNAs result from premature termination, processing and regulatory events such as cis-acting regulation. Furthermore, the polyamine spermidine globally influences the generation of truncated mRNAs. Collectively, our findings provide insights into transcription termination and uncover an abundance of potential RNA regulators in B. burgdorferi.

Premature termination of transcription is shaped by Rho and translated uORFS in Mycobacterium tuberculosis

Little is known about the decisions behind transcription elongation versus termination in the human pathogen Mycobacterium tuberculosis (M.TB). By applying Term-seq to M.TB we found that the majority of transcription termination is premature and associated with translated regions, i.e., within previously annotated or newly identified open reading frames. Computational predictions and Term-seq analysis, upon depletion of termination factor Rho, suggests that Rho-dependent transcription termination dominates all transcription termination sites (TTS), including those associated with regulatory 5' leaders. Moreover, our results suggest that tightly coupled translation, in the form of overlapping stop and start codons, may suppress Rho-dependent termination. This study provides detailed insights into novel M.TB cis-regulatory elements, where Rho-dependent, conditional termination of transcription and translational coupling together play major roles in gene expression control. Our findings contribute to a deeper understanding of the fundamental regulatory mechanisms that enable M.TB adaptation to the host environment offering novel potential points of intervention.

Extensive diversity in RNA termination and regulation revealed by transcriptome mapping for the Lyme pathogen B. burgdorferi

Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. We employed complementary RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. By systematically mapping B. burgdorferi RNA ends at single nucleotide resolution, we delineated complex gene arrangements and operons and mapped untranslated regions (UTRs) and small RNAs (sRNAs). We experimentally tested modes of B. burgdorferi transcription termination and compared our findings to observations in E. coli , P. aeruginosa , and B. subtilis . We discovered 63% of B. burgdorferi RNA 3' ends map upstream or internal to open reading frames (ORFs), suggesting novel mechanisms of regulation. Northern analysis confirmed the presence of stable 5' derived RNAs from mRNAs encoding gene products involved in the unique infectious cycle of B. burgdorferi . We suggest these RNAs resulted from premature termination and regulatory events, including forms of cis- acting regulation. For example, we documented that the polyamine spermidine globally influences the generation of truncated mRNAs. In one case, we showed that high spermidine concentrations increased levels of RNA fragments derived from an mRNA encoding a spermidine import system, with a concomitant decrease in levels of the full- length mRNA. Collectively, our findings revealed new insight into transcription termination and uncovered an abundance of potential RNA regulators.

Global transcriptional responses of pneumococcus to human blood components and cerebrospinal fluid

Streptococcus pneumoniae (pneumococcus) is a leading cause of severe invasive infectious diseases such as sepsis and meningitis. Understanding how pneumococcus adapts and survive in the human bloodstream environment and cerebrospinal fluid (CSF) is important for development of future treatment strategies. This study investigates the global transcriptional response of pneumococcus to human blood components and CSF acquired from discarded and anonymized patient samples. Extensive transcriptional changes to human blood components were observed during early stages of interaction. Plasma-specific responses were primarily related to metabolic components and include strong downregulation of fatty acid biosynthesis genes, and upregulation of nucleotide biosynthesis genes. No transcriptional responses specific to the active plasma proteins (e.g., complement proteins) were observed during early stages of interaction as demonstrated by a differential expression analysis between plasma and heat-inactivated plasma. The red blood cell (RBC)-specific response was far more complex, and included activation of the competence system, differential expression of several two-component systems, phosphotransferase systems and transition metal transporter genes. Interestingly, most of the changes observed for CSF were also observed for plasma. One of the few CSF-specific responses, not observed for plasma, was a strong downregulation of the iron acquisition system piuBCDA. Intriguingly, this transcriptomic analysis also uncovers significant differential expression of more than 20 small non-coding RNAs, most of them in response to RBCs, including small RNAs from uncharacterized type I toxin-antitoxin systems. In summary, this transcriptomic study identifies key pneumococcal metabolic pathways and regulatory genes involved with adaptation to human blood and CSF. Future studies should uncover the potential involvement of these factors with virulence in-vivo.

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

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

High-Resolution Small RNAs Landscape Provides Insights into Alkane Adaptation in the Marine Alkane-Degrader Alcanivorax dieselolei B-5

Alkanes are widespread in the ocean, and Alcanivorax is one of the most ubiquitous alkane-degrading bacteria in the marine ecosystem. Small RNAs (sRNAs) are usually at the heart of regulatory pathways, but sRNA-mediated alkane metabolic adaptability still remains largely unknown due to the difficulties of identification. Here, differential RNA sequencing (dRNA-seq) modified with a size selection (~50-nt to 500-nt) strategy was used to generate high-resolution sRNAs profiling in the model species Alcanivorax dieselolei B-5 under alkane (n-hexadecane) and non-alkane (acetate) conditions. As a result, we identified 549 sRNA candidates at single-nucleotide resolution of 5'-ends, 63.4% of which are with transcription start sites (TSSs), and 36.6% of which are with processing sites (PSSs) at the 5'-ends. These sRNAs originate from almost any location in the genome, regardless of intragenic (65.8%), antisense (20.6%) and intergenic (6.2%) regions, and RNase E may function in the maturation of sRNAs. Most sRNAs locally distribute across the 15 reference genomes of Alcanivorax, and only 7.5% of sRNAs are broadly conserved in this genus. Expression responses to the alkane of several core conserved sRNAs, including 6S RNA, M1 RNA and tmRNA, indicate that they may participate in alkane metabolisms and result in more actively global transcription, RNA processing and stresses mitigation. Two novel CsrA-related sRNAs are identified, which may be involved in the translational activation of alkane metabolism-related genes by sequestering the global repressor CsrA. The relationships of sRNAs with the characterized genes of alkane sensing (ompS), chemotaxis (mcp, cheR, cheW2), transporting (ompT1, ompT2, ompT3) and hydroxylation (alkB1, alkB2, almA) were created based on the genome-wide predicted sRNA-mRNA interactions. Overall, the sRNA landscape lays the ground for uncovering cryptic regulations in critical marine bacterium, among which both the core and species-specific sRNAs are implicated in the alkane adaptive metabolisms.

Transcription-Translation Coupling in Bacteria

In bacteria, transcription and translation take place in the same cellular compartment. Therefore, a messenger RNA can be translated as it is being transcribed, a process known as transcription-translation coupling. This process was already recognized at the dawn of molecular biology, yet the interplay between the two key players, the RNA polymerase and ribosome, remains elusive. Genetic data indicate that an RNA sequence can be translated shortly after it has been transcribed. The closer both processes are in time, the less accessible the RNA sequence is between the RNA polymerase and ribosome. This temporal coupling has important consequences for gene regulation. Biochemical and structural studies have detailed several complexes between the RNA polymerase and ribosome. The in vivo relevance of this physical coupling has not been formally demonstrated. We discuss how both temporal and physical coupling may mesh to produce the phenomenon we know as transcription-translation coupling.

Bacterial Transcription Factors Bind to Coding Regions and Regulate Internal Cryptic Promoters

Transcription factors (TFs) regulate transcription by binding to the specific sequences at the promoter region. However, the mechanisms and functions of TFs binding within the coding sequences (CDS) remain largely elusive in prokaryotes. To this end, we collected 409 data sets for bacterial TFs, including 104 chromatin immunoprecipitation sequencing (ChIP-seq) assays and 305 data sets from the systematic evolution of ligands by exponential enrichment (SELEX) in seven model bacteria. Interestingly, these TFs displayed the same binding capabilities for both coding and intergenic regions. Subsequent biochemical and genetic experiments demonstrated that several TFs bound to the coding regions and regulated the transcription of the binding or adjacent genes. Strand-specific RNA sequencing revealed that these CDS-binding TFs regulated the activity of the cryptic promoters, resulting in the altered transcription of the corresponding antisense RNA. TF RhpR hindered the transcriptional elongation of a subgenic transcript within a CDS. A ChIP-seq and Ribo-seq coanalysis revealed that RhpR influenced the translational efficiency of binding genes. Taken together, the present study reveals three regulatory mechanisms of CDS-bound TFs within individual genes, operons, and antisense RNAs, which demonstrate the variability of the regulatory mechanisms of TFs and expand upon the complexity of bacterial transcriptomes.

A scalable framework for the discovery of functional helicase substrates and helicase-driven regulatory switches

Helicases are ubiquitous motor enzymes that remodel nucleic acids (NA) and NA-protein complexes in key cellular processes. To explore the functional repertoire and specificity landscape of helicases, we devised a screening scheme-Helicase-SELEX (Systematic Evolution of Ligands by EXponential enrichment)-that enzymatically probes substrate and cofactor requirements at global scale. Using the transcription termination Rho helicase of Escherichia coli as a prototype for Helicase-SELEX, we generated a genome-wide map of Rho utilization (Rut) sites. The map reveals many features, including promoter- and intrinsic terminator-associated Rut sites, bidirectional Rut tandems, and cofactor-dependent Rut sites with inverted G > C skewed compositions. We also implemented an H-SELEX variant where we used a model ligand, serotonin, to evolve synthetic Rut sites operating in vitro and in vivo in a ligand-dependent manner. Altogether, our data illustrate the power and flexibility of Helicase-SELEX to seek constitutive or conditional helicase substrates in natural or synthetic NA libraries for fundamental or synthetic biology discovery.

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.

RNase III-CLASH of multi-drug resistant Staphylococcus aureus reveals a regulatory mRNA 3'UTR required for intermediate vancomycin resistance

Treatment of methicillin-resistant Staphylococcus aureus infections is dependent on the efficacy of last-line antibiotics including vancomycin. Treatment failure is commonly linked to isolates with intermediate vancomycin resistance (termed VISA). These isolates have accumulated point mutations that collectively reduce vancomycin sensitivity, often by thickening the cell wall. Changes in regulatory small RNA expression have been correlated with antibiotic stress in VISA isolates however the functions of most RNA regulators is unknown. Here we capture RNA-RNA interactions associated with RNase III using CLASH. RNase III-CLASH uncovers hundreds of novel RNA-RNA interactions in vivo allowing functional characterisation of many sRNAs for the first time. Surprisingly, many mRNA-mRNA interactions are recovered and we find that an mRNA encoding a long 3' untranslated region (UTR) (termed vigR 3'UTR) functions as a regulatory 'hub' within the RNA-RNA interaction network. We demonstrate that the vigR 3'UTR promotes expression of folD and the cell wall lytic transglycosylase isaA through direct mRNA-mRNA base-pairing. Deletion of the vigR 3'UTR re-sensitised VISA to glycopeptide treatment and both isaA and vigR 3'UTR deletions impact cell wall thickness. Our results demonstrate the utility of RNase III-CLASH and indicate that S. aureus uses mRNA-mRNA interactions to co-ordinate gene expression more widely than previously appreciated.

The Small RNA CyaR Activates Translation of the Outer Membrane Haem Receptor chuA in Enterohemorrhagic Escherichia coli

To sense the transition from environment to host, bacteria use a range of environmental cues to control expression of virulence genes. Iron is tightly sequestered in host tissues and in the human pathogen enterohemorrhagic Escherichia coli (EHEC) iron-limitation induces transcription of the outer membrane haem transporter encoded by chuAS. ChuA expression is post-transcriptionally activated at 37°C by a FourU RNA thermometer ensuring that the haem receptor is only expressed under low iron, high temperature conditions that indicate the host. Here we demonstrate that expression of chuA is also independently regulated by the cAMP-responsive small RNA (sRNA) CyaR and transcriptional terminator Rho. These results indicate that chuAS expression is regulated at the transcription initiation, transcript elongation, and translational level. We speculate that additional sensing of the gluconeogenic environment allows further precision in determining when EHEC is at the gastrointestinal epithelium of the host. With previous studies, it appears that the chuAS transcript is controlled by eight regulatory inputs that control expression through six different transcriptional and post-transcriptional mechanisms. The results highlight the ability of regulatory sRNAs to integrate multiple environmental signals into a layered hierarchy of signal input.

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

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

Reassembling a cannon in the DNA defense arsenal: Genetics of StySA, a BREX phage exclusion system in Salmonella lab strains

Understanding mechanisms that shape horizontal exchange in prokaryotes is a key problem in biology. A major limit on DNA entry is imposed by restriction-modification (RM) processes that depend on the pattern of DNA modification at host-specified sites. In classical RM, endonucleolytic DNA cleavage follows detection of unprotected sites on entering DNA. Recent investigation has uncovered BREX (BacteRiophage EXclusion) systems. These RM-like activities employ host protection by DNA modification, but immediate replication arrest occurs without evident of nuclease action on unmodified phage DNA. Here we show that the historical stySA RM locus of Salmonella enterica sv Typhimurium is a variant BREX system. A laboratory strain disabled for both the restriction and methylation activity of StySA nevertheless has wild type sequence in pglX, the modification gene homolog. Instead, flanking genes pglZ and brxC each carry multiple mutations (μ) in their C-terminal domains. We further investigate this system in situ, replacing the mutated pglZμ and brxCμ genes with the WT counterpart. PglZ-WT supports methylation in the presence of either BrxCμ or BrxC-WT but not in the presence of a deletion/insertion allele, ΔbrxC::cat. Restriction requires both BrxC-WT and PglZ-WT, implicating the BrxC C-terminus specifically in restriction activity. These results suggests that while BrxC, PglZ and PglX are principal components of the BREX modification activity, BrxL is required for restriction only. Furthermore, we show that a partial disruption of brxL disrupts transcription globally.

Horizontal gene transfer is a major driver of evolution and adaptation in bacteria. Genes from outside may be beneficial or dangerous to the receiving cell. Benefits include new food sources such as sugars, or new homes by adhesion, or new resistances, as to antibiotics. Dangers are posed by bacteriophages--viruses that take over the cell machinery, multiply, and release progeny to kill sister cells. Host-dependent restriction-modification systems enable defense that distinguishes relatives from strangers: using a modification pattern (M) carried by DNA bases added by the host cell to prevent restriction (R). Sisters and cousin cells will have the same protective pattern on DNA, while DNA of foreign origin will have the wrong M pattern and be restricted (R, rejected). Typically, restriction involves nuclease digestion. Here we address the enigmatic StySA RM system, one of the earliest to be genetically characterized. It is a variant of the newly recognized defense mechanism, BREX. BREX systems also track DNA history via modification pattern, but restrict by a novel, uncharacterized mechanism. Like other BREX family systems, StySA-BREX modification requires multiple components. When StySA-BREX transcription is unbalanced, we find global disruption of gene transcription. The disruption pattern does not suggest SOS-inducing damage to DNA.

Synthetic 3'-UTR valves for optimal metabolic flux control in Escherichia coli

As the design of genetic circuitry for synthetic biology becomes more sophisticated, diverse regulatory bioparts are required. Despite their importance, well-characterized 3′-untranslated region (3′-UTR) bioparts are limited. Thus, transcript 3′-ends require further investigation to understand the underlying regulatory role and applications of the 3′-UTR. Here, we revisited the use of Term-Seq in the Escherichia coli strain K-12 MG1655 to enhance our understanding of 3′-UTR regulatory functions and to provide a diverse collection of tunable 3′-UTR bioparts with a wide termination strength range. Comprehensive analysis of 1,629 transcript 3′-end positions revealed multiple 3′-termini classes generated through transcription termination and RNA processing. The examination of individual Rho-independent terminators revealed a reduction in downstream gene expression over a wide range, which led to the design of novel synthetic metabolic valves that control metabolic fluxes in branched pathways. These synthetic metabolic valves determine the optimal balance of heterologous pathways for maximum target biochemical productivity. The regulatory strategy using 3′-UTR bioparts is advantageous over promoter- or 5′-UTR-based transcriptional control as it modulates gene expression at transcription levels without trans-acting element requirements (e.g. transcription factors). Our results provide a foundational platform for 3′-UTR engineering in synthetic biology applications.

Nanopore sequencing of RNA and cDNA molecules in Escherichia coli

High-throughput sequencing dramatically changed our view of transcriptome architectures and allowed for ground-breaking discoveries in RNA biology. Recently, sequencing of full-length transcripts based on the single-molecule sequencing platform from Oxford Nanopore Technologies (ONT) was introduced and is widely used to sequence eukaryotic and viral RNAs. However, experimental approaches implementing this technique for prokaryotic transcriptomes remain scarce. Here, we present an experimental and bioinformatic workflow for ONT RNA-seq in the bacterial model organism Escherichia coli, which can be applied to any microorganism. Our study highlights critical steps of library preparation and computational analysis and compares the results to gold standards in the field. Furthermore, we comprehensively evaluate the applicability and advantages of different ONT-based RNA sequencing protocols, including direct RNA, direct cDNA, and PCR-cDNA. We find that (PCR)-cDNA-seq offers improved yield and accuracy compared to direct RNA sequencing. Notably, (PCR)-cDNA-seq is suitable for quantitative measurements and can be readily used for simultaneous and accurate detection of transcript 5' and 3' boundaries, analysis of transcriptional units, and transcriptional heterogeneity. In summary, based on our comprehensive study, we show nanopore RNA-seq to be a ready-to-use tool allowing rapid, cost-effective, and accurate annotation of multiple transcriptomic features. Thereby nanopore RNA-seq holds the potential to become a valuable alternative method for RNA analysis in prokaryotes.

Expression of the DeaD RNA helicase is regulated at multiple levels through its long mRNA 5' untranslated region

DEAD-box proteins (DBPs) are a prominent class of RNA remodeling proteins that alter RNA structure, a process they typically perform through an ATP-dependent RNA helicase activity. Although many DBPs have been characterized at the structural and functional level in detail, much less is known about how they are regulated. We previously showed that the messenger RNA (mRNA) for the Escherichia coli (E. coli) DeaD DBP contains an unusually long 5' untranslated region (5' UTR) of 838 nucleotides (nts) and that it is the primary RNA determinant of DeaD autoregulation. We speculated that such a long and complex 5' UTR might regulate deaD expression in additional ways. Here we show that the deaD mRNA 5' UTR regulates deaD expression at two additional levels: temperature dependent expression and through a stem-loop structure overlapping the start codon. These results support the hypothesis that a long 5' UTR can regulate gene expression through multiple mechanisms. Importance The expression of genes is frequently regulated by determinants with the 5' UTR. Although many different regulatory mechanisms that operate via the 5' UTR have been described, the functional relevance of genes with long UTRs is less clear. Here, we show that the 838 nt long 5' UTR in the deaD mRNA regulates the expression of DeaD at multiple levels. We propose that long UTRs originate to provide precise control of gene expression through multiple regulatory mechanisms, and they are indicators of the importance of their associated gene products for cellular adaptation to different environments.

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.

Overlapping genes in natural and engineered genomes

Modern genome-scale methods that identify new genes, such as proteogenomics and ribosome profiling, have revealed, to the surprise of many, that overlap in genes, open reading frames and even coding sequences is widespread and functionally integrated into prokaryotic, eukaryotic and viral genomes. In parallel, the constraints that overlapping regions place on genome sequences and their evolution can be harnessed in bioengineering to build more robust synthetic strains and constructs. With a focus on overlapping protein-coding and RNA-coding genes, this Review examines their discovery, topology and biogenesis in the context of their genome biology. We highlight exciting new uses for sequence overlap to control translation, compress synthetic genetic constructs, and protect against mutation.

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.

Minigene as a Novel Regulatory Element in Toxin-Antitoxin Systems

The axe-txe type II toxin-antitoxin (TA) system is characterized by a complex and multilayered mode of gene expression regulation. Precise and tight control of this process is crucial to keep the toxin in an appropriate balance with the cognate antitoxin until its activation is needed for the cell. In this report, we provide evidence that a minigene encoded within the axe-txe operon influences translation of the Txe toxin. This is the first example to date of such a regulatory mechanism identified in the TA modules. Here, in a series of genetic studies, we employed translational reporter gene fusions to establish the molecular basis of this phenomenon. Our results show that translation of the two-codon mini-ORF displays an in cis mode of action, and positively affects the expression of txe, possibly by increasing its mRNA stability through protection from an endonuclease attack. Moreover, we established that the reading frame in which the two cistrons are encoded, as well as the distance between them, are critical parameters that affect the level of such regulation. In addition, by searching for two-codon ORFs we found sequences of several potential minigenes in the leader sequences of several other toxins belonging to the type II TA family. These findings suggest that this type of gene regulation may not only apply for the axe-txe cassette, but could be more widespread among other TA systems.

RNase III-mediated processing of a trans-acting bacterial sRNA and its cis-encoded antagonist

Bacterial small RNAs (sRNAs) are important post-transcriptional regulators in stress responses and virulence. They can be derived from an expanding list of genomic contexts, such as processing from parental transcripts by RNase E. The role of RNase III in sRNA biogenesis is less well understood despite its well-known roles in rRNA processing, RNA decay, and cleavage of sRNA-mRNA duplexes. Here, we show that RNase III processes a pair of cis-encoded sRNAs (CJnc190 and CJnc180) of the food-borne pathogen Campylobacter jejuni. While CJnc180 processing by RNase III requires CJnc190, RNase III processes CJnc190 independent of CJnc180 via cleavage of an intramolecular duplex. We also show that CJnc190 directly represses translation of the colonization factor PtmG by targeting a G-rich ribosome-binding site, and uncover that CJnc180 is a cis-acting antagonist of CJnc190, indirectly affecting ptmG regulation. Our study highlights a role for RNase III in sRNA biogenesis and adds cis-encoded RNAs to the expanding diversity of transcripts that can antagonize bacterial sRNAs.

Bacterial Chaperone Protein Hfq Facilitates the Annealing of Sponge RNAs to Small Regulatory RNAs

Bacterial small RNAs (sRNAs) in association with the chaperone protein Hfq regulate the expression of many target mRNAs. Since sRNAs' action is crucial to engendering a response to changing environmental conditions, their activity needs to be regulated. One such mechanism occurs at the post-transcriptional level and involves sponge RNAs, which sequester sRNAs affecting their regulatory output. Both types of RNAs were identified on Hfq, but it is not known how Hfq interacts with RNA sponges and stimulates their base-pairing with sRNAs. Here, we used biochemical methods to demonstrate that sponge RNAs resemble sRNAs by their structure and their modes of Hfq binding. Hfq facilitates the annealing of AgvB and 3'ETS^leuZ sponge RNAs to targeted sRNAs: GcvB and RybB, respectively, and each surface of the protein plays a particular role in the process. Moreover, we found that the efficiency of sponge RNA interactions with sRNAs can be improved; therefore, we propose that natural RNA sponges might not sequester sRNAs optimally.

Rho-dependent transcription termination: a revisionist view

Rho is a hexameric bacterial RNA helicase, which became a paradigm of factor-dependent transcription termination. The broadly accepted ("textbook") model posits a series of steps, wherein Rho first binds C-rich Rho utilization (rut) sites on nascent RNA, uses its ATP-dependent translocase activity to catch up with RNA polymerase (RNAP), and either pulls the transcript from the elongation complex or pushes RNAP forward, thus terminating transcription. However, this appealingly simple mechano-chemical model lacks a biological realism and is increasingly at odds with genetic and biochemical data. Here, we summarize recent structural and biochemical studies that have advanced our understanding of molecular details of RNA recognition, termination signaling, and RNAP inactivation in Rho-dependent transcription termination, rebalancing the view in favor of an alternative "allosteric" mechanism. In the revised model, Rho binds RNAP early in elongation assisted by the cofactors NusA and NusG, forming a pre-termination complex (PTC). The formation of PTC allows Rho to continuously sample nascent transcripts for a termination signal, which subsequently traps the elongation complex in an inactive state prior to its dissociation.

Riboregulation in bacteria: From general principles to novel mechanisms of the trp attenuator and its sRNA and peptide products

Gene expression strategies ensuring bacterial survival and competitiveness rely on cis- and trans-acting RNA-regulators (riboregulators). Among the cis-acting riboregulators are transcriptional and translational attenuators, and antisense RNAs (asRNAs). The trans-acting riboregulators are small RNAs (sRNAs) that bind proteins or base pairs with other RNAs. This classification is artificial since some regulatory RNAs act both in cis and in trans, or function in addition as small mRNAs. A prominent example is the archetypical, ribosome-dependent attenuator of tryptophan (Trp) biosynthesis genes. It responds by transcription attenuation to two signals, Trp availability and inhibition of translation, and gives rise to two trans-acting products, the attenuator sRNA rnTrpL and the leader peptide peTrpL. In Escherichia coli, rnTrpL links Trp availability to initiation of chromosome replication and in Sinorhizobium meliloti, it coordinates regulation of split tryptophan biosynthesis operons. Furthermore, in S. meliloti, peTrpL is involved in mRNA destabilization in response to antibiotic exposure. It forms two types of asRNA-containing, antibiotic-dependent ribonucleoprotein complexes (ARNPs), one of them changing the target specificity of rnTrpL. The posttranscriptional role of peTrpL indicates two emerging paradigms: (1) sRNA reprograming by small molecules and (2) direct involvement of antibiotics in regulatory RNPs. They broaden our view on RNA-based mechanisms and may inspire new approaches for studying, detecting, and using antibacterial compounds. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Kinetic modeling reveals additional regulation at co-transcriptional level by post-transcriptional sRNA regulators

Small RNAs (sRNAs) are important gene regulators in bacteria. Many sRNAs act post-transcriptionally by affecting translation and degradation of the target mRNAs upon base-pairing interactions. Here we present a general approach combining imaging and mathematical modeling to determine kinetic parameters at different levels of sRNA-mediated gene regulation that contribute to overall regulation efficacy. Our data reveal that certain sRNAs previously characterized as post-transcriptional regulators can regulate some targets co-transcriptionally, leading to a revised model that sRNA-mediated regulation can occur early in an mRNA’s lifetime, as soon as the sRNA binding site is transcribed. This co-transcriptional regulation is likely mediated by Rho-dependent termination when transcription-coupled translation is reduced upon sRNA binding. Our data also reveal several important kinetic steps that contribute to the differential regulation of mRNA targets by an sRNA. Particularly, binding of sRNA to the target mRNA may dictate the regulation hierarchy observed within an sRNA regulon.

Reyer et al. use fluorescent microscopy and kinetic modeling to find that two sRNAs canonically described as post-transcriptional regulators can regulate their targets co-transcriptionally and determine the in vivo kinetic parameters that dictate differential regulation efficiency.

Diversity and Versatility in Small RNA-Mediated Regulation in Bacterial Pathogens

Bacterial gene expression is under the control of a large set of molecules acting at multiple levels. In addition to the transcription factors (TFs) already known to be involved in global regulation of gene expression, small regulatory RNAs (sRNAs) are emerging as major players in gene regulatory networks, where they allow environmental adaptation and fitness. Developments in high-throughput screening have enabled their detection in the entire bacterial kingdom. These sRNAs influence a plethora of biological processes, including but not limited to outer membrane synthesis, metabolism, TF regulation, transcription termination, virulence, and antibiotic resistance and persistence. Almost always noncoding, they regulate target genes at the post-transcriptional level, usually through base-pair interactions with mRNAs, alone or with the help of dedicated chaperones. There is growing evidence that sRNA-mediated mechanisms of actions are far more diverse than initially thought, and that they go beyond the so-called cis- and trans-encoded classifications. These molecules can be derived and processed from 5' untranslated regions (UTRs), coding or non-coding sequences, and even from 3' UTRs. They usually act within the bacterial cytoplasm, but recent studies showed sRNAs in extracellular vesicles, where they influence host cell interactions. In this review, we highlight the various functions of sRNAs in bacterial pathogens, and focus on the increasing examples of widely diverse regulatory mechanisms that might compel us to reconsider what constitute the sRNA.

BB0562 is a nutritional virulence determinant with lipase activity important for Borrelia burgdorferi infection and survival in fatty acid deficient environments

The Lyme disease spirochete Borrelia burgdorferi relies on uptake of essential nutrients from its host environments for survival and infection. Therefore, nutrient acquisition mechanisms constitute key virulence properties of the pathogen, yet these mechanisms remain largely unknown. In vivo expression technology applied to B. burgdorferi (BbIVET) during mammalian infection identified gene bb0562, which encodes a hypothetical protein comprised of a conserved domain of unknown function, DUF3996. DUF3996 is also found across adjacent encoded hypothetical proteins BB0563 and BB0564, suggesting the possibility that the three proteins could be functionally related. Deletion of bb0562, bb0563 and bb0564 individually and together demonstrated that bb0562 alone was important for optimal disseminated infection in immunocompetent and immunocompromised mice by needle inoculation and tick bite transmission. Moreover, bb0562 promoted spirochete survival during the blood dissemination phase of infection. Gene bb0562 was also found to be important for spirochete growth in low serum media and the growth defect of Δbb0562 B. burgdorferi was rescued with the addition of various long chain fatty acids, particularly oleic acid. In mammals, fatty acids are primarily stored in fat droplets in the form of triglycerides. Strikingly, addition of glyceryl trioleate, the triglyceride form of oleic acid, to the low serum media did not rescue the growth defect of the mutant, suggesting bb0562 may be important for the release of fatty acids from triglycerides. Therefore, we searched for and identified two canonical GXSXG lipase motifs within BB0562, despite the lack of homology to known bacterial lipases. Purified BB0562 demonstrated lipolytic activity dependent on the catalytic serine residues within the two motifs. In sum, we have established that bb0562 is a novel nutritional virulence determinant, encoding a lipase that contributes to fatty acid scavenge for spirochete survival in environments deficient in free fatty acids including the mammalian host.

Borrelia burgdorferi, the causative agent of Lyme disease, has a small genome and lacks the ability to synthesize essential nutrients on its own as well as many of the virulence properties typical of bacterial pathogens that contribute to disease. The clinical manifestations of Lyme disease predominantly result from inflammation in response to the B. burgdorferi infection. Therefore, nutrient acquisition functions constitute key virulence factors for the pathogen. Fatty acids are critical components of B. burgdorferi membranes and lipoproteins, which the spirochete must scavenge from the host environment. Previously, through a genetic screen for B. burgdorferi genes that are expressed during mammalian infection we identified gene of unknown function, bb0562. Herein, we demonstrate that bb0562 encodes a lipase that plays a role in the release of free fatty acids from triglycerides. Furthermore, bb0562 contributes to B. burgdorferi survival and dissemination in the mammalian host. BB0562 is important for spirochete survival in environments low in free fatty acids thereby adding to B. burgdorferi’s arsenal of nutritional virulence determinants necessary for the pathogen to be maintained in the tick-mouse enzootic cycle and to cause disseminated disease.

Burning the Candle at Both Ends: Have Exoribonucleases Driven Divergence of Regulatory RNA Mechanisms in Bacteria?

Regulatory RNAs have emerged as ubiquitous gene regulators in all bacterial species studied to date. The combination of sequence-specific RNA interactions and malleable RNA structure has allowed regulatory RNA to adopt different mechanisms of gene regulation in a diversity of genetic backgrounds. In the model GammaproteobacteriaEscherichia coli and Salmonella, the regulatory RNA chaperone Hfq appears to play a global role in gene regulation, directly controlling ∼20 to 25% of the entire transcriptome. While the model FirmicutesBacillus subtilis and Staphylococcus aureus encode a Hfq homologue, its role has been significantly depreciated. These bacteria also have marked differences in RNA turnover. E. coli and Salmonella degrade RNA through internal endonucleolytic and 3′→5′ exonucleolytic cleavage that appears to allow transient accumulation of mRNA 3′ UTR cleavage fragments that contain stabilizing 3′ structures. In contrast, B. subtilis and S. aureus are able to exonucleolytically attack internally cleaved RNA from both the 5′ and 3′ ends, efficiently degrading mRNA 3′ UTR fragments. Here, we propose that the lack of 5′→3′ exoribonuclease activity in Gammaproteobacteria has allowed the accumulation of mRNA 3′ UTR ends as the “default” setting. This in turn may have provided a larger pool of unconstrained RNA sequences that has fueled the expansion of Hfq function and small RNA (sRNA) regulation in E. coli and Salmonella. Conversely, the exoribonuclease RNase J may be a significant barrier to the evolution of 3′ UTR sRNAs in B. subtilis and S. aureus that has limited the pool of RNA ligands available to Hfq and other sRNA chaperones, depreciating their function in these model Firmicutes.

Methyltransferase-directed orthogonal tagging and sequencing of miRNAs and bacterial small RNAs

Background

Targeted installation of designer chemical moieties on biopolymers provides an orthogonal means for their visualisation, manipulation and sequence analysis. Although high-throughput RNA sequencing is a widely used method for transcriptome analysis, certain steps, such as 3′ adapter ligation in strand-specific RNA sequencing, remain challenging due to structure- and sequence-related biases introduced by RNA ligases, leading to misrepresentation of particular RNA species. Here, we remedy this limitation by adapting two RNA 2′-O-methyltransferases from the Hen1 family for orthogonal chemo-enzymatic click tethering of a 3′ sequencing adapter that supports cDNA production by reverse transcription of the tagged RNA.

Results

We showed that the ssRNA-specific DmHen1 and dsRNA-specific AtHEN1 can be used to efficiently append an oligonucleotide adapter to the 3′ end of target RNA for sequencing library preparation. Using this new chemo-enzymatic approach, we identified miRNAs and prokaryotic small non-coding sRNAs in probiotic Lactobacillus casei BL23. We found that compared to a reference conventional RNA library preparation, methyltransferase-Directed Orthogonal Tagging and RNA sequencing, mDOT-seq, avoids misdetection of unspecific highly-structured RNA species, thus providing better accuracy in identifying the groups of transcripts analysed. Our results suggest that mDOT-seq has the potential to advance analysis of eukaryotic and prokaryotic ssRNAs.

Conclusions

Our findings provide a valuable resource for studies of the RNA-centred regulatory networks in Lactobacilli and pave the way to developing novel transcriptome and epitranscriptome profiling approaches in vitro and inside living cells. As RNA methyltransferases share the structure of the AdoMet-binding domain and several specific cofactor binding features, the basic principles of our approach could be easily translated to other AdoMet-dependent enzymes for the development of modification-specific RNA-seq techniques.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01053-w.

smORFer: a modular algorithm to detect small ORFs in prokaryotes

Emerging evidence places small proteins (≤50 amino acids) more centrally in physiological processes. Yet, their functional identification and the systematic genome annotation of their cognate small open-reading frames (smORFs) remains challenging both experimentally and computationally. Ribosome profiling or Ribo-Seq (that is a deep sequencing of ribosome-protected fragments) enables detecting of actively translated open-reading frames (ORFs) and empirical annotation of coding sequences (CDSs) using the in-register translation pattern that is characteristic for genuinely translating ribosomes. Multiple identifiers of ORFs that use the 3-nt periodicity in Ribo-Seq data sets have been successful in eukaryotic smORF annotation. They have difficulties evaluating prokaryotic genomes due to the unique architecture (e.g. polycistronic messages, overlapping ORFs, leaderless translation, non-canonical initiation etc.). Here, we present a new algorithm, smORFer, which performs with high accuracy in prokaryotic organisms in detecting putative smORFs. The unique feature of smORFer is that it uses an integrated approach and considers structural features of the genetic sequence along with in-frame translation and uses Fourier transform to convert these parameters into a measurable score to faithfully select smORFs. The algorithm is executed in a modular way, and dependent on the data available for a particular organism, different modules can be selected for smORF search.