RNA structures regulating ribosomal protein biosynthesis in bacilli

Regulatory RNAs in Bacillus subtilis: A review on regulatory mechanism and applications in synthetic biology

Bacteria exhibit a rich repertoire of RNA molecules that intricately regulate gene expression at multiple hierarchical levels, including small RNAs (sRNAs), riboswitches, and antisense RNAs. Notably, the majority of these regulatory RNAs lack or have limited protein-coding capacity but play pivotal roles in orchestrating gene expression by modulating transcription, post-transcription or translation processes. Leveraging and redesigning these regulatory RNA elements have emerged as pivotal strategies in the domains of metabolic engineering and synthetic biology. While previous investigations predominantly focused on delineating the roles of regulatory RNA in Gram-negative bacterial models such as Escherichia coli and Salmonella enterica, this review aims to summarize the mechanisms and functionalities of endogenous regulatory RNAs inherent to typical Gram-positive bacteria, notably Bacillus subtilis. Furthermore, we explore the engineering and practical applications of these regulatory RNA elements in the arena of synthetic biology, employing B. subtilis as a foundational chassis.

Unraveling the adaptive strategies of Mycoplasma hominis through proteogenomic profiling of clinical isolates

Introduction

Mycoplasma hominis (M. hominis) belongs to the class Mollicutes, characterized by a very small genome size, reduction of metabolic pathways, including transcription factors, and the absence of a cell wall. Despite this, they adapt well not only to specific niches within the host organism but can also spread throughout the body, colonizing various organs and tissues. The adaptation mechanisms of M. hominis, as well as their regulatory pathways, are poorly understood. It is known that, when adapting to adverse conditions, Mycoplasmas can undergo phenotypic switches that may persist for several generations.

Methods

To investigate the adaptive properties of M. hominis related to survival in the host, we conducted a comparative phenotypic and proteogenomic analysis of eight clinical isolates of M. hominis obtained from patients with urogenital infections and the laboratory strain H-34.

Results

We have shown that clinical isolates differ in phenotypic features from the laboratory strain, form biofilms more effectively and show resistance to ofloxacin. The comparative proteogenomic analysis revealed that, unlike the laboratory strain, the clinical isolates possess several features related to stress survival: they switch carbon metabolism, activating the energetically least advantageous pathway of nucleoside utilization, which allows slowing down cellular processes and transitioning to a starvation state; they reconfigure the repertoire of membrane proteins; they have integrative conjugative elements in their genomes, which are key mediators of horizontal gene transfer. The upregulation of the methylating subunit of the restriction-modification (RM) system type I and the additional components of RM systems found in clinical isolates suggest that DNA methylation may play a role in regulating the adaptation mechanisms of M. hominis in the host organism. It has been shown that based on the proteogenomic profile, namely the genome sequence, protein content, composition of the RM systems and additional subunits HsdM, HsdS and HsdR, composition and number of transposable elements, as well as the sequence of the main variable antigen Vaa, we can divide clinical isolates into two phenotypes: typical colonies (TC), which have a high growth rate, and atypical (aTC) mini-colonies, which have a slow growth rate and exhibit properties similar to persisters.

Discussion

We believe that the key mechanism of adaptation of M. hominis in the host is phenotypic restructuring, leading to a slowing down cellular processes and the formation of small atypical colonies. This is due to a switch in carbon metabolism and activation the pathway of nucleoside utilization. We hypothesize that DNA methylation may play a role in regulating this switch.

Evaluation of suitable reference genes for expression profile analyses of target genes in the coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Curculionidae)

The coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Curculionidae), is a major destructive insect pest of coffee, which impacts the coffee crops negatively. As a draft genome has been completed for this insect, most molecular studies on gene transcriptional levels under different experimental conditions will be conducted using real-time reverse-transcription quantitative polymerase chain reactions (RT-qPCR). However, the lack of suitable internal reference genes will affect the accuracy of RT-qPCR results. In this study, the expression stability of nine candidate reference genes was evaluated under different developmental stages, temperature stress, and Beauveria bassiana infection. Data analyses were completed by four commonly used programs, BestKeeper, NormFinder, geNorm, and RefFinder. The result showed that RPL3 and EF1α combination were recommended as the most stable reference genes for developmental stages. EF1α and RPS3a combination were the top two stable reference genes for B. bassiana infection. RPS3a and RPL3 combination performed as the optimal reference genes both in temperature stress and all samples. Our results should provide a good foundation for the expression profile analyses of target genes in the future, especially for molecular studies on insect genetic development, temperature adaptability, and immune mechanism to entomogenous fungi in H. hampei.

Actin remodeling driven by circLIMA1: sperm cell as an intriguing cellular model

CircRNA cargo in spermatozoa (SPZ) participates in setting cell quality, in terms of morphology and motility. Cannabinoid receptor CB1 activity is correlated with a proper spermatogenesis and epididymal sperm maturation. Despite CB1 promotes endogenous skill to circularize mRNAs in SPZ, few notions are reported regarding the functional link between endocannabinoids and spermatic circRNA cargo. In CB1 knock-out male mice, we performed a complete dataset of spermatic circRNA content by microarray strategy.

Differentially expressed (DE)-circRNAs, as a function of genotype, were identified. Within DE-circRNAs, we focused the attention on circLIMA1, as putative actin-cytoskeleton architecture regulator. The validation of circLIMA1 dependent-competitive endogenous RNA (ceRNA) network (ceRNET) in in vitro cell line confirmed its activity in the regulation of the cytoskeletal actin. Interestingly, a dynamic actin regulation in SPZ nuclei was found during their epididymal maturation.

In this scenario, we showed for the first time an intriguing sperm nuclear actin remodeling, regulated via a ceRNET-independent pathway, consisting in the nuclear shuttling of circLIMA1-QKI interactome and downstream in Gelsolin regulation. In particular, the increased levels of circLIMA1 in CB1 knock-out SPZ, associated with an inefficient depolymerization of nuclear actin, specifically illustrate how endocannabinoids, by regulating circRNA cargo, may contribute to sperm morpho-cellular maturation.

The Dynamics of Mycoplasma gallisepticum Nucleoid Structure at the Exponential and Stationary Growth Phases

The structure and dynamics of bacterial nucleoids play important roles in regulating gene expression. Bacteria of class Mollicutes and, in particular, mycoplasmas feature extremely reduced genomes. They lack multiple structural proteins of the nucleoid, as well as regulators of gene expression. We studied the organization of Mycoplasma gallisepticum nucleoids in the stationary and exponential growth phases at the structural and protein levels. The growth phase transition results in the structural reorganization of M. gallisepticum nucleoid. In particular, it undergoes condensation and changes in the protein content. The observed changes corroborate with the previously identified global rearrangement of the transcriptional landscape in this bacterium during the growth phase transition. In addition, we identified that the glycolytic enzyme enolase functions as a nucleoid structural protein in this bacterium. It is capable of non-specific DNA binding and can form fibril-like complexes with DNA.

Discovery of 17 conserved structural RNAs in fungi

Many non-coding RNAs with known functions are structurally conserved: their intramolecular secondary and tertiary interactions are maintained across evolutionary time. Consequently, the presence of conserved structure in multiple sequence alignments can be used to identify candidate functional non-coding RNAs. Here, we present a bioinformatics method that couples iterative homology search with covariation analysis to assess whether a genomic region has evidence of conserved RNA structure. We used this method to examine all unannotated regions of five well-studied fungal genomes (Saccharomyces cerevisiae, Candida albicans, Neurospora crassa, Aspergillus fumigatus, and Schizosaccharomyces pombe). We identified 17 novel structurally conserved non-coding RNA candidates, which include four H/ACA box small nucleolar RNAs, four intergenic RNAs and nine RNA structures located within the introns and untranslated regions (UTRs) of mRNAs. For the two structures in the 3' UTRs of the metabolic genes GLY1 and MET13, we performed experiments that provide evidence against them being eukaryotic riboswitches.

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.

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.

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.

Genome-wide discovery of structured noncoding RNAs in bacteria

Background

Structured noncoding RNAs (ncRNAs) play essential roles in many biological processes such as gene regulation, signaling, RNA processing, and protein synthesis. Among the most common groups of ncRNAs in bacteria are riboswitches. These cis-regulatory, metabolite-binding RNAs are present in many species where they regulate various metabolic and signaling pathways. Collectively, there are likely to be hundreds of novel riboswitch classes that remain hidden in the bacterial genomes that have already been sequenced, and potentially thousands of classes distributed among various other species in the biosphere. The vast majority of these undiscovered classes are proposed to be exceedingly rare, and so current bioinformatics search techniques are reaching their limits for differentiating between true riboswitch candidates and false positives.

Results

Herein, we exploit a computational search pipeline that can efficiently identify intergenic regions most likely to encode structured ncRNAs. Application of this method to five bacterial genomes yielded nearly 70 novel genetic elements including 30 novel candidate ncRNA motifs. Among the riboswitch candidates identified is an RNA motif involved in the regulation of thiamin biosynthesis.

Conclusions

Analysis of other genomes will undoubtedly lead to the discovery of many additional novel structured ncRNAs, and provide insight into the range of riboswitches and other kinds of ncRNAs remaining to be discovered in bacteria and archaea.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1433-7) contains supplementary material, which is available to authorized users.

Conserved and specific features of Streptococcus pyogenes and Streptococcus agalactiae transcriptional landscapes

BACKGROUND

The human pathogen Streptococcus pyogenes, or group A Streptococcus, is responsible for mild infections to life-threatening diseases. To facilitate the characterization of regulatory networks involved in the adaptation of this pathogen to its different environments and their evolution, we have determined the primary transcriptome of a serotype M1 S. pyogenes strain at single-nucleotide resolution and compared it with that of Streptococcus agalactiae, also from the pyogenic group of streptococci.

RESULTS

By using a combination of differential RNA-sequencing and oriented RNA-sequencing we have identified 892 transcription start sites (TSS) and 885 promoters in the S. pyogenes M1 strain S119. 8.6% of S. pyogenes mRNAs were leaderless, among which 81% were also classified as leaderless in S. agalactiae. 26% of S. pyogenes transcript 5' untranslated regions (UTRs) were longer than 60 nt. Conservation of long 5' UTRs with S. agalactiae allowed us to predict new potential regulatory sequences. In addition, based on the mapping of 643 transcript ends in the S. pyogenes strain S119, we constructed an operon map of 401 monocistrons and 349 operons covering 81.5% of the genome. One hundred fifty-six operons and 254 monocistrons retained the same organization, despite multiple genomic reorganizations between S. pyogenes and S. agalactiae. Genomic reorganization was found to more often go along with variable promoter sequences and 5' UTR lengths. Finally, we identified 117 putative regulatory RNAs, among which nine were regulated in response to magnesium concentration.

CONCLUSIONS

Our data provide insights into transcriptome evolution in pyogenic streptococci and will facilitate the analysis of genetic polymorphisms identified by comparative genomics in S. pyogenes.

Selection of Reference Genes for Optimal Normalization of Quantitative Real-Time Polymerase Chain Reaction Results for Diaphorina citri Adults

The Asian citrus psyllid, Diaphorina citri Kuwayama (Hemiptera: Liviidae), can cause direct damage to citrus trees and is the main vector for the devastating disease, citrus greening disease or huanglongbing. Most molecular studies on this important insect pest use real-time reverse-transcription quantitative polymerase chain reaction (RT-qPCR) to quantify gene expression, including analyzing molecular basis for insecticide resistance in field populations. One critical factor to cause inaccuracy in RT-qPCR results is the lack of appropriate internal reference genes for optimal data normalization. In this study, the expression levels of 10 selected reference genes were evaluated in different tissue samples of psyllid adults and in the insects treated with different temperatures and insecticides. Data were analyzed using different computational algorithms, including Delta Ct, BestKeeper, NormFinder, geNorm, and RefFinder. According to our results, at least two reference genes should be used for the normalization of RT-qPCR data in this insect. The best choices of reference genes for different samples are as follows: ACT1 and Ferritin for different tissue samples, RPS20 and Ferritin for samples treated with different temperatures, TBP and EF1α for samples treated with imidacloprid, and Ferritin and TBP for samples treated with beta-cypermethrin. The reference genes identified in this study should be useful for future studies to analyze the expression patterns of target genes, especially for genes linked with temperature adaptability and insecticide resistance in this insect species in the future.

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.

Disruption of a stem-loop structure located upstream of pseudoknot domain in Tobacco mosaic virus enhanced its infectivity and viral RNA accumulation

A predicted stem-loop structure of 25 nucleotides, located in the coat protein (CP) gene and 3'-UTR sequences of Tobacco mosaic virus (TMV), was validated previously (Guo et al., 2015). In this study, both disrupted stem-loop and nucleotide deletion mutants of TMV replicated more rapidly in Nicotiana benthamiana protoplasts. The TMV mutant with a complete mirrored stem-loop structure showed similar level of viral RNA accumulation as TMV. Recovering the stem-loop structure also resulted in a similar replication level as TMV. All these mutants induced necrosis in N. benthamiana and assembled into typical rigid rod-shaped virions. TMV mutant without the stem-loop structure induced more local lesions in Chenopodium quinoa. When the putative stem-loop structure in Tomato mosaic virus (ToMV) was disrupted, the mutant also showed an enhanced virus replication. This suggests that the stem-loop structure of TMV is a new cis-acting element with a role in virus replication.

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.

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.

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.

Regulatory RNAs in Bacillus subtilis: a Gram-Positive Perspective on Bacterial RNA-Mediated Regulation of Gene Expression

Bacteria can employ widely diverse RNA molecules to regulate their gene expression. Such molecules include trans-acting small regulatory RNAs, antisense RNAs, and a variety of transcriptional attenuation mechanisms in the 5' untranslated region. Thus far, most regulatory RNA research has focused on Gram-negative bacteria, such as Escherichia coli and Salmonella. Hence, there is uncertainty about whether the resulting insights can be extrapolated directly to other bacteria, such as the Gram-positive soil bacterium Bacillus subtilis. A recent study identified 1,583 putative regulatory RNAs in B. subtilis, whose expression was assessed across 104 conditions. Here, we review the current understanding of RNA-based regulation in B. subtilis, and we categorize the newly identified putative regulatory RNAs on the basis of their conservation in other bacilli and the stability of their predicted secondary structures. Our present evaluation of the publicly available data indicates that RNA-mediated gene regulation in B. subtilis mostly involves elements at the 5' ends of mRNA molecules. These can include 5' secondary structure elements and metabolite-, tRNA-, or protein-binding sites. Importantly, sense-independent segments are identified as the most conserved and structured potential regulatory RNAs in B. subtilis. Altogether, the present survey provides many leads for the identification of new regulatory RNA functions in B. subtilis.

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.

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.

The reduction in small ribosomal subunit abundance in ethanol-stressed cells of Bacillus subtilis is mediated by a SigB-dependent antisense RNA

One of the best-characterized general stress responses in bacteria is the σB-mediated stress response of the Gram-positive soil bacterium Bacillus subtilis. The σB regulon contains approximately 200 protein-encoding genes and 136 putative regulatory RNAs. One of these σB-dependent RNAs, named S1136-S1134, was recently mapped as being transcribed from the S1136 promoter on the opposite strand of the essential rpsD gene, which encodes the ribosomal primary-binding protein S4. Accordingly, S1136-S1134 transcription results in an rpsD-overlapping antisense RNA (asRNA). Upon exposure of B. subtilis to ethanol, the S1136 promoter was found to be induced, while rpsD transcription was downregulated. By quantitative PCR, we show that the activation of transcription from the S1136 promoter is directly responsible for the downregulation of rpsD upon ethanol exposure. We also show that this downregulation of rpsD leads to a reduced level of the small (30S) ribosomal subunit upon ethanol stress. The activation of the S1136 promoter thus represents the first example of antisense transcription-mediated regulation in the general stress response of B. subtilis and implicates the reduction of ribosomal protein abundance as a new aspect in the σB-dependent stress response. We propose that the observed reduction in the level of the small ribosomal subunit, which contains the ribosome-decoding center, may protect B. subtilis cells against misreading and spurious translation of possibly toxic aberrant peptides under conditions of ethanol stress.

Ribosomal protein L10(L12)4 autoregulates expression of the Bacillus subtilis rplJL operon by a transcription attenuation mechanism

Ribosomal protein genes are often controlled by autoregulatory mechanisms in which a protein encoded in the operon can either bind to newly synthesized rRNA during rapid growth or to a similar target in its mRNA during poor growth conditions. The rplJL operon encodes the ribosomal L10(L12)₄ complex. In Escherichia coli L10(L12)₄ represses its translation by binding to the rplJL leader transcript. We identified three RNA structures in the Bacillus subtilis rplJL leader transcript that function as an anti-antiterminator, antiterminator or intrinsic terminator. Expression studies with transcriptional and translational fusions indicated that L10(L12)₄ represses rplJL expression at the transcriptional level. RNA binding studies demonstrated that L10(L12)₄ stabilizes the anti-antiterminator structure, while in vitro transcription results indicated that L10(L12)₄ promotes termination. Disruption of anti-antiterminator, antiterminator or terminator function by competitor oligonucleotides in vitro and by mutations in vivo demonstrated that each structure functions as predicted. Thus, rplJL expression is regulated by an autogenous transcription attenuation mechanism in which L10(L12)₄ binding to the anti-antiterminator structure promotes termination. We also found that translation of a leader peptide increases rplJL expression, presumably by inhibiting Rho-dependent termination. Thus, the rplJL operon of B. subtilis is regulated by transcription attenuation and antitermination mechanisms.

Extreme features of the Galdieria sulphuraria organellar genomes: a consequence of polyextremophily?

Nuclear genome sequencing from extremophilic eukaryotes has revealed clues about the mechanisms of adaptation to extreme environments, but the functional consequences of extremophily on organellar genomes are unknown. To address this issue, we assembled the mitochondrial and plastid genomes from a polyextremophilic red alga, Galdieria sulphuraria strain 074 W, and performed a comparative genomic analysis with other red algae and more broadly across eukaryotes. The mitogenome is highly reduced in size and genetic content and exhibits the highest guanine–cytosine skew of any known genome and the fastest substitution rate among all red algae. The plastid genome contains a large number of intergenic stem-loop structures but is otherwise rather typical in size, structure, and content in comparison with other red algae. We suggest that these unique genomic modifications result not only from the harsh conditions in which Galdieria lives but also from its unusual capability to grow heterotrophically, endolithically, and in the dark. These conditions place additional mutational pressures on the mitogenome due to the increased reliance on the mitochondrion for energy production, whereas the decreased reliance on photosynthesis and the presence of numerous stem-loop structures may shield the plastome from similar genomic stress.

The expanding view of RNA and DNA function

RNA and DNA are simple linear polymers consisting of only four major types of subunits, and yet these molecules carry out a remarkable diversity of functions in cells and in the laboratory. Each newly discovered function of natural or engineered nucleic acids enforces the view that prior assessments of nucleic acid function were far too narrow and that many more exciting findings are yet to come. This Perspective highlights just a few of the numerous discoveries over the past 20 years pertaining to nucleic acid function, focusing on those that have been of particular interest to chemical biologists. History suggests that there will continue to be many opportunities to engage chemical biologists in the discovery, creation, and manipulation of nucleic acid function in the years to come.

Discovery and validation of novel and distinct RNA regulators for ribosomal protein S15 in diverse bacterial phyla

Background

Autogenous cis-regulators of ribosomal protein synthesis play a critical role in maintaining the stoichiometry of ribosome components. Structured portions within an mRNA transcript typically interact with specific ribosomal proteins to prevent expression of the entire operon, thus balancing levels of ribosomal proteins across transcriptional units. Three distinct RNA structures from different bacterial phyla have demonstrated interactions with S15 to regulate gene expression; however, these RNAs are distributed across a small fraction of bacterial diversity.

Results

We used comparative genomics in combination with analysis of existing transcriptomic data to identify three novel putative RNA structures associated with the S15 coding region in microbial genomes. These structures are completely distinct from those previously published and encompass potential regulatory regions including ribosome-binding sites. To validate the biological relevance of our findings, we demonstrate that an example of the Alphaproteobacterial RNA from Rhizobium radiobacter specifically interacts with S15 in vitro, and allows in vivo regulation of gene expression in an E. coli reporter system. In addition, structural probing and nuclease protection assays confirm the predicted secondary structure and indicate nucleotides required for protein interaction.

Conclusions

This work illustrates the importance of integrating comparative genomic and transcriptomic approaches during de novo ncRNA identification and reveals a diversity of distinct natural RNA regulators that support analogous biological functions. Furthermore, this work indicates that many additional uncharacterized RNA regulators likely exist within bacterial genomes and that the plasticity of RNA structure allows unique, and likely independently derived, solutions to the same biological problem.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-657) contains supplementary material, which is available to authorized users.

Bacterial RNA motif in the 5' UTR of rpsF interacts with an S6:S18 complex

Approximately half the transcripts encoding ribosomal proteins in Escherichia coli include a structured RNA motif that interacts with a specific ribosomal protein to inhibit gene expression, thus allowing stoichiometric production of ribosome components. However, many of these RNA structures are not widely distributed across bacterial phyla. It is increasingly common for RNA motifs associated with ribosomal protein genes to be identified using comparative genomic methods, yet these are rarely experimentally validated. In this work, we characterize one such motif that precedes operons containing rpsF and rpsR, which encode ribosomal proteins S6 and S18. This RNA structure is widely distributed across many phyla of bacteria despite differences within the downstream operon, and examples are present in both E. coli and Bacillus subtilis. We demonstrate a direct interaction between an example of the RNA from B. subtilis and an S6:S18 complex using in vitro binding assays, verify our predicted secondary structure, and identify a putative protein-binding site. The proposed binding site bears a strong resemblance to the S18 binding site within the 16S rRNA, suggesting molecular mimicry. This interaction is a valuable addition to the canon of ribosomal protein mRNA interactions. This work shows how experimental verification translates computational results into concrete knowledge of biological systems.

Comparative genomics of metabolic capacities of regulons controlled by cis-regulatory RNA motifs in bacteria

Background

In silico comparative genomics approaches have been efficiently used for functional prediction and reconstruction of metabolic and regulatory networks. Riboswitches are metabolite-sensing structures often found in bacterial mRNA leaders controlling gene expression on transcriptional or translational levels.

An increasing number of riboswitches and other cis-regulatory RNAs have been recently classified into numerous RNA families in the Rfam database. High conservation of these RNA motifs provides a unique advantage for their genomic identification and comparative analysis.

Results

A comparative genomics approach implemented in the RegPredict tool was used for reconstruction and functional annotation of regulons controlled by RNAs from 43 Rfam families in diverse taxonomic groups of Bacteria. The inferred regulons include ~5200 cis-regulatory RNAs and more than 12000 target genes in 255 microbial genomes. All predicted RNA-regulated genes were classified into specific and overall functional categories. Analysis of taxonomic distribution of these categories allowed us to establish major functional preferences for each analyzed cis-regulatory RNA motif family. Overall, most RNA motif regulons showed predictable functional content in accordance with their experimentally established effector ligands. Our results suggest that some RNA motifs (including thiamin pyrophosphate and cobalamin riboswitches that control the cofactor metabolism) are widespread and likely originated from the last common ancestor of all bacteria. However, many more analyzed RNA motifs are restricted to a narrow taxonomic group of bacteria and likely represent more recent evolutionary innovations.

Conclusions

The reconstructed regulatory networks for major known RNA motifs substantially expand the existing knowledge of transcriptional regulation in bacteria. The inferred regulons can be used for genetic experiments, functional annotations of genes, metabolic reconstruction and evolutionary analysis. The obtained genome-wide collection of reference RNA motif regulons is available in the RegPrecise database (http://regprecise.lbl.gov/).