Biocomputational prediction of non-coding RNAs in model cyanobacteria

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

[1]  Ilka M. Axmann,et al.  Two distinct types of 6S RNA in Prochlorococcus. , 2007, Gene.

[2]  M. Hagemann,et al.  Construction of promoter probe vectors for Synechocystis sp. PCC 6803 using the light-emitting reporter systems Gfp and LuxAB. , 2000, Journal of microbiological methods.

[3]  S. Kurtz The Vmatch large scale sequence analysis software , 2003 .

[4]  Feng Chen,et al.  Patterns and Implications of Gene Gain and Loss in the Evolution of Prochlorococcus , 2007, PLoS genetics.

[5]  Ilka M. Axmann,et al.  An internal antisense RNA regulates expression of the photosynthesis gene isiA. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[6]  J. Vogel,et al.  Small non-coding RNAs and the bacterial outer membrane. , 2006, Current opinion in microbiology.

[7]  M. Gelfand,et al.  Comparative Genomics of Thiamin Biosynthesis in Procaryotes , 2002, The Journal of Biological Chemistry.

[8]  P. Babitzke,et al.  CsrB sRNA family: sequestration of RNA-binding regulatory proteins. , 2007, Current opinion in microbiology.

[9]  H. Margalit,et al.  A survey of small RNA-encoding genes in Escherichia coli. , 2003, Nucleic acids research.

[10]  David L. Wheeler,et al.  GenBank , 2015, Nucleic Acids Res..

[11]  A. Glatz,et al.  Chaperonin genes of the Synechocystis PCC 6803 are differentially regulated under light-dark transition during heat stress. , 1997, Biochemical and biophysical research communications.

[12]  J. Waterbury,et al.  Generic assignments, strain histories, and properties of pure cultures of cyanobacteria , 1979 .

[13]  P. Stadler,et al.  Secondary structure prediction for aligned RNA sequences. , 2002, Journal of molecular biology.

[14]  J. Vogel,et al.  The RNA chaperone Hfq is essential for the virulence of Salmonella typhimurium , 2007, Molecular microbiology.

[15]  N. Wingreen,et al.  The Small RNA Chaperone Hfq and Multiple Small RNAs Control Quorum Sensing in Vibrio harveyi and Vibrio cholerae , 2004, Cell.

[16]  M. Bes,et al.  Identification of a furA cis antisense RNA in the cyanobacterium Anabaena sp. PCC 7120. , 2006, Journal of molecular biology.

[17]  D. Patel,et al.  RNA bulges as architectural and recognition motifs. , 2000, Structure.

[18]  S. Salzberg,et al.  Rapid, accurate, computational discovery of Rho-independent transcription terminators illuminates their relationship to DNA uptake , 2007, Genome Biology.

[19]  Stan J. J. Brouns,et al.  Small CRISPR RNAs Guide Antiviral Defense in Prokaryotes , 2008, Science.

[20]  Jörg Vogel,et al.  How to find small non-coding RNAs in bacteria , 2005, Biological chemistry.

[21]  J. Vogel,et al.  The cyanobacterial homologue of the RNA chaperone Hfq is essential for motility of Synechocystis sp. PCC 6803. , 2008, Microbiology.

[22]  E. Rivas,et al.  Identification of differentially expressed small non-coding RNAs in the legume endosymbiont Sinorhizobium meliloti by comparative genomics , 2007, Molecular microbiology.

[23]  Sean R. Eddy,et al.  Rfam: annotating non-coding RNAs in complete genomes , 2004, Nucleic Acids Res..

[24]  Hitoshi Nakamoto,et al.  Targeted inactivation of the hrcA repressor gene in cyanobacteria , 2003, FEBS letters.

[25]  Andrew C. Tolonen,et al.  Genetic Manipulation of Prochlorococcus Strain MIT9313: Green Fluorescent Protein Expression from an RSF1010 Plasmid and Tn5 Transposition , 2006, Applied and Environmental Microbiology.

[26]  Jeffrey E. Barrick,et al.  New RNA motifs suggest an expanded scope for riboswitches in bacterial genetic control. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Björn Voß,et al.  Structural analysis of aligned RNAs , 2006, Nucleic acids research.

[28]  Younghoon Lee,et al.  Regulation of 6S RNA biogenesis by switching utilization of both sigma factors and endoribonucleases. , 2004, Nucleic acids research.

[29]  J. Thompson,et al.  CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. , 1994, Nucleic acids research.

[30]  Ilka M. Axmann,et al.  A motif-based search in bacterial genomes identifies the ortholog of the small RNA Yfr1 in all lineages of cyanobacteria , 2007, BMC Genomics.

[31]  Pascale Cossart,et al.  Identification of new noncoding RNAs in Listeria monocytogenes and prediction of mRNA targets , 2007, Nucleic acids research.

[32]  Sam Griffiths-Jones,et al.  RALEE--RNA ALignment Editor in Emacs , 2005, Bioinform..

[33]  Hanah Margalit,et al.  Small RNAs encoded within genetic islands of Salmonella typhimurium show host-induced expression and role in virulence , 2008, Nucleic acids research.

[34]  Peter F Stadler,et al.  Fast and reliable prediction of noncoding RNAs , 2005, Proc. Natl. Acad. Sci. USA.

[35]  W. L. Ruzzo,et al.  6S RNA is a widespread regulator of eubacterial RNA polymerase that resembles an open promoter. , 2005, RNA.

[36]  Ilka M. Axmann,et al.  Identification of cyanobacterial non-coding RNAs by comparative genome analysis , 2005, Genome Biology.

[37]  P. Wincker,et al.  Unraveling the genomic mosaic of a ubiquitous genus of marine cyanobacteria , 2008, Genome Biology.

[38]  S. Tabata,et al.  Complete Genomic Structure of the Bloom-forming Toxic Cyanobacterium Microcystis aeruginosa NIES-843 , 2008, DNA research : an international journal for rapid publication of reports on genes and genomes.

[39]  L. Lindahl,et al.  Diverse mechanisms for regulating ribosomal protein synthesis in Escherichia coli. , 1994, Progress in nucleic acid research and molecular biology.

[40]  R. Simons,et al.  Antisense RNA control in bacteria, phages, and plasmids. , 1994, Annual review of microbiology.

[41]  S. Heeb,et al.  Genome-wide search reveals a novel GacA-regulated small RNA in Pseudomonas species , 2008, BMC Genomics.

[42]  M. Sugita,et al.  A cyanobacterial non-coding RNA, Yfr1, is required for growth under multiple stress conditions. , 2007, Plant & cell physiology.

[43]  S. Gottesman Micros for microbes: non-coding regulatory RNAs in bacteria. , 2005, Trends in genetics : TIG.

[44]  E. Wagner,et al.  A Repeated GGA Motif Is Critical for the Activity and Stability of the Riboregulator RsmY of Pseudomonas fluorescens* , 2004, Journal of Biological Chemistry.

[45]  Z. Wang,et al.  Escherichia coli 6S RNA gene is part of a dual-function transcription unit , 1985, Journal of bacteriology.

[46]  F. Barloy-Hubler,et al.  Identification of chromosomal alpha-proteobacterial small RNAs by comparative genome analysis and detection in Sinorhizobium meliloti strain 1021 , 2007, BMC Genomics.