<title>Introduction</title> <p>Cyanobacteria are a diverse group of photoautotrophic bacteria that occupy a broad range of habitats, including the oceans, lakes, and soil, and are also found as symbionts in many different types of organisms. <italic>Prochlorococcus</italic>, a member of the cyanobacterial lineage, often accounts for up to 50% of the photosynthetic biomass in the open oceans between 40°N and 40°S ##UREF##0##[1]##,##REF##17843668##[2]##. In these areas <italic>Prochlorococcus</italic> numerically dominates the phytoplankton with cell numbers reaching 10<sup>5</sup> cells per mL ##UREF##1##[3]##. Two major ecotypes can be differentiated within the <italic>Prochlorococcus</italic> group, which are relatively adapted to high or low light. They are genetically and physiologically distinct ##REF##9624000##[4]## and are distributed differently in the water column ##REF##12871235##[5]##–##REF##10347047##[7]##, with the high light adapted cells dominating the surface waters, and the low light adapted cells abundant in deep waters.</p> <p>The genomes of 12 <italic>Prochlorococcus</italic> strains, spanning the known microdiversity within the group, have been sequenced (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/genomes/MICROBES/microbial_taxtree.html">http://www.ncbi.nlm.nih.gov/genomes/MICROBES/microbial_taxtree.html</ext-link>). The cells posses the most streamlined genome of a free-living photoautotroph with genome sizes ranging from 1.6 Mbp to 2.7 Mbp ##REF##12917486##[8]##–##REF##12917642##[10]##. The number of modelled protein-coding genes in these genomes is 1,855–3,022 ##REF##18159947##[9]## and the core genome shared by all <italic>Prochlorococcus</italic> strains has been estimated at 1,273 genes ##REF##18159947##[9]##. Several hundred additional genes are specific for one or only a few strains, and they are frequently clustered in genomic islands ##REF##16556843##[11]##,##REF##18159947##[9]##.</p> <p>Genome reduction in this genus has particularly affected the number of regulatory genes. Many otherwise widely distributed two-component systems and DNA-binding proteins are not present in <italic>Prochlorococcus</italic>. This has been linked to the fitness gain conferred by a streamlined genome to organisms existing in a relatively stable environment ##REF##12917486##[8]##. Although the ocean environment may be relatively stable, it does fluctuate, making one wonder how <italic>Prochlorococcus</italic> cells respond to these changing conditions. Perhaps each protein regulator performs multiple regulatory functions in this cell. Alternatively, non-coding RNAs (ncRNAs) may play a major regulatory role compensating the lack of regulatory proteins.</p> <p>ncRNAs are functional RNA molecules, mostly without a protein-coding function, and their genes are normally located in intergenic regions. They frequently play a crucial role in bacterial regulatory networks particularly in response to environmental stress ##REF##9230301##[12]##,##REF##8670904##[13]## and are also known to control plasmid and viral replication ##REF##7826024##[14]##, bacterial virulence ##REF##17163975##[15]## and quorum sensing ##REF##15242645##[16]##. However the function of many ncRNAs remains unknown. <italic>Escherichia coli</italic> has over 70 ncRNAs most of which have been detected by computational prediction ##REF##11448770##[17]##–##REF##11445539##[20]## and “experimental RNomics” ##REF##15718303##[21]##–##REF##14602901##[23]##. These regulators were overlooked by traditional genome annotation due to their short length (50–400 nt in size), the lack of algorithms to search for sequences that are frequently more conserved in secondary structure rather than sequence, and the absence of a protein coding function.</p> <p>Another class of functional RNAs – chromosomally encoded antisense RNAs (asRNAs) also plays a role in the regulation of gene expression. There are no systematic approaches to screen for asRNAs, but RNomics approaches have inadvertently revealed the presence of asRNAs in <italic>Escherichia coli</italic>
##REF##15718303##[21]##–##REF##14602901##[23]##. These <italic>cis</italic>-encoded asRNAs are transcribed from the opposite strand of the same genomic locus as the target (m)RNA and feature 100% base complementarity. In contrast, most ncRNAs studied so far act in <italic>trans</italic> in a different genomic locus having only a short and imperfect base complementarity with the target transcripts (for a detailed review see ##REF##15952886##[24]##).</p> <p>Although a considerable number of <italic>Prochlorococcus</italic> strains have been fully sequenced, only a small number of ncRNAs have been identified in this group. In addition to the ubiquitous signal recognition particle RNA, RNAse P RNA and the tmRNA, encoded by <italic>ffs</italic>, <italic>rnpB</italic> and <italic>ssrA</italic>, seven ncRNAs have been identified in cyanobacteria, all of which were first described in <italic>Prochlorococcus</italic> MED4 and were denoted as Yfr1–Yfr7, for c<underline>Y</underline>anobacterial <underline>F</underline>unctional <underline>R</underline>NA ##REF##16168080##[25]##. Amongst them is Yfr7, which is homologous to 6S RNA ##REF##17640832##[26]## and known to have global regulatory functions in <italic>Escherichia coli</italic>. Another is Yfr1, which has homologues in other cyanobacteria ##REF##17941988##[27]## and in <italic>Synechococcus elongatus</italic> PCC6301, is required for growth under multiple stress conditions ##REF##17664182##[28]##. These two ncRNAs were classified as such after experimental verification of the expression of candidate ncRNAs initially identified from secondary structure conservation using a comparative genomics approach. Little is known regarding <italic>cis</italic>-acting asRNAs in cyanobacteria. Only 3 chromosomally <italic>cis</italic>-encoded asRNAs have been identified so far ##REF##2450053##[29]##–##REF##16324715##[31]##, none of which occur in <italic>Prochlorococcus</italic>.</p> <p>Despite the presence of ncRNAs in <italic>Prochlorococcus</italic> the gene encoding the Hfq RNA chaperone is absent from 10 of the 12 sequenced <italic>Prochlorococcus</italic> strains, including MED4. This is in contrast to other completely sequenced cyanobacteria that all contain an Hfq homologue. Hfq belongs to the eukaryotic and archaeal family of Sm and Sm-like (Lsm) proteins and is found in all domains of life. It facilitates the interaction of ncRNAs with their target mRNAs and is thus involved in many essential regulatory processes including ncRNA-mediated translational regulation ##REF##11804583##[32]##–##REF##11804582##[35]##. Its loss during evolution of the <italic>Prochlorococcus</italic> group may be taken as evidence for a general decay in RNA-dependent gene regulation or as an indication that novel mechanisms for RNA - RNA interactions may exist in this group.</p> <p>In the past few years, new experimental strategies such as ‘experimental RNomics’ and mining microarray expression data in intergenic regions have demonstrated that the number of ncRNAs in microbial genomes is much greater than previously thought (for reviews see ##REF##16436800##[36]##,##REF##16336117##[37]##). In the light of the small number of ncRNAs detected thus far in <italic>Prochlorococcus</italic> we were curious to see whether more ncRNAs are present in <italic>Prochlorococcus</italic> than were detected in the comparative genomics analysis used by Axman et al. ##REF##16168080##[25]##.</p> <p>Using an alternative approach based on microarray expression profiling, we investigated the presence of ncRNAs in <italic>Prochlorococcus</italic> MED4, which has the most compact genome of all sequenced <italic>Prochlorococcus</italic> strains, has few protein coding regulators, and is Hfq-deficient.</p>
<title>Materials and Methods</title> <title>Extraction of Microarray Expression Signals</title> <p>Three independent microarray experiments investigating global changes of gene expression under different light quantities and qualities (light experiment, ##REF##16980454##[38]##), under phage infection (phage experiment, ##REF##17805294##[39]##) and under phosphorus starvation (phosphorus experiment, ##REF##16895994##[40]##) were analyzed. The custom Affymetrix high-density array MED4-9313 that was used features 25-base oligomers identical to the target sequence that are spread over the complete genome comprising all gene coding regions as well as all intergenic regions on both forward and reverse strands with a coverage of every 45 bases in intergenic regions, a special feature that offers the detection of unknown transcripts. The Affymetrix array also contains probes for another <italic>Prochlorococcus</italic> genome (MIT9313) and two cyanophage genomes P-SSP7 and P-SSM4, whose average signal intensities were used to calculate threshold expression signals. For each set of experiments the threshold value used was re-evaluated to ensure high specificity of candidate probes. Specifically, we extracted probes with an expression signal of ≥200 in 18 of 21 arrays from the light experiment. Because of different experimental designs and thus resultant variations in overall expression signals, the threshold filter was adapted for the phage and phosphorus experiment extracting probes with expression signals ≥100 in 4 of 14 or 4 of 10 time points respectively (corresponding to the average of biological triplicates), respectively, with a 2-fold change in at least one time point between control and stress condition. The distribution of probe intensities was adjusted by quantile normalization across different arrays within the same experiment. This procedure minimized array-specific effects and allowed us to determine fold changes of single probes targeting ncRNAs, asRNAs and overlapping transcripts. Rather strict criteria for transcript identification were chosen to ensure a high true positive rate for transcript detection. To minimize the number of 5′ and 3′UTRs detected, probes within 100 nt of the adjacent gene in the same orientation were excluded. Remaining probes were grouped in transcriptional units and further characterized to categories: ORF, asRNA, ncRNA, 5′/3′UTR, pseudogenes and operon elements. The grouping and characterization is based on the localization in the genome and on BLAST searches against 11 <italic>Prochlorococcus</italic> genomes (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/genomes/MICROBES/microbial_taxtree.html">http://www.ncbi.nlm.nih.gov/genomes/MICROBES/microbial_taxtree.html</ext-link>) to identify conserved regions. Genes classified as ORFs encode for a peptide sequence with a start and stop codon without a frame shift and were present in at least two genomes. All PMED_xxxxx ORF notations (including new ORFs) follow that of Kettler et al. ##REF##18159947##[9]## and are available at <ext-link ext-link-type="uri" xlink:href="http://www.microbesonline.org">www.microbesonline.org</ext-link>. ncRNAs and asRNAs were defined as genes without peptide-coding potential localized in intergenic regions and opposite protein-coding genes, respectively. In two special cases ncRNAs with a regulatory RNA component as well as a peptide-encoded component were allowed. For detailed information about grouping see ##SUPPL##4##Table S2##.</p> <title>Culture Conditions</title> <p>
<italic>Prochlorococcus</italic> MED4 was grown at 21°C in AMP1 medium ##UREF##4##[61]## under 30 µmol quanta m<sup>−2</sup> s<sup>−1</sup> continuous white cool light. Culture conditions for microarray experiments are provided elsewhere ##REF##17805294##[39]##,##REF##16895994##[40]##,##REF##16980454##[38]##.</p> <title>5′- and 3′-RACE</title> <p>Total RNA was isolated as previously described ##REF##16980454##[38]## with the following modifications. Cells were harvested by centrifugation at 10,000×g for 10 min at 20°C. The pellet was resuspended in RNA resuspension buffer (10 mM sodium acetate [pH 5.2], 200 mM sucrose, 5 mM EDTA), snap frozen in liquid nitrogen and subsequently stored at −80°C. Total nucleic acids were DNase-treated with Turbo DNA-free (1 U/8 µg RNA, Ambion, USA) for 15 min at 37°C. RNA was precipitated with 1/10 volume 3 M sodium acetate (pH 5.2) and 3 volumes ethanol by centrifugation at 13,000×g for 30 min at 4°C and subsequently resuspended in water. Transcriptional start sites were determined by 5′-RACE following the method of Bensing et al. ##REF##8755555##[62]##. Briefly, RNA was treated with tobacco acid pyrophosphorylase (1 U/1 µg RNA; Epicentre, USA) for 1 h at 37°C followed by phenol/chloroform extraction and ethanol precipitation. A synthetic RNA oligonucleotide (0.5 µl oligonucleotide [10 µM]/ 4 µg RNA; <named-content content-type="gene">AUA UGC GCG AAU UCC UGU AGA ACG AAC ACU AGA AGA AA</named-content>, Invitrogen, Germany) was ligated to RNA using T4 RNA ligase (3 U/1 µg RNA; Fermentas, Germany) for 1 h at 37°C followed by phenol/chloroform extraction and ethanol precipitation. Three control reactions were performed: i) omitting tobacco acid pyrophosphorylase, ii) omitting tobacco acid pyrophosphorylase and RNA oligonucleotide and iii) dephosphorylating RNA prior to ligation with calf intestine alkaline phosphatase (0.1 U/1 µg RNA; Fermentas, Germany) at 37°C for 1 h, followed by phenol/chloroform extraction and ethanol precipitation. For reverse transcription 250 ng linked RNA per gene was incubated with 0.8 U of the Omniscript reverse transcriptase (Qiagen, Germany) in the provided reaction buffer containing 0.08 µM gene specific primer and 1 mM dNTPs. Incubation was carried out at 42°C for 2 h with a final inactivation step at 95°C for 5 min. All reactions were performed in the presence of 40 U Ribolock RNase Inhibitor (Fermentas, Germany). cDNA was amplified by PCR using a gene-specific primer (0.2 µM) and an RNA oligonucleotide-specific primer (0.2 µM) with following the cycling conditions: 93°C/3 min; 35 cycles of 93°C/30 s; 50°C/30 or 55°C/30 or 60°C/30 s, 72°C/45 s; 72°C/5 min in GoTaq reaction buffer containing 1 U GoTaq polymerase (Promega, Germany), 0.2 mM dNTPs and 3.5 mM MgCl<sub>2</sub>. A complete list with all primers used is provided in ##SUPPL##9##Table S7##. Amplified PCR fragments were gel-excised and purified on Nucleospin columns (Macherey &amp; Nagel, Germany) and then cloned into plasmid pGEMT (Promega, Germany). After transformation into <italic>E. coli</italic> XL1-Blue, plasmid inserts were amplified by colony PCR, purified on Nucleospin columns (Macherey &amp; Nagel, Germany) and sequenced using an ABI 3130XL automatic DNA sequencer (Applied Biosystems, USA). To determine the 3′ end of RNAs, 3′RACE was performed following the method described previously ##REF##11448770##[17]##. Briefly, RNA was treated as described above followed by a dephosphorylation with calf intestine alkaline phosphatase (0.2 U/1 µg RNA; Fermentas, Germany) at 37°C for 1 h and a subsequent phenol/chloroform extraction and ethanol precipitation. RNA 3′ ends were linked to a 3′ end blocked RNA oligonucleotide (0.5 µl oligonucleotide [10 µM]/4 µg RNA, <named-content content-type="gene">pAAG AUG AAU GCA ACA CUU CUG UAC GAC UAG AGC AC</named-content>, Metabion, Germany) using 0.8 U/1 µg RNA T4 RNA Ligase (Fermentas, Germany) followed by phenol/chloroform extraction and ethanol precipitation. Reverse transcription was performed as described above with the following modifications: 0.2 µM 3′ RNA oligonucleotide-specific primer and 2.5 mM dNTPs. Subsequent PCR, cloning and sequencing was performed as described above. Determined 5′ and 3′ ends are given in ##SUPPL##4##Tables S2## and ##SUPPL##10##S8##.</p> <title>Northern Analysis</title> <p>RNA samples (50 µg) were denatured for 5 min at 65°C in loading buffer (Fermentas, Germany), separated on 10% urea-polyacrylamide gels for 16 h at 100 V and transferred to Hybond-N nylon membranes (Amersham, Germany) by electroblotting for 1 h at 400 mA. The membranes were hybridized with specific [γ -<sup>32</sup>P]ATP end-labelled oligonucleotides or [α-<sup>32</sup>P]UTP-incorporated transcripts. Hybridization in 50% deionized formamide, 7% SDS, 250 mM NaCl and 120 mM Na(PO4) pH 7.2 was performed over night at 42°C or at 62°C with labelled oligonucleotide probes or labelled transcript probes, respectively. The membranes were washed in 2×SSC (3 M NaCl, 0.3 M sodium citrate, pH 7.0) ##REF##12944965##[55]##, 1% SDS for 10 minutes; 1×SSC, 0.5% SDS for 10 min; and briefly in 0.1×SSC, 0.1% SDS. All wash steps were performed 5°C below hybridization temperature. Signals were detected and analyzed on a Personal Molecular Imager FX system with Quantity One software (BIO-RAD, Germany).</p> <title>Oligonucleotide End Labelling</title> <p>Gene-specific oligonucleotides were labelled with [γ-<sup>32</sup>P]ATP by the exchange reaction of T4 polynucleotide kinase (Fermentas, Germany) using 0.5 U of enzyme, 1.25 µM oligonucleotide, 15 µCi [γ-<sup>32</sup>P]ATP in reaction buffer A for 30 min at 37°C followed by inactivation for 5 min at 95°C.</p> <title>In Vitro Transcription</title> <p>The MAXIscript Kit (Ambion, USA) was used for transcription of probes for use in Northern analyses containing 100 ng PCR-generated DNA template, 500 µM each of ATP, CTP, GTP, 20 µM UTP, 50 µCi [α-<sup>32</sup>P]UTP, 1 µl T7 enzyme mix in reaction buffer amended with SUPERase In RNase inhibitor (Ambion, USA). Transcription was carried out at 37°C for 10 min. Thereafter, the reactions were treated with 2 U of Turbo DNase-free (Ambion, USA) at 37°C for 15 min. The enzyme was heat inactivated for 10 min in the presence of 23 mM EDTA.</p>
<title>Results/Discussion</title> <title>Identification of Novel Transcripts</title> <p>The design of the <italic>Prochlorococcus</italic> custom Affymetrix microarray, which contains probes not only in gene-coding regions but also in intergenic regions (on both strands), and the availability of diverse data sets describing changes in gene expression in response to environmental stresses, allowed us to undertake a focused study of transcriptionally active intergenic regions. Three independent data sets were used: experiments investigating global changes in gene expression under different light quantities and qualities (from here on referred to as the “light experiment”, ##REF##16980454##[38]##), under phage infection (the “phage experiment”, ##REF##17805294##[39]##), and under phosphorus starvation (the “phosphorus experiment”, ##REF##16895994##[40]##), encompassing a total of 95 microarrays. To designate expression signals as novel transcripts, probes had to be above a threshold expression level, and be further than 100 nt from flanking genes (see <xref ref-type="sec" rid="s3">Materials and Methods</xref> for details). After identifying 553 probes that met these criteria from the light, phage and phosphorus experiments, we combined adjacent probes, yielding 276 unique transcriptional units (##FIG##0##Figure 1##). These transcripts were classified as 5′-UTRs, 3′–UTRs, operon elements, pseudogenes and “other” transcripts based on their genome location and experimental information. The “other” transcripts were then classified: as i) ORFs if they had a protein-coding reading frame with a start and stop codon without a frame shift in <italic>Prochlorococcus</italic> MED4 and in the genome of at least one other <italic>Prochlorococcus</italic> strain; ii) ncRNAs if they lacked an ORF, but had structural features typical of ncRNAs, such as compensatory mutations; and iii) asRNAs if they were located on the opposite strand of mRNAs. No assignment was made for 89 transcripts which either could not be verified in independent experiments (see below) or did not have homologues in other genomes. While some of these unclassified transcripts may represent unverified small ORFs or ncRNAs, they may also be the result of artificial expression signals. These may have occurred due to: i) cross-hybridization with duplicated regions that have only small sequence differences; or ii) artificial antisense signals caused by self-priming through hairpin loop extension of the first-strand cDNA, re-priming either from RNA fragments formed during degradation of the RNA templates, or from primers present in the reaction ##REF##17897965##[41]##.</p> <title>Identification of New ncRNAs</title> <p>Twelve novel ncRNAs were identified through microarray analyses. They were verified in independent experiments by rapid analysis of cDNA 5′ ends (5′ RACEs) and Northern hybridizations, which also served to map their first nucleotide and estimate their lengths (##TAB##0##Table 1##, ##FIG##1##Figure 2##). In addition, 5 of the 7 previously described ncRNAs (Yfr2 and Yfr4–Yfr7, ##REF##16168080##[25]##) were also detected. Yfr1 may not have been detected due to its extraordinary small size (54 nt) which may have resulted in its removal during the cDNA clean-up process. Regardless, Yfr1 would be excluded from our analysis because of its close proximity to an annotated protein-coding <italic>trxA</italic> gene (within 100 nt) that is transcribed in the same direction. Expression levels of Yfr3 were the lowest of the previously reported ncRNAs ##REF##16168080##[25]##, which likely explains why we did not detect this ncRNA in the microarray analysis. The internal consistency of these findings provide confidence in our approach, and suggest that even more ncRNAs may exist in <italic>Prochlorococcus</italic>, especially if they are very short, minimally expressed or close to protein-coding genes. Indeed two additional ncRNAs (Yfr12 and Yfr18) are described below that were identified by sequence homology, and that were not found through the microarray analysis because their expression signal was below the set threshold value.</p> <title>Distribution of ncRNAs within the MED4 Genome and Those of Other <italic>Prochlorococcus</italic> Strains</title> <p>Unlike Yfr1, which a sequence motif-based approach ##REF##17941988##[27]## revealed has homologues throughout the cyanobacterial lineage, none of our newly detected ncRNAs were universally present among the cyanobacteria. Indeed BlastN analyses yielded no evidence for their existence outside of the <italic>Prochlorococcus</italic> genus. With the exception of Yfr13, homologues of newly identified ncRNAs were only found in other high light-adapted <italic>Prochlorococcus</italic> strains (##TAB##0##Table 1##). Yfr13 has homologues in seven different <italic>Prochlorococcus</italic> strains, including the two low light-adapted isolates NATL1A and NATL2A, although the genome location is variable in the different strains (##FIG##2##Figure 3##).</p> <p>ncRNAs are non-randomly distributed within the MED4 genome. They are often associated with hypervariable genomic islands, thought to arise by horizontal gene transfer ##REF##16556843##[11]##. MED4 has 5 genomic islands that constitute only about one tenth of the total genome, whereas 9 of the 21 ncRNAs, described here or by Axmann et al. ##REF##16168080##[25]##, are in one of these islands (##TAB##0##Table 1##). The majority of island-associated ncRNAs are located in island 1 (Yfr8–Yfr11 and Yfr2). Three additional ncRNAs occur in island 2 (Yfr15, Yfr16 and Yfr3) and one is found in island 5 (Yfr20). Interestingly, the homologous ncRNAs in other <italic>Prochlorococcus</italic> strains are not always located in the corresponding island but occur somewhere else in the genome. The reverse is also true: some ncRNAs that are not island-associated in MED4 are located in an island region in other strains, indicative of recombination events. No ncRNAs were detected in island 4 even though this is the largest (74.5 kb long) of all islands present in MED4. Island 4 mainly encodes cell surface-relevant proteins such as glycosyltransferases or lipopolysaccharide-forming enzymes ##REF##16556843##[11]##, suggesting these functions are not controlled through ncRNAs. Yfr11 and Yfr16 are highly similar to each other. Based on their sequence identity (74%) and their highly similar secondary structures (##SUPPL##0##Figure S1##) both ncRNAs may regulate the same targets as has been shown for PrrF1 and PrrF2 in <italic>Pseudomonas aeruginosa</italic>
##REF##15210934##[42]## and for Qrr1, Qrr2, Qrr3, and Qrr4 in several <italic>Vibrio</italic> species ##REF##15242645##[16]##. Alternatively, they might act in a related context but with non-identical functions as has recently been described for GlmY and GlmZ ncRNAs of <italic>Escherichia coli</italic>
##REF##18351803##[43]##. Alternatively, Yfr11 and Yfr16 could be functionally equivalent but expressed in a different regulatory context, as is frequently the case for protein-coding genes that occur in multiple copies in a single genome.</p> <title>Differential Expression Suggests the Involvement of ncRNAs in Regulation during Stress Conditions</title> <p>Since the expression of many regulatory RNAs is coupled to the process they help regulate ##REF##9230301##[12]##,##REF##11395466##[44]##,##REF##15522088##[45]##, we explored the differential expression of the ncRNAs we identified as a function of different environmental stresses. The expression levels of several ncRNAs were influenced by light and phage induced stress, but not by phosphorus stress. Two ncRNAs – Yfr19 and Yfr11 were more than twofold downregulated after transfer from darkness to high white light, normal white light or blue light, but were upregulated when DCMU (an inhibitor of the photosynthetic electron transport chain) was added to cells grown in normal white light conditions (##FIG##3##Figure 4##, ##SUPPL##3##Table S1##). Expression of Yfr16, the homolog of Yfr11, followed the same trends, but was less pronounced than for Yfr11. Both the reduced transcript levels during light exposure and the increased amount upon DCMU treatment indicate a link between the redox status of the photosynthetic electron transport chain and these three ncRNAs. High light induced differential expression in the largest number of ncRNAs, and of the highest magnitude, as has been observed for the response of protein-coding genes ##REF##16980454##[38]##. However, only a single ncRNA, Yfr20, was upregulated when cells were transferred from darkness to high light (##FIG##3##Figure 4##, ##SUPPL##3##Table S1##), whereas all other ncRNAs responsive to light stress decreased in their transcript levels. Yfr20 accumulates in high absolute amounts (##FIG##1##Figure 2##). According to 5′RACE, the major accumulating transcript of 89 nt results from a specific initiation of transcription at position 1336435 (accession number BX548174.1). In addition 5′RACE results show that Yfr20 is transcribed together with the upstream located ORF PMED4_15791 as a dicistronic element. PMED4_15791 showed constitutive expression. Thus, the light-dependent expression of Yfr20 is under control of its own promoter. The dicistronic gene arrangement with an upstream located ORF is split in <italic>Prochlorococcus</italic> strains MIT9515 and MIT9312 and contains an additional <italic>hli</italic> gene in between the ncRNA and the ORF homologous to PMED4_15791 (##FIG##4##Figure 5##) providing further evidence for a possible light-regulatory function of Yfr20. Intriguingly, Yfr20 is the only ncRNA encoded in genomic island 5. This island has been characterized as a “phosphorus” island in MED4 since nine genes (nearly all of unknown function) responded when MED4 cells were starved for phosphorus ##REF##16895994##[40]##. However, high light stress caused an additional 15 genes to respond in genomic island 5 ##REF##16556843##[11]##, among them <italic>hli11</italic> and <italic>hli12</italic>, which are located at a distance of less than 2 kb from <italic>yfr20</italic> (##FIG##4##Figure 5##). Although <italic>hli</italic> (high light inducible) proteins can be factors in other stress responses as well, their mode of regulation here indeed suggests that this island plays a role not only in the adaptation to phosphorus starvation but also to stress caused by high light.</p> <p>Two distinct stress responses of the cell in response to phage infection have been identified. Lindell and co-workers ##REF##17805294##[39]## observed an overall reduction in expression of host genes as the major response to phage infection. However, 41 protein-coding host genes were upregulated in the initial or the mid-to-late phases of phage infection. It is hypothesized that genes belonging to group 1 (the first wave of upregulation) constitute a direct defence to phage infection whereas group 2 genes (the second wave) may be induced by the phage ##REF##17805294##[39]##. Two ncRNAs – Yfr9 and Yfr14 - were upregulated in the initial phase of infection (from 1–3 hours after infection, corresponding to group 1 upregulated protein-coding genes ##REF##17805294##[39]##, ##FIG##3##Figure 4##, ##SUPPL##3##Table S1##). Interestingly, both ncRNAs have an antisense-located ncRNA – Yfr8 and Yfr6 – that are constitutively expressed during that time (##FIG##3##Figure 4##, ##SUPPL##3##Table S1##). The two pairs of overlapping ncRNAs are characterized in more detail below (see section on overlapping ncRNAs). An additional ncRNA – Yfr15 - was upregulated during the mid to late phases of infection (from 3 to 8 hours corresponding to group 2 upregulated protein-coding genes ##REF##17805294##[39]##, ##FIG##3##Figure 4##, ##SUPPL##3##Table S1##). Yfr15 is located in genomic island 2 in the vicinity of PMED4_07441 (PMM0686), the most highly upregulated host mRNA during phage infection, although the two genes are located on opposite strands. Also PMED4_07401 (PMM0684) and PMED4_07421 (PMM0685), two further genes that belong to group 2 phage-induced host genes, are located nearby in genomic island 2. This region and Yfr15 may therefore be of prime importance for phage-host interactions.</p> <p>We did not detect a single ncRNA that was significantly differentially expressed under phosphorus limitation. This was very surprising, in light of the 34 protein encoding-genes that are differentially expressed under P-stress in MED4 ##REF##16895994##[40]##, and because in <italic>Escherichia coli</italic> the existence of such ncRNAs was hypothesized based on the observation of Hfq-dependent regulation of <italic>rpoS</italic> in response to this stress ##REF##14526009##[46]##.</p> <title>An Ultraconserved Sequence Motif is Present in Three ncRNAs</title> <p>The ncRNA Yfr10 contains the conserved unadecanucleotide motif <named-content content-type="gene">5′-ACUCCUCACAC-3′</named-content> (##FIG##5##Figure 6##). This motif occurs 3 times in the MED4 genome sequence, which is more frequent than expected by chance: One would expect approximately 0.5 instances of a specific 11 nt motif in a 2 MB genome at equal base distribution. The second occurrence has already been described as belonging to another ncRNA in MED4, Yfr1 ##REF##17941988##[27]##, which is also found throughout the cyanobacterial radiation.##REF##17941988##[27]##. Using Northern analysis and 5′RACE, we showed that also the third copy of this motif is expressed, revealing another ncRNA - Yfr18 (##FIG##1##Figure 2##, ##TAB##0##Table 1##). This one was not detected from our microarray analyses nor the comparative genomics approach ##REF##16168080##[25]##.</p> <p>If two base transitions are allowed there is even a fourth member of this ‘unadecanucleotide-containing’ class of ncRNAs in MED4. This ncRNA - Yfr12 - was identified by sequence similarity to Yfr10 and was verified as an ncRNA as described above (##FIG##1##Figure 2##, ##TAB##0##Table 1##). Yfr12 appears to be a mutated variant of the other three, since the processed 5′ end of the major accumulating RNA species was mapped to the middle of the unadecanucleotide, and two mutations change the sequence at the 3′ end of the motif from <named-content content-type="gene">CACACAC</named-content> to <named-content content-type="gene">CAUAUAC</named-content> (##FIG##5##Figure 6##). Furthermore, the motif can be extended in all 4 ncRNAs by another AC dinucleotide (##FIG##5##Figure 6##) in comparison to the published cyanobacterial consensus ##REF##17941988##[27]##, probably a peculiarity of these RNAs in <italic>Prochlorococcus</italic>.</p> <p>The functions of Yfr10, Yfr12 and Yfr18 in MED4 remain unknown at present. However, a hint about their potential function may be found from their genome context and the fact that the vast majority of functional interactions between ncRNAs and their targets is exerted through base pairing. The genes for Yfr10, Yfr12 and Yfr18 are each directly adjacent to those for ncRNAs Yfr2, Yfr4 and Yfr5 respectively (##TAB##0##Table 1##), the 5′ ends of which may basepair to the 13-nucleotide consensus of Yfr1, Yfr10, Yfr12 and Yfr18 if a single bulging C and one mismatch is allowed (i.e. <named-content content-type="gene">5′-aCUCCUcACACAC-3′</named-content> pairs with <named-content content-type="gene">5′-GUGUGUAGGAG-3′</named-content>). Moreover, one may note that also the two C – to - U transitions in Yfr12 are compatible with this suggested base pairing, and that secondary structure predictions suggest that the conserved motifs in Yfr1, Yfr10 and Yfr18 are exposed as single stranded elements in an otherwise folded region (##FIG##5##Figure 6##). The same is true for the complementary motif in Yfr2–Yfr5, making physical interactions very likely. Evidence for two regulatory RNAs acting upon each other has recently been reported for the first time for GlmY and GlmZ of <italic>Escherichia coli</italic>, and cascades of hierarchically acting regulatory RNAs have been hypothesized for other bacteria as well ##REF##18351803##[43]##. The ncRNAs described in this section are candidates for such interactions in <italic>Prochlorococcus</italic>.</p> <title>Two Overlapping Pairs of ncRNAs</title> <p>The difference between <italic>trans</italic>- and <italic>cis</italic>-encoded ncRNAs and asRNAs is frequently considered fuzzy since both act through base complementarity. However, depending on the length of the overlap, interactions between transcripts from the forward and the reverse DNA strand can be very strong due to the extended perfect sequence complementarity. asRNAs may act as the antidote in toxin-antitoxin systems ##REF##12123448##[47]##,##REF##16385044##[48]## or in gene regulation ##REF##16636284##[30]##,##REF##16324715##[31]## as has been reported for some bacteria. We detected two regions with probable sense/antisense pairing between ncRNAs. One of these regions contains Yfr6 with Yfr14 on the opposite strand (##FIG##6##Figure 7##, ##TAB##0##Table 1##). The second region is located in genomic island 1 containing Yfr8 and Yfr9, each of which are 290 nt in size (##FIG##1##Figures 2## and ##FIG##6##7##, ##TAB##0##Table 1##). One of the ncRNAs in both pairs is upregulated during phage infection (##FIG##3##Figure 4##) and both ncRNA pairs contain a potential peptide-coding open reading frame within the sequence of one of the RNAs (##FIG##6##Figure 7##). The peptide sequence associated with Yfr6 is 33 amino acids long and is highly conserved and widely distributed among high light- and low light <italic>Prochlorococcus</italic> strains. The potential 44 amino acids peptide-coding frame within Yfr9 was found in three other <italic>Prochlorococcus</italic> isolates (MIT9515, MIT9301, AS9601), but has been lost in <italic>Prochlorococcus</italic> strains MIT9312 and MIT9215 due to a frame shift. Homologues of Yfr6 and Yfr9, respectively, have high sequence conservation over their complete ncRNA genes – including the upstream and downstream regions of the potential peptides.</p> <p>Yfr6 and Yfr9 resemble RNAIII from <italic>Staphylococcus aureus</italic> being both relatively long and consisting of a small peptide-coding unit as well as a regulatory RNA. RNAIII is a 510 nt long riboregulator from which the 26 amino acid δ-hemolysin peptide is also translated ##REF##9620969##[49]##. Another bifunctional ncRNA (SgrS) has been described in <italic>Escherichia coli</italic> that contains a conserved ORF (SgrT) in the 5′ region of SgrS, both of which promote recovery from glucose stress in mechanistically distinct fashions ##REF##18042713##[50]##. Moreover, secondary structure predictions of Yfr6 ##REF##16168080##[25]## and of Yfr14, Yfr8, and Yfr9 (##SUPPL##1##Figures S2## and ##SUPPL##2##S3##) support the potential role of these transcripts as functional RNAs as they contain many G - C base pairings and compensatory mutations, which conserve the structure rather than the sequence – a feature of many ncRNAs. On the other hand, highly structured transcript regions are also found in certain mRNAs where they serve as platforms for sophisticated ncRNA-mediated translational control. In the case of the <italic>Escherichia coli tisB</italic> mRNA, for example, this transcript encodes a peptide as short as 29 amino acids towards its 3′ end, yet there is no evidence that <italic>tisB</italic> would act as a riboregulator ##REF##11448770##[17]##,##REF##17499044##[51]##.</p> <p>The fact that Yfr6 and Yfr9 overlap with other ncRNAs (Yfr14 and Yfr8 respectively), suggests that this could be a toxin-antitoxin system – i.e. pairs of genes that code for a stable toxin and an unstable antitoxin. These are well-characterized in other bacteria, where the toxin is usually a toxic peptide that is neutralized or whose synthesis is prevented by the action of the product of the second gene, the antitoxin, which is either protein or RNA. Toxin-antitoxin systems such as the <italic>hok/sok</italic> system of <italic>Escherichia coli</italic> can serve as a natural genetic selection system to ensure presence of a plasmid ##UREF##3##[52]##. Alternatively, chromosomally encoded toxin-antitoxin systems can be beneficial to cell survival under unfavorable growth conditions, sometimes in very sophisticated ways, for instance by transiently curtailing the consumption of nutrients during starvation or by temporarily inhibiting growth and thereby evading the killing effects of certain antibiotics ##REF##15864262##[53]##. Systematic searches for toxin-antitoxin systems have revealed a high abundance in free-living prokaryotes ##REF##15718296##[54]## but none of the seven known toxin-antitoxin families could be identified in <italic>Prochlorococcus</italic>. There is a growing number of examples of chromosomal toxin-antitoxin systems that use a <italic>cis</italic>-encoded asRNA as an antitoxin. Our data suggest that Yfr6/Yfr14 and Yfr8/Yfr9 may be potential candidates for toxin-antitoxin systems in <italic>Prochlorococcus</italic> MED4.</p> <title>A High Number of <italic>cis</italic>-Encoded Antisense RNAs</title> <p>Little attention has been given to chromosomally <italic>cis</italic>-encoded asRNAs until recently, and only a few have been described for cyanobacteria ##REF##2450053##[29]##–##REF##16324715##[31]##. Surprisingly, we detected 24 asRNAs in our analyses, which vary between 100 to 600 nt in size (##SUPPL##4##Table S2##). Some are differentially expressed under different light conditions and under phage infection (##SUPPL##5##Table S3##).</p> <p>High light treatment caused one asRNA to be upregulated and one to be downregulated (##SUPPL##5##Table S3##). asRNA asMED4_15721 was upregulated twofold when cells were transferred from darkness to high light (##SUPPL##5##Table S3##). This behavior is similar to that of Yfr20 (see above), and like Yfr20, this asRNA is located in genomic island 5, lending additional support for a function of this island in light stress adaptation. Notably, five of the 24 asRNAs are complementary to mRNAs that code for photosystem I subunits (<italic>psaB</italic> and <italic>psaC</italic>) or for photosystem II subunits (<italic>psbB</italic> and <italic>psbO</italic>, <italic>psbX</italic>), respectively. The concentrations of asRNAs of photosystem II genes did not change when cultures were shifted from darkness to different light quantities and qualities, whereas transcript levels of their target mRNA decreased slightly (##SUPPL##5##Table S3##). In contrast, levels of the photosystem I asRNAs asMED4_17331 (antisense of <italic>psaB</italic>) and asMED4_18171 (antisense of <italic>psaC</italic>) decreased when cells were shifted from darkness to light, following the same trend as their mRNA counterparts (##SUPPL##5##Table S3##). Surprisingly, however, the latter asRNAs decreased in amount when transferred from darkness to medium white light whereas their target mRNAs did not (##SUPPL##5##Table S3##), which might indicate a light-dosage specific regulation of these asRNAs.</p> <p>We also found asRNAs that are differentially expressed in cells infected by phage. asMED4_04601 is upregulated during the initial stages of phage infection, whereas its target mRNA (PMED4_04601) is constitutively expressed throughout the infection process. Interestingly, PMED4_04601 shows 67% amino acid identity to the central region of the potential Yfr6 peptide. In the case of PMED4_07401 (PMM0684) both the respective asRNA (asMED4_07401) and its target mRNA are upregulated from mid-to-late phase of phage infection.</p> <p>The number of 24 asRNAs detected in our analyses appears high, especially as the microarrays used for this study did not contain probes for antisense regions of protein coding genes (see <xref ref-type="sec" rid="s3">methods</xref>). Therefore we could only detect those asRNAs found in intergenic regions whose corresponding ORF was not originally annotated and asRNAs located in 5′ and 3′ UTRs. However, our mapping results revealed that asRNAs located opposite of 5′/3′ UTRs frequently overlap major parts of the adjacent coding sequences. While it is not possible to infer the functions of these asRNAs in <italic>Prochlorococus</italic>, in the cyanobacterium <italic>Synechocystis</italic> PCC6803 the asRNA IsrR occurs in higher quantities than the <italic>cis</italic>-encoded mRNA <italic>isiA</italic> under normal growth conditions, leading to degradation of RNA duplexes by RNase III ##REF##16636284##[30]##. Under stress conditions the expression of the mRNA is increased leading to free mRNA molecules that can be translated. This mode of action is highly unlikely in the case of the PMED4_07401 mRNA/asRNA pair, because they are co-upregulated during phage infection, pointing towards a protective rather than a degradative role. Interestingly, RNAse E is also among the genes upregulated during phage infection. It is hypothesized that ribonuclease activity could be utilized by the phage to degrade host RNA to generate nucleotides for phage replication ##REF##17805294##[39]##.</p> <title>Other Overlapping Transcripts</title> <p>In addition to non-protein coding asRNAs, we found several pairs of transcripts that are transcribed from complementary strands and that potentially code for proteins. These complementary transcripts overlap entirely with each other or with a major part of their 5′ or 3′ UTR. We confirmed these overlapping regions experimentally and found that they span between 74 nt to 333 nt at least (##SUPPL##4##Table S2##). Eight out of twelve of these overlapping regions are found in the same position in other <italic>Prochlorococcus</italic> genomes, whereas the other 4 are found in different regions of the genomes. At this point it is not clear whether the overlaps between these protein-encoding transcripts would interfere with their transcription, transcript accumulation or translation. We found evidence for both scenarios. Whereas PMED4_14671 (located in the opposite 3′UTR region of PMED4_14661) is upregulated 14 fold when light intensity is increased, PMED4_14661 (PMM1300) remains at basal transcript levels. Contrary to the above, PMED4_11211 (PMM0997) and PMED4_11201 (located in the opposite 3′UTR region of PMED4_11211) are inversely regulated under different light conditions and DCMU treatment (##SUPPL##6##Table S4##) indicative of either interference during transcription, or coupled degradation, as observed for the asRNA and mRNA IsrR/<italic>isiA</italic> in the cyanobacterium <italic>Synechocystis</italic> PCC6803 ##REF##16636284##[30]##.</p> <p>In cases where only 3′-UTRs overlap this might not be of relevance because the transcriptional machinery should not be constricted. However, most of the transcripts we found overlap the 5′-UTRs and/or complete protein-coding regions and therefore are highly likely to be of regulatory relevance.</p> <title>Identification of New ORFs</title> <p>It is difficult to identify short genes in bacterial genomes using annotation algorithms, because the number of possible reading frames increases the shorter the search window becomes. Therefore, many of the widely used annotation programs (e.g. GLIMMER, GeneMark and CRITICA) constrain the minimum length of an ORF and thus a considerable number of small ORFs remain unannotated. In a recent study of the <italic>Prochlorococcus</italic> pan genome, for example, hypothetical ORFs shorter than 50 amino acids were excluded unless they were found in more than one genome ##REF##18159947##[9]##.</p> <p>Our microarray analyses lead to the observation of 113 new ORFs (not including asRNAs with potential protein coding sequences; ##SUPPL##4##Table S2##) ranging between 33 to 130 amino acids in size that were not annotated in the first published <italic>Prochlorococcus</italic> MED4 genome version (accession number BX548174). The new annotation of Kettler et al. ##REF##18159947##[9]## (accession number BX548174.1, for new ORF IDs refer to: <ext-link ext-link-type="uri" xlink:href="http://www.microbesonline.org">www.microbesonline.org</ext-link>) also found 89 of the new ORFs. BlastP searches against the non-redundant NCBI database revealed that 10 of the 24 remaining ORFs have an annotated counterpart in at least one other genome whereas 14 ORFs represent short proteins that have not been detected thus far (##SUPPL##7##Table S5##). Using TblastN, all of the additional 14 novel ORFs were found in other genomes, even though they had gone undetected by computational annotation tools (##SUPPL##8##Table S6##).</p> <title>Phage-Related Open Reading Frames</title> <p>Amongst the newly discovered protein-coding genes are three ORFs that have homologues in cyanophage genomes. PMED4_16122, which is located in genomic island 5 in MED4, is homologous to PSSM4_095 in the MED4-infecting cyanophage PSSM4 ##REF##12944965##[55]##, suggesting gene transfer between a PSSM4-like phage and this island, in particular since no other homologues were found in the nr database. PMED4_15491 is in vicinity of genomic island 5 in the host genome and has one homolog in each of two cyanophage genomes – P-SSM4 and P-SSM2 (PSSM4_181 and PSSM2_278). According to TblastN results this ORF is present in numerous <italic>Prochlorococcus</italic> genomes but not in other cyanobacteria. The third newly discovered ORF with a homolog in a cyanophage is PMED4_10681, which is found in the genome of cyanophage P-SSM2. Unlike PMED4_16122, however, this ORF is not in a genomic island in the host genome, and furthermore, is widely distributed over the cyanobacterial radiation with homologues in all <italic>Prochlorococcus</italic> strains, <italic>Synechococcus elongatus</italic> strains PCC 6301 and PCC 7942, <italic>Fremyella diplosiphon</italic>, <italic>Nostoc</italic> sp. PCC 7120, <italic>Anabaena variabilis</italic> and <italic>Synechocystis</italic> sp. PCC 6803 (##SUPPL##8##Table S6##). The broad distribution of PMED4_10681 suggests that it plays an important function in cyanobacteria, and emphasizes the importance of better annotation of small ORFs.</p> <title>Conclusions</title> <p>Here we have described 14 novel ncRNAs, which increases the total number of ncRNAs in this organism to 24 (including Yfr1-7, ffs, tmRNA, RNase P RNA). One sixth of the 24 ncRNAs (Yfr1, Yfr3, Yfr12 and Yfr18) were undetectable from microarray analyses under the conditions tested. Therefore it is likely that even more ncRNAs are present in <italic>Prochlorococcus</italic> MED4. The proportion of ncRNAs in the <italic>Prochlorococcus</italic> MED4 genome is comparable with those found in enterobacteria like <italic>Escherichia coli</italic>, i.e. 1–2% of the genes encode ncRNAs. In comparison, the 6 identified protein regulators in <italic>Prochlorococcus</italic>
##REF##12917642##[10]## is a small number relative to the 32 two-component response regulators present in <italic>Escherichia coli</italic>
##REF##9205844##[56]##. This suggests that regulation of gene expression through ncRNAs plays an important role in <italic>Prochlorococcus</italic>' response to environmental cues. The relatively high number of ncRNAs is intriguing as it may represent a mode of adaptation to the extremely low nutrient conditions of the open oceans. Regulation by ncRNAs may require fewer resources than would be required for the synthesis of protein regulators. Furthermore, in the course of genome reduction there might have been a positive selection pressure for keeping small regulators, e.g. ncRNAs rather than large protein regulators. How ncRNAs function in <italic>Prochlorococcus</italic> is at present unclear. The absence of Hfq in MED4 suggests that the ncRNAs found in this strain represent a core-set of ncRNAs that function without the support of a chaperone, or with a novel chaperone yet to be identified.</p> <p>The genomic islands of <italic>Prochlorococcus</italic> are disproportionately connected to ecological functions in this group of cyanobacteria ##REF##17693088##[57]##,##REF##16556843##[11]##. Here we have shown that approximately half of the <italic>Prochlorococcus</italic> ncRNAs are located in genomic islands suggesting that the function of these molecules is relevant for determining the relative fitness of ecotypes within <italic>Prochlorococcus</italic>. This is analogous to the accumulation of genes coding for ncRNA in pathogenicity islands in <italic>Staphylococcus aureus</italic>
##REF##16183745##[58]## and <italic>Salmonella typhimurium</italic>
##REF##17971080##[59]## as well as in genomic islands of <italic>Sinorhizobium meliloti</italic>
##REF##18093320##[60]##, and suggests that this phenomenon could be wide-spread for finely tuned specialization within microbial groups.</p>
<title>Results/Discussion</title> <title>Identification of Novel Transcripts</title> <p>The design of the <italic>Prochlorococcus</italic> custom Affymetrix microarray, which contains probes not only in gene-coding regions but also in intergenic regions (on both strands), and the availability of diverse data sets describing changes in gene expression in response to environmental stresses, allowed us to undertake a focused study of transcriptionally active intergenic regions. Three independent data sets were used: experiments investigating global changes in gene expression under different light quantities and qualities (from here on referred to as the “light experiment”, ##REF##16980454##[38]##), under phage infection (the “phage experiment”, ##REF##17805294##[39]##), and under phosphorus starvation (the “phosphorus experiment”, ##REF##16895994##[40]##), encompassing a total of 95 microarrays. To designate expression signals as novel transcripts, probes had to be above a threshold expression level, and be further than 100 nt from flanking genes (see <xref ref-type="sec" rid="s3">Materials and Methods</xref> for details). After identifying 553 probes that met these criteria from the light, phage and phosphorus experiments, we combined adjacent probes, yielding 276 unique transcriptional units (##FIG##0##Figure 1##). These transcripts were classified as 5′-UTRs, 3′–UTRs, operon elements, pseudogenes and “other” transcripts based on their genome location and experimental information. The “other” transcripts were then classified: as i) ORFs if they had a protein-coding reading frame with a start and stop codon without a frame shift in <italic>Prochlorococcus</italic> MED4 and in the genome of at least one other <italic>Prochlorococcus</italic> strain; ii) ncRNAs if they lacked an ORF, but had structural features typical of ncRNAs, such as compensatory mutations; and iii) asRNAs if they were located on the opposite strand of mRNAs. No assignment was made for 89 transcripts which either could not be verified in independent experiments (see below) or did not have homologues in other genomes. While some of these unclassified transcripts may represent unverified small ORFs or ncRNAs, they may also be the result of artificial expression signals. These may have occurred due to: i) cross-hybridization with duplicated regions that have only small sequence differences; or ii) artificial antisense signals caused by self-priming through hairpin loop extension of the first-strand cDNA, re-priming either from RNA fragments formed during degradation of the RNA templates, or from primers present in the reaction ##REF##17897965##[41]##.</p> <title>Identification of New ncRNAs</title> <p>Twelve novel ncRNAs were identified through microarray analyses. They were verified in independent experiments by rapid analysis of cDNA 5′ ends (5′ RACEs) and Northern hybridizations, which also served to map their first nucleotide and estimate their lengths (##TAB##0##Table 1##, ##FIG##1##Figure 2##). In addition, 5 of the 7 previously described ncRNAs (Yfr2 and Yfr4–Yfr7, ##REF##16168080##[25]##) were also detected. Yfr1 may not have been detected due to its extraordinary small size (54 nt) which may have resulted in its removal during the cDNA clean-up process. Regardless, Yfr1 would be excluded from our analysis because of its close proximity to an annotated protein-coding <italic>trxA</italic> gene (within 100 nt) that is transcribed in the same direction. Expression levels of Yfr3 were the lowest of the previously reported ncRNAs ##REF##16168080##[25]##, which likely explains why we did not detect this ncRNA in the microarray analysis. The internal consistency of these findings provide confidence in our approach, and suggest that even more ncRNAs may exist in <italic>Prochlorococcus</italic>, especially if they are very short, minimally expressed or close to protein-coding genes. Indeed two additional ncRNAs (Yfr12 and Yfr18) are described below that were identified by sequence homology, and that were not found through the microarray analysis because their expression signal was below the set threshold value.</p> <title>Distribution of ncRNAs within the MED4 Genome and Those of Other <italic>Prochlorococcus</italic> Strains</title> <p>Unlike Yfr1, which a sequence motif-based approach ##REF##17941988##[27]## revealed has homologues throughout the cyanobacterial lineage, none of our newly detected ncRNAs were universally present among the cyanobacteria. Indeed BlastN analyses yielded no evidence for their existence outside of the <italic>Prochlorococcus</italic> genus. With the exception of Yfr13, homologues of newly identified ncRNAs were only found in other high light-adapted <italic>Prochlorococcus</italic> strains (##TAB##0##Table 1##). Yfr13 has homologues in seven different <italic>Prochlorococcus</italic> strains, including the two low light-adapted isolates NATL1A and NATL2A, although the genome location is variable in the different strains (##FIG##2##Figure 3##).</p> <p>ncRNAs are non-randomly distributed within the MED4 genome. They are often associated with hypervariable genomic islands, thought to arise by horizontal gene transfer ##REF##16556843##[11]##. MED4 has 5 genomic islands that constitute only about one tenth of the total genome, whereas 9 of the 21 ncRNAs, described here or by Axmann et al. ##REF##16168080##[25]##, are in one of these islands (##TAB##0##Table 1##). The majority of island-associated ncRNAs are located in island 1 (Yfr8–Yfr11 and Yfr2). Three additional ncRNAs occur in island 2 (Yfr15, Yfr16 and Yfr3) and one is found in island 5 (Yfr20). Interestingly, the homologous ncRNAs in other <italic>Prochlorococcus</italic> strains are not always located in the corresponding island but occur somewhere else in the genome. The reverse is also true: some ncRNAs that are not island-associated in MED4 are located in an island region in other strains, indicative of recombination events. No ncRNAs were detected in island 4 even though this is the largest (74.5 kb long) of all islands present in MED4. Island 4 mainly encodes cell surface-relevant proteins such as glycosyltransferases or lipopolysaccharide-forming enzymes ##REF##16556843##[11]##, suggesting these functions are not controlled through ncRNAs. Yfr11 and Yfr16 are highly similar to each other. Based on their sequence identity (74%) and their highly similar secondary structures (##SUPPL##0##Figure S1##) both ncRNAs may regulate the same targets as has been shown for PrrF1 and PrrF2 in <italic>Pseudomonas aeruginosa</italic>
##REF##15210934##[42]## and for Qrr1, Qrr2, Qrr3, and Qrr4 in several <italic>Vibrio</italic> species ##REF##15242645##[16]##. Alternatively, they might act in a related context but with non-identical functions as has recently been described for GlmY and GlmZ ncRNAs of <italic>Escherichia coli</italic>
##REF##18351803##[43]##. Alternatively, Yfr11 and Yfr16 could be functionally equivalent but expressed in a different regulatory context, as is frequently the case for protein-coding genes that occur in multiple copies in a single genome.</p> <title>Differential Expression Suggests the Involvement of ncRNAs in Regulation during Stress Conditions</title> <p>Since the expression of many regulatory RNAs is coupled to the process they help regulate ##REF##9230301##[12]##,##REF##11395466##[44]##,##REF##15522088##[45]##, we explored the differential expression of the ncRNAs we identified as a function of different environmental stresses. The expression levels of several ncRNAs were influenced by light and phage induced stress, but not by phosphorus stress. Two ncRNAs – Yfr19 and Yfr11 were more than twofold downregulated after transfer from darkness to high white light, normal white light or blue light, but were upregulated when DCMU (an inhibitor of the photosynthetic electron transport chain) was added to cells grown in normal white light conditions (##FIG##3##Figure 4##, ##SUPPL##3##Table S1##). Expression of Yfr16, the homolog of Yfr11, followed the same trends, but was less pronounced than for Yfr11. Both the reduced transcript levels during light exposure and the increased amount upon DCMU treatment indicate a link between the redox status of the photosynthetic electron transport chain and these three ncRNAs. High light induced differential expression in the largest number of ncRNAs, and of the highest magnitude, as has been observed for the response of protein-coding genes ##REF##16980454##[38]##. However, only a single ncRNA, Yfr20, was upregulated when cells were transferred from darkness to high light (##FIG##3##Figure 4##, ##SUPPL##3##Table S1##), whereas all other ncRNAs responsive to light stress decreased in their transcript levels. Yfr20 accumulates in high absolute amounts (##FIG##1##Figure 2##). According to 5′RACE, the major accumulating transcript of 89 nt results from a specific initiation of transcription at position 1336435 (accession number BX548174.1). In addition 5′RACE results show that Yfr20 is transcribed together with the upstream located ORF PMED4_15791 as a dicistronic element. PMED4_15791 showed constitutive expression. Thus, the light-dependent expression of Yfr20 is under control of its own promoter. The dicistronic gene arrangement with an upstream located ORF is split in <italic>Prochlorococcus</italic> strains MIT9515 and MIT9312 and contains an additional <italic>hli</italic> gene in between the ncRNA and the ORF homologous to PMED4_15791 (##FIG##4##Figure 5##) providing further evidence for a possible light-regulatory function of Yfr20. Intriguingly, Yfr20 is the only ncRNA encoded in genomic island 5. This island has been characterized as a “phosphorus” island in MED4 since nine genes (nearly all of unknown function) responded when MED4 cells were starved for phosphorus ##REF##16895994##[40]##. However, high light stress caused an additional 15 genes to respond in genomic island 5 ##REF##16556843##[11]##, among them <italic>hli11</italic> and <italic>hli12</italic>, which are located at a distance of less than 2 kb from <italic>yfr20</italic> (##FIG##4##Figure 5##). Although <italic>hli</italic> (high light inducible) proteins can be factors in other stress responses as well, their mode of regulation here indeed suggests that this island plays a role not only in the adaptation to phosphorus starvation but also to stress caused by high light.</p> <p>Two distinct stress responses of the cell in response to phage infection have been identified. Lindell and co-workers ##REF##17805294##[39]## observed an overall reduction in expression of host genes as the major response to phage infection. However, 41 protein-coding host genes were upregulated in the initial or the mid-to-late phases of phage infection. It is hypothesized that genes belonging to group 1 (the first wave of upregulation) constitute a direct defence to phage infection whereas group 2 genes (the second wave) may be induced by the phage ##REF##17805294##[39]##. Two ncRNAs – Yfr9 and Yfr14 - were upregulated in the initial phase of infection (from 1–3 hours after infection, corresponding to group 1 upregulated protein-coding genes ##REF##17805294##[39]##, ##FIG##3##Figure 4##, ##SUPPL##3##Table S1##). Interestingly, both ncRNAs have an antisense-located ncRNA – Yfr8 and Yfr6 – that are constitutively expressed during that time (##FIG##3##Figure 4##, ##SUPPL##3##Table S1##). The two pairs of overlapping ncRNAs are characterized in more detail below (see section on overlapping ncRNAs). An additional ncRNA – Yfr15 - was upregulated during the mid to late phases of infection (from 3 to 8 hours corresponding to group 2 upregulated protein-coding genes ##REF##17805294##[39]##, ##FIG##3##Figure 4##, ##SUPPL##3##Table S1##). Yfr15 is located in genomic island 2 in the vicinity of PMED4_07441 (PMM0686), the most highly upregulated host mRNA during phage infection, although the two genes are located on opposite strands. Also PMED4_07401 (PMM0684) and PMED4_07421 (PMM0685), two further genes that belong to group 2 phage-induced host genes, are located nearby in genomic island 2. This region and Yfr15 may therefore be of prime importance for phage-host interactions.</p> <p>We did not detect a single ncRNA that was significantly differentially expressed under phosphorus limitation. This was very surprising, in light of the 34 protein encoding-genes that are differentially expressed under P-stress in MED4 ##REF##16895994##[40]##, and because in <italic>Escherichia coli</italic> the existence of such ncRNAs was hypothesized based on the observation of Hfq-dependent regulation of <italic>rpoS</italic> in response to this stress ##REF##14526009##[46]##.</p> <title>An Ultraconserved Sequence Motif is Present in Three ncRNAs</title> <p>The ncRNA Yfr10 contains the conserved unadecanucleotide motif <named-content content-type="gene">5′-ACUCCUCACAC-3′</named-content> (##FIG##5##Figure 6##). This motif occurs 3 times in the MED4 genome sequence, which is more frequent than expected by chance: One would expect approximately 0.5 instances of a specific 11 nt motif in a 2 MB genome at equal base distribution. The second occurrence has already been described as belonging to another ncRNA in MED4, Yfr1 ##REF##17941988##[27]##, which is also found throughout the cyanobacterial radiation.##REF##17941988##[27]##. Using Northern analysis and 5′RACE, we showed that also the third copy of this motif is expressed, revealing another ncRNA - Yfr18 (##FIG##1##Figure 2##, ##TAB##0##Table 1##). This one was not detected from our microarray analyses nor the comparative genomics approach ##REF##16168080##[25]##.</p> <p>If two base transitions are allowed there is even a fourth member of this ‘unadecanucleotide-containing’ class of ncRNAs in MED4. This ncRNA - Yfr12 - was identified by sequence similarity to Yfr10 and was verified as an ncRNA as described above (##FIG##1##Figure 2##, ##TAB##0##Table 1##). Yfr12 appears to be a mutated variant of the other three, since the processed 5′ end of the major accumulating RNA species was mapped to the middle of the unadecanucleotide, and two mutations change the sequence at the 3′ end of the motif from <named-content content-type="gene">CACACAC</named-content> to <named-content content-type="gene">CAUAUAC</named-content> (##FIG##5##Figure 6##). Furthermore, the motif can be extended in all 4 ncRNAs by another AC dinucleotide (##FIG##5##Figure 6##) in comparison to the published cyanobacterial consensus ##REF##17941988##[27]##, probably a peculiarity of these RNAs in <italic>Prochlorococcus</italic>.</p> <p>The functions of Yfr10, Yfr12 and Yfr18 in MED4 remain unknown at present. However, a hint about their potential function may be found from their genome context and the fact that the vast majority of functional interactions between ncRNAs and their targets is exerted through base pairing. The genes for Yfr10, Yfr12 and Yfr18 are each directly adjacent to those for ncRNAs Yfr2, Yfr4 and Yfr5 respectively (##TAB##0##Table 1##), the 5′ ends of which may basepair to the 13-nucleotide consensus of Yfr1, Yfr10, Yfr12 and Yfr18 if a single bulging C and one mismatch is allowed (i.e. <named-content content-type="gene">5′-aCUCCUcACACAC-3′</named-content> pairs with <named-content content-type="gene">5′-GUGUGUAGGAG-3′</named-content>). Moreover, one may note that also the two C – to - U transitions in Yfr12 are compatible with this suggested base pairing, and that secondary structure predictions suggest that the conserved motifs in Yfr1, Yfr10 and Yfr18 are exposed as single stranded elements in an otherwise folded region (##FIG##5##Figure 6##). The same is true for the complementary motif in Yfr2–Yfr5, making physical interactions very likely. Evidence for two regulatory RNAs acting upon each other has recently been reported for the first time for GlmY and GlmZ of <italic>Escherichia coli</italic>, and cascades of hierarchically acting regulatory RNAs have been hypothesized for other bacteria as well ##REF##18351803##[43]##. The ncRNAs described in this section are candidates for such interactions in <italic>Prochlorococcus</italic>.</p> <title>Two Overlapping Pairs of ncRNAs</title> <p>The difference between <italic>trans</italic>- and <italic>cis</italic>-encoded ncRNAs and asRNAs is frequently considered fuzzy since both act through base complementarity. However, depending on the length of the overlap, interactions between transcripts from the forward and the reverse DNA strand can be very strong due to the extended perfect sequence complementarity. asRNAs may act as the antidote in toxin-antitoxin systems ##REF##12123448##[47]##,##REF##16385044##[48]## or in gene regulation ##REF##16636284##[30]##,##REF##16324715##[31]## as has been reported for some bacteria. We detected two regions with probable sense/antisense pairing between ncRNAs. One of these regions contains Yfr6 with Yfr14 on the opposite strand (##FIG##6##Figure 7##, ##TAB##0##Table 1##). The second region is located in genomic island 1 containing Yfr8 and Yfr9, each of which are 290 nt in size (##FIG##1##Figures 2## and ##FIG##6##7##, ##TAB##0##Table 1##). One of the ncRNAs in both pairs is upregulated during phage infection (##FIG##3##Figure 4##) and both ncRNA pairs contain a potential peptide-coding open reading frame within the sequence of one of the RNAs (##FIG##6##Figure 7##). The peptide sequence associated with Yfr6 is 33 amino acids long and is highly conserved and widely distributed among high light- and low light <italic>Prochlorococcus</italic> strains. The potential 44 amino acids peptide-coding frame within Yfr9 was found in three other <italic>Prochlorococcus</italic> isolates (MIT9515, MIT9301, AS9601), but has been lost in <italic>Prochlorococcus</italic> strains MIT9312 and MIT9215 due to a frame shift. Homologues of Yfr6 and Yfr9, respectively, have high sequence conservation over their complete ncRNA genes – including the upstream and downstream regions of the potential peptides.</p> <p>Yfr6 and Yfr9 resemble RNAIII from <italic>Staphylococcus aureus</italic> being both relatively long and consisting of a small peptide-coding unit as well as a regulatory RNA. RNAIII is a 510 nt long riboregulator from which the 26 amino acid δ-hemolysin peptide is also translated ##REF##9620969##[49]##. Another bifunctional ncRNA (SgrS) has been described in <italic>Escherichia coli</italic> that contains a conserved ORF (SgrT) in the 5′ region of SgrS, both of which promote recovery from glucose stress in mechanistically distinct fashions ##REF##18042713##[50]##. Moreover, secondary structure predictions of Yfr6 ##REF##16168080##[25]## and of Yfr14, Yfr8, and Yfr9 (##SUPPL##1##Figures S2## and ##SUPPL##2##S3##) support the potential role of these transcripts as functional RNAs as they contain many G - C base pairings and compensatory mutations, which conserve the structure rather than the sequence – a feature of many ncRNAs. On the other hand, highly structured transcript regions are also found in certain mRNAs where they serve as platforms for sophisticated ncRNA-mediated translational control. In the case of the <italic>Escherichia coli tisB</italic> mRNA, for example, this transcript encodes a peptide as short as 29 amino acids towards its 3′ end, yet there is no evidence that <italic>tisB</italic> would act as a riboregulator ##REF##11448770##[17]##,##REF##17499044##[51]##.</p> <p>The fact that Yfr6 and Yfr9 overlap with other ncRNAs (Yfr14 and Yfr8 respectively), suggests that this could be a toxin-antitoxin system – i.e. pairs of genes that code for a stable toxin and an unstable antitoxin. These are well-characterized in other bacteria, where the toxin is usually a toxic peptide that is neutralized or whose synthesis is prevented by the action of the product of the second gene, the antitoxin, which is either protein or RNA. Toxin-antitoxin systems such as the <italic>hok/sok</italic> system of <italic>Escherichia coli</italic> can serve as a natural genetic selection system to ensure presence of a plasmid ##UREF##3##[52]##. Alternatively, chromosomally encoded toxin-antitoxin systems can be beneficial to cell survival under unfavorable growth conditions, sometimes in very sophisticated ways, for instance by transiently curtailing the consumption of nutrients during starvation or by temporarily inhibiting growth and thereby evading the killing effects of certain antibiotics ##REF##15864262##[53]##. Systematic searches for toxin-antitoxin systems have revealed a high abundance in free-living prokaryotes ##REF##15718296##[54]## but none of the seven known toxin-antitoxin families could be identified in <italic>Prochlorococcus</italic>. There is a growing number of examples of chromosomal toxin-antitoxin systems that use a <italic>cis</italic>-encoded asRNA as an antitoxin. Our data suggest that Yfr6/Yfr14 and Yfr8/Yfr9 may be potential candidates for toxin-antitoxin systems in <italic>Prochlorococcus</italic> MED4.</p> <title>A High Number of <italic>cis</italic>-Encoded Antisense RNAs</title> <p>Little attention has been given to chromosomally <italic>cis</italic>-encoded asRNAs until recently, and only a few have been described for cyanobacteria ##REF##2450053##[29]##–##REF##16324715##[31]##. Surprisingly, we detected 24 asRNAs in our analyses, which vary between 100 to 600 nt in size (##SUPPL##4##Table S2##). Some are differentially expressed under different light conditions and under phage infection (##SUPPL##5##Table S3##).</p> <p>High light treatment caused one asRNA to be upregulated and one to be downregulated (##SUPPL##5##Table S3##). asRNA asMED4_15721 was upregulated twofold when cells were transferred from darkness to high light (##SUPPL##5##Table S3##). This behavior is similar to that of Yfr20 (see above), and like Yfr20, this asRNA is located in genomic island 5, lending additional support for a function of this island in light stress adaptation. Notably, five of the 24 asRNAs are complementary to mRNAs that code for photosystem I subunits (<italic>psaB</italic> and <italic>psaC</italic>) or for photosystem II subunits (<italic>psbB</italic> and <italic>psbO</italic>, <italic>psbX</italic>), respectively. The concentrations of asRNAs of photosystem II genes did not change when cultures were shifted from darkness to different light quantities and qualities, whereas transcript levels of their target mRNA decreased slightly (##SUPPL##5##Table S3##). In contrast, levels of the photosystem I asRNAs asMED4_17331 (antisense of <italic>psaB</italic>) and asMED4_18171 (antisense of <italic>psaC</italic>) decreased when cells were shifted from darkness to light, following the same trend as their mRNA counterparts (##SUPPL##5##Table S3##). Surprisingly, however, the latter asRNAs decreased in amount when transferred from darkness to medium white light whereas their target mRNAs did not (##SUPPL##5##Table S3##), which might indicate a light-dosage specific regulation of these asRNAs.</p> <p>We also found asRNAs that are differentially expressed in cells infected by phage. asMED4_04601 is upregulated during the initial stages of phage infection, whereas its target mRNA (PMED4_04601) is constitutively expressed throughout the infection process. Interestingly, PMED4_04601 shows 67% amino acid identity to the central region of the potential Yfr6 peptide. In the case of PMED4_07401 (PMM0684) both the respective asRNA (asMED4_07401) and its target mRNA are upregulated from mid-to-late phase of phage infection.</p> <p>The number of 24 asRNAs detected in our analyses appears high, especially as the microarrays used for this study did not contain probes for antisense regions of protein coding genes (see <xref ref-type="sec" rid="s3">methods</xref>). Therefore we could only detect those asRNAs found in intergenic regions whose corresponding ORF was not originally annotated and asRNAs located in 5′ and 3′ UTRs. However, our mapping results revealed that asRNAs located opposite of 5′/3′ UTRs frequently overlap major parts of the adjacent coding sequences. While it is not possible to infer the functions of these asRNAs in <italic>Prochlorococus</italic>, in the cyanobacterium <italic>Synechocystis</italic> PCC6803 the asRNA IsrR occurs in higher quantities than the <italic>cis</italic>-encoded mRNA <italic>isiA</italic> under normal growth conditions, leading to degradation of RNA duplexes by RNase III ##REF##16636284##[30]##. Under stress conditions the expression of the mRNA is increased leading to free mRNA molecules that can be translated. This mode of action is highly unlikely in the case of the PMED4_07401 mRNA/asRNA pair, because they are co-upregulated during phage infection, pointing towards a protective rather than a degradative role. Interestingly, RNAse E is also among the genes upregulated during phage infection. It is hypothesized that ribonuclease activity could be utilized by the phage to degrade host RNA to generate nucleotides for phage replication ##REF##17805294##[39]##.</p> <title>Other Overlapping Transcripts</title> <p>In addition to non-protein coding asRNAs, we found several pairs of transcripts that are transcribed from complementary strands and that potentially code for proteins. These complementary transcripts overlap entirely with each other or with a major part of their 5′ or 3′ UTR. We confirmed these overlapping regions experimentally and found that they span between 74 nt to 333 nt at least (##SUPPL##4##Table S2##). Eight out of twelve of these overlapping regions are found in the same position in other <italic>Prochlorococcus</italic> genomes, whereas the other 4 are found in different regions of the genomes. At this point it is not clear whether the overlaps between these protein-encoding transcripts would interfere with their transcription, transcript accumulation or translation. We found evidence for both scenarios. Whereas PMED4_14671 (located in the opposite 3′UTR region of PMED4_14661) is upregulated 14 fold when light intensity is increased, PMED4_14661 (PMM1300) remains at basal transcript levels. Contrary to the above, PMED4_11211 (PMM0997) and PMED4_11201 (located in the opposite 3′UTR region of PMED4_11211) are inversely regulated under different light conditions and DCMU treatment (##SUPPL##6##Table S4##) indicative of either interference during transcription, or coupled degradation, as observed for the asRNA and mRNA IsrR/<italic>isiA</italic> in the cyanobacterium <italic>Synechocystis</italic> PCC6803 ##REF##16636284##[30]##.</p> <p>In cases where only 3′-UTRs overlap this might not be of relevance because the transcriptional machinery should not be constricted. However, most of the transcripts we found overlap the 5′-UTRs and/or complete protein-coding regions and therefore are highly likely to be of regulatory relevance.</p> <title>Identification of New ORFs</title> <p>It is difficult to identify short genes in bacterial genomes using annotation algorithms, because the number of possible reading frames increases the shorter the search window becomes. Therefore, many of the widely used annotation programs (e.g. GLIMMER, GeneMark and CRITICA) constrain the minimum length of an ORF and thus a considerable number of small ORFs remain unannotated. In a recent study of the <italic>Prochlorococcus</italic> pan genome, for example, hypothetical ORFs shorter than 50 amino acids were excluded unless they were found in more than one genome ##REF##18159947##[9]##.</p> <p>Our microarray analyses lead to the observation of 113 new ORFs (not including asRNAs with potential protein coding sequences; ##SUPPL##4##Table S2##) ranging between 33 to 130 amino acids in size that were not annotated in the first published <italic>Prochlorococcus</italic> MED4 genome version (accession number BX548174). The new annotation of Kettler et al. ##REF##18159947##[9]## (accession number BX548174.1, for new ORF IDs refer to: <ext-link ext-link-type="uri" xlink:href="http://www.microbesonline.org">www.microbesonline.org</ext-link>) also found 89 of the new ORFs. BlastP searches against the non-redundant NCBI database revealed that 10 of the 24 remaining ORFs have an annotated counterpart in at least one other genome whereas 14 ORFs represent short proteins that have not been detected thus far (##SUPPL##7##Table S5##). Using TblastN, all of the additional 14 novel ORFs were found in other genomes, even though they had gone undetected by computational annotation tools (##SUPPL##8##Table S6##).</p> <title>Phage-Related Open Reading Frames</title> <p>Amongst the newly discovered protein-coding genes are three ORFs that have homologues in cyanophage genomes. PMED4_16122, which is located in genomic island 5 in MED4, is homologous to PSSM4_095 in the MED4-infecting cyanophage PSSM4 ##REF##12944965##[55]##, suggesting gene transfer between a PSSM4-like phage and this island, in particular since no other homologues were found in the nr database. PMED4_15491 is in vicinity of genomic island 5 in the host genome and has one homolog in each of two cyanophage genomes – P-SSM4 and P-SSM2 (PSSM4_181 and PSSM2_278). According to TblastN results this ORF is present in numerous <italic>Prochlorococcus</italic> genomes but not in other cyanobacteria. The third newly discovered ORF with a homolog in a cyanophage is PMED4_10681, which is found in the genome of cyanophage P-SSM2. Unlike PMED4_16122, however, this ORF is not in a genomic island in the host genome, and furthermore, is widely distributed over the cyanobacterial radiation with homologues in all <italic>Prochlorococcus</italic> strains, <italic>Synechococcus elongatus</italic> strains PCC 6301 and PCC 7942, <italic>Fremyella diplosiphon</italic>, <italic>Nostoc</italic> sp. PCC 7120, <italic>Anabaena variabilis</italic> and <italic>Synechocystis</italic> sp. PCC 6803 (##SUPPL##8##Table S6##). The broad distribution of PMED4_10681 suggests that it plays an important function in cyanobacteria, and emphasizes the importance of better annotation of small ORFs.</p> <title>Conclusions</title> <p>Here we have described 14 novel ncRNAs, which increases the total number of ncRNAs in this organism to 24 (including Yfr1-7, ffs, tmRNA, RNase P RNA). One sixth of the 24 ncRNAs (Yfr1, Yfr3, Yfr12 and Yfr18) were undetectable from microarray analyses under the conditions tested. Therefore it is likely that even more ncRNAs are present in <italic>Prochlorococcus</italic> MED4. The proportion of ncRNAs in the <italic>Prochlorococcus</italic> MED4 genome is comparable with those found in enterobacteria like <italic>Escherichia coli</italic>, i.e. 1–2% of the genes encode ncRNAs. In comparison, the 6 identified protein regulators in <italic>Prochlorococcus</italic>
##REF##12917642##[10]## is a small number relative to the 32 two-component response regulators present in <italic>Escherichia coli</italic>
##REF##9205844##[56]##. This suggests that regulation of gene expression through ncRNAs plays an important role in <italic>Prochlorococcus</italic>' response to environmental cues. The relatively high number of ncRNAs is intriguing as it may represent a mode of adaptation to the extremely low nutrient conditions of the open oceans. Regulation by ncRNAs may require fewer resources than would be required for the synthesis of protein regulators. Furthermore, in the course of genome reduction there might have been a positive selection pressure for keeping small regulators, e.g. ncRNAs rather than large protein regulators. How ncRNAs function in <italic>Prochlorococcus</italic> is at present unclear. The absence of Hfq in MED4 suggests that the ncRNAs found in this strain represent a core-set of ncRNAs that function without the support of a chaperone, or with a novel chaperone yet to be identified.</p> <p>The genomic islands of <italic>Prochlorococcus</italic> are disproportionately connected to ecological functions in this group of cyanobacteria ##REF##17693088##[57]##,##REF##16556843##[11]##. Here we have shown that approximately half of the <italic>Prochlorococcus</italic> ncRNAs are located in genomic islands suggesting that the function of these molecules is relevant for determining the relative fitness of ecotypes within <italic>Prochlorococcus</italic>. This is analogous to the accumulation of genes coding for ncRNA in pathogenicity islands in <italic>Staphylococcus aureus</italic>
##REF##16183745##[58]## and <italic>Salmonella typhimurium</italic>
##REF##17971080##[59]## as well as in genomic islands of <italic>Sinorhizobium meliloti</italic>
##REF##18093320##[60]##, and suggests that this phenomenon could be wide-spread for finely tuned specialization within microbial groups.</p>
<title>Conclusions</title> <p>Here we have described 14 novel ncRNAs, which increases the total number of ncRNAs in this organism to 24 (including Yfr1-7, ffs, tmRNA, RNase P RNA). One sixth of the 24 ncRNAs (Yfr1, Yfr3, Yfr12 and Yfr18) were undetectable from microarray analyses under the conditions tested. Therefore it is likely that even more ncRNAs are present in <italic>Prochlorococcus</italic> MED4. The proportion of ncRNAs in the <italic>Prochlorococcus</italic> MED4 genome is comparable with those found in enterobacteria like <italic>Escherichia coli</italic>, i.e. 1–2% of the genes encode ncRNAs. In comparison, the 6 identified protein regulators in <italic>Prochlorococcus</italic>
##REF##12917642##[10]## is a small number relative to the 32 two-component response regulators present in <italic>Escherichia coli</italic>
##REF##9205844##[56]##. This suggests that regulation of gene expression through ncRNAs plays an important role in <italic>Prochlorococcus</italic>' response to environmental cues. The relatively high number of ncRNAs is intriguing as it may represent a mode of adaptation to the extremely low nutrient conditions of the open oceans. Regulation by ncRNAs may require fewer resources than would be required for the synthesis of protein regulators. Furthermore, in the course of genome reduction there might have been a positive selection pressure for keeping small regulators, e.g. ncRNAs rather than large protein regulators. How ncRNAs function in <italic>Prochlorococcus</italic> is at present unclear. The absence of Hfq in MED4 suggests that the ncRNAs found in this strain represent a core-set of ncRNAs that function without the support of a chaperone, or with a novel chaperone yet to be identified.</p> <p>The genomic islands of <italic>Prochlorococcus</italic> are disproportionately connected to ecological functions in this group of cyanobacteria ##REF##17693088##[57]##,##REF##16556843##[11]##. Here we have shown that approximately half of the <italic>Prochlorococcus</italic> ncRNAs are located in genomic islands suggesting that the function of these molecules is relevant for determining the relative fitness of ecotypes within <italic>Prochlorococcus</italic>. This is analogous to the accumulation of genes coding for ncRNA in pathogenicity islands in <italic>Staphylococcus aureus</italic>
##REF##16183745##[58]## and <italic>Salmonella typhimurium</italic>
##REF##17971080##[59]## as well as in genomic islands of <italic>Sinorhizobium meliloti</italic>
##REF##18093320##[60]##, and suggests that this phenomenon could be wide-spread for finely tuned specialization within microbial groups.</p>
<p>Conceived and designed the experiments: CS. Performed the experiments: CS. Analyzed the data: CS MEF BV. Contributed reagents/materials/analysis tools: MEF WRH. Wrote the paper: CS DL SWC WRH.</p> <p>
<italic>Prochlorococcus</italic>, an extremely small cyanobacterium that is very abundant in the world's oceans, has a very streamlined genome. On average, these cells have about 2,000 genes and very few regulatory proteins. The limited capability of regulation is thought to be a result of selection imposed by a relatively stable environment in combination with a very small genome. Furthermore, only ten non-coding RNAs (ncRNAs), which play crucial regulatory roles in all forms of life, have been described in <italic>Prochlorococcus</italic>. Most strains also lack the RNA chaperone Hfq, raising the question of how important this mode of regulation is for these cells. To explore this question, we examined the transcription of intergenic regions of <italic>Prochlorococcus</italic> MED4 cells subjected to a number of different stress conditions: changes in light qualities and quantities, phage infection, or phosphorus starvation. Analysis of Affymetrix microarray expression data from intergenic regions revealed 276 novel transcriptional units. Among these were 12 new ncRNAs, 24 antisense RNAs (asRNAs), as well as 113 short mRNAs. Two additional ncRNAs were identified by homology, and all 14 new ncRNAs were independently verified by Northern hybridization and 5′RACE. Unlike its reduced suite of regulatory proteins, the number of ncRNAs relative to genome size in <italic>Prochlorococcus</italic> is comparable to that found in other bacteria, suggesting that RNA regulators likely play a major role in regulation in this group. Moreover, the ncRNAs are concentrated in previously identified genomic islands, which carry genes of significance to the ecology of this organism, many of which are not of cyanobacterial origin. Expression profiles of some of these ncRNAs suggest involvement in light stress adaptation and/or the response to phage infection consistent with their location in the hypervariable genomic islands.</p> <title>Author Summary</title> <p>
<italic>Prochlorococcus</italic> 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 <italic>Prochlorococcus</italic>—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 <italic>Prochlorococcus</italic> to particular environmental challenges.</p>
<title>Supporting Information</title>
<fig id="pgen-1000173-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000173.g001</object-id><label>Figure 1</label><caption><title>Computational pipeline for the identification of novel genetic elements from microarray data.</title><p>The data were derived from experiments studying stress of <italic>Prochlorococcus</italic> MED4 induced by light changes, phosphorus starvation, and phage infection. Note that the sum of classified transcripts is unequal to the number of transcriptional units due to the occurrence of dicistronic elements. In addition to the 12 ncRNAs newly detected by microarray analysis, two ncRNAs Yfr12 and Yfr18 were found based on similarity.</p></caption></fig> <fig id="pgen-1000173-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000173.g002</object-id><label>Figure 2</label><caption><title>Detection of novel ncRNAs by Northern hybridization.</title><p>The lengths of hybridizing fragments were calculated by comparison to two different size markers (M) – size in nt is shown to the left of the markers. 50 µg of total RNA was loaded per lane on a 10% polyacrylamide gel that was electro-blotted on Hybond-N nylon membranes. All ncRNAs were detected with single strand RNA probes except for those labeled with an asterisk, which were hybridized to an oligonucleotide probe.</p></caption></fig> <fig id="pgen-1000173-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000173.g003</object-id><label>Figure 3</label><caption><title>Genomic location of Yfr13 (red boxes) in seven different <italic>Prochlorococcus</italic> strains.</title><p>Homologous genes in different strains are indicated by blue boxes and those without homology in this genomic region are shown in white. Gene designations are given according to published genome sequences for MED4 (BX548174.1), MIT9312 (CP000111.1), MIT9515 (CP000552.1), MIT9215 (CP000825), MIT9301 (CP000576.1), NATL1A (CP000553.1), NATL2A (CP000095.2).</p></caption></fig> <fig id="pgen-1000173-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000173.g004</object-id><label>Figure 4</label><caption><title>Regulation of ncRNAs under light and phage stress (no changes were found under P-stress).</title><p>Fold changes (log2) of triplicate experiments are shown. For the light experiment grey vertical lines indicate the border of ≥2-fold change after 45 min of light treatment. Only ncRNAs with significant differential expression and their oppositely located ncRNA pair (when present) are illustrated for the phage experiment. Standard errors derive from experimental repetitions of all probes that gave a signal above threshold.</p></caption></fig> <fig id="pgen-1000173-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000173.g005</object-id><label>Figure 5</label><caption><title>Genomic location of Yfr20 (red boxes) in three high-light adapted <italic>Prochlorococcus</italic> strains.</title><p>Homologous genes are indicated by the same colors in the different strains and those without homology in this genomic region are shown in white. Gene designations are given according to published genome sequences for MED4 (BX548174.1), MIT9312 (CP000111.1) and MIT9515 (CP000552.1).</p></caption></fig> <fig id="pgen-1000173-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000173.g006</object-id><label>Figure 6</label><caption><title>A highly conserved RNA motif.</title><p>The conserved motif (<named-content content-type="gene">5′-ACUCCUCACAC-3′</named-content>) is highlighted in blue and the <italic>Prochlorococcus</italic>–specific extension of the motif is indicated by green letters. The two base mutations in Yfr12 are colored in purple. The color code for shaded bases is as follows: red, all sequences have the same two nucleotides; ochre, two types of base pairs occur; green, three types of base pairs occur; turquoise, four types of base pairs occur. The saturation decreases with the number of sequences unable to form a base pair at this position. Circles around letters indicate compensatory mutations. Consensus structures for Yfr1, Yfr10, Yfr18 were predicted with RNAlishapes ##REF##17020924##[63]## and drawn with RNAplot ##UREF##5##[64]##. Homologous ncRNAs found in other <italic>Prochlorococcus</italic> strains that were used for structure predictions are listed in parentheses. The structure of Yfr12 was folded using RNAshapes and relies on the single sequence of Yfr12 from MED4.</p></caption></fig> <fig id="pgen-1000173-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000173.g007</object-id><label>Figure 7</label><caption><title>Gene arrangement of ncRNA pairs Yfr8/Yfr9 and Yfr6/Yfr14.</title><p>Potential peptide-coding open reading frames within the sequence of the RNAs are shown as white boxes. Positions in the genome are indicated by numbers for the forward and reverse (c) strand.</p></caption></fig>
<table-wrap id="pgen-1000173-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pgen.1000173.t001</object-id><label>Table 1</label><caption><title>Synopsis of ncRNAs in <italic>Prochlorococcus</italic> MED4. Newly identified ncRNAs are based on microarray expression data or on sequence similarity (Yfr12 and Yfr18) according to the scheme shown in ##FIG##0##Figure 1##.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">probe names</td><td align="left" rowspan="1" colspan="1">new name</td><td align="left" rowspan="1" colspan="1">genomic location</td><td align="left" rowspan="1" colspan="1">length in nt<xref ref-type="table-fn" rid="nt102">a</xref>
</td><td align="left" rowspan="1" colspan="1">adjacent genes</td><td align="left" rowspan="1" colspan="1">genomic context</td><td align="left" rowspan="1" colspan="1">occurence in <italic>Prochlorococcus</italic>
</td><td align="left" rowspan="1" colspan="1">Genomic island</td><td align="left" rowspan="1" colspan="1">mean MFE z score<xref ref-type="table-fn" rid="nt103">b</xref>
</td><td align="left" rowspan="1" colspan="1">SVM RNA-class probability<xref ref-type="table-fn" rid="nt103">b</xref>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">IG0367f11-17</td><td align="left" rowspan="1" colspan="1">Yfr8</td><td align="left" rowspan="1" colspan="1">339952..340241</td><td align="left" rowspan="1" colspan="1">180, 220, 290</td><td align="left" rowspan="1" colspan="1">PMED4_03781 and PMED4_03791 (PMM0356)</td><td align="left" rowspan="1" colspan="1">→ → →</td><td align="left" rowspan="1" colspan="1">MIT9301, MIT9515, MIT9312, AS9601, MIT9215</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">−4.64</td><td align="left" rowspan="1" colspan="1">0.999616</td></tr><tr><td align="left" rowspan="1" colspan="1">IG0367r2-5</td><td align="left" rowspan="1" colspan="1">Yfr9</td><td align="left" rowspan="1" colspan="1">complement(339887..340176)</td><td align="left" rowspan="1" colspan="1">180, 220, 290</td><td align="left" rowspan="1" colspan="1">PMED4_03781 and PMED4_03791 (PMM0356)</td><td align="left" rowspan="1" colspan="1">→ ← →</td><td align="left" rowspan="1" colspan="1">MIT9301, MIT9515, MIT9312, AS9601, MIT9215</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">−2.54</td><td align="left" rowspan="1" colspan="1">0.676448</td></tr><tr><td align="left" rowspan="1" colspan="1">IG0375f18-20</td><td align="left" rowspan="1" colspan="1">Yfr10</td><td align="left" rowspan="1" colspan="1">346661..346738</td><td align="left" rowspan="1" colspan="1">78, 270<xref ref-type="table-fn" rid="nt106">e</xref>
</td><td align="left" rowspan="1" colspan="1">PMED4_03891 and Yfr2</td><td align="left" rowspan="1" colspan="1">← → →</td><td align="left" rowspan="1" colspan="1">MIT9301, MIT9312</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">−4.69<xref ref-type="table-fn" rid="nt111">j</xref>
</td><td align="left" rowspan="1" colspan="1">0.991539<xref ref-type="table-fn" rid="nt111">j</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">IG0390r20-22</td><td align="left" rowspan="1" colspan="1">Yfr11</td><td align="left" rowspan="1" colspan="1">complement(358606..358754)</td><td align="left" rowspan="1" colspan="1">149</td><td align="left" rowspan="1" colspan="1">PMED4_04111 and PMED4_04121</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">(PMM0378)</td><td align="left" rowspan="1" colspan="1">← ← ←</td><td align="left" rowspan="1" colspan="1">MIT9515</td><td align="left" rowspan="1" colspan="1">1</td><td align="left" rowspan="1" colspan="1">−2.38<xref ref-type="table-fn" rid="nt107">f</xref>
</td><td align="left" rowspan="1" colspan="1">0.998347<xref ref-type="table-fn" rid="nt107">f</xref>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">Yfr12</td><td align="left" rowspan="1" colspan="1">383230..383312</td><td align="left" rowspan="1" colspan="1">83</td><td align="left" rowspan="1" colspan="1">PMED4_04471 and Yfr4</td><td align="left" rowspan="1" colspan="1">→ → →</td><td align="left" rowspan="1" colspan="1">MIT9301, MIT9312</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">−3.09<xref ref-type="table-fn" rid="nt111">j</xref>
</td><td align="left" rowspan="1" colspan="1">0.996117<xref ref-type="table-fn" rid="nt111">j</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">IG0668r5</td><td align="left" rowspan="1" colspan="1">Yfr13</td><td align="left" rowspan="1" colspan="1">complement(621601..621701)</td><td align="left" rowspan="1" colspan="1">101</td><td align="left" rowspan="1" colspan="1">PMED4_07021 and PMED4_07031</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">(PMM0652)</td><td align="left" rowspan="1" colspan="1">← ← ←</td><td align="left" rowspan="1" colspan="1">MIT9301, MIT9515, MIT9312, MIT9215, NATL1a, NATL2a</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">−1.78<xref ref-type="table-fn" rid="nt108">g</xref>
</td><td align="left" rowspan="1" colspan="1">0.600913<xref ref-type="table-fn" rid="nt108">g</xref>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">IG0676f22-24</td><td align="left" rowspan="1" colspan="1">Yfr14</td><td align="left" rowspan="1" colspan="1">627811..628062</td><td align="left" rowspan="1" colspan="1">252</td><td align="left" rowspan="1" colspan="1">PMED4_07142 and PMED4_07151</td><td align="left" rowspan="1" colspan="1">→ → →</td><td align="left" rowspan="1" colspan="1">n.d.</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">n.d.</td><td align="left" rowspan="1" colspan="1">n.d.</td></tr><tr><td align="left" rowspan="1" colspan="1">IG0702r18-21</td><td align="left" rowspan="1" colspan="1">Yfr15</td><td align="left" rowspan="1" colspan="1">complement(654301..654562)</td><td align="left" rowspan="1" colspan="1">86, 262</td><td align="left" rowspan="1" colspan="1">PMED4_07441 and Yfr3</td><td align="left" rowspan="1" colspan="1">→ ← →</td><td align="left" rowspan="1" colspan="1">n.d.</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">n.d.</td><td align="left" rowspan="1" colspan="1">n.d.</td></tr><tr><td align="left" rowspan="1" colspan="1">IG0704r24-25</td><td align="left" rowspan="1" colspan="1">Yfr16</td><td align="left" rowspan="1" colspan="1">complement(657392..657529)</td><td align="left" rowspan="1" colspan="1">138</td><td align="left" rowspan="1" colspan="1">PMED4_07521 and PMED4_07522</td><td align="left" rowspan="1" colspan="1">← ← →</td><td align="left" rowspan="1" colspan="1">MIT9515</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">−2.38<xref ref-type="table-fn" rid="nt107">f</xref>
</td><td align="left" rowspan="1" colspan="1">0.998347<xref ref-type="table-fn" rid="nt107">f</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">IG0714f13</td><td align="left" rowspan="1" colspan="1">Yfr17</td><td align="left" rowspan="1" colspan="1">664667..664930<xref ref-type="table-fn" rid="nt105">d</xref>
</td><td align="left" rowspan="1" colspan="1">76, 285, 450</td><td align="left" rowspan="1" colspan="1">PMED4_07691 and PMED4_07701 (PMM0699)</td><td align="left" rowspan="1" colspan="1">← → →</td><td align="left" rowspan="1" colspan="1">MIT9301, MIT9515, MIT9312,AS9601, MIT9215</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">−1.05<xref ref-type="table-fn" rid="nt109">h</xref>
</td><td align="left" rowspan="1" colspan="1">0.000000<xref ref-type="table-fn" rid="nt109">h</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">Yfr18</td><td align="left" rowspan="1" colspan="1">complement(972284..972388)</td><td align="left" rowspan="1" colspan="1">29, 76</td><td align="left" rowspan="1" colspan="1">Yfr5 and PMED4_11681</td><td align="left" rowspan="1" colspan="1">← ← →</td><td align="left" rowspan="1" colspan="1">MIT9301, MIT9515, MIT9312, AS9601, MIT9215</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">−0.92<xref ref-type="table-fn" rid="nt110">i</xref>
</td><td align="left" rowspan="1" colspan="1">0.996639<xref ref-type="table-fn" rid="nt110">i</xref>
</td></tr><tr><td align="left" rowspan="1" colspan="1">IG1142r19-25</td><td align="left" rowspan="1" colspan="1">Yfr21<xref ref-type="table-fn" rid="nt104">c</xref>
</td><td align="left" rowspan="1" colspan="1">complement(1069553..1069757)</td><td align="left" rowspan="1" colspan="1">205</td><td align="left" rowspan="1" colspan="1">PMED4_12671 and PMED4_12681</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">(PMM1121)</td><td align="left" rowspan="1" colspan="1">← ← ←</td><td align="left" rowspan="1" colspan="1">no homologues</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">n.d.</td><td align="left" rowspan="1" colspan="1">n.d.</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">IG1157r8</td><td align="left" rowspan="1" colspan="1">Yfr19</td><td align="left" rowspan="1" colspan="1">complement(1086081..1086162)</td><td align="left" rowspan="1" colspan="1">82</td><td align="left" rowspan="1" colspan="1">PMED4_12891 and PMED4_12901</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">(PMM1135)</td><td align="left" rowspan="1" colspan="1">→ ← →</td><td align="left" rowspan="1" colspan="1">MIT9301, MIT9312, AS9601, MIT9215</td><td align="left" rowspan="1" colspan="1">-</td><td align="left" rowspan="1" colspan="1">−3.89<xref ref-type="table-fn" rid="nt108">g</xref>
</td><td align="left" rowspan="1" colspan="1">0.999504<xref ref-type="table-fn" rid="nt108">g</xref>
</td><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1"/></tr><tr><td align="left" rowspan="1" colspan="1">IG1412f30</td><td align="left" rowspan="1" colspan="1">Yfr20</td><td align="left" rowspan="1" colspan="1">1336435..1336523</td><td align="left" rowspan="1" colspan="1">89</td><td align="left" rowspan="1" colspan="1">PMED4_15791 and PMED4_15792</td><td align="left" rowspan="1" colspan="1">→ → →</td><td align="left" rowspan="1" colspan="1">MIT9515, MIT9312</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">−0.77</td><td align="left" rowspan="1" colspan="1">0.142982</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pgen.1000173.s001"><label>Figure S1</label><caption><p>Structures of Yfr11 and Yfr16. The structures were predicted in RNAlishapes ##UREF##4##[61]## and were drawn with RNAViz ##UREF##5##[64]##. Conserved bases are shaded in grey.</p><p>(1.41 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s002"><label>Figure S2</label><caption><p>Consensus structures of Yfr6 and Yfr14. The structures were predicted with RNAlishapes ##UREF##4##[61]## and drawn with RNAplot ##REF##8755555##[62]## and are based on sequence alignments including MED4, MIT9312, MIT9515, MIT9301, AS9601, NATL1A and SS120 (in parenthesis). The color code is the same as for ##FIG##5##Figure 6##.</p><p>(1.43 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s003"><label>Figure S3</label><caption><p>Consensus structures of Yfr8 and Yfr9. The structures were predicted with RNAlishapes ##UREF##4##[61]## and drawn with RNAplot ##REF##8755555##[62]## and are based on sequence alignments including MED4, MIT9312, MIT9515, MIT9301, MIT9215 and AS9601 (in parenthesis). The color code is the same as for ##FIG##5##Figure 6##.</p><p>(2.66 MB EPS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s004"><label>Table S1</label><caption><p>Differential expression of ncRNAs during light, phage and phosphorus stress. Microarray ratios are average values ± standard errors of all probes (with expression values above threshold) targeting the respective ncRNA for triplicate biological repeats.</p><p>(0.07 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s005"><label>Table S2</label><caption><p>Transcripts detected from microarray expression analysis classified as new ncRNAs, asRNAs, small ORFs, pseudogenes and UTR regions of adjacent genes. Capital letters and underlined letters in columns “Sequence 5′end” and “Sequence 3′end” indicate mapped ends and TATA boxes, respectively. Transcript detection included microarray analysis (M), 5′ and 3′ RACE (5′R and 3′R), and Northern hybridization (N). For ORFs that were verified by 5′RACE numbers given in brackets correspond to coordinates of the CDS and numbers without brackets indicate mapped 5′ and 3′ ends.</p><p>(0.07 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s006"><label>Table S3</label><caption><p>Differential expression of asRNAs and their sense transcripts during light, phage and phosphorus stress. asRNAs were named with the prefix “as” according to the gene ID found at <ext-link ext-link-type="uri" xlink:href="http://www.microbesonline.org">www.microbesonline.org</ext-link> for the respective gene they are overlapping with. Microarray ratios are average values ± standard errors of all probes targeting the respective asRNA and its sense transcript for biological triplicates.</p><p>(0.10 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s007"><label>Table S4</label><caption><p>Differential expression of overlapping transcripts during light, phage and phosphorus stress. Microarray ratios are average values ± standard errors of all probes targeting the respective transcript for biological triplicates.</p><p>(0.06 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s008"><label>Table S5</label><caption><p>BlastP results for new protein-coding genes against the non-redundant NCBI database. Blast search criteria were set to a cut-off E-value below 0.05 using default algorithm parameters except for compositional adjustment that was disabled. Blasts were performed in February 2008 against the NCBI non-redundant database.</p><p>(0.20 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s009"><label>Table S6</label><caption><p>TblastN results for new protein-coding genes against the non-redundant NCBI database. Blast search criteria were set to a cut-off E-value below 0.05 using default algorithm parameters except for compositional adjustment that was disabled. Blasts were performed in February 2008 against the NCBI non-redundant database.</p><p>(0.08 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s010"><label>Table S7</label><caption><p>List of used oligonucleotides.</p><p>(0.04 MB XLS)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pgen.1000173.s011"><label>Table S8</label><caption><p>Mapped 5′ends of ORFs opposite of asRNAs that were not detected by the microarray screen.</p><p>(0.02 MB XLS)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><p>Microarray probe IDs are given together with the respective ncRNA name, location and length. Computational support values (mean MFE z score and SVM RNA class probability) have been calculated according to Washietl et al. ##REF##15665081##[65]##.</p></fn><fn id="nt102"><label>a</label><p>length determined from Northern results and verified by 5′RACE.</p></fn><fn id="nt103"><label>b</label><p>values derived from RNAz.</p></fn><fn id="nt104"><label>c</label><p>not detected by Northern blot.</p></fn><fn id="nt105"><label>d</label><p>length and location determined from 5′RACE results.</p></fn><fn id="nt106"><label>e</label><p>precursor of Yfr10 and Yfr2.</p></fn><fn id="nt107"><label>f</label><p>values derived from consensus structure of Yfr11, Yfr16 and a single homologue in MIT9515.</p></fn><fn id="nt108"><label>g</label><p>values without MIT9301, MIT9312 and MIT9215.</p></fn><fn id="nt109"><label>h</label><p>values derived from 264 nt long RNA as determined by 5′RACE; without AS9601, MIT9301 and MIT9215.</p></fn><fn id="nt110"><label>i</label><p>values without MIT9301 and MIT9215.</p></fn><fn id="nt111"><label>j</label><p>values derived using the same MIT9301 and MIT9312 homologues.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p>This work was supported in part by grants from the Deutsche Forschungsgemeinschaft SPP 1258 “Sensory and regulatory RNAs in prokaryotes” (CS), the Morasha program of the Israel Science Foundation Grant no. 1504/06 (DL), the NSF, DOE- GTL, the Gordon and Betty Moore Foundation (SWC), and the BMBF - Freiburg Initiative in Systems Biology (WRH). DL is a Shillman Fellow.</p></fn></fn-group>
<graphic xlink:href="pgen.1000173.g001"/> <graphic xlink:href="pgen.1000173.g002"/> <graphic id="pgen-1000173-t001-1" xlink:href="pgen.1000173.t001"/> <graphic xlink:href="pgen.1000173.g003"/> <graphic xlink:href="pgen.1000173.g004"/> <graphic xlink:href="pgen.1000173.g005"/> <graphic xlink:href="pgen.1000173.g006"/> <graphic xlink:href="pgen.1000173.g007"/>
<media xlink:href="pgen.1000173.s001.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s002.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s003.eps"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s004.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s005.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s006.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s007.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s008.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s009.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s010.xls"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pgen.1000173.s011.xls"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["1"], "element-citation": ["\n"], "surname": ["Goericke", "Welschmeyer"], "given-names": ["R", "NA"], "year": ["1993"], "article-title": ["The marine prochlorophyte "], "italic": ["Prochlorococcus"], "source": ["Deep-Sea Res"], "volume": ["40"], "fpage": ["2283"], "lpage": ["2294"]}, {"label": ["3"], "element-citation": ["\n"], "surname": ["Partensky", "Hess", "Vaulot"], "given-names": ["F", "WR", "D"], "year": ["1999"], "article-title": ["\n"], "italic": ["Prochlorococcus"], "source": ["Microbiol Molr Biol Rev"], "volume": ["63(1)"], "fpage": ["106"], "lpage": ["127"]}, {"label": ["19"], "element-citation": ["\n"], "surname": ["Storz", "Gottesman", "Gesteland", "Cech", "Atkins"], "given-names": ["G", "S", "RF", "TR", "JF"], "year": ["2006"], "article-title": ["Versatile roles of small RNA regulators in bacteria."], "source": ["The RNA World"], "publisher-loc": ["Cold Spring Harbor"], "publisher-name": ["Cold Spring Harbor Laboratory Press"], "fpage": ["567"], "lpage": ["594"]}, {"label": ["52"], "element-citation": ["\n"], "surname": ["Gerdes", "Jacobsen", "Franch"], "given-names": ["K", "JS", "T"], "year": ["1997"], "article-title": ["Plasmid stabilization by post-segregational killing."], "source": ["Genet Eng"], "volume": ["19"], "fpage": ["49"], "lpage": ["61"]}, {"label": ["61"], "element-citation": ["\n"], "surname": ["Moore", "Coe", "Zinser", "Saito", "Sullivan"], "given-names": ["LR", "A", "ER", "MA", "MB"], "year": ["2007"], "article-title": ["Culturing the marine cyanobacterium "], "italic": ["Prochlorococcus"], "source": ["Limnol Oceanogr Methods"], "volume": ["5"], "fpage": ["353"], "lpage": ["362"]}, {"label": ["64"], "element-citation": ["\n"], "surname": ["Hofacker", "Fontana", "Stadler", "Bonhoeffer", "Tacker"], "given-names": ["IL", "W", "PF", "M", "M"], "year": ["1994"], "article-title": ["Fast Folding and Comparison of RNA Secondary Structures."], "source": ["Monatshefte f Chemie"], "volume": ["125"], "fpage": ["167"], "lpage": ["188"]}]
PLoS Genet. 2008 Aug 29; 4(8):e1000173