Using EMOTE to map the exact 5'-ends of processed RNA on a transcriptome-wide scale

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.

Insights into the global effect on Staphylococcus aureus growth arrest by induction of the endoribonuclease MazF toxin

A crucial bacterial strategy to avoid killing by antibiotics is to enter a growth arrested state, yet the molecular mechanisms behind this process remain elusive. The conditional overexpression of mazF, the endoribonuclease toxin of the MazEF toxin–antitoxin system in Staphylococcus aureus, is one approach to induce bacterial growth arrest, but its targets remain largely unknown. We used overexpression of mazF and high-throughput sequence analysis following the exact mapping of non-phosphorylated transcriptome ends (nEMOTE) technique to reveal in vivo toxin cleavage sites on a global scale. We obtained a catalogue of MazF cleavage sites and unearthed an extended MazF cleavage specificity that goes beyond the previously reported one. We correlated transcript cleavage and abundance in a global transcriptomic profiling during mazF overexpression. We observed that MazF affects RNA molecules involved in ribosome biogenesis, cell wall synthesis, cell division and RNA turnover and thus deliver a plausible explanation for how mazF overexpression induces stasis. We hypothesize that autoregulation of MazF occurs by directly modulating the MazEF operon, such as the rsbUVW genes that regulate the sigma factor SigB, including an observed cleavage site on the MazF mRNA that would ultimately play a role in entry and exit from bacterial stasis.

Mapping 5'-Ends and Their Phosphorylation State With EMOTE, TSS-EMOTE, and nEMOTE

Knowing the exact 5'-end of a specific RNA molecule has long been of great interest, since this information can reveal transcription start sites, which subsequently pinpoints promoter regions and potential binding sites for transcriptional regulators. Mapping of the 5'-end can also identify endoribonucleolytic cleavage sites, for the maturation of ribosomal RNA or events which initiate degradation of all or parts of the RNA molecule. If the 5' nucleotide of an RNA is triphosphorylated, then the nucleotide represents the original, unaltered 5'-end of the RNA. Similarly, if an endoribonuclease cleaves the RNA with a mechanism that leaves a monophosphate on the 5'-end downstream cleavage product, then it is useful to search exclusively for monophosphorylated 5'-ends to identify the cleavage sites of the endoribonuclease. This chapter details three variations of the EMOTE (exact mapping of transcriptome ends) protocol, which can identify RNA 5'-ends on a transcriptome-wide scale, and simultaneously provide information about the phosphorylation state of the detected ends. EMOTE, TSS-EMOTE (transcription start specific EMOTE), and nEMOTE (nonphosphorylated EMOTE) reveals mono-, tri-, and nonphosphorylated 5'-ends, respectively, and the protocols can be performed individually or in parallel on the same RNA sample.

Post-transcriptional control of virulence gene expression in Staphylococcus aureus

Opportunistic pathogens have to be ready to change life-style whenever the occasion arises, and therefore need to keep tight control over the expression of their virulence factors. Doubly so for commensal bacteria, such as Staphylococcus aureus, which should avoid harming their hosts when they are in a state of peaceful co-existence. S. aureus carries very few sigma factors to help define the transcriptional programs, but instead uses a plethora of small RNA molecules and RNA-RNA interactions to regulate gene expression post-transcriptionally. The endoribonucleases RNase III and RNase Y contribute to this regulatory diversity, and provide a link to RNA-decay and intra-cellular spatiotemporal control of expression. In this review we describe some of these post-transcriptional mechanisms as well as some of the novel transcriptomic approaches that have been used to find and to study them.

Deep sequencing approaches for the analysis of prokaryotic transcriptional boundaries and dynamics

The identification of the protein-coding regions of a genome is straightforward due to the universality of start and stop codons. However, the boundaries of the transcribed regions, conditional operon structures, non-coding RNAs and the dynamics of transcription, such as pausing of elongation, are non-trivial to identify, even in the comparatively simple genomes of prokaryotes. Traditional methods for the study of these areas, such as tiling arrays, are noisy, labour-intensive and lack the resolution required for densely-packed bacterial genomes. Recently, deep sequencing has become increasingly popular for the study of the transcriptome due to its lower costs, higher accuracy and single nucleotide resolution. These methods have revolutionised our understanding of prokaryotic transcriptional dynamics. Here, we review the deep sequencing and data analysis techniques that are available for the study of transcription in prokaryotes, and discuss the bioinformatic considerations of these analyses.

Cell Cycle Constraints and Environmental Control of Local DNA Hypomethylation in α-Proteobacteria

Heritable DNA methylation imprints are ubiquitous and underlie genetic variability from bacteria to humans. In microbial genomes, DNA methylation has been implicated in gene transcription, DNA replication and repair, nucleoid segregation, transposition and virulence of pathogenic strains. Despite the importance of local (hypo)methylation at specific loci, how and when these patterns are established during the cell cycle remains poorly characterized. Taking advantage of the small genomes and the synchronizability of α-proteobacteria, we discovered that conserved determinants of the cell cycle transcriptional circuitry establish specific hypomethylation patterns in the cell cycle model system Caulobacter crescentus. We used genome-wide methyl-N6-adenine (m6A-) analyses by restriction-enzyme-cleavage sequencing (REC-Seq) and single-molecule real-time (SMRT) sequencing to show that MucR, a transcriptional regulator that represses virulence and cell cycle genes in S-phase but no longer in G1-phase, occludes 5'-GANTC-3' sequence motifs that are methylated by the DNA adenine methyltransferase CcrM. Constitutive expression of CcrM or heterologous methylases in at least two different α-proteobacteria homogenizes m6A patterns even when MucR is present and affects promoter activity. Environmental stress (phosphate limitation) can override and reconfigure local hypomethylation patterns imposed by the cell cycle circuitry that dictate when and where local hypomethylation is instated.

TSS-EMOTE, a refined protocol for a more complete and less biased global mapping of transcription start sites in bacterial pathogens

Background

Bacteria rely on efficient gene regulatory mechanisms to switch between genetic programs when they are facing new environments. Although this regulation can occur at many different levels, one of the key steps is the initiation of transcription. Identification of the first nucleotide transcribed by the RNA polymerase is therefore essential to understand the underlying regulatory processes, since this provides insight on promoter strength and binding sites for transcriptional regulators, and additionally reveals the exact 5’ untranslated region of the transcripts, which often contains elements that regulate translation.

Results

Here we present data from a novel TSS-EMOTE assay (Transcription Start Specific Exact Mapping Of Transcriptome Ends) to precisely map the transcription initiation sites of four entire transcriptomes. TSS-EMOTE is a variation of the dRNA-seq method, which has been combined with the EMOTE protocol, in order to increase detection of longer transcripts and limit biases introduced by PCR amplification of the Illumina sequencing library. Using TSS-EMOTE, 2018 promoters were detected in the opportunistic pathogen Staphylococcus aureus, and detailed consensus sequences could be obtained for the RNA polymerase recognition elements (e.g. sigma factor binding sites). The data also revealed a 94 nt median length of the 5’ untranslated region in S. aureus, giving important insights for future work on translational regulation. Additionally, the transcriptomes of three other opportunistic pathogens, Staphylococcus epidermidis, Acinetobacter baumannii and Enterobacter aerogenes, were examined, and the identified promoter locations were then used to generate a map of the operon structure for each of the four organisms.

Conclusions

Mapping transcription start sites, and subsequent correlation with the genomic sequence, provides a multitude of important information about the regulation of gene expression, both at the transcriptional and translational level, by defining 5’ untranslated regions and sigma-factor binding sites. We have here mapped transcription start sites in four important pathogens using TSS-EMOTE, a method that works with both long and 3’-phosphorylated RNA molecules, and which incorporates Unique Molecular Identifiers (UMIs) to allow unbiased quantification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3211-3) contains supplementary material, which is available to authorized users.

Rifampin Resistance rpoB Alleles or Multicopy Thioredoxin/Thioredoxin Reductase Suppresses the Lethality of Disruption of the Global Stress Regulator spx in Staphylococcus aureus

Staphylococcus aureus is capable of causing a remarkable spectrum of disease, ranging from mild skin eruptions to life-threatening infections. The survival and pathogenic potential of S. aureus depend partly on its ability to sense and respond to changes in its environment. Spx is a thiol/oxidative stress sensor that interacts with the C-terminal domain of the RNA polymerase RpoA subunit, leading to changes in gene expression that help sustain viability under various conditions. Using genetic and deep-sequencing methods, we show that spx is essential in S. aureus and that a previously reported Δspx strain harbored suppressor mutations that allowed it to grow without spx One of these mutations is a single missense mutation in rpoB (a P-to-L change at position 519 encoded by rpoB [rpoB-P519L]) that conferred high-level resistance to rifampin. This mutation alone was found to be sufficient to bypass the requirement for spx The generation of rifampin resistance libraries led to the discovery of an additional rpoB mutation, R484H, which supported strains with the spx disruption. Other rifampin resistance mutations either failed to support the Δspx mutant or were recovered at unexpectedly low frequencies in genetic transduction experiments. The amino acid residues encoded by rpoB-P519L and -R484H map in close spatial proximity and comprise a highly conserved region of RpoB. We also discovered that multicopy expression of either trxA (encoding thioredoxin) or trxB (encoding thioredoxin reductase) supports strains with the deletion of spx Our results reveal intriguing properties, especially of RNA polymerase, that compensate for the loss of an essential gene that is a key mediator of diverse processes in S. aureus, including redox and thiol homeostasis, antibiotic resistance, growth, and metabolism.

IMPORTANCE

The survival and pathogenicity of S. aureus depend on complex genetic programs. An objective for combating this insidious organism entails dissecting genetic regulatory circuits and discovering promising new targets for therapeutic intervention. In this study, we discovered that Spx, an RNA polymerase-interacting stress regulator implicated in many stress responses in S. aureus, including responses to oxidative and cell wall antibiotics, is essential. We describe two mechanisms that suppress the lethality of spx disruption. One mechanism highlights how only certain rifampin resistance-encoding alleles of RpoB confer new properties on RNA polymerase, with important mechanistic implications. We describe additional stress conditions where the loss of spx is deleterious, thereby highlighting Spx as a multifaceted regulator and attractive drug discovery target.

How does sub-cellular localization affect the fate of bacterial mRNA?

Recently a number of seminal studies have revealed that both sequence and spatio-temporal factors govern RNA decay in bacteria, which is crucial for regulation of gene expression. Ribonucleases have been described that not only exhibit sequence preferences, but also are sub-cellularly localised. Furthermore, the RNA itself is distributed in an organised manner and does not diffuse freely or randomly within the bacterial cells. Thus, even within the sub-micrometer distances of the bacterial intra-cellular space, the positions of the enzymes and their substrates are kept in check. Adding to this complexity is the secondary structure and sequence specificity that many, perhaps all, ribonucleases exhibit, including those that are responsible for "general" RNA degradation. In this review, the implications of these novel findings are discussed and specific examples from Staphylococcus aureus are analysed.

Growth control switch by a DNA-damage-inducible toxin-antitoxin system in Caulobacter crescentus

Bacterial toxin-antitoxin systems (TASs) are thought to respond to various stresses, often inducing growth-arrested (persistent) sub-populations of cells whose housekeeping functions are inhibited. Many such TASs induce this effect through the translation-dependent RNA cleavage (RNase) activity of their toxins, which are held in check by their cognate antitoxins in the absence of stress. However, it is not always clear whether specific mRNA targets of orthologous RNase toxins are responsible for their phenotypic effect, which has made it difficult to accurately place the multitude of TASs within cellular and adaptive regulatory networks. Here, we show that the TAS HigBA of Caulobacter crescentus can promote and inhibit bacterial growth dependent on the dosage of HigB, a toxin regulated by the DNA damage (SOS) repressor LexA in addition to its antitoxin HigA, and the target selectivity of HigB's mRNA cleavage activity. HigB reduced the expression of an efflux pump that is toxic to a polarity control mutant, cripples the growth of cells lacking LexA, and targets the cell cycle circuitry. Thus, TASs can have outcome switching activity in bacterial adaptive (stress) and systemic (cell cycle) networks.

Decay-Initiating Endoribonucleolytic Cleavage by RNase Y Is Kept under Tight Control via Sequence Preference and Sub-cellular Localisation

Bacteria depend on efficient RNA turnover, both during homeostasis and when rapidly altering gene expression in response to changes. Nevertheless, remarkably few details are known about the rate-limiting steps in targeting and decay of RNA. The membrane-anchored endoribonuclease RNase Y is a virulence factor in Gram-positive pathogens. We have obtained a global picture of Staphylococcus aureus RNase Y sequence specificity using RNA-seq and the novel transcriptome-wide EMOTE method. Ninety-nine endoribonucleolytic sites produced in vivo were precisely mapped, notably inside six out of seven genes whose half-lives increase the most in an RNase Y deletion mutant, and additionally in three separate transcripts encoding degradation ribonucleases, including RNase Y itself, suggesting a regulatory network. We show that RNase Y is required to initiate the major degradation pathway of about a hundred transcripts that are inaccessible to other ribonucleases, but is prevented from promiscuous activity by membrane confinement and sequence preference for guanosines.