Archaeal homologs of eukaryotic methylation guide small nucleolar RNAs: lessons from the Pyrococcus genomes

Ribose methylation is a prevalent type of nucleotide modification in rRNA. Eukaryotic rRNAs display a complex pattern of ribose methylations, amounting to 55 in yeast Saccharomyces cerevisiae and about 100 in vertebrates. Ribose methylations of eukaryotic rRNAs are each guided by a cognate small RNA, belonging to the family of box C/D antisense snoRNAs, through transient formation of a specific base-pairing at the rRNA modification site. In prokaryotes, the pattern of rRNA ribose methylations has been fully characterized in a single species so far, Escherichia coli, which contains only four ribose methylated rRNA nucleotides. However, the hyperthermophile archaeon Sulfolobus solfataricus contains, like eukaryotes, a large number of (yet unmapped) rRNA ribose methylations and homologs of eukaryotic box C/D small nucleolar ribonuclear proteins have been identified in archaeal genomes. We have therefore searched archaeal genomes for potential homologs of eukaryotic methylation guide small nucleolar RNAs, by combining searches for structured motifs with homology searches. We have identified a family of 46 small RNAs, conserved in the genomes of three hyperthermophile Pyrococcus species, which we have experimentally characterized in Pyrococcus abyssi. The Pyrococcus small RNAs, the first reported homologs of methylation guide small nucleolar RNAs in organisms devoid of a nucleus, appear as a paradigm of minimalist box C/D antisense RNAs. They differ from their eukaryotic homologs by their outstanding structural homogeneity, extended consensus box motifs and the quasi-systematic presence of two (instead of one) rRNA antisense elements. Remarkably, for each small RNA the two antisense elements always match rRNA sequences close to each other in rRNA structure, suggesting an important role in rRNA folding. Only a few of the predicted P. abyssi rRNA ribose methylations have been detected so far. Further analysis of these archaeal small RNAs could provide new insights into the origin and functions of methylation guide small nucleolar RNAs and illuminate the still elusive role of rRNA ribose methylations.     

A computational approach to identify genes for functional RNAs in genomic sequences

Currently there is no successful computational approach for identification of genes encoding novel functional RNAs (fRNAs) in genomic sequences. We have developed a machine learning approach using neural networks and support vector machines to extract common features among known RNAs for prediction of new RNA genes in the unannotated regions of prokaryotic and archaeal genomes. The Escherichia coli genome was used for development, but we have applied this method to several other bacterial and archaeal genomes. Networks based on nucleotide composition were 80–90% accurate in jackknife testing experiments for bacteria and 90–99% for hyperthermophilic archaea. We also achieved a significant improvement in accuracy by combining these predictions with those obtained using a second set of parameters consisting of known RNA sequence motifs and the calculated free energy of folding. Several known fRNAs not included in the training datasets were identified as well as several hundred predicted novel RNAs. These studies indicate that there are many unidentified RNAs in simple genomes that can be predicted computationally as a precursor to experimental study. Public access to our RNA gene predictions and an interface for user predictions is available via the web.     

sRNAbase: A Web Resource for Small RNAs in Enterobacteriaceae\m{1}

Bacterial small RNAs (sRNAs) have gained considerable attention due to their multivalent roles in the survival and pathogenesis of bacteria and mostly identified through bio-computational methods. A manually curated web-resource, sRNAbase has been constructed to give comprehensive and exhaustive information on non-coding small RNAs excluding tRNAs and rRNAs in Enterobacteriaceae family. The sRNA entries curated in sRNAbase contain experimentally verified small RNAs available in the literature and their partial/non-homologs reported within the related genomes from our earlier studies. The sRNAbase aims to facilitate the scientific community by providing information on the physical genomic location of the non-coding small RNAs, its alias names, sequences, strand orientation, gene identification numbers of the conserved genes that sandwiches the particular sRNA with possible functional role and a link to the PubMed literatures. Currently, sRNAbase holding information on 1986 entries belongs to 80 sRNA families spread over 45 Enterobacteriaceae genomes. The sRNAbase is accessible on the web at http://bicmku.in:8081/srnabase/.     

Comprehensive alanine-scanning mutagenesis of Escherichia coli CsrA defines two subdomains of critical functional importance

The RNA-binding protein CsrA (carbon storage regulator) of Escherichia coli is a global regulator of gene expression and is representative of the CsrA/RsmA family of bacterial proteins. These proteins act by regulating mRNA translation and stability and are antagonized by binding to small noncoding RNAs. Although the RNA target sequence and structure for CsrA binding have been well defined, little information exists concerning the protein requirements for RNA recognition. The three-dimensional structures of three CsrA/RsmA proteins were recently solved, revealing a novel protein fold consisting of two interdigitated monomers. Here, we performed comprehensive alanine-scanning mutagenesis on csrA of E. coli and tested the 58 resulting mutants for regulation of glycogen accumulation, motility, and biofilm formation. Quantitative effects of these mutations on expression of glgCA'-'lacZ, flhDC'-'lacZ, and pgaA'-'lacZ translational fusions were also examined, and eight of the mutant proteins were purified and tested for RNA binding. These studies identified two regions of the amino acid sequence that were critical for regulation and RNA binding, located within the first (beta1, residues 2-7) and containing the last (beta5, residues 40-47) beta-strands of CsrA. The beta1 and beta5 strands of opposite monomers lie adjacent and parallel to each other in the three-dimensional structure of this protein. Given the symmetry of the CsrA dimer, these findings imply that two distinct RNA binding surfaces or functional subdomains lie on opposite sides of the protein.     

细菌全局性调控因子CsrA的结构与功能

碳存储调控因子A (carbon storage regulator, CsrA) 是一种RNA结合蛋白,在细菌的碳代谢、生物被膜形成、运动性、病原菌毒力、群体感应、环二鸟苷酸信号合成、应激感应等多种生理过程中具有重要调节功能,是全局性调控蛋白.它通过与靶标mRNA的特异结合,抑制其翻译或增强其稳定性来调控下游基因的表达,属于转录后调控因子的范畴.CsrA蛋白的表达与活性受碳存储调控(Csr)系统本身多个自主调节回路的精密控制：　一些小的非编码RNA (snmRNAs,如CsrB/C）作为拮抗因子与CsrA二聚体结合并抑制其活性;而这些snmRNAs在体内又可在CsrD的辅助下被核糖核酸内切酶E和多核苷酸磷酸化酶降解,释放CsrA的活性.当前,对于Csr系统的调节作用、调控通路与机制的研究是细菌学研究的热点,本文综述了该蛋白及Csr系统的结构、功能和作用机制的最新研究进展.     

Global Assessment of Genomic Regions Required for Growth in Mycobacterium tuberculosis

Identifying genomic elements required for viability is central to our understanding of the basic physiology of bacterial pathogens. Recently, the combination of high-density mutagenesis and deep sequencing has allowed for the identification of required and conditionally required genes in many bacteria. Genes, however, make up only a part of the complex genomes of important bacterial pathogens. Here, we use an unbiased analysis to comprehensively identify genomic regions, including genes, domains, and intergenic elements, required for the optimal growth of Mycobacterium tuberculosis, a major global health pathogen. We found that several proteins jointly contain both domains required for optimal growth and domains that are dispensable. In addition, many non-coding regions, including regulatory elements and non-coding RNAs, are critical for mycobacterial growth. Our analysis shows that the genetic requirements for growth are more complex than can be appreciated using gene-centric analysis.     

Solution structure of the carbon storage regulator protein CsrA from Escherichia coli

The carbon storage regulator A (CsrA) is a protein responsible for the repression of a variety of stationary-phase genes in bacteria. In this work, we describe the nuclear magnetic resonance (NMR)-based structure of the CsrA dimer and its RNA-binding properties. CsrA is a dimer of two identical subunits, each composed of five strands, a small alpha-helix and a flexible C terminus. NMR titration experiments suggest that the beta1-beta2 and beta3-beta4 loops and the C-terminal helix are important elements in RNA binding. Even though the beta3-beta4 loop contains a highly conserved RNA-binding motif, GxxG, typical of KH domains, our structure excludes CsrA from being a member of this protein family, as previously suggested. A mechanism for the recognition of mRNAs downregulated by CsrA is proposed.     

Molecular basis for the wide range of affinity found in Csr/Rsm protein–RNA recognition

The carbon storage regulator/regulator of secondary metabolism (Csr/Rsm) type of small non-coding RNAs (sRNAs) is widespread throughout bacteria and acts by sequestering the global translation repressor protein CsrA/RsmE from the ribosome binding site of a subset of mRNAs. Although we have previously described the molecular basis of a high affinity RNA target bound to RsmE, it remains unknown how other lower affinity targets are recognized by the same protein. Here, we have determined the nuclear magnetic resonance solution structures of five separate GGA binding motifs of the sRNA RsmZ of Pseudomonas fluorescens in complex with RsmE. The structures explain how the variation of sequence and structural context of the GGA binding motifs modulate the binding affinity for RsmE by five orders of magnitude (∼10 nM to ∼3 mM, K_d). Furthermore, we see that conformational adaptation of protein side-chains and RNA enable recognition of different RNA sequences by the same protein contributing to binding affinity without conferring specificity. Overall, our findings illustrate how the variability in the Csr/Rsm protein–RNA recognition allows a fine-tuning of the competition between mRNAs and sRNAs for the CsrA/RsmE protein.     

Prediction of structured non-coding RNAs in the genomes of the nematodes Caenorhabditis elegans and Caenorhabditis briggsae

We present a survey for non-coding RNAs and other structured RNA motifs in the genomes of Caenorhabditis elegans and Caenorhabditis briggsae using the RNAz program. This approach explicitly evaluates comparative sequence information to detect stabilizing selection acting on RNA secondary structure. We detect 3,672 structured RNA motifs, of which only 678 are known non-translated RNAs (ncRNAs) or clear homologs of known C. elegans ncRNAs. Most of these signals are located in introns or at a distance from known protein-coding genes. With an estimated false positive rate of about 50% and a sensitivity on the order of 50%, we estimate that the nematode genomes contain between 3,000 and 4,000 RNAs with evolutionary conserved secondary structures. Only a small fraction of these belongs to the known RNA classes, including tRNAs, snoRNAs, snRNAs, or microRNAs. A relatively small class of ncRNA candidates is associated with previously observed RNA-specific upstream elements.