Identification of small RNAs in diverse bacterial species

Genome-wide annotation of transcript boundaries using bacterial Rend-seq datasets

Accurate annotation to single-nucleotide resolution of the transcribed regions in genomes is key to optimally analyse RNA-seq data, understand regulatory events and for the design of experiments. However, currently most genome annotations provided by GenBank generally lack information about untranslated regions. Additionally, information regarding genomic locations of non-coding RNAs, such as sRNAs, or anti-sense RNAs is frequently missing. To provide such information, diverse RNA-seq technologies, such as Rend-seq, have been developed and applied to many bacterial species. However, incorporating this vast amount of information into annotation files has been limited and is bioinformatically challenging, resulting in UTRs and other non-coding elements being overlooked or misrepresented. To overcome this problem, we present pyRAP (python Rend-seq Annotation Pipeline), a software package that analyses Rend-seq datasets to accurately resolve transcript boundaries genome-wide. We report the use of pyRAP to find novel transcripts, transcript isoforms, and RNase-dependent sRNA processing events. In Bacillus subtilis we uncovered 63 novel transcripts and provide genomic coordinates with single-nucleotide resolution for 2218 5'UTRs, 1864 3'UTRs and 161 non-coding RNAs. In Escherichia coli, we report 117 novel transcripts, 2429 5'UTRs, 1619 3'UTRs and 91 non-coding RNAs, and in Staphylococcus aureus, 16 novel transcripts, 664 5'UTRs, 696 3'UTRs, and 81 non-coding RNAs. Finally, we use pyRAP to produce updated annotation files for B. subtilis 168, E. coli K-12 MG1655, and S. aureus 8325 for use in the wider microbial genomics research community.

Identification and initial characterization of Hfq-associated sRNAs in Histophilus somni strain 2336

Small RNAs (sRNA), in association with the global chaperone regulator Hfq, positively or negatively regulate gene expression in bacteria. For this study, Histophilus somni sRNAs that bind to Hfq were identified and then partially characterized. The Hfq-associated sRNAs in H. somni were isolated and identified by co-immunoprecipitation using anti-Hfq antibody, followed by sRNA sequencing. Sequence analysis of the sRNA samples identified 100 putative sRNAs, out of which 16 were present in pathogenic strain 2336, but not in non-pathogenic strain 129Pt. Bioinformatic analyses suggested that the sRNAs HS9, HS79, and HS97 could bind to many genes putatively involved in virulence/biofilm formation. Furthermore, multi-sequence alignment of the sRNA regions in the genome revealed that HS9 and HS97 could interact with sigma 54, which is a transcription factor linked to important bacterial traits, including motility, virulence, and biofilm formation. Northern blotting was used to determine the approximate size, abundance and any processing events attributed to the sRNAs. Selected sRNA candidates were confirmed to bind Hfq, as determined by electrophoretic mobility shift assays using sRNAs synthesized by in vitro transcription and recombinant Hfq. The exact transcriptional start site of the sRNA candidates was determined by RNA ligase-mediated rapid amplification of cDNA ends, followed by cloning and sequencing. This is the first investigation of H. somni sRNAs that show they may have important regulatory roles in virulence and biofilm formation.

Transcription start site mapping and small RNA profiling of Leptospira biflexa serovar Patoc

Leptospirosis is an emerging zoonotic disease caused by bacterial species of the genus Leptospira. However, the regulatory mechanisms and pathways underlying the adaptation of pathogenic and non-pathogenic Leptospira spp. in different environmental conditions remain elusive. Leptospira biflexa is a non-pathogenic species of Leptospira that lives exclusively in a natural environment. It is an ideal model not only for exploring molecular mechanisms underlying the environmental survival of Leptospira species but also for identifying virulence factors unique to Leptospira's pathogenic species. In this study, we aim to establish the transcription start site (TSS) landscape and the small RNA (sRNA) profile of L. biflexa serovar Patoc grown to exponential and stationary phases via differential RNA-seq (dRNA-seq) and small RNA-seq (sRNA-seq) analyses, respectively. Our dRNA-seq analysis uncovered a total of 2726 TSSs, which are also used to identify other elements, e.g., promoter and untranslated regions (UTRs). Besides, our sRNA-seq analysis revealed a total of 603 sRNA candidates, comprising 16 promoter-associated sRNAs, 184 5'UTR-derived sRNAs, 230 true intergenic sRNAs, 136 5'UTR-antisense sRNAs, and 130 open reading frame (ORF)-antisense sRNAs. In summary, these findings reflect the transcriptional complexity of L. biflexa serovar Patoc under different growth conditions and help to facilitate our understanding of regulatory networks in L. biflexa. To the best of our knowledge, this is the first study reporting the TSS landscape of L. biflexa. The TSS and sRNA landscapes of L. biflexa can also be compared with its pathogenic counterparts, e.g., L. borgpetersenii and L. interrogans, to identify features contributing to their environmental survival and virulence.

Identification and Regulatory Roles of a New Csr Small RNA from Arctic Pseudoalteromonas fuliginea BSW20308 in Temperature Responses

Small RNAs (sRNAs) play a very important role in gene regulation at the posttranscriptional level. However, sRNAs from nonmodel microorganisms, extremophiles in particular, have been rarely explored. We discovered a putative sRNA, termed Pf1 sRNA, in Pseudoalteromonas fuliginea BSW20308 isolated from the polar regions in our previous work. In this study, we identified the sRNA and investigated its regulatory role in gene expression under different temperatures. Pf1 sRNA was confirmed to be a new member of the CsrB family but has little sequence similarity with Escherichia coli CsrB. However, Pf1 sRNA was able to bind to CsrA from E. coli and P. fuliginea BSW20308 to regulate glycogen synthesis. The Pf1 sRNA knockout strain (ΔPf1) affected motility, fitness, and global gene expression in transcriptomes and proteomes at 4°C and 32°C. Genes related to carbon metabolism, amino acid metabolism, salinity tolerance, antibiotic resistance, oxidative stress, motility, chemotaxis, biofilm, and secretion systems were differentially expressed in the wild-type strain and the ΔPf1 mutant. Our study suggested that Pf1 sRNA might play an important role in response to environmental changes by regulating global gene expression. Specific targets of the Pf1 sRNA-CsrA system were tentatively proposed, such as genes involved in the type VI secretion system, TonB-dependent receptors, and response regulators, but most of them have an unknown function. Since this is the first study of CsrB family sRNA in Pseudoalteromonas and microbes from the polar regions, it provides a novel insight at the posttranscriptional level into the responses and adaptation to temperature changes in bacteria from extreme environments. This study also sheds light on the evolution of sRNA in extreme environments and expands the bacterial sRNA database. IMPORTANCE Previous research on microbial temperature adaptation has focused primarily on functional genes, with little attention paid to posttranscriptional regulation. Small RNAs, the major posttranscriptional modulators of gene expression, are greatly underexplored, especially in nonpathogenic and nonmodel microorganisms. In this study, we verified the first Csr sRNA, named Pf1 sRNA, from Pseudoalteromonas, a model genus for studying cold adaptation. We revealed that Pf1 sRNA played an important role in global regulation and was indispensable in improving fitness. This study provided us a comprehensive view of sRNAs from Pseudoalteromonas and expanded our understanding of bacterial sRNAs from extreme environments.

High-Resolution Small RNAs Landscape Provides Insights into Alkane Adaptation in the Marine Alkane-Degrader Alcanivorax dieselolei B-5

Alkanes are widespread in the ocean, and Alcanivorax is one of the most ubiquitous alkane-degrading bacteria in the marine ecosystem. Small RNAs (sRNAs) are usually at the heart of regulatory pathways, but sRNA-mediated alkane metabolic adaptability still remains largely unknown due to the difficulties of identification. Here, differential RNA sequencing (dRNA-seq) modified with a size selection (~50-nt to 500-nt) strategy was used to generate high-resolution sRNAs profiling in the model species Alcanivorax dieselolei B-5 under alkane (n-hexadecane) and non-alkane (acetate) conditions. As a result, we identified 549 sRNA candidates at single-nucleotide resolution of 5'-ends, 63.4% of which are with transcription start sites (TSSs), and 36.6% of which are with processing sites (PSSs) at the 5'-ends. These sRNAs originate from almost any location in the genome, regardless of intragenic (65.8%), antisense (20.6%) and intergenic (6.2%) regions, and RNase E may function in the maturation of sRNAs. Most sRNAs locally distribute across the 15 reference genomes of Alcanivorax, and only 7.5% of sRNAs are broadly conserved in this genus. Expression responses to the alkane of several core conserved sRNAs, including 6S RNA, M1 RNA and tmRNA, indicate that they may participate in alkane metabolisms and result in more actively global transcription, RNA processing and stresses mitigation. Two novel CsrA-related sRNAs are identified, which may be involved in the translational activation of alkane metabolism-related genes by sequestering the global repressor CsrA. The relationships of sRNAs with the characterized genes of alkane sensing (ompS), chemotaxis (mcp, cheR, cheW2), transporting (ompT1, ompT2, ompT3) and hydroxylation (alkB1, alkB2, almA) were created based on the genome-wide predicted sRNA-mRNA interactions. Overall, the sRNA landscape lays the ground for uncovering cryptic regulations in critical marine bacterium, among which both the core and species-specific sRNAs are implicated in the alkane adaptive metabolisms.

Identification of sigma factor 54-regulated small non-coding RNAs by employing genome-wide and transcriptome-based methods in rhizobium strains

Rhizobium-legume symbiosis is considered as the major contributor of biological nitrogen fixation. Bacterial small non-coding RNAs are crucial regulators in several cellular adaptation processes that occur due to the changes in metabolism, physiology, or the external environment. Identifying and analysing the conditional specific/sigma factor-54 regulated sRNAs provides a better understanding of sRNA regulation/mechanism in symbiotic association. In the present study, we have identified sigma factor 54-regulated sRNAs from the genome of six rhizobium strains and from the RNA-seq data of free-living and symbiotic conditions of Bradyrhizobium diazoefficiens USDA 110 to identify the novel putative sRNAs that are over expressed during the regulation of nitrogen fixation. A total of 1351 sRNAs were predicted from the genome of six rhizobium strains and 1375 sRNAs were predicted from the transcriptome data of B. diazoefficiens USDA 110. Analysis of target mRNA for these novel sRNAs was inferred to target several nodulation and nitrogen fixation genes including nodC, nodJ, nodY, nodJ, nodM, nodW, nodZ, nifD, nifN, nifQ, fixK, fixL, fdx, nolB, and several cytochrome proteins. In addition, sRNAs of B. diazoefficiens USDA 110 which targeted the regulatory genes of nitrogen fixation were confirmed by wet-lab experiments with semi-quantitative reverse transcription polymerase chain reaction. Predicted target mRNAs were functionally classified based on the COG analysis and GO annotations. The genome-wide and transcriptome-based integrated methods have led to the identification of several sRNAs involved in the nodulation and symbiosis. Further validation of the functional role of these sRNAs can help in exploring the role of sRNAs in nitrogen metabolism during free-living and symbiotic association with legumes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03394-x.

CRISPRtracrRNA: robust approach for CRISPR tracrRNA detection

MOTIVATION

The CRISPR-Cas9 system is a Type II CRISPR system that has rapidly become the most versatile and widespread tool for genome engineering. It consists of two components, the Cas9 effector protein, and a single guide RNA that combines the spacer (for identifying the target) with the tracrRNA, a trans-activating small RNA required for both crRNA maturation and interference. While there are well-established methods for screening Cas effector proteins and CRISPR arrays, the detection of tracrRNA remains the bottleneck in detecting Class 2 CRISPR systems.

RESULTS

We introduce a new pipeline CRISPRtracrRNA for screening and evaluation of tracrRNA candidates in genomes. This pipeline combines evidence from different components of the Cas9-sgRNA complex. The core is a newly developed structural model via covariance models from a sequence-structure alignment of experimentally validated tracrRNAs. As additional evidence, we determine the terminator signal (required for the tracrRNA transcription) and the RNA-RNA interaction between the CRISPR array repeat and the 5'-part of the tracrRNA. Repeats are detected via an ML-based approach (CRISPRidenify). Providing further evidence, we detect the cassette containing the Cas9 (Type II CRISPR systems) and Cas12 (Type V CRISPR systems) effector protein. Our tool is the first for detecting tracrRNA for Type V systems.

AVAILABILITY AND IMPLEMENTATION

The implementation of the CRISPRtracrRNA is available on GitHub upon requesting the access permission, (https://github.com/BackofenLab/CRISPRtracrRNA). Data generated in this study can be obtained upon request to the corresponding person: Rolf Backofen (backofen@informatik.uni-freiburg.de).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The "guiding" principles of noncoding RNA function

The human genome is pervasively transcribed and yet only a small fraction of these RNAs (less than 2%) are known to code for proteins. The vast majority of the RNAs are classified as noncoding RNAs (ncRNAs) and are further subgrouped as small (shorter than 200 bases) and long noncoding RNAs. The ncRNAs have been identified in all three domains of life and regulate diverse cellular processes through transcriptional and posttranscriptional gene regulation. Most of these RNAs work in conjunction with proteins forming a wide array of base pairing interactions. The determinants of these base pairing interactions are now becoming more evident and show striking similarities among the diverse group of ncRNAs. Here we present a mechanistic overview of pairing between RNA-RNA or RNA-DNA that dictates the function of ncRNAs; we provide examples to illustrate that ncRNAs work through shared evolutionary mechanisms that encompasses a guide-target interaction, involving not only classical Watson-Crick but also noncanonical Wobble and Hoogsteen base pairing. We also highlight the similarities in target selection, proofreading, and the ruler mechanism of ncRNA-protein complexes that confers target specificity and target site selection. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs RNA-Based Catalysis > RNA-Mediated Cleavage RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution.

Small RNAs that target G-rich sequences are generated by diverse biogenesis pathways in Epsilonproteobacteria

Bacterial small RNAs (sRNAs) are widespread post-transcriptional regulators controlling bacterial stress responses and virulence. Nevertheless, little is known about how they arise and evolve. Homologues can be difficult to identify beyond the strain level using sequence-based approaches, and similar functionalities can arise by convergent evolution. Here, we found that the virulence-associated CJnc190 sRNA of the foodborne pathogen Campylobacter jejuni resembles the RepG sRNA from the gastric pathogen Helicobacter pylori. However, while both sRNAs bind G-rich sites in their target mRNAs using a C/U-rich loop, they largely differ in their biogenesis. RepG is transcribed from a stand-alone gene and does not require processing, whereas CJnc190 is transcribed from two promoters as precursors that are processed by RNase III and also has a cis-encoded antagonist, CJnc180. By comparing CJnc190 homologues in diverse Campylobacter species, we show that RNase III-dependent processing of CJnc190 appears to be a conserved feature even outside of C. jejuni. We also demonstrate the CJnc180 antisense partner is expressed in C. coli, yet here might be derived from the 3'UTR of the upstream flagella-related gene. Our analysis of G-tract targeting sRNAs in Epsilonproteobacteria demonstrates that similar sRNAs can have markedly different biogenesis pathways.

Function of Non-coding RNA in Helicobacter pylori-Infected Gastric Cancer

Gastric cancer is a common malignant tumor of the digestive system. Its occurrence and development are the result of a combination of genetic, environmental, and microbial factors. Helicobacter pylori infection is a chronic infection that is closely related to the occurrence of gastric tumorigenesis. Non-coding RNA has been demonstrated to play a very important role in the organism, exerting a prominent role in the carcinogenesis, proliferation, apoptosis, invasion, metastasis, and chemoresistance of tumor progression. H. pylori infection affects the expression of non-coding RNA at multiple levels such as genetic polymorphisms and signaling pathways, thereby promoting or inhibiting tumor progression or chemoresistance. This paper mainly introduces the relationship between H. pylori-infected gastric cancer and non-coding RNA, providing a new perspective for gastric cancer treatment.

Identification of a non-coding RNA and its putative involvement in the regulation of tetanus toxin synthesis in Clostridium tetani

Clostridium tetani produces the tetanus toxin (TeNT), one of the most powerful bacterial toxins known to humankind and responsible for tetanus. The regulation of toxin expression is complex and involves the alternative sigma factor TetR as well as other regulators. Here, a transcriptional analysis of the TeNT-encoding large plasmid of C. tetani identified a putative non-coding small RNA (sRNA), located in close vicinity of the 3′ untranslated region of the tent gene. A northern blot experiment could identify a respective sRNA with a size of approx. 140 nucleotides. Sequence analysis showed that the sRNA contains a 14-nucleotide region that is complementary to a 5′ located region of tent. In order to investigate the function of the sRNA, we applied a RNA interference approach targeting the sRNA in two C. tetani wild-type strains; the constructed antisense C. tetani strains showed an approx. threefold increase in both extracellular and total TeNT production compared to the respective wild-type strains. In addition, recombinant C. tetani strains were constructed that contained tent-locus harboring plasmids with and without the sRNA. However, the introduction of the tent-locus without the sRNA in a C. tetani strain lacking the wild-type TeNT-encoding large plasmid resulted in a lower TeNT production compared to the same strain with recombinant tent-locus with the sRNA. This suggests that the expression or the effect of the sRNA is modulated by the C. tetani genetic background, notably that of the wild-type TeNT-encoding large plasmid. In addition, some recombinant strains exhibited modulated growth patterns, characterized by premature bacterial cell lysis. Taken together, our data indicate that the sRNA acts as a negative regulator of TeNT synthesis, with a possible impact on the growth of C. tetani. We hypothesize that the role of this sRNA is to limit toxin levels in the exponential growth phase in order to prevent premature bacterial lysis.

Optimization and Molecular Characterization of Lipase Producing Bacillus subtilis Strain Rcps3 and Bacillus fumarioli Strain Rcps4 from Oil-contaminated Soils of Warangal

Microbial sources are regularly used as reliable biocatalysts sources which are often used in the process and production industry. Demands for such organisms with greater capacity of intended enzyme production are on the rise. Lipase is important enzyme used in the biotechnological process of hydrolysis of fats in almost all the relevant industries We have utilized the local oil-contaminated soil resources to search for efficacious bacterial strains that have excellent lipase activity. We were successful in identifying two such bacterial sources, namely, Bacillus subtilis strain RCPS3 and Bacillus fumarioli strain RCPS4, responsible for lipase production from oil effluent contaminated soil of Telangana. This is the first report of these two strains from this part of India that are involved in lipase production. The strains were isolated, optimized, and purified using standard microbiology protocols and were characterized at the molecular level using the biomarker 16s ribosomal RNA genes of the strains. The identified and isolated bacterial strains were confirmed as Bacillus subtilis strain RCPS3, and Bacillus fumarioli strain RCPS4 through molecular and computational characterization.

OCCAM: prediction of small ORFs in bacterial genomes by means of a target-decoy database approach and machine learning techniques

Small open reading frames (ORFs) have been systematically disregarded by automatic genome annotation. The difficulty in finding patterns in tiny sequences is the main reason that makes small ORFs to be overlooked by computational procedures. However, advances in experimental methods show that small proteins can play vital roles in cellular activities. Hence, it is urgent to make progress in the development of computational approaches to speed up the identification of potential small ORFs. In this work, our focus is on bacterial genomes. We improve a previous approach to identify small ORFs in bacteria. Our method uses machine learning techniques and decoy subject sequences to filter out spurious ORF alignments. We show that an advanced multivariate analysis can be more effective in terms of sensitivity than applying the simplistic and widely used e-value cutoff. This is particularly important in the case of small ORFs for which alignments present higher e-values than usual. Experiments with control datasets show that the machine learning algorithms used in our method to curate significant alignments can achieve average sensitivity and specificity of 97.06% and 99.61%, respectively. Therefore, an important step is provided here toward the construction of more accurate computational tools for the identification of small ORFs in bacteria.

Whole genome sequence of moderate halophilic marine bacterium Marinobacter litoralis SW-45: Abundance of non-coding RNAs

A report on the de novo Whole Genome Sequence (WGS) of Marinobacter litoralis SW-45, a moderately salt-tolerant bacterium isolated from the seawater in Malaysia is presented. The strain has a genome size of 3.45 Mb and is capable of producing halophilic lipase, protease and esterase enzymes. Computational prediction of non-coding RNA (ncRNA) genes in M. litoralis SW-45 was performed using standalone software known as the non-coding RNA characterization (nocoRNAc). In addition, a phylogenetic tree showing the evolutionary relationship between the strain and other members of the genus Marinobacter was constructed using 16SrRNA sequence information. A total of 385 ncRNA transcripts, 1124 terminator region, and 2350 Stress Induced Duplex Destabilization sites were predicted. The current WGS shotgun project has provided the relevant genetic information that may be useful for the strain's improvement studies. This manuscript gives the first description of M. litoralis with a complete genome.

Regulation of Transcription Termination of Small RNAs and by Small RNAs: Molecular Mechanisms and Biological Functions

Accurate and efficient transcription termination is an important step for cells to generate functional RNA transcripts. In bacteria, two mechanisms are responsible for terminating transcription: intrinsic (Rho-independent) termination and Rho-dependent termination. Growing examples suggest that neither type of transcription termination is static, but instead are highly dynamic and regulated. Regulatory small RNAs (sRNAs) are key players in bacterial stress responses, are frequently expressed under specific growth conditions, and are predominantly terminated through the intrinsic termination mechanism. Once made, sRNAs can base-pair with mRNA targets and regulate mRNA translation and stability. Recent findings suggest that alterations in the efficiency of intrinsic termination for sRNAs under various growth conditions may affect the availability of a given sRNA and the ability of the sRNA to function properly. Moreover, alterations of mRNA structure, translation, and accessibility by sRNAs have the potential to impact the access of Rho factor to mRNAs and thus termination of the mRNA. Indeed, recent studies have revealed that some sRNAs can modulate target gene expression by stimulating or inhibiting Rho-dependent termination, thus expanding the regulatory power of bacterial sRNAs. Here we review the current knowledge on intrinsic termination of sRNAs and sRNA-mediated Rho-dependent termination of protein coding genes in bacteria.

Computational prediction of microRNAs in marine bacteria of the genus Thalassospira

MicroRNAs (miRNAs) are key players in regulation of gene expression at post-transcription level in eukaryotic cells. MiRNAs have been intensively studied in plants, animals and viruses. The investigations of bacterial miRNAs have gained less attention, except for the recent studies on miRNAs derived from Streptococcus mutans ATCC 25175 and Escherichia coli DH10B. In this study, high-throughput sequencing approach was employed to investigate the miRNA population in bacteria of the genus Thalassospira using both the miRDeep2 and CID-miRNA methods. A total of 984 putative miRNAs were identified in 9 species of the genus Thalassospira using both miRDeep and CID-miRNA analyses. Fifty seven conserved putative miRNAs were found in different species of the genus Thalassospira, and up to 6 miRNAs were found to be present at different locations in the T. alkalitolerans JCM 18968^T, T. lucentensis QMT2^T and T. xianhensis P-4^T. None of the putative miRNAs was found to share sequence to the reported miRNAs in E. coli DH10B and S. mutans ATCC 25175. The findings provide a comprehensive list of computationally identified miRNAs in 9 bacterial species of the genus Thalassospira and addressed the existing knowledge gap on the presence of miRNAs in the Thalassospira genomes.

Synthetic Biology of Small RNAs and Riboswitches

In bacteria and archaea, small RNAs (sRNAs) regulate complex networks through antisense interactions with target mRNAs in trans, and riboswitches regulate gene expression in cis based on the ability to bind small-molecule ligands. Although our understanding and characterization of these two important regulatory RNA classes is far from complete, these RNA-based mechanisms have proven useful for a wide variety of synthetic biology applications. Besides classic and contemporary applications in the realm of metabolic engineering and orthogonal gene control, this review also covers newer applications of regulatory RNAs as biosensors, logic gates, and tools to determine RNA-RNA interactions. A separate section focuses on critical insights gained and challenges posed by fundamental studies of sRNAs and riboswitches that should aid future development of synthetic regulatory RNAs.

Bacterial Small RNAs in the Genus Herbaspirillum spp

The genus Herbaspirillum includes several strains isolated from different grasses. The identification of non-coding RNAs (ncRNAs) in the genus Herbaspirillum is an important stage studying the interaction of these molecules and the way they modulate physiological responses of different mechanisms, through RNA⁻RNA interaction or RNA⁻protein interaction. This interaction with their target occurs through the perfect pairing of short sequences (cis-encoded ncRNAs) or by the partial pairing of short sequences (trans-encoded ncRNAs). However, the companion Hfq can stabilize interactions in the trans-acting class. In addition, there are Riboswitches, located at the 5' end of mRNA and less often at the 3' end, which respond to environmental signals, high temperatures, or small binder molecules. Recently, CRISPR (clustered regularly interspaced palindromic repeats), in prokaryotes, have been described that consist of serial repeats of base sequences (spacer DNA) resulting from a previous exposure to exogenous plasmids or bacteriophages. We identified 285 ncRNAs in Herbaspirillum seropedicae (H. seropedicae) SmR1, expressed in different experimental conditions of RNA-seq material, classified as cis-encoded ncRNAs or trans-encoded ncRNAs and detected RNA riboswitch domains and CRISPR sequences. The results provide a better understanding of the participation of this type of RNA in the regulation of the metabolism of bacteria of the genus Herbaspirillum spp.

Algae induce siderophore biosynthesis in the freshwater bacterium Cupriavidus necator H16

Cupriachelin is a photoreactive lipopeptide siderophore produced by the freshwater bacterium Cupriavidus necator H16. In the presence of sunlight, the iron-loaded siderophore undergoes photolytic cleavage, thereby releasing solubilized iron into the environment. This iron is not only available to the siderophore producer, but also to the surrounding microbial community. In this study, the cupriachelin-based interaction between C. necator H16 and the freshwater diatom Navicula pelliculosa was investigated. A reporter strain of the bacterium was constructed to study differential expression levels of the cupriachelin biosynthesis gene cucJ in response to varying environmental conditions. Not only iron starvation, but also culture supernatants of N. pelliculosa were found to induce cupriachelin biosynthesis. The transcription factors involved in this differential gene expression were identified using DNA-protein pulldown assays. Besides the well-characterized ferric uptake regulator, a two-component system was found to tune the expression of cupriachelin biosynthesis genes in the presence of diatom supernatants.

Regulatory sRNAs in Cyanobacteria

As the transcriptional and post-transcriptional regulators of gene expression, small RNAs (sRNAs) play important roles in every domain of life in organisms. It has been discovered gradually that bacteria possess multiple means of gene regulation using RNAs. They have been continuously used as model organisms for photosynthesis, metabolism, biotechnology, evolution, and nitrogen fixation for many decades. Cyanobacteria, one of the most ancient life forms, constitute all kinds of photoautotrophic bacteria and exist in almost any environment on this planet. It is believed that a complex RNA-based regulatory mechanism functions in cyanobacteria to help them adapt to changes and stresses in diverse environments. Although lagging far behind other model microorganisms, such as yeast and Escherichia coli, more and more non-coding regulatory sRNAs have been recognized in cyanobacteria during the past decades. In this article, by focusing on cyanobacterial sRNAs, the approaches for detection and targeting of sRNAs will be summarized, four major mechanisms and regulatory functions will be generalized, eight types of cis-encoded sRNA and four types of trans-encoded sRNAs will be reviewed in detail, and their possible physiological functions will be further discussed.

Two virulent sRNAs identified by genomic sequencing target the type III secretion system in rice bacterial blight pathogen

BACKGROUND

Small non-coding RNA (sRNA) short sequences regulate various biological processes in all organisms, including bacteria that are animal or plant pathogens. Virulent or pathogenicity-associated sRNAs have been increasingly elucidated in animal pathogens but little is known about similar category of sRNAs in plant-pathogenic bacteria. This is particularly true regarding rice bacterial blight pathogen Xanthomonas oryzae pathovar oryzae (Xoo) as studies on the virulent role of Xoo sRNAs is very limited at present.

RESULTS

The number and genomic distribution of sRNAs in Xoo were determined by bioinformatics analysis based on high throughput sequencing (sRNA-Seq) of the bacterial cultures from virulence-inducing and standard growth media, respectively. A total of 601 sRNAs were identified in the Xoo genome and ten virulent sRNA candidates were screened out based on significant differences of their expression levels between the culture conditions. In addition, trans3287 and trans3288 were also selected as candidates due to high expression levels in both media. The differential expression of 12 sRNAs evidenced by the sRNA-Seq data was confirmed by a convincing quantitative method. Based on genetic analysis of Xoo ΔsRNA mutants generated by deletion of the 12 single sRNAs, trans217 and trans3287 were characterized as virulent sRNAs. They are essential not only for the formation of bacterial blight in a susceptible rice variety Nipponbare but also for the induction of hypersensitive response (HR) in nonhost plant tobacco. Xoo Δtrans217 and Δtrans3287 mutants fail to induce bacterial blight in Nipponbare and also fail to induce the HR in tobacco, whereas, genetic complementation restores both mutants to the wild type in the virulent performance and HR induction. Similar effects of gene knockout and complementation were found in the expression of hrpG and hrpX genes, which encode regulatory proteins of the type III secretion system. Consistently, secretion of a type III effector, PthXo1, is blocked in Δtrans217 or Δtrans3287 bacterial cultures but retrieved by genetic complementation to both mutants.

CONCLUSIONS

The genetic analysis characterizes trans217 and trans3287 as pathogenicity-associated sRNAs essential for the bacterial virulence on the susceptible rice variety and for the HR elicitation in the nonhost plant. The molecular evidence suggests that both virulent sRNAs regulate the bacterial virulence by targeting the type III secretion system.

Identification of a Novel Small RNA srvg23535 in Vibrio alginolyticus ZJ-T and Its Characterization With Phenotype MicroArray Technology

Small non-coding RNAs (sRNAs) are important modulators of gene expression and are involved in the pathogenesis and survival of prokaryotes. However, few studies have been conducted with Vibrio alginolyticus, which limits our ability to probe the regulation of virulence and environmental adaptation by sRNAs in this opportunistic pathogen. In this study, the sRNA candidate srvg23535 was identified in V. alginolyticus ZJ-T. The precise transcript end, secondary structure, and sequence conservation were determined. A srvg23535 null mutant was constructed and characterized by using Phenotype MicroArray (PM) technology. In silico target prediction was conducted by IntaRNA and TargetRNA2. Subsequently, a 107 nt transcript was validated with a sigma70 promoter at the 5' end and a Rho-independent terminator at the 3' end. The sRNA srvg23535 had four stem-loop structures and was conserved among Vibrio harveyi, Vibrio parahaemolyticus, and Vibrio splendidus. Deletion of srvg23535 in V. alginolyticus ZJ-T led to a weaker utilization of D-mannose, D-melibiose, lactulose, and inosine as carbon sources but stronger utilization of L-cysteine as nitrogen source. Moreover, the srvg2353 mutant showed stronger resistance to osmotic stress but weaker resistance to pH stress. Additionally, a total of 22 common targets were identified and several were related to the observed phenotype of the mutant. This study indicated that the novel sRNA, srvg23535, is conserved and restricted to Vibrio spp., affecting the utilization of several carbon and nitrogen sources and the response to osmotic and pH stress. These results extend our understanding of sRNA regulation in V. alginolyticus and provide a significant resource for the further study of the precise target mRNAs of srvg23535, which may provide targets for antibacterial therapeutic or attenuated vaccines against Vibrio spp.

Advances and prospects in metabolic engineering of Zymomonas mobilis

Biorefinery of biomass-based biofuels and biochemicals by microorganisms is a competitive alternative of traditional petroleum refineries. Zymomonas mobilis is a natural ethanologen with many desirable characteristics, which makes it an ideal industrial microbial biocatalyst for commercial production of desirable bioproducts through metabolic engineering. In this review, we summarize the metabolic engineering progress achieved in Z. mobilis to expand its substrate and product ranges as well as to enhance its robustness against stressful conditions such as inhibitory compounds within the lignocellulosic hydrolysates and slurries. We also discuss a few metabolic engineering strategies that can be applied in Z. mobilis to further develop it as a robust workhorse for economic lignocellulosic bioproducts. In addition, we briefly review the progress of metabolic engineering in Z. mobilis related to the classical synthetic biology cycle of "Design-Build-Test-Learn", as well as the progress and potential to develop Z. mobilis as a model chassis for biorefinery practices in the synthetic biology era.

Identification of Plasmid-Encoded sRNAs in a blaNDM-1-Harboring Multidrug-Resistance Plasmid pNDM-HK in Enterobacteriaceae

Small RNAs (sRNAs) play significant roles in regulating gene expression post-transcriptionally in response to environmental changes in bacteria. In this work, we identified and characterized six novel sRNAs from an emerging multidrug-resistance (MDR) plasmid pNDM-HK, a New Delhi metallo-β-lactamase 1 gene (bla_NDM−1)-carrying IncL/M plasmid that has caused worldwide threat in recent years. These sRNAs are located at different regions of pNDM-HK, such as replication, stability, and variable regions. Moreover, one of the plasmid-encoded sRNAs (NDM-sR3) functions in an Hfq-dependent manner and possibly plays roles in the fitness of pNDM-HK carrying bacteria. In addition, we attempted to construct the phylogenetic tree based on these novel sRNAs and surprisingly, the sRNA-phylogenetic tree provided significant information about the evolutionary pathway of pNDM-HK, including possible gene acquisition and insertion from relevant plasmids. Moreover, the sRNA-phylogenetic tree can specifically cluster the IncM2 type and distinguish it from other IncL/M subtypes. In summary, this is the first study to systematically identify and characterize sRNAs from clinically-isolated MDR plasmids. We believe that these newly found sRNAs could lead to further understanding and new directions to study the evolution and dissemination of the clinically MDR bacterial plasmids.

Small RNA-mediated regulation in bacteria: A growing palette of diverse mechanisms

Small RNAs (sRNAs) in bacteria have evolved with diverse mechanisms to balance their target gene expression in response to changes in the environment. Accumulating studies on bacterial regulatory processes firmly established that sRNAs modulate their target gene expression generally at the posttranscriptional level. Identification of large number of sRNAs by advanced technologies, like deep sequencing, tilling microarray, indicates the existence of a plethora of distinctive sRNA-mediated regulatory mechanisms in bacteria. Types of the novel mechanisms are increasing with the discovery of new sRNAs. Complementary base pairing between sRNAs and target RNAs assisted by RNA chaperones like Hfq and ProQ, in many occasions, to regulate the cognate gene expression is prevalent in sRNA mechanisms. sRNAs, in most studied cases, can directly base pair with target mRNA to remodel its expression. Base pairing can happen either in the untranslated regions or in the coding regions of mRNA to activate/repress its translation. sRNAs also act as target mimic to titrate away different regulatory RNAs from its target. Other mechanism includes the sequestration of regulatory proteins, especially transcription factors, by sRNAs. Numerous sRNAs, following analogous mechanism, are widespread in bacteria, and thus, has drawn immense attention for the development of RNA-based technologies. Nevertheless, typical sRNA mechanisms are also discovered to be confined in some bacteria. Analysis of the sRNA mechanisms unravels their existence in both the single step processes and the complex regulatory networks with a global effect on cell physiology. This review deals with the diverse array of mechanisms, which sRNAs follow to maintain bacterial lifestyle.

A Genome-Wide Prediction and Identification of Intergenic Small RNAs by Comparative Analysis in Mesorhizobium huakuii 7653R

In bacteria, small non-coding RNAs (sRNAs) are critical regulators of cellular adaptation to changes in metabolism, physiology, or the external environment. In the last decade, more than 2000 of sRNA families have been reported in the Rfam database and have been shown to exert various regulatory functions in bacterial transcription and translation. However, little is known about sRNAs and their functions in Mesorhizobium. Here, we predicted putative sRNAs in the intergenic regions (IGRs) of M. huakuii 7653R by genome-wide comparisons with four related Mesorhizobial strains. The expression and transcribed regions of candidate sRNAs were analyzed using a set of high-throughput RNA deep sequencing data. In all, 39 candidate sRNAs were found, with 5 located in the symbiotic megaplasmids and 34 in the chromosome of M. huakuii 7653R. Of these, 24 were annotated as functional sRNAs in the Rfam database and 15 were recognized as putative novel sRNAs. The expression of nine selected sRNAs was confirmed by Northern blotting, and most of the nine selected sRNAs were highly expressed in 28 dpi nodules and under symbiosis-mimicking conditions. For those putative novel sRNAs, functional categorizations of their target genes were performed by analyzing the enriched GO terms. In addition, MH_s15 was shown to be an abundant and conserved sRNA.

Genome-Wide Detection of Small Regulatory RNAs in Deep-Sea Bacterium Shewanella piezotolerans WP3

Shewanella are one of the most abundant Proteobacteria in the deep-sea and are renowned for their versatile electron accepting capacities. The molecular mechanisms involved in their adaptation to diverse and extreme environments are not well understood. Small non-coding RNAs (sRNAs) are known for modulating the gene expression at transcriptional and posttranscriptional levels, subsequently playing a key role in microbial adaptation. To understand the potential roles of sRNAs in the adaptation of Shewanella toward deep-sea environments, here an in silico approach was utilized to detect the sRNAs in the genome of Shewanella piezotolerans WP3, a piezotolerant and psychrotolerant deep-sea iron reducing bacterium. After scanning 3673 sets of 5' and 3' UTRs of orthologous genes, 209 sRNA candidates were identified with high confidence in S. piezotolerans WP3. About 92% (193 out of 209) of these putative sRNAs belong to the class trans-encoded RNAs, suggesting that trans-regulatory RNAs are the dominant class of sRNAs in S. piezotolerans WP3. The remaining 16 cis-regulatory RNAs were validated through quantitative polymerase chain reaction. Five cis-sRNAs were further shown to act as cold regulated sRNAs. Our study provided additional evidence at the transcriptional level to decipher the microbial adaptation mechanisms to extreme environmental conditions.

The malS-5′UTR regulates hisG, a key gene in the histidine biosynthetic pathway in Salmonella enterica serovar Typhi

Bacterial noncoding RNAs (ncRNA) regulate diverse cellular processes, including virulence and environmental fitness. The malS 5′ untranslated region (named malS-5′UTR) was identified as a regulatory ncRNA that increases the invasive capacity of Salmonella enterica serovar Typhi. An IntaRNA search suggested base pairing between malS-5′UTR and hisG mRNA, a key gene in the histidine biosynthetic pathway. Overexpression of malS-5′UTR markedly reduced bacterial growth in minimal medium without histidine. Overexpression of malS-5′UTR increased mRNA from his operon genes, independently of the bax gene, and decreased HisG protein in Salmonella Typhi. RNA structure analysis showed base pairing of the malS-5′UTR RNA with the hisG mRNA across the ribosome binding site. Thus, we propose that malS-5′UTR inhibited hisG translation, probably by base pairing to the Shine–Dalgarno sequence.

Several Hfq-dependent alterations in physiology of Yersinia enterocolitica O:3 are mediated by derepression of the transcriptional regulator RovM

In bacteria, the RNA chaperone Hfq enables pairing of small regulatory RNAs with their target mRNAs and therefore is a key player of post-transcriptional regulation network. As a global regulator, Hfq is engaged in the adaptation to external environment, regulation of metabolism and bacterial virulence. In this study we used RNA-sequencing and quantitative proteomics (LC-MS/MS) to elucidate the role of this chaperone in the physiology and virulence of Yersinia enterocolitica serotype O:3. This global approach revealed the profound impact of Hfq on gene and protein expression. Furthermore, the role of Hfq in the cell morphology, metabolism, cell wall integrity, resistance to external stresses and pathogenicity was evaluated. Importantly, our results revealed that several alterations typical for the hfq-negative phenotype were due to derepression of the transcriptional factor RovM. The overexpression of RovM caused by the loss of Hfq chaperone resulted in extended growth defect, alterations in the lipid A structure, motility and biofilm formation defects, as well as changes in mannitol utilization. Furthermore, in Y. enterocolitica RovM only in the presence of Hfq affected the abundance of RpoS. Finally, the impact of hfq and rovM mutations on the virulence was assessed in the mouse infection model.

Modeling sRNA-Regulated Plasmid Maintenance

We study a theoretical model for the toxin-antitoxin (hok/sok) mechanism for plasmid maintenance in bacteria. Toxin-antitoxin systems enforce the maintenance of a plasmid through post-segregational killing of cells that have lost the plasmid. Key to their function is the tight regulation of expression of a protein toxin by an sRNA antitoxin. Here, we focus on the nonlinear nature of the regulatory circuit dynamics of the toxin-antitoxin mechanism. The mechanism relies on a transient increase in protein concentration rather than on the steady state of the genetic circuit. Through a systematic analysis of the parameter dependence of this transient increase, we confirm some known design features of this system and identify new ones: for an efficient toxin-antitoxin mechanism, the synthesis rate of the toxin’s mRNA template should be lower that of the sRNA antitoxin, the mRNA template should be more stable than the sRNA antitoxin, and the mRNA-sRNA complex should be more stable than the sRNA antitoxin. Moreover, a short half-life of the protein toxin is also beneficial to the function of the toxin-antitoxin system. In addition, we study a therapeutic scenario in which a competitor mRNA is introduced to sequester the sRNA antitoxin, causing the toxic protein to be expressed.

The Small RNA GcvB Promotes Mutagenic Break Repair by Opposing the Membrane Stress Response

Microbes and human cells possess mechanisms of mutagenesis activated by stress responses. Stress-inducible mutagenesis mechanisms may provide important models for mutagenesis that drives host-pathogen interactions, antibiotic resistance, and possibly much of evolution generally. In Escherichia coli, repair of DNA double-strand breaks is switched to a mutagenic mode, using error-prone DNA polymerases, via the SOS DNA damage and general (σ^S) stress responses. We investigated small RNA (sRNA) clients of Hfq, an RNA chaperone that promotes mutagenic break repair (MBR), and found that GcvB promotes MBR by allowing a robust σ^S response, achieved via opposing the membrane stress (σ^E) response. Cells that lack gcvB were MBR deficient and displayed reduced σ^S-dependent transcription but not reduced σ^S protein levels. The defects in MBR and σ^S-dependent transcription in ΔgcvB cells were alleviated by artificially increasing σ^S levels, implying that GcvB promotes mutagenesis by allowing a normal σ^S response. ΔgcvB cells were highly induced for the σ^E response, and blocking σ^E response induction restored both mutagenesis and σ^S-promoted transcription. We suggest that GcvB may promote the σ^S response and mutagenesis indirectly, by promoting membrane integrity, which keeps σ^E levels lower. At high levels, σ^E might outcompete σ^S for binding RNA polymerase and so reduce the σ^S response and mutagenesis. The data show the delicate balance of stress response modulation of mutagenesis.

IMPORTANCE

Mutagenesis mechanisms upregulated by stress responses promote de novo antibiotic resistance and cross-resistance in bacteria, antifungal drug resistance in yeasts, and genome instability in cancer cells under hypoxic stress. This paper describes the role of a small RNA (sRNA) in promoting a stress-inducible-mutagenesis mechanism, mutagenic DNA break repair in Escherichia coli The roles of many sRNAs in E. coli remain unknown. This study shows that ΔgcvB cells, which lack the GcvB sRNA, display a hyperactivated membrane stress response and reduced general stress response, possibly because of sigma factor competition for RNA polymerase. This results in a mutagenic break repair defect. The data illuminate a function of GcvB sRNA in opposing the membrane stress response, and thus indirectly upregulating mutagenesis.

Contribution of bacterial promoter elements to transcription start site detection accuracy

Accurately detecting transcription start sites (TSS) is a starting point for understanding gene transcription, and an important ingredient in a number of applications necessary for functional gene annotation, such as gene and operon predictions. Available methods for TSS detection in bacteria use very different description of the bacterial promoter structure and all of them show low accuracy. It is therefore unclear which promoter features should be included in TSS recognition, and how their accuracy impacts the search detection. We here address this question for [Formula: see text] and [Formula: see text] (an alternative [Formula: see text] factor) promoters in E. coli. We find that [Formula: see text]35 element, which is considered exchangeable, and is often not included in TSS search, contributes to the search accuracy equally (for [Formula: see text], or more (for [Formula: see text] than the ubiquitous [Formula: see text]10 element. Surprisingly, the sequence of the spacer between [Formula: see text]35 and [Formula: see text]10 promoter elements, which is commonly included in TSS detection, significantly decreases the search accuracy for [Formula: see text] promoters. However, the spacer sequence improves the search accuracy for [Formula: see text] promoters, which we attribute to a presence of sequence conservation. Overall, there is as much as [Formula: see text]50% false positive reduction for optimally implemented promoter features in [Formula: see text], underlying necessity for accurate promoter element alignments.

Bacterial Small Regulatory RNAs and Hfq Protein

Small regulatory RNA (sRNA) is a unique noncoding RNA involved in regulation of gene expression in both eukaryotic and bacterial cells. This short review discusses examples of positive and negative translation regulation by sRNAs in bacteria and participation of Hfq in these processes. The importance of structure investigation of nucleotide-protein and RNA-protein complexes for designing a model of Hfq interaction with both mRNA and sRNA simultaneously is demonstrated.

Regulatory RNAs in Bacillus subtilis: a Gram-Positive Perspective on Bacterial RNA-Mediated Regulation of Gene Expression

Bacteria can employ widely diverse RNA molecules to regulate their gene expression. Such molecules include trans-acting small regulatory RNAs, antisense RNAs, and a variety of transcriptional attenuation mechanisms in the 5' untranslated region. Thus far, most regulatory RNA research has focused on Gram-negative bacteria, such as Escherichia coli and Salmonella. Hence, there is uncertainty about whether the resulting insights can be extrapolated directly to other bacteria, such as the Gram-positive soil bacterium Bacillus subtilis. A recent study identified 1,583 putative regulatory RNAs in B. subtilis, whose expression was assessed across 104 conditions. Here, we review the current understanding of RNA-based regulation in B. subtilis, and we categorize the newly identified putative regulatory RNAs on the basis of their conservation in other bacilli and the stability of their predicted secondary structures. Our present evaluation of the publicly available data indicates that RNA-mediated gene regulation in B. subtilis mostly involves elements at the 5' ends of mRNA molecules. These can include 5' secondary structure elements and metabolite-, tRNA-, or protein-binding sites. Importantly, sense-independent segments are identified as the most conserved and structured potential regulatory RNAs in B. subtilis. Altogether, the present survey provides many leads for the identification of new regulatory RNA functions in B. subtilis.

Use of Metatranscriptomics in Microbiome Research

The human intestinal microbiome is a microbial ecosystem that expresses as many as 100 times more genes than the human host, thereby constituting an important component of the human holobiome, which contributes to multiple health and disease processes. As most commensal species are difficult or impossible to culture, genomic characterization of microbiome composition and function, under various environmental conditions, comprises a central tool in understanding its roles in health and disease. The first decade of microbiome research was mainly characterized by usage of DNA sequencing-based 16S rDNA and shotgun metagenome sequencing, allowing for the elucidation of microbial composition and genome structure. Technological advances in RNA-seq have recently provided us with an ability to gain insight into the genes that are actively expressed in complex bacterial communities, enabling the elucidation of the functional changes that dictate the microbiome functions at given contexts, its interactions with the host, and functional alterations that accompany the conversion of a healthy microbiome toward a disease-driving configuration. Here, we highlight some of the key metatranscriptomics strategies that are implemented to determine microbiota gene expression and its regulation and discuss the advantages and potential challenges associated with these approaches.

A Transcriptome Map of Actinobacillus pleuropneumoniae at Single-Nucleotide Resolution Using Deep RNA-Seq

Actinobacillus pleuropneumoniae is the pathogen of porcine contagious pleuropneumoniae, a highly contagious respiratory disease of swine. Although the genome of A. pleuropneumoniae was sequenced several years ago, limited information is available on the genome-wide transcriptional analysis to accurately annotate the gene structures and regulatory elements. High-throughput RNA sequencing (RNA-seq) has been applied to study the transcriptional landscape of bacteria, which can efficiently and accurately identify gene expression regions and unknown transcriptional units, especially small non-coding RNAs (sRNAs), UTRs and regulatory regions. The aim of this study is to comprehensively analyze the transcriptome of A. pleuropneumoniae by RNA-seq in order to improve the existing genome annotation and promote our understanding of A. pleuropneumoniae gene structures and RNA-based regulation. In this study, we utilized RNA-seq to construct a single nucleotide resolution transcriptome map of A. pleuropneumoniae. More than 3.8 million high-quality reads (average length ~90 bp) from a cDNA library were generated and aligned to the reference genome. We identified 32 open reading frames encoding novel proteins that were mis-annotated in the previous genome annotations. The start sites for 35 genes based on the current genome annotation were corrected. Furthermore, 51 sRNAs in the A. pleuropneumoniae genome were discovered, of which 40 sRNAs were never reported in previous studies. The transcriptome map also enabled visualization of 5'- and 3'-UTR regions, in which contained 11 sRNAs. In addition, 351 operons covering 1230 genes throughout the whole genome were identified. The RNA-Seq based transcriptome map validated annotated genes and corrected annotations of open reading frames in the genome, and led to the identification of many functional elements (e.g. regions encoding novel proteins, non-coding sRNAs and operon structures). The transcriptional units described in this study provide a foundation for future studies concerning the gene functions and the transcriptional regulatory architectures of this pathogen.

Towards the complete small RNome of Acinetobacter baumannii

In recent years, the Gram-negative bacterium Acinetobacter baumannii has garnered considerable attention for its unprecedented capacity to rapidly develop resistance to antibacterial therapeutics. This is coupled with the seemingly epidemic emergence of new hyper-virulent strains. Although strain-specific differences for A. baumannii isolates have been well described, these studies have primarily focused on proteinaceous factors. At present, only limited publications have investigated the presence and role of small regulatory RNA (sRNA) transcripts. Herein, we perform such an analysis, describing the RNA-seq-based identification of 78 A. baumannii sRNAs in the AB5075 background. Together with six previously identified elements, we include each of these in a new genome annotation file, which will serve as a tool to investigate regulatory events in this organism. Our work reveals that the sRNAs display high expression, accounting for >50 % of the 20 most strongly expressed genes. Through conservation analysis we identified six classes of similar sRNAs, with one found to be particularly abundant and homologous to regulatory, C4 antisense RNAs found in bacteriophages. These elements appear to be processed from larger transcripts in an analogous manner to the phage C4 molecule and are putatively controlled by two further sRNAs that are strongly antisense to them. Collectively, this study offers a detailed view of the sRNA content of A. baumannii, exposing sequence and structural conservation amongst these elements, and provides novel insight into the potential evolution, and role, of these understudied regulatory molecules.

Discovery and profiling of small RNAs responsive to stress conditions in the plant pathogen Pectobacterium atrosepticum

Background

Small RNAs (sRNAs) have emerged as important regulatory molecules and have been studied in several bacteria. However, to date, there have been no whole-transcriptome studies on sRNAs in any of the Soft Rot Enterobacteriaceae (SRE) group of pathogens. Although the main ecological niches for these pathogens are plants, a significant part of their life cycle is undertaken outside their host within adverse soil environment. However, the mechanisms of SRE adaptation to this harsh nutrient-deficient environment are poorly understood.

Results

In the study reported herein, by using strand-specific RNA-seq analysis and in silico sRNA predictions, we describe the sRNA pool of Pectobacterium atrosepticum and reveal numerous sRNA candidates, including those that are induced during starvation-activated stress responses. Consequently, strand-specific RNA-seq enabled detection of 137 sRNAs and sRNA candidates under starvation conditions; 25 of these sRNAs were predicted for this bacterium in silico. Functional annotations were computationally assigned to 68 sRNAs. The expression of sRNAs in P. atrosepticum was compared under growth-promoting and starvation conditions: 68 sRNAs were differentially expressed with 47 sRNAs up-regulated under nutrient-deficient conditions. Conservation analysis using BLAST showed that most of the identified sRNAs are conserved within the SRE. Subsequently, we identified 9 novel sRNAs within the P. atrosepticum genome.

Conclusions

Since many of the identified sRNAs are starvation-induced, the results of our study suggests that sRNAs play key roles in bacterial adaptive response. Finally, this work provides a basis for future experimental characterization and validation of sRNAs in plant pathogens.

Electronic supplementary material

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

Accelerating Discovery and Functional Analysis of Small RNAs with New Technologies

Over the past decade, bacterial small RNAs (sRNAs) have gone from a biological curiosity to being recognized as a major class of regulatory molecules. High-throughput methods for sampling the transcriptional output of bacterial cells demonstrate that sRNAs are universal features of bacterial transcriptomes, are plentiful, and appear to vary extensively over evolutionary time. With ever more bacteria coming under study, the question becomes how can we accelerate the discovery and functional characterization of sRNAs in diverse organisms. New technologies built on high-throughput sequencing are emerging that can rapidly provide global insight into the numbers and functions of sRNAs in bacteria of interest, providing information that can shape hypotheses and guide research. In this review, we describe recent developments in transcriptomics (RNA-seq) and functional genomics that we expect to help us develop an integrated, systems-level view of sRNA biology in bacteria.

5′-UTR of malS increases the invasive capacity of Salmonella enterica serovar Typhi by influencing the expression of bax

ABSTRACT 

Aim: An RNA-seq analysis recently identified a 236-nucleotide transcript upstream from malS in Salmonella enterica serovar Typhi. Here, we investigated its molecular characteristics and function. Materials & methods: RACE and northern blotting were used to determine the molecular characteristics of the sequence, and mutagenesis, microarray, immunoblotting and an invasion assay were used to investigate the functions of the transcript. Results: The transcript was identified as the malS 5′-untranslated region (UTR), which could influence the expression of the flagellar and SPI-1 genes and the invasion of HeLa cells by S. Typhi. Deletion of bax increased the expression of the invasion genes and the invasive capacity of S. Typhi, whereas the expression of the malS 5′-UTR reduced the expression of bax. Conclusion: The malS 5′-UTR reduces the expression of bax and increases the invasive capacity of S. Typhi.

Rapid and high-throughput construction of microbial cell-factories with regulatory noncoding RNAs

Due to global crises such as pollution and depletion of fossil fuels, sustainable technologies based on microbial cell-factories have been garnering great interest as an alternative to chemical factories. The development of microbial cell-factories is imperative in cutting down the overall manufacturing cost. Thus, diverse metabolic engineering strategies and engineering tools have been established to obtain a preferred genotype and phenotype displaying superior productivity. However, these tools are limited to only a handful of genes with permanent modification of a genome and significant labor costs, and this is one of the bottlenecks associated with biofactory construction. Therefore, a groundbreaking rapid and high-throughput engineering tool is needed for efficient construction of microbial cell-factories. During the last decade, copious small noncoding RNAs (ncRNAs) have been discovered in bacteria. These are involved in substantial regulatory roles like transcriptional and post-transcriptional gene regulation by modulating mRNA elongation, stability, or translational efficiency. Because of their vulnerability, ncRNAs can be used as another layer of conditional control over gene expression without modifying chromosomal sequences, and hence would be a promising high-throughput tool for metabolic engineering. Here, we review successful design principles and applications of ncRNAs for high-throughput metabolic engineering or physiological studies of diverse industrially important microorganisms.

Identification of a Novel Small Non-Coding RNA Modulating the Intracellular Survival of Brucella melitensis

Bacterial small non-coding RNAs (sRNAs) are gene expression modulators respond to environmental changes, stressful conditions, and pathogenesis. In this study, by using a combined bioinformatic and experimental approach, eight novel sRNA genes were identified in intracellular pathogen Brucella melitensis. BSR0602, one sRNA that was highly induced in stationary phase, was further examined and found to modulate the intracellular survival of B. melitensis. BSR0602 was present at very high levels in vitro under stresses similar to those encountered during infection in host macrophages. Furthermore, BSR0602 was found to be highly expressed in the spleens of infected mice, suggesting its potential role in the control of pathogenesis. BSR0602 targets the mRNAs coding for gntR, a global transcriptional regulator, which is required for B. melitensis virulence. Overexpression of BSR0602 results in distinct reduction in the gntR mRNA level. B. melitensis with high level of BSR0602 is defective in bacteria intracellular survival in macrophages and defective in growth in the spleens of infected mice. Therefore, BSR0602 may directly inhibit the expression of gntR, which then impairs Brucellae intracellular survival and contributes to Brucella infection. Our findings suggest that BSR0602 is responsible for bacterial adaptation to stress conditions and thus modulate B. melitensis intracellular survival.

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.

Small transcriptome analysis indicates that the enzyme RppH influences both the quality and quantity of sRNAs in Neisseria gonorrhoeae

Prokaryotic mRNA turnover can be initiated by the removal of pyrophosphate from the 5' end of a transcript using the RNA pyrophosphohydrolase enzyme RppH. Following the initial dephosphorylation step, RNaseE then degrades the message into small oligonucleotide segments. This study assessed the small RNA transcriptome of Neisseria gonorrhoeae strain MS11 in two genetic backgrounds; using wild type cells as well as cells carrying a rppH insertional mutation. It was found that the presence of the RppH enzyme affected both the quantity and length of small RNAs (sRNAs) in various chromosomal locations and involved sense transcripts (seRNAs), transcripts originating from the opposite strand (asRNAs) as well as inter-genic-derived RNAs (IGRs). In comparing the two transcriptomes, we found that not all small RNAs were expressed in both genetic backgrounds, suggesting that RppH apparently targets only a subset of transcripts. Overall, this study shows that small RNAs can be detected from the majority of genes within the chromosome, as well as from inter-genic regions, and that more sRNA transcripts are detected in the absence of the RppH enzyme.

Cell wall stress activates expression of a novel stress response facilitator (SrfA) under σ22 (AlgT/U) control in Pseudomonas aeruginosa

The ECF (extracytoplasmic function) alternative sigma factor, σ(22) (AlgT/U), is required for expression of the algD promoter of the operon for alginate biosynthesis in Pseudomonas aeruginosa. Alginate production promotes chronic pulmonary infections by this opportunistic pathogen in patients with cystic fibrosis and chronic obstructive pulmonary disease. σ(22) is normally sequestered, but its deregulation and activation occur either by mutation in mucA (encoding an anti-sigma factor) or in response to envelope stress, such as inhibition of peptidoglycan synthesis. The σ(22) stress response system includes many genes in addition to those for alginate. In the present study, we characterized an intergenic region between ORFs PA2559 and PA2560 in PAO1 for a σ(22)-dependent, stress-responsive transcript, described here as PA2559.1. Northern analysis and transcript end-mapping indicated the PA2559.1 transcript was ~310 nt in length. Examination of the DNA sequence upstream of +1 revealed a σ(22) core promoter motif, GAATTT-N16-TCTGT, and site-directed mutagenesis confirmed this to be a σ(22)-dependent promoter that was highly activated during cell wall stress. PA2559.1 also contained an ORF that demonstrated increased translational activity upon cell wall stress. As determined by mutant analysis, the protein encoded by PA2559.1 was shown to play a positive role in the σ(22)-dependent activation of the algD promoter under stress in both sessile (i.e. biofilm) and planktonic conditions. Thus, it appeared to act as a stress response facilitator and so was named SrfA. The sequence of SrfA was found to be novel in nature and extremely well conserved only in P. aeruginosa, suggesting that it is of high evolutionary importance in this species.

Discovery of ethanol-responsive small RNAs in Zymomonas mobilis

Zymomonas mobilis is a bacterium that can produce ethanol by fermentation. Due to its unique metabolism and efficient ethanol production, Z. mobilis has attracted special interest for biofuel energy applications; an important area of study is the regulation of those specific metabolic pathways. Small RNAs (sRNAs) have been studied as molecules that function as transcriptional regulators in response to cellular stresses. While sRNAs have been discovered in various organisms by computational prediction and experimental approaches, their discovery in Z. mobilis has not yet been reported. In this study, we have applied transcriptome analysis and computational predictions to facilitate identification and validation of 15 novel sRNAs in Z. mobilis. We furthermore characterized their expression in the context of high and low levels of intracellular ethanol. Here, we report that 3 of the sRNAs (Zms2, Zms4, and Zms6) are differentially expressed under aerobic and anaerobic conditions, when low and high ethanol productions are observed, respectively. Importantly, when we tested the effect of ethanol stress on the expression of sRNAs in Z. mobilis, Zms2, Zms6, and Zms18 showed differential expression under 5% ethanol stress conditions. These data suggest that in this organism regulatory RNAs can be associated with metabolic functions involved in ethanol stress responses.

Bacillus subtilis Systems Biology: Applications of -Omics Techniques to the Study of Endospore Formation

Endospore-forming bacteria, with Bacillus subtilis being the prevalent model organism, belong to the phylum Firmicutes. Although the last common ancestor of all Firmicutes is likely to have been an endospore-forming species, not every lineage in the phylum has maintained the ability to produce endospores (hereafter, spores). In 1997, the release of the full genome sequence for B. subtilis strain 168 marked the beginning of the genomic era for the study of spore formation (sporulation). In this original genome sequence, 139 of the 4,100 protein-coding genes were annotated as sporulation genes. By the time a revised genome sequence with updated annotations was published in 2009, that number had increased significantly, especially since transcriptional profiling studies (transcriptomics) led to the identification of several genes expressed under the control of known sporulation transcription factors. Over the past decade, genome sequences for multiple spore-forming species have been released (including several strains in the Bacillus anthracis / Bacillus cereus group and many Clostridium species), and phylogenomic analyses have revealed many conserved sporulation genes. Parallel advances in transcriptomics led to the identification of small untranslated regulatory RNAs (sRNAs), including some that are expressed during sporulation. An extended array of -omics techniques, i.e., techniques designed to probe gene function on a genome-wide scale, such as proteomics, metabolomics, and high-throughput protein localization studies, have been implemented in microbiology. Combined with the use of new computational methods for predicting gene function and inferring regulatory relationships on a global scale, these -omics approaches are uncovering novel information about sporulation and a variety of other bacterial cell processes.