RNA sequencing uncovers antisense RNAs and novel small RNAs in Streptococcus pyogenes

Dynamics of diversified A-to-I editing in Streptococcus pyogenes is governed by changes in mRNA stability

Adenosine-to-inosine (A-to-I) RNA editing plays an important role in the post-transcriptional regulation of eukaryotic cell physiology. However, our understanding of the occurrence, function and regulation of A-to-I editing in bacteria remains limited. Bacterial mRNA editing is catalysed by the deaminase TadA, which was originally described to modify a single tRNA in Escherichia coli. Intriguingly, several bacterial species appear to perform A-to-I editing on more than one tRNA. Here, we provide evidence that in the human pathogen Streptococcus pyogenes, tRNA editing has expanded to an additional tRNA substrate. Using RNA sequencing, we identified more than 27 editing sites in the transcriptome of S. pyogenes SF370 and demonstrate that the adaptation of S. pyogenes TadA to a second tRNA substrate has also diversified the sequence context and recoding scope of mRNA editing. Based on the observation that editing is dynamically regulated in response to several infection-relevant stimuli, such as oxidative stress, we further investigated the underlying determinants of editing dynamics and identified mRNA stability as a key modulator of A-to-I editing. Overall, our findings reveal the presence and diversification of A-to-I editing in S. pyogenes and provide novel insights into the plasticity of the editome and its regulation in bacteria.

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.

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

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

Insights into 6S RNA in lactic acid bacteria (LAB)

BACKGROUND

6S RNA is a regulator of cellular transcription that tunes the metabolism of cells. This small non-coding RNA is found in nearly all bacteria and among the most abundant transcripts. Lactic acid bacteria (LAB) constitute a group of microorganisms with strong biotechnological relevance, often exploited as starter cultures for industrial products through fermentation. Some strains are used as probiotics while others represent potential pathogens. Occasional reports of 6S RNA within this group already indicate striking metabolic implications. A conceivable idea is that LAB with 6S RNA defects may metabolize nutrients faster, as inferred from studies of Echerichia coli. This may accelerate fermentation processes with the potential to reduce production costs. Similarly, elevated levels of secondary metabolites might be produced. Evidence for this possibility comes from preliminary findings regarding the production of surfactin in Bacillus subtilis, which has functions similar to those of bacteriocins. The prerequisite for its potential biotechnological utility is a general characterization of 6S RNA in LAB.

RESULTS

We provide a genomic annotation of 6S RNA throughout the Lactobacillales order. It laid the foundation for a bioinformatic characterization of common 6S RNA features. This covers secondary structures, synteny, phylogeny, and product RNA start sites. The canonical 6S RNA structure is formed by a central bulge flanked by helical arms and a template site for product RNA synthesis. 6S RNA exhibits strong syntenic conservation. It is usually flanked by the replication-associated recombination protein A and the universal stress protein A. A catabolite responsive element was identified in over a third of all 6S RNA genes. It is known to modulate gene expression based on the available carbon sources. The presence of antisense transcripts could not be verified as a general trait of LAB 6S RNAs.

CONCLUSIONS

Despite a large number of species and the heterogeneity of LAB, the stress regulator 6S RNA is well-conserved both from a structural as well as a syntenic perspective. This is the first approach to describe 6S RNAs and short 6S RNA-derived transcripts beyond a single species, spanning a large taxonomic group covering multiple families. It yields universal insights into this regulator and complements the findings derived from other bacterial model organisms.

An Inventory of CiaR-Dependent Small Regulatory RNAs in Streptococci

Bacteria adapt to the different environments encountered by rapid and tightly controlled regulations involving complex networks. A first line of control is transcriptional with regulators such as two-component systems (TCSs) that respond to physical and chemical perturbations. It is followed by posttranscriptional regulations in which small regulatory RNAs (sRNAs) may affect RNA translation. Streptococci are opportunistic pathogens for humans and farm animals. The TCS CiaRH is highly conserved among this genus and crucial in bacterial survival under stressful conditions. In several streptococcal species, some sRNAs belong to the CiaRH regulon and are called csRNAs for cia-dependent sRNAs. In this review, we start by focusing on the Streptococcus species harboring a CiaRH TCS. Then the role of CiaRH in streptococcal pathogenesis is discussed in the context of recent studies. Finally, we give an overview of csRNAs and their functions in Streptococci with a focus on their importance in bacterial adaptation and virulence.

Gene Regulation and Transcriptomics

Borrelia (Borreliella) burgdorferi, along with closely related species, is the etiologic agent of Lyme disease. The spirochete subsists in an enzootic cycle that encompasses acquisition from a vertebrate host to a tick vector and transmission from a tick vector to a vertebrate host. To adapt to its environment and persist in each phase of its enzootic cycle, B. burgdorferi wields three systems to regulate the expression of genes: the RpoN-RpoS alternative sigma factor cascade, the Hk1/Rrp1 two-component system and its product c-di-GMP, and the stringent response mediated by Rel_Bbu and DksA. These regulatory systems respond to enzootic phase-specific signals and are controlled or fine- tuned by transcription factors, including BosR and BadR, as well as small RNAs, including DsrABb and Bb6S RNA. In addition, several other DNA-binding and RNA-binding proteins have been identified, although their functions have not all been defined. Global changes in gene expression revealed by high-throughput transcriptomic studies have elucidated various regulons, albeit technical obstacles have mostly limited this experimental approach to cultivated spirochetes. Regardless, we know that the spirochete, which carries a relatively small genome, regulates the expression of a considerable number of genes required for the transitions between the tick vector and the vertebrate host as well as the adaptation to each.

An RNA-seq based comparative approach reveals the transcriptome-wide interplay between 3'-to-5' exoRNases and RNase Y

RNA degradation is an essential process that allows bacteria to control gene expression and adapt to various environmental conditions. It is usually initiated by endoribonucleases (endoRNases), which produce intermediate fragments that are subsequently degraded by exoribonucleases (exoRNases). However, global studies of the coordinated action of these enzymes are lacking. Here, we compare the targetome of endoRNase Y with the targetomes of 3′-to-5′ exoRNases from Streptococcus pyogenes, namely, PNPase, YhaM, and RNase R. We observe that RNase Y preferentially cleaves after guanosine, generating substrate RNAs for the 3′-to-5′ exoRNases. We demonstrate that RNase Y processing is followed by trimming of the newly generated 3′ ends by PNPase and YhaM. Conversely, the RNA 5′ ends produced by RNase Y are rarely further trimmed. Our strategy enables the identification of processing events that are otherwise undetectable. Importantly, this approach allows investigation of the intricate interplay between endo- and exoRNases on a genome-wide scale.

Bacterial RNA degradation is typically initiated by endoribonucleases and followed by exoribonucleases. Here the authors report the targetome of endoRNase Y in Streptococcus pyogenes, revealing the interplay between RNase Y and 3′-to-5′ exoribonuclease PNPase and YhaM.

New Pathogenesis Mechanisms and Translational Leads Identified by Multidimensional Analysis of Necrotizing Myositis in Primates

A fundamental goal of contemporary biomedical research is to understand the molecular basis of disease pathogenesis and exploit this information to develop targeted and more-effective therapies. Necrotizing myositis caused by the bacterial pathogen Streptococcus pyogenes is a devastating human infection with a high mortality rate and few successful therapeutic options. We used dual transcriptome sequencing (RNA-seq) to analyze the transcriptomes of S. pyogenes and host skeletal muscle recovered contemporaneously from infected nonhuman primates. The in vivo bacterial transcriptome was strikingly remodeled compared to organisms grown in vitro, with significant upregulation of genes contributing to virulence and altered regulation of metabolic genes. The transcriptome of muscle tissue from infected nonhuman primates (NHPs) differed significantly from that of mock-infected animals, due in part to substantial changes in genes contributing to inflammation and host defense processes. We discovered significant positive correlations between group A streptococcus (GAS) virulence factor transcripts and genes involved in the host immune response and inflammation. We also discovered significant correlations between the magnitude of bacterial virulence gene expression in vivo and pathogen fitness, as assessed by previously conducted genome-wide transposon-directed insertion site sequencing (TraDIS). By integrating the bacterial RNA-seq data with the fitness data generated by TraDIS, we discovered five new pathogen genes, namely, S. pyogenes 0281 (Spy0281 [dahA]), ihk-irr, slr, isp, and ciaH, that contribute to necrotizing myositis and confirmed these findings using isogenic deletion-mutant strains. Taken together, our study results provide rich new information about the molecular events occurring in severe invasive infection of primate skeletal muscle that has extensive translational research implications.IMPORTANCE Necrotizing myositis caused by Streptococcus pyogenes has high morbidity and mortality rates and relatively few successful therapeutic options. In addition, there is no licensed human S. pyogenes vaccine. To gain enhanced understanding of the molecular basis of this infection, we employed a multidimensional analysis strategy that included dual RNA-seq and other data derived from experimental infection of nonhuman primates. The data were used to target five streptococcal genes for pathogenesis research, resulting in the unambiguous demonstration that these genes contribute to pathogen-host molecular interactions in necrotizing infections. We exploited fitness data derived from a recently conducted genome-wide transposon mutagenesis study to discover significant correlation between the magnitude of bacterial virulence gene expression in vivo and pathogen fitness. Collectively, our findings have significant implications for translational research, potentially including vaccine efforts.

Large-scale and significant expression from pseudogenes in Sodalis glossinidius - a facultative bacterial endosymbiont

The majority of bacterial genomes have high coding efficiencies, but there are some genomes of intracellular bacteria that have low gene density. The genome of the endosymbiont Sodalis glossinidius contains almost 50 % pseudogenes containing mutations that putatively silence them at the genomic level. We have applied multiple 'omic' strategies, combining Illumina and Pacific Biosciences Single-Molecule Real-Time DNA sequencing and annotation, stranded RNA sequencing and proteome analysis to better understand the transcriptional and translational landscape of Sodalis pseudogenes, and potential mechanisms for their control. Between 53 and 74 % of the Sodalis transcriptome remains active in cell-free culture. The mean sense transcription from coding domain sequences (CDSs) is four times greater than that from pseudogenes. Comparative genomic analysis of six Illumina-sequenced Sodalis isolates from different host Glossina species shows pseudogenes make up ~40 % of the 2729 genes in the core genome, suggesting that they are stable and/or that Sodalis is a recent introduction across the genus Glossina as a facultative symbiont. These data shed further light on the importance of transcriptional and translational control in deciphering host-microbe interactions. The combination of genomics, transcriptomics and proteomics gives a multidimensional perspective for studying prokaryotic genomes with a view to elucidating evolutionary adaptation to novel environmental niches.

The roles of small RNAs: insights from bacterial quorum sensing

Ribonucleic acids (RNAs) mainly played auxiliary roles in regulations of genetic processes while recent explorations into small non-coding RNAs (sRNAs) in bacteria have broadened the scope of RNAs studies in these processes. sRNAs have been demonstrated to be involved in various genetic processes and to regulate a variety of bacterial physiologies. Comparatively, quorum sensing (QS) is a mature bacterial cell signaling system which regulates bacteria physiologies as well. Prokaryotic sRNAs studies in the status quo have revealed an emerging picture of trans-kingdom signaling regulation and increasing investigations have demonstrated the feasibility of inter-kingdom signaling as the consequence of QS. We therefore review such phenomena and their similarities to investigate the potential of prokaryote-sourced interkingdom signaling and regulation.

Conserved and specific features of Streptococcus pyogenes and Streptococcus agalactiae transcriptional landscapes

BACKGROUND

The human pathogen Streptococcus pyogenes, or group A Streptococcus, is responsible for mild infections to life-threatening diseases. To facilitate the characterization of regulatory networks involved in the adaptation of this pathogen to its different environments and their evolution, we have determined the primary transcriptome of a serotype M1 S. pyogenes strain at single-nucleotide resolution and compared it with that of Streptococcus agalactiae, also from the pyogenic group of streptococci.

RESULTS

By using a combination of differential RNA-sequencing and oriented RNA-sequencing we have identified 892 transcription start sites (TSS) and 885 promoters in the S. pyogenes M1 strain S119. 8.6% of S. pyogenes mRNAs were leaderless, among which 81% were also classified as leaderless in S. agalactiae. 26% of S. pyogenes transcript 5' untranslated regions (UTRs) were longer than 60 nt. Conservation of long 5' UTRs with S. agalactiae allowed us to predict new potential regulatory sequences. In addition, based on the mapping of 643 transcript ends in the S. pyogenes strain S119, we constructed an operon map of 401 monocistrons and 349 operons covering 81.5% of the genome. One hundred fifty-six operons and 254 monocistrons retained the same organization, despite multiple genomic reorganizations between S. pyogenes and S. agalactiae. Genomic reorganization was found to more often go along with variable promoter sequences and 5' UTR lengths. Finally, we identified 117 putative regulatory RNAs, among which nine were regulated in response to magnesium concentration.

CONCLUSIONS

Our data provide insights into transcriptome evolution in pyogenic streptococci and will facilitate the analysis of genetic polymorphisms identified by comparative genomics in S. pyogenes.

Prioritizing bona fide bacterial small RNAs with machine learning classifiers

Bacterial small (sRNAs) are involved in the control of several cellular processes. Hundreds of putative sRNAs have been identified in many bacterial species through RNA sequencing. The existence of putative sRNAs is usually validated by Northern blot analysis. However, the large amount of novel putative sRNAs reported in the literature makes it impractical to validate each of them in the wet lab. In this work, we applied five machine learning approaches to construct twenty models to discriminate bona fide sRNAs from random genomic sequences in five bacterial species. Sequences were represented using seven features including free energy of their predicted secondary structure, their distances to the closest predicted promoter site and Rho-independent terminator, and their distance to the closest open reading frames (ORFs). To automatically calculate these features, we developed an sRNA Characterization Pipeline (sRNACharP). All seven features used in the classification task contributed positively to the performance of the predictive models. The best performing model obtained a median precision of 100% at 10% recall and of 64% at 40% recall across all five bacterial species, and it outperformed previous published approaches on two benchmark datasets in terms of precision and recall. Our results indicate that even though there is limited sRNA sequence conservation across different bacterial species, there are intrinsic features in the genomic context of sRNAs that are conserved across taxa. We show that these features are utilized by machine learning approaches to learn a species-independent model to prioritize bona fide bacterial sRNAs.

In vivo 3'-to-5' exoribonuclease targetomes of Streptococcus pyogenes

To cope with harsh environments and cause infection, bacteria need to constantly adjust gene expression. Ribonucleases (RNases) control the abundance of regulatory and protein-coding RNA through degradation and maturation. The current characterization of 3′-to-5′ exoribonucleases (exoRNases), processing RNAs from their 3′ end, is solely based on the description of a limited number of targets processed by these RNases. Here, we characterized bacterial 3′-to-5′ exoRNase targetomes. We show that YhaM, polynucleotide phosphorylase (PNPase), and RNase R have exoribonucleolytic activities in the human pathogen Streptococcus pyogenes. We demonstrate that PNPase is the main 3′-to-5′ exoRNase participating in RNA decay, we show that RNase R has a limited processing activity, and we describe an intriguing RNA processing behavior for YhaM.

mRNA decay plays an essential role in the control of gene expression in bacteria. Exoribonucleases (exoRNases), which trim transcripts starting from the 5′ or 3′ end, are particularly important to fully degrade unwanted transcripts and renew the pool of nucleotides available in the cell. While recent techniques have allowed genome-wide identification of ribonuclease (RNase) targets in bacteria in vivo, none of the 3′-to-5′ exoRNase targetomes (i.e., global processing sites) have been studied so far. Here, we report the targetomes of YhaM, polynucleotide phosphorylase (PNPase), and RNase R of the human pathogen Streptococcus pyogenes. We determined that YhaM is an unspecific enzyme that trims a few nucleotides and targets the majority of transcript ends, generated either by transcription termination or by endonucleolytic activity. The molecular determinants for YhaM-limited processivity are yet to be deciphered. We showed that PNPase clears the cell from mRNA decay fragments produced by endoribonucleases (endoRNases) and is the major 3′-to-5′ exoRNase for RNA turnover in S. pyogenes. In particular, PNPase is responsible for the degradation of regulatory elements from 5′ untranslated regions. However, we observed little RNase R activity in standard culture conditions. Overall, our study sheds light on the very distinct features of S. pyogenes 3′-to-5′ exoRNases.

RNase Y-mediated regulation of the streptococcal pyrogenic exotoxin B

Endoribonuclease Y (RNase Y) is a crucial regulator of virulence in Gram-positive bacteria. In the human pathogen Streptococcus pyogenes, RNase Y is required for the expression of the major secreted virulence factor streptococcal pyrogenic exotoxin B (SpeB), but the mechanism involved in this regulation remains elusive. Here, we demonstrate that the 5′ untranslated region of speB mRNA is processed by several RNases including RNase Y. In particular, we identify two RNase Y cleavage sites located downstream of a guanosine (G) residue. To assess whether this nucleotide is required for RNase Y activity in vivo, we mutated it and demonstrate that the presence of this G residue is essential for the processing of the speB mRNA 5′ UTR by RNase Y. Although RNase Y directly targets and processes speB, we show that RNase Y-mediated regulation of speB expression occurs primarily at the transcriptional level and independently of the processing in the speB mRNA 5′ UTR. To conclude, we demonstrate for the first time that RNase Y processing of an mRNA target requires the presence of a G. We also provide new insights on the speB 5′ UTR and on the role of RNase Y in speB regulation.

Functional annotation and distribution overview of RNA families in 27 Streptococcus agalactiae genomes

Background

Streptococcus agalactiae, also known as Group B Streptococcus (GBS), is a Gram-positive bacterium that colonizes the gastrointestinal and genitourinary tract of humans. This bacterium has also been isolated from various animals, such as fish and cattle. Non-coding RNAs (ncRNAs) can act as regulators of gene expression in bacteria, such as Streptococcus pneumoniae and Streptococcus pyogenes. However, little is known about the genomic distribution of ncRNAs and RNA families in S. agalactiae.

Results

Comparative genome analysis of 27 S. agalactiae strains showed more than 5 thousand genomic regions identified and classified as Core, Exclusive, and Shared genome sequences. We identified 27 to 89 RNA families per genome distributed over these regions, from these, 25 were in Core regions while Shared and Exclusive regions showed variations amongst strains. We propose that the amount and type of ncRNA present in each genome can provide a pattern to contribute in the identification of the clonal types.

Conclusions

The identification of RNA families provides an insight over ncRNAs, sRNAs and ribozymes function, that can be further explored as targets for antibiotic development or studied in gene regulation of cellular processes. RNA families could be considered as markers to determine infection capabilities of different strains. Lastly, pan-genome analysis of GBS including the full range of functional transcripts provides a broader approach in the understanding of this pathogen.

Electronic supplementary material

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

A Glycine Riboswitch in Streptococcus pyogenes Controls Expression of a Sodium:Alanine Symporter Family Protein Gene

Regulatory RNAs play important roles in the control of bacterial gene expression. In this study, we investigated gene expression regulation by a putative glycine riboswitch located in the 5'-untranslated region of a sodium:alanine symporter family (SAF) protein gene in the group A Streptococcus pyogenes serotype M49 strain 591. Glycine-dependent gene expression mediated by riboswitch activity was studied using a luciferase reporter gene system. Maximal reporter gene expression was observed in the absence of glycine and in the presence of low glycine concentrations. Differences in glycine-dependent gene expression were not based on differential promoter activity. Expression of the SAF protein gene and the downstream putative cation efflux protein gene was investigated in wild-type bacteria by RT-qPCR transcript analyses. During growth in the presence of glycine (≥1 mM), expression of the genes were downregulated. Northern blot analyses revealed premature transcription termination in the presence of high glycine concentrations. Growth in the presence of 0.1 mM glycine led to the production of a full-length transcript. Furthermore, stability of the SAF protein gene transcript was drastically reduced in the presence of glycine. We conclude that the putative glycine riboswitch in S. pyogenes serotype M49 strain 591 represses expression of the SAF protein gene and the downstream putative cation efflux protein gene in the presence of high glycine concentrations. Sequence and secondary structure comparisons indicated that the streptococcal riboswitch belongs to the class of tandem aptamer glycine riboswitches.

RNase III-Binding-mRNAs Revealed Novel Complementary Transcripts in Streptomyces

cis-Antisense RNAs (asRNAs) provide very simple and effective gene expression control due to the perfect complementarity between regulated and regulatory transcripts. In Streptomyces, the antibiotic-producing clade, the antisense control system is not yet understood, although it might direct the organism's complex development. Initial studies in Streptomyces have found a number of asRNAs. Apart from this, hundreds of mRNAs have been shown to bind RNase III, the double strand-specific endoribonuclease. In this study, we tested 17 mRNAs that have been previously co-precipitated with RNase III for antisense expression. Our RACE mapping showed that all of these mRNAs possess cognate asRNA. Additional tests for antisense expression uncovered as-adpA, as-rnc, as3983, as-sigB, as-sigH, and as-sigR RNAs. Northern blots detected the expression profiles of 18 novel transcripts. Noteworthy, we also found that only a minority of asRNAs respond to the absence of RNase III enzyme by increasing their cellular levels. Our findings suggest that antisense expression is widespread in Streptomyces, including genes of such important developmental regulators, as AdpA, RNase III, and sigma factors.

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.