High-resolution RNA 3'-ends mapping of bacterial Rho-dependent transcripts

Physiological roles of an Acinetobacter-specific σ factor

The Gram-negative pathogen Acinetobacter baumannii is considered an "urgent threat" to human health due to its propensity to become antibiotic resistant. Understanding the distinct regulatory paradigms used by A. baumannii to mitigate cellular stresses may uncover new therapeutic targets. Many γ-proteobacteria use the extracytoplasmic function (ECF) σ factor, RpoE, to invoke envelope homeostasis networks in response to stress. Acinetobacter species contain the poorly characterized ECF "SigAb"; however, it is unclear if SigAb has the same physiological role as RpoE. Here, we show that SigAb is a metal stress-responsive ECF that appears unique to Acinetobacter species and distinct from RpoE-like ECFs. We combine promoter mutagenesis, motif scanning, and chromatin immunoprecipitation-sequencing (ChIP-seq) to define the direct SigAb regulon, which consists of genes encoding SigAb itself, the stringent response mediator, RelA, and the uncharacterized small RNA, "SabS." However, RNA-seq of strains overexpressing SigAb revealed a large, indirect regulon containing hundreds of genes. Metal resistance genes are key elements of the indirect regulon, as CRISPRi knockdown of sigAb or sabS resulted in increased copper sensitivity and excess copper-induced SigAb-dependent transcription. Furthermore, we found that two uncharacterized genes in the sigAb operon, "aabA" and "aabB," have anti-SigAb activity. Finally, employing a targeted Tn-seq approach that uses CRISPR-associated transposons, we show that sigAb, aabA, and aabB are important for fitness even during optimal growth conditions. Our work reveals new physiological roles for SigAb and SabS, provides a novel approach for assessing gene fitness, and highlights the distinct regulatory architecture of A. baumannii.

IMPORTANCE

Acinetobacter baumannii is a hospital-acquired pathogen, and many strains are resistant to multiple antibiotics. Understanding how A. baumannii senses and responds to stress may uncover novel routes to treat infections. Here, we examine how the Acinetobacter-specific transcription factor, SigAb, mitigates stress. We find that SigAb directly regulates only a small number of genes, but indirectly controls hundreds of genes that have substantial impacts on cell physiology. We show that SigAb is required for maximal growth, even during optimal conditions, and is acutely required during growth in the presence of elevated copper. Given that copper toxicity plays roles in pathogenesis and on copper-containing surfaces in hospitals, we speculate that SigAb function may be important in clinically relevant contexts.

Rho-dependent termination and RNase E-mediated cleavage: dual pathways for RNA 3' end processing in polycistronic mRNA

"Pre-full-length" transcripts are produced at the end of the polycistronic galactose (gal) operon, 5' galE-galT-galK-galM 3', via Rho-dependent transcription termination (RDT) and -independent transcription termination. The 3' end of the full-length galETKM mRNA is acquired by exonucleolytic processing of the 3'-OH ends of the pre-full-length transcripts. However, the gal operon produces an mRNA termed galE whose 3' end forms approximately 120 nucleotides downstream of the galE stop codon, within the subsequent gene, galT, thereby establishing polarity in gene expression. In this study, we investigated the molecular processes that generate the 3' end of galE mRNA. We discovered that the 3' ends of pre-galE mRNA are produced in the middle of galT as a result of the combination of two separate molecular processes-one previously reported as RDT and the other as unreported RNase E-mediated transcript cleavage. The 3' ends of pre-galE mRNA undergo exonucleolytic processing to the 3' end of galE mRNA observed in vivo. A hairpin structure containing an 8 bp stem and a 4-nucleotide loop, located 5-10 nucleotides upstream of the 3' ends of galE mRNA, blocks exoribonuclease digestion and renders transcript stability. These findings demonstrate that RNase E-contrary to its general role in mRNA degradation-produces RNA 3' ends that regulate polarity in gene expression.IMPORTANCEThis study reports the findings of two molecular mechanisms that generate the 3' ends of pre-galE mRNA in the gal operon, viz., Rho-dependent transcription termination and RNase E-mediated cleavage. These 3' ends are subsequently processed to produce stable galE mRNA with a hairpin structure that prevents exoribonuclease degradation. This mechanism establishes gene expression polarity by generating the 3' end of galE mRNA within galT in contrast to the usual mRNA degradation role of RNase E. The study reveals a unique role of RNase E in mRNA processing and stability.

Rho-dependent termination: a bacterial evolutionary capacitor for stress resistance

Since the Modern Synthesis, interest has grown in resolving the "black box" between genotype and phenotype. Contained within this black box are highly plastic RNA and proteins with global effects on chromosome integrity and gene expression that serve as evolutionary capacitors - elements that enable the accumulation and buffering of genetic variation in normal conditions and reveal hidden genetic variation when induced by environmental stress. Discussion of evolutionary capacitors has primarily focused on eukaryotic translation factors and chaperones, such as Hsp90 and PSI+ prion. However, due to the coupling of transcription and translation in prokaryotes, transcription factors can be equally impactful in the modulation of gene expression and phenotypes. In this review, we discuss the prokaryotic transcription terminator Rho and how mutagenesis and plasticity of Rho influence epistasis, evolvability, and adaptation to stress in bacteria. We discuss the effects of variation in Rho generated by nature, laboratory mutagenesis, and experimental evolution; and how this variation is constrained or encouraged by Rho's extensive network of protein interactors. Exploring Rho's role as an evolutionary capacitor, along with identifying additional elements that can serve this function, can significantly advance our understanding of how organisms adapt to thrive in diverse environments.

BacTermFinder: a comprehensive and general bacterial terminator finder using a CNN ensemble

A terminator is a DNA region that ends the transcription process. Currently, multiple computational tools are available for predicting bacterial terminators. However, these methods are specialized for certain bacteria or terminator type (i.e. intrinsic or factor-dependent). In this work, we developed BacTermFinder using an ensemble of convolutional neural networks (CNNs) receiving as input four different representations of terminator sequences. To develop BacTermFinder, we collected roughly 41 000 bacterial terminators (intrinsic and factor-dependent) of 22 species with varying GC-content (from 28% to 71%) from published studies that used RNA-seq technologies. We evaluated BacTermFinder's performance on terminators of five bacterial species (not used for training BacTermFinder) and two archaeal species. BacTermFinder's performance was compared with that of four other bacterial terminator prediction tools. Based on our results, BacTermFinder outperforms all other four approaches in terms of average recall without increasing the number of false positives. Moreover, BacTermFinder identifies both types of terminators (intrinsic and factor-dependent) and generalizes to archaeal terminators. Additionally, we visualized the saliency map of the CNNs to gain insights on terminator motif per species. BacTermFinder is publicly available at https://github.com/BioinformaticsLabAtMUN/BacTermFinder.

Deciphering the Coupling State-Dependent Transcription Termination in the Escherichia coli Galactose Operon

The distance between the ribosome and the RNA polymerase active centers, known as the mRNA loop length, is crucial for transcription-translation coupling. Despite the existence of multiple expressomes with varying mRNA loop lengths, their in vivo roles remain largely unexplored. This study examines the mechanisms governing transcription termination in the Escherichia coli galactose operon, revealing a crucial role in the transcription and translation coupling state. The operon utilizes both Rho-independent and Rho-dependent terminators. Our findings demonstrate that long-loop coupled transcription-translation complexes preferentially terminate at the upstream Rho-independent terminator, while short-loop complexes bypass it, terminating at the downstream Rho-dependent terminator. The efficiency of the Rho-independent terminator is enhanced by an extended U-track, suggesting a novel mechanism to overcome ribosome inhibition. These results uncover a new regulatory layer in transcription termination, challenging the traditional view of this process as random and highlighting a predetermined mechanism based on the coupling state. We propose that tandem terminators may function as regulatory checkpoints under fluctuating ribosome-RNAP coupling conditions, which can occur due to specific cellular states or factors affecting ribosome or RNAP binding efficiency. This suggests a previously overlooked mechanism that could refine transcription termination choices and expand our understanding of transcription regulation.

The conserved noncoding RNA ModT coordinates growth and virulence in Clostridioides difficile

Bacterial noncoding RNAs fulfill a variety of cellular functions as catalysts, as scaffolds in protein complexes or as regulators of gene expression. They often exhibit complex tertiary structures that are a key determinant of their biochemical function. Here, we characterize the structured "raiA motif" RNA from Clostridioides difficile, which is conserved in more than 2,500 bacterial species from the phyla Bacillota and Actinomycetota. We show that its transcript abundance and stability in exponentially growing bacteria rivals that of ribosomal RNAs. Deletion of the "raiA motif" RNA is associated with delayed transition into stationary phase, and changes in stationary phase pathways such as spore formation, hence we rename it ModT (modulator of transition phase). Mechanistically, we show that ModT-mediated changes in cellular cyclic di-GMP levels are linked to the pronounced sporulation defect in the modT mutant. Importantly, we show that expression profiles and isoform patterns of ModT are conserved in Clostridium perfringens and Paeniclostridium sordellii, and that these orthologs can functionally complement ModT in C. difficile. Chemical structure probing of ModT in vivo reveals dynamic refolding and provides initial evidence for a potential association of ModT with proteins. In summary, our findings indicate that ModT fulfills a conserved role in regulating growth transitions in bacteria and provide a crucial step towards delineating its molecular mechanism.

PIPETS: a statistically informed, gene-annotation agnostic analysis method to study bacterial termination using 3'-end sequencing

BACKGROUND

Over the last decade the drop in short-read sequencing costs has allowed experimental techniques utilizing sequencing to address specific biological questions to proliferate, oftentimes outpacing standardized or effective analysis approaches for the data generated. There are growing amounts of bacterial 3'-end sequencing data, yet there is currently no commonly accepted analysis methodology for this datatype. Most data analysis approaches are somewhat ad hoc and, despite the presence of substantial signal within annotated genes, focus on genomic regions outside the annotated genes (e.g. 3' or 5' UTRs). Furthermore, the lack of consistent systematic analysis approaches, as well as the absence of genome-wide ground truth data, make it impossible to compare conclusions generated by different labs, using different organisms.

RESULTS

We present PIPETS, (Poisson Identification of PEaks from Term-Seq data), an R package available on Bioconductor that provides a novel analysis method for 3'-end sequencing data. PIPETS is a statistically informed, gene-annotation agnostic methodology. Across two different datasets from two different organisms, PIPETS identified significant 3'-end termination signal across a wider range of annotated genomic contexts than existing analysis approaches, suggesting that existing approaches may miss biologically relevant signal. Furthermore, assessment of the previously called 3'-end positions not captured by PIPETS showed that they were uniformly very low coverage.

CONCLUSIONS

PIPETS provides a broadly applicable platform to explore and analyze 3'-end sequencing data sets from across different organisms. It requires only the 3'-end sequencing data, and is broadly accessible to non-expert users.

PIPETS: A statistically informed, gene-annotation agnostic analysis method to study bacterial termination using 3'-end sequencing

Background

Over the last decade the drop in short-read sequencing costs has allowed experimental techniques utilizing sequencing to address specific biological questions to proliferate, oVentimes outpacing standardized or effective analysis approaches for the data generated. There are growing amounts of bacterial 3'-end sequencing data, yet there is currently no commonly accepted analysis methodology for this datatype. Most data analysis approaches are somewhat ad hoc and, despite the presence of substantial signal within annotated genes, focus on genomic regions outside the annotated genes (e.g. 3' or 5' UTRs). Furthermore, the lack of consistent systematic analysis approaches, as well as the absence of genome-wide ground truth data, make it impossible to compare conclusions generated by different labs, using different organisms.

Results

We present PIPETS, (Poisson Identification of PEaks from Term-Seq data), an R package available on Bioconductor that provides a novel analysis method for 3'-end sequencing data. PIPETS is a statistically informed, gene-annotation agnostic methodology. Across two different datasets from two different organisms, PIPETS identified significant 3'-end termination signal across a wider range of annotated genomic contexts than existing analysis approaches, suggesting that existing approaches may miss biologically relevant signal. Furthermore, assessment of the previously called 3'-end positions not captured by PIPETS showed that they were uniformly very low coverage.

Conclusions

PIPETS provides a broadly applicable placorm to explore and analyze 3'-end sequencing data sets from across different organisms. It requires only the 3'-end sequencing data, and is broadly accessible to non-expert users.

Rho-dependent transcriptional switches regulate the bacterial response to cold shock

Binding of the bacterial Rho helicase to nascent transcripts triggers Rho-dependent transcription termination (RDTT) in response to cellular signals that modulate mRNA structure and accessibility of Rho utilization (Rut) sites. Despite the impact of temperature on RNA structure, RDTT was never linked to the bacterial response to temperature shifts. We show that Rho is a central player in the cold-shock response (CSR), challenging the current view that CSR is primarily a posttranscriptional program. We identify Rut sites in 5'-untranslated regions of key CSR genes/operons (cspA, cspB, cspG, and nsrR-rnr-yjfHI) that trigger premature RDTT at 37°C but not at 15°C. High concentrations of RNA chaperone CspA or nucleotide changes in the cspA mRNA leader reduce RDTT efficiency, revealing how RNA restructuring directs Rho to activate CSR genes during the cold shock and to silence them during cold acclimation. These findings establish a paradigm for how RNA thermosensors can modulate gene expression.

Physiological Roles of an Acinetobacter-specific σ Factor

The Gram-negative pathogen Acinetobacter baumannii is considered an "urgent threat" to human health due to its propensity to become antibiotic resistant. Understanding the distinct regulatory paradigms used by A. baumannii to mitigate cellular stresses may uncover new therapeutic targets. Many γ-proteobacteria use the extracytoplasmic function (ECF) σ factor, RpoE, to invoke envelope homeostasis networks in response to stress. Acinetobacter species contain the poorly characterized ECF "SigAb;" however, it is unclear if SigAb has the same physiological role as RpoE. Here, we show that SigAb is a metal stress-responsive ECF that appears unique to Acinetobacter species and distinct from RpoE. We combine promoter mutagenesis, motif scanning, and ChIP-seq to define the direct SigAb regulon, which consists of sigAb itself, the stringent response mediator, relA, and the uncharacterized small RNA, "sabS." However, RNA-seq of strains overexpressing SigAb revealed a large, indirect regulon containing hundreds of genes. Metal resistance genes are key elements of the indirect regulon, as CRISPRi knockdown of sigAb or sabS resulted in increased copper sensitivity and excess copper induced SigAb-dependent transcription. Further, we found that two uncharacterized genes in the sigAb operon, "aabA" and "aabB", have anti-SigAb activity. Finally, employing a targeted Tn-seq approach that uses CRISPR-associated transposons, we show that sigAb, aabA, and aabB are important for fitness even during optimal growth conditions. Our work reveals new physiological roles for SigAb and SabS, provides a novel approach for assessing gene fitness, and highlights the distinct regulatory architecture of A. baumannii.

Co-transcriptional gene regulation in eukaryotes and prokaryotes

Many steps of RNA processing occur during transcription by RNA polymerases. Co-transcriptional activities are deemed commonplace in prokaryotes, in which the lack of membrane barriers allows mixing of all gene expression steps, from transcription to translation. In the past decade, an extraordinary level of coordination between transcription and RNA processing has emerged in eukaryotes. In this Review, we discuss recent developments in our understanding of co-transcriptional gene regulation in both eukaryotes and prokaryotes, comparing methodologies and mechanisms, and highlight striking parallels in how RNA polymerases interact with the machineries that act on nascent RNA. The development of RNA sequencing and imaging techniques that detect transient transcription and RNA processing intermediates has facilitated discoveries of transcription coordination with splicing, 3'-end cleavage and dynamic RNA folding and revealed physical contacts between processing machineries and RNA polymerases. Such studies indicate that intron retention in a given nascent transcript can prevent 3'-end cleavage and cause transcriptional readthrough, which is a hallmark of eukaryotic cellular stress responses. We also discuss how coordination between nascent RNA biogenesis and transcription drives fundamental aspects of gene expression in both prokaryotes and eukaryotes.

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.

Reporter Gene-Based qRT-PCR Assay for Rho-Dependent Termination In Vivo

In bacteria, the Rho protein mediates Rho-dependent termination (RDT) by identifying a non-specific cytosine-rich Rho utilization site on the newly synthesized RNA. As a result of RDT, downstream RNA transcription is reduced. Due to the bias in reverse transcription and PCR amplification, we could not identify the RDT site by directly measuring the amount of mRNA upstream and downstream of RDT sites. To overcome this difficulty, we employed a 77 bp reporter gene argX, (coding tRNAarg) from Brevibacterium albidum, and we transcriptionally fused it to the sequences to be assayed. We constructed a series of plasmids by combining a segment of the galactose (gal) operon sequences, both with and without the RDT regions at the ends of cistrons (galE, galT, and galM) upstream of argX. The RNA polymerase will transcribe the gal operon sequence and argX unless it encounters the RDT encoded by the inserted sequence. Since the quantitative real-time PCR (qRT-PCR) method detects the steady state following mRNA synthesis and degradation, we observed that tRNAarg is degraded at the same rate in these transcriptional fusion plasmids. Therefore, the amount of tRNAarg can directly reflect the mRNA synthesis. Using this approach, we were able to effectively assay the RDTs and Rho-independent termination (RIT) in the gal operon by quantifying the relative amount of tRNAarg using qRT-PCR analyses. The resultant RDT% for galET, galTK, and at the end of galM were 36, 26, and 63, individually. The resultant RIT% at the end of the gal operon is 33%. Our findings demonstrate that combining tRNAarg with qRT-PCR can directly measure RIT, RDT, or any other signal that attenuates transcription efficiencies in vivo, making it a useful tool for gene expression research.

Selective RNA Processing and Stabilization are Multi-Layer and Stoichiometric Regulators of Gene Expression in Escherichia coli

Selective RNA processing and stabilization (SRPS) facilitates the differential expression of multiple genes in polycistronic operons. However, how the coordinated actions of SRPS-related enzymes affect stoichiometric regulation remains unclear. In the present study, the first genome-wide targetome analysis is reported of these enzymes in Escherichia coli, at a single-nucleotide resolution. A strictly linear relationship is observed between the RNA pyrophosphohydrolase processing ratio and scores assigned to the first three nucleotides of the primary transcript. Stem-loops associated with PNPase targetomes exhibit a folding free energy that is negatively correlated with the termination ratio of PNPase at the 3' end. More than one-tenth of the RNase E processing sites in the 5'-untranslated regions（UTR） form different stem-loops that affect ribosome-binding and translation efficiency. The effectiveness of the SRPS elements is validated using a dual-fluorescence reporter system. The findings highlight a multi-layer and quantitative regulatory method for optimizing the stoichiometric expression of genes in bacteria and promoting the application of SRPS in synthetic biology.

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.

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.

Rho-dependent transcription termination is the dominant mechanism in Mycobacterium tuberculosis

Intrinsic and Rho-dependent transcription termination mechanisms regulate gene expression and recycle RNA polymerase in bacteria. Both the modes are well studied in Escherichia coli, and a few other organisms. The understanding of Rho function is limited in most other bacteria including mycobacteria. Here, we highlight the dominance of Rho-dependent termination in mycobacteria and validate Rho as a key regulatory factor. The lower abundance of intrinsic terminators, high cellular levels of Rho, and its genome-wide association with a majority of transcriptionally active genes indicate the pronounced role of Rho-mediated termination in Mycobacterium tuberculosis (Mtb). Rho modulates the termination of RNA synthesis for both protein-coding and stable RNA genes in Mtb. Concordantly, the depletion of Rho in mycobacteria impact its growth and enhances the transcription read-through at 3' ends of the transcription units. We demonstrate that MtbRho is catalytically active in the presence of RNA with varied secondary structures. These properties suggest an evolutionary adaptation of Rho as the efficient and preponderant mode of transcription termination in mycobacteria.

Identification of a Rho-Dependent Termination Site In Vivo Using Synthetic Small RNA

Rho promotes Rho-dependent termination (RDT) at the Rho-dependent terminator, producing a variable-length region without secondary structure at the 3' end of mRNA. Determining the exact RDT site in vivo is challenging, because the 3' end of mRNA is rapidly removed after RDT by 3'-to-5' exonuclease processing. Here, we applied synthetic small RNA (sysRNA) to identify the RDT region in vivo by exploiting its complementary base-pairing ability to target mRNA. Through the combined analyses of rapid amplification of cDNA 3' ends, primer extension, and capillary electrophoresis, we could precisely map and quantify mRNA 3' ends. We found that complementary double-stranded RNA (dsRNA) formed between sysRNA and mRNA was efficiently cleaved by RNase III in the middle of the dsRNA region. The formation of dsRNA appeared to protect the cleaved RNA 3' ends from rapid degradation by 3'-to-5' exonuclease, thereby stabilizing the mRNA 3' end. We further verified that the signal intensity at the 3' end was positively correlated with the amount of mRNA. By constructing a series of sysRNAs with close target sites and comparing the difference in signal intensity at the 3' end of wild-type and Rho-impaired strains, we finally identified a region of increased mRNA expression within the 21-bp range, which was determined as the RDT region. Our results demonstrated the ability to use sysRNA as a novel tool to identify RDT regions in vivo and expand the range of applications of sysRNA. IMPORTANCE sysRNA, which was formerly widely employed, has steadily lost popularity as more novel techniques for suppressing gene expression come into existence because of issues such as unstable inhibition effect and low inhibition efficiency. However, it remains an interesting topic as a regulatory tool due to its ease of design and low metabolic burden on cells. Here, for the first time, we discovered a new method to identify RDT regions in vivo using sysRNA. This new feature is important because since the discovery of the Rho protein in 1969, specific identification of RDT sites in vivo has been difficult due to the rapid processing of RNA 3' ends by exonucleases, and sysRNA might provide a new approach to address this challenge.

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.

The recognition of structured elements by a conserved groove distant from domains associated with catalysis is an essential determinant of RNase E

RNase E is an endoribonuclease found in many bacteria, including important human pathogens. Within Escherichia coli, it has been shown to have a major role in both the maturation of all classes of RNA involved in translation and the initiation of mRNA degradation. Thus, knowledge of the major determinants of RNase E cleavage is central to our understanding and manipulation of bacterial gene expression. We show here that the binding of RNase E to structured RNA elements is crucial for the processing of tRNA, can activate catalysis and may be important in mRNA degradation. The recognition of structured elements by RNase E is mediated by a recently discovered groove that is distant from the domains associated with catalysis. The functioning of this groove is shown here to be essential for E. coli cell viability and may represent a key point of evolutionary divergence from the paralogous RNase G family, which we show lack amino acid residues conserved within the RNA-binding groove of members of the RNase E family. Overall, this work provides new insights into the recognition and cleavage of RNA by RNase E and provides further understanding of the basis of RNase E essentiality in E. coli.

Comprehensive transcription terminator atlas for Bacillus subtilis

The transcriptome-wide contributions of Rho-dependent and intrinsic (Rho-independent) transcription termination mechanisms in bacteria are unclear. By sequencing released transcripts in a wild-type strain and strains containing deficiencies in NusA, NusG and/or Rho (10 strains), we produced an atlas of terminators for the model Gram-positive bacterium Bacillus subtilis. We found that NusA and NusG stimulate 77% and 19% of all intrinsic terminators, respectively, and that both proteins participate in Rho-dependent termination. We also show that Rho stimulates termination at 10% of the intrinsic terminators in vivo. We recapitulated Rho-stimulated intrinsic termination at 5 terminators in vitro and found that Rho requires the KOW domain of NusG to stimulate this process at one of these terminators. Computational analyses of our atlas using RNAstructure, MEME suite and DiffLogo, combined with in vitro transcription experiments, revealed that Rho stimulates intrinsic terminators with weak hairpins and/or U-rich tracts by remodelling the RNA upstream of the intrinsic terminator to prevent the formation of RNA structures that could otherwise compete with the terminator hairpin. We also identified 56 putative examples of 'hybrid Rho-dependent termination', wherein classical Rho-dependent termination occurs after readthrough of a Rho-stimulated intrinsic terminator.

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.

PNPase and RhlB Interact and Reduce the Cellular Availability of Oxidized RNA in Deinococcus radiodurans

8-Oxo-7,8-dihydroguanine (8-oxoG) is a major RNA modification caused by oxidative stresses and has been implicated in carcinogenesis, neurodegeneration, and aging. Several RNA-binding proteins have been shown to have a binding preference for 8-oxoG-modified RNA in eukaryotes and protect cells from oxidative stress. To date, polynucleotide phosphorylase (PNPase) is one of the most well-characterized proteins in bacteria that recognize 8-oxoG-modified RNA, but how PNPase cooperates with other proteins to process oxidized RNA is still unclear. Here, we use RNA affinity chromatography and mass spectrometry to search for proteins that preferably bind 8-oxoG-modified RNA in Deinococcus radiodurans, an extremophilic bacterium with extraordinary resistance to oxidative stresses. We identified four proteins that preferably bind to oxidized RNA: PNPase (DR_2063), DEAD box RNA helicase (DR_0335/RhlB), ribosomal protein S1 (DR_1983/RpsA), and transcriptional termination factor (DR_1338/Rho). Among these proteins, PNPase and RhlB exhibit high-affinity binding to 8-oxoG-modified RNA in a dose-independent manner. Deletions of PNPase and RhlB caused increased sensitivity of D. radiodurans to oxidative stress. We further showed that PNPase and RhlB specifically reduce the cellular availability of 8-oxoG-modified RNA but have no effect on oxidized DNA. Importantly, PNPase directly interacts with RhlB in D. radiodurans; however, no additional phenotypic effect was observed for the double deletion of pnp and rhlB compared to the single deletions. Overall, our findings suggest the roles of PNPase and RhlB in targeting 8-oxoG-modified RNAs and thereby constitute an important component of D. radiodurans resistance to oxidative stress. IMPORTANCE Oxidative RNA damage can be caused by oxidative stress, such as hydrogen peroxide, ionizing radiation, and antibiotic treatment. 8-oxo-7,8-dihydroguanine (8-oxoG), a major type of oxidized RNA, is highly mutagenic and participates in a variety of disease occurrences and development. Although several proteins have been identified to recognize 8-oxoG-modified RNA, the knowledge of how RNA oxidative damage is controlled largely remains unclear, especially in nonmodel organisms. In this study, we identified four RNA binding proteins that show higher binding affinity to 8-oxoG-modified RNA compared to unmodified RNA in the extremophilic bacterium Deinococcus radiodurans, which can endure high levels of oxidative stress. Two of the proteins, polynucleotide phosphorylase (PNPase) and DEAD-box RNA helicase (RhlB), interact with each other and reduce the cellular availability of 8-oxoG-modified RNA under oxidative stress. As such, this work contributes to our understanding of how RNA oxidation is influenced by RNA binding proteins in bacteria.

Factor-stimulated intrinsic termination: getting by with a little help from some friends

Transcription termination is known to occur via two mechanisms in bacteria, intrinsic termination (also frequently referred to as Rho-independent or factor-independent termination) and Rho-dependent termination. Based primarily on in vitro studies using Escherichia coli RNA polymerase, it was generally assumed that intrinsic termination and Rho-dependent termination are distinct mechanisms and that the signals required for intrinsic termination are present primarily within the nucleic acids. In this review, we detail recent findings from studies in Bacillus subtilis showing that intrinsic termination in this organism is highly stimulated by NusA, NusG, and even Rho. In NusA-stimulated intrinsic termination, NusA facilitates the formation of weak terminator hairpins and compensates for distal U-rich tract interruptions. In NusG-stimulated intrinsic termination, NusG stabilizes a sequence-dependent pause at the point of termination, which extends the time frame for RNA hairpins with weak terminal base pairs to form in either a NusA-stimulated or a NusA-independent fashion. In Rho-stimulated intrinsic termination, Rho prevents the formation of antiterminator-like RNA structures that could otherwise compete with the terminator hairpin. Combined, NusA, NusG, and Rho stimulate approximately 97% of all intrinsic terminators in B. subtilis. Thus, the general view that intrinsic termination is primarily a factor-independent process needs to be revised to account for recent findings. Moreover, the historical distinction between Rho-dependent and intrinsic termination is overly simplistic and needs to be modernized.

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.

In vivo regulation of bacterial Rho-dependent transcription termination by the nascent RNA

Bacterial Rho is a RNA-dependent ATPase that functions in the termination of DNA transcription. However, the in vivo nature of the bacterial Rho-dependent terminators, as well as the mechanism of the Rho-dependent termination process, are not fully understood. Here, we measured the in vivo termination efficiencies of 72 Rho-dependent terminators in E. coli by systematically performing qRT-PCR analyses of cDNA prepared from mid-log phase bacterial cultures. We found that these terminators exhibited a wide range of efficiencies, and many behaved differently in vivo compared to the predicted or experimentally determined efficiencies in vitro. Rho-utilization sites (rut sites) present in the RNA terminator sequences are characterized by the presence of C-rich/G-poor sequences, or C>G bubbles. We found that weaker terminators exhibited a robust correlation with the properties (size, length, density, etc.) of these C>G bubbles of their respective rut sites, while stronger terminators lack this correlation, suggesting a limited role of rut sequences in controlling in vivo termination efficiencies. We also found that in vivo termination efficiencies are dependent on the rates of ATP hydrolysis as well as Rho-translocation on the nascent RNA. We demonstrate that weaker terminators, in addition to having rut sites with diminished C>G bubble sizes, are dependent on the Rho-auxiliary factor, NusG, in vivo. From these results, we concluded that in vivo Rho-dependent termination follows a nascent RNA-dependent pathway, where Rho-translocation along the RNA is essential and rut sequences may recruit Rho in vivo, but Rho-rut binding strengths do not regulate termination efficiencies.

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.

Pervasive Transcription-coupled DNA repair in E. coli

Global Genomic Repair (GGR) and Transcription-Coupled Repair (TCR) have been viewed, respectively, as major and minor sub-pathways of the nucleotide excision repair (NER) process that removes bulky lesions from the genome. Here we applied a next generation sequencing assay, CPD-seq, in E. coli to measure the levels of cyclobutane pyrimidine dimer (CPD) lesions before, during, and after UV-induced genotoxic stress, and, therefore, to determine the rate of genomic recovery by NER at a single nucleotide resolution. We find that active transcription is necessary for the repair of not only the template strand (TS), but also the non-template strand (NTS), and that the bulk of TCR is independent of Mfd - a DNA translocase that is thought to be necessary and sufficient for TCR in bacteria. We further show that repair of both TS and NTS is enhanced by increased readthrough past Rho-dependent terminators. We demonstrate that UV-induced genotoxic stress promotes global antitermination so that TCR is more accessible to the antisense, intergenic, and other low transcribed regions. Overall, our data suggest that GGR and TCR are essentially the same process required for complete repair of the bacterial genome.

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.

System-Level Analysis of Transcriptional and Translational Regulatory Elements in Streptomyces griseus

Bacteria belonging to Streptomyces have the ability to produce a wide range of secondary metabolites through a shift from primary to secondary metabolism regulated by complex networks activated after vegetative growth terminates. Despite considerable effort to understand the regulatory elements governing gene expression related to primary and secondary metabolism in Streptomyces, system-level information remains limited. In this study, we integrated four multi-omics datasets from Streptomyces griseus NBRC 13350: RNA-seq, ribosome profiling, dRNA-seq, and Term-Seq, to analyze the regulatory elements of transcription and translation of differentially expressed genes during cell growth. With the functional enrichment of gene expression in different growth phases, one sigma factor regulon and four transcription factor regulons governing differential gene transcription patterns were found. In addition, the regulatory elements of transcription termination and post-transcriptional processing at transcript 3'-end positions were elucidated, including their conserved motifs, stem-loop RNA structures, and non-terminal locations within the polycistronic operons, and the potential regulatory elements of translation initiation and elongation such as 5'-UTR length, RNA structures at ribosome-bound sites, and codon usage were investigated. This comprehensive genetic information provides a foundational genetic resource for strain engineering to enhance secondary metabolite production in Streptomyces.

Genome-scale analysis of genetic regulatory elements in Streptomyces avermitilis MA-4680 using transcript boundary information

Background

The gram-positive bacterium, Streptomyces avermitilis, holds industrial importance as the producer of avermectin, a widely used anthelmintic agent, and a heterologous expression host of secondary metabolite-biosynthetic gene clusters. Despite its industrial importance, S. avermitilis’ genome organization and regulation of gene expression remain poorly understood. In this study, four different types of Next-Generation Sequencing techniques, including dRNA-Seq, Term-Seq, RNA-Seq and ribosome profiling, were applied to S. avermitilis to determine transcription units of S. avermitilis at a genome-wide level and elucidate regulatory elements for transcriptional and translational control of individual transcription units.

Result

By applying dRNA-Seq and Term-Seq to S. avermitilis MA-4680, a total of 2361 transcription start sites and 2017 transcript 3′-end positions were identified, respectively, leading to determination of 1601 transcription units encoded in S. avermitilis’ genome. Cataloguing the transcription units and integrated analysis of multiple high-throughput data types revealed the presence of diverse regulatory elements for gene expression, such as promoters, 5′-UTRs, terminators, 3′-UTRs and riboswitches. The conserved promoter motifs were identified from 2361 transcription start sites as 5′-TANNNT and 5′-BTGACN for the − 10 and − 35 elements, respectively. The − 35 element and spacer lengths between − 10 and − 35 elements were critical for transcriptional regulation of functionally distinct genes, suggesting the involvement of unique sigma factors. In addition, regulatory sequences recognized by antibiotic regulatory proteins were identified from the transcription start site information. Analysis of the 3′-end of RNA transcript revealed that stem structure formation is a major determinant for transcription termination of most transcription units.

Conclusions

The transcription unit architecture elucidated from the transcripts’ boundary information provides insights for unique genetic regulatory mechanisms of S. avermitilis. Our findings will elevate S. avermitilis’ potential as a production host for a diverse set of secondary metabolites.

Different Regulatory Modes of Synechocystis sp. PCC 6803 in Response to Photosynthesis Inhibitory Conditions

Cyanobacteria are promising industrial platforms owing to their ability to produce diverse natural secondary metabolites and nonnative value-added biochemicals from CO₂ and light. To fully utilize their industrial potency, it is critical to understand their photosynthetic efficiency under various environmental conditions. In this study, we elucidated the inhibitory mechanisms of photosynthesis under high-light and low-temperature stress conditions in the model cyanobacterium Synechocystis sp. PCC 6803. Under each stress condition, the transcript abundance and translation efficiency were measured using transcriptome sequencing (RNA-seq) and ribosome profiling, and the genome-wide transcription unit architecture was constructed by data integration of transcription start sites and transcript 3′-end positions obtained from differential RNA-seq and sequencing of 3′-ends (Term-seq), respectively. Our results suggested that the mode of photosynthesis inhibition differed between the two stress conditions; high light stress induced photodamage responses, while low temperature stress impaired the translation efficiency of photosynthesis-associated genes. In particular, poor translation of photosystem I resulted from ribosome stalling at the untranslated regions, affecting the overall photosynthetic yield under low temperature stress. Our comprehensive multiomics analysis with transcription unit architecture provides foundational information on photosynthesis for future industrial strain development.

IMPORTANCE Cyanobacteria are a compelling biochemical production platform for their ability to propagate using light and atmospheric CO₂ via photosynthesis. However, the engineering of strains is hampered by limited understanding of photosynthesis under diverse environmental conditions such as high-light and low-temperature stresses. Herein, we decipher the transcriptomic and translatomic responses of the photosynthetic efficiency to stress conditions using the integrative analysis of multiomic data generated by RNA-seq and ribosome profiling, respectively. Through the generated massive data, along with the guide of the genome-wide transcription unit architecture constructed by transcription start sites and transcript 3′-end positions, we identified the factors affecting photosynthesis at transcription, posttranscription, and translation levels. Importantly, the high-light stress induces photodamage responses, and the low-temperature stress cripples the translation efficiency of photosynthesis-associated genes. The resulting insights provide pivotal information for future cyanobacterial cell factories powered by the engineering toward robust photosynthesis ability.

Analysing the fitness cost of antibiotic resistance to identify targets for combination antimicrobials

Mutations in the rifampicin (Rif)-binding site of RNA polymerase (RNAP) confer antibiotic resistance and often have global effects on transcription that compromise fitness and stress tolerance of resistant mutants. We suggested that the non-essential genome, through its impact on the bacterial transcription cycle, may represent an untapped source of targets for combination antimicrobial therapies. Using transposon sequencing, we carried out a genome-wide analysis of fitness cost in a clinically common rpoB H526Y mutant. We find that genes whose products enable increased transcription elongation rates compound the fitness costs of resistance whereas genes whose products function in cell wall synthesis and division mitigate it. We validate our findings by showing that the cell wall synthesis and division defects of rpoB H526Y result from an increased transcription elongation rate that is further exacerbated by the activity of the uracil salvage pathway and unresponsiveness of the mutant RNAP to the alarmone ppGpp. We applied our findings to identify drugs that inhibit more readily rpoB H526Y and other Rif^R alleles from the same phenotypic class. Thus, genome-wide analysis of fitness cost of antibiotic-resistant mutants should expedite the discovery of new combination therapies and delineate cellular pathways that underlie the molecular mechanisms of cost.

The Complex Transcriptional Landscape of Magnetosome Gene Clusters in Magnetospirillum gryphiswaldense

Magnetosomes are complex membrane organelles synthesized by magnetotactic bacteria (MTB) for navigation in the Earth’s magnetic field. In the alphaproteobacterium Magnetospirillum gryphiswaldense, all steps of magnetosome formation are tightly controlled by >30 specific genes arranged in several gene clusters. However, the transcriptional organization of the magnetosome gene clusters has remained poorly understood. Here, by applying Cappable-seq and whole-transcriptome shotgun RNA sequencing, we show that mamGFDCop and feoAB1op are transcribed as single transcriptional units, whereas multiple transcription start sites (TSS) are present in mms6op, mamXYop, and the long (>16 kb) mamABop. Using a bioluminescence reporter assay and promoter knockouts, we demonstrate that most of the identified TSS originate from biologically meaningful promoters which mediate production of multiple transcripts and are functionally relevant for proper magnetosome biosynthesis. In addition, we identified a strong promoter in a large intergenic region within mamXYop, which likely drives transcription of a noncoding RNA important for gene expression in this operon. In summary, our data suggest a more complex transcriptional architecture of the magnetosome operons than previously recognized, which is largely conserved in other magnetotactic Magnetospirillum species and, thus, is likely fundamental for magnetosome biosynthesis in these organisms.

IMPORTANCE Magnetosomes have emerged as a model system to study prokaryotic organelles and a source of biocompatible magnetic nanoparticles for various biomedical applications. However, the lack of knowledge about the transcriptional organization of magnetosome gene clusters has severely impeded the engineering, manipulation, and transfer of this highly complex biosynthetic pathway into other organisms. Here, we provide a high-resolution image of the previously unappreciated transcriptional landscape of the magnetosome operons. Our findings are important for further unraveling the complex genetic framework of magnetosome biosynthesis. In addition, they will facilitate the rational reengineering of magnetic bacteria for improved bioproduction of tunable magnetic nanoparticles, as well as transplantation of magnetosome biosynthesis into foreign hosts by synthetic biology approaches. Overall, our study exemplifies how a genetically complex pathway is orchestrated at the transcriptional level to ensure the balanced expression of the numerous constituents required for the proper assembly of one of the most intricate prokaryotic organelles.

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.

Siblings or doppelgängers? Deciphering the evolution of structured cis-regulatory RNAs beyond homology

Structured cis-regulatory RNAs have evolved across all domains of life, highlighting the utility and plasticity of RNA as a regulatory molecule. Homologous RNA sequences and structures often have similar functions, but homology may also be deceiving. The challenges that derive from trying to assign function to structure and vice versa are not trivial. Bacterial riboswitches, viral and eukaryotic IRESes, CITEs, and 3′ UTR elements employ an array of mechanisms to exert their effects. Bioinformatic searches coupled with biochemical and functional validation have elucidated some shared and many unique ways cis-regulators are employed in mRNA transcripts. As cis-regulatory RNAs are resolved in greater detail, it is increasingly apparent that shared homology can mask the full spectrum of mRNA cis-regulator functional diversity. Furthermore, similar functions may be obscured by lack of obvious sequence similarity. Thus looking beyond homology is crucial for furthering our understanding of RNA-based regulation.

Multi-Omic Analyses Reveal Habitat Adaptation of Marine Cyanobacterium Synechocystis sp. PCC 7338

Cyanobacteria are considered as promising microbial cell factories producing a wide array of bio-products. Among them, Synechocystis sp. PCC 7338 has the advantage of growing in seawater, rather than requiring arable land or freshwater. Nonetheless, how this marine cyanobacterium grows under the high salt stress condition remains unknown. Here, we determined its complete genome sequence with the embedded regulatory elements and analyzed the transcriptional changes in response to a high-salt environment. Complete genome sequencing revealed a 3.70 mega base pair genome and three plasmids with a total of 3,589 genes annotated. Differential RNA-seq and Term-seq data aligned to the complete genome provided genome-wide information on genetic regulatory elements, including promoters, ribosome-binding sites, 5'- and 3'-untranslated regions, and terminators. Comparison with freshwater Synechocystis species revealed Synechocystis sp. PCC 7338 genome encodes additional genes, whose functions are related to ion channels to facilitate the adaptation to high salt and high osmotic pressure. Furthermore, a ferric uptake regulator binding motif was found in regulatory regions of various genes including SigF and the genes involved in energy metabolism, suggesting the iron-regulatory network is connected to not only the iron acquisition, but also response to high salt stress and photosynthesis. In addition, the transcriptomics analysis demonstrated a cyclic electron transport through photosystem I was actively used by the strain to satisfy the demand for ATP under high-salt environment. Our comprehensive analyses provide pivotal information to elucidate the genomic functions and regulations in Synechocystis sp. PCC 7338.

Elucidating the Regulatory Elements for Transcription Termination and Posttranscriptional Processing in the Streptomyces clavuligerus Genome

Identification of transcriptional regulatory elements in the GC-rich Streptomyces genome is essential for the production of novel biochemicals from secondary metabolite biosynthetic gene clusters (smBGCs). Despite many efforts to understand the regulation of transcription initiation in smBGCs, information on the regulation of transcription termination and posttranscriptional processing remains scarce. In this study, we identified the transcriptional regulatory elements in β-lactam antibiotic-producing Streptomyces clavuligerus ATCC 27064 by determining a total of 1,427 transcript 3'-end positions (TEPs) using the term-seq method. Termination of transcription was governed by three classes of TEPs, of which each displayed unique sequence features. The data integration with transcription start sites and transcriptome data generated 1,648 transcription units (TUs) and 610 transcription unit clusters (TUCs). TU architecture showed that the transcript abundance in TU isoforms of a TUC was potentially affected by the sequence context of their TEPs, suggesting that the regulatory elements of TEPs could control the transcription level in additional layers. We also identified TU features of a xenobiotic response element (XRE) family regulator and DUF397 domain-containing protein, particularly showing the abundance of bidirectional TEPs. Finally, we found that 189 noncoding TUs contained potential cis- and trans-regulatory elements that played a major role in regulating the 5' and 3' UTR. These findings highlight the role of transcriptional regulatory elements in transcription termination and posttranscriptional processing in Streptomyces sp.IMPORTANCE Streptomyces sp. is a great source of bioactive secondary metabolites, including antibiotics, antifungal agents, antiparasitic agents, immunosuppressant compounds, and other drugs. Secondary metabolites are synthesized via multistep conversions of the precursor molecules from primary metabolism, governed by multicomplex enzymes from secondary metabolite biosynthetic gene clusters. As their production is closely related with the growth phase and dynamic cellular status in response to various intra- and extracellular signals, complex regulatory systems tightly control the gene expressions related to secondary metabolism. In this study, we determined genome-wide transcript 3'-end positions and transcription units in the β-lactam antibiotic producer Streptomyces clavuligerus ATCC 27064 to elucidate the transcriptional regulatory elements in transcription termination and posttranscriptional processing by integration of multiomics data. These unique features, such as transcript 3'-end sequence, potential riboregulators, and potential 3'-untranslated region (UTR) cis-regulatory elements, can be potentially used to design engineering tools that can regulate the transcript abundance of genes for enhancing secondary metabolite production.

Engineering Transcriptional Interference through RNA Polymerase Processivity Control

Antisense transcription is widespread in all kingdoms of life and has been shown to influence gene expression through transcriptional interference (TI), a phenomenon in which one transcriptional process negatively influences another in cis. The processivity, or uninterrupted transcription, of an RNA polymerase (RNAP) is closely tied to levels of antisense transcription in bacterial genomes, but its influence on TI, while likely important, is not well-characterized. Here, we show that TI can be tuned through processivity control via three distinct antitermination strategies: the antibiotic bicyclomycin, phage protein Psu, and ribosome-RNAP coupling. We apply these methods toward TI and tune ribosome-RNAP coupling to produce 38-fold transcription-level gene repression due to both RNAP collisions and antisense RNA interference. We then couple protein roadblock and TI to design minimal genetic NAND and NOR logic gates. Together, these results show the importance of processivity control for strong TI and demonstrate TI's potential for synthetic biology.

NusG is an intrinsic transcription termination factor that stimulates motility and coordinates gene expression with NusA

NusA and NusG are transcription factors that stimulate RNA polymerase pausing in Bacillus subtilis. While NusA was known to function as an intrinsic termination factor in B. subtilis, the role of NusG in this process was unknown. To examine the individual and combinatorial roles that NusA and NusG play in intrinsic termination, Term-seq was conducted in wild type, NusA depletion, ΔnusG, and NusA depletion ΔnusG strains. We determined that NusG functions as an intrinsic termination factor that works alone and cooperatively with NusA to facilitate termination at 88% of the 1400 identified intrinsic terminators. Our results indicate that NusG stimulates a sequence-specific pause that assists in the completion of suboptimal terminator hairpins with weak terminal A-U and G-U base pairs at the bottom of the stem. Loss of NusA and NusG leads to global misregulation of gene expression and loss of NusG results in flagella and swimming motility defects.

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

Many bacterial genes are regulated by RNA elements in their 5´ untranslated regions (UTRs). However, the full complement of these elements is not known even in the model bacterium Escherichia coli. Using complementary RNA-sequencing approaches, we detected large numbers of 3´ ends in 5´ UTRs and open reading frames (ORFs), suggesting extensive regulation by premature transcription termination. We documented regulation for multiple transcripts, including spermidine induction involving Rho and translation of an upstream ORF for an mRNA encoding a spermidine efflux pump. In addition to discovering novel sites of regulation, we detected short, stable RNA fragments derived from 5´ UTRs and sequences internal to ORFs. Characterization of three of these transcripts, including an RNA internal to an essential cell division gene, revealed that they have independent functions as sRNA sponges. Thus, these data uncover an abundance of cis- and trans-acting RNA regulators in bacterial 5´ UTRs and internal to ORFs.

In most organisms, specific segments of a cell’s genetic information are copied to form single-stranded molecules of various sizes and purposes. Each of these RNA molecules, as they are known, is constructed as a chain that starts at the 5´ end and terminates at the 3´ end.

Certain RNAs carry the information present in a gene, which provides the instructions that a cell needs to build proteins. Some, however, are ‘non-coding’ and instead act to fine-tune the activity of other RNAs. These regulatory RNAs can be separate from the RNAs they control, or they can be embedded in the very sequences they regulate; new evidence also shows that certain regulatory RNAs can act in both ways.

Many regulatory RNAs are yet to be catalogued, even in simple, well-studied species such as the bacterium Escherichia coli. Here, Adams et al. aimed to better characterize the regulatory RNAs present in E. coli by mapping out the 3´ ends of every RNA molecule in the bacterium.

This revealed many new regulatory RNAs and offered insights into where these sequences are located. For instance, the results show that several of these RNAs were embedded within RNA produced from larger genes. Some were nested in coding RNAs, and were parts of a longer RNA sequence that is adjacent to the protein coding segment. Others, however, were present within the instructions that code for a protein.

The work by Adams et al. reveals that regulatory RNAs can be located in unexpected places, and provides a method for identifying them. This can be applied to other types of bacteria, in particular in species with few known RNA regulators.

Bacterial 3'UTRs: A Useful Resource in Post-transcriptional Regulation

Bacterial messenger RNAs (mRNAs) are composed of 5′ and 3′ untranslated regions (UTRs) that flank the coding sequences (CDSs). In eukaryotes, 3′UTRs play key roles in post-transcriptional regulatory mechanisms. Shortening or deregulation of these regions is associated with diseases such as cancer and metabolic disorders. Comparatively, little is known about the functions of 3′UTRs in bacteria. Over the past few years, 3′UTRs have emerged as important players in the regulation of relevant bacterial processes such as virulence, iron metabolism, and biofilm formation. This MiniReview is an update for the different 3′UTR-mediated mechanisms that regulate gene expression in bacteria. Some of these include 3′UTRs that interact with the 5′UTR of the same transcript to modulate translation, 3′UTRs that are targeted by specific ribonucleases, RNA-binding proteins and small RNAs (sRNAs), and 3′UTRs that act as reservoirs of trans-acting sRNAs, among others. In addition, recent findings regarding a differential evolution of bacterial 3′UTRs and its impact in the species-specific expression of orthologous genes are also discussed.

Systems and synthetic biology to elucidate secondary metabolite biosynthetic gene clusters encoded in Streptomyces genomes

Covering: 2010 to 2020 Over the last few decades, Streptomyces have been extensively investigated for their ability to produce diverse bioactive secondary metabolites. Recent advances in Streptomyces research have been largely supported by improvements in high-throughput technology 'omics'. From genomics, numerous secondary metabolite biosynthetic gene clusters were predicted, increasing their genomic potential for novel bioactive compound discovery. Additional omics, including transcriptomics, translatomics, interactomics, proteomics and metabolomics, have been applied to obtain a system-level understanding spanning entire bioprocesses of Streptomyces, revealing highly interconnected and multi-layered regulatory networks for secondary metabolism. The comprehensive understanding derived from this systematic information accelerates the rational engineering of Streptomyces to enhance secondary metabolite production, integrated with the exploitation of the highly efficient 'Design-Build-Test-Learn' cycle in synthetic biology. In this review, we describe the current status of omics applications in Streptomyces research to better understand the organism and exploit its genetic potential for higher production of valuable secondary metabolites and novel secondary metabolite discovery.

Genome-scale determination of 5´ and 3´ boundaries of RNA transcripts in Streptomyces genomes

Streptomyces species are gram-positive bacteria with GC-rich linear genomes and they serve as dominant reservoirs for producing clinically and industrially important secondary metabolites. Genome mining of Streptomyces revealed that each Streptomyces species typically encodes 20–50 secondary metabolite biosynthetic gene clusters (smBGCs), emphasizing their potential for novel compound discovery. Unfortunately, most of smBGCs are uncharacterized in terms of their products and regulation since they are silent under laboratory culture conditions. To translate the genomic potential of Streptomyces to practical applications, it is essential to understand the complex regulation of smBGC expression and to identify the underlying regulatory elements. To progress towards these goals, we applied two Next-Generation Sequencing methods, dRNA-Seq and Term-Seq, to industrially relevant Streptomyces species to reveal the 5´ and 3´ boundaries of RNA transcripts on a genome scale. This data provides a fundamental resource to aid our understanding of Streptomyces’ regulation of smBGC expression and to enhance their potential for secondary metabolite synthesis.

Measurement(s)	5´-ends of transcripts • 3´-ends of transcripts • RNA • TSS • transcription_termination_signal
Technology Type(s)	dRNA-Seq • Term-Seq • RNA sequencing
Factor Type(s)	Streptomyces growth phase
Sample Characteristic - Organism	Streptomyces avermitilis • Streptomyces clavuligerus • Streptomyces coelicolor • Streptomyces griseus • Streptomyces lividans • Streptomyces tsukubensis • Streptomyces venezuelae

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.13259393

Pre-termination Transcription Complex: Structure and Function

Rho is a general transcription termination factor playing essential roles in RNA polymerase (RNAP) recycling, gene regulation, and genomic stability in most bacteria. Traditional models of transcription termination postulate that hexameric Rho loads onto RNA prior to contacting RNAP and then translocates along the transcript in pursuit of the moving RNAP to pull RNA from it. Here, we report the cryoelectron microscopy (cryo-EM) structures of two termination process intermediates. Prior to interacting with RNA, Rho forms a specific "pre-termination complex" (PTC) with RNAP and elongation factors NusA and NusG, which stabilize the PTC. RNA exiting RNAP interacts with NusA before entering the central channel of Rho from the distal C-terminal side of the ring. We map the principal interactions in the PTC and demonstrate their critical role in termination. Our results support a mechanism in which the formation of a persistent PTC is a prerequisite for termination.

Interplay between Regulatory RNAs and Signal Transduction Systems during Bacterial Infection

The ability of pathogenic bacteria to stably infect the host depends on their capacity to respond and adapt to the host environment and on the efficiency of their defensive mechanisms. Bacterial envelope provides a physical barrier protecting against environmental threats. It also constitutes an important sensory interface where numerous sensing systems are located. Signal transduction systems include Two-Component Systems (TCSs) and alternative sigma factors. These systems are able to sense and respond to the ever-changing environment inside the host, altering the bacterial transcriptome to mitigate the impact of the stress. The regulatory networks associated with signal transduction systems comprise small regulatory RNAs (sRNAs) that can be directly involved in the expression of virulence factors. The aim of this review is to describe the importance of TCS- and alternative sigma factor-associated sRNAs in human pathogens during infection. The currently available genome-wide approaches for studies of TCS-regulated sRNAs will be discussed. The differences in the signal transduction mediated by TCSs between bacteria and higher eukaryotes and the specificity of regulatory RNAs for their targets make them appealing targets for discovery of new strategies to fight against multi-resistant bacteria.