Identification of 88 regulatory small RNAs in the TIGR4 strain of the human pathogen Streptococcus pneumoniae

An antisense RNA regulates production of DnaA and affects sporulation in Bacillus subtilis

DnaA is the replication initiator and a transcription factor in virtually all bacteria. Although the synthesis and activity of DnaA are highly regulated, the mechanisms of regulation vary between organisms. We found that production of DnaA in Bacillus subtilis is regulated by an antisense RNA that overlaps with the 5' untranslated region upstream of the dnaA open reading frame. We initially observed this RNA in in vitro transcription experiments and found that its production was inhibited by DnaA. This RNA, now called ArrA for antisense RNA repressor of dnaA, is made in vivo. We identified the arrA promoter and made a mutation that greatly reduced (or eliminated) production of ArrA RNA in vitro and in vivo. In vivo, this arrA promoter mutation caused an increase in the amount of mRNA and protein from dnaA and dnaN, indicating that arrA expression normally inhibits expression of the dnaA-dnaN operon. The arrA mutation also caused a delay in sporulation that was alleviated by loss of sda, a sporulation-inhibitory gene that is directly activated by DnaA. arrA appears to be conserved in some members of the Bacillus genus, indicating that arrA has evolved in at least some endospore-forming bacteria to modulate production of DnaA and enable timely and robust sporulation.

Assessing the conservation and targets of putative sRNAs in Streptococcus pneumoniae

RNA regulators are often found in regulatory networks and mediate growth and virulence in bacteria. Small RNAs (sRNAs) are non-coding RNAs that modulate translation initiation and mRNA degradation by base-pairing. To better understand the role of sRNAs in pathogenicity several studies identified sRNAs in Streptococcus pneumoniae; however, little functional characterization has followed. This study's goals are: 1) survey putative sRNAs in S. pneumoniae; 2) assess the conservation of these sRNAs; and 3) examine their predicted targets. Three previous studies in S. pneumoniae identified 287 putative sRNAs by high-throughput sequencing. This study narrows the candidates to a list of 58 putative sRNAs. BLAST analysis indicates that the 58 sequences are highly conserved across the S. pneumoniae pangenome, and 25 of them are identified sporadically in other Streptococcus species. However, only 2 have corresponding sequences identified across several Streptococcus species. We used four RNA-target prediction programs to predict targets for each of the 58 putative sRNAs. Across all probable predictions, six sRNAs have overlapping targets predicted by multiple programs, four of which target numerous transposase encoding transcripts. sRNAs targeting transposase genes display nearly identical and perfect base-pairing. One sRNA, M63 (Spd_sr37), has several probable targets in the CcpA regulon, a network responsible for global catabolite repression, suggesting a possible biological function in carbon metabolism control. Each M63-target interaction exhibits unique base-pairing increasing confidence in the biological relevance of the result. This study produces a list of S. pneumoniae putative sRNAs whose predicted targets suggest functional significance in transposon and carbon metabolism regulation.

The five homologous CiaR-controlled Ccn sRNAs of Streptococcus pneumoniae modulate Zn-resistance

Zinc is a vital transition metal for all bacteria; however, elevated intracellular free Zn levels can result in mis-metalation of Mn-dependent enzymes. For Mn-centric bacteria such as Streptococcus pneumoniae that primarily use Mn instead of Fe as an enzyme cofactor, Zn is particularly toxic at high concentrations. Here, we report our identification and characterization of the function of the five homologous, CiaRH-regulated Ccn sRNAs in controlling S. pneumoniae virulence and metal homeostasis. We show that deletion of all five ccn genes (ccnA, ccnB, ccnC, ccnD, and ccnE) from S. pneumoniae strains D39 (serotype 2) and TIGR4 (serotype 4) causes Zn hypersensitivity and an attenuation of virulence in a murine invasive pneumonia model. We provide evidence that bioavailable Zn disproportionately increases in S. pneumoniae strains lacking the five ccn genes. Consistent with a response to Zn intoxication or relatively high intracellular free Zn levels, expression of genes encoding the CzcD Zn exporter and the Mn-independent ribonucleotide reductase, NrdD-NrdG, were increased in the ΔccnABCDE mutant relative to its isogenic ccn+ parent strain. The growth inhibition by Zn that occurs as the result of loss of the ccn genes is rescued by supplementation with Mn or Oxyrase, a reagent that removes dissolved oxygen. Lastly, we found that the Zn-dependent growth inhibition of the ΔccnABCDE strain was not altered by deletion of sodA, whereas the ccn+ ΔsodA strain phenocopied the ΔccnABCDE strain. Overall, our results indicate that the Ccn sRNAs have a crucial role in preventing Zn intoxication in S. pneumoniae.

A 3'UTR-derived small RNA represses pneumolysin synthesis and facilitates pneumococcal brain invasion

Pneumolysin (Ply) of Streptococcus pneumoniae (pneumococcus) at relatively high and low levels facilitates pneumococcal invasion into the lung and brain, respectively; however, the regulatory mechanisms of Ply expression are poorly understood. Here, we find that a small RNA plyT, processed from the 3'UTR of the ply operon, is expressed higher in anaerobically- than in statically-cultured pneumococcus D39. Using bioinformatic, biochemical and genetic approaches, we reveal that PlyT inhibits Ply synthesis and hemolytic activities by pairing with an RBS-embedded intergenic region of the ply operon. The RNA-binding protein SPD_1558 facilitates the pairing. Importantly, PlyT inhibition of Ply synthesis is stronger in anaerobic culture and leads to lower Ply abundance. Deletion of plyT decreases the number of pneumococci in the infected mouse brain and reduces the virulence, demonstrating that PlyT-regulated lower Ply in oxygen-void microenvironments, such as the blood, is important for pneumococcus to cross the blood-brain barrier and invade the brain. PlyT-mediated repression of Ply synthesis at anoxic niches is also verified in pneumococcal serotype 4 and 14 strains; moreover, the ply operon with a 3'UTR-embedded plyT, and the pairing sequences of IGR and plyT are highly conserved among pneumococcal strains, implying PlyT-regulated Ply synthesis might be widely employed by pneumococcus.

The five homologous CiaR-controlled Ccn sRNAs of Streptococcus pneumoniae modulate Zn-resistance

Zinc is a vital transition metal for Streptococcus pneumoniae, but is deadly at high concentrations. In S. pneumoniae, elevated intracellular free Zn levels result in mis-metallation of key Mn-dependent metabolic and superoxide detoxifying enzymes resulting in Zn intoxication. Here, we report our identification and characterization of the function of the five homologous, CiaRH-regulated Ccn sRNAs in controlling S. pneumoniae virulence and metal homeostasis. We show that deletion of all five ccn genes (ccnA, ccnB, ccnC, ccnD, and ccnE) from S. pneumoniae strains D39 (serotype 2) and TIGR4 (serotype 4) causes Zn hypersensitivity and an attenuation of virulence in a murine invasive pneumonia model. We provide evidence that bioavailable Zn disproportionately increases in S. pneumoniae strains lacking the five ccn genes. Consistent with a response to Zn intoxication or relatively high intracellular free Zn levels, expression of genes encoding the CzcD Zn exporter and the Mn-independent ribonucleotide reductase, NrdD-NrdG, were increased in the ΔccnABCDE mutant relative to its isogenic ccn+ parent strain. The growth inhibition by Zn that occurs as the result of loss of the ccn genes is rescued by supplementation with Mn or Oxyrase™, a reagent that removes dissolved oxygen. Lastly, we found that the Zn-dependent growth inhibition of the ΔccnABCDE strain was not altered by deletion of sodA, whereas the ccn+ ΔsodA strain phenocopied the ΔccnABCDE strain. Overall, our results indicate that the Ccn sRNAs have a crucial role in preventing Zn intoxication in S. pneumoniae.

Diverse regulatory pathways modulate bet hedging of competence induction in epigenetically-differentiated phase variants of Streptococcus pneumoniae

Despite enabling Streptococcus pneumoniae to acquire antibiotic resistance and evade vaccine-induced immunity, transformation occurs at variable rates across pneumococci. Phase variants of isolate RMV7, distinguished by altered methylation patterns driven by the translocating variable restriction-modification (tvr) locus, differed significantly in their transformation efficiencies and biofilm thicknesses. These differences were replicated when the corresponding tvr alleles were introduced into an RMV7 derivative lacking the locus. RNA-seq identified differential expression of the type 1 pilus, causing the variation in biofilm formation, and inhibition of competence induction in the less transformable variant, RMV7domi. This was partly attributable to RMV7domi's lower expression of ManLMN, which promoted competence induction through importing N-acetylglucosamine. This effect was potentiated by analogues of some proteobacterial competence regulatory machinery. Additionally, one of RMV7domi's phage-related chromosomal island was relatively active, which inhibited transformation by increasing expression of the stress response proteins ClpP and HrcA. However, HrcA increased competence induction in the other variant, with its effects depending on Ca2+ supplementation and heat shock. Hence the heterogeneity in transformation efficiency likely reflects the diverse signalling pathways by which it is affected. This regulatory complexity will modulate population-wide responses to synchronising quorum sensing signals to produce co-ordinated yet stochastic bet hedging behaviour.

Post-vaccine epidemiology of serotype 3 pneumococci identifies transformation inhibition through prophage-driven alteration of a non-coding RNA

BACKGROUND

The respiratory pathogen Streptococcus pneumoniae (the pneumococcus) is a genetically diverse bacterium associated with over 101 immunologically distinct polysaccharide capsules (serotypes). Polysaccharide conjugate vaccines (PCVs) have successfully eliminated multiple targeted serotypes, yet the mucoid serotype 3 has persisted despite its inclusion in PCV13. This capsule type is predominantly associated with a single globally disseminated strain, GPSC12 (clonal complex 180).

METHODS

A genomic epidemiology study combined previous surveillance datasets of serotype 3 pneumococci to analyse the population structure, dynamics, and differences in rates of diversification within GPSC12 during the period of PCV introductions. Transcriptomic analyses, whole genome sequencing, mutagenesis, and electron microscopy were used to characterise the phenotypic impact of loci hypothesised to affect this strain's evolution.

RESULTS

GPSC12 was split into clades by a genomic analysis. Clade I, the most common, rarely underwent transformation, but was typically infected with the prophage ϕOXC141. Prior to the introduction of PCV13, this clade's composition shifted towards a ϕOXC141-negative subpopulation in a systematically sampled UK collection. In the post-PCV13 era, more rapidly recombining non-Clade I isolates, also ϕOXC141-negative, have risen in prevalence. The low in vitro transformation efficiency of a Clade I isolate could not be fully explained by the ~100-fold reduction attributable to the serotype 3 capsule. Accordingly, prophage ϕOXC141 was found to modify csRNA3, a non-coding RNA that inhibits the induction of transformation. This alteration was identified in ~30% of all pneumococci and was particularly common in the unusually clonal serotype 1 GPSC2 strain. RNA-seq and quantitative reverse transcriptase PCR experiments using a genetically tractable pneumococcus demonstrated the altered csRNA3 was more effective at inhibiting production of the competence-stimulating peptide pheromone. This resulted in a reduction in the induction of competence for transformation.

CONCLUSION

This interference with the quorum sensing needed to induce competence reduces the risk of the prophage being deleted by homologous recombination. Hence the selfish prophage-driven alteration of a regulatory RNA limits cell-cell communication and horizontal gene transfer, complicating the interpretation of post-vaccine population dynamics.

Role of PatAB Transporter in Efflux of Levofloxacin in Streptococcus pneumoniae

PatAB is an ABC bacterial transporter that facilitates the export of antibiotics and dyes. The overexpression of patAB genes conferring efflux-mediated fluoroquinolone resistance has been observed in several laboratory strains and clinical isolates of Streptococcus pneumoniae. Using transformation and whole-genome sequencing, we characterized the fluoroquinolone-resistance mechanism of one S. pneumoniae clinical isolate without mutations in the DNA topoisomerase genes. We identified the PatAB fluoroquinolone efflux-pump as the mechanism conferring a low-level resistance to ciprofloxacin (8 µg/mL) and levofloxacin (4 µg/mL). Genetic transformation experiments with different amplimers revealed that the entire patA plus the 5'-terminus of patB are required for levofloxacin-efflux. By contrast, only the upstream region of the patAB operon, plus the region coding the N-terminus of PatA containing the G39D, T43A, V48A and D100N amino acid changes, are sufficient to confer a ciprofloxacin-efflux phenotype, thus suggesting differences between fluoroquinolones in their binding and/or translocation pathways. In addition, we identified a novel single mutation responsible for the constitutive and ciprofloxacin-inducible upregulation of patAB. This mutation is predicted to destabilize the putative rho-independent transcriptional terminator located upstream of patA, increasing transcription of downstream genes. This is the first report demonstrating the role of the PatAB transporter in levofloxacin-efflux in a pneumoccocal clinical isolate.

Developing New Tools to Fight Human Pathogens: A Journey through the Advances in RNA Technologies

A long scientific journey has led to prominent technological advances in the RNA field, and several new types of molecules have been discovered, from non-coding RNAs (ncRNAs) to riboswitches, small interfering RNAs (siRNAs) and CRISPR systems. Such findings, together with the recognition of the advantages of RNA in terms of its functional performance, have attracted the attention of synthetic biologists to create potent RNA-based tools for biotechnological and medical applications. In this review, we have gathered the knowledge on the connection between RNA metabolism and pathogenesis in Gram-positive and Gram-negative bacteria. We further discuss how RNA techniques have contributed to the building of this knowledge and the development of new tools in synthetic biology for the diagnosis and treatment of diseases caused by pathogenic microorganisms. Infectious diseases are still a world-leading cause of death and morbidity, and RNA-based therapeutics have arisen as an alternative way to achieve success. There are still obstacles to overcome in its application, but much progress has been made in a fast and effective manner, paving the way for the solid establishment of RNA-based therapies in the future.

Comparative Genomic Analysis of Streptococcus pneumoniae Strains: Penicillin Non-susceptible Multi-drug-Resistant Serotype 19A Isolates

Streptococcus pneumoniae can cause several diseases including otitis media, sinusitis, pneumonia, sepsis and meningitis. The introduction of pneumococcal vaccines has changed the molecular epidemiological and antibiotic resistance profiles of related diseases. Analysis of molecular patterns and genome sequences of clinical strains may facilitate the identification of novel drug resistance mechanism. Three multidrug resistance 19A isolates were verified, serotyped and the complete genomes were sequenced combining the Pacific Biosciences and the Illumina Miseq platform. Genomic annotation revealed that similar central networks were found in the clinical isolates, and Mauve alignments indicated high similarity between different strains. The pan-genome analysis showed the shared and unique cluster in the strains. Mobile elements were predicted in the isolates including prophages and CRISPER systems, which may participate in the virulence and antibiotic resistance of the strains. The presence of 31 virulence factor genes was predicted from other pathogens for PRSP 19339 and 19343, while 30 for PRSP 19087. Meanwhile, 33 genes antibiotic resistance genes were predicted including antibiotic resistance genes, antibiotic-target genes and antibiotic biosynthesis genes. Further analysis of the antibiotic resistance genes revealed new mutations in the isolates. By comparative genomic analysis, we contributed to the understanding of resistance mechanism of the clinical isolates with other serotype strains, which could facilitate the concrete drug resistance mechanism study.

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.

A Small Non-Coding RNA Modulates Expression of Pilus-1 Type in Streptococcus pneumoniae

Streptococcus pneumoniae is a major cause of morbidity and mortality worldwide, and about 30% of the pneumococcal clinical isolates show type I pili-like structures. These long proteinaceous polymers extending from the bacterial surface are encoded by pilus islet 1 and play major roles in adhesion and host colonization. Pili expression is bistable and is controlled by the transcriptional activator RlrA. In this work, we demonstrate that the previously identified small noncoding RNA srn135 also participates in pilus regulation. Our findings show that srn135 is generated upon processing of the 5′-UTR region of rrgA messenger and its deletion prevents the synthesis of RrgA, the main pili adhesin. Moreover, overexpression of srn135 increases the expression of all pili genes and rises the percentage of piliated bacteria within a clonal population. This regulation is mediated by the stabilization of rlrA mRNA since higher levels of srn135 increase its half-life to 165%. Our findings suggest that srn135 has a dual role in pilus expression acting both in cis- (on the RrgA levels) and in trans- (modulating the levels of RlrA) and contributes to the delicate balance between pili expressing and non-expressing bacteria.

Involvement of non-coding RNAs during infection of rice by Acidovorax oryzae

The expression of non-coding RNAs (ncRNAs) has been observed in a variety of bacteria. However, the function of ncRNAs and their regulatory targets are largely unknown, and few ncRNAs are found to be associated with bacterial virulence. The bacterial brown stripe pathogen Acidovorax oryzae (Ao) RS-1 shows a high level of condition-dependent differential expression of ncRNA, which we identified in a genome wide screen. We experimentally validated 66 differentially expressed ncRNAs using an integrative analysis of conservative genome sequences and transcriptomic data during in vivo interaction of the bacterial pathogen with the rice plant. To test the relevance of the differentially expressed ncRNAs, we chose four with different positions within the genome, and with different secondary structures and promoter activities. The results show that the overexpression of the four ncRNAs caused a significant change in virulence-related phenotypes, resistance to various environmental stresses, expression of secretion systems and effector proteins, while changing the expression of ncRNA putative target genes. We conclude that these ncRNAs are examples for the inherent regulatory roles for many of the observed ncRNAs in response to changing conditions such as host interaction or environmental adaption.

Recent Research Advances in Small Regulatory RNAs in Streptococcus

Small non-coding RNAs (sRNAs) are a class of regulatory RNAs 20-500 nucleotides in length, which have recently been discovered in prokaryotic organisms. sRNAs are key regulators in many biological processes, such as sensing various environmental changes and regulating intracellular gene expression through binding target mRNAs or proteins. Bacterial sRNAs have recently been rapidly mined, thus providing new insights into the regulatory network of biological functions in prokaryotes. Although most bacterial sRNAs have been discovered and studied in Escherichia coli and other Gram-negative bacteria, sRNAs have increasingly been predicted and verified in Gram-positive bacteria in the past decade. The genus Streptococcus includes many commensal and pathogenic Gram-positive bacteria. However, current understanding of sRNA-mediated regulation in Streptococcus is limited. Most known sRNAs in Streptococcus are associated with the regulation of virulence. In this review, we summarize recent advances in understanding of the functions and mechanisms of sRNAs in Streptococcus, and we discuss the RNA chaperone protein and synthetic sRNA-mediated gene regulation, with the aim of providing a reference for the study of microbial sRNAs.

Interplay between Regulatory RNAs and Signal Transduction Systems during Bacterial Infection

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

S1 Domain RNA-Binding Protein CvfD Is a New Posttranscriptional Regulator That Mediates Cold Sensitivity, Phosphate Transport, and Virulence in Streptococcus pneumoniae D39

Posttranscriptional gene regulation often involves RNA-binding proteins that modulate mRNA translation and/or stability either directly through protein-RNA interactions or indirectly by facilitating the annealing of small regulatory RNAs (sRNAs). The human pathogen Streptococcus pneumoniae D39 (pneumococcus) does not encode homologs to RNA-binding proteins known to be involved in promoting sRNA stability and function, such as Hfq or ProQ, even though it contains genes for at least 112 sRNAs. However, the pneumococcal genome contains genes for other RNA-binding proteins, including at least six S1 domain proteins: ribosomal protein S1 (rpsA), polynucleotide phosphorylase (pnpA), RNase R (rnr), and three proteins with unknown functions. Here, we characterize the function of one of these conserved, yet uncharacterized, S1 domain proteins, SPD_1366, which we have renamed CvfD (conserved virulence factor D), since loss of the protein results in attenuation of virulence in a murine pneumonia model. We report that deletion of cvfD impacts the expression of 144 transcripts, including the pst1 operon, encoding phosphate transport system 1 in S. pneumoniae We further show that CvfD posttranscriptionally regulates the PhoU2 master regulator of the pneumococcal dual-phosphate transport system by binding phoU2 mRNA and impacting PhoU2 translation. CvfD not only controls expression of phosphate transporter genes but also functions as a pleiotropic regulator that impacts cold sensitivity and the expression of sRNAs and genes involved in diverse cellular functions, including manganese uptake and zinc efflux. Together, our data show that CvfD exerts a broad impact on pneumococcal physiology and virulence, partly by posttranscriptional gene regulation.IMPORTANCE Recent advances have led to the identification of numerous sRNAs in the major human respiratory pathogen S. pneumoniae However, little is known about the functions of most sRNAs or RNA-binding proteins involved in RNA biology in pneumococcus. In this paper, we characterize the phenotypes and one target of the S1 domain RNA-binding protein CvfD, a homolog of general stress protein 13 identified, but not extensively characterized, in other Firmicutes species. Pneumococcal CvfD is a broadly pleiotropic regulator, whose absence results in misregulation of divalent cation homeostasis, reduced translation of the PhoU2 master regulator of phosphate uptake, altered metabolism and sRNA amounts, cold sensitivity, and attenuation of virulence. These findings underscore the critical roles of RNA biology in pneumococcal physiology and virulence.

Grad-seq in a Gram-positive bacterium reveals exonucleolytic sRNA activation in competence control

RNA–protein interactions are the crucial basis for many steps of bacterial gene expression, including post‐transcriptional control by small regulatory RNAs (sRNAs). In stark contrast to recent progress in the analysis of Gram‐negative bacteria, knowledge about RNA–protein complexes in Gram‐positive species remains scarce. Here, we used the Grad‐seq approach to draft a comprehensive landscape of such complexes in Streptococcus pneumoniae, in total determining the sedimentation profiles of ~ 88% of the transcripts and ~ 62% of the proteins of this important human pathogen. Analysis of in‐gradient distributions and subsequent tag‐based protein capture identified interactions of the exoribonuclease Cbf1/YhaM with sRNAs that control bacterial competence for DNA uptake. Unexpectedly, the nucleolytic activity of Cbf1 stabilizes these sRNAs, thereby promoting their function as repressors of competence. Overall, these results provide the first RNA/protein complexome resource of a Gram‐positive species and illustrate how this can be utilized to identify new molecular factors with functions in RNA‐based regulation of virulence‐relevant pathways.

Targeting Hsp20 Using the Novel Small Non-coding RNA DnrH Regulates Heat Tolerance in Deinococcus radiodurans

Small non-coding RNAs (ncRNAs) are a class of regulatory molecules, which remain understudied in bacteria. In the extremophilic bacterium Deinococcus radiodurans, although hundreds of ncRNAs have been identified, few have been characterized in detail. In this study, we report the identification and characterization of a novel heat-inducible ncRNA named DnrH. Heat tolerance analysis showed that deleting DnrH significantly inhibited viability in response to high temperature conditions. Comparative phenotypic and qRT-PCR analyses of a DnrH mutant (ΔDnrH) and wild-type (WT) D. radiodurans suggested that DnrH is potentially involved in regulating the expression of the heat shock-related gene Hsp20. Microscale thermophoresis and genetic complementation showed that a 28-nucleotide (nt) sequence in the stem-loop structure of DnrH (143–170 nt) pairs with its counterpart in the coding region of Hsp20 mRNA (91–117 nt) via a 22 nt region. In vivo, mutation of the 22-nt region in the D. radiodurans genome led to a reduction in heat tolerance similar to that observed in the DnrH-mutant. Our results show that DnrH positively influences heat tolerance by increasing the transcription of Hsp20 mRNA, demonstrating, for the first time, a ncRNA that directly controls the expression of a heat stress-resistance gene. This work provides new insight into the heat stress response mechanism of D. radiodurans as well as other extremophiles that express similar Hsp20 proteins.

Dual RNA-seq in Streptococcus pneumoniae Infection Reveals Compartmentalized Neutrophil Responses in Lung and Pleural Space

The factors that regulate the passage of bacteria between different anatomical compartments are unclear. We have used an experimental model of infection with Streptococcus pneumoniae to examine the host and bacterial factors involved in the passage of bacteria from the lung to the pleural space. The transcriptional profile of host and bacterial cells within the pleural space and lung was analyzed using deep sequencing of the entire transcriptome using the technique of dual RNA-seq. We found significant differences in the host and bacterial RNA profiles in infection, which shed light on the key factors that allow passage of this bacterium into the pleural space.

Streptococcus pneumoniae is the dominant cause of community-acquired pneumonia worldwide. Invasion of the pleural space is common and results in increased mortality. We set out to determine the bacterial and host factors that influence invasion of the pleural space. In a murine model of pneumococcal infection, we isolated neutrophil-dominated samples of bronchoalveolar and pleural fluid containing bacteria 48 hours after infection. Using dual RNA sequencing (RNA-seq), we characterized bacterial and host transcripts that were differentially regulated between these compartments and bacteria in broth and resting neutrophils, respectively. Pleural and lung samples showed upregulation of genes involved in the positive regulation of neutrophil extravasation but downregulation of genes mediating bacterial killing. Compared to the lung samples, cells within the pleural space showed marked upregulation of many genes induced by type I interferons, which are cytokines implicated in preventing bacterial transmigration across epithelial barriers. Differences in the bacterial transcripts between the infected samples and bacteria grown in broth showed the upregulation of genes in the bacteriocin locus, the pneumococcal surface adhesin PsaA, and the glycopeptide resistance gene vanZ; the gene encoding the ClpP protease was downregulated in infection. One hundred sixty-nine intergenic putative small bacterial RNAs were also identified, of which 43 (25.4%) small RNAs had been previously described. Forty-two of the small RNAs were upregulated in pleura compared to broth, including many previously identified as being important in virulence. Our results have identified key host and bacterial responses to invasion of the pleural space that can be potentially exploited to develop alternative antimicrobial strategies for the prevention and treatment of pneumococcal pleural disease.

IMPORTANCE The factors that regulate the passage of bacteria between different anatomical compartments are unclear. We have used an experimental model of infection with Streptococcus pneumoniae to examine the host and bacterial factors involved in the passage of bacteria from the lung to the pleural space. The transcriptional profile of host and bacterial cells within the pleural space and lung was analyzed using deep sequencing of the entire transcriptome using the technique of dual RNA-seq. We found significant differences in the host and bacterial RNA profiles in infection, which shed light on the key factors that allow passage of this bacterium into the pleural space.

A Mn-sensing riboswitch activates expression of a Mn2+/Ca2+ ATPase transporter in Streptococcus

Maintaining manganese (Mn) homeostasis is important for the virulence of numerous bacteria. In the human respiratory pathogen Streptococcus pneumoniae, the Mn-specific importer PsaBCA, exporter MntE, and transcriptional regulator PsaR establish Mn homeostasis. In other bacteria, Mn homeostasis is controlled by yybP-ykoY family riboswitches. Here, we characterize a yybP-ykoY family riboswitch upstream of the mgtA gene encoding a PII-type ATPase in S. pneumoniae, suggested previously to function in Ca2+ efflux. We show that the mgtA riboswitch aptamer domain adopts a canonical yybP-ykoY structure containing a three-way junction that is compacted in the presence of Ca2+ or Mn2+ at a physiological Mg2+ concentration. Although Ca2+ binds to the RNA aptamer with higher affinity than Mn2+, in vitro activation of transcription read-through of mgtA by Mn2+ is much greater than by Ca2+. Consistent with this result, mgtA mRNA and protein levels increase ≈5-fold during cellular Mn stress, but only in genetic backgrounds of S. pneumoniae and Bacillus subtilis that exhibit Mn2+ sensitivity, revealing that this riboswitch functions as a failsafe 'on' signal to prevent Mn2+ toxicity in the presence of high cellular Mn2+. In addition, our results suggest that the S. pneumoniae yybP-ykoY riboswitch functions to regulate Ca2+ efflux under these conditions.

Deep genome annotation of the opportunistic human pathogen Streptococcus pneumoniae D39

A precise understanding of the genomic organization into transcriptional units and their regulation is essential for our comprehension of opportunistic human pathogens and how they cause disease. Using single-molecule real-time (PacBio) sequencing we unambiguously determined the genome sequence of Streptococcus pneumoniae strain D39 and revealed several inversions previously undetected by short-read sequencing. Significantly, a chromosomal inversion results in antigenic variation of PhtD, an important surface-exposed virulence factor. We generated a new genome annotation using automated tools, followed by manual curation, reflecting the current knowledge in the field. By combining sequence-driven terminator prediction, deep paired-end transcriptome sequencing and enrichment of primary transcripts by Cappable-Seq, we mapped 1015 transcriptional start sites and 748 termination sites. We show that the pneumococcal transcriptional landscape is complex and includes many secondary, antisense and internal promoters. Using this new genomic map, we identified several new small RNAs (sRNAs), RNA switches (including sixteen previously misidentified as sRNAs), and antisense RNAs. In total, we annotated 89 new protein-encoding genes, 34 sRNAs and 165 pseudogenes, bringing the S. pneumoniae D39 repertoire to 2146 genetic elements. We report operon structures and observed that 9% of operons are leaderless. The genome data are accessible in an online resource called PneumoBrowse (https://veeninglab.com/pneumobrowse) providing one of the most complete inventories of a bacterial genome to date. PneumoBrowse will accelerate pneumococcal research and the development of new prevention and treatment strategies.

Redefining the Small Regulatory RNA Transcriptome in Streptococcus pneumoniae Serotype 2 Strain D39

Streptococcus pneumoniae (pneumococcus) is a major human respiratory pathogen and a leading cause of bacterial pneumonia worldwide. Small regulatory RNAs (sRNAs), which often act by posttranscriptionally regulating gene expression, have been shown to be crucial for the virulence of S. pneumoniae and other bacterial pathogens. Over 170 putative sRNAs have been identified in the S. pneumoniae TIGR4 strain (serotype 4) through transcriptomic studies, and a subset of these sRNAs has been further implicated in regulating pneumococcal pathogenesis. However, there is little overlap in the sRNAs identified among these studies, which indicates that the approaches used for sRNA identification were not sufficiently sensitive and robust and that there are likely many more undiscovered sRNAs encoded in the S. pneumoniae genome. Here, we sought to comprehensively identify sRNAs in Avery's virulent S. pneumoniae strain D39 using two independent RNA sequencing (RNA-seq)-based approaches. We developed an unbiased method for identifying novel sRNAs from bacterial RNA-seq data and have further tested the specificity of our analysis program toward identifying sRNAs encoded by both strains D39 and TIGR4. Interestingly, the genes for 15% of the putative sRNAs identified in strain TIGR4, including ones previously implicated in virulence, are not present in the strain D39 genome, suggesting that the differences in sRNA repertoires between these two serotypes may contribute to their strain-specific virulence properties. Finally, this study has identified 66 new sRNA candidates in strain D39, 30 of which have been further validated, raising the total number of sRNAs that have been identified in strain D39 to 112.IMPORTANCE Recent work has shown that sRNAs play crucial roles in S. pneumoniae pathogenesis, as inactivation of nearly one-third of the putative sRNA genes identified in one study led to reduced fitness or virulence in a murine model. Yet our understanding of sRNA-mediated gene regulation in S. pneumoniae has been hindered by limited knowledge about these regulatory RNAs, including which sRNAs are synthesized by different S. pneumoniae strains. We sought to address this problem by developing a sensitive sRNA detection technique to identify sRNAs in S. pneumoniae D39. A comparison of our data set reported here to those of other RNA-seq studies for S. pneumoniae strain D39 and TIGR4 has provided new insights into the S. pneumoniae sRNA transcriptome.

The Transcriptional landscape of Streptococcus pneumoniae TIGR4 reveals a complex operon architecture and abundant riboregulation critical for growth and virulence

Efficient and highly organized regulation of transcription is fundamental to an organism’s ability to survive, proliferate, and quickly respond to its environment. Therefore, precise mapping of transcriptional units and understanding their regulation is crucial to determining how pathogenic bacteria cause disease and how they may be inhibited. In this study, we map the transcriptional landscape of the bacterial pathogen Streptococcus pneumoniae TIGR4 by applying a combination of high-throughput RNA-sequencing techniques. We successfully map 1864 high confidence transcription termination sites (TTSs), 790 high confidence transcription start sites (TSSs) (742 primary, and 48 secondary), and 1360 low confidence TSSs (74 secondary and 1286 primary) to yield a total of 2150 TSSs. Furthermore, our study reveals a complex transcriptome wherein environment-respondent alternate transcriptional units are observed within operons stemming from internal TSSs and TTSs. Additionally, we identify many putative cis-regulatory RNA elements and riboswitches within 5’-untranslated regions (5’-UTR). By integrating TSSs and TTSs with independently collected RNA-Seq datasets from a variety of conditions, we establish the response of these regulators to changes in growth conditions and validate several of them. Furthermore, to demonstrate the importance of ribo-regulation by 5’-UTR elements for in vivo virulence, we show that the pyrR regulatory element is essential for survival, successful colonization and infection in mice suggesting that such RNA elements are potential drug targets. Importantly, we show that our approach of combining high-throughput sequencing with in vivo experiments can reconstruct a global understanding of regulation, but also pave the way for discovery of compounds that target (ribo-)regulators to mitigate virulence and antibiotic resistance.

The canonical relationship between a bacterial operon and the mRNA transcript produced from the operon has become significantly more complex as numerous regulatory mechanisms that impact the stability, translational efficiency, and early termination rates for mRNA transcripts have been described. With the rise of antibiotic resistance, these mechanisms offer new potential targets for antibiotic development. In this study we used a combination of high-throughput sequencing technologies to assess genome-wide transcription start and stop sites, as well as determine condition specific global transcription patterns in the human pathogen Streptococcus pneumoniae. We find that the majority of multi-gene operons have alternative start and stop sites enabling condition specific regulation of genes within the same operon. Furthermore, we identified many putative RNA regulators that are widespread in the S. pneumoniae pan-genome. Finally, we show that separately collected RNA-Seq data enables identification of conditional triggers for regulatory RNAs, and experimentally demonstrate that our approach may be used to identify drug-able RNA targets by establishing that pyrR RNA functionality is critical for successful S. pneumoniae mouse colonization and infection. Thus, our study not only uses genome-wide high-throughput approaches to identify putative RNA regulators, but also establishes the importance of such regulators in S. pneumoniae virulence.

Bacterial Small RNAs in Mixed Regulatory Networks

Small regulatory RNAs are now recognized as key regulators of gene expression in bacteria. They accumulate under specific conditions, most often because their synthesis is directly controlled by transcriptional regulators, including but not limited to alternative sigma factors and response regulators of two-component systems. In turn, small RNAs regulate, mostly at the posttranscriptional level, expression of multiple genes, among which are genes encoding transcriptional regulators. Small RNAs are thus embedded in mixed regulatory circuits combining transcriptional and posttranscriptional controls, and whose properties are discussed here.

Genome-wide identification and characterization of small RNAs in Rhodobacter capsulatus and identification of small RNAs affected by loss of the response regulator CtrA

Small non-coding RNAs (sRNAs) are involved in the control of numerous cellular processes through various regulatory mechanisms, and in the past decade many studies have identified sRNAs in a multitude of bacterial species using RNA sequencing (RNA-seq). Here, we present the first genome-wide analysis of sRNA sequencing data in Rhodobacter capsulatus, a purple nonsulfur photosynthetic alphaproteobacterium. Using a recently developed bioinformatics approach, sRNA-Detect, we detected 422 putative sRNAs from R. capsulatus RNA-seq data. Based on their sequence similarity to sRNAs in a sRNA collection, consisting of published putative sRNAs from 23 additional bacterial species, and RNA databases, the sequences of 124 putative sRNAs were conserved in at least one other bacterial species; and, 19 putative sRNAs were assigned a predicted function. We bioinformatically characterized all putative sRNAs and applied machine learning approaches to calculate the probability of a nucleotide sequence to be a bona fide sRNA. The resulting quantitative model was able to correctly classify 95.2% of sequences in a validation set. We found that putative cis-targets for antisense and partially overlapping sRNAs were enriched with protein-coding genes involved in primary metabolic processes, photosynthesis, compound binding, and with genes forming part of macromolecular complexes. We performed differential expression analysis to compare the wild type strain to a mutant lacking the response regulator CtrA, an important regulator of gene expression in R. capsulatus, and identified 18 putative sRNAs with differing levels in the two strains. Finally, we validated the existence and expression patterns of four novel sRNAs by Northern blot analysis.

Comparative Analyses of Selection Operating on Nontranslated Intergenic Regions of Diverse Bacterial Species

Nontranslated intergenic regions (IGRs) compose 10-15% of bacterial genomes, and contain many regulatory elements with key functions. Despite this, there are few systematic studies on the strength and direction of selection operating on IGRs in bacteria using whole-genome sequence data sets. Here we exploit representative whole-genome data sets from six diverse bacterial species: Staphylococcus aureus, Streptococcus pneumoniae, Mycobacterium tuberculosis, Salmonella enterica, Klebsiella pneumoniae, and Escherichia coli We compare patterns of selection operating on IGRs using two independent methods: the proportion of singleton mutations and the d_I/d_S ratio, where d_I is the number of intergenic SNPs per intergenic site. We find that the strength of purifying selection operating over all intergenic sites is consistently intermediate between that operating on synonymous and nonsynonymous sites. Ribosome binding sites and noncoding RNAs tend to be under stronger selective constraint than promoters and Rho-independent terminators. Strikingly, a clear signal of purifying selection remains even when all these major categories of regulatory elements are excluded, and this constraint is highest immediately upstream of genes. While a paucity of variation means that the data for M. tuberculosis are more equivocal than for the other species, we find strong evidence for positive selection within promoters of this species. This points to a key adaptive role for regulatory changes in this important pathogen. Our study underlines the feasibility and utility of gauging the selective forces operating on bacterial IGRs from whole-genome sequence data, and suggests that our current understanding of the functionality of these sequences is far from complete.

Identification and validation of sRNAs in Edwardsiella tarda S08

Bacterial small non-coding RNAs (sRNAs) are known as novel regulators involved in virulence, stress responsibility, and so on. Recently, a lot of new researches have highlighted the critical roles of sRNAs in fine-tune gene regulation in both prokaryotes and eukaryotes. Edwardsiella tarda (E. tarda) is a gram-negative, intracellular pathogen that causes edwardsiellosis in fish. Thus far, no sRNA has been reported in E. tarda. The present study represents the first attempt to identify sRNAs in E. tarda S08. Ten sRNAs were validated by RNA sequencing and quantitative PCR (qPCR). ET_sRNA_1 and ET_sRNA_2 were homolous to tmRNA and GcvB, respectively. However, the other candidate sRNAs have not been reported till now. The cellular abundance of 10 validated sRNA was detected by qPCR at different growth phases to monitor their biosynthesis. Nine candidate sRNAs were expressed in the late-stage of exponential growth and stationary stages of growth (36~60 h). And the expression of the nine sRNAs was growth phase-dependent. But ET_sRNA_10 was almost expressed all the time and reached the highest peak at 48 h. Their targets were predicted by TargetRNA2 and each sRNA target contains some genes that directly or indirectly relate to virulence. These results preliminary showed that sRNAs probably play a regulatory role of virulence in E. tarda.

Absence of tmRNA Has a Protective Effect against Fluoroquinolones in Streptococcus pneumoniae

The transfer messenger RNA (tmRNA), encoded by the ssrA gene, is a small non-coding RNA involved in trans-translation that contributes to the recycling of ribosomes stalled on aberrant mRNAs. In most bacteria, its inactivation has been related to a decreased ability to respond to and recover from a variety of stress conditions. In this report, we investigated the role of tmRNA in stress adaptation in the human pathogen Streptococcus pneumoniae. We constructed a tmRNA deletion mutant and analyzed its response to several lethal stresses. The ΔssrA strain grew slower than the wild type, indicating that, although not essential, tmRNA is important for normal pneumococcal growth. Moreover, deletion of tmRNA increased susceptibility to UV irradiation, to exogenous hydrogen peroxide and to antibiotics that inhibit protein synthesis and transcription. However, the ΔssrA strain was more resistant to fluoroquinolones, showing twofold higher MIC values and up to 1000-fold higher survival rates than the wild type. Deletion of SmpB, the other partner in trans-translation, also reduced survival to levofloxacin in a similar extent. Accumulation of intracellular reactive oxygen species associated to moxifloxacin and levofloxacin treatment was also highly reduced (∼100-fold). Nevertheless, the ΔssrA strain showed higher intracellular accumulation of ethidium bromide and levofloxacin than the wild type, suggesting that tmRNA deficiency protects pneumococcal cells from fluoroquinolone-mediated killing. In fact, analysis of chromosome integrity revealed that deletion of tmRNA prevented the fragmentation of the chromosome associated to levofloxacin treatment. Moreover, such protective effect appears to relay mainly on inhibition of protein synthesis, since a similar effect was observed with antibiotics that inhibit that process. The emergence and spread of drug-resistant pneumococci is a matter of concern and these results contribute to a better comprehension of the mechanisms underlying fluoroquinolones action.

Non-coding RNAs as antibiotic targets

Antibiotics inhibit a wide range of essential processes in the bacterial cell, including replication, transcription, translation and cell wall synthesis. In many instances, these antibiotics exert their effects through association with non-coding RNAs. This review highlights many classical antibiotic targets (e.g. rRNAs and the ribosome), explores a number of emerging targets (e.g. tRNAs, RNase P, riboswitches and small RNAs), and discusses the future directions and challenges associated with non-coding RNAs as antibiotic targets.

Analysis of Small RNAs in Streptococcus mutans under Acid Stress-A New Insight for Caries Research

Streptococcus mutans (S. mutans) is the major clinical pathogen responsible for dental caries. Its acid tolerance has been identified as a significant virulence factor for its survival and cariogenicity in acidic conditions. Small RNAs (sRNAs) are recognized as key regulators of virulence and stress adaptation. Here, we constructed three libraries of sRNAs with small size exposed to acidic conditions for the first time, followed by verification using qRT-PCR. The levels of two sRNAs and target genes predicted to be bioinformatically related to acid tolerance were further evaluated under different acid stress conditions (pH 7.5, 6.5, 5.5, and 4.5) at three time points (0.5, 1, and 2 h). Meanwhile, bacterial growth characteristics and vitality were assessed. We obtained 1879 sRNAs with read counts of at least 100. One hundred and ten sRNAs were perfectly mapped to reported msRNAs in S. mutans. Ten out of 18 sRNAs were validated by qRT-PCR. The survival of bacteria declined as the acid was increased from pH 7.5 to 4.5 at each time point. The bacteria can proliferate under each pH except pH 4.5 with time. The levels of sRNAs gradually decreased from pH 7.5 to 5.5, and slightly increased in pH 4.5; however, the expression levels of target mRNAs were up-regulated in acidic conditions than in pH 7.5. These results indicate that some sRNAs are specially induced at acid stress conditions, involving acid adaptation, and provide a new insight into exploring the complex acid tolerance for S. mutans.

Stress Physiology of Lactic Acid Bacteria

Lactic acid bacteria (LAB) are important starter, commensal, or pathogenic microorganisms. The stress physiology of LAB has been studied in depth for over 2 decades, fueled mostly by the technological implications of LAB robustness in the food industry. Survival of probiotic LAB in the host and the potential relatedness of LAB virulence to their stress resilience have intensified interest in the field. Thus, a wealth of information concerning stress responses exists today for strains as diverse as starter (e.g., Lactococcus lactis), probiotic (e.g., several Lactobacillus spp.), and pathogenic (e.g., Enterococcus and Streptococcus spp.) LAB. Here we present the state of the art for LAB stress behavior. We describe the multitude of stresses that LAB are confronted with, and we present the experimental context used to study the stress responses of LAB, focusing on adaptation, habituation, and cross-protection as well as on self-induced multistress resistance in stationary phase, biofilms, and dormancy. We also consider stress responses at the population and single-cell levels. Subsequently, we concentrate on the stress defense mechanisms that have been reported to date, grouping them according to their direct participation in preserving cell energy, defending macromolecules, and protecting the cell envelope. Stress-induced responses of probiotic LAB and commensal/pathogenic LAB are highlighted separately due to the complexity of the peculiar multistress conditions to which these bacteria are subjected in their hosts. Induction of prophages under environmental stresses is then discussed. Finally, we present systems-based strategies to characterize the "stressome" of LAB and to engineer new food-related and probiotic LAB with improved stress tolerance.

DiiA is a novel dimorphic cell wall protein of Streptococcus pneumoniae involved in invasive disease

OBJECTIVES

Many outer multidomain proteins play fundamental virulent roles in an allele-dependent manner. We aimed to investigate the influence of the outer SP1992 protein, here renamed DiiA (Dimorphic invasion-involved A), in pneumococcal disease.

METHODS

The presence and type of diiA allele was screened by PCR in 560 clinical isolates. Isogenic mutants carrying progressive diiA deletions were constructed and checked in mouse models of infection. DiiA binding to human molecules was carried out by surface plasmon resonance.

RESULTS

The diiA gene is exclusive of Streptococcus pneumoniae and included in the core genome. DiiA variants contain one or two imperfect repeats (R1 and R2), an unstructured region and a cell-wall anchor domain. Clonal complexes carrying both repeats were associated with invasive disease, while those carrying R2 preferentially caused non-invasive syndromes in patients with underlying risk factors. Mutants lacking both repeats were less efficient in nasopharyngeal colonization and dissemination from lungs. Moreover, the ΔdiiA defective strain suffered a severe impairment in bacterial proliferation in blood. Purified DiiA bound to collagen and lactoferrin with high affinity.

CONCLUSIONS

DiiA is a distinctive pneumococcal virulence factor contributing to colonization and long-term invasion in this pathogen.

Transcriptome landscape of Lactococcus lactis reveals many novel RNAs including a small regulatory RNA involved in carbon uptake and metabolism

RNA sequencing has revolutionized genome-wide transcriptome analyses, and the identification of non-coding regulatory RNAs in bacteria has thus increased concurrently. Here we reveal the transcriptome map of the lactic acid bacterial paradigm Lactococcus lactis MG1363 by employing differential RNA sequencing (dRNA-seq) and a combination of manual and automated transcriptome mining. This resulted in a high-resolution genome annotation of L. lactis and the identification of 60 cis-encoded antisense RNAs (asRNAs), 186 trans-encoded putative regulatory RNAs (sRNAs) and 134 novel small ORFs. Based on the putative targets of asRNAs, a novel classification is proposed. Several transcription factor DNA binding motifs were identified in the promoter sequences of (a)sRNAs, providing insight in the interplay between lactococcal regulatory RNAs and transcription factors. The presence and lengths of 14 putative sRNAs were experimentally confirmed by differential Northern hybridization, including the abundant RNA 6S that is differentially expressed depending on the available carbon source. For another sRNA, LLMGnc_147, functional analysis revealed that it is involved in carbon uptake and metabolism. L. lactis contains 13% leaderless mRNAs (lmRNAs) that, from an analysis of overrepresentation in GO classes, seem predominantly involved in nucleotide metabolism and DNA/RNA binding. Moreover, an A-rich sequence motif immediately following the start codon was uncovered, which could provide novel insight in the translation of lmRNAs. Altogether, this first experimental genome-wide assessment of the transcriptome landscape of L. lactis and subsequent sRNA studies provide an extensive basis for the investigation of regulatory RNAs in L. lactis and related lactococcal species.

The Mechanisms of Virulence Regulation by Small Noncoding RNAs in Low GC Gram-Positive Pathogens

The discovery of small noncoding regulatory RNAs (sRNAs) in bacteria has grown tremendously recently, giving new insights into gene regulation. The implementation of computational analysis and RNA sequencing has provided new tools to discover and analyze potential sRNAs. Small regulatory RNAs that act by base-pairing to target mRNAs have been found to be ubiquitous and are the most abundant class of post-transcriptional regulators in bacteria. The majority of sRNA studies has been limited to E. coli and other gram-negative bacteria. However, examples of sRNAs in gram-positive bacteria are still plentiful although the detailed gene regulation mechanisms behind them are not as well understood. Strict virulence control is critical for a pathogen’s survival and many sRNAs have been found to be involved in that process. This review outlines the targets and currently known mechanisms of trans-acting sRNAs involved in virulence regulation in various gram-positive pathogens. In addition, their shared characteristics such as CU interaction motifs, the role of Hfq, and involvement in two-component regulators, riboswitches, quorum sensing, or toxin/antitoxin systems are described.

Genome-wide analyses of small non-coding RNAs in streptococci

Streptococci represent a diverse group of Gram-positive bacteria, which colonize a wide range of hosts among animals and humans. Streptococcal species occur as commensal as well as pathogenic organisms. Many of the pathogenic species can cause severe, invasive infections in their hosts leading to a high morbidity and mortality. The consequence is a tremendous suffering on the part of men and livestock besides the significant financial burden in the agricultural and healthcare sectors. An environmentally stimulated and tightly controlled expression of virulence factor genes is of fundamental importance for streptococcal pathogenicity. Bacterial small non-coding RNAs (sRNAs) modulate the expression of genes involved in stress response, sugar metabolism, surface composition, and other properties that are related to bacterial virulence. Even though the regulatory character is shared by this class of RNAs, variation on the molecular level results in a high diversity of functional mechanisms. The knowledge about the role of sRNAs in streptococci is still limited, but in recent years, genome-wide screens for sRNAs have been conducted in an increasing number of species. Bioinformatics prediction approaches have been employed as well as expression analyses by classical array techniques or next generation sequencing. This review will give an overview of whole genome screens for sRNAs in streptococci with a focus on describing the different methods and comparing their outcome considering sRNA conservation among species, functional similarities, and relevance for streptococcal infection.

Small regulatory RNAs in Streptococcus pneumoniae: discovery and biological functions

Streptococcus pneumoniae is a prominent human pathogen responsible for many severe diseases and the leading cause of childhood mortality worldwide. The pneumococcus is remarkably adept at colonizing and infecting different niches in the human body, and its adaptation to dynamic host environment is a central aspect of its pathogenesis. In the last decade, increasing findings have evidenced small RNAs (sRNAs) as vital regulators in a number of important processes in bacteria. In S. pneumoniae, a small antisense RNA was first discovered in the pMV158 plasmid as a copy number regulator. More recently, genome-wide screens revealed that the pneumococcal genome also encodes multiple sRNAs, many of which have important roles in virulence while some are implicated in competence control. The knowledge of the sRNA-mediated regulation in pneumococcus remains very limited, and future research is needed for better understanding of functions and mechanisms. Here, we provide a comprehensive summary of the current knowledge on sRNAs from S. pneumoniae, focusing mainly on the trans-encoded sRNAs.

Genome-wide identification of small RNAs in Bifidobacterium animalis subsp. lactis KLDS 2.0603 and their regulation role in the adaption to gastrointestinal environment

OBJECTIVE

Bifidobacteria are one of the predominant bacterial species in the human gastrointestinal tract (GIT) and play a vital role in the host's health by acting as probiotics. However, how they regulate themselves to adapt to GIT of their host remains unknown.

METHODS

Eighteen bifidobacterial strains were used to analyze their adaptive capacities towards simulated GIT environment. The strain with highest survival rate and adhesion ability was selected for comparative genome as well as transcriptomic analysis.

RESULTS

The Bifidobacterium animalis subsp. lactis KLDS 2.0603 strain was demonstrated to have the highest survival rate and adhesion ability in simulated GIT treatments. The comparative genome analysis revealed that the KLDS 2.0603 has most similar whole genome sequence compared with BB-12 strain. Eleven intergenic sRNAs were identified after genomes prediction and transcriptomic analysis of KLDS 2.0603. Transcriptomic analysis also showed that genes (mainly sRNAs targeted genes) and sRNAs were differentially expressed in different stress conditions, suggesting that sRNAs might play a crucial role in regulating genes involved in the stress resistance of this strain towards environmental changes.

CONCLUSIONS

This study first provided deep and comprehensive insights into the regulation of KLDS 2.0603 strain at transcription and post-transcription level towards environmental.

Genome-wide analyses in bacteria show small-RNA enrichment for long and conserved intergenic regions

Interest in finding small RNAs (sRNAs) in bacteria has significantly increased in recent years due to their regulatory functions. Development of high-throughput methods and more sophisticated computational algorithms has allowed rapid identification of sRNA candidates in different species. However, given their various sizes (50 to 500 nucleotides [nt]) and their potential genomic locations in the 5' and 3' untranslated regions as well as in intergenic regions, identification and validation of true sRNAs have been challenging. In addition, the evolution of bacterial sRNAs across different species continues to be puzzling, given that they can exert similar functions with various sequences and structures. In this study, we analyzed the enrichment patterns of sRNAs in 13 well-annotated bacterial species using existing transcriptome and experimental data. All intergenic regions were analyzed by WU-BLAST to examine conservation levels relative to species within or outside their genus. In total, more than 900 validated bacterial sRNAs and 23,000 intergenic regions were analyzed. The results indicate that sRNAs are enriched in intergenic regions, which are longer and more conserved than the average intergenic regions in the corresponding bacterial genome. We also found that sRNA-coding regions have different conservation levels relative to their flanking regions. This work provides a way to analyze how noncoding RNAs are distributed in bacterial genomes and also shows conserved features of intergenic regions that encode sRNAs. These results also provide insight into the functions of regions surrounding sRNAs and into optimization of RNA search algorithms.

Identification of 23 novel conserved microRNAs in three rice cultivars

Plant microRNAs (miRNAs) play important roles as modulators of gene expression at the post-transcriptional level. Previous studies have shown that high-throughput sequencing is a powerful tool for the identification of miRNAs, and it is believed that many more miRNAs remain to be discovered. Here, we found 23 novel conserved miRNAs from three rice cultivars by high-throughput sequencing and further identified these through subsequent cloning and quantitative real-time polymerase chain reaction (qPCR). Eight of these novel miRNAs were detected with significant signals in the three rice cultivars by northern blotting assays. The quantitative analysis of their expression profiles showed that most of these miRNAs were perfectly or imperfectly negatively correlated with their target genes, which suggests that these miRNAs may play important roles during rice development. This is the first genome-wide investigation of miRNAs from different rice cultivars, and the data obtained expand the known rice miRNA inventory and provide further information about the regulatory roles played by miRNAs in rice development.

A multicopy sRNA of Listeria monocytogenes regulates expression of the virulence adhesin LapB

The multicopy sRNA LhrC of the intracellular pathogen Listeria monocytogenes has been shown to be induced under infection-relevant conditions, but its physiological role and mechanism of action is not understood. In an attempt to pinpoint the exact terms of LhrC expression, cell envelope stress could be defined as a specific inducer of LhrC. In this process, the two-component system LisRK was shown to be indispensable for expression of all five copies of LhrC. lapB mRNA, encoding a cell wall associated protein that was recently identified as an important virulence factor, was disclosed to be directly bound by LhrC leading to an impediment of its translation. Although LhrC binds to Hfq, it does not require the RNA chaperone for stability or lapB mRNA interaction. The mechanism of LhrC-lapB mRNA binding was shown to involve three redundant CU-rich sites and a structural rearrangement in the sRNA. This study represents an extensive depiction of a so far uncharacterized multicopy sRNA and reveals interesting new aspects concerning its regulation, virulence association and mechanism of target binding.

Impact of Hfq on the Bacillus subtilis transcriptome

The RNA chaperone Hfq acts as a central player in post-transcriptional gene regulation in several Gram-negative Bacteria, whereas comparatively little is known about its role in Gram-positive Bacteria. Here, we studied the function of Hfq in Bacillus subtilis, and show that it confers a survival advantage. A comparative transcriptome analysis revealed mRNAs with a differential abundance that are governed by the ResD-ResE system required for aerobic and anaerobic respiration. Expression of resD was found to be up-regulated in the hfq- strain. Furthermore, several genes of the GerE and ComK regulons were de-regulated in the hfq- background. Surprisingly, only six out of >100 known and predicted small RNAs (sRNAs) showed altered abundance in the absence of Hfq. Moreover, Hfq positively affected the transcript abundance of genes encoding type I toxin-antitoxin systems. Taken the moderate effect on sRNA levels and mRNAs together, it seems rather unlikely that Hfq plays a central role in RNA transactions in Bacillus subtilis.

Deep sequencing-based identification of small regulatory RNAs in Synechocystis sp. PCC 6803

Synechocystis sp. PCC 6803 is a genetically tractable model organism for photosynthesis research. The genome of Synechocystis sp. PCC 6803 consists of a circular chromosome and seven plasmids. The importance of small regulatory RNAs (sRNAs) as mediators of a number of cellular processes in bacteria has begun to be recognized. However, little is known regarding sRNAs in Synechocystis sp. PCC 6803. To provide a comprehensive overview of sRNAs in this model organism, the sRNAs of Synechocystis sp. PCC 6803 were analyzed using deep sequencing, and 7,951,189 reads were obtained. High quality mapping reads (6,127,890) were mapped onto the genome and assembled into 16,192 transcribed regions (clusters) based on read overlap. A total number of 5211 putative sRNAs were revealed from the genome and the 4 megaplasmids, and 27 of these molecules, including four from plasmids, were confirmed by RT-PCR. In addition, possible target genes regulated by all of the putative sRNAs identified in this study were predicted by IntaRNA and analyzed for functional categorization and biological pathways, which provided evidence that sRNAs are indeed involved in many different metabolic pathways, including basic metabolic pathways, such as glycolysis/gluconeogenesis, the citrate cycle, fatty acid metabolism and adaptations to environmentally stress-induced changes. The information from this study provides a valuable reservoir for understanding the sRNA-mediated regulation of the complex physiology and metabolic processes of cyanobacteria.

Small regulatory RNAs from low-GC Gram-positive bacteria

Small regulatory RNAs (sRNAs) that act by base-pairing were first discovered in so-called accessory DNA elements—plasmids, phages, and transposons—where they control replication, maintenance, and transposition. Since 2001, a huge body of work has been performed to predict and identify sRNAs in a multitude of bacterial genomes. The majority of chromosome-encoded sRNAs have been investigated in E. coli and other Gram-negative bacteria. However, during the past five years an increasing number of sRNAs were found in Gram-positive bacteria. Here, we outline our current knowledge on chromosome-encoded sRNAs from low-GC Gram-positive species that act by base-pairing, i.e., an antisense mechanism. We will focus on sRNAs with known targets and defined regulatory mechanisms with special emphasis on Bacillus subtilis.

Determination of sRNA expressions by RNA-seq in Yersinia pestis grown in vitro and during infection

Background

Small non-coding RNAs (sRNAs) facilitate host-microbe interactions. They have a central function in the post-transcriptional regulation during pathogenic lifestyles. Hfq, an RNA-binding protein that many sRNAs act in conjunction with, is required for Y. pestis pathogenesis. However, information on how Yersinia pestis modulates the expression of sRNAs during infection is largely unknown.

Methodology and Principal Findings

We used RNA-seq technology to identify the sRNA candidates expressed from Y. pestis grown invitro and in the infected lungs of mice. A total of 104 sRNAs were found, including 26 previously annotated sRNAs, by searching against the Rfam database with 78 novel sRNA candidates. Approximately 89% (93/104) of these sRNAs from Y. pestis are shared with its ancestor Y. pseudotuberculosis. Ninety-seven percent of these sRNAs (101/104) are shared among more than 80 sequenced genomes of 135 Y. pestis strains. These 78 novel sRNAs include 62 intergenic and 16 antisense sRNAs. Fourteen sRNAs were selected for verification by independent Northern blot analysis. Results showed that nine selected sRNA transcripts were Hfq-dependent. Interestingly, three novel sRNAs were identified as new members of the transcription factor CRP regulon. Semi-quantitative analysis revealed that Y. pestis from the infected lungs induced the expressions of six sRNAs including RyhB1, RyhB2, CyaR/RyeE, 6S RNA, RybB and sR039 and repressed the expressions of four sRNAs, including CsrB, CsrC, 4.5S RNA and sR027.

Conclusions and Significance

This study is the first attempt to subject RNA from Y. pestis-infected samples to direct high-throughput sequencing. Many novel sRNAs were identified and the expression patterns of relevant sRNAs in Y. pestis during invitro growth and invivo infection were revealed. The annotated sRNAs accounted for the most abundant sRNAs either expressed in bacteria grown invitro or differentially expressed in the infected lungs. These findings suggested these sRNAs may have important functions in Y. pestis physiology or pathogenesis.

Target evaluation of the non-coding csRNAs reveals a link of the two-component regulatory system CiaRH to competence control in Streptococcus pneumoniae R6

The two-component regulatory system CiaRH of Streptococcus pneumoniae controls 25 genes, five of which specify homologous small non-coding csRNAs (cia-dependent small RNAs). The csRNAs were predicted to act regulatory as base-pairing sRNAs, but their targets have not been identified. By csRNA gene inactivations we established that the major phenotypes associated with a hyperactive CiaRH system, enhanced β-lactam resistance and prevention of genetic competence, are dependent on the csRNAs. Computational target predictions and evaluations by translational fusions identified six genes to be under csRNA control: spr0081, spr0371, spr0159, spr0551, spr1097 and spr2043(comC). Measuring the effect of single csRNAs on three targets indicated that they acted additively. One of the targets, comC(spr2043), encoding the precursor of the competence stimulating pheromone CSP, constitutes a link of CiaRH to competence control. Partially disrupting predicted csRNA-comC complementarity led to strongly diminished repression by the csRNAs and to transformability in a strain with a hyperactive CiaRH. Thus, a hyperactive CiaRH system prevents competence development by csRNA-dependent post-transcriptional repression of CSP production. The csRNAs are also involved in competence regulation in the wild-type strain R6, but their activity is only apparent in the absence of the protease gene htrA, another CiaRH regulon member.

Genome-wide identification of regulatory RNAs in the human pathogen Clostridium difficile

Clostridium difficile is an emergent pathogen, and the most common cause of nosocomial diarrhea. In an effort to understand the role of small noncoding RNAs (sRNAs) in C. difficile physiology and pathogenesis, we used an in silico approach to identify 511 sRNA candidates in both intergenic and coding regions. In parallel, RNA–seq and differential 5′-end RNA–seq were used for global identification of C. difficile sRNAs and their transcriptional start sites at three different growth conditions (exponential growth phase, stationary phase, and starvation). This global experimental approach identified 251 putative regulatory sRNAs including 94 potential trans riboregulators located in intergenic regions, 91 cis-antisense RNAs, and 66 riboswitches. Expression of 35 sRNAs was confirmed by gene-specific experimental approaches. Some sRNAs, including an antisense RNA that may be involved in control of C. difficile autolytic activity, showed growth phase-dependent expression profiles. Expression of each of 16 predicted c-di-GMP-responsive riboswitches was observed, and experimental evidence for their regulatory role in coordinated control of motility and biofilm formation was obtained. Finally, we detected abundant sRNAs encoded by multiple C. difficile CRISPR loci. These RNAs may be important for C. difficile survival in bacteriophage-rich gut communities. Altogether, this first experimental genome-wide identification of C. difficile sRNAs provides a firm basis for future RNome characterization and identification of molecular mechanisms of sRNA–based regulation of gene expression in this emergent enteropathogen.

The emergent human pathogen Clostridium difficile is a major cause of nosocomial diarrhea associated with antibiotic therapy. During the last few years, severe forms of C. difficile infections became more frequent due to the emergence of hypervirulent isolates. Despite intensive studies, many questions regarding the mechanisms controlling C. difficile virulence remain unanswered. We hypothesized that C. difficile, a member of an ancient group of bacteria, might widely use ancestral RNA–based mechanisms to control its gene expression for better adaptation to host conditions. Indeed, using next-generation sequencing technology, we identified a great number and a large diversity of potential RNA regulators in this pathogen. We obtained experimental evidence for regulatory roles of a particular class of regulatory RNAs responding to c-di-GMP, a universal bacterial signaling molecule regulating motility, biofilm formation, and virulence. We also detected abundant small RNA products of recently discovered adaptive prokaryotic immunity CRISPR-Cas systems that might be important for C. difficile survival in gut communities. Our findings suggest that small RNA molecules may play a major role in regulatory processes during C. difficile infection cycle and as such are promising targets of new therapeutic strategies.

Synergies between RNA degradation and trans-translation in Streptococcus pneumoniae: cross regulation and co-transcription of RNase R and SmpB

BACKGROUND

Ribonuclease R (RNase R) is an exoribonuclease that recognizes and degrades a wide range of RNA molecules. It is a stress-induced protein shown to be important for the establishment of virulence in several pathogenic bacteria. RNase R has also been implicated in the trans-translation process. Transfer-messenger RNA (tmRNA/SsrA RNA) and SmpB are the main effectors of trans-translation, an RNA and protein quality control system that resolves challenges associated with stalled ribosomes on non-stop mRNAs. Trans-translation has also been associated with deficiencies in stress-response mechanisms and pathogenicity.

RESULTS

In this work we study the expression of RNase R in the human pathogen Streptococcus pneumoniae and analyse the interplay of this enzyme with the main components of the trans-translation machinery (SmpB and tmRNA/SsrA). We show that RNase R is induced after a 37°C to 15°C temperature downshift and that its levels are dependent on SmpB. On the other hand, our results revealed a strong accumulation of the smpB transcript in the absence of RNase R at 15°C. Transcriptional analysis of the S. pneumoniae rnr gene demonstrated that it is co-transcribed with the flanking genes, secG and smpB. Transcription of these genes is driven from a promoter upstream of secG and the transcript is processed to yield mature independent mRNAs. This genetic organization seems to be a common feature of Gram positive bacteria, and the biological significance of this gene cluster is further discussed.

CONCLUSIONS

This study unravels an additional contribution of RNase R to the trans-translation system by demonstrating that smpB is regulated by this exoribonuclease. RNase R in turn, is shown to be under the control of SmpB. These proteins are therefore mutually dependent and cross-regulated. The data presented here shed light on the interactions between RNase R, trans-translation and cold-shock response in an important human pathogen.

Predicting sRNAs and their targets in bacteria

Bacterial small RNAs (sRNAs) are an emerging class of regulatory RNAs of about 40–500 nucleotides in length and, by binding to their target mRNAs or proteins, get involved in many biological processes such as sensing environmental changes and regulating gene expression. Thus, identification of bacterial sRNAs and their targets has become an important part of sRNA biology. Current strategies for discovery of sRNAs and their targets usually involve bioinformatics prediction followed by experimental validation, emphasizing a key role for bioinformatics prediction. Here, therefore, we provided an overview on prediction methods, focusing on the merits and limitations of each class of models. Finally, we will present our thinking on developing related bioinformatics models in future.