Global analysis of mRNA decay intermediates in Bacillus subtilis wild-type and polynucleotide phosphorylase-deletion strains

Strategies for improving fengycin production: a review

Fengycin is an important member of the lipopeptide family with a wide range of applications in the agricultural, food, medical and cosmetic industries. However, its commercial application is severely hindered by low productivity and high cost. Therefore, numerous studies have been devoted to improving the production of fengycin. We summarize these studies in this review with the aim of providing a reference and guidance for future researchers. This review begins with an overview of the synthesis mechanism of fengycin via the non-ribosomal peptide synthetases (NRPS), and then delves into the strategies for improving the fengycin production in recent years. These strategies mainly include fermentation optimization and metabolic engineering, and the metabolic engineering encompasses enhancement of precursor supply, application of regulatory factors, promoter engineering, and application of genome-engineering (genome shuffling and genome-scale metabolic network model). Finally, we conclude this review with a prospect of fengycin production.

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.

The previously uncharacterized RnpM (YlxR) protein modulates the activity of ribonuclease P in Bacillus subtilis in vitro

Even though Bacillus subtilis is one of the most studied organisms, no function has been identified for about 20% of its proteins. Among these unknown proteins are several RNA- and ribosome-binding proteins suggesting that they exert functions in cellular information processing. In this work, we have investigated the RNA-binding protein YlxR. This protein is widely conserved in bacteria and strongly constitutively expressed in B. subtilis suggesting an important function. We have identified the RNA subunit of the essential RNase P as the binding partner of YlxR. The main activity of RNase P is the processing of 5' ends of pre-tRNAs. In vitro processing assays demonstrated that the presence of YlxR results in reduced RNase P activity. Chemical cross-linking studies followed by in silico docking analysis and experiments with site-directed mutant proteins suggest that YlxR binds to the region of the RNase P RNA that is important for binding and cleavage of the pre-tRNA substrate. We conclude that the YlxR protein is a novel interaction partner of the RNA subunit of RNase P that serves to finetune RNase P activity to ensure appropriate amounts of mature tRNAs for translation. We rename the YlxR protein RnpM for RNase P modulator.

A comprehensive study of the interactions in the B. subtilis degradosome with special emphasis on the role of the small proteins SR1P and SR7P

Here, we employ coelution experiments and far-western blotting to identify stable interactions between the main components of the B. subtilis degradosome and the small proteins SR1P and SR7P. Our data indicate that B. subtilis has a degradosome comprising at least RNases Y and PnpA, enolase, phosphofructokinase, glycerol-3-phosphate dehydrogenase GapA, and helicase CshA that can be co-purified without cross-linking. All interactions were corroborated by far-western blotting with proteins purified from E. coli. Previously, we discovered that stress-induced SR7P binds enolase to enhance its interaction with and activity of enolase-bound RNase Y (RnY), while SR1P transcribed under gluconeogenic conditions interacts with GapA to stimulate its interaction with and the activity of RnjA (RnjA). We show that SR1P can directly bind RnjA, RnY, and PnpA independently of GapA, whereas SR7P only interacts with enolase. Northern blotting suggests that the degradation of individual RNAs in B. subtilis under gluconeogenic or stress conditions depends on either RnjA or RnY alone or on RnjA-SR1P, RnY-SR1P, or RnY-Eno. In vitro degradation assays with RnY or RnjA substrates corroborate the in vivo role of SR1P. Currently, it is unknown which substrate property is decisive for the utilization of one of the complexes.

Bacillus subtilis NrnB is expressed during sporulation and acts as a unique 3'-5' exonuclease

All cells employ a combination of endo- and exoribonucleases to degrade long RNA polymers to fragments 2-5 nucleotides in length. These short RNA fragments are processed to monoribonucleotides by nanoRNases. Genetic depletion of nanoRNases has been shown to increase abundance of short RNAs. This deleteriously affects viability, virulence, and fitness, indicating that short RNAs are a metabolic burden. Previously, we provided evidence that NrnA is the housekeeping nanoRNase for Bacillus subtilis. Herein, we investigate the biological and biochemical functions of the evolutionarily related protein, B. subtilis NrnB (NrnBBs). These experiments show that NrnB is surprisingly different from NrnA. While NrnA acts at the 5' terminus of RNA substrates, NrnB acts at the 3' terminus. Additionally, NrnA is expressed constitutively under standard growth conditions, yet NrnB is selectively expressed during endospore formation. Furthermore, NrnA processes only short RNAs, while NrnB unexpectedly processes both short RNAs and longer RNAs. Indeed, inducible expression of NrnB can even complement the loss of the known global 3'-5' exoribonucleases, indicating that it acts as a general exonuclease. Together, these data demonstrate that NrnB proteins, which are widely found in Firmicutes, Epsilonproteobacteria and Archaea, are fundamentally different than NrnA proteins and may be used for specialized purposes.

Processing and decay of 6S-1 and 6S-2 RNAs in Bacillus subtilis

Noncoding 6S RNAs regulate transcription by binding to the active site of bacterial RNA polymerase holoenzymes. Processing and decay of 6S-1 and 6S-2 RNA were investigated in Bacillus subtilis by northern blot and RNA-seq analyses using different RNase knockout strains, as well as by in vitro processing assays. For both 6S RNA paralogs, we identified a key-but mechanistically different-role of RNase J1. RNase J1 catalyzes 5'-end maturation of 6S-1 RNA, yet relatively inefficient and possibly via the enzyme's "sliding endonuclease" activity. 5'-end maturation has no detectable effect on 6S-1 RNA function, but rather regulates its decay: The generated 5'-monophosphate on matured 6S-1 RNA propels endonucleolytic cleavage in its apical loop region. The major 6S-2 RNA degradation pathway is initiated by endonucleolytic cleavage in the 5'-central bubble to trigger 5'-to-3'-exoribonucleolytic degradation of the downstream fragment by RNase J1. The four 3'-exonucleases of B. subtilis-RNase R, PNPase, YhaM, and particularly RNase PH-are involved in 3'-end trimming of both 6S RNAs, degradation of 6S-1 RNA fragments, and decay of abortive transcripts (so-called product RNAs, ∼14 nt in length) synthesized on 6S-1 RNA during outgrowth from stationary phase. In the case of the growth-retarded RNase Y deletion strain, we were unable to infer a specific role of RNase Y in 6S RNA decay. Yet, a participation of RNase Y in 6S RNA decay still remains possible, as evidence for such a function may have been obscured by overlapping substrate specificities of RNase Y, RNase J1, and RNase J2.

Identification of Genes Required for Swarming Motility in Bacillus subtilis Using Transposon Mutagenesis and High-Throughput Sequencing (TnSeq)

Bacillus subtilis exhibits swarming motility, a flagellar-mediated form of surface motility. Here, we use transposon mutagenesis and sequencing (TnSeq) to perform a high-throughput screen for candidate genes required for swarming. The TnSeq approach identified all of the known genes required for flagellar biosynthesis and nearly all of the previously reported regulators that promote swarming. Moreover, we identified an additional 36 genes that improve swarming and validated them individually. Among these, two mutants with severe defects were recovered, including fliT, required for flagellar biosynthesis, and a gene of unknown function, yolB, whose defect could not be attributed to a lack of flagella. In addition to discovering additional genes required for B. subtilis swarming, our work validates TnSeq as a powerful approach for comprehensively identifying genes important for nonessential processes such as colony expansion on plates. IMPORTANCE In TnSeq, transposons are randomly inserted throughout the chromosome at a population level, but insertions that disrupt genes of essential function cause strains that carry them to fall out of the population and appear underrepresented at the sequence level. Here, we apply TnSeq to the nonessential phenotype of motility in B. subtilis and spatially select for cells proficient in swarming. We find that insertions in nearly all genes previously identified as required for swarming are underrepresented in TnSeq analysis, and we identify 36 additional genes that enhance swarming. We demonstrate that TnSeq is a powerful tool for the genetic analysis of motility and likely other nonlethal screens for which enrichment is strong.

Analysis of mRNA Decay Intermediates in Bacillus subtilis 3' Exoribonuclease and RNA Helicase Mutant Strains

The Bacillus subtilis genome encodes four 3' exoribonucleases: polynucleotide phosphorylase (PNPase), RNase R, RNase PH, and YhaM. Previous work showed that PNPase, encoded by the pnpA gene, is the major 3' exonuclease involved in mRNA turnover; in a pnpA deletion strain, numerous mRNA decay intermediates accumulate. Whether B. subtilis mRNA decay occurs in the context of a degradosome complex is controversial. In this study, global mapping of mRNA decay intermediate 3' ends within coding sequences was performed in strains that were either deleted for or had an inactivating point mutation in the pnpA gene. The patterns of 3'-end accumulation in these strains were highly similar, which may have implications for the role of a degradosome in mRNA decay. A comparison with mapped 3' ends in a strain lacking CshA, the major RNA helicase, indicated that many mRNAs require both PNPase and CshA for efficient decay. Transcriptome sequencing (RNA-seq) analysis of strains lacking RNase R suggested that this enzyme did not play a major role in mRNA turnover in the wild-type strain. Strains were constructed that contained only one of the four known 3' exoribonucleases. When RNase R was the only 3' exonuclease present, it was able to degrade a model mRNA efficiently, showing processive decay even through a strong stem-loop structure that inhibits PNPase processivity. Strains containing only RNase PH or only YhaM were also insensitive to this RNA secondary structure, suggesting the existence of another, as-yet-unidentified, 3' exoribonuclease. IMPORTANCE The ability to rapidly change bacterial gene expression programs in response to environmental conditions is highly dependent on the efficient turnover of mRNA. While much is known about the regulation of gene expression at the transcriptional and translational levels, much less is known about the intermediate step of mRNA decay. Here, we mapped the 3' ends of mRNA decay intermediates in strains that were missing the major 3' exoribonuclease PNPase or the RNA helicase CshA. We also assessed the roles of three other B. subtilis 3' exonucleases in the mRNA decay process. The data confirm the primary role of PNPase in mRNA turnover and suggest the involvement of one or more unknown RNases.

Discovery and initial characterization of YloC, a novel endoribonuclease in Bacillus subtilis

The Bacillus subtilis genome is predicted to encode numerous ribonucleases, including four 3' exoribonucleases that have been characterized to some extent. A strain containing gene knockouts of all four known 3' exoribonucleases is viable, suggesting that one or more additional RNases remain to be discovered. A protein extract from the quadruple RNase mutant strain was fractionated and RNase activity was followed, resulting in the identification of an enzyme activity catalyzed by the YloC protein. YloC is an endoribonuclease and is a member of the highly conserved "YicC family" of proteins that is widespread in bacteria. YloC is a metal-dependent enzyme that catalyzes the cleavage of single-stranded RNA, preferentially at U residues, and exists in an oligomeric form, most likely a hexamer. As such, YloC shares some characteristics with the SARS-CoV Nsp15 endoribonuclease. While the in vivo function of YloC in B. subtilis is yet to be determined, YloC was found to act similarly to YicC in an Escherichia coli in vivo assay that assesses decay of the small RNA, RyhB. Thus, YloC may play a role in small RNA regulation.

Pivotal Roles for Ribonucleases in Streptococcus pneumoniae Pathogenesis

RNases perform indispensable functions in regulating gene expression in many bacterial pathogens by processing and/or degrading RNAs. Despite the pivotal role of RNases in regulating bacterial virulence factors, the functions of RNases have not yet been studied in the major human respiratory pathogen Streptococcus pneumoniae (pneumococcus). Here, we sought to determine the impact of two conserved RNases, the endoribonuclease RNase Y and exoribonuclease polynucleotide phosphorylase (PNPase), on the physiology and virulence of S. pneumoniae serotype 2 strain D39. We report that RNase Y and PNPase are essential for pneumococcal pathogenesis, as both deletion mutants showed strong attenuation of virulence in murine models of invasive pneumonia. Genome-wide transcriptomic analysis revealed that the abundances of nearly 200 mRNA transcripts were significantly increased, whereas those of several pneumococcal small regulatory RNAs (sRNAs), including the Ccn (CiaR-controlled noncoding RNA) sRNAs, were altered in the Δrny mutant relative to the wild-type strain. Additionally, lack of RNase Y resulted in pleiotropic phenotypes that included defects in pneumococcal cell morphology and growth in vitro. In contrast, Δpnp mutants showed no growth defect in vitro but differentially expressed a total of 40 transcripts, including the tryptophan biosynthesis operon genes and numerous 5' cis-acting regulatory RNAs, a majority of which were previously shown to impact pneumococcal disease progression in mice using the serotype 4 strain TIGR4. Together, our data suggest that RNase Y exerts a global impact on pneumococcal physiology, while PNPase mediates virulence phenotypes, likely through sRNA regulation. IMPORTANCE Streptococcus pneumoniae is a notorious human pathogen that adapts to conditions in distinct host tissues and responds to host cell interactions by adjusting gene expression. RNases are key players that modulate gene expression by mediating the turnover of regulatory and protein-coding transcripts. Here, we characterized two highly conserved RNases, RNase Y and PNPase, and evaluated their impact on the S. pneumoniae transcriptome for the first time. We show that PNPase influences the levels of a narrow set of mRNAs but a large number of regulatory RNAs primarily implicated in virulence control, whereas RNase Y has a more sweeping effect on gene expression, altering levels of transcripts involved in diverse cellular processes, including cell division, metabolism, stress response, and virulence. This study further reveals that RNase Y regulates expression of genes governing competence by mediating the turnover of CiaR-controlled noncoding (Ccn) sRNAs.

Gene expression profile analysis and target gene discovery of Mycobacterium tuberculosis biofilm

Tuberculosis (TB) caused by Mycobacterium tuberculosis (M. tuberculosis) is a fatal infectious disease to human health, and the drug tolerance and immune evasion of M. tuberculosis were reported to be related to its biofilm formation; however, the difficulty of M. tuberculosis biofilm culture and its unknown global mechanism impede its further research. Here, we developed a modified in vitro M. tuberculosis biofilm model with shorter culture time. Then we used Illumina RNA-seq technology to determine the global gene expression profile of M. tuberculosis H37Rv biofilms. Over 437 genes are expressed at significantly different levels in biofilm cells than in planktonic cells; among them, 153 were downregulated and 284 were upregulated. Go enrichment analysis and KEGG pathway analysis showed that genes involved in biosynthesis and metabolism of sulfur metabolism, steroid degradation, atrazine degradation, mammalian cell entry protein complex, etc. are involved in M. tuberculosis biofilm cells. Especially, ATP-binding cassette (ABC) transporters Rv1217c and Rv1218c were significantly upregulated in biofilm, whereas efflux pump inhibitors (EPIs) piperine and 1-(1-naphthylmethyl)-piperazine (NMP) inhibited biofilm formation and the expression of the Rv1217c and Rv1218c genes in a concentration-dependent manner, respectively, indicating Rv1217c and Rv1218c are potential target genes of M. tuberculosis biofilm. This study is the first RNA-Seq-based transcriptome profiling of M. tuberculosis biofilms and provides insights into a potential strategy for M. tuberculosis biofilm inhibition. KEY POINTS: • Characterize M. tuberculosis transcriptomes in biofilm cells by RNA-seq. • Inhibit the expression of Rv1217c and Rv1218c repressed biofilm formation.

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.

Polynucleotide phosphorylase and RNA helicase CshA cooperate in Bacillus subtilis mRNA decay

Polynucleotide phosphorylase (PNPase), a 3' exoribonuclease that degrades RNA in the 3'-to-5' direction, is the major mRNA decay activity in Bacillus subtilis. PNPase is known to be inhibited in vitro by strong RNA secondary structure, and rapid mRNA turnover in vivo is thought to require an RNA helicase activity working in conjunction with PNPase. The most abundant RNA helicase in B. subtilis is CshA. We found for three small, monocistronic mRNAs that, for some RNA sequences, PNPase processivity was unimpeded even without CshA, whereas others required CshA for efficient degradation. A novel colour screen for decay of mRNA in B. subtilis was created, using mRNA encoded by the slrA gene, which is degraded from its 3' end by PNPase. A significant correlation between the predicted strength of a stem-loop structure, located in the body of the message, and PNPase processivity was observed. Northern blot analysis confirmed that PNPase processivity was greatly hindered by the internal RNA structure, and even more so in the absence of CshA. Three other B. subtilis RNA helicases did not appear to be involved in mRNA decay during vegetative growth. The results confirm the hypothesis that efficient 3' exonucleolytic decay of B. subtilis RNA depends on the combined activity of PNPase and CshA.

Physiological RNA dynamics in RNA-Seq analysis

Physiological RNA dynamics cause problems in transcriptome analysis. Physiological RNA accumulation affects the analysis of RNA quantification, and physiological RNA degradation affects the analysis of the RNA sequence length, feature site and quantification. In the present article, we review the effects of physiological degradation and accumulation of RNA on analysing RNA sequencing data. Physiological RNA accumulation and degradation probably led to such phenomena as incorrect estimations of transcription quantification, differential expressions, co-expressions, RNA decay rates, alternative splicing, boundaries of transcription, novel genes, new single-nucleotide polymorphisms, small RNAs and gene fusion. Thus, the transcriptomic data obtained up to date warrant further scrutiny. New and improved techniques and bioinformatics software are needed to produce accurate data in transcriptome research.

Bacterial ribonucleases and their roles in RNA metabolism

Ribonucleases (RNases) are mediators in most reactions of RNA metabolism. In recent years, there has been a surge of new information about RNases and the roles they play in cell physiology. In this review, a detailed description of bacterial RNases is presented, focusing primarily on those from Escherichia coli and Bacillus subtilis, the model Gram-negative and Gram-positive organisms, from which most of our current knowledge has been derived. Information from other organisms is also included, where relevant. In an extensive catalog of the known bacterial RNases, their structure, mechanism of action, physiological roles, genetics, and possible regulation are described. The RNase complement of E. coli and B. subtilis is compared, emphasizing the similarities, but especially the differences, between the two. Included are figures showing the three major RNA metabolic pathways in E. coli and B. subtilis and highlighting specific steps in each of the pathways catalyzed by the different RNases. This compilation of the currently available knowledge about bacterial RNases will be a useful tool for workers in the RNA field and for others interested in learning about this area.

RNases and Helicases in Gram-Positive Bacteria

RNases are key enzymes involved in RNA maturation and degradation. Although they play a crucial role in all domains of life, bacteria, archaea, and eukaryotes have evolved with their own sets of RNases and proteins modulating their activities. In bacteria, these enzymes allow modulation of gene expression to adapt to rapidly changing environments. Today, >20 RNases have been identified in both Escherichia coli and Bacillus subtilis, the paradigms of the Gram-negative and Gram-positive bacteria, respectively. However, only a handful of these enzymes are common to these two organisms and some of them are essential to only one. Moreover, although sets of RNases can be very similar in closely related bacteria such as the Firmicutes Staphylococcus aureus and B. subtilis, the relative importance of individual enzymes in posttranscriptional regulation in these organisms varies. In this review, we detail the role of the main RNases involved in RNA maturation and degradation in Gram-positive bacteria, with an emphasis on the roles of RNase J1, RNase III, and RNase Y. We also discuss how other proteins such as helicases can modulate the RNA-degradation activities of these enzymes.

PNPase is involved in the coordination of mRNA degradation and expression in stationary phase cells of Escherichia coli

Background

Exoribonucleases are crucial for RNA degradation in Escherichia coli but the roles of RNase R and PNPase and their potential overlap in stationary phase are not well characterized. Here, we used a genome-wide approach to determine how RNase R and PNPase affect the mRNA half-lives in the stationary phase. The genome-wide mRNA half-lives were determined by a dynamic analysis of transcriptomes after transcription arrest. We have combined the analysis of mRNA half-lives with the steady-state concentrations (transcriptome) to provide an integrated overview of the in vivo activity of these exoribonucleases at the genome-scale.

Results

The values of mRNA half-lives demonstrated that the mRNAs are very stable in the stationary phase and that the deletion of RNase R or PNPase caused only a limited mRNA stabilization. Intriguingly the absence of PNPase provoked also the destabilization of many mRNAs. These changes in mRNA half-lives in the PNPase deletion strain were associated with a massive reorganization of mRNA levels and also variation in several ncRNA concentrations. Finally, the in vivo activity of the degradation machinery was found frequently saturated by mRNAs in the PNPase mutant unlike in the RNase R mutant, suggesting that the degradation activity is limited by the deletion of PNPase but not by the deletion of RNase R.

Conclusions

This work had identified PNPase as a central player associated with mRNA degradation in stationary phase.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5259-8) contains supplementary material, which is available to authorized users.

Major 3'-5' Exoribonucleases in the Metabolism of Coding and Non-coding RNA

3'-5' exoribonucleases are key enzymes in the degradation of superfluous or aberrant RNAs and in the maturation of precursor RNAs into their functional forms. The major bacterial 3'-5' exoribonucleases responsible for both these activities are PNPase, RNase II and RNase R. These enzymes are of ancient nature with widespread distribution. In eukaryotes, PNPase and RNase II/RNase R enzymes can be found in the cytosol and in mitochondria and chloroplasts; RNase II/RNase R-like enzymes are also found in the nucleus. Humans express one PNPase (PNPT1) and three RNase II/RNase R family members (Dis3, Dis3L and Dis3L2). These enzymes take part in a multitude of RNA surveillance mechanisms that are critical for translation accuracy. Although active against a wide range of both coding and non-coding RNAs, the different 3'-5' exoribonucleases exhibit distinct substrate affinities. The latest studies on these RNA degradative enzymes have contributed to the identification of additional homologue proteins, the uncovering of novel RNA degradation pathways, and to a better comprehension of several disease-related processes and response to stress, amongst many other exciting findings. Here, we provide a comprehensive and up-to-date overview on the function, structure, regulation and substrate preference of the key 3'-5' exoribonucleases involved in RNA metabolism.

Maturation of polycistronic mRNAs by the endoribonuclease RNase Y and its associated Y-complex in Bacillus subtilis

Endonucleolytic cleavage within polycistronic mRNAs can lead to differential stability, and thus discordant abundance, among cotranscribed genes. RNase Y, the major endonuclease for mRNA decay in Bacillus subtilis, was originally identified for its cleavage activity toward the cggR-gapA operon, an event that differentiates the synthesis of a glycolytic enzyme from its transcriptional regulator. A three-protein Y-complex (YlbF, YmcA, and YaaT) was recently identified as also being required for this cleavage in vivo, raising the possibility that it is an accessory factor acting to regulate RNase Y. However, whether the Y-complex is broadly required for RNase Y activity is unknown. Here, we used end-enrichment RNA sequencing (Rend-seq) to globally identify operon mRNAs that undergo maturation posttranscriptionally by RNase Y and the Y-complex. We found that the Y-complex is required for the majority of RNase Y-mediated mRNA maturation events and also affects riboswitch abundance in B. subtilis In contrast, noncoding RNA maturation by RNase Y often does not require the Y-complex. Furthermore, deletion of RNase Y has more pleiotropic effects on the transcriptome and cell growth than deletions of the Y-complex. We propose that the Y-complex is a specificity factor for RNase Y, with evidence that its role is conserved in Staphylococcus aureus.

Polynucleotide phosphorylase is involved in the control of lipopeptide fengycin production in Bacillus subtilis

Bacillus subtilis is a wealth source of lipopeptide molecules such as iturins, surfactins and fengycins or plipastatins endowed with a range of biological activities. These molecules, designated secondary metabolites, are synthesized via non-ribosomal peptides synthesis (NRPS) machinery and are most often subjected to a complex regulation with involvement of several regulatory factors. To gain novel insights on mechanism regulating fengycin production, we investigated the effect of the fascinating polynucleotide phosphorylase (PNPase), as well as the effect of lipopeptide surfactin. Compared to the wild type, the production of fengycin in the mutant strains B. subtilis BBG235 and BBG236 altered for PNPase has not only decreased to about 70 and 40%, respectively, but also hampered its antifungal activity towards the plant pathogen Botrytis cinerea. On the other hand, mutant strains BBG231 (srfAA^-) and BBG232 (srfAC^-) displayed different levels of fengycin production. BBG231 had registered an important decrease in fengycin production, comparable to that observed for BBG235 or BBG236. This study permitted to establish that the products of pnpA gene (PNPase), and srfAA^- (surfactin synthetase) are involved in fengycin production.

Bacterial RNA Biology on a Genome Scale

Bacteria are an exceedingly diverse group of organisms whose molecular exploration is experiencing a renaissance. While the classical view of bacterial gene expression was relatively simple, the emerging view is more complex, encompassing extensive post-transcriptional control involving riboswitches, RNA thermometers, and regulatory small RNAs (sRNAs) associated with the RNA-binding proteins CsrA, Hfq, and ProQ, as well as CRISPR/Cas systems that are programmed by RNAs. Moreover, increasing interest in members of the human microbiota and environmental microbial communities has highlighted the importance of understudied bacterial species with largely unknown transcriptome structures and RNA-based control mechanisms. Collectively, this creates a need for global RNA biology approaches that can rapidly and comprehensively analyze the RNA composition of a bacterium of interest. We review such approaches with a focus on RNA-seq as a versatile tool to investigate the different layers of gene expression in which RNA is made, processed, regulated, modified, translated, and turned over.

Genome-wide analysis shows that RNase G plays a global role in the stability of mRNAs in Stenotrophomonas maltophilia

Gene expression is determined by critical processes such as RNA synthesis and degradation. Ribonucleases participate in the coordinated and differential decay of messenger RNAs. We describe a suitable method of normalization and calculation of mRNAs half-life values quantified by RNA-Seq. We determined the mRNA half-lives of more than 2000 genes in Stenotrophomonas maltophilia D457 and in an isogenic RNase G deficient mutant. Median half-lives were 2,74 and 3 min in the wild-type and the rng-deficient strain, respectively. The absence of RNase G resulted in an overall enhancement of mRNA half-life times, showing that many RNAs are targets of RNase G in S. maltophilia. Around 40 genes are likely to be regulated directly by RNase G since their half-lives were more than two-fold higher in the rng-deficient mutant. Gene length, GC content or expression levels did not correlate with mRNAs lifetimes, although groups of genes with different functions showed different RNA half-lives. Further, we predicted 1542 gene pairs to be part of the same operons in S. maltophilia. In contrast to what was described for other bacteria, our data indicate that RNase G has a global role in mRNA stability and consequently in the regulation of S. maltophilia gene expression.

Global survey of mRNA levels and decay rates of Chlamydia trachomatis trachoma and lymphogranuloma venereum biovars

Interpreting the intricate bacterial transcriptomics implies understanding the dynamic relationship established between de novo transcription and the degradation of transcripts. Here, we performed a comparative overview of gene expression levels and mRNA decay rates for different-biovar (trachoma and lymphogranuloma venereum) strains of the obligate intracellular bacterium Chlamydia trachomatis. By using RNA-sequencing to measure gene expression levels at mid developmental stage and mRNA decay rates upon rifampicin-based transcription blockage, we observed that: i) 60-70% of the top-50 expressed genes encode proteins with unknown function and proteins involved in "Translation, ribosomal structure and biogenesis" for all strains; ii) the expression ranking by genes' functional categories was in general concordant among different-biovar strains; iii) the median of the half-life time (t_1/2) values of transcripts were 15-17 min, indicating that the degree of transcripts' stability seems to correlate with the bacterial intracellular life-style, as these values are considerably higher than the ones observed in other studies for facultative intracellular and free-living bacteria; iv) transcript decay rates were highly heterogeneous within each C. trachomatis strain and did not correlate with steady-state expression levels; v) only at very few instances (essentially at gene functional category level) was possible to unveil dissimilarities potentially underlying phenotypic differences between biovars. In summary, the unveiled transcriptomic scenario, marked by a general lack of correlation between transcript production and degradation and a huge inter-transcript heterogeneity in decay rates, likely reflects the challenges underlying the unique biphasic developmental cycle of C. trachomatis and its intricate interactions with the human host, which probably exacerbate the complexity of the bacterial transcription regulation.

Identification of endoribonuclease specific cleavage positions reveals novel targets of RNase III in Streptococcus pyogenes

A better understanding of transcriptional and post-transcriptional regulation of gene expression in bacteria relies on studying their transcriptome. RNA sequencing methods are used not only to assess RNA abundance but also the exact boundaries of primary and processed transcripts. Here, we developed a method, called identification of specific cleavage position (ISCP), which enables the identification of direct endoribonuclease targets in vivo by comparing the 5΄ and 3΄ ends of processed transcripts between wild type and RNase deficient strains. To demonstrate the ISCP method, we used as a model the double-stranded specific RNase III in the human pathogen Streptococcus pyogenes. We mapped 92 specific cleavage positions (SCPs) among which, 48 were previously described and 44 are new, with the characteristic 2 nucleotides 3΄ overhang of RNase III. Most SCPs were located in untranslated regions of RNAs. We screened for RNase III targets using transcriptomic differential expression analysis (DEA) and compared those with the RNase III targets identified using the ISCP method. Our study shows that in S. pyogenes, under standard growth conditions, RNase III has a limited impact both on antisense transcripts and on global gene expression with the expression of most of the affected genes being downregulated in an RNase III deletion mutant.

Mini review: Recombinant production of tailored bio-pharmaceuticals in different Bacillus strains and future perspectives

Bio-pharmaceuticals like antibodies, hormones and growth factors represent about one-fifth of commercial pharmaceuticals. Host candidates of growing interest for recombinant production of these proteins are strains of the genus Bacillus, long being established for biotechnological production of homologous and heterologous proteins. Bacillus strains benefit from development of efficient expression systems in the last decades and emerge as major industrial workhorses for recombinant proteins due to easy cultivation, non-pathogenicity and their ability to secrete recombinant proteins directly into extracellular medium allowing cost-effective downstream processing. Their broad product portfolio of pharmaceutically relevant recombinant proteins described in research include antibody fragments, growth factors, interferons and interleukins, insulin, penicillin G acylase, streptavidin and different kinases produced in various cultivation systems like microtiter plates, shake flasks and bioreactor systems in batch, fed-batch and continuous mode. To further improve production and secretion performance of Bacillus, bottlenecks and limiting factors concerning proteases, chaperones, secretion machinery or feedback mechanisms can be identified on different cell levels from genomics and transcriptomics via proteomics to metabolomics and fluxomics. For systematical identification of recurring patterns characteristic of given regulatory systems and key genetic targets, systems biology and omics-technology provide suitable and promising approaches, pushing Bacillus further towards industrial application for recombinant pharmaceutical protein production.

Role of RNase Y in Clostridium perfringens mRNA Decay and Processing

RNase Y is a major endoribonuclease that plays a crucial role in mRNA degradation and processing. We study the role of RNase Y in the Gram-positive anaerobic pathogen Clostridium perfringens, which until now has not been well understood. Our study implies an important role for RNase Y-mediated RNA degradation and processing in virulence gene expression and the physiological development of the organism. We began by constructing an RNase Y conditional knockdown strain in order to observe the importance of RNase Y on growth and virulence. Our resulting transcriptome analysis shows that RNase Y affects the expression of many genes, including toxin-producing genes. We provide data to show that RNase Y depletion repressed several toxin genes in C. perfringens and involved the virR-virS two-component system. We also observe evidence that RNase Y is indispensable for processing and stabilizing the transcripts of colA (encoding a major toxin collagenase) and pilA2 (encoding a major pilin component of the type IV pili). Posttranscriptional regulation of colA is known to be mediated by cleavage in the 5' untranslated region (5'UTR), and we observe that RNase Y depletion diminishes colA 5'UTR processing. We show that RNase Y is also involved in the posttranscriptional stabilization of pilA2 mRNA, which is thought to be important for host cell adherence and biofilm formation.

IMPORTANCE

RNases have important roles in RNA degradation and turnover in all organisms. C. perfringens is a Gram-positive anaerobic spore-forming bacterial pathogen that produces numerous extracellular enzymes and toxins, and it is linked to digestive disorders and disease. A highly conserved endoribonuclease, RNase Y, affects the expression of hundreds of genes, including toxin genes, and studying these effects is useful for understanding C. perfringens specifically and RNases generally. Moreover, RNase Y is involved in processing specific transcripts, and we observed that this processing in C. perfringens results in the stabilization of mRNAs encoding a toxin and bacterial extracellular apparatus pili. Our study shows that RNase activity is associated with gene expression, helping to determine the growth, proliferation, and virulence of C. perfringens.

A Novel Strategy for Exploitation of Host RNase E Activity by a Marine Cyanophage

Previous studies have shown that infection of Prochlorococcus MED4 by the cyanophage P-SSP7 leads to increased transcript levels of host endoribonuclease (RNase) E. However, it has remained enigmatic whether this is part of a host defense mechanism to degrade phage messenger RNA (mRNA) or whether this single-strand RNA-specific RNase is utilized by the phage. Here we describe a hitherto unknown means through which this cyanophage increases expression of RNase E during phage infection and concomitantly protects its own RNA from degradation. We identified two functionally different RNase E mRNA variants, one of which is significantly induced during phage infection. This transcript lacks the 5' UTR, is considerably more stable than the other transcript, and is likely responsible for increased RNase E protein levels during infection. Furthermore, selective enrichment and in vivo analysis of double-stranded RNA (dsRNA) during infection revealed that phage antisense RNAs (asRNAs) sequester complementary mRNAs to form dsRNAs, such that the phage protein-coding transcriptome is nearly completely covered by asRNAs. In contrast, the host protein-coding transcriptome is only partially covered by asRNAs. These data suggest that P-SSP7 orchestrates degradation of host RNA by increasing RNase E expression while masking its own transcriptome from RNase E degradation in dsRNA complexes. We propose that this combination of strategies contributes significantly to phage progeny production.

A multistress responsive type I toxin-antitoxin system: bsrE/SR5 from the B. subtilis chromosome

bsrE/SR5 is a type I TA system from prophage-like element P6 of the B. subtilis chromosome. The 256 nt bsrE RNA encodes a 30 aa toxin. The antitoxin SR5 is a 163 nt antisense RNA. Both genes overlap at their 3' ends. Overexpression of bsrE causes cell lysis on agar plates, which can be neutralized by sr5 overexpression, whereas deletion of the chromosomal sr5 copy has no effect. SR5 is short-lived with a half-life of ≈7 min, whereas bsrE RNA is stable with a half-life of >80 min. The sr5 promoter is 10-fold stronger than the bsrE promoter. SR5 interacts with the 3' UTR of bsrE RNA, thereby promoting its degradation by recruiting RNase III. RNase J1 is the main RNase responsible for SR5 and bsrE RNA degradation, and PnpA processes an SR5 precursor to the mature RNA. Hfq stabilizes SR5, but is not required for its inhibitory function. While bsrE RNA is affected by temperature shock and alkaline stress, the amount of SR5 is significantly influenced by various stresses, among them pH, anoxia and iron limitation. Only the latter one is dependent on sigB. Both RNAs are extremely unstable upon ethanol stress due to rapid degradation by RNase Y.

Mapping of internal monophosphate 5' ends of Bacillus subtilis messenger RNAs and ribosomal RNAs in wild-type and ribonuclease-mutant strains

The recent findings that the narrow-specificity endoribonuclease RNase III and the 5′ exonuclease RNase J1 are not essential in the Gram-positive model organism, Bacillus subtilis, facilitated a global analysis of internal 5′ ends that are generated or acted upon by these enzymes. An RNA-Seq protocol known as PARE (Parallel Analysis of RNA Ends) was used to capture 5′ monophosphorylated RNA ends in ribonuclease wild-type and mutant strains. Comparison of PARE peaks in strains with RNase III present or absent showed that, in addition to its well-known role in ribosomal (rRNA) processing, many coding sequences and intergenic regions appeared to be direct targets of RNase III. These target sites were, in most cases, not associated with a known antisense RNA. The PARE analysis also revealed an accumulation of 3′-proximal peaks that correlated with the absence of RNase J1, confirming the importance of RNase J1 in degrading RNA fragments that contain the transcription terminator structure. A significant result from the PARE analysis was the discovery of an endonuclease cleavage just 2 nts downstream of the 16S rRNA 3′ end. This latter observation begins to answer, at least for B. subtilis, a long-standing question on the exonucleolytic versus endonucleolytic nature of 16S rRNA maturation.

Expression of multiple Bacillus subtilis genes is controlled by decay of slrA mRNA from Rho-dependent 3' ends

Timely turnover of RNA is an important element in the control of bacterial gene expression, but relatively few specific targets of RNA turnover regulation are known. Deletion of the Bacillus subtilis pnpA gene, encoding the major 3' exonuclease turnover enzyme, polynucleotide phosphorylase (PNPase), was shown previously to cause a motility defect correlated with a reduced level of the 32-gene fla/che flagellar biosynthesis operon transcript.fla/che operon transcript abundance has been shown to be inhibited by an excess of the small regulatory protein, SlrA, and here we find that slrA mRNA accumulated in the pnpA-deletion mutant. Mutation of slrA was epistatic to mutation of pnpA for the motility-related phenotype. Further, Rho-dependent termination was required for PNPase turnover of slrA mRNA. When the slrA gene was provided with a Rho-independent transcription terminator, gene regulation was no longer PNPase-dependent. Thus we show that the slrA transcript is a direct target of PNPase and that regulation of RNA turnover is a major determinant of motility gene expression. The interplay of specific transcription termination and mRNA decay mechanisms suggests selection for fine-tuning of gene expression.

Interaction of Bacillus subtilis Polynucleotide Phosphorylase and RNase Y: STRUCTURAL MAPPING AND EFFECT ON mRNA TURNOVER

Polynucleotide phosphorylase (PNPase), a 3'-to-5' phosphorolytic exoribonuclease, is thought to be the primary enzyme responsible for turnover ofBacillus subtilismRNA. The role of PNPase inB. subtilismRNA decay has been analyzed previously by comparison of mRNA profiles in a wild-type strainversusa strain that is deleted forpnpA, the gene encoding PNPase. Recent studies have provided evidence for a degradosome-like complex inB. subtilisthat is built around the major decay-initiating endonuclease, RNase Y, and there is ample evidence for a strong interaction between PNPase and RNase Y. The role of the PNPase-RNase Y interaction in the exonucleolytic function of PNPase needs to be clarified. We sought to construct aB. subtilisstrain containing a catalytically active PNPase that could not interact with RNase Y. Mapping studies of the PNPase-RNase Y interaction were guided by a homology model ofB. subtilisPNPase based on the known structure of theEscherichia coliPNPase in complex with an RNase E peptide. Mutations inB. subtilisresidues predicted to be involved in RNase Y binding showed a loss of PNPase-RNase Y interaction. Two mRNAs whose decay is dependent on RNase Y and PNPase were examined in strains containing full-length PNPase that was either catalytically active but unable to interact with RNase Y, or catalytically inactive but able to interact with RNase Y. At least for these two mRNAs, disruption of the PNPase-RNase Y interaction did not appear to affect mRNA turnover.

Airpnp: Auto- and Integrated Regulation of Polynucleotide Phosphorylase

The properties and expression of polynucleotide phosphorylase (PNPase), capable of both RNA degradation and polymerization, have been studied for 60 years. In this issue of the Journal of Bacteriology, Park et al. (H. Park, H. Yakhnin, M. Connolly, T. Romeo, and P. Babitzke, J Bacteriol 197:3751-3759, 2015, http://dx.doi.org/10.1128/JB.00721-15) write the latest chapter on the complex regulation of pnp gene expression involving CsrA. I describe how this new piece of the puzzle fits into the global scheme of PNPase autoregulation and how this is influenced by central carbon metabolism at both the posttranscriptional level and that of enzyme activity.

sRNA and mRNA turnover in Gram-positive bacteria

It is widely recognized that RNA degradation plays a critical role in gene regulation when fast adaptation of cell growth is required to respond to stress and changing environmental conditions. Bacterial ribonucleases acting alone or in concert with various trans-acting regulatory factors are important mediators of RNA degradation. Here, we will give an overview of what is known about ribonucleases in several Gram-positive bacteria, their specificities and mechanisms of action. In addition, we will illustrate how sRNAs act in a coordinated manner with ribonucleases to regulate the turnover of particular mRNA targets, and the complex interplay existing between the ribosome, the ribonucleases and RNAs.

Genome-wide study of mRNA degradation and transcript elongation in Escherichia coli

An essential part of gene expression is the coordination of RNA synthesis and degradation, which occurs in the same cellular compartment in bacteria. Here, we report a genome-wide RNA degradation study in Escherichia coli using RNA-seq, and present evidence that the stereotypical exponential RNA decay curve obtained using initiation inhibitor, rifampicin, consists of two phases: residual RNA synthesis, a delay in the interruption of steady state that is dependent on distance relative to the mRNA's 5′ end, and the exponential decay. This gives a more accurate RNA lifetime and RNA polymerase elongation rate simultaneously genome-wide. Transcripts typically have a single RNA decay constant along all positions, which is distinct between different operons, indicating that RNA stability is unlikely determined by local sequences. These measurements allowed us to establish a model for RNA processing involving co-transcriptional degradation, providing quantitative description of the macromolecular coordination in gene expression in bacteria on a system-wide level.

Caught in the act: RNA-Seq provides novel insights into mRNA degradation

RNA degradation is a major mechanism of post-transcriptional control of gene expression. Moreover, the rapid turnover of bacterial mRNAs is central to the fast adaptation of these organisms to changing environmental conditions by the regulation of transcription initiation. In most bacteria, RNA processing and degradation require the concerted action of endo- and exoribonucleases. In Molecular Microbiology, Liu et al. (2014) have analysed RNA processing by polynucleotide phosphorylase, the major 3'-5' exonuclease in Bacillus subtilis. For the first time, they were able to study RNA degradation by this enzyme at single nucleotide resolution. The work provides novel insights into the mechanism by which RNA degradation acts in B. subtilis. Moreover, Liu et al. demonstrate that the post-transcriptional control of central regulators affects the expression of whole regulons.