Björn Voß

Zusammenfassung

RNA–RNA inter- and intramolecular interactions are fundamental for numerous biological processes. While there are reasonable approaches to map RNA secondary structures genome-wide, understanding how different RNAs interact to carry out their regulatory functions requires mapping of intermolecular base pairs. Recently, different strategies to detect RNA–RNA duplexes in living cells, so called direct duplex detection (DDD) methods, have been developed. Common to all is the Psoralen-mediated in vivo RNA crosslinking followed by RNA Proximity Ligation to join the two interacting RNA strands. Sequencing of the RNA via classical RNA-seq and subsequent specialised bioinformatic analyses the result in the prediction of inter- and intramolecular RNA–RNA interactions. Existing approaches adapt standard RNA-seq analysis pipelines, but often neglect inherent features of RNA–RNA interactions that are useful for filtering and statistical assessment. Here we present RNAnue, a general pipeline for the inference of RNA–RNA interactions from DDD experiments that takes into account hybridisation potential and statistical significance to improve prediction accuracy. We applied RNAnue to data from different DDD studies and compared our results to those of the original methods. This showed that RNAnue performs better in terms of quantity and quality of predictions.

Zusammenfassung

The function and mode of action of small regulatory RNAs is currently still understudied in archaea. In the halophilic archaeon H. volcanii a plethora of sRNAs have been identified, however, in-depth functional analysis is missing for most of them. We selected a small RNA (s479) from H. volcanii for detailed characterization. The sRNA gene is encoded between a CRISPR RNA locus and the Cas protein gene cluster, the s479 deletion strain is viable and was characterized in detail. Transcriptome studies of wild type Haloferax cells and the deletion mutant revealed up-regulation of six genes in the deletion strain, showing that the sRNA has a clearly defined function. Three of the six up-regulated genes encode potential zinc transporter proteins (ZnuA1, ZnuB1, ZnuC1) suggesting involvement of s479 in regulation of zinc transport. Upregulation of these genes in the deletion strain was confirmed by northern blot and proteome analyses. Furthermore, electrophoretic mobility shift assays demonstrate a direct interaction of s479 with the target znuC1 mRNA. Proteome comparison of wild type and deletion strains further expanded the regulon of s479 deeply rooting this sRNA within the metabolism of H. volcanii especially the regulation of transporter abundance. Interestingly, s479 is not only encoded next to CRISPR-cas genes but the mature s479 contains a crRNA-like 5´ handle and experiments with Cas protein deletion strains indicate maturation by Cas6 and interaction with Cas proteins. Together this might suggest that the CRISPR-Cas system is involved in s479 function.

Zusammenfassung

The correct prediction of bacterial sRNA homologs is a prerequisite for many downstream analyses based on comparative genomics, but it is frequently challenging due to the short length and distinct heterogeneity of such homologs. GLASSGo is an efficient tool for the prediction of sRNA homologs from a single input query. To make the algorithm available to a broader community, we offer a Docker container along with a free-access web service. For non-computer scientists, the web service provides a user-friendly interface. However, capabilities were lacking so far for batch processing, version control, and direct interaction with compatible software applications as a workflow management system can provide.Here we present GLASSGo 1.5.2, an updated version that is fully incorporated into the workflow management system Galaxy. The improved version contains a new feature for extracting the upstream regions, allowing the search for conserved promoter elements. Additionally, it supports the use of accession numbers instead of the outdated GI numbers, which widens the applicability of the tool.GLASSGo is available at https://github.com/lotts/GLASSgo/ under the MIT license and is accompanied by instruction and application data. Furthermore, it can be installed into any Galaxy instance using the Galaxy ToolShed.

Zusammenfassung

Pseudomonas putida is recognized as a very promising strain for industrial application due to its high redox capacity and frequently observed tolerance towards organic solvents. In this research, we studied the metabolic and transcriptional response of P. putida KT2440 exposed to large-scale heterogeneous mixing conditions in the form of repeated glucose shortage. Cellular responses were mimicked in an experimental setup comprising a stirred tank reactor and a connected plug flow reactor. We deciphered that a stringent response-like transcriptional regulation programme is frequently induced, which seems to be linked to the intracellular pool of 3-hydroxyalkanoates (3-HA) that are known to serve as precursors for polyhydroxyalkanoates (PHA). To be precise, P. putida is endowed with a survival strategy likely to access cellular PHA, amino acids and glycogen in few seconds under glucose starvation to obtain ATP from respiration, thereby replenishing the reduced ATP levels and the adenylate energy charge. Notably, cells only need 0.4\% of glucose uptake to build those 3-HA-based energy buffers. Concomitantly, genes that are related to amino acid catabolism and $\beta$-oxidation are upregulated during the transient absence of glucose. Furthermore, we provide a detailed list of transcriptional short- and long-term responses that increase the cellular maintenance by about 17\% under the industrial-like conditions tested.

Zusammenfassung

Clustered regularly interspaced short palindromic repeats (CRISPR) and their associated proteins (Cas) are essential genetic elements in many archaeal and bacterial genomes, playing a key role in a prokaryote adaptive immune system against invasive foreign elements. In recent years, the CRISPR-Cas system has also been engineered to facilitate target gene editing in eukaryotic genomes. Bioinformatics played an essential role in the detection and analysis of CRISPR systems and here we review the bioinformatics-based efforts that pushed the field of CRISPR-Cas research further. We discuss the bioinformatics tools that have been published over the last few years and, finally, present the most popular tools for the design of CRISPR-Cas9 guides.

Zusammenfassung

Bacterial small RNAs (sRNAs) are important post-transcriptional regulators of gene expression. The functional and evolutionary characterization of sRNAs requires the identification of homologs, which is frequently challenging due to their heterogeneity, short length and partly, little sequence conservation. We developed the GLobal Automatic Small RNA Search go (GLASSgo) algorithm to identify sRNA homologs in complex genomic databases starting from a single sequence. GLASSgo combines an iterative BLAST strategy with pairwise identity filtering and a graph-based clustering method that utilizes RNA secondary structure information. We tested the specificity, sensitivity and runtime of GLASSgo, BLAST and the combination RNAlien/cmsearch in a typical use case scenario on 40 bacterial sRNA families. The sensitivity of the tested methods was similar, while the specificity of GLASSgo and RNAlien/cmsearch was significantly higher than that of BLAST. GLASSgo was on average ~87 times faster than RNAlien/cmsearch, and only ~7.5 times slower than BLAST, which shows that GLASSgo optimizes the trade-off between speed and accuracy in the task of finding sRNA homologs. GLASSgo is fully automated, whereas BLAST often recovers only parts of homologs and RNAlien/cmsearch requires extensive additional bioinformatic work to get a comprehensive set of homologs. GLASSgo is available as an easy-to-use web server to find homologous sRNAs in large databases.

Zusammenfassung

De-novo assembly of a metagenomic dataset obtained from a decaying cyanobacterial Trichodesmium bloom from the New Caledonian lagoon resulted in a complete giant phage genome of 257,908bp, obtained independently with multiple assembly tools. Noteworthy, gammaproteobacteria were an abundant fraction in the sequenced samples. Mapping of the raw reads with 99% accuracy to the giant phage genome resulted in an average coverage of 262X. The closest sequenced relatives, albeit still distant, are the Pseudomonas phages PaBG from Lake Baikal and Lu11 isolated from a soil sample from the Philippines. The phage reported here might belong to the same family within the Myoviridae as PaBG and Lu11 and would thus be its first marine member, indicating a more widespread occurrence of this group. We named this phage NCTB (New Caledonia Trichodesmium Bloom) after its origin.

Zusammenfassung

For eukaryotes there seems to be no doubt that differences on the trancriptomic level substantially contribute to the process of species diversification, whereas for bacteria this is thought to be less important. Recent years saw a significant increase in full transcriptome studies for bacteria, which provided deep insight into the architecture of bacterial transcriptomes. Most notably, it became evident that, in contrast to previous scientific consensus, bacterial transcriptomes are quite complex. There exist a large number of cis-antisense RNAs, non-coding RNAs, overlapping transcripts and RNA elements that regulate transcription, such as riboswitches. Furthermore, processing and degradation of RNA has gained interest, because it has a significant impact on the composition of the transcriptome. In this review, we summarize recent findings and put them into a broader context with respect to the complexity of bacterial transcriptomes and its putative biological meanings.

Zusammenfassung

Motivation: We developed VisualGraphX, a web-based, interactive visualization tool for large-scale graphs. Current graph visualization tools that follow the rich-internet paradigm lack an interactive and scalable visualization of graph-based data. VisualGraphX aims to provide a universal graph visualization tool that empowers the users to efficiently explore the data for themselves at a large scale. It is available as a visualization plugin for the Galaxy platform, such that VisualGraphX can be integrated into custom analysis pipelines. Availability and Implementation: VisualGraphX has been released as a visualization plugin for the Galaxy platform under AFL 3.0 and is available with instructions and application data at http://gitlab.com/comptrans/VisualGraphX/. Contact: bjoern.voss@ibvt.uni-stuttgart.de

Zusammenfassung

The visualization of massive datasets, such as those resulting from comparative metatranscriptome analyses or the analysis of microbial population structures using ribosomal RNA sequences, is a challenging task. We developed a new method called CoVennTree (Comparative weighted Venn Tree) that simultaneously compares up to three multifarious datasets by aggregating and propagating information from the bottom to the top level and produces a graphical output in Cytoscape. With the introduction of weighted Venn structures, the contents and relationships of various datasets can be correlated and simultaneously aggregated without losing information. We demonstrate the suitability of this approach using a dataset of 16S rDNA sequences obtained from microbial populations at three different depths of the Gulf of Aqaba in the Red Sea. CoVennTree has been integrated into the Galaxy ToolShed and can be directly downloaded and integrated into the user instance.

Zusammenfassung

Synechocystis sp. PCC 6803 is the most popular cyanobacterial model for prokaryotic photosynthesis and for metabolic engineering to produce biofuels. Genomic and transcriptomic comparisons between closely related bacteria are powerful approaches to infer insights into their metabolic potentials and regulatory networks. To enable a comparative approach, we generated the draft genome sequence of Synechocystis sp. PCC 6714, a closely related strain of 6803 (16S rDNA identity 99.4\%) that also is amenable to genetic manipulation. Both strains share 2838 protein-coding genes, leaving 845 unique genes in Synechocystis sp. PCC 6803 and 895 genes in Synechocystis sp. PCC 6714. The genetic differences include a prophage in the genome of strain 6714, a different composition of the pool of transposable elements, and a ∼40 kb genomic island encoding several glycosyltransferases and transport proteins. We verified several physiological differences that were predicted on the basis of the respective genome sequence. Strain 6714 exhibited a lower tolerance to Zn2+ ions, associated with the lack of a corresponding export system and a lowered potential of salt acclimation due to the absence of a transport system for the re-uptake of the compatible solute glucosylglycerol. These new data will support the detailed comparative analyses of this important cyanobacterial group than has been possible thus far. Genome information for Synechocystis sp. PCC 6714 has been deposited in Genbank (accession no AMZV01000000).

Zusammenfassung

Synechocystis sp. strain PCC 6714 is a unicellular cyanobacterium closely related to the popular model organism Synechocystis sp. strain PCC 6803. A combination of PacBio SMRT and Illumina GAIIx data results in a highly accurate finished genome sequence that provides a reliable resource for further comparative analyses.

Zusammenfassung

Marine diatoms constitute a major component of eukaryotic phytoplankton and stand at the crossroads of several evolutionary lineages. These microalgae possess peculiar genomic features and novel combinations of genes acquired from bacterial, animal and plant ancestors. Furthermore, they display both DNA methylation and gene silencing activities. Yet, the biogenesis and regulatory function of small RNAs (sRNAs) remain ill defined in diatoms. PMID: 25142710

Zusammenfassung

RNA-seq and its variant differential RNA-seq (dRNA-seq) are today routine methods for transcriptome analysis in bacteria. While expression profiling and transcriptional start site prediction are standard tasks today, the problem of identifying transcriptional units in a genome-wide fashion is still not solved for prokaryotic systems. PMID: 24780064

Zusammenfassung

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 4091 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.

Zusammenfassung

RNA molecules, especially non-coding RNAs, play vital roles in the cell and their biological functions are mostly determined by structural properties. Often, these properties are related to dynamic changes in the structure, as in the case of riboswitches, and thus the analysis of RNA folding kinetics is crucial for their study. Exact approaches to kinetic folding are computationally expensive and, thus, limited to short sequences. In a previous study, we introduced a position-specific abstraction based on helices which we termed helix index shapes (hishapes) and a hishape-based algorithm for near-optimal folding pathway computation, called HiPath. The combination of these approaches provides an abstract view of the folding space that offers information about the global features. PMID: 24575751

Zusammenfassung

Prochlorococcus is a genus of abundant and ecologically important marine cyanobacteria. Here, we present a comprehensive comparison of the structure and composition of the transcriptomes of two Prochlorococcus strains, which, despite their similarities, have adapted their gene pool to specific environmental constraints. We present genome-wide maps of transcriptional start sites (TSS) for both organisms, which are representatives of the two most diverse clades within the two major ecotypes adapted to high- and low-light conditions, respectively. Our data suggest antisense transcription for three-quarters of all genes, which is substantially more than that observed in other bacteria. We discovered hundreds of TSS within genes, most notably within 16 of the 29 prochlorosin genes, in strain MIT9313. A direct comparison revealed very little conservation in the location of TSS and the nature of non-coding transcripts between both strains. We detected extremely short 5′ untranslated regions with a median length of only 27 and 29 nt for MED4 and MIT9313, respectively, and for 8\% of all protein-coding genes the median distance to the start codon is only 10 nt or even shorter. These findings and the absence of an obvious Shine–Dalgarno motif suggest that leaderless translation and ribosomal protein S1-dependent translation constitute alternative mechanisms for translation initiation in Prochlorococcus. We conclude that genome-wide antisense transcription is a major component of the transcriptional output from these relatively small genomes and that a hitherto unrecognized high degree of complexity and variability of gene expression exists in their transcriptional architecture.

Zusammenfassung

Nodularia spumigena is a filamentous diazotrophic cyanobacterium that dominates the annual late summer cyanobacterial blooms in the Baltic Sea. But N. spumigena also is common in brackish water bodies worldwide, suggesting special adaptation allowing it to thrive at moderate salinities. A draft genome analysis of N. spumigena sp. CCY9414 yielded a single scaffold of 5,462,271 nucleotides in length on which genes for 5,294 proteins were annotated. A subsequent strand-specific transcriptome analysis identified more than 6,000 putative transcriptional start sites (TSS). Orphan TSSs located in intergenic regions led us to predict 764 non-coding RNAs, among them 70 copies of a possible retrotransposon and several potential RNA regulators, some of which are also present in other N2-fixing cyanobacteria. Approximately 4\% of the total coding capacity is devoted to the production of secondary metabolites, among them the potent hepatotoxin nodularin, the linear spumigin and the cyclic nodulapeptin. The transcriptional complexity associated with genes involved in nitrogen fixation and heterocyst differentiation is considerably smaller compared to other Nostocales. In contrast, sophisticated systems exist for the uptake and assimilation of iron and phosphorus compounds, for the synthesis of compatible solutes, and for the formation of gas vesicles, required for the active control of buoyancy. Hence, the annotation and interpretation of this sequence provides a vast array of clues into the genomic underpinnings of the physiology of this cyanobacterium and indicates in particular a competitive edge of N. spumigena in nutrient-limited brackish water ecosystems.

Zusammenfassung

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.

Zusammenfassung

Synechocystis sp. PCC 6803 is a widely used model cyanobacterium for studying photosynthesis, phototaxis, the production of biofuels and many other aspects. Here we present a re-sequencing study of the genome and seven plasmids of one of the most widely used Synechocystis sp. PCC 6803 substrains, the glucose tolerant and motile Moscow or ‘PCC-M’ strain, revealing considerable evidence for recent microevolution. Seven single nucleotide polymorphisms (SNPs) specifically shared between ‘PCC-M’ and the ‘PCC-N and PCC-P’ substrains indicate that ‘PCC-M’ belongs to the ‘PCC’ group of motile strains. The identified indels and SNPs in ‘PCC-M’ are likely to affect glucose tolerance, motility, phage resistance, certain stress responses as well as functions in the primary metabolism, potentially relevant for the synthesis of alkanes. Three SNPs in intergenic regions could affect the promoter activities of two protein-coding genes and one cis-antisense RNA. Two deletions in ‘PCC-M’ affect parts of clustered regularly interspaced short palindrome repeats-associated spacer-repeat regions on plasmid pSYSA, in one case by an unusual recombination between spacer sequences.

Zusammenfassung

RNA has many pivotal functions especially in the regulation of gene expression by ncRNAs. Identification of their structure is an important requirement for understanding their function. Structure prediction alone is often insufficient for this task, due to algorithmic problems, parameter inaccuracies, and biological peculiarities. Among the latter, there are base modifications, cotranscriptional folding leading to folding traps, and conformational switching as in the case of riboswitches. All these require more in-depth analysis of the folding space. The major drawback, which all methods have to cope with, is the exponential growth of the folding space. Therefore, methods are often limited in the sequence length they can analyze, or they make use of heuristics, sampling, or abstraction. Our approach adopts the abstraction strategy and remedies some problems of existing methods. We introduce a position-specific abstraction based on helices that we term helix index shapes, or hishapes for short. Utilizing a dynamic programming framework, we have implemented this abstraction in the program RNAHeliCes. Furthermore, we developed two hishape-based methods, one for energy barrier estimation, called HiPath, and one for abstract structure comparison, termed HiTed. We demonstrate the superior performance of HiPath compared to other existing methods and the competitive accuracy of HiTed. RNAHeliCes, together with HiPath and HiTed, are available for download at http://www.cyanolab.de/software/RNAHeliCes.htm.

Zusammenfassung

There has been an increasing interest in cyanobacteria because these photosynthetic organisms convert solar energy into biomass and because of their potential for the production of biofuels. However, the exploitation of cyanobacteria for bioengineering requires knowledge of their transcriptional organization. Using differential RNA sequencing, we have established a genome-wide map of 3,527 transcriptional start sites (TSS) of the model organism Synechocystis sp. PCC6803. One-third of all TSS were located upstream of an annotated gene; another third were on the reverse complementary strand of 866 genes, suggesting massive antisense transcription. Orphan TSS located in intergenic regions led us to predict 314 noncoding RNAs (ncRNAs). Complementary microarray-based RNA profiling verified a high number of noncoding transcripts and identified strong ncRNA regulations. Thus, ∼64\% of all TSS give rise to antisense or ncRNAs in a genome that is to 87\% protein coding. Our data enhance the information on promoters by a factor of 40, suggest the existence of additional small peptide-encoding mRNAs, and provide corrected 5′ annotations for many genes of this cyanobacterium. The global TSS map will facilitate the use of Synechocystis sp. PCC6803 as a model organism for further research on photosynthesis and energy research.

Zusammenfassung

RNA has many pivotal functions especially in the regulation of gene expression by ncRNAs. Identification of their structure is an important requirement for understanding their function. Structure prediction alone is often insufficient for this task, due to algorithmic problems, parameter inaccuracies and also due to biological peculiarities. Among the latter we have base modifications, co-transcriptional folding leading to folding traps and conformational switching as in the case of riboswitches. All these require more in-depth analysis of the folding space. The major drawback, which all methods have to cope with, is the exponential growth of the folding space. Therefore, methods are often limited in the sequence length they can analyze or make use of heuristics, sampling or abstraction. Our approach adopts the abstraction strategy and remedies some problems of existing methods. We introduce a position-specific abstraction based on helices which we term helix index shapes or hishapes for short. Utilizing a dynamic programming framework, we have implemented this abstraction in the program RNAHeliCes. Furthermore, we present results of empirical studies on the size of the search space, showing that it is a reasonable extension to existing methods. We further show the application of RNAHeliCes to some well-studied classes of RNAs and discuss possible use cases of our method. RNAHeliCes is available for download at http://www.cyanolab.de/software/RNAHeliCes.htm.

Zusammenfassung

The ability of some RNA molecules to switch between different metastable conformations plays an important role in cellular processes. In order to identify such molecules and to predict their conformational changes one has to investigate the refolding pathways. As a qualitative measure of these transitions, the barrier height marks the energy peak along such refolding paths. We introduce a meta-heuristic to estimate such barriers, which is an NP-complete problem. To guide an arbitrary path heuristic, the method uses RNA shape representative structures as intermediate checkpoints for detours. This enables a broad but estimationcient search for refolding pathways. The resulting Shape Triples meta-heuristic enables a close to optimal estimation of the barrier height that outperforms the precision of the employed path heuristic.

Zusammenfassung

In response to nitrogen deficiency, some cyanobacteria develop heterocysts, a terminally differentiated cell type, specialized for the fixation of atmospheric nitrogen. In Nostocales, this differentiation process is controlled by two major regulators, NtcA and HetR, but additional unknown factors are likely to be involved as well. In the context of a genome-wide search for potential non-coding RNAs, we identified an array of 12 tandem repeats that is transcribed in large amounts when cells enter conditions that trigger cell differentiation and switch to nitrogen fixation. The main accumulating transcript, which we suggest designating nitrogen stress-induced RNA 1 (NsiR1), has properties similar to regulatory non-coding RNAs. In Anabaena sp. PCC 7120, it is about 60 nt in length, has a very distinct predicted secondary structure, and is expressed very early and transiently after nitrogen step-down. Moreover, its expression requires HetR and NtcA and is restricted to cells that are differentiating into heterocysts, clearly placing NsiR1 within the regulon that controls the switch to nitrogen fixation and heterocyst formation. The genomic arrangement of NsiR1, located upstream of hetF, a gene whose product is involved in heterocyst formation, is conserved in all five Nostocales whose genomes are completely sequenced. Additionally, we detected NsiR1 expression in 19 different heterocyst-forming cyanobacteria. Our data suggest that every repeat is a complete transcriptional unit furnished with a cell-type-specific promoter and a Rho-independent terminator, which gives rise to a very high NsiR1 transcript level. NsiR1 is the first known bacterial non-coding RNA that is specifically upregulated in response to nitrogen step-down.

Zusammenfassung

Retrograde signaling is a pathway of communication from mitochondria and plastids to the nucleus in the context of cell differentiation, development and stress response. In Chlamydomonas reinhardtii, the tetrapyrroles Mg-protoporphyrin IX (MgProto) and heme are only synthesized within the chloroplast, and they have been implicated in the retrograde control of nuclear gene expression in this unicellular green alga. Feeding the two tetrapyrroles to Chlamydomonas cultures was previously shown to transiently induce five nuclear genes, three of which encode the heat shock proteins HSP70A, B, and E. In contrast, controversial results exist on the possible role of MgProto in the repression of genes for light harvesting proteins in higher plants, raising the question of how important this mode of regulation is. Here we used genome-wide transcriptional profiling to measure the global impact of these tetrapyrroles on gene regulation and the scope of the response. We identified almost 1000 genes whose expression level changed transiently but significantly. Among them were only a few genes for photosynthetic proteins but several encoding enzymes of the TCA cycle, heme-binding proteins, stress-response proteins, as well as proteins involved in protein folding and degradation. More than \textgreater50\% of the latter class of genes were also regulated by heat shock. The observed drastic fold changes at RNA level did not correlate with similar changes of protein concentrations under the tested experimental conditions. Phylogenetic profiling revealed that genes of putative endosymbiontic origin are not overrepresented among the responding genes. This and the transient nature of changes in gene expression suggest a signaling role of both tetrapyrroles as secondary messengers for adaptive responses affecting the entire cell and not only organellar proteins.

Zusammenfassung

Based on computational prediction of RNA secondary structures, a long antisense RNA (asRNA) was found in chloroplasts of Arabidopsis, Nicotiana tabacum and poplar, which occurs in two to three major transcripts. Mapping of primary 5' ends, northern hybridizations and quantitative real-time reverse transcription polymerase chain reaction (qPCR) experiments demonstrated that these transcripts originate from a promoter that is typical for the plastid-encoded RNA polymerase and are over their full length in antisense orientation to the gene ndhB and therefore were designated asRNA\_ndhB. The asRNA\_ndhB transcripts predominantly accumulate in young leaves and at physiological growth temperatures. Two nucleotide positions in the mRNA that are subject to C-to-U RNA editing and which were previously found to be sensitive to elevated temperatures are covered by asRNA\_ndhB. Nevertheless, the correlation between the accumulation of asRNA\_ndhB and RNA editing appeared weak in a temperature shift experiment. With asRNA\_ndhB, we describe the first asRNA of plant chloroplasts that covers RNA editing sites, as well as a group II intron splice acceptor site, and that is under developmental control, raising the possibility that long asRNAs could be involved in RNA maturation or the control of RNA stability.

Zusammenfassung

Regulatory RNA plays a pivotal role in the regulation of bacterial gene expression. Here, five small RNAs were studied in Sinorhizobium meliloti - SmrC15, SmrC16, Sra33, 6S and the signal recognition particle (SRP) RNA, which are conserved among at least seven different Rhizobium and Sinorhizobium species. The amount of SmrC16 decreased in stationary phase, while the other RNAs were up-regulated. The smallest changes, maximally 2-fold, were observed for 6S RNA. In the distantly related Bradyrhizobium japonicum, the amount of 6S RNA was not increased in stationary phase, suggesting some functional divergence in the roles of this molecule in Rhizobiales in comparison to Escherichia coli. Different decay rates were observed for SmrC15 and SmrC16 of S. meliloti upon rifampicin treatment, revealing posttranscriptional regulation during growth. The use of a Deltahfq mutant showed that Hfq protects full-length SmrC16 from degradation and stabilises its specific degradation products.

Zusammenfassung

BACKGROUND:In bacteria, non-coding RNAs (ncRNA) are crucial regulators of gene expression, controlling various stress responses, virulence, and motility. Previous work revealed a relatively high number of ncRNAs in some marine cyanobacteria. However, for efficient genetic and biochemical analysis it would be desirable to identify a set of ncRNA candidate genes in model cyanobacteria that are easy to manipulate and for which extended mutant, transcriptomic and proteomic data sets are available.RESULTS:Here we have used comparative genome analysis for the biocomputational prediction of ncRNA genes and other sequence/structure-conserved elements in intergenic regions of the three unicellular model cyanobacteria Synechocystis PCC6803, Synechococcus elongatus PCC6301 and Thermosynechococcus elongatus BP1 plus the toxic Microcystis aeruginosa NIES843. The unfiltered numbers of predicted elements in these strains is 383, 168, 168, and 809, respectively, combined into 443 sequence clusters, whereas the numbers of individual elements with high support are 94, 56, 64, and 406, respectively. Removing also transposon-associated repeats, finally 78, 53, 42 and 168 sequences, respectively, are left belonging to 109 different clusters in the data set. Experimental analysis of selected ncRNA candidates in Synechocystis PCC6803 validated new ncRNAs originating from the fabF-hoxH and apcC-prmA intergenic spacers and three highly expressed ncRNAs belonging to the Yfr2 family of ncRNAs. Yfr2a promoter-luxAB fusions confirmed a very strong activity of this promoter and indicated a stimulation of expression if the cultures were exposed to elevated light intensities.CONCLUSION:Comparison to entries in Rfam and experimental testing of selected ncRNA candidates in Synechocystis PCC6803 indicate a high reliability of the current prediction, despite some contamination by the high number of repetitive sequences in some of these species. In particular, we identified in the four species altogether 8 new ncRNA homologs belonging to the Yfr2 family of ncRNAs. Modelling of RNA secondary structures indicated two conserved single-stranded sequence motifs that might be involved in RNA-protein interactions or in the recognition of target RNAs. Since our analysis has been restricted to find ncRNA candidates with a reasonable high degree of conservation among these four cyanobacteria, there might be many more, requiring direct experimental approaches for their identification.

Zusammenfassung

In recent years, sRNAs (small non-coding RNAs) have been found to be abundant in eukaryotes and bacteria and have been recognized as a novel class of gene expression regulators. In contrast, much less is known about sRNAs in archaea, except for snoRNAs (small nucleolar RNAs) that are involved in the modification of bases in stable RNAs. Therefore bioinformatic and experimental RNomics approaches were undertaken to search for the presence of sRNAs in the model archaeon Haloferax volcanii, resulting in more than 150 putative sRNA genes being identified. Northern blot analyses were used to study (differential) expression of sRNA genes. Several chromosomal deletion mutants of sRNA genes were generated and compared with the wild-type. It turned out that two sRNAs are essential for growth at low salt concentrations and high temperatures respectively, and one is involved in the regulation of carbon metabolism. Taken together, it could be shown that sRNAs are as abundant in H. volcanii as they are in well-studied bacterial species and that they fulfil important biological roles under specific conditions.

Zusammenfassung

Prochlorococcus is the most abundant phototroph in the vast, nutrient-poor areas of the ocean. It plays an important role in the ocean carbon cycle, and is a key component of the base of the food web. All cells share a core set of about 1,200 genes, augmented with a variable number of “flexible” genes. Many of the latter are located in genomic islands—hypervariable regions of the genome that encode functions important in differentiating the niches of “ecotypes.” Of major interest is how cells with such a small genome regulate cellular processes, as they lack many of the regulatory proteins commonly found in bacteria. We show here that contrary to the regulatory proteins, ncRNAs are present at levels typical of bacteria, revealing that they might have a disproportional regulatory role in Prochlorococcus—likely an adaptation to the extremely low-nutrient conditions of the open oceans, combined with the constraints of a small genome. Some of the ncRNAs were differentially expressed under stress conditions, and a high number of them were found to be associated with genomic islands, suggesting functional links between these RNAs and the response of Prochlorococcus to particular environmental challenges.

Zusammenfassung

BACKGROUND:MicroRNAs (miRNAs) are regulatory RNA molecules that are specified by their mode of action, the structure of primary transcripts, and their typical size of 20-24 nucleotides. Frequently, not only single miRNAs but whole families of closely related miRNAs have been found in animals and plants. Some families are widely conserved among different plant taxa. Hence, it is evident that these conserved miRNAs are of ancient origin and indicate essential functions that have been preserved over long evolutionary time scales. In contrast, other miRNAs seem to be species-specific and consequently must possess very distinct functions. Thus, the analysis of an early-branching species provides a window into the early evolution of fundamental regulatory processes in plants.RESULTS:Based on a combined experimental-computational approach, we report on the identification of 48 novel miRNAs and their putative targets in the moss Physcomitrella patens. From these, 18 miRNAs and two targets were verified in independent experiments. As a result of our study, the number of known miRNAs in Physcomitrella has been raised to 78. Functional assignments to mRNAs targeted by these miRNAs revealed a bias towards genes that are involved in regulation, cell wall biosynthesis and defense. Eight miRNAs were detected with different expression in protonema and gametophore tissue. The miRNAs 1-50 and 2-51 are located on a shared precursor that are separated by only one nucleotide and become processed in a tissue-specific way.CONCLUSION:Our data provide evidence for a surprisingly diverse and complex miRNA population in Physcomitrella. Thus, the number and function of miRNAs must have significantly expanded during the evolution of early land plants. As we have described here within, the coupled maturation of two miRNAs from a shared precursor has not been previously identified in plants.

Zusammenfassung

Different forms of the 6S non-coding RNA (ncRNA) exist in enterobacteria and in B. subtilis but there is only limited information about this RNA from other groups of bacteria. Prochlorococcus is an oceanic, ecologically important, cyanobacterium. It possesses the most streamlined genome within the cyanobacterial phylum, lacking many regulatory proteins and mechanisms well-known from other bacteria. Here we show the accumulation of two distinct types of 6S RNA in Prochlorococcus MED4. One of these RNAs is transcribed from a specific promoter located 23 nucleotides downstream the terminal codon of the purK gene, whereas the longer transcript is produced by processing from a purK-6S RNA precursor. The expression of both 6S transcripts is under diel control, reaching maxima during the day and minima coinciding with the S- and G2-like phases which are typical for synchronized cultures of this prokaryote. Based on data from four closely related Prochlorococcus strains and 11 environmental sequences from the Sargasso Sea, a previously unknown structural element is predicted within the 6S RNA 5' domain by comparative computational analysis. The divergent expression in synchronized cultures and unusual structural domains that were detected based on metagenomic data sets indicate that 6S RNA is an extremely important global regulator in these marine cyanobacteria.

Zusammenfassung

BACKGROUND:Non-coding RNAs (ncRNA) are regulators of gene expression in all domains of life. They control growth and differentiation, virulence, motility and various stress responses. The identification of ncRNAs can be a tedious process due to the heterogeneous nature of this molecule class and the missing sequence similarity of orthologs, even among closely related species. The small ncRNA Yfr1 has previously been found in the Prochlorococcus/Synechococcus group of marine cyanobacteria.RESULTS:Here we show that screening available genome sequences based on an RNA motif and followed by experimental analysis works successfully in detecting this RNA in all lineages of cyanobacteria. Yfr1 is an abundant ncRNA between 54 and 69 nt in size that is ubiquitous for cyanobacteria except for two low light-adapted strains of Prochlorococcus, MIT 9211 and SS120, in which it must have been lost secondarily. Yfr1 consists of two predicted stem-loop elements separated by an unpaired sequence of 16-20 nucleotides containing the ultraconserved undecanucleotide 5'-ACUCCUCACAC-3'.CONCLUSION:Starting with an ncRNA previously found in a narrow group of cyanobacteria only, we show here the highly specific and sensitive identification of its homologs within all lineages of cyanobacteria, whereas it was not detected within the genome sequences of E. coli and of 7 other eubacteria belonging to the alpha-proteobacteria, chlorobiaceae and spirochaete. The integration of RNA motif prediction into computational pipelines for the detection of ncRNAs in bacteria appears as a promising step to improve the quality of such predictions.

Zusammenfassung

Micro-RNAs (miRNAs) are one class of small non-coding RNAs that have important regulatory roles in higher plants. Much less is known about their prevalence and function in lower land plants. Previously we cloned 100 non-structural small RNAs from the moss Physcomitrella patens but could annotate only 11 as miRNAs. To identify additional moss miRNAs among cloned small RNAs we have analyzed their genomic sequences for a characteristic miRNA precursor-like structure. This analysis revealed 19 new moss miRNAs that are predicted to be encoded by 22 putative foldbacks. Northern blot analysis confirmed the expression of 14 new miRNA representatives. Half of these were gametophore specific, the rest were detected at low levels in the protonema. We predicted 12 genes as targets of nine new miRNAs. Three of these show homology to transcription factors and the others appear to play roles in diverse physiological processes including light and cytokine signaling, which have not to date been shown to be regulated by a miRNA in flowering plants. Four target genes, which show homology to ATN1-like protein kinase, NAC transcription factors and a cytokinin receptor, have been validated by miRNA-mediated mRNA cleavage. In addition, our analysis revealed that seven small RNAs represent miRNA* and three represent intermediates of pre-miRNA processing, providing evidence for specific DICER-like cleavage steps during miRNA biogenesis in moss. Our findings suggest that miRNAs are common in mosses and set the stage for the elucidation of their varied biological functions.

Zusammenfassung

BACKGROUND:Soon after the first algorithms for RNA folding became available, it was recognised that the prediction of only one energetically optimal structure is insufficient to achieve reliable results. An in-depth analysis of the folding space as a whole appeared necessary to deduce the structural properties of a given RNA molecule reliably. Folding space analysis comprises various methods such as suboptimal folding, computation of base pair probabilities, sampling procedures and abstract shape analysis. Common to many approaches is the idea of partitioning the folding space into classes of structures, for which certain properties can be derived.RESULTS:In this paper we extend the approach of abstract shape analysis. We show how to compute the accumulated probabilities of all structures that share the same shape. While this implies a complete (non-heuristic) analysis of the folding space, the computational effort depends only on the size of the shape space, which is much smaller. This approach has been integrated into the tool RNAshapes, and we apply it to various RNAs.CONCLUSION:Analyses of conformational switches show the existence of two shapes with probabilities approximately 23 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaiabikdaYaqaaiabiodaZaaaaaa@2EA2@ vs. 13 MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacH8akY=wiFfYdH8Gipec8Eeeu0xXdbba9frFj0=OqFfea0dXdd9vqai=hGuQ8kuc9pgc9s8qqaq=dirpe0xb9q8qiLsFr0=vr0=vr0dc8meaabaqaciaacaGaaeqabaqabeGadaaakeaadaWcaaqaaiabigdaXaqaaiabiodaZaaaaaa@2EA0@, whereas the analysis of a microRNA precursor reveals one shape with a probability near to 1.0. Furthermore, it is shown that a shape can outperform an energetically more favourable one by achieving a higher probability. From these results, and the fact that we use a complete and exact analysis of the folding space, we conclude that this approach opens up new and promising routes for investigating and understanding RNA secondary structure.

Zusammenfassung

Computational analysis of RNA secondary structure is a classical field of biosequence analysis, which has recently gained momentum due to the manyfold regulatory functions of RNA that have become apparent. We present five recent computational approaches that address the problems of synoptic folding space analysis, pseudoknot prediction, structure alignment, comparative structure prediction, and miRNA target prediction. All these programs are in current use and are available via the Bielefeld Bioinformatics Server at http://bibiserv.techfak.uni-bielefeld.de.

Zusammenfassung

Summary: We introduce RNAshapes, a new software package that integrates three RNA analysis tools based on the abstract shapes approach: the analysis of shape representatives, the calculation of shape probabilities and the consensus shapes approach. This new package is completely reimplemented in C and outruns the original implementations significantly in runtime and memory requirements. Additionally, we added a number of useful features like suboptimal folding with correct dangling energies, structure graph output, shape matching and a sliding window approach. Availability: RNAshapes is freely available at http://bibiserv.techfak.uni-bielefeld.de/rnashapes/ as C source code, and as compiled binaries for the most common computer architectures. For Microsoft Windows, we also offer a graphical user interface with convenient access to the complete functionality of the package. Contact: psteffen@techfak.uni-bielefeld.de

Zusammenfassung

Recent years have seen a rapid increase in the need for highperformance computing. These demands come from disciplines such as particle physics traditionally relying on High Performance Computing (HPC) but lately also from the various branches of life science that have matured into quantitative disciplines. The classical infrastructure of university computer centres results to be unsuited to cope with the new requirements for a multitude of reasons. Here we discuss the causes of this failure and present a solution developed at the University of Freiburg in a collaborative effort of several faculties. We demonstrate that using state of the art grid computing technology the problem can now be addressed in a bottom-up approach. The organizational, technical, and financial components of our framework, the Black Forest Grid Initiative (BFG) are described and results of its implementation are presented. In the process, a number of new questions have emerged which the next phase of our project needs to address.

Zusammenfassung

The knowledge about classes of non-coding RNAs (ncRNAs) is growing very fast and it is mainly the structure which is the common characteristic property shared by members of the same class. For correct characterization of such classes it is therefore of great importance to analyse the structural features in great detail. In this manuscript I present RNAlishapes which combines various secondary structure analysis methods, such as suboptimal folding and shape abstraction, with a comparative approach known as RNA alignment folding. RNAlishapes makes use of an extended thermodynamic model and covariance scoring, which allows to reward covariation of paired bases. Applying the algorithm to a set of bacterial trp-operon leaders using shape abstraction it was able to identify the two alternating conformations of this attenuator. Besides providing in-depth analysis methods for aligned RNAs, the tool also shows a fairly well prediction accuracy. Therefore, RNAlishapes provides the community with a powerful tool for structural analysis of classes of RNAs and is also a reasonable method for consensus structure prediction based on sequence alignments. RNAlishapes is available for online use and download at http://rna.cyanolab.de.

Zusammenfassung

The analysis of RNA secondary structure has become more and more important throughout the last decades after it was recognised that RNA does not only serve as a passive messenger (mRNA), but also as a functional compound of the cell. Furthermore, it was elucidated that mainly the structure rather than the sequence determines the function of such non-protein-coding RNA. This means that two RNA molecules which have low sequence similarity but high structure similarity are likely to have a similar function. The prediction of RNA secondary structure is based on parameters that have been measured in vitro. This results in rather static parameters, that do not incorporate the dynamic change of environment occurring in living organisms. Nevertheless, the use of these parameters, that are summarised in the energy model, gave valuable results, especially for short sequences. Several refinements throughout the years improved the predictions, but still the calculated optimal structure is not guaranteed to correspond to the native one. In this case, and due to the fact that the native structure is feasible under the energy model, it is common practice to additionally calculate suboptimal structures and incorporate these in the study. The set of all suboptimal structures is referred to as the structure space, which actually holds the information needed to answer questions such as: Is the optimal structure also the native one? Are there more than one structure an RNA molecule can adopt? How well-defined is the optimal structure? Major problems in the analysis of the structure space are its size and its shape. The number of suboptimal structures is exponential in the sequence length, which means that for sequences of moderate length the size quickly exceeds several billion. Besides the size, the appearance of the structure space complicates its study. The structure space can be imagined as a rough landscape with valleys, holding local optimal structures, separated by mountains and saddles. This landscape is not smooth but cliffy and complex, which prevents the development of a practical and still intuitive visualisation. In general, the intention of structure space analysis is not its visualisation, but its complexity also hampers approaches to derive specific features hidden in the structure space. Despite these problems, several tools exist that analyse the complete structure space or at least a part of it to answer the aforementioned questions. Among these are MFOLD which produces a subset of all possible structures according to a threshold of structural similarity, SFOLD which samples the structures in a probabilistic fashion and provides a method to identify alternating structures, RNAsubopt to produce all suboptimal structures within a given energy threshold, barriers to identify valleys, mountains and saddles of the structure landscape, and others. My contribution to this area of research is twofold: First, I present paRNAss (prediction of alternating RNA secondary structures) which focuses on the detection of conformational switches and analyses the structure space based on pairwise comparisons. paRNAss has been available since 1997 and I could improve its predictive power as well as its speed which made possible a systematic evaluation. During this evaluation it turned out that paRNAss can even be used to identify more than two competing structures and hence get a deeper insight into the structure space. The second tool I introduce is RNAshapes which facilitates different kinds of analyses. The algorithm makes use of abstract representations of the secondary structure to compute only those that are morphologically dissimilar, i.e. are composed of different structural elements. Structures being morphologically similar are pooled in a class of structures and each class is represented by its best member. The list of these representatives gives a general overview of what is there in the structure space. In addition to this, I introduce an algorithm to compute probabilities of the aforementioned classes of structures. This gives hints to properties such as alternating secondary structures (two classes with similar probabilities) and structural well-definedness (one class with very high probability).

Zusammenfassung

The function of a non-protein-coding RNA is often determined by its structure. Since experimental determination of RNA structure is time-consuming and expensive, its computational prediction is of great interest, and efficient solutions based on thermodynamic parameters are known. Frequently, however, the predicted minimum free energy structures are not the native ones, leading to the necessity of generating suboptimal solutions. While this can be accomplished by a number of programs, the user is often confronted with large outputs of similar structures, although he or she is interested in structures with more fundamental differences, or, in other words, with different abstract shapes. Here, we formalize the concept of abstract shapes and introduce their efficient computation. Each shape of an RNA molecule comprises a class of similar structures and has a representative structure of minimal free energy within the class. Shape analysis is implemented in the program RNAshapes. We applied RNAshapes to the prediction of optimal and suboptimal abstract shapes of several RNAs. For a given energy range, the number of shapes is considerably smaller than the number of structures, and in all cases, the native structures were among the top shape representatives. This demonstrates that the researcher can quickly focus on the structures of interest, without processing up to thousands of near-optimal solutions. We complement this study with a large-scale analysis of the growth behaviour of structure and shape spaces. RNAshapes is available for download and as an online version on the Bielefeld Bioinformatics Server.

Zusammenfassung

Motivation: There are various cases where the biological function of an RNA molecule involves a reversible change of conformation. paRNAss is a software approach to the prediction of such structural switching in RNA. It is based on three hypotheses about the secondary structure space of a switching RNA molecule that can be evaluated by RNA folding and structure comparison. In the positive case, the predicted structural switching must be verified experimentally. Results: After reviewing the strategy used in paRNAss, we present recent improvements on the algorithmic level of the approach, and the results of an evaluation procedure, comprising 1500 RNA sequences. It could be shown that the paRNAss approach performs well on known examples for conformational switching in RNA. The overall number of positive predictions was small, whereas for human 3′ UTRs, representing regulatory important regions, it was substantially higher than for arbitrary natural and random sequences. Availability: paRNAss is available as a Web service at http://bibiserv.techfak.uni-bielefeld.de/parnass Supplementary information: Detailed information on the analyses summarized in Table 1 can be found at http://bibiserv.techfak.uni-bielefeld.de/parnass/examples.html