Biocomputational prediction of non-coding RNAs in model cyanobacteria

Small antisense RNA ThfR positively regulates Thf1 in Synechocystis sp. PCC 6803

Thylakoid formation1 (Thf1), encoded by sll1414 (thf1), is a multifunctional protein conserved in all photosynthetic organisms. thf1 expression is highly induced by high light in Synechocystis during photosynthesis-related stress. In this study, differential RNA sequencing analysis of the Synechocystis sp. PCC 6803 revealed a small antisense RNA (asRNA) gene located on the reverse-complementary strand of the thf1 gene. The full length of this asRNA (designated ThfR) was determined by 5' and 3' RACE analysis. The accumulation of thf1 mRNA was up-regulated synchronously with the ThfR level during survival after high-light stress or nitrogen starvation. Under nitrogen starvation or high-light stress, compared with the wild type, a ThfR overexpression mutant demonstrated relatively more Thf1 protein content, while a ThfR reduced-expression mutant accumulated less Thf1 protein. Furthermore, the overexpression of ThfR enhanced the electron transport rate and the proliferation of cyanobacteria under high-light stress. These results, which we confirmed further using an Escherichia coli sRNA expression platform, suggest that the thf1 gene is positively regulated by ThfR, possibly through protection of the RAUUW element at the RNase E cleavage site. This study represents the first report, to our knowledge, of a cis-transcript antisense RNA that targets thf1 in Synechocystis sp. PCC 6803 and provides evidence that ThfR regulates photosynthesis by positively modulating thf1 under high-light conditions.

Genetic Responses of Metabolically Active Limnospira indica Strain PCC 8005 Exposed to γ-Radiation during Its Lifecycle

Two morphotypes of the cyanobacterial Limnospira indica (formerly Arthrospira sp.) strain PCC 8005, denoted as P2 (straight trichomes) and P6 (helical trichomes), were subjected to chronic gamma radiation from spent nuclear fuel (SNF) rods at a dose rate of ca. 80 Gy·h-1 for one mass doubling period (approximately 3 days) under continuous light with photoautotrophic metabolism fully active. Samples were taken for post-irradiation growth recovery and RNA-Seq transcriptional analysis at time intervals of 15, 40, and 71.5 h corresponding to cumulative doses of ca. 1450, 3200, and 5700 Gy, respectively. Both morphotypes, which were previously reported by us to display different antioxidant capacities and differ at the genomic level in 168 SNPs, 48 indels and 4 large insertions, recovered equally well from 1450 and 3200 Gy. However, while the P2 straight type recovered from 5700 Gy by regaining normal growth within 6 days, the P6 helical type took about 13 days to recover from this dose, indicating differences in their radiation tolerance and response. To investigate these differences, P2 and P6 cells exposed to the intermediate dose of gamma radiation (3200 Gy) were analyzed for differential gene expression by RNA-Seq analysis. Prior to batch normalization, a total of 1553 genes (887 and 666 of P2 and P6, respectively, with 352 genes in common) were selected based on a two-fold change in expression and a false discovery rate FDR smaller or equal to 0.05. About 85% of these 1553 genes encoded products of yet unknown function. Of the 229 remaining genes, 171 had a defined function while 58 genes were transcribed into non-coding RNA including 21 tRNAs (all downregulated). Batch normalization resulted in 660 differentially expressed genes with 98 having a function and 32 encoding RNA. From PCC 8005-P2 and PCC 8005-P6 expression patterns, it emerges that although the cellular routes used by the two substrains to cope with ionizing radiation do overlap to a large extent, both strains displayed a distinct preference of priorities.

Discovery of 17 conserved structural RNAs in fungi

Many non-coding RNAs with known functions are structurally conserved: their intramolecular secondary and tertiary interactions are maintained across evolutionary time. Consequently, the presence of conserved structure in multiple sequence alignments can be used to identify candidate functional non-coding RNAs. Here, we present a bioinformatics method that couples iterative homology search with covariation analysis to assess whether a genomic region has evidence of conserved RNA structure. We used this method to examine all unannotated regions of five well-studied fungal genomes (Saccharomyces cerevisiae, Candida albicans, Neurospora crassa, Aspergillus fumigatus, and Schizosaccharomyces pombe). We identified 17 novel structurally conserved non-coding RNA candidates, which include four H/ACA box small nucleolar RNAs, four intergenic RNAs and nine RNA structures located within the introns and untranslated regions (UTRs) of mRNAs. For the two structures in the 3' UTRs of the metabolic genes GLY1 and MET13, we performed experiments that provide evidence against them being eukaryotic riboswitches.

Inverse regulation of light harvesting and photoprotection is mediated by a 3'-end-derived sRNA in cyanobacteria

Phycobilisomes (PBSs), the principal cyanobacterial antenna, are among the most efficient macromolecular structures in nature, and are used for both light harvesting and directed energy transfer to the photosynthetic reaction center. However, under unfavorable conditions, excess excitation energy needs to be rapidly dissipated to avoid photodamage. The orange carotenoid protein (OCP) senses light intensity and induces thermal energy dissipation under stress conditions. Hence, its expression must be tightly controlled; however, the molecular mechanism of this regulation remains to be elucidated. Here, we describe the discovery of a posttranscriptional regulatory mechanism in Synechocystis sp. PCC 6803 in which the expression of the operon encoding the allophycocyanin subunits of the PBS is directly and in an inverse fashion linked to the expression of OCP. This regulation is mediated by ApcZ, a small regulatory RNA that is derived from the 3'-end of the tetracistronic apcABC-apcZ operon. ApcZ inhibits ocp translation under stress-free conditions. Under most stress conditions, apc operon transcription decreases and ocp translation increases. Thus, a key operon involved in the collection of light energy is functionally connected to the expression of a protein involved in energy dissipation. Our findings support the view that regulatory RNA networks in bacteria evolve through the functionalization of mRNA 3'-UTRs.

Identification of small RNAs involved in nitrogen fixation in Anabaena sp. PCC 7120 based on RNA-seq under steady state conditions

Purpose

Anabaena sp. PCC7120 is a genetically tractable model organism for nitrogen fixation and photosynthesis research. The importance of small regulatory RNAs (sRNAs) as mediators of a number of cellular processes in bacteria has begun to be recognized. Bacterial sRNA binds to target genes through base pairing, and play a regulatory role. Many studies have shown that bacterial sRNA can regulate cell stress response, carbon and nitrogen fixation, and so on. However, little is known about sRNAs in Anabaena sp. PCC 7120 regarded to nitrogen fixation under later steady state.

Methods

To provide a comprehensive study of sRNAs in this model organism, the sRNA (< 200 nt) extracted from Anabaena sp. PCC 7120 under nitrogen step-down treatment of 12 days, together with the sRNA from the control, was analyzed using deep RNA sequencing. Possible target genes regulated by all identified putative sRNAs were predicted by IntaRNA and further analyzed for functional categorizations for biological pathways.

Result

Totally, 14,132 transcripts were produced from the de novo assembly. Among them, transcripts that are located either in the intergenic region or antisense strand were kept, which resulted in 1219 sRNA candidates, for further analysis. RPKM-based differential expression analysis showed that 418 sRNAs were significantly differentially expressed between the samples from control (nitrogen addition, N+) and nitrogen depletion, (N−). Among them, 303 sRNAs were significantly upregulated, whereas 115 sRNAs were significantly downregulated. RT-PCR of 18 randomly chosen sRNAs showed a similar pattern as RNA-seq result, which confirmed the reliability of the RNA-seq data. In addition, the possible target genes regulated by unique sRNAs of Anabaena sp. PCC 7120 under nitrogen addition (N+) condition or that under nitrogen depletion (N−) condition were analyzed for functional categorization and biological pathways, which provided the evidences that sRNAs were indeed involved in many different metabolic pathways.

Conclusion

The information from the present study provides a valuable reference for understanding the sRNA-mediated regulation of the nitrogen fixation in Anabaena PCC 7120 under steady state conditions.

Regulatory RNA at the crossroads of carbon and nitrogen metabolism in photosynthetic cyanobacteria

Cyanobacteria are photosynthetic bacteria that populate widely different habitats. Accordingly, cyanobacteria exhibit a wide spectrum of lifestyles, physiologies, and morphologies and possess genome sizes and gene numbers which may vary by up to a factor of ten within the phylum. Consequently, large differences exist between individual species in the size and complexity of their regulatory networks. Several non-coding RNAs have been identified that play crucial roles in the acclimation responses of cyanobacteria to changes in the environment. Some of these regulatory RNAs are conserved throughout the cyanobacterial phylum, while others exist only in a few taxa. Here we give an overview on characterized regulatory RNAs in cyanobacteria, with a focus on regulators of photosynthesis, carbon and nitrogen metabolism. However, chances are high that these regulators represent just the tip of the iceberg.

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.

Inactivation of the RNA helicase CrhR impacts a specific subset of the transcriptome in the cyanobacterium Synechocystis sp. PCC 6803

DEAD-box RNA-helicases catalyze the reorganization of structured RNAs and the formation of RNP complexes. The cyanobacterium Synechocystis sp. PCC 6803 encodes a single DEAD-box RNA helicase, CrhR (Slr0083), whose expression is regulated by abiotic stresses that alter the redox potential of the photosynthetic electron transport chain, including temperature downshift. Despite its proposed effect on RNA metabolism and its known relevance in cold-stress adaptation, the reported impact of a CrhR knockout on the cold adaption of the transcriptome only identified eight affected genes. Here, we utilized a custom designed microarray to assess the impact of the absence of CrhR RNA helicase activity on the transcriptome, independent of cold stress. CrhR truncation impacts an RNA subset comprising ~10% of the ncRNA and also ~10% of the mRNA transcripts. While equal numbers of mRNAs showed increased as well as decreased abundance, more than 90% of the ncRNAs showed enhanced expression in the absence of CrhR, indicative of a negative effect on ncRNA transcription or stability. We further tested the effect of CrhR on the stability of strongly responding RNAs that identify examples of post-transcriptional and transcriptional regulation. The data suggest that CrhR impacts multiple aspects of RNA metabolism in Synechocystis.

The Distinctive Regulation of Cyanobacterial Glutamine Synthetase

Glutamine synthetase (GS) features prominently in bacterial nitrogen assimilation as it catalyzes the entry of bioavailable nitrogen in form of ammonium into cellular metabolism. The classic example, the comprehensively characterized GS of enterobacteria, is subject to exquisite regulation at multiple levels, among them gene expression regulation to control GS abundance, as well as feedback inhibition and covalent modifications to control enzyme activity. Intriguingly, the GS of the ecologically important clade of cyanobacteria features fundamentally different regulatory systems to those of most prokaryotes. These include the interaction with small proteins, the so-called inactivating factors (IFs) that inhibit GS linearly with their abundance. In addition to this protein interaction-based regulation of GS activity, cyanobacteria use alternative elements to control the synthesis of GS and IFs at the transcriptional level. Moreover, cyanobacteria evolved unique RNA-based regulatory mechanisms such as glutamine riboswitches to tightly tune IF abundance. In this review, we aim to outline the current knowledge on the distinctive features of the cyanobacterial GS encompassing the overall control of its activity, sensing the nitrogen status, transcriptional and post-transcriptional regulation, as well as strain-specific differences.

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.

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.

Systematic and functional identification of small non-coding RNAs associated with exogenous biofuel stress in cyanobacterium Synechocystis sp. PCC 6803

BACKGROUND

The unicellular model cyanobacterium Synechocystis sp. PCC 6803 is considered a promising microbial chassis for biofuel production. However, its low tolerance to biofuel toxicity limits its potential application. Although recent studies showed that bacterial small RNAs (sRNAs) play important roles in regulating cellular processes in response to various stresses, the role of sRNAs in resisting exogenous biofuels is yet to be determined.

RESULTS

Based on genome-wide sRNA sequencing combined with systematic analysis of previous transcriptomic and proteomic data under the same biofuel or environmental perturbations, we report the identification of 133 trans-encoded sRNA transcripts with high-resolution mapping of sRNAs in Synechocystis, including 23 novel sRNAs identified for the first time. In addition, according to quantitative expression analysis and sRNA regulatory network prediction, sRNAs potentially involved in biofuel tolerance were identified and functionally confirmed by constructing sRNA overexpression or suppression strains of Synechocystis. Notably, overexpression of sRNA Nc117 revealed an improved tolerance to ethanol and butanol, while suppression of Nc117 led to increased sensitivity.

CONCLUSIONS

The study provided the first comprehensive responses to exogenous biofuels at the sRNA level in Synechocystis and opens an avenue to engineering sRNA regulatory elements for improved biofuel tolerance in the cyanobacterium Synechocystis.

Inactivation of iron-sulfur cluster biogenesis regulator SufR in Synechocystis sp. PCC 6803 induces unique iron-dependent protein-level responses

BACKGROUND

Iron-sulfur (Fe-S) clusters are protein-bound cofactors associated with cellular electron transport and redox sensing, with multiple specific functions in oxygen-evolving photosynthetic cyanobacteria. The aim here was to elucidate protein-level effects of the transcriptional repressor SufR involved in the regulation of Fe-S cluster biogenesis in the cyanobacterium Synechocystis sp. PCC 6803.

METHODS

The approach was to quantitate 94 pre-selected target proteins associated with various metabolic functions using SRM in Synechocystis. The evaluation was conducted in response to sufR deletion under different iron conditions, and complemented with EPR analysis on the functionality of the photosystems I and II as well as with RT-qPCR to verify the effects of SufR also on transcript level.

RESULTS

The results on both protein and transcript levels show that SufR acts not only as a repressor of the suf operon when iron is available but also has other direct and indirect functions in the cell, including maintenance of the expression of pyruvate:ferredoxin oxidoreductase NifJ and other Fe-S cluster proteins under iron sufficient conditions. Furthermore, the results imply that in the absence of iron the suf operon is repressed by some additional regulatory mechanism independent of SufR.

CONCLUSIONS

The study demonstrates that Fe-S cluster metabolism in Synechocystis is stringently regulated, and has complex interactions with multiple primary functions in the cell, including photosynthesis and central carbon metabolism.

GENERAL SIGNIFICANCE

The study provides new insight into the regulation of Fe-S cluster biogenesis via suf operon, and the associated wide-ranging protein-level changes in photosynthetic cyanobacteria.

Small Antisense RNA RblR Positively Regulates RuBisCo in Synechocystis sp. PCC 6803

Small regulatory RNAs (sRNAs) function as transcriptional and post-transcriptional regulators of gene expression in organisms from all domains of life. Cyanobacteria are thought to have developed a complex RNA-based regulatory mechanism. In the current study, by genome-wide analysis of differentially expressed small RNAs in Synechocystis sp. PCC 6803 under high light conditions, we discovered an asRNA (RblR) that is 113nt in length and completely complementary to its target gene rbcL, which encodes the large chain of RuBisCO, the enzyme that catalyzes carbon fixation. Further analysis of the RblR(+)/(−) mutants revealed that RblR acts as a positive regulator of rbcL under various stress conditions; Suppressing RblR adversely affects carbon assimilation and thus the yield, and those phenotypes of both the wild type and the overexpressor could be downgraded to the suppressor level by carbonate depletion, indicated a regulatory role of RblR in CO₂ assimilation. In addition, a real-time expression platform in Escherichia coli was setup and which confirmed that RblR promoted the translation of the rbcL mRNA into the RbcL protein. The present study is the first report of a regulatory RNA that targets RbcL in Synechocystis sp. PCC 6803, and provides strong evidence that RblR regulates photosynthesis by positively modulating rbcL expression in Synechocystis.

The Complex Transcriptional Response of Acaryochloris marina to Different Oxygen Levels

Ancient oxygenic photosynthetic prokaryotes produced oxygen as a waste product, but existed for a long time under an oxygen-free (anoxic) atmosphere, before an oxic atmosphere emerged. The change in oxygen levels in the atmosphere influenced the chemistry and structure of many enzymes that contained prosthetic groups that were inactivated by oxygen. In the genome of Acaryochloris marina, multiple gene copies exist for proteins that are normally encoded by a single gene copy in other cyanobacteria. Using high throughput RNA sequencing to profile transcriptome responses from cells grown under microoxic and hyperoxic conditions, we detected 8446 transcripts out of the 8462 annotated genes in the Cyanobase database. Two-thirds of the 50 most abundant transcripts are key proteins in photosynthesis. Microoxic conditions negatively affected the levels of expression of genes encoding photosynthetic complexes, with the exception of some subunits. In addition to the known regulation of the multiple copies of psbA, we detected a similar transcriptional pattern for psbJ and psbU, which might play a key role in the altered components of photosystem II. Furthermore, regulation of genes encoding proteins important for reactive oxygen species-scavenging is discussed at genome level, including, for the first time, specific small RNAs having possible regulatory roles under varying oxygen levels.

Small proteins in cyanobacteria provide a paradigm for the functional analysis of the bacterial micro-proteome

Background

Despite their versatile functions in multimeric protein complexes, in the modification of enzymatic activities, intercellular communication or regulatory processes, proteins shorter than 80 amino acids (μ-proteins) are a systematically underestimated class of gene products in bacteria. Photosynthetic cyanobacteria provide a paradigm for small protein functions due to extensive work on the photosynthetic apparatus that led to the functional characterization of 19 small proteins of less than 50 amino acids. In analogy, previously unstudied small ORFs with similar degrees of conservation might encode small proteins of high relevance also in other functional contexts.

Results

Here we used comparative transcriptomic information available for two model cyanobacteria, Synechocystis sp. PCC 6803 and Synechocystis sp. PCC 6714 for the prediction of small ORFs. We found 293 transcriptional units containing candidate small ORFs ≤80 codons in Synechocystis sp. PCC 6803, also including the known mRNAs encoding small proteins of the photosynthetic apparatus. From these transcriptional units, 146 are shared between the two strains, 42 are shared with the higher plant Arabidopsis thaliana and 25 with E. coli. To verify the existence of the respective μ-proteins in vivo, we selected five genes as examples to which a FLAG tag sequence was added and re-introduced them into Synechocystis sp. PCC 6803. These were the previously annotated gene ssr1169, two newly defined genes norf1 and norf4, as well as nsiR6(nitrogen stress-induced RNA 6) and hliR1(high light-inducible RNA 1) , which originally were considered non-coding. Upon activation of expression via the Cu^2+.responsive petE promoter or from the native promoters, all five proteins were detected in Western blot experiments.

Conclusions

The distribution and conservation of these five genes as well as their regulation of expression and the physico-chemical properties of the encoded proteins underline the likely great bandwidth of small protein functions in bacteria and makes them attractive candidates for functional studies.

Experimental Identification of Small Non-Coding RNAs in the Model Marine Bacterium Ruegeria pomeroyi DSS-3

In oligotrophic ocean waters where bacteria are often subjected to chronic nutrient limitation, community transcriptome sequencing has pointed to the presence of highly abundant small RNAs (sRNAs). The role of sRNAs in regulating response to nutrient stress was investigated in a model heterotrophic marine bacterium Ruegeria pomeroyi grown in continuous culture under carbon (C) and nitrogen (N) limitation. RNAseq analysis identified 99 putative sRNAs. Sixty-nine were cis-encoded and located antisense to a presumed target gene. Thirty were trans-encoded and initial target prediction was performed computationally. The most prevalent functional roles of genes anti-sense to the cis-sRNAs were transport, cell-cell interactions, signal transduction, and transcriptional regulation. Most sRNAs were transcribed equally under both C and N limitation, and may be involved in a general stress response. However, 14 were regulated differentially between the C and N treatments and may respond to specific nutrient limitations. A network analysis of the predicted target genes of the R. pomeroyi cis-sRNAs indicated that they average fewer connections than typical protein-encoding genes, and appear to be more important in peripheral or niche-defining functions encoded in the pan genome.

Identification of Conserved and Potentially Regulatory Small RNAs in Heterocystous Cyanobacteria

Small RNAs (sRNAs) are a growing class of non-protein-coding transcripts that participate in the regulation of virtually every aspect of bacterial physiology. Heterocystous cyanobacteria are a group of photosynthetic organisms that exhibit multicellular behavior and developmental alternatives involving specific transcriptomes exclusive of a given physiological condition or even a cell type. In the context of our ongoing effort to understand developmental decisions in these organisms we have undertaken an approach to the global identification of sRNAs. Using differential RNA-Seq we have previously identified transcriptional start sites for the model heterocystous cyanobacterium Nostoc sp. PCC 7120. Here we combine this dataset with a prediction of Rho-independent transcriptional terminators and an analysis of phylogenetic conservation of potential sRNAs among 89 available cyanobacterial genomes. In contrast to predictive genome-wide approaches, the use of an experimental dataset comprising all active transcriptional start sites (differential RNA-Seq) facilitates the identification of bona fide sRNAs. The output of our approach is a dataset of predicted potential sRNAs in Nostoc sp. PCC 7120, with different degrees of phylogenetic conservation across the 89 cyanobacterial genomes analyzed. Previously described sRNAs appear among the predicted sRNAs, demonstrating the performance of the algorithm. In addition, new predicted sRNAs are now identified that can be involved in regulation of different aspects of cyanobacterial physiology, including adaptation to nitrogen stress, the condition that triggers differentiation of heterocysts (specialized nitrogen-fixing cells). Transcription of several predicted sRNAs that appear exclusively in the genomes of heterocystous cyanobacteria is experimentally verified by Northern blot. Cell-specific transcription of one of these sRNAs, NsiR8 (nitrogen stress-induced RNA 8), in developing heterocysts is also demonstrated.

The sRNA NsiR4 is involved in nitrogen assimilation control in cyanobacteria by targeting glutamine synthetase inactivating factor IF7

Glutamine synthetase (GS), a key enzyme in biological nitrogen assimilation, is regulated in multiple ways in response to varying nitrogen sources and levels. Here we show a small regulatory RNA, NsiR4 (nitrogen stress-induced RNA 4), which plays an important role in the regulation of GS in cyanobacteria. NsiR4 expression in the unicellular Synechocystis sp. PCC 6803 and in the filamentous, nitrogen-fixing Anabaena sp. PCC 7120 is stimulated through nitrogen limitation via NtcA, the global transcriptional regulator of genes involved in nitrogen metabolism. NsiR4 is widely conserved throughout the cyanobacterial phylum, suggesting a conserved function. In silico target prediction, transcriptome profiling on pulse overexpression, and site-directed mutagenesis experiments using a heterologous reporter system showed that NsiR4 interacts with the 5'UTR of gifA mRNA, which encodes glutamine synthetase inactivating factor (IF)7. In Synechocystis, we observed an inverse relationship between the levels of NsiR4 and the accumulation of IF7 in vivo. This NsiR4-dependent modulation of gifA (IF7) mRNA accumulation influenced the glutamine pool and thus [Formula: see text] assimilation via GS. As a second target, we identified ssr1528, a hitherto uncharacterized nitrogen-regulated gene. Competition experiments between WT and an ΔnsiR4 KO mutant showed that the lack of NsiR4 led to decreased acclimation capabilities of Synechocystis toward oscillating nitrogen levels. These results suggest a role for NsiR4 in the regulation of nitrogen metabolism in cyanobacteria, especially for the adaptation to rapid changes in available nitrogen sources and concentrations. NsiR4 is, to our knowledge, the first identified bacterial sRNA regulating the primary assimilation of a macronutrient.

Regulatory RNAs in photosynthetic cyanobacteria

Regulatory RNAs play versatile roles in bacteria in the coordination of gene expression during various physiological processes, especially during stress adaptation. Photosynthetic bacteria use sunlight as their major energy source. Therefore, they are particularly vulnerable to the damaging effects of excess light or UV irradiation. In addition, like all bacteria, photosynthetic bacteria must adapt to limiting nutrient concentrations and abiotic and biotic stress factors. Transcriptome analyses have identified hundreds of potential regulatory small RNAs (sRNAs) in model cyanobacteria such as Synechocystis sp. PCC 6803 or Anabaena sp. PCC 7120, and in environmentally relevant genera such as Trichodesmium, Synechococcus and Prochlorococcus. Some sRNAs have been shown to actually contain μORFs and encode short proteins. Examples include the 40-amino-acid product of the sml0013 gene, which encodes the NdhP subunit of the NDH1 complex. In contrast, the functional characterization of the non-coding sRNA PsrR1 revealed that the 131 nt long sRNA controls photosynthetic functions by targeting multiple mRNAs, providing a paradigm for sRNA functions in photosynthetic bacteria. We suggest that actuatons comprise a new class of genetic elements in which an sRNA gene is inserted upstream of a coding region to modify or enable transcription of that region.

Variations in the non-coding transcriptome as a driver of inter-strain divergence and physiological adaptation in bacteria

In all studied organisms, a substantial portion of the transcriptome consists of non-coding RNAs that frequently execute regulatory functions. Here, we have compared the primary transcriptomes of the cyanobacteria Synechocystis sp. PCC 6714 and PCC 6803 under 10 different conditions. These strains share 2854 protein-coding genes and a 16S rRNA identity of 99.4%, indicating their close relatedness. Conserved major transcriptional start sites (TSSs) give rise to non-coding transcripts within the sigB gene, from the 5′UTRs of cmpA and isiA, and 168 loci in antisense orientation. Distinct differences include single nucleotide polymorphisms rendering promoters inactive in one of the strains, e.g., for cmpR and for the asRNA PsbA2R. Based on the genome-wide mapped location, regulation and classification of TSSs, non-coding transcripts were identified as the most dynamic component of the transcriptome. We identified a class of mRNAs that originate by read-through from an sRNA that accumulates as a discrete and abundant transcript while also serving as the 5′UTR. Such an sRNA/mRNA structure, which we name ‘actuaton’, represents another way for bacteria to remodel their transcriptional network. Our findings support the hypothesis that variations in the non-coding transcriptome constitute a major evolutionary element of inter-strain divergence and capability for physiological adaptation.

Trichodesmium genome maintains abundant, widespread noncoding DNA in situ, despite oligotrophic lifestyle

Understanding the evolution of the free-living, cyanobacterial, diazotroph Trichodesmium is of great importance because of its critical role in oceanic biogeochemistry and primary production. Unlike the other >150 available genomes of free-living cyanobacteria, only 63.8% of the Trichodesmium erythraeum (strain IMS101) genome is predicted to encode protein, which is 20-25% less than the average for other cyanobacteria and nonpathogenic, free-living bacteria. We use distinctive isolates and metagenomic data to show that low coding density observed in IMS101 is a common feature of the Trichodesmium genus, both in culture and in situ. Transcriptome analysis indicates that 86% of the noncoding space is expressed, although the function of these transcripts is unclear. The density of noncoding, possible regulatory elements predicted in Trichodesmium, when normalized per intergenic kilobase, was comparable and twofold higher than that found in the gene-dense genomes of the sympatric cyanobacterial genera Synechococcus and Prochlorococcus, respectively. Conserved Trichodesmium noncoding RNA secondary structures were predicted between most culture and metagenomic sequences, lending support to the structural conservation. Conservation of these intergenic regions in spatiotemporally separated Trichodesmium populations suggests possible genus-wide selection for their maintenance. These large intergenic spacers may have developed during intervals of strong genetic drift caused by periodic blooms of a subset of genotypes, which may have reduced effective population size. Our data suggest that transposition of selfish DNA, low effective population size, and high-fidelity replication allowed the unusual "inflation" of noncoding sequence observed in Trichodesmium despite its oligotrophic lifestyle.

Exploring the secrets of long noncoding RNAs

High-throughput sequencing has revealed that the majority of RNAs have no capacity to encode protein. Among these non-coding transcripts, recent work has focused on the roles of long noncoding RNAs (lncRNAs) of >200 nucleotides. Although many of their attributes, such as patterns of expression, remain largely unknown, lncRNAs have key functions in transcriptional, post-transcriptional, and epigenetic gene regulation; Also, new work indicates their functions in scaffolding ribonuclear protein complexes. In plants, genome-wide identification of lncRNAs has been conducted in several species, including Zea mays, and recent research showed that lncRNAs regulate flowering time in the photoperiod pathway, and function in nodulation. In this review, we discuss the basic mechanisms by which lncRNAs regulate key cellular processes, using the large body of knowledge on animal and yeast lncRNAs to illustrate the significance of emerging work on lncRNAs in plants.

The hypothetical protein Ycf46 is involved in regulation of CO2 utilization in the cyanobacterium Synechocystis sp. PCC 6803

The Ycf46 mutant of Synechocystis showed growth inhibition under low dissolved CO 2 conditions, suggesting a role for the Ycf46 protein in the process of photosynthetic CO 2 uptake and utilization. Hypothetical chloroplast open reading frame Ycf46 proteins are highly conserved in all cyanobacterial lineages and most algal chloroplast genomes, but their exact function is still unknown. In the cyanobacterium Synechocystis sp. PCC 6803, the Ycf46 encoding gene slr0374 is part of an operon (with slr0373 and slr0376) and responds to many environmental stresses. Transcript levels of the slr0373, slr0374 and slr0376 genes were increased under a low concentration of dissolved inorganic carbon (Ci). Compared with the wild type, the mutant lacking slr0374 showed growth arrest under Ci-deficient conditions but not under iron-deficient or low-light conditions. In addition, the mutant grew more slowly than the wild type under pH 6.0 conditions in which CO2 was the dominant Ci source, indicating the mutant cells had weak CO2 uptake and/or utilization ability. Supplying a high concentration of CO2 (5 %, v/v) to the mutant restored its phenotype to the wild type level. The photosynthetic activity of the mutant was inhibited to a lesser extent by a carbonic anhydrase inhibitor than that of the wild type, which specifically blocked CO2 uptake. Inactivation of slr0374 decreased expression of the ecaB gene and reduced carbonic anhydrase activity. A subcellular localization assay indicated that the Ycf46 protein was soluble. By co-immunoprecipitation assay using Slr0374 as a bait-protein, potential interacting proteins in the size range of 30 kDa were identified. These results suggest that the Ycf46 protein plays a role in the regulation of photosynthesis in cyanobacteria, especially in CO2 uptake and utilization.

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.

Mesoscopic model and free energy landscape for protein-DNA binding sites: analysis of cyanobacterial promoters

The identification of protein binding sites in promoter sequences is a key problem to understand and control regulation in biochemistry and biotechnological processes. We use a computational method to analyze promoters from a given genome. Our approach is based on a physical model at the mesoscopic level of protein-DNA interaction based on the influence of DNA local conformation on the dynamics of a general particle along the chain. Following the proposed model, the joined dynamics of the protein particle and the DNA portion of interest, only characterized by its base pair sequence, is simulated. The simulation output is analyzed by generating and analyzing the Free Energy Landscape of the system. In order to prove the capacity of prediction of our computational method we have analyzed nine promoters of Anabaena PCC 7120. We are able to identify the transcription starting site of each of the promoters as the most populated macrostate in the dynamics. The developed procedure allows also to characterize promoter macrostates in terms of thermo-statistical magnitudes (free energy and entropy), with valuable biological implications. Our results agree with independent previous experimental results. Thus, our methods appear as a powerful complementary tool for identifying protein binding sites in promoter sequences.

Binding of specific proteins to particular sites in the DNA sequence is a fundamental issue for gene regulation in molecular biology and genetic engineering. A deep understanding of cell physiology requires the analysis of a plethora of genes involving characterization of their promoter architectures that determine their regulation and gene transcription. In order to locate the promoter elements of a given gene, experimental determination of its transcription start site (TSS) is required. This is an expensive, time-consuming task that, depending on our requirements, could be simplified using computational analysis as a first approach. Nevertheless, most computational methods lack a physical basis on the protein-DNA interaction mechanism. We adopt here this strategy, by using a simple model for protein-DNA interaction to find TSS in a bunch of cyanobacteria promoters. We make use of physical tools to characterize these TSS and to relate them with biological properties as the relative strength of the promoter. Our study shows how a model based on a coarse-grained description of a biomolecule can give valuable insight on its biological function.

The small regulatory RNA SyR1/PsrR1 controls photosynthetic functions in cyanobacteria

Little is known so far about RNA regulators of photosynthesis in plants, algae, or cyanobacteria. The small RNA PsrR1 (formerly SyR1) has been discovered in Synechocystis sp PCC 6803 and appears to be widely conserved within the cyanobacterial phylum. Expression of PsrR1 is induced shortly after a shift from moderate to high-light conditions. Artificial overexpression of PsrR1 led to a bleaching phenotype under moderate light growth conditions. Advanced computational target prediction suggested that several photosynthesis-related mRNAs could be controlled by PsrR1, a finding supported by the results of transcriptome profiling experiments upon pulsed overexpression of this small RNA in Synechocystis sp PCC 6803. We confirmed the interaction between PsrR1 and the ribosome binding regions of the psaL, psaJ, chlN, and cpcA mRNAs by mutational analysis in a heterologous reporter system. Focusing on psaL as a specific target, we show that the psaL mRNA is processed by RNase E only in the presence of PsrR1. Furthermore, we provide evidence for a posttranscriptional regulation of psaL by PsrR1 in the wild type at various environmental conditions and analyzed the consequences of PsrR1-based regulation on photosystem I. In summary, computational and experimental data consistently establish the small RNA PsrR1 as a regulatory factor controlling photosynthetic functions.

Alkane Biosynthesis Genes in Cyanobacteria and Their Transcriptional Organization

In cyanobacteria, alkanes are synthesized from a fatty acyl-ACP by two enzymes, acyl–acyl carrier protein reductase and aldehyde deformylating oxygenase. Despite the great interest in the exploitation for biofuel production, nothing is known about the transcriptional organization of their genes or the physiological function of alkane synthesis. The comparison of 115 microarray datasets indicates the relatively constitutive expression of aar and ado genes. The analysis of 181 available genomes showed that in 90% of the genomes both genes are present, likely indicating their physiological relevance. In 61% of them they cluster together with genes encoding acetyl-CoA carboxyl transferase and a short-chain dehydrogenase, strengthening the link to fatty acid metabolism and in 76% of the genomes they are located in tandem, suggesting constraints on the gene arrangement. However, contrary to the expectations for an operon, we found in Synechocystis sp. PCC 6803 specific promoters for the two genes, sll0208 (ado) and sll0209 (aar), which give rise to monocistronic transcripts. Moreover, the upstream located ado gene is driven by a proximal as well as a second, distal, promoter, from which a third transcript, the ~160 nt sRNA SyR9 is transcribed. Thus, the transcriptional organization of the alkane biosynthesis genes in Synechocystis sp. PCC 6803 is of substantial complexity. We verified all three promoters to function independently from each other and show a similar promoter arrangement also in the more distant Nodularia spumigena, Trichodesmium erythraeum, Anabaena sp. PCC 7120, Prochlorococcus MIT9313, and MED4. The presence of separate regulatory elements and the dominance of monocistronic mRNAs suggest the possible autonomous regulation of ado and aar. The complex transcriptional organization of the alkane synthesis gene cluster has possible metabolic implications and should be considered when manipulating the expression of these genes in cyanobacteria.

Comparative analysis of the primary transcriptome of Synechocystis sp. PCC 6803

RNA-seq and especially differential RNA-seq-type transcriptomic analyses (dRNA-seq) are powerful analytical tools, as they not only provide insights into gene expression changes but also provide detailed information about all promoters active at a given moment, effectively giving a deep insight into the transcriptional landscape. Synechocystis sp. PCC 6803 (Synechocystis 6803) is a unicellular model cyanobacterium that is widely used in research fields from ecology, photophysiology to systems biology, modelling and biotechnology. Here, we analysed the response of the Synechocystis 6803 primary transcriptome to different, environmentally relevant stimuli. We established genome-wide maps of the transcriptional start sites active under 10 different conditions relevant for photosynthetic growth and identified 4,091 transcriptional units, which provide information about operons, 5′ and 3′ untranslated regions (UTRs). Based on a unique expression factor, we describe regulons and relevant promoter sequences at single-nucleotide resolution. Finally, we report several sRNAs with an intriguing expression pattern and therefore likely function, specific for carbon depletion (CsiR1), nitrogen depletion (NsiR4), phosphate depletion (PsiR1), iron stress (IsaR1) or photosynthesis (PsrR1). This dataset is accompanied by comprehensive information providing extensive visualization and data access to allow an easy-to-use approach for the design of experiments, the incorporation into modelling studies of the regulatory system and for comparative analyses.

CopraRNA and IntaRNA: predicting small RNA targets, networks and interaction domains

CopraRNA (Comparative prediction algorithm for small RNA targets) is the most recent asset to the Freiburg RNA Tools webserver. It incorporates and extends the functionality of the existing tool IntaRNA (Interacting RNAs) in order to predict targets, interaction domains and consequently the regulatory networks of bacterial small RNA molecules. The CopraRNA prediction results are accompanied by extensive postprocessing methods such as functional enrichment analysis and visualization of interacting regions. Here, we introduce the functionality of the CopraRNA and IntaRNA webservers and give detailed explanations on their postprocessing functionalities. Both tools are freely accessible at http://rna.informatik.uni-freiburg.de.

Biocomputional construction of a gene network under acid stress in Synechocystis sp. PCC 6803

Acid stress is one of the most serious threats that cyanobacteria have to face, and it has an impact at all levels from genome to phenotype. However, very little is known about the detailed response mechanism to acid stress in this species. We present here a general analysis of the gene regulatory network of Synechocystis sp. PCC 6803 in response to acid stress using comparative genome analysis and biocomputational prediction. In this study, we collected 85 genes and used them as an initial template to predict new genes through co-regulation, protein-protein interactions and the phylogenetic profile, and 179 new genes were obtained to form a complete template. In addition, we found that 11 enriched pathways such as glycolysis are closely related to the acid stress response. Finally, we constructed a regulatory network for the intricate relationship of these genes and summarize the key steps in response to acid stress. This is the first time a bioinformatic approach has been taken systematically to gene interactions in cyanobacteria and the elaboration of their cell metabolism and regulatory pathways under acid stress, which is more efficient than a traditional experimental study. The results also provide theoretical support for similar research into environmental stresses in cyanobacteria and possible industrial applications.

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.

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

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

Riboregulators and the role of Hfq in photosynthetic bacteria

Anoxygenic and oxygenic bacteria directly convert solar energy into biomass using photosynthesis. The formation and composition of photosynthetic complexes has to be tightly controlled in response to environmental conditions, as exposure to sunlight can be harmful due to the generation of reactive oxygen species and the damaging effects of UV irradiation. Therefore, photosynthetic bacteria are exposed to a particular set of regulatory challenges in addition to those that also affect other bacteria, requiring sophisticated regulatory systems. Indeed, hundreds of potential regulatory RNAs have been identified in photosynthetic model bacteria as well as antisense RNAs (asRNAs) of up to several kb in length that protect certain mRNAs from degradation. The trans-acting small non-coding RNAs (sRNAs), PcrZ and PsrR1, control pigment and photosystem biogenesis in Rhodobacter sphaeroides and cyanobacteria, respectively. The asRNAs IsrR and As1_flv4 act as negative regulators and the asRNAs PsbA2R and PsbA3R as positive effectors of photosynthesis gene expression in Synechocystis 6803.

Transcriptomic response to prolonged ethanol production in the cyanobacterium Synechocystis sp. PCC6803

BACKGROUND

The production of biofuels in photosynthetic microalgae and cyanobacteria is a promising alternative to the generation of fuels from fossil resources. To be economically competitive, producer strains need to be established that synthesize the targeted product at high yield and over a long time. Engineering cyanobacteria into forced fuel producers should considerably interfere with overall cell homeostasis, which in turn might counteract productivity and sustainability of the process. Therefore, in-depth characterization of the cellular response upon long-term production is of high interest for the targeted improvement of a desired strain.

RESULTS

The transcriptome-wide response to continuous ethanol production was examined in Synechocystis sp. PCC6803 using high resolution microarrays. In two independent experiments, ethanol production rates of 0.0338% (v/v) ethanol d-1 and 0.0303% (v/v) ethanol d-1 were obtained over 18 consecutive days, measuring two sets of biological triplicates in fully automated photobioreactors. Ethanol production caused a significant (~40%) delay in biomass accumulation, the development of a bleaching phenotype and a down-regulation of light harvesting capacity. However, microarray analyses performed at day 4, 7, 11 and 18 of the experiment revealed only three mRNAs with a strongly modified accumulation level throughout the course of the experiment. In addition to the overexpressed adhA (slr1192) gene, this was an approximately 4 fold reduction in cpcB (sll1577) and 3 to 6 fold increase in rps8 (sll1809) mRNA levels. Much weaker modifications of expression level or modifications restricted to day 18 of the experiment were observed for genes involved in carbon assimilation (Ribulose bisphosphate carboxylase and Glutamate decarboxylase). Molecular analysis of the reduced cpcB levels revealed a post-transcriptional processing of the cpcBA operon mRNA leaving a truncated mRNA cpcA* likely not competent for translation. Moreover, western blots and zinc-enhanced bilin fluorescence blots confirmed a severe reduction in the amounts of both phycocyanin subunits, explaining the cause of the bleaching phenotype.

CONCLUSIONS

Changes in gene expression upon induction of long-term ethanol production in Synechocystis sp. PCC6803 are highly specific. In particular, we did not observe a comprehensive stress response as might have been expected.

A complex of Cas proteins 5, 6, and 7 is required for the biogenesis and stability of clustered regularly interspaced short palindromic repeats (crispr)-derived rnas (crrnas) in Haloferax volcanii

The clustered regularly interspaced short palindromic repeats/CRISPR-associated (CRISPR-Cas) system is a prokaryotic defense mechanism against foreign genetic elements. A plethora of CRISPR-Cas versions exist, with more than 40 different Cas protein families and several different molecular approaches to fight the invading DNA. One of the key players in the system is the CRISPR-derived RNA (crRNA), which directs the invader-degrading Cas protein complex to the invader. The CRISPR-Cas types I and III use the Cas6 protein to generate mature crRNAs. Here, we show that the Cas6 protein is necessary for crRNA production but that additional Cas proteins that form a CRISPR-associated complex for antiviral defense (Cascade)-like complex are needed for crRNA stability in the CRISPR-Cas type I-B system in Haloferax volcanii in vivo. Deletion of the cas6 gene results in the loss of mature crRNAs and interference. However, cells that have the complete cas gene cluster (cas1–8b) removed and are transformed with the cas6 gene are not able to produce and stably maintain mature crRNAs. crRNA production and stability is rescued only if cas5, -6, and -7 are present. Mutational analysis of the cas6 gene reveals three amino acids (His-41, Gly-256, and Gly-258) that are essential for pre-crRNA cleavage, whereas the mutation of two amino acids (Ser-115 and Ser-224) leads to an increase of crRNA amounts. This is the first systematic in vivo analysis of Cas6 protein variants. In addition, we show that the H. volcanii I-B system contains a Cascade-like complex with a Cas7, Cas5, and Cas6 core that protects the crRNA.

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.

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.

Predicting sRNAs and their targets in bacteria

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

Computational identification and functional predictions of long noncoding RNA in Zea mays

BACKGROUND

Computational analysis of cDNA sequences from multiple organisms suggests that a large portion of transcribed DNA does not code for a functional protein. In mammals, noncoding transcription is abundant, and often results in functional RNA molecules that do not appear to encode proteins. Many long noncoding RNAs (lncRNAs) appear to have epigenetic regulatory function in humans, including HOTAIR and XIST. While epigenetic gene regulation is clearly an essential mechanism in plants, relatively little is known about the presence or function of lncRNAs in plants.

METHODOLOGY/PRINCIPAL FINDINGS

To explore the connection between lncRNA and epigenetic regulation of gene expression in plants, a computational pipeline using the programming language Python has been developed and applied to maize full length cDNA sequences to identify, classify, and localize potential lncRNAs. The pipeline was used in parallel with an SVM tool for identifying ncRNAs to identify the maximal number of ncRNAs in the dataset. Although the available library of sequences was small and potentially biased toward protein coding transcripts, 15% of the sequences were predicted to be noncoding. Approximately 60% of these sequences appear to act as precursors for small RNA molecules and may function to regulate gene expression via a small RNA dependent mechanism. ncRNAs were predicted to originate from both genic and intergenic loci. Of the lncRNAs that originated from genic loci, ∼20% were antisense to the host gene loci.

CONCLUSIONS/SIGNIFICANCE

Consistent with similar studies in other organisms, noncoding transcription appears to be widespread in the maize genome. Computational predictions indicate that maize lncRNAs may function to regulate expression of other genes through multiple RNA mediated mechanisms.

Non-coding RNAs in marine Synechococcus and their regulation under environmentally relevant stress conditions

Regulatory small RNAs (sRNAs) have crucial roles in the adaptive responses of bacteria to changes in the environment. Thus far, potential regulatory RNAs have been studied mainly in marine picocyanobacteria in genetically intractable Prochlorococcus, rendering their molecular analysis difficult. Synechococcus sp. WH7803 is a model cyanobacterium, representative of the picocyanobacteria from the mesotrophic areas of the ocean. Similar to the closely related Prochlorococcus it possesses a relatively streamlined genome and a small number of genes, but is genetically tractable. Here, a comparative genome analysis was performed for this and four additional marine Synechococcus to identify the suite of possible sRNAs and other RNA elements. Based on the prediction and on complementary microarray profiling, we have identified several known as well as 32 novel sRNAs. Some sRNAs overlap adjacent coding regions, for instance for the central photosynthetic gene psbA. Several of these novel sRNAs responded specifically to environmentally relevant stress conditions. Among them are six sRNAs changing their accumulation level under cold stress, six responding to high light and two to iron limitation. Target predictions suggested genes encoding components of the light-harvesting apparatus as targets of sRNAs originating from genomic islands and that one of the iron-regulated sRNAs might be a functional homolog of RyhB. These data suggest that marine Synechococcus mount adaptive responses to these different stresses involving regulatory sRNAs.

Comprehensive transcriptome analysis reveals novel genes involved in cardiac glycoside biosynthesis and mlncRNAs associated with secondary metabolism and stress response in Digitalis purpurea

Abstract

Conclusions

Through comprehensive transcriptome analysis, we not only identified 29 novel gene families potentially involved in the biosynthesis of cardiac glycosides but also characterized a large number of mlncRNAs. Our results suggest the importance of mlncRNAs in secondary metabolism and stress response in D. purpurea.

Small RNAs of the Bradyrhizobium/Rhodopseudomonas lineage and their analysis

Small RNAs (sRNAs) play a pivotal role in bacterial gene regulation. However, the sRNAs of the vast majority of bacteria with sequenced genomes still remain unknown since sRNA genes are usually difficult to recognize and thus not annotated. Here, expression of seven sRNAs (BjrC2a, BjrC2b, BjrC2c, BjrC68, BjrC80, BjrC174 and BjrC1505) predicted by genome comparison of Bradyrhizobium and Rhodopseudomonas members, was verified by RNA gel blot hybridization, microarray and deep sequencing analyses of RNA from the soybean symbiont Bradyrhizobium japonicum USDA 110. BjrC2a, BjrC2b and BjrC2c belong to the RNA family RF00519, while the other sRNAs are novel. For some of the sRNAs we observed expression differences between free-living bacteria and bacteroids in root nodules. The amount of BjrC1505 was decreased in nodules. By contrast, the amount of BjrC2a, BjrC68, BjrC80, BjrC174 and the previously described 6S RNA was increased in nodules, and accumulation of truncated forms of these sRNAs was observed. Comparative genomics and deep sequencing suggest that BjrC2a is an antisense RNA regulating the expression of inositol-monophosphatase. The analyzed sRNAs show a different degree of conservation in Rhizobiales, and expression of homologs of BjrC2, BjrC68, BjrC1505, and 6S RNA was confirmed in the free-living purple bacterium Rhodopseudomonas palustris 5D.

Identification of three novel antisense RNAs in the fur locus from unicellular cyanobacteria

The interplay between Fur (ferric uptake regulator) proteins and small, non-coding RNAs has been described as a key regulatory loop in several bacteria. In the filamentous cyanobacterium Anabaena sp. PCC 7120, a large dicistronic transcript encoding the putative membrane protein Alr1690 and an α-furA RNA is involved in the modulation of the global regulator FurA. In this work we report the existence of three novel antisense RNAs in cyanobacteria and show that a cis α-furA RNA is conserved in very different genomic contexts, namely in the unicellular cyanobacteria Microcystis aeruginosa PCC 7806 and Synechocystis sp. PCC 6803. Syα-fur RNA covers only part of the coding sequence of the fur orthologue sll0567, whose flanking genes encode two hypothetical proteins. Transcriptional analysis of fur and its adjacent genes in Microcystis unravels a highly compact organization of this locus involving overlapping transcripts. Maα-fur RNA spans the whole Mafur CDS and part of the flanking dnaJ and sufE sequences. In addition, Mafur seems to be part of a dicistronic operon encoding this regulator and an α-sufE RNA. These results allow new insights into the transcriptomes of two unicellular cyanobacteria and suggest that in M. aeruginosa PCC 7806, the α-fur and α-sufE RNAs might participate in a regulatory connection between the genes of the dnaJ–fur–sufE locus.

Conservation and Occurrence of Trans-Encoded sRNAs in the Rhizobiales

Post-transcriptional regulation by trans-encoded sRNAs, for example via base-pairing with target mRNAs, is a common feature in bacteria and influences various cell processes, e.g., response to stress factors. Several studies based on computational and RNA-seq approaches identified approximately 180 trans-encoded sRNAs in Sinorhizobium meliloti. The initial point of this report is a set of 52 trans-encoded sRNAs derived from the former studies. Sequence homology combined with structural conservation analyses were applied to elucidate the occurrence and distribution of conserved trans-encoded sRNAs in the order of Rhizobiales. This approach resulted in 39 RNA family models (RFMs) which showed various taxonomic distribution patterns. Whereas the majority of RFMs was restricted to Sinorhizobium species or the Rhizobiaceae, members of a few RFMs were more widely distributed in the Rhizobiales. Access to this data is provided via the RhizoGATE portal [1,2].

Cyanobacterial genomics for ecology and biotechnology

Cyanobacteria are the only prokaryotes that directly convert solar energy and CO(2) into organic matter by oxygenic photosynthesis, explaining their relevance for primary production in many ecosystems and the increasing interest for biotechnology. At present, there are more than 60 cyanobacteria for which a total genome sequence is publicly available. These cyanobacteria belong to different lifestyles and origins, coming from marine and freshwater aquatic environments, as well as terrestrial and symbiotic habitats. Genome sizes vary by a factor of six, from 1.44 Mb to 9.05 Mb, with the number of reported genes ranging from 1241 to 8462. Several studies have demonstrated how these sequences could be used to successfully infer important ecological, physiological and biotechnologically relevant characteristics. However, sequences of cyanobacterial origin also comprise a significant portion of certain metagenomes. Moreover, genome analysis has been employed for culture-independent approaches and for resequencing mutant strains, a very recent tool in cyanobacterial research.

Small RNA-mediated control of the Agrobacterium tumefaciens GABA binding protein

Wounded plants activate a complex defence programme in response to Agrobacterium tumefaciens. They synthesize the non-proteinogenic amino acid γ-aminobutyric acid (GABA), which stimulates degradation of the quorum sensing signal N-(3-oxo-octanoyl) homoserine lactone. GABA is transported into A. tumefaciens via an ABC transporter dependent on the periplasmic binding protein Atu2422. We demonstrate that expression of atu2422 and two other ABC transporter genes is downregulated by the conserved small RNA (sRNA) AbcR1 (for ABC regulator). AbcR1 is encoded in tandem with another sRNA, which is similar in sequence and structure. Both sRNAs accumulate during stationary phase but only the absence of AbcR1 resulted in significant accumulation of Atu2422 and increased GABA import. AbcR1 inhibits initiation of atu2422 translation by masking its Shine-Dalgarno sequence and thereby reduces stability of the atu2422 transcript. It is the first described bacterial sRNA that controls uptake of a plant-generated signalling molecule. Given that similar sRNAs and ABC transporter genes are present in various Rhizobiaceae and in Brucella, it is likely that such sRNA-mediated control impacts a number of host-microbe interactions.

Regulatory RNAs in cyanobacteria: developmental decisions, stress responses and a plethora of chromosomally encoded cis-antisense RNAs

Cyanobacteria are the only prokaryotes which directly convert solar energy into biomass using oxygenic photosynthesis. Therefore, these bacteria are of interest for the production of biofuels, biotechnology and are of tremendous relevance for primary carbon fixation in many ecosystems. Mechanisms controlling gene expression cannot be understood entirely without information on the numbers and functions of regulatory RNAs. In cyanobacteria, non-coding RNAs have been characterized from simple unicellular species such as Prochlorococcus up to complex species such as Anabaena. Several of these RNAs function in the control of stress responses, photosynthesis, outer cell membrane protein biosynthesis and the differentiation of cells.