Most RNAs regulating ribosomal protein biosynthesis in Escherichia coli are narrowly distributed to Gammaproteobacteria

MnoSR removal in Mycobacterium smegmatis triggers broad transcriptional response to 1,3-propanediol and glucose as sole carbon sources

Introduction

The two-component signal transduction systems play an essential role in the adaptation of bacteria to changing environmental conditions. One of them is the MnoSR system involved in the regulation of methylotrophic metabolism in M. smegmatis.

Methods

Mycobacterium smegmatis mutant strains ΔmnoS, ΔmnoR and ΔmnoS/R lacking functional mnoS, mnoR and both genes were generated using a homologous recombination approach. MnoR recombinant protein was purified by affinity column chromatography. The present study employs molecular biology techniques: cloning strategies, global RNA sequencing, qRT-PCR, EMSA, Microscale thermophoresis, and bioinformatics analysis.

Results and discussion

The ∆mnoS, ∆mnoR, and ∆mnoS/R mutant strains were generated and cultured in the presence of defined carbon sources. Growth curve analysis confirmed that inactivation of the MnoSR impairs the ability of M. smegmatis cells to use alcohols such as 1,3-propanediol and ethanol but improves the bacterial growth on ethylene glycol, xylitol, and glycerol. The total RNA sequencing method was employed to understand the importance of MnoSR in the global responses of mycobacteria to limited carbon access and in carbon-rich conditions. The loss of MnoSR significantly affected carbon utilization in the case of mycobacteria cultured on glucose or 1,3-propanediol as sole carbon sources as it influenced the expression of multiple metabolic pathways. The numerous transcriptional changes could not be linked to the presence of evident MnoR DNA-binding sites within the promotor regions for the genes outside of the mno operon. This was confirmed by EMSA and microscale thermophoresis with mutated MnoR binding consensus region. Our comprehensive analysis highlights the system's vital role in metabolic adaptability, providing insights into its potential impact on the environmental survival of mycobacteria.

Lemon essential oil nanoemulsions: Potential natural inhibitors against Escherichia coli

Lemon essential oil (LEO) is a common natural antibacterial substance, and encapsulating LEO into nanoemulsions (NEs) can improve their stability and broaden its application. Our study aimed to investigate the bacterial inhibitory effect of LEO-NEs against Escherichia coli (E. coli). Results showed that the minimum inhibitory concentration (MIC) of LEO-NEs was 6.25 mg/mL, and the time-kill curve showed that E. coli were significantly killed by LEO-NEs after 5 h of treatment at 1MIC. Flow-cytometry analysis showed that LEO-NEs adversely affected the cell-membrane depolarisation, cell-membrane integrity, and efflux pump function of E. coli. Confocal laser scanning microscopy demonstrated that 8MIC of LEO-NEs induced changes in the cell-membrane permeability and cell-wall integrity of E. coli. Proteomic results suggested that the mode of action LEO-NEs against E. coli was to enhance bacterial chemotaxis and significantly inhibit ribosomal assembly. They may also affect butyric acid, ascorbic acid and aldehyde metabolism, and sulphur-relay system pathways. In conclusion, LEO-NEs had potential application as a natural antibacterial agent for the control of E. coli in the food industry.

Extraribosomal Functions of Bacterial Ribosomal Proteins-An Update, 2023

Ribosomal proteins (r-proteins) are abundant, highly conserved, and multifaceted cellular proteins in all domains of life. Most r-proteins have RNA-binding properties and can form protein-protein contacts. Bacterial r-proteins govern the co-transcriptional rRNA folding during ribosome assembly and participate in the formation of the ribosome functional sites, such as the mRNA-binding site, tRNA-binding sites, the peptidyl transferase center, and the protein exit tunnel. In addition to their primary role in a cell as integral components of the protein synthesis machinery, many r-proteins can function beyond the ribosome (the phenomenon known as moonlighting), acting either as individual regulatory proteins or in complexes with various cellular components. The extraribosomal activities of r-proteins have been studied over the decades. In the past decade, our understanding of r-protein functions has advanced significantly due to intensive studies on ribosomes and gene expression mechanisms not only in model bacteria like Escherichia coli or Bacillus subtilis but also in little-explored bacterial species from various phyla. The aim of this review is to update information on the multiple functions of r-proteins in bacteria.

Single-Molecule Correlated Chemical Probing: A Revolution in RNA Structure Analysis

RNA molecules convey biological information both in their linear sequence and in their base-paired secondary and tertiary structures. Chemical probing experiments, which involve treating an RNA with a reagent that modifies conformationally dynamic nucleotides, have broadly enabled examination of short- and long-range RNA structure in diverse contexts, including in living cells. For decades, chemical probing experiments have been interpreted in a per-nucleotide way, such that the reactivity measured at each nucleotide reports the average structure at a position over all RNA molecules within a sample. However, there are numerous important cases where per-nucleotide chemical probing falls short, including for RNAs that are bound by proteins, RNAs that form complex higher order structures, and RNAs that sample multiple conformations.Recent experimental and computational innovations have started a revolution in RNA structure analysis by transforming chemical probing into a massively parallel, single-molecule experiment. Enabled by a specialized reverse transcription strategy called mutational profiling (MaP), multiple chemical modification events can be measured within individual RNA molecules. Nucleotides that communicate structurally through direct base pairing or large-scale folding-unfolding transitions will react with chemical probes in a correlated manner, thereby revealing structural complexity hidden to conventional approaches. These single-molecule correlated chemical probing (smCCP) experiments can be interpreted to directly identify nucleotides that base pair (the PAIR-MaP strategy) and to reveal long-range, through-space structural communication (RING-MaP). Correlated probing can also define the thermodynamic populations of complex RNA ensembles (DANCE-MaP). Complex RNA-protein networks can be interrogated by cross-linking proteins to RNA and measuring correlations between cross-linked positions (RNP-MaP).smCCP thus visualizes RNA secondary and higher-order structure with unprecedented accuracy, defining novel structures, RNA-protein interaction networks, time-resolved dynamics, and allosteric structural switches. These strategies are not mutually exclusive; in favorable cases, multiple levels of RNA structure ─ base pairing, through-space structural communication, and equilibrium ensembles ─ can be resolved concurrently. The physical experimentation required for smCCP is profoundly simple, and experiments are readily performed in cells on RNAs of any size, including large noncoding RNAs and mRNAs. Single-molecule correlated chemical probing is paving the way for a new generation of biophysical studies on RNA in living systems.

At the Crossroad of Nucleotide Dynamics and Protein Synthesis in Bacteria

Nucleotides are at the heart of the most essential biological processes in the cell, be it as key protagonists in the dogma of molecular biology or by regulating multiple metabolic pathways. The dynamic nature of nucleotides, the cross talk between them, and their constant feedback to and from the cell's metabolic state position them as a hallmark of adaption toward environmental and growth challenges. It has become increasingly clear how the activity of RNA polymerase, the synthesis and maintenance of tRNAs, mRNA translation at all stages, and the biogenesis and assembly of ribosomes are fine-tuned by the pools of intracellular nucleotides. With all aspects composing protein synthesis involved, the ribosome emerges as the molecular hub in which many of these nucleotides encounter each other and regulate the state of the cell. In this review, we aim to highlight intracellular nucleotides in bacteria as dynamic characters permanently cross talking with each other and ultimately regulating protein synthesis at various stages in which the ribosome is mainly the principal character.

Unexpected absence of ribosomal protein genes from metagenome-assembled genomes

Metagenome-assembled genomes (MAGs) have revealed the hidden diversity and functions of uncultivated microbes, but their reconstruction from metagenomes remains a computationally difficult task. Repetitive or exogenous sequences, such as ribosomal RNA and horizontally transferred genes, are frequently absent from MAGs because of misassembly and binning errors. Here, we report that ribosomal protein genes are also often absent from MAGs, although they are neither repetitive nor exogenous. Comprehensive analyses of more than 190,000 MAGs revealed that these genes could be missing in more than 20-40% of near-complete (i.e., with completeness of 90% or higher) MAGs. While some uncultivated environmental microbes intrinsically lack some ribosomal protein genes, we found that this unexpected absence is largely due to special evolutionary patterns of codon usage bias in ribosomal protein genes and algorithmic characteristics of metagenomic binning, which is dependent on tetranucleotide frequencies of contigs. This problem reflects the microbial life-history strategy. Fast-growing microbes tend to have this difficulty, likely because of strong evolutionary pressures on ribosomal protein genes toward the efficient assembly of ribosomes. Our observations caution those who study genomics and phylogeny of uncultivated microbes, the diversity and evolution of microbial genes in the central dogma, and bioinformatics in metagenomics.

Differences in gene expression patterns between cultured and natural Haloquadratum walsbyi ecotypes

Solar crystallizer ponds are characterized by high population density with a relatively simple community structure in terms of species composition. The microbial community in the solar saltern of Santa Pola (Alicante, Spain), is largely dominated by the hyperhalophilic square archaeon Haloquadratum walsbyi. Here we studied metatranscriptomes retrieved from a crystallizer pond during the winter of 2012 and summer of 2014 and compared Hqr. walsbyi's transcription patterns with that of the cultured strain Hqr. walsbyi HBSQ001. Significant differences were found between natural and the cultured grown strain in the distribution of transcript levels per gene. This likely reflects the adaptation of the cultured strain to the relative homogeneous growth conditions while the natural species, which is represented by multiple ecotypes, is adapted to heterogeneous environmental conditions and challenges of nutrient competition, viral attack, and other stressors. An important consequence of this study is that expression patterns obtained under artificial cultivation conditions cannot be directly extrapolated to gene expression under natural conditions. Moreover, we found 195 significantly differential expressed genes between the seasons, with 140 genes being higher expressed in winter and mainly encode proteins involved in energy and carbon source acquiring processes, and in stress responses.

Discovery of a large-scale, cell-state-responsive allosteric switch in the 7SK RNA using DANCE-MaP

7SK is a conserved noncoding RNA that regulates transcription by sequestering the transcription factor P-TEFb. 7SK function entails complex changes in RNA structure, but characterizing RNA dynamics in cells remains an unsolved challenge. We developed a single-molecule chemical probing strategy, DANCE-MaP (deconvolution and annotation of ribonucleic conformational ensembles), that defines per-nucleotide reactivity, direct base pairing interactions, tertiary interactions, and thermodynamic populations for each state in RNA structural ensembles from a single experiment. DANCE-MaP reveals that 7SK RNA encodes a large-scale structural switch that couples dissolution of the P-TEFb binding site to structural remodeling at distal release factor binding sites. The 7SK structural equilibrium shifts in response to cell growth and stress and can be targeted to modulate expression of P-TEFbresponsive genes. Our study reveals that RNA structural dynamics underlie 7SK function as an integrator of diverse cellular signals to control transcription and establishes the power of DANCE-MaP to define RNA dynamics in cells.

Interplay between non-coding RNA transcription, stringent phenotype and antibiotic production in Streptomyces

While in recent years the key role of non-coding RNAs (ncRNAs) in regulation of gene expression has become increasingly evident, their interaction with the global regulatory circuits is still obscure. Here we analyzed the structure and organization of the transcriptome of Streptomyces ambofaciens, the producer of spiramycin. We identified ncRNAs including 45 small-RNAs (sRNAs) and 119 antisense-RNAs (asRNAs I) that appear transcribed from dedicated promoters. Some sRNAs and asRNAs are unprecedented in Streptomyces, and were predicted to target mRNAs encoding proteins involved in transcription, translation, ribosomal structure and biogenesis, and regulation of morphological and biochemical differentiation. We then compared ncRNA expression in three strains: i.) the wild type strain; ii.) an isogenic pirA-defective mutant with central carbon metabolism imbalance, "relaxed" phenotype, and repression of antibiotic production; iii.) a pirA-derivative strain harboring a "stringent" RNA polymerase that suppresses pirA-associated phenotypes. Data indicated that expression of most ncRNAs was correlated to the stringent/relaxed phenotype suggesting novel effector mechanisms of the stringent response.

Yeast Crf1p: An activator in need is an activator indeed

Ribosome biogenesis is an energetically costly process, and tight regulation is required for stoichiometric balance between components. This requires coordination of RNA polymerases I, II, and III. Lack of nutrients or the presence of stress leads to downregulation of ribosome biogenesis, a process for which mechanistic target of rapamycin complex I (mTORC1) is key. mTORC1 activity is communicated by means of specific transcription factors, and in yeast, which is a primary model system in which transcriptional coordination has been delineated, transcription factors involved in regulation of ribosomal protein genes include Fhl1p and its cofactors, Ifh1p and Crf1p. Ifh1p is an activator, whereas Crf1p has been implicated in maintaining the repressed state upon mTORC1 inhibition. Computational analyses of evolutionary relationships have indicated that Ifh1p and Crf1p descend from a common ancestor. Here, we discuss recent evidence, which suggests that Crf1p also functions as an activator. We propose a model that consolidates available experimental evidence, which posits that Crf1p functions as an alternate activator to prevent the stronger activator Ifh1p from re-binding gene promoters upon mTORC1 inhibition. The correlation between retention of Crf1p in related yeast strains and duplication of ribosomal protein genes suggests that this backup activation may be important to ensure gene expression when Ifh1p is limiting. With ribosome biogenesis as a hallmark of cell growth, failure to control assembly of ribosomal components leads to several human pathologies. A comprehensive understanding of mechanisms underlying this process is therefore of the essence.

Proteomic profiling of the endogenous peptides of MRSA and MSSA

Staphylococcus aureus is a Gram-positive bacterium that can cause diverse skin and soft tissue infections. Methicillin-resistant Staphylococcus aureus (MRSA) can cause more severe infections than methicillin-susceptible Staphylococcus aureus (MSSA). Nevertheless, the physiological and metabolic regulation of MSSA and MRSA has not been well studied. In light of the increased interest in endogenous peptides and recognition of the important roles that they play, we studied the endogenous peptidome of MSSA and MRSA. We identified 1,065 endogenous peptides, among which 435 were differentially expressed (DE), with 292 MSSA-abundant endogenous peptides and 35 MRSA-abundant endogenous peptides. MSSA-abundant endogenous peptides have significantly enriched “VXXXK” motif of at the C-terminus. MSSA-abundant endogenous peptides are involved in penicillin-binding and immune responses, whereas MRSA-abundant endogenous peptides are associated with antibiotic resistance and increased toxicity. Our characterization of the peptidome of MSSA and MRSA provides a rich resource for future studies to explore the functional regulation of drug resistance in S. aureus and may also help elucidate the mechanisms of its pathogenicity and the development of treatments.

Regulation of Ribosomal Protein Synthesis in Mycobacteria: The Autogenous Control of rpsO

The autogenous regulation of ribosomal protein (r-protein) synthesis plays a key role in maintaining the stoichiometry of ribosomal components in bacteria. In this work, taking the rpsO gene as a classic example, we addressed for the first time the in vivo regulation of r-protein synthesis in the mycobacteria M. smegmatis (Msm) and M. tuberculosis (Mtb). We used a strategy based on chromosomally integrated reporters under the control of the rpsO regulatory regions and the ectopic expression of Msm S15 to measure its impact on the reporter expression. Because the use of E. coli as a host appeared inefficient, a fluorescent reporter system was developed by inserting Msm or Mtb rpsO-egfp fusions into the Msm chromosome and expressing Msm S15 or E. coli S15 in trans from a novel replicative shuttle vector, pAMYC. The results of the eGFP expression measurements in Msm cells provided evidence that the rpsO gene in Msm and Mtb was feedback-regulated at the translation level. The mutagenic analysis showed that the folding of Msm rpsO 5'UTR in a pseudoknot appeared crucial for repression by both Msm S15 and E. coli S15, thus indicating a striking resemblance of the rpsO feedback control in mycobacteria and in E. coli.

Comparative genomics provides structural and functional insights into Bacteroides RNA biology

Bacteria employ noncoding RNA molecules for a wide range of biological processes, including scaffolding large molecular complexes, catalyzing chemical reactions, defending against phages, and controlling gene expression. Secondary structures, binding partners, and molecular mechanisms have been determined for numerous small noncoding RNAs (sRNAs) in model aerobic bacteria. However, technical hurdles have largely prevented analogous analyses in the anaerobic gut microbiota. While experimental techniques are being developed to investigate the sRNAs of gut commensals, computational tools and comparative genomics can provide immediate functional insight. Here, using Bacteroides thetaiotaomicron as a representative microbiota member, we illustrate how comparative genomics improves our understanding of the RNA biology in an understudied gut bacterium. We investigate putative RNA-binding proteins and predict a Bacteroides cold-shock protein homologue to have an RNA-related function. We apply an in-silico protocol incorporating both sequence and structural analysis to determine the consensus structures and conservation of nine Bacteroides noncoding RNA families. Using structure probing, we validate and refine these predictions, and deposit them in the Rfam database. Through synteny analyses, we illustrate how genomic co-conservation can serve as a predictor of sRNA function. Altogether, this work showcases the power of RNA informatics for investigating the RNA biology of anaerobic microbiota members.

Interplay between Non-Coding RNA Transcription, Stringent/Relaxed Phenotype and Antibiotic Production in Streptomyces ambofaciens

While in recent years the key role of non-coding RNAs (ncRNAs) in the regulation of gene expression has become increasingly evident, their interaction with the global regulatory circuits is still obscure. Here we analyzed the structure and organization of the transcriptome of Streptomyces ambofaciens, the producer of spiramycin. We identified ncRNAs including 45 small-RNAs (sRNAs) and 119 antisense-RNAs (asRNAs I) that appear transcribed from dedicated promoters. Some sRNAs and asRNAs are unprecedented in Streptomyces and were predicted to target mRNAs encoding proteins involved in transcription, translation, ribosomal structure and biogenesis, and regulation of morphological and biochemical differentiation. We then compared ncRNA expression in three strains: (i) the wild-type strain; (ii) an isogenic pirA-defective mutant with central carbon metabolism imbalance, “relaxed” phenotype, and repression of antibiotic production; and (iii) a pirA-derivative strain harboring a “stringent” RNA polymerase that suppresses pirA-associated phenotypes. Data indicated that the expression of most ncRNAs was correlated to the stringent/relaxed phenotype suggesting novel effector mechanisms of the stringent response.

Quantitative mapping of mRNA 3' ends in Pseudomonas aeruginosa reveals a pervasive role for premature 3' end formation in response to azithromycin

Pseudomonas aeruginosa produces serious chronic infections in hospitalized patients and immunocompromised individuals, including patients with cystic fibrosis. The molecular mechanisms by which P. aeruginosa responds to antibiotics and other stresses to promote persistent infections may provide new avenues for therapeutic intervention. Azithromycin (AZM), an antibiotic frequently used in cystic fibrosis treatment, is thought to improve clinical outcomes through a number of mechanisms including impaired biofilm growth and quorum sensing (QS). The mechanisms underlying the transcriptional response to AZM remain unclear. Here, we interrogated the P. aeruginosa transcriptional response to AZM using a fast, cost-effective genome-wide approach to quantitate RNA 3’ ends (3pMap). We also identified hundreds of P. aeruginosa genes with high incidence of premature 3’ end formation indicative of riboregulation in their transcript leaders using 3pMap. AZM treatment of planktonic and biofilm cultures alters the expression of hundreds of genes, including those involved in QS, biofilm formation, and virulence. Strikingly, most genes downregulated by AZM in biofilms had increased levels of intragenic 3’ ends indicating premature transcription termination, transcriptional pausing, or accumulation of stable intermediates resulting from the action of nucleases. Reciprocally, AZM reduced premature intragenic 3’ end termini in many upregulated genes. Most notably, reduced termination accompanied robust induction of obgE, a GTPase involved in persister formation in P. aeruginosa. Our results support a model in which AZM-induced changes in 3’ end formation alter the expression of central regulators which in turn impairs the expression of QS, biofilm formation and stress response genes, while upregulating genes associated with persistence.

Pseudomonas aeruginosa is a common source of hospital-acquired infections and causes prolonged illness in patients with cystic fibrosis. P. aeruginosa infections are often treated with the macrolide antibiotic azithromycin, which changes the expression of many genes involved in infection. By examining such expression changes at nucleotide resolution, we found azithromycin treatment alters the locations of mRNA 3’ ends suggesting most downregulated genes are subject to premature 3’ end formation. We further identified candidate RNA regulatory elements that P. aeruginosa may use to control gene expression. Our work provides new insights in P. aeruginosa gene regulation and its response to antibiotics.

Mfd regulates RNA polymerase association with hard-to-transcribe regions in vivo, especially those with structured RNAs

RNA polymerase (RNAP) encounters various roadblocks during transcription. These obstacles can impede RNAP movement and influence transcription, ultimately necessitating the activity of RNAP-associated factors. One such factor is the bacterial protein Mfd, a highly conserved DNA translocase and evolvability factor that interacts with RNAP. Although Mfd is thought to function primarily in the repair of DNA lesions that stall RNAP, increasing evidence suggests that it may also be important for transcription regulation. However, this is yet to be fully characterized. To shed light on Mfd's in vivo functions, we identified the chromosomal regions where it associates. We analyzed Mfd's impact on RNAP association and transcription regulation genome-wide. We found that Mfd represses RNAP association at many chromosomal regions. We found that these regions show increased RNAP pausing, suggesting that they are hard to transcribe. Interestingly, we noticed that the majority of the regions where Mfd regulates transcription contain highly structured regulatory RNAs. The RNAs identified regulate a myriad of biological processes, ranging from metabolism to transfer RNA regulation to toxin-antitoxin (TA) functions. We found that cells lacking Mfd are highly sensitive to toxin overexpression. Finally, we found that Mfd promotes mutagenesis in at least one toxin gene, suggesting that its function in regulating transcription may promote evolution of certain TA systems and other regions containing strong RNA secondary structures. We conclude that Mfd is an RNAP cofactor that is important, and at times critical, for transcription regulation at hard-to-transcribe regions, especially those that express structured regulatory RNAs.

Characterization of Regulatory Elements of L11 and L1 Operons in Thermophilic Bacteria and Archaea

Ribosomal protein L1 is a conserved two-domain protein that is involved in formation of the L1 stalk of the large ribosomal subunit. When there are no free binding sites available on the ribosomal 23S RNA, the protein binds to the specific site on the mRNA of its own operon (L11 operon in bacteria and L1 operon in archaea) preventing translation. Here we show that the regulatory properties of the r-protein L1 and its domain I are conserved in the thermophilic bacteria Thermus and Thermotoga and in the halophilic archaeon Haloarcula marismortui. At the same time the revealed features of the operon regulation in thermophilic bacteria suggest presence of two regulatory regions.

Comparative genomics identifies thousands of candidate structured RNAs in human microbiomes

BACKGROUND

Structured RNAs play varied bioregulatory roles within microbes. To date, hundreds of candidate structured RNAs have been predicted using informatic approaches that search for motif structures in genomic sequence data. The human microbiome contains thousands of species and strains of microbes. Yet, much of the metagenomic data from the human microbiome remains unmined for structured RNA motifs primarily due to computational limitations.

RESULTS

We sought to apply a large-scale, comparative genomics approach to these organisms to identify candidate structured RNAs. With a carefully constructed, though computationally intensive automated analysis, we identify 3161 conserved candidate structured RNAs in intergenic regions, as well as 2022 additional candidate structured RNAs that may overlap coding regions. We validate the RNA expression of 177 of these candidate structures by analyzing small fragment RNA-seq data from four human fecal samples.

CONCLUSIONS

This approach identifies a wide variety of candidate structured RNAs, including tmRNAs, antitoxins, and likely ribosome protein leaders, from a wide variety of taxa. Overall, our pipeline enables conservative predictions of thousands of novel candidate structured RNAs from human microbiomes.

Autoregulation of yeast ribosomal proteins discovered by efficient search for feedback regulation

Post-transcriptional autoregulation of gene expression is common in bacteria but many fewer examples are known in eukaryotes. We used the yeast collection of genes fused to GFP as a rapid screen for examples of feedback regulation in ribosomal proteins by overexpressing a non-regulatable version of a gene and observing the effects on the expression of the GFP-fused version. We tested 95 ribosomal protein genes and found a wide continuum of effects, with 30% showing at least a 3-fold reduction in expression. Two genes, RPS22B and RPL1B, showed over a 10-fold repression. In both cases the cis-regulatory segment resides in the 5’ UTR of the gene as shown by placing that segment of the mRNA upstream of GFP alone and demonstrating it is sufficient to cause repression of GFP when the protein is over-expressed. Further analyses showed that the intron in the 5’ UTR of RPS22B is required for regulation, presumably because the protein inhibits splicing that is necessary for translation. The 5’ UTR of RPL1B contains a sequence and structure motif that is conserved in the binding sites of Rpl1 orthologs from bacteria to mammals, and mutations within the motif eliminate repression.

Here, the authors screen for feedback regulation of ribosomal proteins by overexpressing a non- regulatable version of a gene and observing its effects on the expression of the GFP-fused version. They find that 30% show at least a 3-fold reduction in expression and two genes show a 10-fold reduction with the regulatory site being in the 5’ untranslated region of the gene.

Single-Nucleotide RNA Maps for the Two Major Nosocomial Pathogens Enterococcus faecalis and Enterococcus faecium

Enterococcus faecalis and faecium are two major representative clinical strains of the Enterococcus genus and are sadly notorious to be part of the top agents responsible for nosocomial infections. Despite their critical implication in worldwide public healthcare, essential and available resources such as deep transcriptome annotations remain poor, which also limits our understanding of post-transcriptional control small regulatory RNA (sRNA) functions in these bacteria. Here, using the dRNA-seq technique in combination with ANNOgesic analysis, we successfully mapped and annotated transcription start sites (TSS) of both E. faecalis V583 and E. faecium AUS0004 at single nucleotide resolution. Analyzing bacteria in late exponential phase, we capture ~40% (E. faecalis) and 43% (E. faecium) of the annotated protein-coding genes, determine 5′ and 3′ UTR (untranslated region) length, and detect instances of leaderless mRNAs. The transcriptome maps revealed sRNA candidates in both bacteria, some found in previous studies and new ones. Expression of candidate sRNAs is being confirmed under biologically relevant environmental conditions. This comprehensive global TSS mapping atlas provides a valuable resource for RNA biology and gene expression analysis in the Enterococci. It can be accessed online at www.helmholtz-hiri.de/en/datasets/enterococcus through an instance of the genomic viewer JBrowse.

Regulation of ribosomal protein genes: An ordered anarchy

Ribosomal protein genes are among the most highly expressed genes in most cell types. Their products are generally essential for ribosome synthesis, which is the cornerstone for cell growth and proliferation. Many cellular resources are dedicated to producing ribosomal proteins and thus this process needs to be regulated in ways that carefully balance the supply of nascent ribosomal proteins with the demand for new ribosomes. Ribosomal protein genes have classically been viewed as a uniform interconnected regulon regulated in eukaryotic cells by target of rapamycin and protein kinase A pathway in response to changes in growth conditions and/or cellular status. However, recent literature depicts a more complex picture in which the amount of ribosomal proteins produced varies between genes in response to two overlapping regulatory circuits. The first includes the classical general ribosome‐producing program and the second is a gene‐specific feature responsible for fine‐tuning the amount of ribosomal proteins produced from each individual ribosomal gene. Unlike the general pathway that is mainly controlled at the level of transcription and translation, this specific regulation of ribosomal protein genes is largely achieved through changes in pre‐mRNA splicing efficiency and mRNA stability. By combining general and specific regulation, the cell can coordinate ribosome production, while allowing functional specialization and diversity. Here we review the many ways ribosomal protein genes are regulated, with special focus on the emerging role of posttranscriptional regulatory events in fine‐tuning the expression of ribosomal protein genes and its role in controlling the potential variation in ribosome functions.

This article is categorized under:

Translation > Ribosome Biogenesis

Translation > Ribosome Structure/Function

Translation > Translation Regulation

Autogenous regulation in vivo of the rpmE gene encoding ribosomal protein L31 (bL31), a key component of the protein-protein intersubunit bridge B1b

Bacterial ribosomal proteins (r-proteins) encoded by nonessential genes often carry out very important tasks in translation. In particular, this is the case of a small basic bacteria-specific r-protein L31 (bL31). Recent studies revealed a crucial role of bL31 in formation of the protein-protein intersubunit bridge B1b and hence its contribution to ribosome dynamics. Our goal was to study in vivo regulation of the rpmE operon encoding bL31. We used a previously developed approach based on chromosomally integrated fusions with the lacZ reporter. E. coli rpmE is transcribed from two promoter regions, and translation of both mRNA transcripts was shown to be feedback regulated by bL31, indicating that the autogenous operator is located within the shorter transcript. The bL31-mediated control of rpmE is gene-specific, as no regulation was found for rpmE-unrelated reporters. Thus, bL31, as many other r-proteins, possesses dual activity in living cells, acting both as an integral ribosome component and an autogenous repressor. Phylogenetic studies revealed the presence of a highly conserved stem-loop structure in the rpmE 5'UTR, a presumable translational operator targeted by bL31, which was further confirmed by site-directed mutagenesis. This stable operator stem-loop separates an AU-rich translational enhancer from a Shine-Dalgarno element, which is a rare case of a noncontiguous translation initiation region. Sequence/structure computational approaches classify bL31 as an RNA-binding protein, consistent with its repressor function discovered here. Mutational analysis of bL31 showed that its unstructured amino-terminal part enriched in lysine is necessary for the repressor activity.

Discovery of 20 novel ribosomal leader candidates in bacteria and archaea

BACKGROUND

RNAs perform many functions in addition to supplying coding templates, such as binding proteins. RNA-protein interactions are important in multiple processes in all domains of life, and the discovery of additional protein-binding RNAs expands the scope for studying such interactions. To find such RNAs, we exploited a form of ribosomal regulation. Ribosome biosynthesis must be tightly regulated to ensure that concentrations of rRNAs and ribosomal proteins (r-proteins) match. One regulatory mechanism is a ribosomal leader (r-leader), which is a domain in the 5' UTR of an mRNA whose genes encode r-proteins. When the concentration of one of these r-proteins is high, the protein binds the r-leader in its own mRNA, reducing gene expression and thus protein concentrations. To date, 35 types of r-leaders have been validated or predicted.

RESULTS

By analyzing additional conserved RNA structures on a multi-genome scale, we identified 20 novel r-leader structures. Surprisingly, these included new r-leaders in the highly studied organisms Escherichia coli and Bacillus subtilis. Our results reveal several cases where multiple unrelated RNA structures likely bind the same r-protein ligand, and uncover previously unknown r-protein ligands. Each r-leader consistently occurs upstream of r-protein genes, suggesting a regulatory function. That the predicted r-leaders function as RNAs is supported by evolutionary correlations in the nucleotide sequences that are characteristic of a conserved RNA secondary structure. The r-leader predictions are also consistent with the locations of experimentally determined transcription start sites.

CONCLUSIONS

This work increases the number of known or predicted r-leader structures by more than 50%, providing additional opportunities to study structural and evolutionary aspects of RNA-protein interactions. These results provide a starting point for detailed experimental studies.

Recent advances in understanding mitochondrial genome diversity

Ever since its discovery, the double-stranded DNA contained in the mitochondria of eukaryotes has fascinated researchers because of its bacterial endosymbiotic origin, crucial role in encoding subunits of the respiratory complexes, compact nature, and specific inheritance mechanisms. In the last few years, high-throughput sequencing techniques have accelerated the sequencing of mitochondrial genomes (mitogenomes) and uncovered the great diversity of organizations, gene contents, and modes of replication and transcription found in living eukaryotes. Some early divergent lineages of unicellular eukaryotes retain certain synteny and gene content resembling those observed in the genomes of alphaproteobacteria (the inferred closest living group of mitochondria), whereas others adapted to anaerobic environments have drastically reduced or even lost the mitogenome. In the three main multicellular lineages of eukaryotes, mitogenomes have pursued diverse evolutionary trajectories in which different types of molecules (circular versus linear and single versus multipartite), gene structures (with or without self-splicing introns), gene contents, gene orders, genetic codes, and transfer RNA editing mechanisms have been selected. Whereas animals have evolved a rather compact mitochondrial genome between 11 and 50 Kb in length with a highly conserved gene content in bilaterians, plants exhibit large mitochondrial genomes of 66 Kb to 11.3 Mb with large intergenic repetitions prone to recombination, and fungal mitogenomes have intermediate sizes of 12 to 236 Kb.

Comparative Metatranscriptomics of Periodontitis Supports a Common Polymicrobial Shift in Metabolic Function and Identifies Novel Putative Disease-Associated ncRNAs

Periodontitis is an inflammatory disease that deteriorates bone supporting teeth afflicting ∼743 million people worldwide. Bacterial communities associated with disease have been classified into red, orange, purple, blue, green, and yellow complexes based on their roles in the periodontal pocket. Previous metagenomic and metatranscriptomics analyses suggest a common shift in metabolic signatures in disease vs. healthy communities with up-regulated processes including pyruvate fermentation, histidine degradation, amino acid metabolism, TonB-dependent receptors. In this work, we examine existing metatranscriptome datasets to identify the commonly differentially expressed transcripts and potential underlying RNA regulatory mechanisms behind the metabolic shifts. Raw RNA-seq reads from three studies (including 49 healthy and 48 periodontitis samples) were assembled into transcripts de novo. Analyses revealed 859 differentially expressed (DE) transcripts, 675 more- and 174 less-expressed. Only ∼20% of the DE transcripts originate from the pathogenic red/orange complexes, and ∼50% originate from organisms unaffiliated with a complex. Comparison of expression profiles revealed variations among disease samples; while specific metabolic processes are commonly up-regulated, the underlying organisms are diverse both within and across disease associated communities. Surveying DE transcripts for known ncRNAs from the Rfam database identified a large number of tRNAs and tmRNAs as well as riboswitches (FMN, glycine, lysine, and SAM) in more prevalent transcripts and the cobalamin riboswitch in both more and less prevalent transcripts. In silico discovery identified many putative ncRNAs in DE transcripts. We report 15 such putative ncRNAs having promising covariation in the predicted secondary structure and interesting genomic context. Seven of these are antisense of ribosomal proteins that are novel and may involve maintaining ribosomal protein stoichiometry during the disease associated metabolic shift. Our findings describe the role of organisms previously unaffiliated with disease and identify the commonality in progression of disease across three metatranscriptomic studies. We find that although the communities are diverse between individuals, the switch in metabolic signatures characteristic of disease is typically achieved through the contributions of several community members. Furthermore, we identify many ncRNAs (both known and putative) which may facilitate the metabolic shifts associated with periodontitis.

RNA base-pairing complexity in living cells visualized by correlated chemical probing

RNA structure and dynamics are critical to biological function. However, strategies for determining RNA structure in vivo are limited, with established chemical probing and newer duplex detection methods each having deficiencies. Here we convert the common reagent dimethyl sulfate into a useful probe of all 4 RNA nucleotides. Building on this advance, we introduce PAIR-MaP, which uses single-molecule correlated chemical probing to directly detect base-pairing interactions in cells. PAIR-MaP has superior resolution compared to alternative experiments, can resolve multiple sets of pairing interactions for structurally dynamic RNAs, and enables highly accurate structure modeling, including of RNAs containing multiple pseudoknots and extensively bound by proteins. Application of PAIR-MaP to human RNase MRP and 2 bacterial messenger RNA 5' untranslated regions reveals functionally important and complex structures undetected by prior analyses. PAIR-MaP is a powerful, experimentally concise, and broadly applicable strategy for directly visualizing RNA base pairs and dynamics in cells.

Deciphering global gene expression and regulation strategy in Escherichia coli during carbon limitation

Despite decades of studies meant to analyse the bacterial response to carbon limitation, we still miss a high-resolution overview of the situation. All gene expression changes observed in such conditions cannot solely be accounted for by the global regulator Crp either free or bound to its effector, cyclic AMP. Here, for the first time, we evaluated the response of both CDS (protein-coding sequence) and ncRNA (non-coding RNA) genes to carbon limitation, revealed cellular functions of differentially expressed genes systematically, quantified the contribution of Crp-cAMP and other factors to regulation and deciphered regulation strategies at a genomewide scale. Approximately one-third of the differentially expressed genes we identified responded to Crp-cAMP via its direct or indirect control, while the remaining genes were subject to growth rate-dependent control or were controlled by other regulators, especially RpoS. Importantly, gene regulation mechanisms can be established by expression pattern studies. Here, we propose a comprehensive picture of how cells respond to carbon scarcity. The global regulation strategies thus exposed illustrate that the response of cell to carbon scarcity is not limited to maintaining sufficient carbon metabolism via cAMP signalling while the main response is to adjust metabolism to cope with a slow growth rate.

rRNA Mimicry in RNA Regulation of Gene Expression

The rRNA is the largest and most abundant RNA in bacterial and archaeal cells. It is also one of the best-characterized RNAs in terms of its structural motifs and sequence variation. Production of ribosome components including >50 ribosomal proteins (r-proteins) consumes significant cellular resources. Thus, RNA cis-regulatory structures that interact with r-proteins to repress further r-protein synthesis play an important role in maintaining appropriate stoichiometry between r-proteins and rRNA. Classically, such mRNA structures were thought to directly mimic the rRNA. However, more than 30 years of research has demonstrated that a variety of different recognition and regulatory paradigms are present. This review will demonstrate how structural mimicry between the rRNA and mRNA cis-regulatory structures may take many different forms. The collection of mRNA structures that interact with r-proteins to regulate r-protein operons are best characterized in Escherichia coli, but are increasingly found within species from nearly all phyla of bacteria and several archaea. Furthermore, they represent a unique opportunity to assess the plasticity of RNA structure in the context of RNA-protein interactions. The binding determinants imposed by r-proteins to allow regulation can be fulfilled in many ways. Some r-protein-interacting mRNAs are immediately obvious as rRNA mimics from primary sequence similarity, others are identifiable only after secondary or tertiary structure determination, and some show no obvious similarity. In addition, across different bacterial species a host of different mechanisms of action have been characterized, showing that there is no simple one-size-fits-all solution.

The Transcriptional landscape of Streptococcus pneumoniae TIGR4 reveals a complex operon architecture and abundant riboregulation critical for growth and virulence

Efficient and highly organized regulation of transcription is fundamental to an organism’s ability to survive, proliferate, and quickly respond to its environment. Therefore, precise mapping of transcriptional units and understanding their regulation is crucial to determining how pathogenic bacteria cause disease and how they may be inhibited. In this study, we map the transcriptional landscape of the bacterial pathogen Streptococcus pneumoniae TIGR4 by applying a combination of high-throughput RNA-sequencing techniques. We successfully map 1864 high confidence transcription termination sites (TTSs), 790 high confidence transcription start sites (TSSs) (742 primary, and 48 secondary), and 1360 low confidence TSSs (74 secondary and 1286 primary) to yield a total of 2150 TSSs. Furthermore, our study reveals a complex transcriptome wherein environment-respondent alternate transcriptional units are observed within operons stemming from internal TSSs and TTSs. Additionally, we identify many putative cis-regulatory RNA elements and riboswitches within 5’-untranslated regions (5’-UTR). By integrating TSSs and TTSs with independently collected RNA-Seq datasets from a variety of conditions, we establish the response of these regulators to changes in growth conditions and validate several of them. Furthermore, to demonstrate the importance of ribo-regulation by 5’-UTR elements for in vivo virulence, we show that the pyrR regulatory element is essential for survival, successful colonization and infection in mice suggesting that such RNA elements are potential drug targets. Importantly, we show that our approach of combining high-throughput sequencing with in vivo experiments can reconstruct a global understanding of regulation, but also pave the way for discovery of compounds that target (ribo-)regulators to mitigate virulence and antibiotic resistance.

The canonical relationship between a bacterial operon and the mRNA transcript produced from the operon has become significantly more complex as numerous regulatory mechanisms that impact the stability, translational efficiency, and early termination rates for mRNA transcripts have been described. With the rise of antibiotic resistance, these mechanisms offer new potential targets for antibiotic development. In this study we used a combination of high-throughput sequencing technologies to assess genome-wide transcription start and stop sites, as well as determine condition specific global transcription patterns in the human pathogen Streptococcus pneumoniae. We find that the majority of multi-gene operons have alternative start and stop sites enabling condition specific regulation of genes within the same operon. Furthermore, we identified many putative RNA regulators that are widespread in the S. pneumoniae pan-genome. Finally, we show that separately collected RNA-Seq data enables identification of conditional triggers for regulatory RNAs, and experimentally demonstrate that our approach may be used to identify drug-able RNA targets by establishing that pyrR RNA functionality is critical for successful S. pneumoniae mouse colonization and infection. Thus, our study not only uses genome-wide high-throughput approaches to identify putative RNA regulators, but also establishes the importance of such regulators in S. pneumoniae virulence.

Ribosome Profiling of Synechocystis Reveals Altered Ribosome Allocation at Carbon Starvation

Ribosome profiling accesses the translational step of gene expression via deep sequencing of ribosome-protected mRNA footprints. Pairing of ribosome profiling and transcriptomics data provides a translational efficiency for each gene. Here, the translatome and transcriptome of the model cyanobacterium Synechocystis were compared under carbon-replete and carbon starvation conditions. The latter may be experienced when cyanobacteria are cultivated in poorly mixed bioreactors or engineered to be product-secreting cell factories. A small fraction of genes (<200), including stress response genes, showed changes in translational efficiency during carbon starvation, indicating condition-dependent translation-level regulation. We observed ribosome occupancy in untranslated regions, possibly due to an alternative translation initiation mechanism in Synechocystis. The higher proportion of ribosomes residing in untranslated regions during carbon starvation may be a mechanism to quickly inactivate superfluous ribosomes. This work provides the first ribosome profiling data for cyanobacteria and reveals new regulation strategies for coping with nutrient limitation.

Fitness advantages conferred by the L20-interacting RNA cis-regulator of ribosomal protein synthesis in Bacillus subtilis

In many bacteria, ribosomal proteins autogenously repress their own expression by interacting with RNA structures typically located in the 5'-UTRs of their mRNA transcripts. This regulation is necessary to maintain a balance between ribosomal proteins and rRNA to ensure proper ribosome production. Despite advances in noncoding RNA discovery and validation of RNA-protein regulatory interactions, the selective pressures that govern the formation and maintenance of such RNA cis-regulators in the context of an organism remain largely undetermined. To examine the impact disruptions to this regulation have on bacterial fitness, we introduced point mutations that abolish ribosomal protein binding and regulation into the RNA structure that controls expression of ribosomal proteins L20 and L35 within the Bacillus subtilis genome. Our studies indicate that removing this regulation results in reduced log phase growth, improper rRNA maturation, and the accumulation of a kinetically trapped or misassembled ribosomal particle at low temperatures, suggesting defects in ribosome synthesis. Such work emphasizes the important role regulatory RNAs play in the stoichiometric production of ribosomal components for proper ribosome composition and overall organism viability and reinforces the potential of targeting ribosomal protein production and ribosome assembly with novel antimicrobials.

Pervasive Regulatory Functions of mRNA Structure Revealed by High-Resolution SHAPE Probing

mRNAs can fold into complex structures that regulate gene expression. Resolving such structures de novo has remained challenging and has limited our understanding of the prevalence and functions of mRNA structure. We use SHAPE-MaP experiments in living E. coli cells to derive quantitative, nucleotide-resolution structure models for 194 endogenous transcripts encompassing approximately 400 genes. Individual mRNAs have exceptionally diverse architectures, and most contain well-defined structures. Active translation destabilizes mRNA structure in cells. Nevertheless, mRNA structure remains similar between in-cell and cell-free environments, indicating broad potential for structure-mediated gene regulation. We find that the translation efficiency of endogenous genes is regulated by unfolding kinetics of structures overlapping the ribosome binding site. We discover conserved structured elements in 35% of UTRs, several of which we validate as novel protein binding motifs. RNA structure regulates every gene studied here in a meaningful way, implying that most functional structures remain to be discovered.

Challenges of ligand identification for the second wave of orphan riboswitch candidates

Orphan riboswitch candidates are noncoding RNA motifs whose representatives are believed to function as genetic regulatory elements, but whose target ligands have yet to be identified. The study of certain orphans, particularly classes that have resisted experimental validation for many years, has led to the discovery of important biological pathways and processes once their ligands were identified. Previously, we highlighted details for four of the most common and intriguing orphan riboswitch candidates. This facilitated the validation of riboswitches for the signaling molecules c-di-AMP, ZTP, and ppGpp, the metal ion Mn2+, and the metabolites guanidine and PRPP. Such studies also yield useful linkages between the ligands sensed by the riboswitches and numerous biochemical pathways. In the current report, we describe the known characteristics of 30 distinct classes of orphan riboswitch candidates - some of which have remained unsolved for over a decade. We also discuss the prospects for uncovering novel biological insights via focused studies on these RNAs. Lastly, we make recommendations for experimental objectives along the path to finding ligands for these mysterious RNAs.

Identification of sRNA mediated responses to nutrient depletion in Burkholderia pseudomallei

The Burkholderia genus includes many species that are known to survive in diverse environmental conditions including low nutrient environments. One species, Burkholderia pseudomallei is a versatile pathogen that can survive in a wide range of hosts and environmental conditions. In this study, we investigated how a nutrient depleted growth environment evokes sRNA mediated responses by B. pseudomallei. Computationally predicted B. pseudomallei D286 sRNAs were mapped to RNA-sequencing data for cultures grown under two conditions: (1) BHIB as a nutrient rich media reference environment and (2) M9 media as a nutrient depleted stress environment. The sRNAs were further selected to identify potentially cis-encoded systems by investigating their possible interactions with their flanking genes. The mappings of predicted sRNA genes and interactions analysis to their flanking genes identified 12 sRNA candidates that may possibly have cis-acting regulatory roles that are associated to a nutrient depleted growth environment. Our approach can be used for identifying novel sRNA genes and their possible role as cis-mediated regulatory systems.

Genome-wide mRNA processing in methanogenic archaea reveals post-transcriptional regulation of ribosomal protein synthesis

Unlike stable RNAs that require processing for maturation, prokaryotic cellular mRNAs generally follow an 'all-or-none' pattern. Herein, we used a 5΄ monophosphate transcript sequencing (5΄P-seq) that specifically captured the 5΄-end of processed transcripts and mapped the genome-wide RNA processing sites (PSSs) in a methanogenic archaeon. Following statistical analysis and stringent filtration, we identified 1429 PSSs, among which 23.5% and 5.4% were located in 5΄ untranslated region (uPSS) and intergenic region (iPSS), respectively. A predominant uridine downstream PSSs served as a processing signature. Remarkably, 5΄P-seq detected overrepresented uPSS and iPSS in the polycistronic operons encoding ribosomal proteins, and the majority upstream and proximal ribosome binding sites, suggesting a regulatory role of processing on translation initiation. The processed transcripts showed increased stability and translation efficiency. Particularly, processing within the tricistronic transcript of rplA-rplJ-rplL enhanced the translation of rplL, which can provide a driving force for the 1:4 stoichiometry of L10 to L12 in the ribosome. Growth-associated mRNA processing intensities were also correlated with the cellular ribosomal protein levels, thereby suggesting that mRNA processing is involved in tuning growth-dependent ribosome synthesis. In conclusion, our findings suggest that mRNA processing-mediated post-transcriptional regulation is a potential mechanism of ribosomal protein synthesis and stoichiometry.

Identification of 15 candidate structured noncoding RNA motifs in fungi by comparative genomics

BACKGROUND

With the development of rapid and inexpensive DNA sequencing, the genome sequences of more than 100 fungal species have been made available. This dataset provides an excellent resource for comparative genomics analyses, which can be used to discover genetic elements, including noncoding RNAs (ncRNAs). Bioinformatics tools similar to those used to uncover novel ncRNAs in bacteria, likewise, should be useful for searching fungal genomic sequences, and the relative ease of genetic experiments with some model fungal species could facilitate experimental validation studies.

RESULTS

We have adapted a bioinformatics pipeline for discovering bacterial ncRNAs to systematically analyze many fungal genomes. This comparative genomics pipeline integrates information on conserved RNA sequence and structural features with alternative splicing information to reveal fungal RNA motifs that are candidate regulatory domains, or that might have other possible functions. A total of 15 prominent classes of structured ncRNA candidates were identified, including variant HDV self-cleaving ribozyme representatives, atypical snoRNA candidates, and possible structured antisense RNA motifs. Candidate regulatory motifs were also found associated with genes for ribosomal proteins, S-adenosylmethionine decarboxylase (SDC), amidase, and HexA protein involved in Woronin body formation. We experimentally confirm that the variant HDV ribozymes undergo rapid self-cleavage, and we demonstrate that the SDC RNA motif reduces the expression of SAM decarboxylase by translational repression. Furthermore, we provide evidence that several other motifs discovered in this study are likely to be functional ncRNA elements.

CONCLUSIONS

Systematic screening of fungal genomes using a computational discovery pipeline has revealed the existence of a variety of novel structured ncRNAs. Genome contexts and similarities to known ncRNA motifs provide strong evidence for the biological and biochemical functions of some newly found ncRNA motifs. Although initial examinations of several motifs provide evidence for their likely functions, other motifs will require more in-depth analysis to reveal their functions.

Community-like genome in single cells of the sulfur bacterium Achromatium oxaliferum

Polyploid bacteria are common, but the genetic and functional diversity resulting from polyploidy is unknown. Here we use single-cell genomics, metagenomics, single-cell amplicon sequencing, and fluorescence in situ hybridization, to show that individual cells of Achromatium oxaliferum, the world’s biggest known freshwater bacterium, harbor genetic diversity typical of whole bacterial communities. The cells contain tens of transposable elements, which likely cause the unprecedented diversity that we observe in the sequence and synteny of genes. Given the high within-cell diversity of the usually conserved 16S ribosomal RNA gene, we suggest that gene conversion occurs in multiple, separated genomic hotspots. The ribosomal RNA distribution inside the cells hints to spatially differential gene expression. We also suggest that intracellular gene transfer may lead to extensive gene reshuffling and increased diversity.

The cells of Achromatium bacteria are remarkably large and contain multiple chromosome copies. Here, Ionescu et al. show that chromosome copies within individual cells display high diversity, similar to that of bacterial communities, and contain tens of transposable elements.

Recognizing RNA structural motifs in HT-SELEX data for ribosomal protein S15

BACKGROUND

Proteins recognize many different aspects of RNA ranging from single stranded regions to discrete secondary or tertiary structures. High-throughput sequencing (HTS) of in vitro selected populations offers a large scale method to study RNA-proteins interactions. However, most existing analysis methods require that the binding motifs are enriched in the population relative to earlier rounds, and that motifs are found in a loop or single stranded region of the potential RNA secondary structure. Such methods do not generalize to all RNA-protein interaction as some RNA binding proteins specifically recognize more complex structures such as double stranded RNA.

RESULTS

In this study, we use HT-SELEX derived populations to study the landscape of RNAs that interact with Geobacillus kaustophilus ribosomal protein S15. Our data show high sequence and structure diversity and proved intractable to existing methods. Conventional programs identified some sequence motifs, but these are found in less than 5-10% of the total sequence pool. Therefore, we developed a novel framework to analyze HT-SELEX data. Our process accounts for both sequence and structure components by abstracting the overall secondary structure into smaller substructures composed of a single base-pair stack, which allows us to leverage existing approaches already used in k-mer analysis to identify enriched motifs. By focusing on secondary structure motifs composed of specific two base-pair stacks, we identified significantly enriched or depleted structure motifs relative to earlier rounds.

CONCLUSIONS

Discrete substructures are likely to be important to RNA-protein interactions, but they are difficult to elucidate. Substructures can help make highly diverse sequence data more tractable. The structure motifs provide limited accuracy in predicting enrichment suggesting that G. kaustophilus S15 can either recognize many different secondary structure motifs or some aspects of the interaction are not captured by the analysis. This highlights the importance of considering secondary and tertiary structure elements and their role in RNA-protein interactions.

A statistical test for conserved RNA structure shows lack of evidence for structure in lncRNAs

Many functional RNAs have an evolutionarily conserved secondary structure. Conservation of RNA base pairing induces pairwise covariations in sequence alignments. We developed a computational method, R-scape (RNA Structural Covariation Above Phylogenetic Expectation), that quantitatively tests whether covariation analysis supports the presence of a conserved RNA secondary structure. R-scape analysis finds no statistically significant support for proposed secondary structures of the long noncoding RNAs HOTAIR, SRA, and Xist.

Regulation of Ribosomal Protein Operons rplM-rpsI, rpmB-rpmG, and rplU-rpmA at the Transcriptional and Translational Levels

It is widely assumed that in the best-characterized model bacterium Escherichia coli, transcription units encoding ribosomal proteins (r-proteins) and regulation of their expression have been already well defined. However, transcription start sites for several E. coli r-protein operons have been established only very recently, so that information concerning the regulation of these operons at the transcriptional or posttranscriptional level is still missing. This paper describes for the first time the in vivo regulation of three r-protein operons, rplM-rpsI, rpmB-rpmG, and rplU-rpmA The results demonstrate that transcription of all three operons is subject to ppGpp/DksA-dependent negative stringent control under amino acid starvation, in parallel with the rRNA operons. By using single-copy translational fusions with the chromosomal lacZ gene, we show here that at the translation level only one of these operons, rplM-rpsI, is regulated by the mechanism of autogenous repression involving the 5' untranslated region (UTR) of the operon mRNA, while rpmB-rpmG and rplU-rpmA are not subject to this type of regulation. This may imply that translational feedback control is not a general rule for modulating the expression of E. coli r-protein operons. Finally, we report that L13, a primary protein in 50S ribosomal subunit assembly, serves as a repressor of rplM-rpsI expression in vivo, acting at a target within the rplM translation initiation region. Thus, L13 represents a novel example of regulatory r-proteins in bacteria.

IMPORTANCE

It is important to obtain a deeper understanding of the regulatory mechanisms responsible for coordinated and balanced synthesis of ribosomal components. In this paper, we highlight the major role of a stringent response in regulating transcription of three previously unexplored r-protein operons, and we show that only one of them is subject to feedback regulation at the translational level. Improved knowledge of the regulatory pathways controlling ribosome biogenesis may promote the development of novel antibacterial agents.

RNA regulators responding to ribosomal protein S15 are frequent in sequence space

There are several natural examples of distinct RNA structures that interact with the same ligand to regulate the expression of homologous genes in different organisms. One essential question regarding this phenomenon is whether such RNA regulators are the result of convergent or divergent evolution. Are the RNAs derived from some common ancestor and diverged to the point where we cannot identify the similarity, or have multiple solutions to the same biological problem arisen independently? A key variable in assessing these alternatives is how frequently such regulators arise within sequence space. Ribosomal protein S15 is autogenously regulated via an RNA regulator in many bacterial species; four apparently distinct regulators have been functionally validated in different bacterial phyla. Here, we explore how frequently such regulators arise within a partially randomized sequence population. We find many RNAs that interact specifically with ribosomal protein S15 from Geobacillus kaustophilus with biologically relevant dissociation constants. Furthermore, of the six sequences we characterize, four show regulatory activity in an Escherichia coli reporter assay. Subsequent footprinting and mutagenesis analysis indicates that protein binding proximal to regulatory features such as the Shine–Dalgarno sequence is sufficient to enable regulation, suggesting that regulation in response to S15 is relatively easily acquired.

Studying RNA Homology and Conservation with Infernal: From Single Sequences to RNA Families

Emerging high-throughput technologies have led to a deluge of putative non-coding RNA (ncRNA) sequences identified in a wide variety of organisms. Systematic characterization of these transcripts will be a tremendous challenge. Homology detection is critical to making maximal use of functional information gathered about ncRNAs: identifying homologous sequence allows us to transfer information gathered in one organism to another quickly and with a high degree of confidence. ncRNA presents a challenge for homology detection, as the primary sequence is often poorly conserved and de novo secondary structure prediction and search remain difficult. This unit introduces methods developed by the Rfam database for identifying "families" of homologous ncRNAs starting from single "seed" sequences, using manually curated sequence alignments to build powerful statistical models of sequence and structure conservation known as covariance models (CMs), implemented in the Infernal software package. We provide a step-by-step iterative protocol for identifying ncRNA homologs and then constructing an alignment and corresponding CM. We also work through an example for the bacterial small RNA MicA, discovering a previously unreported family of divergent MicA homologs in genus Xenorhabdus in the process. © 2016 by John Wiley & Sons, Inc.

The role of mRNA structure in bacterial translational regulation

The characteristics of bacterial messenger RNAs (mRNAs) that influence translation efficiency provide many convenient handles for regulation of gene expression, especially when coupled with the processes of transcription termination and mRNA degradation. An mRNA's structure, especially near the site of initiation, has profound consequences for how readily it is translated. This property allows bacterial gene expression to be altered by changes to mRNA structure induced by temperature, or interactions with a wide variety of cellular components including small molecules, other RNAs (such as sRNAs and tRNAs), and RNA-binding proteins. This review discusses the links between mRNA structure and translation efficiency, and how mRNA structure is manipulated by conditions and signals within the cell to regulate gene expression. The range of RNA regulators discussed follows a continuum from very complex tertiary structures such as riboswitch aptamers and ribosomal protein-binding sites to thermosensors and mRNA:sRNA interactions that involve only base-pairing interactions. Furthermore, the high degrees of diversity observed for both mRNA structures and the mechanisms by which inhibition of translation occur have significant consequences for understanding the evolution of bacterial translational regulation. WIREs RNA 2017, 8:e1370. doi: 10.1002/wrna.1370 For further resources related to this article, please visit the WIREs website.

Involvement of ribosomal protein L6 in assembly of functional 50S ribosomal subunit in Escherichia coli cells

Ribosomal protein L6, an essential component of the large (50S) subunit, primarily binds to helix 97 of 23S rRNA and locates near the sarcin/ricin loop of helix 95 that directly interacts with GTPase translation factors. Although L6 is believed to play important roles in factor-dependent ribosomal function, crucial biochemical evidence for this hypothesis has not been obtained. We constructed and characterized an Escherichia coli mutant bearing a chromosomal L6 gene (rplF) disruption and carrying a plasmid with an arabinose-inducible L6 gene. Although this ΔL6 mutant grew more slowly than its wild-type parent, it proliferated in the presence of arabinose. Interestingly, cell growth in the absence of arabinose was biphasic. Early growth lasted only a few generations (LI-phase) and was followed by a suspension of growth for several hours (S-phase). This suspension was followed by a second growth phase (LII-phase). Cells harvested at both LI- and S-phases contained ribosomes with reduced factor-dependent GTPase activity and accumulated 50S subunit precursors (45S particles). The 45S particles completely lacked L6. Complete 50S subunits containing L6 were observed in all growth phases regardless of the L6-depleted condition, implying that the ΔL6 mutant escaped death because of a leaky expression of L6 from the complementing plasmid. We conclude that L6 is essential for the assembly of functional 50S subunits at the late stage. We thus established conditions for the isolation of L6-depleted 50S subunits, which are essential to study the role of L6 in translation.

Loss of Conserved Noncoding RNAs in Genomes of Bacterial Endosymbionts

The genomes of intracellular symbiotic or pathogenic bacteria, such as of Buchnera, Mycoplasma, and Rickettsia, are typically smaller compared with their free-living counterparts. Here we showed that noncoding RNA (ncRNA) families, which are conserved in free-living bacteria, frequently could not be detected by computational methods in the small genomes. Statistical tests demonstrated that their absence is not an artifact of low GC content or small deletions in these small genomes, and thus it was indicative of an independent loss of ncRNAs in different endosymbiotic lineages. By analyzing the synteny (conservation of gene order) between the reduced and nonreduced genomes, we revealed instances of protein-coding genes that were preserved in the reduced genomes but lost cis-regulatory elements. We found that the loss of cis-regulatory ncRNA sequences, which regulate the expression of cognate protein-coding genes, is characterized by the reduction of secondary structure formation propensity, GC content, and length of the corresponding genomic regions.

Abundant Intergenic TAACTGA Direct Repeats and Putative Alternate RNA Polymerase β' Subunits in Marine Beggiatoaceae Genomes: Possible Regulatory Roles and Origins

The genome sequences of several giant marine sulfur-oxidizing bacteria present evidence of a possible post-transcriptional regulatory network that may have been transmitted to or from two distantly related bacteria lineages. The draft genome of a Cand. “Maribeggiatoa” filament from the Guaymas Basin (Gulf of California, Mexico) seafloor contains 169 sets of TAACTGA direct repeats and one indirect repeat, with two to six copies per set. Related heptamers are rarely or never found as direct repeats. TAACTGA direct repeats are also found in some other Beggiatoaceae, Thiocystis violascens, a range of Cyanobacteria, and five Bacteroidetes. This phylogenetic distribution suggests they may have been transmitted horizontally, but no mechanism is evident. There is no correlation between total TAACTGA occurrences and repeats per genome. In most species the repeat units are relatively short, but longer arrays of up to 43 copies are found in several Bacteroidetes and Cyanobacteria. The majority of TAACTGA repeats in the Cand. “Maribeggiatoa” Orange Guaymas (BOGUAY) genome are within several nucleotides upstream of a putative start codon, suggesting they may be binding sites for a post-transcriptional regulator. Candidates include members of the ribosomal protein S1, Csp (cold shock protein), and Csr (carbon storage regulator) families. No pattern was evident in the predicted functions of the open reading frames (ORFs) downstream of repeats, but some encode presumably essential products such as ribosomal proteins. Among these is an ORF encoding a possible alternate or modified RNA polymerase beta prime subunit, predicted to have the expected subunit interaction domains but lacking most catalytic residues. A similar ORF was found in the Thioploca ingrica draft genome, but in no others. In both species they are immediately upstream of putative sensor kinase genes with nearly identical domain structures. In the marine Beggiatoaceae, a role for the TAACTGA repeats in translational regulation is suggested. More speculatively, the putative alternate RNA polymerase subunit could be a negative transcriptional regulator.

Co-evolution of Bacterial Ribosomal Protein S15 with Diverse mRNA Regulatory Structures

RNA-protein interactions are critical in many biological processes, yet how such interactions affect the evolution of both partners is still unknown. RNA and protein structures are impacted very differently by mechanisms of genomic change. While most protein families are identifiable at the nucleotide level across large phylogenetic distances, RNA families display far less nucleotide similarity and are often only shared by closely related bacterial species. Ribosomal protein S15 has two RNA binding functions. First, it is a ribosomal protein responsible for organizing the rRNA during ribosome assembly. Second, in many bacterial species S15 also interacts with a structured portion of its own transcript to negatively regulate gene expression. While the first interaction is conserved in most bacteria, the second is not. Four distinct mRNA structures interact with S15 to enable regulation, each of which appears to be independently derived in different groups of bacteria. With the goal of understanding how protein-binding specificity may influence the evolution of such RNA regulatory structures, we examine whether examples of these mRNA structures are able to interact with, and regulate in response to, S15 homologs from organisms containing distinct mRNA structures. We find that despite their shared RNA binding function in the rRNA, S15 homologs have distinct RNA recognition profiles. We present a model to explain the specificity patterns observed, and support this model by with further mutagenesis. After analyzing the patterns of conservation for the S15 protein coding sequences, we also identified amino acid changes that alter the binding specificity of an S15 homolog. In this work we demonstrate that homologous RNA-binding proteins have different specificity profiles, and minor changes to amino acid sequences, or to RNA structural motifs, can have large impacts on RNA-protein recognition.

An S6:S18 complex inhibits translation of E. coli rpsF

More than half of the ribosomal protein operons in Escherichia coli are regulated by structures within the mRNA transcripts that interact with specific ribosomal proteins to inhibit further protein expression. This regulation is accomplished using a variety of mechanisms and the RNA structures responsible for regulation are often not conserved across bacterial phyla. A widely conserved mRNA structure preceding the ribosomal protein operon containing rpsF and rpsR (encoding S6 and S18) was recently identified through comparative genomics. Examples of this RNA from both E. coli and Bacillus subtilis were shown to interact in vitro with an S6:S18 complex. In this work, we demonstrate that in E. coli, this RNA structure regulates gene expression in response to the S6:S18 complex. β-galactosidase activity from a lacZ reporter translationally fused to the 5' UTR and first nine codons of E. coli rpsF is reduced fourfold by overexpression of a genomic fragment encoding both S6 and S18 but not by overexpression of either protein individually. Mutations to the mRNA structure, as well as to the RNA-binding site of S18 and the S6-S18 interaction surfaces of S6 and S18, are sufficient to derepress β-galactosidase activity, indicating that the S6:S18 complex is the biologically active effector. Measurement of transcript levels shows that although reporter levels do not change upon protein overexpression, levels of the native transcript are reduced fourfold, suggesting that the mRNA regulator prevents translation and this effect is amplified on the native transcript by other mechanisms.