Biocomputational prediction of small non-coding RNAs in Streptomyces

Biocomputational Identification of sRNAs in Leptospira interrogans Serovar Lai

Regulatory small RNAs (sRNA) are RNA transcripts that are not translated into proteins but act as functional RNAs. Pathogenic Leptospira cause an epidemic spirochaetal zoonosis, Leptospirosis. It is speculated that Leptospiral sRNAs are involved in orchestrating their pathogenicity. In this study, biocomputational approach was adopted to identify Leptospiral sRNAs. In this study, two sRNA prediction programs, i.e., RNAz and nocoRNAc, were employed to screen the reference genome of Leptospira interrogans serovar Lai. Out of 126 predicted sRNAs, there are 96 cis-antisense sRNAs, 28 trans-encoded sRNAs and 2 sRNAs that partially overlap with protein-coding genes in a sense orientation. To determine whether these candidates are expressed in the pathogen, they were compared with the coverage files generated from our RNA-seq datasets. It was found out that 7 predicted sRNAs are expressed in mid-log phase, stationary phase, serum stress, temperature stress and iron stress while 2 sRNAs are expressed in mid-log phase, stationary phase, serum stress, and temperature stress. Besides, their expressions were also confirmed experimentally via RT-PCR. These experimentally validated candidates were also subjected to mRNA target prediction using TargetRNA2. Taken together, our study demonstrated that biocomputational strategy can serve as an alternative or as a complementary strategy to the laborious and expensive deep sequencing methods not only to uncover putative sRNAs but also to predict their targets in bacteria. In fact, this is the first study that integrates computational approach to predict putative sRNAs in L. interrogans serovar Lai.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-022-01050-9.

Comparative genome analysis of mycobacteria focusing on tRNA and non-coding RNA

BACKGROUND

The Mycobacterium genus encompasses at least 192 named species, many of which cause severe diseases such as tuberculosis. Non-tuberculosis mycobacteria (NTM) can also infect humans and animals. Some are of emerging concern because they show high resistance to commonly used antibiotics while others are used and evaluated in bioremediation or included in anticancer vaccines.

RESULTS

We provide the genome sequences for 114 mycobacterial type strains and together with 130 available mycobacterial genomes we generated a phylogenetic tree based on 387 core genes and supported by average nucleotide identity (ANI) data. The 244 genome sequences cover most of the species constituting the Mycobacterium genus. The genome sizes ranged from 3.2 to 8.1 Mb with an average of 5.7 Mb, and we identified 14 new plasmids. Moreover, mycobacterial genomes consisted of phage-like sequences ranging between 0 and 4.64% dependent on mycobacteria while the number of IS elements varied between 1 and 290. Our data also revealed that, depending on the mycobacteria, the number of tRNA and non-coding (nc) RNA genes differ and that their positions on the chromosome varied. We identified a conserved core set of 12 ncRNAs, 43 tRNAs and 18 aminoacyl-tRNA synthetases among mycobacteria.

CONCLUSIONS

Phages, IS elements, tRNA and ncRNAs appear to have contributed to the evolution of the Mycobacterium genus where several tRNA and ncRNA genes have been horizontally transferred. On the basis of our phylogenetic analysis, we identified several isolates of unnamed species as new mycobacterial species or strains of known mycobacteria. The predicted number of coding sequences correlates with genome size while the number of tRNA, rRNA and ncRNA genes does not. Together these findings expand our insight into the evolution of the Mycobacterium genus and as such they establish a platform to understand mycobacterial pathogenicity, their evolution, antibiotic resistance/tolerance as well as the function and evolution of ncRNA among mycobacteria.

An extended catalogue of ncRNAs in Streptomyces coelicolor reporting abundant tmRNA, RNase-P RNA and RNA fragments derived from pre-ribosomal RNA leader sequences

Streptomyces coelicolor is a model organism for studying streptomycetes. This genus possesses relevant medical and economical roles, because it produces many biologically active metabolites of pharmaceutical interest, including the majority of commercialized antibiotics. In this bioinformatic study, the transcriptome of S. coelicolor has been analyzed to identify novel RNA species and quantify the expression of both annotated and novel transcripts in solid and liquid growth medium cultures at different times. The major characteristics disclosed in this study are: (i) the diffuse antisense transcription; (ii) the great abundance of transfer-messenger RNAs (tmRNA); (iii) the abundance of rnpB transcripts, paramount for the RNase-P complex; and (iv) the presence of abundant fragments derived from pre-ribosomal RNA leader sequences of unknown biological function. Overall, this study extends the catalogue of ncRNAs in S. coelicolor and suggests an important role of non-coding transcription in the regulation of biologically active molecule production.

Ms1 RNA Interacts With the RNA Polymerase Core in Streptomyces coelicolor and Was Identified in Majority of Actinobacteria Using a Linguistic Gene Synteny Search

Bacteria employ small non-coding RNAs (sRNAs) to regulate gene expression. Ms1 is an sRNA that binds to the RNA polymerase (RNAP) core and affects the intracellular level of this essential enzyme. Ms1 is structurally related to 6S RNA that binds to a different form of RNAP, the holoenzyme bearing the primary sigma factor. 6S RNAs are widespread in the bacterial kingdom except for the industrially and medicinally important Actinobacteria. While Ms1 RNA was identified in Mycobacterium, it is not clear whether Ms1 RNA is present also in other Actinobacteria species. Here, using a computational search based on secondary structure similarities combined with a linguistic gene synteny approach, we identified Ms1 RNA in Streptomyces. In S. coelicolor, Ms1 RNA overlaps with the previously annotated scr3559 sRNA with an unknown function. We experimentally confirmed that Ms1 RNA/scr3559 associates with the RNAP core without the primary sigma factor HrdB in vivo. Subsequently, we applied the computational approach to other Actinobacteria and identified Ms1 RNA candidates in 824 Actinobacteria species, revealing Ms1 RNA as a widespread class of RNAP binding sRNAs, and demonstrating the ability of our multifactorial computational approach to identify weakly conserved sRNAs in evolutionarily distant genomes.

Interplay between non-coding RNA transcription, stringent phenotype and antibiotic production in Streptomyces

While in recent years the key role of non-coding RNAs (ncRNAs) in regulation of gene expression has become increasingly evident, their interaction with the global regulatory circuits is still obscure. Here we analyzed the structure and organization of the transcriptome of Streptomyces ambofaciens, the producer of spiramycin. We identified ncRNAs including 45 small-RNAs (sRNAs) and 119 antisense-RNAs (asRNAs I) that appear transcribed from dedicated promoters. Some sRNAs and asRNAs are unprecedented in Streptomyces, and were predicted to target mRNAs encoding proteins involved in transcription, translation, ribosomal structure and biogenesis, and regulation of morphological and biochemical differentiation. We then compared ncRNA expression in three strains: i.) the wild type strain; ii.) an isogenic pirA-defective mutant with central carbon metabolism imbalance, "relaxed" phenotype, and repression of antibiotic production; iii.) a pirA-derivative strain harboring a "stringent" RNA polymerase that suppresses pirA-associated phenotypes. Data indicated that expression of most ncRNAs was correlated to the stringent/relaxed phenotype suggesting novel effector mechanisms of the stringent response.

Genome-scale analysis of genetic regulatory elements in Streptomyces avermitilis MA-4680 using transcript boundary information

Background

The gram-positive bacterium, Streptomyces avermitilis, holds industrial importance as the producer of avermectin, a widely used anthelmintic agent, and a heterologous expression host of secondary metabolite-biosynthetic gene clusters. Despite its industrial importance, S. avermitilis’ genome organization and regulation of gene expression remain poorly understood. In this study, four different types of Next-Generation Sequencing techniques, including dRNA-Seq, Term-Seq, RNA-Seq and ribosome profiling, were applied to S. avermitilis to determine transcription units of S. avermitilis at a genome-wide level and elucidate regulatory elements for transcriptional and translational control of individual transcription units.

Result

By applying dRNA-Seq and Term-Seq to S. avermitilis MA-4680, a total of 2361 transcription start sites and 2017 transcript 3′-end positions were identified, respectively, leading to determination of 1601 transcription units encoded in S. avermitilis’ genome. Cataloguing the transcription units and integrated analysis of multiple high-throughput data types revealed the presence of diverse regulatory elements for gene expression, such as promoters, 5′-UTRs, terminators, 3′-UTRs and riboswitches. The conserved promoter motifs were identified from 2361 transcription start sites as 5′-TANNNT and 5′-BTGACN for the − 10 and − 35 elements, respectively. The − 35 element and spacer lengths between − 10 and − 35 elements were critical for transcriptional regulation of functionally distinct genes, suggesting the involvement of unique sigma factors. In addition, regulatory sequences recognized by antibiotic regulatory proteins were identified from the transcription start site information. Analysis of the 3′-end of RNA transcript revealed that stem structure formation is a major determinant for transcription termination of most transcription units.

Conclusions

The transcription unit architecture elucidated from the transcripts’ boundary information provides insights for unique genetic regulatory mechanisms of S. avermitilis. Our findings will elevate S. avermitilis’ potential as a production host for a diverse set of secondary metabolites.

6S-Like scr3559 RNA Affects Development and Antibiotic Production in Streptomyces coelicolor

Regulatory RNAs control a number of physiological processes in bacterial cells. Here we report on a 6S-like RNA transcript (scr3559) that affects both development and antibiotic production in Streptomyces coelicolor. Its expression is enhanced during the transition to stationary phase. Strains that over-expressed the scr3559 gene region exhibited a shortened exponential growth phase in comparison with a control strain; accelerated aerial mycelium formation and spore maturation; alongside an elevated production of actinorhodin and undecylprodigiosin. These observations were supported by LC-MS analyses of other produced metabolites, including: germicidins, desferrioxamines, and coelimycin. A subsequent microarray differential analysis revealed increased expression of genes associated with the described morphological and physiological changes. Structural and functional similarities between the scr3559 transcript and 6S RNA, and its possible employment in regulating secondary metabolite production are discussed.

Interplay between Non-Coding RNA Transcription, Stringent/Relaxed Phenotype and Antibiotic Production in Streptomyces ambofaciens

While in recent years the key role of non-coding RNAs (ncRNAs) in the regulation of gene expression has become increasingly evident, their interaction with the global regulatory circuits is still obscure. Here we analyzed the structure and organization of the transcriptome of Streptomyces ambofaciens, the producer of spiramycin. We identified ncRNAs including 45 small-RNAs (sRNAs) and 119 antisense-RNAs (asRNAs I) that appear transcribed from dedicated promoters. Some sRNAs and asRNAs are unprecedented in Streptomyces and were predicted to target mRNAs encoding proteins involved in transcription, translation, ribosomal structure and biogenesis, and regulation of morphological and biochemical differentiation. We then compared ncRNA expression in three strains: (i) the wild-type strain; (ii) an isogenic pirA-defective mutant with central carbon metabolism imbalance, “relaxed” phenotype, and repression of antibiotic production; and (iii) a pirA-derivative strain harboring a “stringent” RNA polymerase that suppresses pirA-associated phenotypes. Data indicated that the expression of most ncRNAs was correlated to the stringent/relaxed phenotype suggesting novel effector mechanisms of the stringent response.

sRNA scr5239 Involved in Feedback Loop Regulation of Streptomyces coelicolor Central Metabolism

In contrast to transcriptional regulation, post-transcriptional regulation and the role of small non-coding RNAs (sRNAs) in streptomycetes are not well studied. Here, we focus on the highly conserved sRNA scr5239 in Streptomyces coelicolor. A proteomics approach revealed that the sRNA regulates several metabolic enzymes, among them phosphoenolpyruvate carboxykinase (PEPCK), a key enzyme of the central carbon metabolism. The sRNA scr5239 represses pepck at the post-transcriptional level and thus modulates the intracellular level of phosphoenolpyruvate (PEP). The expression of scr5239 in turn is dependent on the global transcriptional regulator DasR, thus creating a feedback loop regulation of the central carbon metabolism. By post-transcriptional regulation of PEPCK and in all likelihood other targets, scr5239 adds an additional layer to the DasR regulatory network and provides a tool to control the metabolism dependent on the available carbon source.

Mycobacterium tuberculosis 6C sRNA binds multiple mRNA targets via C-rich loops independent of RNA chaperones

Bacterial small regulatory RNAs (sRNAs) are the most abundant class of post-transcriptional regulators and have been well studied in Gram-negative bacteria. Little is known about the functions and mechanisms of sRNAs in high GC Gram-positive bacteria including Mycobacterium and Streptomyces. Here, we performed an in-depth study of 6C sRNA of Mycobacterium tuberculosis, which is conserved among high GC Gram-positive bacteria. Forty-seven genes were identified as possible direct targets of 6C sRNA and 15 of them were validated using an in vivo translational lacZ fusion system. We found that 6C sRNA plays a pleotropic role and regulates genes involved in various cellular processes, including DNA replication and protein secretion. Mapping the interactions of 6C sRNA with mRNA targets panD and dnaB revealed that the C-rich loops of 6C sRNA act as direct binding sites to mRNA targets. Unlike in Gram-negative bacteria where RNA binding proteins Hfq and ProQ are required, the interactions of 6C sRNA with mRNAs appear to be independent of RNA chaperones. Our findings suggest that the multiple G–C pairings between single stranded regions are sufficient to establish stable interactions between 6C sRNA and mRNA targets, providing a mechanism for sRNAs in high GC Gram-positive bacteria.

On a stake-out: Mycobacterial small RNA identification and regulation

Persistence of mycobacteria in the hostile environment of human macrophage is pivotal for its successful pathogenesis. Rapid adaptation to diverse stresses is the key aspect for their survival in the host cells. A range of heterogeneous mechanisms operate in bacteria to retaliate stress conditions. Small RNAs (sRNA) have been implicated in many of those mechanisms in either a single or multiple regulatory networks to post-transcriptionally modulate bacterial gene expression. Although small RNA profiling in mycobacteria by advanced technologies like deep sequencing, tilling microarray etc. have identified hundreds of sRNA, however, a handful of those small RNAs have been unearthed with precise regulatory mechanism. Extensive investigations on sRNA-mediated gene regulations in eubacteria like Escherichia coli revealed the existence of a plethora of distinctive sRNA mechanisms e.g. base pairing, protein sequestration, RNA decoy etc. Increasing studies on mycobacterial sRNA also discovered several eccentric mechanisms where sRNAs act at the posttranscriptional stage to either activate or repress target gene expression that lead to promote mycobacterial survival in stresses. Several intrinsic features like high GC content, absence of any homologue of abundant RNA chaperones, Hfq and ProQ, isolate sRNA mechanisms of mycobacteria from that of other bacteria. An insightful approach has been taken in this review to describe sRNA identification and its regulations in mycobacterial species especially in Mycobacterium tuberculosis.

DNA mapping and kinetic modeling of the HrdB regulon in Streptomyces coelicolor

HrdB in streptomycetes is a principal sigma factor whose deletion is lethal. This is also the reason why its regulon has not been investigated so far. To overcome experimental obstacles, for investigating the HrdB regulon, we constructed a strain whose HrdB protein was tagged by an HA epitope. ChIP-seq experiment, done in 3 repeats, identified 2137 protein-coding genes organized in 337 operons, 75 small RNAs, 62 tRNAs, 6 rRNAs and 3 miscellaneous RNAs. Subsequent kinetic modeling of regulation of protein-coding genes with HrdB alone and with a complex of HrdB and a transcriptional cofactor RbpA, using gene expression time series, identified 1694 genes that were under their direct control. When using the HrdB-RbpA complex in the model, an increase of the model fidelity was found for 322 genes. Functional analysis revealed that HrdB controls the majority of gene groups essential for the primary metabolism and the vegetative growth. Particularly, almost all ribosomal protein-coding genes were found in the HrdB regulon. Analysis of promoter binding sites revealed binding motif at the -10 region and suggested the possible role of mono- or di-nucleotides upstream of the -10 element.

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.

Comparative genomics of Mycobacterium mucogenicum and Mycobacterium neoaurum clade members emphasizing tRNA and non-coding RNA

Background

Mycobacteria occupy various ecological niches and can be isolated from soil, tap water and ground water. Several cause diseases in humans and animals. To get deeper insight into our understanding of mycobacterial evolution focusing on tRNA and non-coding (nc)RNA, we conducted a comparative genome analysis of Mycobacterium mucogenicum (Mmuc) and Mycobacterium neoaurum (Mneo) clade members.

Results

Genome sizes for Mmuc- and Mneo-clade members vary between 5.4 and 6.5 Mbps with the complete Mmuc^T (type strain) genome encompassing 6.1 Mbp. The number of tRNA genes range between 46 and 79 (including one pseudo tRNA gene) with 39 tRNA genes common among the members of these clades, while additional tRNA genes were probably acquired through horizontal gene transfer. Selected tRNAs and ncRNAs (RNase P RNA, tmRNA, 4.5S RNA, Ms1 RNA and 6C RNA) are expressed, and the levels for several of these are higher in stationary phase compared to exponentially growing cells. The rare tRNA^IleTAT isoacceptor and two for mycobacteria novel ncRNAs: the Lactobacillales-derived GOLLD RNA and a homolog to the antisense Salmonella typhimurium phage Sar RNA, were shown to be present and expressed in certain Mmuc-clade members.

Conclusions

Phages, IS elements, horizontally transferred tRNA gene clusters, and phage-derived ncRNAs appears to have influenced the evolution of the Mmuc- and Mneo-clades. While the number of predicted coding sequences correlates with genome size, the number of tRNA coding genes does not. The majority of the tRNA genes in mycobacteria are transcribed mainly from single genes and the levels of certain ncRNAs, including RNase P RNA (essential for the processing of tRNAs), are higher at stationary phase compared to exponentially growing cells. We provide supporting evidence that Ms1 RNA represents a mycobacterial 6S RNA variant. The evolutionary routes for the ncRNAs RNase P RNA, tmRNA and Ms1 RNA are different from that of the core genes.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1447-7) contains supplementary material, which is available to authorized users.

Small RNAs of Haloferax mediterranei: Identification and Potential Involvement in Nitrogen Metabolism

Small RNAs have been studied in detail in domains Bacteria and Eukarya but, in the case of the domain Archaea, the knowledge is scarce and the physiological function of these small RNAs (sRNAs) is still uncertain. To extend the knowledge of sRNAs in the domain Archaea and their possible role in the regulation of the nitrogen assimilation metabolism in haloarchaea, Haloferax mediterranei has been used as a model microorganism. The bioinformatic approach has allowed for the prediction of 295 putative sRNAs genes in the genome of H. mediterranei, 88 of which have been verified by means of RNA-Sequencing (RNA-Seq). The secondary structure of these sRNAs and their possible targets have been identified. Curiously, some of them present as possible target genes relating to nitrogen assimilation, such as glutamate dehydrogenase and the nitrogen regulatory PII protein. Analysis of RNA-Seq data has also revealed differences in the expression pattern of 16 sRNAs according to the nitrogen source. Consequently, RNomic and bioinformatic approaches used in this work have allowed for the identification of new sRNAs in H. mediterranei, some of which show different expression patterns depending on the nitrogen source. This suggests that these sRNAs could be involved in the regulation of nitrogen assimilation and can constitute an important gene regulatory network.

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

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

Genome-Wide Discovery of Putative sRNAs in Paracoccus denitrificans Expressed under Nitrous Oxide Emitting Conditions

Nitrous oxide (N₂O) is a stable, ozone depleting greenhouse gas. Emissions of N₂O into the atmosphere continue to rise, primarily due to the use of nitrogen-containing fertilizers by soil denitrifying microbes. It is clear more effective mitigation strategies are required to reduce emissions. One way to help develop future mitigation strategies is to address the currently poor understanding of transcriptional regulation of the enzymes used to produce and consume N₂O. With this ultimate aim in mind we performed RNA-seq on a model soil denitrifier, Paracoccus denitrificans, cultured anaerobically under high N₂O and low N₂O emitting conditions, and aerobically under zero N₂O emitting conditions to identify small RNAs (sRNAs) with potential regulatory functions transcribed under these conditions. sRNAs are short (∼40–500 nucleotides) non-coding RNAs that regulate a wide range of activities in many bacteria. Hundred and sixty seven sRNAs were identified throughout the P. denitrificans genome which are either present in intergenic regions or located antisense to ORFs. Furthermore, many of these sRNAs are differentially expressed under high N₂O and low N₂O emitting conditions respectively, suggesting they may play a role in production or reduction of N₂O. Expression of 16 of these sRNAs have been confirmed by RT-PCR. Ninety percent of the sRNAs are predicted to form secondary structures. Predicted targets include transporters and a number of transcriptional regulators. A number of sRNAs were conserved in other members of the α-proteobacteria. Better understanding of the sRNA factors which contribute to expression of the machinery required to reduce N₂O will, in turn, help to inform strategies for mitigation of N₂O emissions.

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.

A comparison of key aspects of gene regulation in Streptomyces coelicolor and Escherichia coli using nucleotide-resolution transcription maps produced in parallel by global and differential RNA sequencing

Streptomyces coelicolor is a model for studying bacteria renowned as the foremost source of natural products used clinically. Post-genomic studies have revealed complex patterns of gene expression and links to growth, morphological development and individual genes. However, the underlying regulation remains largely obscure, but undoubtedly involves steps after transcription initiation. Here we identify sites involved in RNA processing and degradation as well as transcription within a nucleotide-resolution map of the transcriptional landscape. This was achieved by combining RNA-sequencing approaches suited to the analysis of GC-rich organisms. Escherichia coli was analysed in parallel to validate the methodology and allow comparison. Previously, sites of RNA processing and degradation had not been mapped on a transcriptome-wide scale for E. coli. Through examples, we show the value of our approach and data sets. This includes the identification of new layers of transcriptional complexity associated with several key regulators of secondary metabolism and morphological development in S. coelicolor and the identification of host-encoded leaderless mRNA and rRNA processing associated with the generation of specialized ribosomes in E. coli. New regulatory small RNAs were identified for both organisms. Overall the results illustrate the diversity in mechanisms used by different bacterial groups to facilitate and regulate gene expression.

Identification of sRNAs expressed by the human pathogen Neisseria gonorrhoeae under disparate growth conditions

In the last several years, bacterial gene regulation via small RNAs (sRNAs) has been recognized as an important mechanism controlling expression of essential proteins that are critical to bacterial growth and metabolism. Technologies such as RNA-seq are rapidly expanding the field of sRNAs and are enabling a global view of the “sRNAome” of several bacterial species. While numerous sRNAs have been identified in a variety of both Gram-negative and Gram-positive bacteria, only a very small number have been fully characterized in the human pathogen Neisseria gonorrhoeae, the etiological agent of the STD gonorrhea. Here we present the first analysis of N. gonorrhoeae specifically focused on the identification of sRNAs through RNA-seq analysis of the organism cultured under different in vitro growth conditions. Using a new computational program, Rockhopper, to analyze prokaryotic RNA-seq data obtained from N. gonorrhoeae we identified several putative sRNAs and confirmed their expression and size through Northern blot analysis. In addition, RNA was collected from four different growth conditions (iron replete and deplete, as well as with and without co-culture with human endocervical cells). Many of the putative sRNAs identified shoed varying expression levels relative to the different growth conditions examine or were detected only under certain conditions but not others. Comparisons of identified sRNAs with the regulatory pattern of putative mRNA targets revealed possible functional roles for these sRNAs. These studies are the first to carry out a global analysis of N. gonorrhoeae specifically focused on sRNAs and show that RNA-mediated regulation may be an important mechanism of gene control in this human pathogen.

Noncoding RNA in Mycobacteria

Efforts to understand the molecular basis of mycobacterial gene regulation are dominated by a protein-centric view. However, there is a growing appreciation that noncoding RNA, i.e., RNA that is not translated, plays a role in a wide variety of molecular mechanisms. Noncoding RNA comprises rRNA, tRNA, 4.5S RNA, RnpB, and transfer-messenger RNA, as well as a vast population of regulatory RNA, often dubbed "the dark matter of gene regulation." The regulatory RNA species comprise 5' and 3' untranslated regions and a rapidly expanding category of transcripts with the ability to base-pair with mRNAs or to interact with proteins. Regulatory RNA plays a central role in the bacterium's response to changes in the environment, and in this article we review emerging information on the presence and abundance of different types of noncoding RNA in mycobacteria.

Small non-coding RNAs in streptomycetes

Streptomycetes are Gram-positive, GC-rich, soil dwelling bacteria, occurring ubiquitary throughout nature. They undergo extensive morphological changes from spores to filamentous mycelia and produce a plethora of secondary metabolites. Owing to their complex life cycle, streptomycetes require efficient regulatory machinery for the control of gene expression. Therefore, they possess a large diversity of regulators. Within this review we summarize the current knowledge about the importance of small non-coding RNA for the control of gene expression in these organisms.

Complex intra-operonic dynamics mediated by a small RNA in Streptomyces coelicolor

Streptomyces are predominantly soil-dwelling bacteria that are best known for their multicellular life cycle and their prodigious metabolic capabilities. They are also renowned for their regulatory capacity and flexibility, with each species encoding >60 sigma factors, a multitude of transcription factors, and an increasing number of small regulatory RNAs. Here, we describe our characterization of a conserved small RNA (sRNA), scr4677. In the model species Streptomyces coelicolor, this sRNA is located in the intergenic region separating SCO4677 (an anti-sigma factor-encoding gene) and SCO4676 (a putative regulatory protein-encoding gene), close to the SCO4676 translation start site in an antisense orientation. There appears to be considerable genetic interplay between these different gene products, with wild type expression of scr4677 requiring function of the anti-sigma factor SCO4677, and scr4677 in turn influencing the abundance of SCO4676-associated transcripts. The scr4677-mediated effects were independent of RNase III (a double stranded RNA-specific nuclease), with RNase III having an unexpectedly positive influence on the level of SCO4676-associated transcripts. We have shown that both SCO4676 and SCO4677 affect the production of the blue-pigmented antibiotic actinorhodin under specific growth conditions, and that this activity appears to be independent of scr4677.

Prediction and characterization of small non-coding RNAs related to secondary metabolites in Saccharopolyspora erythraea

Saccharopolyspora erythraea produces a large number of secondary metabolites with biological activities, including erythromycin. Elucidation of the mechanisms through which the production of these secondary metabolites is regulated may help to identify new strategies for improved biosynthesis of erythromycin. In this paper, we describe the systematic prediction and analysis of small non-coding RNAs (sRNAs) in S. erythraea, with the aim to elucidate sRNA-mediated regulation of secondary metabolite biosynthesis. In silico and deep-sequencing technologies were applied to predict sRNAs in S. erythraea. Six hundred and forty-seven potential sRNA loci were identified, of which 382 cis-encoded antisense RNA are complementary to protein-coding regions and 265 predicted transcripts are located in intergenic regions. Six candidate sRNAs (sernc292, sernc293, sernc350, sernc351, sernc361, and sernc389) belong to four gene clusters (tpc3, pke, pks6, and nrps5) that are involved in secondary metabolite biosynthesis. Deep-sequencing data showed that the expression of all sRNAs in the strain HL3168 E3 (E3) was higher than that in NRRL23338 (M), except for sernc292 and sernc361 expression. The relative expression of six sRNAs in strain M and E3 were validated by qRT-PCR at three different time points (24, 48, and 72 h). The results showed that, at each time point, the transcription levels of sernc293, sernc350, sernc351, and sernc389 were higher in E3 than in M, with the largest difference observed at 72 h, whereas no signals for sernc292 and sernc361 were detected. sernc293, sernc350, sernc351, and sernc389 probably regulate iron transport, terpene metabolism, geosmin synthesis, and polyketide biosynthesis, respectively. The major significance of this study is the successful prediction and identification of sRNAs in genomic regions close to the secondary metabolism-related genes in S. erythraea. A better understanding of the sRNA-target interaction would help to elucidate the complete range of functions of sRNAs in S. erythraea, including sRNA-mediated regulation of erythromycin biosynthesis.

Inverse regulation of Fe- and Ni-containing SOD genes by a Fur family regulator Nur through small RNA processed from 3'UTR of the sodF mRNA

Superoxide dismutases (SODs) are widely distributed enzymes that convert superoxides to hydrogen peroxide and molecular oxygen, using various metals as cofactors. Many actinobacteria contain genes for both Ni-containing (sodN) and Fe-containing (sodF) SODs. In Streptomyces coelicolor, expression of the sodF and sodN genes is inversely regulated by nickel-specific Nur, a Fur-family regulator. With sufficient nickel, Nur directly represses sodF transcription, while inducing sodN indirectly. Bioinformatic search revealed that a conserved 19-nt stretch upstream of sodN matches perfectly with the sodF downstream sequence. We found that the sodF gene produced a stable small-sized RNA species (s-SodF) of ∼ 90 nt that harbors the anti-sodN sequence complementary to sodN mRNA from the 5'-end up to the ribosome binding site. Absence of nearby promoters and sensitivity to 5'-phosphate-specific exonuclease indicated that the s-SodF RNA is a likely processed product of sodF mRNA. The s-SodF RNA caused a significant decrease in the half-life of the sodN mRNA. Therefore, Nur activates sodN expression through inhibiting the synthesis of sodF mRNA, from which inhibitory s-SodF RNA is generated. This reveals a novel mechanism by which antagonistic regulation of one gene is achieved by small RNA processed from the 3'UTR of another gene's mRNA.

Inference of sigma factor controlled networks by using numerical modeling applied to microarray time series data of the germinating prokaryote

A computational model of gene expression was applied to a novel test set of microarray time series measurements to reveal regulatory interactions between transcriptional regulators represented by 45 sigma factors and the genes expressed during germination of a prokaryote Streptomyces coelicolor. Using microarrays, the first 5.5 h of the process was recorded in 13 time points, which provided a database of gene expression time series on genome-wide scale. The computational modeling of the kinetic relations between the sigma factors, individual genes and genes clustered according to the similarity of their expression kinetics identified kinetically plausible sigma factor-controlled networks. Using genome sequence annotations, functional groups of genes that were predominantly controlled by specific sigma factors were identified. Using external binding data complementing the modeling approach, specific genes involved in the control of the studied process were identified and their function suggested.

SdrA, a new DeoR family regulator involved in Streptomyces avermitilis morphological development and antibiotic production

The SAV3339 (SdrA) protein of Streptomyces avermitilis, a member of the DeoR family of regulators, was assessed to determine its in vivo function by gene knockdown through the use of cis-encoded noncoding RNA and knockout of the sdrA gene. These analyses revealed that SdrA represents another class of Streptomyces regulator that controls morphological development and antibiotic production.

Comparative analysis of non-coding RNAs in the antibiotic-producing Streptomyces bacteria

Background

Non-coding RNAs (ncRNAs) are key regulatory elements that control a wide range of cellular processes in all bacteria in which they have been studied. Taking advantage of recent technological innovations, we set out to fully explore the ncRNA potential of the multicellular, antibiotic-producing Streptomyces bacteria.

Results

Using a comparative RNA sequencing analysis of three divergent model streptomycetes (S. coelicolor, S. avermitilis and S. venezuelae), we discovered hundreds of novel cis-antisense RNAs and intergenic small RNAs (sRNAs). We identified a ubiquitous antisense RNA species that arose from the overlapping transcription of convergently-oriented genes; we termed these RNA species ‘cutoRNAs’, for convergent untranslated overlapping RNAs. Conservation between different classes of ncRNAs varied greatly, with sRNAs being more conserved than antisense RNAs. Many species-specific ncRNAs, including many distinct cutoRNA pairs, were located within antibiotic biosynthetic clusters, including the actinorhodin, undecylprodigiosin, and coelimycin clusters of S. coelicolor, the chloramphenicol cluster of S. venezuelae, and the avermectin cluster of S. avermitilis.

Conclusions

These findings indicate that ncRNAs, including a novel class of antisense RNA, may exert a previously unrecognized level of regulatory control over antibiotic production in these bacteria. Collectively, this work has dramatically expanded the ncRNA repertoire of three Streptomyces species and has established a critical foundation from which to investigate ncRNA function in this medically and industrially important bacterial genus.

Synthetic RNA Silencing of Actinorhodin Biosynthesis in Streptomyces coelicolor A3(2)

We demonstrate the first application of synthetic RNA gene silencers in Streptomyces coelicolor A3(2). Peptide nucleic acid and expressed antisense RNA silencers successfully inhibited actinorhodin production. Synthetic RNA silencing was target-specific and is a new tool for gene regulation and metabolic engineering studies in Streptomyces.

Computational small RNA prediction in bacteria

Bacterial, small RNAs were once regarded as potent regulators of gene expression and are now being considered as essential for their diversified roles. Many small RNAs are now reported to have a wide array of regulatory functions, ranging from environmental sensing to pathogenesis. Traditionally, noncoding transcripts were rarely detected by means of genetic screens. However, the availability of approximately 2200 prokaryotic genome sequences in public databases facilitates the efficient computational search of those molecules, followed by experimental validation. In principle, the following four major computational methods were applied for the prediction of sRNA locations from bacterial genome sequences: (1) comparative genomics, (2) secondary structure and thermodynamic stability, (3) ‘Orphan’ transcriptional signals and (4) ab initio methods regardless of sequence or structure similarity; most of these tools were applied to locate the putative genomic sRNA locations followed by experimental validation of those transcripts. Therefore, computational screening has simplified the sRNA identification process in bacteria. In this review, a plethora of small RNA prediction methods and tools that have been reported in the past decade are discussed comprehensively and assessed based on their attributes, compatibility, and their prediction accuracy.

Molecular regulation of antibiotic biosynthesis in streptomyces

Streptomycetes are the most abundant source of antibiotics. Typically, each species produces several antibiotics, with the profile being species specific. Streptomyces coelicolor, the model species, produces at least five different antibiotics. We review the regulation of antibiotic biosynthesis in S. coelicolor and other, nonmodel streptomycetes in the light of recent studies. The biosynthesis of each antibiotic is specified by a large gene cluster, usually including regulatory genes (cluster-situated regulators [CSRs]). These are the main point of connection with a plethora of generally conserved regulatory systems that monitor the organism's physiology, developmental state, population density, and environment to determine the onset and level of production of each antibiotic. Some CSRs may also be sensitive to the levels of different kinds of ligands, including products of the pathway itself, products of other antibiotic pathways in the same organism, and specialized regulatory small molecules such as gamma-butyrolactones. These interactions can result in self-reinforcing feed-forward circuitry and complex cross talk between pathways. The physiological signals and regulatory mechanisms may be of practical importance for the activation of the many cryptic secondary metabolic gene cluster pathways revealed by recent sequencing of numerous Streptomyces genomes.

Computational discovery and RT-PCR validation of novel Burkholderia conserved and Burkholderia pseudomallei unique sRNAs

BACKGROUND

The sRNAs of bacterial pathogens are known to be involved in various cellular roles including environmental adaptation as well as regulation of virulence and pathogenicity. It is expected that sRNAs may also have similar functions for Burkholderia pseudomallei, a soil bacterium that can adapt to diverse environmental conditions, which causes the disease melioidosis and is also able to infect a wide variety of hosts.

RESULTS

By integrating several proven sRNA prediction programs into a computational pipeline, available Burkholderia spp. genomes were screened to identify sRNA gene candidates. Orthologous sRNA candidates were then identified via comparative analysis. From the total prediction, 21 candidates were found to have Rfam homologs. RT-PCR and sequencing of candidate sRNA genes of unknown functions revealed six putative sRNAs which were highly conserved in Burkholderia spp. and two that were unique to B. pseudomallei present in a normal culture conditions transcriptome. The validated sRNAs include potential cis-acting elements associated with the modulation of methionine metabolism and one B. pseudomallei-specific sRNA that is expected to bind to the Hfq protein.

CONCLUSIONS

The use of the pipeline developed in this study and subsequent comparative analysis have successfully aided in the discovery and shortlisting of sRNA gene candidates for validation. This integrated approach identified 29 B. pseudomallei sRNA genes - of which 21 have Rfam homologs and 8 are novel.

Streptomyces coelicolor sRNA scr5239 inhibits agarase expression by direct base pairing to the dagA coding region

Transcriptional regulation of primary and secondary metabolism is well-studied in Streptomyces coelicolor, a model organism for antibiotic production and cell differentiation. In contrast, little is known about post-transcriptional regulation and the potential functions of small non-coding RNAs (sRNAs) in this Gram-positive, GC-rich soil bacterium. Here, we report the identification and characterization of scr5239, an sRNA highly conserved in the genus Streptomyces. The sRNA is 159 nt long, composed of five stem-loops, and encoded in the intergenic region between SCO5238 and SCO5239. scr5239 expression is constitutive under several stress and growth conditions but dependent on the nitrogen supply. scr5239 decreases the production of the antibiotic actinorhodin, and represses expression of the extracellular agarase dagA at the post-transcriptional level by direct base pairing to the coding region 33 nt downstream of the ribosome-binding site.

Convergent transcription in the butyrolactone regulon in Streptomyces coelicolor confers a bistable genetic switch for antibiotic biosynthesis

cis-encoded antisense RNAs (cis asRNA) have been reported to participate in gene expression regulation in both eukaryotic and prokaryotic organisms. Its presence in Streptomyces coelicolor has also been reported recently; however, its role has yet to be fully investigated. Using mathematical modeling we explore the role of cis asRNA produced as a result of convergent transcription in scbA-scbR genetic switch. scbA and scbR gene pair, encoding repressor–amplifier proteins respectively, mediates the synthesis of a signaling molecule, the γ-butyrolactone SCB1 and controls the onset of antibiotic production. Our model considers that transcriptional interference caused by convergent transcription of two opposing RNA polymerases results in fatal collision and transcriptional termination, which suppresses transcription efficiency. Additionally, convergent transcription causes sense and antisense interactions between complementary sequences from opposing strands, rendering the full length transcript inaccessible for translation. We evaluated the role of transcriptional interference and the antisense effect conferred by convergent transcription on the behavior of scbA-scbR system. Stability analysis showed that while transcriptional interference affects the system, it is asRNA that confers scbA-scbR system the characteristics of a bistable switch in response to the signaling molecule SCB1. With its critical role of regulating the onset of antibiotic synthesis the bistable behavior offers this two gene system the needed robustness to be a genetic switch. The convergent two gene system with potential of transcriptional interference is a frequent feature in various genomes. The possibility of asRNA regulation in other such gene-pairs is yet to be examined.

Deep sequencing-based identification of small non-coding RNAs in Streptomyces coelicolor

Streptomyces coelicolor is considered the model organism among Gram positive, GC rich bacteria. Its genome has been sequenced but little is known about the occurrence and distribution of small non-coding RNAs in this biotechnologically relevant organism. Using deep sequencing we analyzed the transcriptome at the end of exponential growth, which corresponds to the onset of secondary metabolism. We mapped 193 transcriptional start sites of mRNA genes and identified putative new and alternative open reading frames. We identified 63 non-coding RNAs including 29 cis encoded antisense RNAs, and confirmed expression for 11, most of them being growth-phase dependent. A comparison between the sequencing results and bioinformatic sRNA predictions using Dynalign and RNAz revealed only a small overlap between the different approaches.

RNIE: genome-wide prediction of bacterial intrinsic terminators

Bacterial Rho-independent terminators (RITs) are important genomic landmarks involved in gene regulation and terminating gene expression. In this investigation we present RNIE, a probabilistic approach for predicting RITs. The method is based upon covariance models which have been known for many years to be the most accurate computational tools for predicting homology in structural non-coding RNAs. We show that RNIE has superior performance in model species from a spectrum of bacterial phyla. Further analysis of species where a low number of RITs were predicted revealed a highly conserved structural sequence motif enriched near the genic termini of the pathogenic Actinobacteria, Mycobacterium tuberculosis. This motif, together with classical RITs, account for up to 90% of all the significantly structured regions from the termini of M. tuberculosis genic elements. The software, predictions and alignments described below are available from http://github.com/ppgardne/RNIE.

nocoRNAc: characterization of non-coding RNAs in prokaryotes

Background

The interest in non-coding RNAs (ncRNAs) constantly rose during the past few years because of the wide spectrum of biological processes in which they are involved. This led to the discovery of numerous ncRNA genes across many species. However, for most organisms the non-coding transcriptome still remains unexplored to a great extent. Various experimental techniques for the identification of ncRNA transcripts are available, but as these methods are costly and time-consuming, there is a need for computational methods that allow the detection of functional RNAs in complete genomes in order to suggest elements for further experiments. Several programs for the genome-wide prediction of functional RNAs have been developed but most of them predict a genomic locus with no indication whether the element is transcribed or not.

Results

We present NOCORNAc, a program for the genome-wide prediction of ncRNA transcripts in bacteria. NOCORNAc incorporates various procedures for the detection of transcriptional features which are then integrated with functional ncRNA loci to determine the transcript coordinates. We applied RNAz and NOCORNAc to the genome of Streptomyces coelicolor and detected more than 800 putative ncRNA transcripts most of them located antisense to protein-coding regions. Using a custom design microarray we profiled the expression of about 400 of these elements and found more than 300 to be transcribed, 38 of them are predicted novel ncRNA genes in intergenic regions. The expression patterns of many ncRNAs are similarly complex as those of the protein-coding genes, in particular many antisense ncRNAs show a high expression correlation with their protein-coding partner.

Conclusions

We have developed NOCORNAc, a framework that facilitates the automated characterization of functional ncRNAs. NOCORNAc increases the confidence of predicted ncRNA loci, especially if they contain transcribed ncRNAs. NOCORNAc is not restricted to intergenic regions, but it is applicable to the prediction of ncRNA transcripts in whole microbial genomes. The software as well as a user guide and example data is available at http://www.zbit.uni-tuebingen.de/pas/nocornac.htm.

Small RNAs in the genus Clostridium

The genus Clostridium includes major human pathogens and species important to cellulose degradation, the carbon cycle, and biotechnology. Small RNAs (sRNAs) are emerging as crucial regulatory molecules in all organisms, but they have not been investigated in clostridia. Research on sRNAs in clostridia is hindered by the absence of a systematic method to identify sRNA candidates, thus delegating clostridial sRNA research to a hit-and-miss process. Thus, we wanted to develop a method to identify potential sRNAs in the Clostridium genus to open up the field of sRNA research in clostridia. Using comparative genomics analyses combined with predictions of rho-independent terminators and promoters, we predicted sRNAs in 21 clostridial genomes: Clostridium acetobutylicum, C. beijerinckii, C. botulinum (eight strains), C. cellulolyticum, C. difficile, C. kluyveri (two strains), C. novyi, C. perfringens (three strains), C. phytofermentans, C. tetani, and C. thermocellum. Although more than one-third of predicted sRNAs have Shine-Dalgarno (SD) sequences, only one-sixth have a start codon downstream of SD sequences; thus, most of the predicted sRNAs are noncoding RNAs. Quantitative reverse transcription-PCR (Q-RT-PCR) and Northern analysis were employed to test the presence of a randomly chosen set of sRNAs in C. acetobutylicum and several C. botulinum strains, leading to the confirmation of a large fraction of the tested sRNAs. We identified a conserved, novel sRNA which, together with the downstream gene coding for an ATP-binding cassette (ABC) transporter gene, responds to the antibiotic clindamycin. The number of predicted sRNAs correlated with the physiological function of the species (high for pathogens, low for cellulolytic, and intermediate for solventogenic), but not with 16S rRNA-based phylogeny.

The suboptimal structures find the optimal RNAs: homology search for bacterial non-coding RNAs using suboptimal RNA structures

Non-coding RNAs (ncRNAs) are regulatory molecules encoded in the intergenic or intragenic regions of the genome. In prokaryotes, biocomputational identification of homologs of known ncRNAs in other species often fails due to weakly evolutionarily conserved sequences, structures, synteny and genome localization, except in the case of evolutionarily closely related species. To eliminate results from weak conservation, we focused on RNA structure, which is the most conserved ncRNA property. Analysis of the structure of one of the few well-studied bacterial ncRNAs, 6S RNA, demonstrated that unlike optimal and consensus structures, suboptimal structures are capable of capturing RNA homology even in divergent bacterial species. A computational procedure for the identification of homologous ncRNAs using suboptimal structures was created. The suggested procedure was applied to strongly divergent bacterial species and was capable of identifying homologous ncRNAs.

Experimental identification and characterization of 97 novel npcRNA candidates in Salmonella enterica serovar Typhi

We experimentally identified and characterized 97 novel, non-protein-coding RNA candidates (npcRNAs) from the human pathogen Salmonella enterica serovar Typhi (hereafter referred to as S. typhi). Three were specific to S. typhi, 22 were restricted to Salmonella species and 33 were differentially expressed during S. typhi growth. We also identified Salmonella Pathogenicity Island-derived npcRNAs that might be involved in regulatory mechanisms of virulence, antibiotic resistance and pathogenic specificity of S. typhi. An in-depth characterization of S. typhi StyR-3 npcRNA showed that it specifically interacts with RamR, the transcriptional repressor of the ramA gene, which is involved in the multidrug resistance (MDR) of Salmonella. StyR-3 interfered with RamR–DNA binding activity and thus potentially plays a role in regulating ramA gene expression, resulting in the MDR phenotype. Our study also revealed a large number of cis-encoded antisense npcRNA candidates, supporting previous observations of global sense–antisense regulatory networks in bacteria. Finally, at least six of the npcRNA candidates interacted with the S. typhi Hfq protein, supporting an important role of Hfq in npcRNA networks. This study points to novel functional npcRNA candidates potentially involved in various regulatory roles including the pathogenicity of S. typhi.

A broadening world of bacterial small RNAs

The ubiquity of small RNAs (sRNAs) in bacteria is now well established. These transcripts are the members of regulatory circuits involved in diverse processes ranging from stress adaptation to virulence to metabolism. Recent large-scale searches suggest that there exist many times more sRNAs than previously predicted even in the best studied bacterial transcriptomes. On the basis of these and other recent findings of regulatory sRNAs that do not function in a 'classical' manner, we propose that the working definition of sRNAs be broadened.

Noncoding RNA of glutamine synthetase I modulates antibiotic production in Streptomyces coelicolor A3(2)

Overexpression of antisense chromosomal cis-encoded noncoding RNAss (ncRNAs) in glutamine synthetase I resulted in a decrease in growth, protein synthesis, and antibiotic production in Streptomyces coelicolor. In addition, we predicted 3,597 cis-encoded ncRNAs and validated 13 of them experimentally, including several ncRNAs that are differentially expressed in bacterial hormone-defective mutants.

Small RNA identification in Enterobacteriaceae using synteny and genomic backbone retention II

Small RNAs are bacterial counterparts of noncoding RNAs. Increasing evidence being added in the literature indicates that these small RNAs play major roles in prokaryotes both at the transcriptome and proteome levels. Based on comparative genomic studies, we present manually curated small RNA regions in 25 recently completed genomes from Enterobacteriaceae. The study is a continuation of our earlier work that uses the presence of small RNAs sandwiched between specific conserved flanking genes retaining genomic backbone and gene synteny. Based on this study, a total of 931 identified sRNA/sRNA regions are reported. This data contains 498 small RNA homologs, 80 putative small RNA regions containing partial stretches of homologous sequences, and 353 putative nonhomologous sRNA regions. This homologs/partial homologs includes, 84 putative small RNA homologous regions retaining at least one of the conserved flanking genes pair which may possibly act as hotspots for genetic pool insertion/deletion in genomes. Nonhomologous CsrB sRNA region reported by us in Yersinia pseudotuberculosis IP32953 has been experimentally confirmed by Kulkarni's group and sraH and ryeE sRNAs from Erwinia carotovora subsp. atroseptica SCRI1043 recently added to the Rfam database are indicative proof of our positive approach.

Identification of small RNAs in Mycobacterium tuberculosis

In spite of being one of our most prominent bacterial pathogens, the presence of small regulatory RNAs (sRNAs) has not previously been investigated in Mycobacterium tuberculosis. Post-transcriptional regulation of gene expression by sRNA molecules has been demonstrated in a wide range of pathogenic bacteria and has been shown to play a significant role in the control of virulence. By screening cDNA libraries prepared from low-molecular weight RNA from M. tuberculosis we have identified nine putative sRNA molecules, including cis-encoded antisense transcripts from within open reading frames and trans-encoded transcripts from intergenic regions. sRNAs displayed differential expression between exponential and stationary phase, and during a variety of stress conditions. Two of the cis-encoded sRNAs were associated with genes encoding enzymes involved in lipid metabolism, desA1 and pks12. These sRNAs showed complementarity to multiple M. tuberculosis genes, suggesting the potential to act as both cis-encoded and trans-encoded sRNAs. Overexpression of selected trans-encoded sRNAs had profound impact on growth of M. tuberculosis and M. smegmatis. This is the first experimental evidence of sRNAs in M. tuberculosis and it will be important to consider the potential influence of sRNA regulation when studying the transcriptome and the proteome of M. tuberculosis during infection.