The challenge of regulation in a minimal photoautotroph: non-coding RNAs in Prochlorococcus

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.

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.

Toward a systems-level understanding of gene regulatory, protein interaction, and metabolic networks in cyanobacteria

Cyanobacteria are essential primary producers in marine ecosystems, playing an important role in both carbon and nitrogen cycles. In the last decade, various genome sequencing and metagenomic projects have generated large amounts of genetic data for cyanobacteria. This wealth of data provides researchers with a new basis for the study of molecular adaptation, ecology and evolution of cyanobacteria, as well as for developing biotechnological applications. It also facilitates the use of multiplex techniques, i.e., expression profiling by high-throughput technologies such as microarrays, RNA-seq, and proteomics. However, exploration and analysis of these data is challenging, and often requires advanced computational methods. Also, they need to be integrated into our existing framework of knowledge to use them to draw reliable biological conclusions. Here, systems biology provides important tools. Especially, the construction and analysis of molecular networks has emerged as a powerful systems-level framework, with which to integrate such data, and to better understand biological relevant processes in these organisms. In this review, we provide an overview of the advances and experimental approaches undertaken using multiplex data from genomic, transcriptomic, proteomic, and metabolomic studies in cyanobacteria. Furthermore, we summarize currently available web-based tools dedicated to cyanobacteria, i.e., CyanoBase, CyanoEXpress, ProPortal, Cyanorak, CyanoBIKE, and CINPER. Finally, we present a case study for the freshwater model cyanobacteria, Synechocystis sp. PCC6803, to show the power of meta-analysis, and the potential to extrapolate acquired knowledge to the ecologically important marine cyanobacteria genus, Prochlorococcus.

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.

Comparative genome analysis of the closely related Synechocystis strains PCC 6714 and PCC 6803

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 Zn(2+) 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).

Comparative transcriptomics of two environmentally relevant cyanobacteria reveals unexpected transcriptome diversity

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.

Implications of streamlining theory for microbial ecology

Whether a small cell, a small genome or a minimal set of chemical reactions with self-replicating properties, simplicity is beguiling. As Leonardo da Vinci reportedly said, 'simplicity is the ultimate sophistication'. Two diverging views of simplicity have emerged in accounts of symbiotic and commensal bacteria and cosmopolitan free-living bacteria with small genomes. The small genomes of obligate insect endosymbionts have been attributed to genetic drift caused by small effective population sizes (Ne). In contrast, streamlining theory attributes small cells and genomes to selection for efficient use of nutrients in populations where Ne is large and nutrients limit growth. Regardless of the cause of genome reduction, lost coding potential eventually dictates loss of function. Consequences of reductive evolution in streamlined organisms include atypical patterns of prototrophy and the absence of common regulatory systems, which have been linked to difficulty in culturing these cells. Recent evidence from metagenomics suggests that streamlining is commonplace, may broadly explain the phenomenon of the uncultured microbial majority, and might also explain the highly interdependent (connected) behavior of many microbial ecosystems. Streamlining theory is belied by the observation that many successful bacteria are large cells with complex genomes. To fully appreciate streamlining, we must look to the life histories and adaptive strategies of cells, which impose minimum requirements for complexity that vary with niche.

Translational control of small heat shock genes in mesophilic and thermophilic cyanobacteria by RNA thermometers

Cyanobacteria constitute a heterogeneous phylum of oxygen-producing, photosynthetic prokaryotes. They are susceptible to various stress conditions like heat, salt, or light stress, all inducing the cyanobacterial heat shock response (HSR). Cyanobacterial small heat shock proteins (sHsps) are known to preserve thylakoid membrane integrity under stress conditions, thereby protecting the photosynthesis machinery. In Synechocystis sp PCC 6803, synthesis of the sHsp Hsp17 is regulated by an RNA thermometer (RNAT) in the 5′-untranslated region (5′-UTR) of the hsp17 mRNA. RNATs are direct temperature sensors that control expression of many bacterial heat shock and virulence genes. They hinder translation at low temperatures by base pairing, thus blocking ribosome access to the mRNA.

 

To explore the temperature range in which RNATs act, we studied various RNAT candidates upstream of sHsp genes from mesophilic and thermophilic cyanobacteria. The mesophilic cyanobacteria Anabaena variabilis and Nostoc sp chromosomally encode two sHsps each. Reporter gene studies suggested RNAT-mediated post-transcriptional regulation of shsp expression in both organisms. Detailed structural analysis of the two A. variabilis candidates revealed two novel RNAT types. The first, avashort, regulates translation primarily by masking of the AUG translational start codon. The second, featuring an extended initial hairpin, thus named avalong, presumably makes use of complex tertiary interaction. The 5′-UTR of the small heat shock gene hspA in the thermophile Thermosynechococcus elongatus is predicted to adopt an extended secondary structure. Structure probing revealed that the ribosome binding site was blocked at temperatures below 55 °C. The results of this study demonstrate that cyanobacteria commonly use RNATs to control expression of their small heat shock genes.

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.

The transcriptome landscape of Prochlorococcus MED4 and the factors for stabilizing the core genome

BACKGROUND

Gene gain and loss frequently occurs in the cyanobacterium Prochlorococcus, a phototroph that numerically dominates tropical and subtropical open oceans. However, little is known about the stabilization of its core genome, which contains approximately 1250 genes, in the context of genome streamlining. Using Prochlorococcus MED4 as a model organism, we investigated the constraints on core genome stabilization using transcriptome profiling.

RESULTS

RNA-Seq technique was used to obtain the transcriptome map of Prochlorococcus MED4, including operons, untranslated regions, non-coding RNAs, and novel genes. Genome-wide expression profiles revealed that three factors contribute to core genome stabilization. First, a negative correlation between gene expression levels and protein evolutionary rates was observed. Highly expressed genes were overrepresented in the core genome but not in the flexible genome. Gene necessity was determined as a second powerful constraint on genome evolution through functional enrichment analysis. Third, quick mRNA turnover may increase corresponding proteins' fidelity among genes that were abundantly expressed. Together, these factors influence core genome stabilization during MED4 genome evolution.

CONCLUSIONS

Gene expression, gene necessity, and mRNA turnover contribute to core genome maintenance during cyanobacterium Prochlorococcus genus evolution.

Proteomic and transcriptomic analyses of "Candidatus Pelagibacter ubique" describe the first PII-independent response to nitrogen limitation in a free-living Alphaproteobacterium

Nitrogen is one of the major nutrients limiting microbial productivity in the ocean, and as a result, most marine microorganisms have evolved systems for responding to nitrogen stress. The highly abundant alphaproteobacterium "Candidatus Pelagibacter ubique," a cultured member of the order Pelagibacterales (SAR11), lacks the canonical GlnB, GlnD, GlnK, and NtrB/NtrC genes for regulating nitrogen assimilation, raising questions about how these organisms respond to nitrogen limitation. A survey of 266 Alphaproteobacteria genomes found these five regulatory genes nearly universally conserved, absent only in intracellular parasites and members of the order Pelagibacterales, including "Ca. Pelagibacter ubique." Global differences in mRNA and protein expression between nitrogen-limited and nitrogen-replete cultures were measured to identify nitrogen stress responses in "Ca. Pelagibacter ubique" strain HTCC1062. Transporters for ammonium (AmtB), taurine (TauA), amino acids (YhdW), and opines (OccT) were all elevated in nitrogen-limited cells, indicating that they devote increased resources to the assimilation of nitrogenous organic compounds. Enzymes for assimilating amine into glutamine (GlnA), glutamate (GltBD), and glycine (AspC) were similarly upregulated. Differential regulation of the transcriptional regulator NtrX in the two-component signaling system NtrY/NtrX was also observed, implicating it in control of the nitrogen starvation response. Comparisons of the transcriptome and proteome supported previous observations of uncoupling between transcription and translation in nutrient-deprived "Ca. Pelagibacter ubique" cells. Overall, these data reveal a streamlined, PII-independent response to nitrogen stress in "Ca. Pelagibacter ubique," and likely other Pelagibacterales, and show that they respond to nitrogen stress by allocating more resources to the assimilation of nitrogen-rich organic compounds.

IMPORTANCE

Pelagibacterales are extraordinarily abundant and play a pivotal role in marine geochemical cycles, as one of the major recyclers of labile dissolved organic matter. They are also models for understanding how streamlining selection can reshape chemoheterotroph metabolism. Streamlining and its broad importance to environmental microbiology are emerging slowly from studies that reveal the complete genomes of uncultured organisms. Here, we report another remarkable example of streamlined metabolism in Pelagibacterales, this time in systems that control nitrogen assimilation. Pelagibacterales are major contributors to metatranscriptomes and metaproteomes from ocean systems, where patterns of gene expression are used to gain insight into ocean conditions and geochemical cycles. The data presented here supply background that is essential to interpreting data from field studies.

RNase III-dependent down-regulation of ftsH by an artificial internal sense RNA in Anabaena sp. PCC 7120

RNase III is a group of dsRNA-specific ribonucleases that play important roles in RNA processing and metabolism. Alr0280 and All4107 in Anabaena sp. PCC7120 are highly similar to RNase III enzymes. In vitro, recombinant Alr0280 showed RNase III activity. In the same cyanobacterium, the expression of ftsH (FtsH protease) could be suppressed by overexpression of an artificial sense RNA (aaftsH) that was complementary to aftsH, an internal antisense RNA. The aaftsH interference was abolished by inactivation of alr0280, the RNase III-encoding gene, and restored by complementation of the mutant. A cyanobacterial homolog to hen1, an RNA methyltransferase gene, may also be required for the aaftsH interference. This is the first report of RNase III-dependent sense RNA interference in cyanobacteria, and the underlying mechanism remains to be elucidated.

Toxin-antitoxin systems on the large defense plasmid pSYSA of Synechocystis sp. PCC 6803

Bacterial toxin-antitoxin (TA) systems are genetic elements, which are encoded by plasmid as well as chromosomal loci and mediate plasmid and genomic island maintenance through post-segregational killing mechanisms. TA systems exist in surprisingly high numbers in all prokaryotes, but cyanobacterial TA systems have been only very poorly experimentally characterized so far. Cyanobacteria are the only prokaryotes that perform oxygenic photosynthesis. As such, cyanobacteria are of high ecological importance and are considered promising for the production of biofuels. Here, we present the molecular characterization of the sll7003/ssl7004 TA system encoded on plasmid pSYSA of the model cyanobacterium Synechocystis sp. PCC 6803 as involving a Mg(2+)-dependent RNA endonuclease activity targeting single-stranded RNA regions and demonstrate the functionality of four more TA systems encoded on this 100,749-bp plasmid. Furthermore, one additional type I, one additional type II, and three freestanding TA system components are predicted on pSYSA, all of which appear active judged by their expression. By harboring at least seven simultaneously active TA systems, pSYSA appears as the plasmid most strongly selected for among all plasmids studied in this respect thus far. These results point to a high biological relevance of pSYSA, whose coding capacity is 75% devoted to three distinct clustered regularly interspaced short palindromic repeats (CRISPR) systems mediating antiviral defense.

The search for new chlorophyll-binding proteins in the cyanobacterium Synechocystis sp. PCC 6803

Light harvesting provides a major challenge in the production of biofuels from microorganisms; while sunlight provides the energy necessary for biomass/biofuel production, at the same time it damages the cells. The genome of Synechocystis sp. PCC 6803 was searched for open reading frames that might code for yet unidentified chlorophyll-binding proteins with low molecular mass that could be involved in stress-adaptation. Amongst 9167 hypothetical ORFs corresponding to potential polypeptides of 100 amino acids or less, two were identified that had the potential to be pigment-binding, because they (i) encoded a potential transmembrane region, (ii) showed sequence similarity with known chlorophyll-binding domains, (iii) were conserved in other cyanobacterial species, and (iv) their codon adaptation index indicated significant translation probability. The two ORFs were located complementary (antisense) and internal to the ferrochelatase (hemH) and the pyruvate dehydrogenase (pdh) genes and therefore were named a-fch and a-pdh, respectively. Transcription of both genes was confirmed; however, no translated proteins could be detected immunologically. Whereas mutations within a-pdh or a-fch did not lead to any obvious phenotype, it is clear that transcripts and proteins over and above the currently known set may play a role in defining the physiology of cyanobacteria and other organisms.

Synechococcus sp. Strain PCC 7002 Transcriptome: Acclimation to Temperature, Salinity, Oxidative Stress, and Mixotrophic Growth Conditions

Synechococcus sp. strain PCC 7002 is a unicellular, euryhaline cyanobacterium. It is a model organism for studies of cyanobacterial metabolism and has great potential for biotechnological applications. It exhibits an exceptional tolerance of high-light irradiation and shows very rapid growth. The habitats from which this and closely related strains were isolated are subject to changes in several environmental factors, including light, nutrient supply, temperature, and salinity. In this study global transcriptome profiling via RNAseq has been used to perform a comparative and integrated study of global changes in cells grown at different temperatures, at different salinities, and under mixotrophic conditions, when a metabolizable organic carbon source was present. Furthermore, the transcriptomes were investigated for cells that were subjected to a heat shock and that were exposed to oxidative stress. Lower growth temperatures caused relatively minor changes of the transcriptome; the most prominent changes affected fatty acid desaturases. A heat shock caused severe changes of the transcriptome pattern; transcripts for genes associated with major metabolic pathways declined and those for different chaperones increased dramatically. Oxidative stress, however, left the transcript pattern almost unaffected. When grown at high salinity, Synechococcus sp. PCC 7002 had increased expression of genes involved in compatible solute biosynthesis and showed increased mRNA levels for several genes involved in electron transport. Transcripts of two adjacent genes dramatically increased upon growth at high salinity; the respective proteins are putatively involved in coping with oxidative stress and in triggering ion channels. Only minor changes were observed when cells were grown at low salinity or when the growth medium was supplemented with glycerol. However, the transcriptome data suggest that cells must acclimate to excess reducing equivalents when a reduced C-source is present.

The antisense RNA As1_flv4 in the Cyanobacterium Synechocystis sp. PCC 6803 prevents premature expression of the flv4-2 operon upon shift in inorganic carbon supply

The functional relevance of natural cis-antisense transcripts is mostly unknown. Here we have characterized the association of three antisense RNAs and one intergenically encoded noncoding RNA with an operon that plays a crucial role in photoprotection of photosystem II under low carbon conditions in the cyanobacterium Synechocystis sp. PCC 6803. Cyanobacteria show strong gene expression dynamics in response to a shift of cells from high carbon to low levels of inorganic carbon (C_i), but the regulatory mechanisms are poorly understood. Among the most up-regulated genes in Synechocystis are flv4, sll0218, and flv2, which are organized in the flv4-2 operon. The flavodiiron proteins encoded by this operon open up an alternative electron transfer route, likely starting from the Q_B site in photosystem II, under photooxidative stress conditions. Our expression analysis of cells shifted from high carbon to low carbon demonstrated an inversely correlated transcript accumulation of the flv4-2 operon mRNA and one antisense RNA to flv4, designated as As1_flv4. Overexpression of As1_flv4 led to a decrease in flv4-2 mRNA. The promoter activity of as1_flv4 was transiently stimulated by C_i limitation and negatively regulated by the AbrB-like transcription regulator Sll0822, whereas the flv4-2 operon was positively regulated by the transcription factor NdhR. The results indicate that the tightly regulated antisense RNA As1_flv4 establishes a transient threshold for flv4-2 expression in the early phase after a change in C_i conditions. Thus, it prevents unfavorable synthesis of the proteins from the flv4-2 operon.

Acclimation of the Global Transcriptome of the Cyanobacterium Synechococcus sp. Strain PCC 7002 to Nutrient Limitations and Different Nitrogen Sources

The unicellular, euryhaline cyanobacterium Synechococcus sp. strain PCC 7002 is a model organism for laboratory-based studies of cyanobacterial metabolism and is a potential platform for biotechnological applications. Two of its most notable properties are its exceptional tolerance of high-light intensity and very rapid growth under optimal conditions. In this study, transcription profiling by RNAseq has been used to perform an integrated study of global changes in transcript levels in cells subjected to limitation for the major nutrients CO₂, nitrogen, sulfate, phosphate, and iron. Transcriptional patterns for cells grown on nitrate, ammonia, and urea were also studied. Nutrient limitation caused strong decreases of transcript levels of the genes encoding major metabolic pathways, especially for components of the photosynthetic apparatus, CO₂ fixation, and protein biosynthesis. Uptake mechanisms for the respective nutrients were strongly up-regulated. The transcription data further suggest that major changes in the composition of the NADH dehydrogenase complex occur upon nutrient limitation. Transcripts for flavoproteins increased strongly when CO₂ was limiting. Genes involved in protection from oxidative stress generally showed high, constitutive transcript levels, which possibly explains the high-light tolerance of this organism. The transcriptomes of cells grown with ammonia or urea as nitrogen source showed increased transcript levels for components of the CO₂ fixation machinery compared to cells grown with nitrate, but in general transcription differences in cells grown on different N-sources exhibited surprisingly minor differences.

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.

Unexpected Diversity of Chloroplast Noncoding RNAs as Revealed by Deep Sequencing of the Arabidopsis Transcriptome

Noncoding RNAs (ncRNA) are widely expressed in both prokaryotes and eukaryotes. Eukaryotic ncRNAs are commonly micro- and small-interfering RNAs (18-25 nt) involved in posttranscriptional gene silencing, whereas prokaryotic ncRNAs vary in size and are involved in various aspects of gene regulation. Given the prokaryotic origin of organelles, the presence of ncRNAs might be expected; however, the full spectrum of organellar ncRNAs has not been determined systematically. Here, strand-specific RNA-Seq analysis was used to identify 107 candidate ncRNAs from Arabidopsis thaliana chloroplasts, primarily encoded opposite protein-coding and tRNA genes. Forty-eight ncRNAs were shown to accumulate by RNA gel blot as discrete transcripts in wild-type (WT) plants and/or the pnp1-1 mutant, which lacks the chloroplast ribonuclease polynucleotide phosphorylase (cpPNPase). Ninety-eight percent of the ncRNAs detected by RNA gel blot had different transcript patterns between WT and pnp1-1, suggesting cpPNPase has a significant role in chloroplast ncRNA biogenesis and accumulation. Analysis of materials deficient for other major chloroplast ribonucleases, RNase R, RNase E, and RNase J, showed differential effects on ncRNA accumulation and/or form, suggesting specificity in RNase-ncRNA interactions. 5' end mapping demonstrates that some ncRNAs are transcribed from dedicated promoters, whereas others result from transcriptional read-through. Finally, correlations between accumulation of some ncRNAs and the symmetrically transcribed sense RNA are consistent with a role in RNA stability. Overall, our data suggest that this extensive population of ncRNAs has the potential to underpin a previously underappreciated regulatory mode in the chloroplast.

Chloroplast RNase J compensates for inefficient transcription termination by removal of antisense RNA

Ribonuclease J is an essential enzyme, and the Bacillus subtilis ortholog possesses both endoribonuclease and 5' → 3' exoribonuclease activities. Chloroplasts also contain RNase J, which has been postulated to participate, as both an exo- and endonuclease, in the maturation of polycistronic mRNAs. Here we have examined recombinant Arabidopsis RNase J and found both 5' → 3' exoribonuclease and endonucleolytic activities. Virus-induced gene silencing was used to reduce RNase J expression in Arabidopsis and Nicotiana benthamiana, leading to chlorosis but surprisingly few disruptions in the cleavage of polycistronic rRNA and mRNA precursors. In contrast, antisense RNAs accumulated massively, suggesting that the failure of chloroplast RNA polymerase to terminate effectively leads to extensive symmetric transcription products that are normally eliminated by RNase J. Mung bean nuclease digestion and polysome analysis revealed that this antisense RNA forms duplexes with sense strand transcripts and prevents their translation. We conclude that a major role of chloroplast RNase J is RNA surveillance to prevent overaccumulation of antisense RNA, which would otherwise exert deleterious effects on chloroplast gene expression.

Dynamics of transcriptional start site selection during nitrogen stress-induced cell differentiation in Anabaena sp. PCC7120

The fixation of atmospheric N(2) by cyanobacteria is a major source of nitrogen in the biosphere. In Nostocales, such as Anabaena, this process is spatially separated from oxygenic photosynthesis and occurs in heterocysts. Upon nitrogen step-down, these specialized cells differentiate from vegetative cells in a process controlled by two major regulators: NtcA and HetR. However, the regulon controlled by these two factors is only partially defined, and several aspects of the differentiation process have remained enigmatic. Using differential RNA-seq, we experimentally define a genome-wide map of >10,000 transcriptional start sites (TSS) of Anabaena sp. PCC7120, a model organism for the study of prokaryotic cell differentiation and N(2) fixation. By analyzing the adaptation to nitrogen stress, our global TSS map provides insight into the dynamic changes that modify the transcriptional organization at a critical step of the differentiation process. We identify >900 TSS with minimum fold change in response to nitrogen deficiency of eight. From these TSS, at least 209 were under control of HetR, whereas at least 158 other TSS were potentially directly controlled by NtcA. Our analysis of the promoters activated during the switch to N(2) fixation adds hundreds of protein-coding genes and noncoding transcripts to the list of potentially involved factors. These data experimentally define the NtcA regulon and the DIF(+) motif, a palindrome at or close to position -35 that seems essential for heterocyst-specific expression of certain genes.

A tiling microarray for global analysis of chloroplast genome expression in cucumber and other plants

Plastids are small organelles equipped with their own genomes (plastomes). Although these organelles are involved in numerous plant metabolic pathways, current knowledge about the transcriptional activity of plastomes is limited. To solve this problem, we constructed a plastid tiling microarray (PlasTi-microarray) consisting of 1629 oligonucleotide probes. The oligonucleotides were designed based on the cucumber chloroplast genomic sequence and targeted both strands of the plastome in a non-contiguous arrangement. Up to 4 specific probes were designed for each gene/exon, and the intergenic regions were covered regularly, with 70-nt intervals. We also developed a protocol for direct chemical labeling and hybridization of as little as 2 micrograms of chloroplast RNA. We used this protocol for profiling the expression of the cucumber chloroplast plastome on the PlasTi-microarray. Owing to the high sequence similarity of plant plastomes, the newly constructed microarray can be used to study plants other than cucumber. Comparative hybridization of chloroplast transcriptomes from cucumber, Arabidopsis, tomato and spinach showed that the PlasTi-microarray is highly versatile.

cis-antisense RNA, another level of gene regulation in bacteria

A substantial amount of antisense transcription is a hallmark of gene expression in eukaryotes. However, antisense transcription was first demonstrated in bacteria almost 50 years ago. The transcriptomes of bacteria as different as Helicobacter pylori, Bacillus subtilis, Escherichia coli, Synechocystis sp. strain PCC6803, Mycoplasma pneumoniae, Sinorhizobium meliloti, Geobacter sulfurreducens, Vibrio cholerae, Chlamydia trachomatis, Pseudomonas syringae, and Staphylococcus aureus have now been reported to contain antisense RNA (asRNA) transcripts for a high percentage of genes. Bacterial asRNAs share functional similarities with trans-acting regulatory RNAs, but in addition, they use their own distinct mechanisms. Among their confirmed functional roles are transcription termination, codegradation, control of translation, transcriptional interference, and enhanced stability of their respective target transcripts. Here, we review recent publications indicating that asRNAs occur as frequently in simple unicellular bacteria as they do in higher organisms, and we provide a comprehensive overview of the experimentally confirmed characteristics of asRNA actions and intimately linked quantitative aspects. Emerging functional data suggest that asRNAs in bacteria mediate a plethora of effects and are involved in far more processes than were previously anticipated. Thus, the functional impact of asRNAs should be considered when developing new strategies against pathogenic bacteria and when optimizing bacterial strains for biotechnology.

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.

Transcriptome response of high- and low-light-adapted Prochlorococcus strains to changing iron availability

Prochlorococcus contributes significantly to ocean primary productivity. The link between primary productivity and iron in specific ocean regions is well established and iron limitation of Prochlorococcus cell division rates in these regions has been shown. However, the extent of ecotypic variation in iron metabolism among Prochlorococcus and the molecular basis for differences is not understood. Here, we examine the growth and transcriptional response of Prochlorococcus strains, MED4 and MIT9313, to changing iron concentrations. During steady state, MIT9313 sustains growth at an order-of-magnitude lower iron concentration than MED4. To explore this difference, we measured the whole-genome transcriptional response of each strain to abrupt iron starvation and rescue. Only four of the 1159 orthologs of MED4 and MIT9313 were differentially expressed in response to iron in both strains. However, in each strain, the expression of over a hundred additional genes changed, many of which are in labile genomic regions, suggesting a role for lateral gene transfer in establishing diversity of iron metabolism among Prochlorococcus. Furthermore, we found that MED4 lacks three genes near the iron-deficiency-induced gene (idiA) that are present and induced by iron stress in MIT9313. These genes are interesting targets for studying the adaptation of natural Prochlorococcus assemblages to local iron conditions as they show more diversity than other genomic regions in environmental metagenomic databases.

Antisense RNA protects mRNA from RNase E degradation by RNA-RNA duplex formation during phage infection

The ecologically important cyanobacterium Prochlorococcus possesses the smallest genome among oxyphototrophs, with a reduced suite of protein regulators and a disproportionately high number of regulatory RNAs. Many of these are asRNAs, raising the question whether they modulate gene expression through the protection of mRNA from RNase E degradation. To address this question, we produced recombinant RNase E from Prochlorococcus sp. MED4, which functions optimally at 12 mM Mg²⁺, pH 9 and 35°C. RNase E cleavage assays were performed with this recombinant protein to assess enzyme activity in the presence of single- or double-stranded RNA substrates. We found that extraordinarily long asRNAs of 3.5 and 7 kb protect a set of mRNAs from RNase E degradation that accumulate during phage infection. These asRNA–mRNA duplex formations mask single-stranded recognition sites of RNase E, leading to increased stability of the mRNAs. Such interactions directly modulate RNA stability and provide an explanation for enhanced transcript abundance of certain mRNAs during phage infection. Protection from RNase E-triggered RNA decay may constitute a hitherto unknown regulatory function of bacterial cis-asRNAs, impacting gene expression.

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.

Bacterial antisense RNAs: how many are there, and what are they doing?

Antisense RNAs encoded on the DNA strand opposite another gene have the potential to form extensive base-pairing interactions with the corresponding sense RNA. Unlike other smaller regulatory RNAs in bacteria, antisense RNAs range in size from tens to thousands of nucleotides. The numbers of antisense RNAs reported for different bacteria vary extensively, but hundreds have been suggested in some species. If all of these reported antisense RNAs are expressed at levels sufficient to regulate the genes encoded opposite them, antisense RNAs could significantly impact gene expression in bacteria. Here, we review the evidence for these RNA regulators and describe what is known about the functions and mechanisms of action for some of these RNAs. Important considerations for future research as well as potential applications are also discussed.

Overexpression of a natural chloroplast-encoded antisense RNA in tobacco destabilizes 5S rRNA and retards plant growth

BACKGROUND

The roles of non-coding RNAs in regulating gene expression have been extensively studied in both prokaryotes and eukaryotes, however few reports exist as to their roles in organellar gene regulation. Evidence for accumulation of natural antisense RNAs (asRNAs) in chloroplasts comes from the expressed sequence tag database and cDNA libraries, while functional data have been largely obtained from artificial asRNAs. In this study, we used Nicotiana tabacum to investigate the effect on sense strand transcripts of overexpressing a natural chloroplast asRNA, AS5, which is complementary to the region which encodes the 5S rRNA and tRNAArg.

RESULTS

AS5-overexpressing (AS5ox) plants obtained by chloroplast transformation exhibited slower growth and slightly pale green leaves. Analysis of AS5 transcripts revealed four distinct species in wild-type (WT) and AS5ox plants, and additional AS5ox-specific products. Of the corresponding sense strand transcripts, tRNAArg overaccumulated several-fold in transgenic plants whereas 5S rRNA was unaffected. However, run-on transcription showed that the 5S-trnR region was transcribed four-fold more in the AS5ox plants compared to WT, indicating that overexpression of AS5 was associated with decreased stability of 5S rRNA. In addition, polysome analysis of the transformants showed less 5S rRNA and rbcL mRNA associated with ribosomes.

CONCLUSIONS

Our results suggest that AS5 can modulate 5S rRNA levels, giving it the potential to affect Chloroplast translation and plant growth. More globally, overexpression of asRNAs via chloroplast transformation may be a useful strategy for defining their functions.

Ultraviolet stress delays chromosome replication in light/dark synchronized cells of the marine cyanobacterium Prochlorococcus marinus PCC9511

BACKGROUND

The marine cyanobacterium Prochlorococcus is very abundant in warm, nutrient-poor oceanic areas. The upper mixed layer of oceans is populated by high light-adapted Prochlorococcus ecotypes, which despite their tiny genome (approximately 1.7 Mb) seem to have developed efficient strategies to cope with stressful levels of photosynthetically active and ultraviolet (UV) radiation. At a molecular level, little is known yet about how such minimalist microorganisms manage to sustain high growth rates and avoid potentially detrimental, UV-induced mutations to their DNA. To address this question, we studied the cell cycle dynamics of P. marinus PCC9511 cells grown under high fluxes of visible light in the presence or absence of UV radiation. Near natural light-dark cycles of both light sources were obtained using a custom-designed illumination system (cyclostat). Expression patterns of key DNA synthesis and repair, cell division, and clock genes were analyzed in order to decipher molecular mechanisms of adaptation to UV radiation.

RESULTS

The cell cycle of P. marinus PCC9511 was strongly synchronized by the day-night cycle. The most conspicuous response of cells to UV radiation was a delay in chromosome replication, with a peak of DNA synthesis shifted about 2 h into the dark period. This delay was seemingly linked to a strong downregulation of genes governing DNA replication (dnaA) and cell division (ftsZ, sepF), whereas most genes involved in DNA repair (such as recA, phrA, uvrA, ruvC, umuC) were already activated under high visible light and their expression levels were only slightly affected by additional UV exposure.

CONCLUSIONS

Prochlorococcus cells modified the timing of the S phase in response to UV exposure, therefore reducing the risk that mutations would occur during this particularly sensitive stage of the cell cycle. We identified several possible explanations for the observed timeshift. Among these, the sharp decrease in transcript levels of the dnaA gene, encoding the DNA replication initiator protein, is sufficient by itself to explain this response, since DNA synthesis starts only when the cellular concentration of DnaA reaches a critical threshold. However, the observed response likely results from a more complex combination of UV-altered biological processes.

The gene ssl3076 encodes a protein mediating the salt-induced expression of ggpS for the biosynthesis of the compatible solute glucosylglycerol in Synechocystis sp. strain PCC 6803

Acclimation to high salt concentrations involves concerted changes in gene expression. For the majority of salt-regulated genes, the mechanism underlying the induction process is not known. The gene ggpS (sll1566), which encodes the glucosylglycerol-phosphate synthase responsible for the synthesis of the compatible solute glucosylglycerol (GG), is specifically induced by salt in the cyanobacterial model strain Synechocystis sp. strain PCC 6803. To identify mechanisms mediating this salt-specific gene regulation, the ggpS promoter was analyzed in more detail. 5' rapid amplification of cDNA ends (5'-RACE) experiments revealed that the adjacent open reading frame (ORF), which is annotated as unknown protein Ssl3076, overlaps with the transcriptional start site of the ggpS gene. Reporter gene expression analyses indicated an essential role for the intact ssl3076 gene in the salt-regulated transcription of a gfp reporter gene. Promoter fragments containing a mutated ssl3076 lost the salt regulation; similarly, a frameshift mutation in ssl3076 resulted in a high level of ggpS expression under low-salt conditions, thereby establishing this small ORF, named ggpR, as a negative regulator of ggpS. Interestingly, small ORFs were also found adjacent to ggpS genes in the genomes of other GG-accumulating cyanobacteria. These results suggest that the GgpR protein represses ggpS expression under low-salt conditions, whereas in salt-shocked and salt-acclimated cells a stress-proportional ggpS expression occurs, leading to GG accumulation.

Short RNA half-lives in the slow-growing marine cyanobacterium Prochlorococcus

BACKGROUND

RNA turnover plays an important role in the gene regulation of microorganisms and influences their speed of acclimation to environmental changes. We investigated whole-genome RNA stability of Prochlorococcus, a relatively slow-growing marine cyanobacterium doubling approximately once a day, which is extremely abundant in the oceans.

RESULTS

Using a combination of microarrays, quantitative RT-PCR and a new fitting method for determining RNA decay rates, we found a median half-life of 2.4 minutes and a median decay rate of 2.6 minutes for expressed genes - twofold faster than that reported for any organism. The shortest transcript half-life (33 seconds) was for a gene of unknown function, while some of the longest (approximately 18 minutes) were for genes with high transcript levels. Genes organized in operons displayed intriguing mRNA decay patterns, such as increased stability, and delayed onset of decay with greater distance from the transcriptional start site. The same phenomenon was observed on a single probe resolution for genes greater than 2 kb.

CONCLUSIONS

We hypothesize that the fast turnover relative to the slow generation time in Prochlorococcus may enable a swift response to environmental changes through rapid recycling of nucleotides, which could be advantageous in nutrient poor oceans. Our growing understanding of RNA half-lives will help us interpret the growing bank of metatranscriptomic studies of wild populations of Prochlorococcus. The surprisingly complex decay patterns of large transcripts reported here, and the method developed to describe them, will open new avenues for the investigation and understanding of RNA decay for all organisms.

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

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

A long antisense RNA in plant chloroplasts

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.

An antisense RNA in a lytic cyanophage links psbA to a gene encoding a homing endonuclease

Cyanophage genomes frequently possess the psbA gene, encoding the D1 polypeptide of photosystem II. This protein is believed to maintain host photosynthetic capacity during infection and enhance phage fitness under high-light conditions. Although the first documented cyanophage-encoded psbA gene contained a group I intron, this feature has not been widely reported since, despite a plethora of new sequences becoming available. In this study, we show that in cyanophage S-PM2, this intron is spliced during the entire infection cycle. Furthermore, we report the widespread occurrence of psbA introns in marine metagenomic libraries, and with psbA often adjacent to a homing endonuclease (HE). Bioinformatic analysis of the intergenic region between psbA and the adjacent HE gene F-CphI in S-PM2 showed the presence of an antisense RNA (asRNA) connecting these two separate genetic elements. The asRNA is co-regulated with psbA and F-CphI, suggesting its involvement with their expression. Analysis of scaffolds from global ocean survey datasets shows this asRNA to be commonly associated with the 3' end of cyanophage psbA genes, implying that this potential mechanism of regulating marine 'viral' photosynthesis is evolutionarily conserved. Although antisense transcription is commonly found in eukaryotic and increasingly also in prokaryotic organisms, there has been no indication for asRNAs in lytic phages so far. We propose that this asRNA also provides a means of preventing the formation of mobile group I introns within cyanophage psbA genes.

Heterocyst-specific transcription of NsiR1, a non-coding RNA encoded in a tandem array of direct repeats in cyanobacteria

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.

Comparative genomics reveals 104 candidate structured RNAs from bacteria, archaea, and their metagenomes

BACKGROUND

Structured noncoding RNAs perform many functions that are essential for protein synthesis, RNA processing, and gene regulation. Structured RNAs can be detected by comparative genomics, in which homologous sequences are identified and inspected for mutations that conserve RNA secondary structure.

RESULTS

By applying a comparative genomics-based approach to genome and metagenome sequences from bacteria and archaea, we identified 104 candidate structured RNAs and inferred putative functions for many of these. Twelve candidate metabolite-binding RNAs were identified, three of which were validated, including one reported herein that binds the coenzyme S-adenosylmethionine. Newly identified cis-regulatory RNAs are implicated in photosynthesis or nitrogen regulation in cyanobacteria, purine and one-carbon metabolism, stomach infection by Helicobacter, and many other physiological processes. A candidate riboswitch termed crcB is represented in both bacteria and archaea. Another RNA motif may control gene expression from 3'-untranslated regions of mRNAs, which is unusual for bacteria. Many noncoding RNAs that likely act in trans are also revealed, and several of the noncoding RNA candidates are found mostly or exclusively in metagenome DNA sequences.

CONCLUSIONS

This work greatly expands the variety of highly structured noncoding RNAs known to exist in bacteria and archaea and provides a starting point for biochemical and genetic studies needed to validate their biologic functions. Given the sustained rate of RNA discovery over several similar projects, we expect that far more structured RNAs remain to be discovered from bacterial and archaeal organisms.

Abundance of type I toxin-antitoxin systems in bacteria: searches for new candidates and discovery of novel families

Small, hydrophobic proteins whose synthesis is repressed by small RNAs (sRNAs), denoted type I toxin-antitoxin modules, were first discovered on plasmids where they regulate plasmid stability, but were subsequently found on a few bacterial chromosomes. We used exhaustive PSI-BLAST and TBLASTN searches across 774 bacterial genomes to identify homologs of known type I toxins. These searches substantially expanded the collection of predicted type I toxins, revealed homology of the Ldr and Fst toxins, and suggested that type I toxin-antitoxin loci are not spread by horizontal gene transfer. To discover novel type I toxin-antitoxin systems, we developed a set of search parameters based on characteristics of known loci including the presence of tandem repeats and clusters of charged and bulky amino acids at the C-termini of short proteins containing predicted transmembrane regions. We detected sRNAs for three predicted toxins from enterohemorrhagic Escherichia coli and Bacillus subtilis, and showed that two of the respective proteins indeed are toxic when overexpressed. We also demonstrated that the local free-energy minima of RNA folding can be used to detect the positions of the sRNA genes. Our results suggest that type I toxin-antitoxin modules are much more widely distributed among bacteria than previously appreciated.

Deep sequencing analysis of the Methanosarcina mazei Gö1 transcriptome in response to nitrogen availability

Methanosarcina mazei and related mesophilic archaea are the only organisms fermenting acetate, methylamines, and methanol to methane and carbon dioxide, contributing significantly to greenhouse gas production. The biochemistry of these metabolic processes is well studied, and genome sequences are available, yet little is known about the overall transcriptional organization and the noncoding regions representing 25% of the 4.01-Mb genome of M. mazei. We present a genome-wide analysis of transcription start sites (TSS) in M. mazei grown under different nitrogen availabilities. Pyrosequencing-based differential analysis of primary vs. processed 5' ends of transcripts discovered 876 TSS across the M. mazei genome. Unlike in other archaea, in which leaderless mRNAs are prevalent, the majority of the detected mRNAs in M. mazei carry long untranslated 5' regions. Our experimental data predict a total of 208 small RNA (sRNA) candidates, mostly from intergenic regions but also antisense to 5' and 3' regions of mRNAs. In addition, 40 new small mRNAs with ORFs of < or = 30 aa were identified, some of which might have dual functions as mRNA and regulatory sRNA. We confirmed differential expression of several sRNA genes in response to nitrogen availability. Inspection of their promoter regions revealed a unique conserved sequence motif associated with nitrogen-responsive regulation, which might serve as a regulator binding site upstream of the common IIB recognition element. Strikingly, several sRNAs antisense to mRNAs encoding transposases indicate nitrogen-dependent transposition events. This global TSS map in archaea will facilitate a better understanding of transcriptional and posttranscriptional control in the third domain of life.

Seed-based INTARNA prediction combined with GFP-reporter system identifies mRNA targets of the small RNA Yfr1

Motivation:Prochlorococcus possesses the smallest genome of all sequenced photoautotrophs. Although the number of regulatory proteins in the genome is very small, the relative number of small regulatory RNAs is comparable with that of other bacteria. The compact genome size of Prochlorococcus offers an ideal system to search for targets of small RNAs (sRNAs) and to refine existing target prediction algorithms.

Results: Target predictions for the cyanobacterial sRNA Yfr1 were carried out with INTARNA in Prochlorococcus MED4. The ultraconserved Yfr1 sequence motif was defined as the putative interaction seed. To study the impact of Yfr1 on its predicted mRNA targets, a reporter system based on green fluorescent protein (GFP) was applied. We show that Yfr1 inhibits the translation of two predicted targets. We used mutation analysis to confirm that Yfr1 directly regulates its targets by an antisense interaction sequestering the ribosome binding site, and to assess the importance of interaction site accessibility.

Contact:backofen@informatik.uni-freiburg.de; claudia.steglich@biologie.uni-freiburg.de

Supplementary information:Supplementary data are available at Bioinformatics online.

Evidence for a major role of antisense RNAs in cyanobacterial gene regulation

Information on the numbers and functions of naturally occurring antisense RNAs (asRNAs) in eubacteria has thus far remained incomplete. Here, we screened the model cyanobacterium Synechocystis sp. PCC 6803 for asRNAs using four different methods. In the final data set, the number of known noncoding RNAs rose from 6 earlier identified to 60 and of asRNAs from 1 to 73 (28 were verified using at least three methods). Among these, there are many asRNAs to housekeeping, regulatory or metabolic genes, as well as to genes encoding electron transport proteins. Transferring cultures to high light, carbon-limited conditions or darkness influenced the expression levels of several asRNAs, suggesting their functional relevance. Examples include the asRNA to rpl1, which accumulates in a light-dependent manner and may be required for processing the L11 r-operon and the SyR7 noncoding RNA, which is antisense to the murF 5′ UTR, possibly modulating murein biosynthesis. Extrapolated to the whole genome, ∼10% of all genes in Synechocystis are influenced by asRNAs. Thus, chromosomally encoded asRNAs may have an important function in eubacterial regulatory networks.

Experimental approaches for the discovery and characterization of regulatory small RNA

Following the pioneering screens for small regulatory RNAs (sRNAs) in Escherichia coli in 2001, sRNAs are now being identified in almost every branch of the eubacterial kingdom. Experimental strategies have become increasingly important for sRNA discovery, thanks to increased availability of tiling arrays and fast progress in the development of high-throughput cDNA sequencing (RNA-Seq). The new technologies also facilitate genome-wide discovery of potential target mRNAs by sRNA pulse-expression coupled to transcriptomics, and immunoprecipitation with RNA-binding proteins such as Hfq. Moreover, the staggering rate of new sRNAs demands mechanistic analysis of target regulation. We will also review the available toolbox for wet lab-based research, including in vivo and in vitro reporter systems, genetic methods and biochemical co-purification of sRNA interaction partners.

On the facultative requirement of the bacterial RNA chaperone, Hfq

The pleiotropic post-transcriptional regulator Hfq is an RNA chaperone that facilitates pairing interactions between small regulatory RNAs (sRNAs) and their mRNA targets in several bacteria. However, this classical pattern, derived from the Escherichia coli model, is not applicable to the whole bacterial kingdom. In this article we discuss the facultative requirement for Hfq for sRNA-mRNA duplex formation among bacteria and the specific features of the Hfq protein and RNA duplexes that might account for the dispensability or requirement of the chaperone. Apparent links between the need for Hfq, the GC content of bacterial genomes and the free energy of experimentally validated sRNA-mRNA pairing interactions are presented.

Metatranscriptomics reveals unique microbial small RNAs in the ocean's water column

Microbial gene expression in the environment has recently been assessed via pyrosequencing of total RNA extracted directly from natural microbial assemblages. Several such 'metatranscriptomic' studies have reported that many complementary DNA sequences shared no significant homology with known peptide sequences, and so might represent transcripts from uncharacterized proteins. Here we report that a large fraction of cDNA sequences detected in microbial metatranscriptomic data sets are comprised of well-known small RNAs (sRNAs), as well as new groups of previously unrecognized putative sRNAs (psRNAs). These psRNAs mapped specifically to intergenic regions of microbial genomes recovered from similar habitats, displayed characteristic conserved secondary structures and were frequently flanked by genes that indicated potential regulatory functions. Depth-dependent variation of psRNAs generally reflected known depth distributions of broad taxonomic groups, but fine-scale differences in the psRNAs within closely related populations indicated potential roles in niche adaptation. Genome-specific mapping of a subset of psRNAs derived from predominant planktonic species such as Pelagibacter revealed recently discovered as well as potentially new regulatory elements. Our analyses show that metatranscriptomic data sets can reveal new information about the diversity, taxonomic distribution and abundance of sRNAs in naturally occurring microbial communities, and indicate their involvement in environmentally relevant processes including carbon metabolism and nutrient acquisition.

Ecological genomics of marine picocyanobacteria

Marine picocyanobacteria of the genera Prochlorococcus and Synechococcus numerically dominate the picophytoplankton of the world ocean, making a key contribution to global primary production. Prochlorococcus was isolated around 20 years ago and is probably the most abundant photosynthetic organism on Earth. The genus comprises specific ecotypes which are phylogenetically distinct and differ markedly in their photophysiology, allowing growth over a broad range of light and nutrient conditions within the 45 degrees N to 40 degrees S latitudinal belt that they occupy. Synechococcus and Prochlorococcus are closely related, together forming a discrete picophytoplankton clade, but are distinguishable by their possession of dissimilar light-harvesting apparatuses and differences in cell size and elemental composition. Synechococcus strains have a ubiquitous oceanic distribution compared to that of Prochlorococcus strains and are characterized by phylogenetically discrete lineages with a wide range of pigmentation. In this review, we put our current knowledge of marine picocyanobacterial genomics into an environmental context and present previously unpublished genomic information arising from extensive genomic comparisons in order to provide insights into the adaptations of these marine microbes to their environment and how they are reflected at the genomic level.