Identification of 22 candidate structured RNAs in bacteria using the CMfinder comparative genomics pipeline

A spermidine riboswitch class in bacteria exploits a close variant of an aptamer for the enzyme cofactor S-adenosylmethionine

Natural polyamines such as spermidine and spermine cations have characteristics that make them highly likely to be sensed by riboswitches, such as their general affinity to polyanionic RNA and their broad contributions to cell physiology. Despite previous claims that polyamine riboswitches exist, evidence of their biological functions has remained unconvincing. Here, we report that rare variants of bacterial S-adenosylmethionine-I (SAM-I) riboswitches reject SAM and have adapted to selectively sense spermidine. These spermidine-sensing riboswitch variants are associated with genes whose protein products are directly involved in the production of spermidine and other polyamines. Biochemical and genetic assays demonstrate that representatives of this riboswitch class robustly function as genetic "off" switches, wherein spermidine binding causes premature transcription termination to suppress the expression of polyamine biosynthetic genes. These findings confirm the existence of natural spermidine-sensing riboswitches in bacteria and expand the list of variant riboswitch classes that have adapted to bind different ligands.

Structure and ion-dependent folding of k-junctions

k-Junctions are elaborated forms of kink turns with an additional helix on the nonbulged strand, thus forming a three-way helical junction. Two were originally identified in the structures of Arabidopsis and Escherichia coli thiamine pyrophosphate (TPP) riboswitches, and another called DUF-3268 was tentatively identified from sequence information. In this work we show that the Arabidopsis and E. coli riboswitch k-junctions fold in response to the addition of magnesium or sodium ions, and that atomic mutations that should disrupt key hydrogen bonding interactions greatly impair folding. Using X-ray crystallography, we have determined the structure of the DUF-3268 RNA and thus confirmed that it is a k-junction. It also folds upon the addition of metal ions, though requiring a 40-fold lower concentration of either divalent or monovalent ions. The key difference between the DUF-3268 and riboswitch k-junctions is the lack of nucleotides inserted between G1b and A2b in the former. We show that this insertion is primarily responsible for the difference in folding properties. Finally, we show that the DUF-3268 can functionally substitute for the k-junction in the E. coli TPP riboswitch such that the chimera can bind the TPP ligand, although less avidly.

An O2-sensing diguanylate cyclase broadly affects the aerobic transcriptome in the phytopathogen Pectobacterium carotovorum

Pectobacterium carotovorum is an important plant pathogen responsible for the destruction of crops through bacterial soft rot, which is modulated by oxygen (O₂) concentration. A soluble globin coupled sensor protein, Pcc DgcO (also referred to as PccGCS) is one way through which P. carotovorum senses oxygen. DgcO contains a diguanylate cyclase output domain producing c-di-GMP. Synthesis of the bacterial second messenger c-di-GMP is increased upon oxygen binding to the sensory globin domain. This work seeks to understand regulation of function by DgcO at the transcript level. RNA sequencing and differential expression analysis revealed that the deletion of DgcO only affects transcript levels in cells grown under aerobic conditions. Differential expression analysis showed that DgcO deletion alters transcript levels for metal transporters. These results, followed by inductively coupled plasma-mass spectrometry showing decreased concentrations of six biologically relevant metals upon DgcO deletion, provide evidence that a globin coupled sensor can affect cellular metal content. These findings improve the understanding of the transcript level control of O₂-dependent phenotypes in an important phytopathogen and establish a basis for further studies on c-di-GMP-dependent functions in P. carotovorum.

The structural features of the ligand-free moaA riboswitch and its ion-dependent folding

Riboswitches are structural elements of mRNA involved in the regulation of gene expression by responding to specific cellular metabolites. To fulfil their regulatory function, riboswitches prefold into an active state, the so-called binding competent form, that guarantees metabolite binding and allows a consecutive refolding of the RNA. Here, we describe the folding pathway to the binding competent form as well as the ligand free structure of the moaA riboswitch of E. coli. This RNA proposedly responds to the molybdenum cofactor (Moco), a highly oxygen-sensitive metabolite, essential in the carbon and sulfur cycles of eukaryotes. K⁺- and Mg²⁺-dependent footprinting assays and spectroscopic investigations show a high degree of structure formation of this RNA already at very low ion-concentrations. Mg²⁺ facilitates additionally a general compaction of the riboswitch towards its proposed active structure. We show that this fold agrees with the earlier suggested secondary structure which included also a long-range tetraloop/tetraloop-receptor like interaction. Metal ion cleavage assays revealed specific Mg²⁺-binding pockets within the moaA riboswitch. These Mg²⁺ binding pockets are good indicators for the potential Moco binding site, since in riboswitches, Mg²⁺ was shown to be necessary to bind phosphate-carrying metabolites. The importance of the phosphate and of other functional groups of Moco is highlighted by binding assays with tetrahydrobiopterin, the reduced and oxygen-sensitive core moiety of Moco. We demonstrate that the general molecular shape of pterin by its own is insufficient for the recognition by the riboswitch.

Transcription start site mapping and small RNA profiling of Leptospira biflexa serovar Patoc

Leptospirosis is an emerging zoonotic disease caused by bacterial species of the genus Leptospira. However, the regulatory mechanisms and pathways underlying the adaptation of pathogenic and non-pathogenic Leptospira spp. in different environmental conditions remain elusive. Leptospira biflexa is a non-pathogenic species of Leptospira that lives exclusively in a natural environment. It is an ideal model not only for exploring molecular mechanisms underlying the environmental survival of Leptospira species but also for identifying virulence factors unique to Leptospira's pathogenic species. In this study, we aim to establish the transcription start site (TSS) landscape and the small RNA (sRNA) profile of L. biflexa serovar Patoc grown to exponential and stationary phases via differential RNA-seq (dRNA-seq) and small RNA-seq (sRNA-seq) analyses, respectively. Our dRNA-seq analysis uncovered a total of 2726 TSSs, which are also used to identify other elements, e.g., promoter and untranslated regions (UTRs). Besides, our sRNA-seq analysis revealed a total of 603 sRNA candidates, comprising 16 promoter-associated sRNAs, 184 5'UTR-derived sRNAs, 230 true intergenic sRNAs, 136 5'UTR-antisense sRNAs, and 130 open reading frame (ORF)-antisense sRNAs. In summary, these findings reflect the transcriptional complexity of L. biflexa serovar Patoc under different growth conditions and help to facilitate our understanding of regulatory networks in L. biflexa. To the best of our knowledge, this is the first study reporting the TSS landscape of L. biflexa. The TSS and sRNA landscapes of L. biflexa can also be compared with its pathogenic counterparts, e.g., L. borgpetersenii and L. interrogans, to identify features contributing to their environmental survival and virulence.

Discovering riboswitches: the past and the future

Riboswitches are structured noncoding RNA domains used by many bacteria to monitor the concentrations of target ligands and regulate gene expression accordingly. In the past 20 years over 55 distinct classes of natural riboswitches have been discovered that selectively sense small molecules or elemental ions, and thousands more are predicted to exist. Evidence suggests that some riboswitches might be direct descendants of the RNA-based sensors and switches that were likely present in ancient organisms before the evolutionary emergence of proteins. We provide an overview of the current state of riboswitch research, focusing primarily on the discovery of riboswitches, and speculate on the major challenges facing researchers in the field.

A variant of guanidine-IV riboswitches exhibits evidence of a distinct ligand specificity

Riboswitches are regulatory RNAs that specifically bind a small molecule or ion. Like metabolite-binding proteins, riboswitches can evolve new ligand specificities, and some examples of this phenomenon have been validated. As part of work based on comparative genomics to discover novel riboswitches, we encountered a candidate riboswitch with striking similarities to the recently identified guanidine-IV riboswitch. This candidate riboswitch, the Gd4v motif, is predicted in four distinct bacterial phyla, thus almost as widespread as the guanidine-IV riboswitch. Bioinformatic and experimental analysis suggest that the Gd4v motif is a riboswitch that binds a ligand other than guanidine. It is found associated with gene classes that differ from genes regulated by confirmed guanidine riboswitches. In inline-probing assays, we showed that free guanidine binds only weakly to one of the tested sequences of the variant. Further tested compounds did not show binding, attenuation of transcription termination, or activation of a genetic reporter construct. We characterized an N-acetyltransferase frequently associated with the Gd4v motif and compared its substrate preference to an N-acetyltransferase that occurs under control of guanidine-IV riboswitches. The substrates of this Gd4v-motif-associated enzyme did not show activity for Gd4v RNA binding or transcription termination. Hence, the ligand of the candidate riboswitch motif remains unidentified. The variant RNA motif is predominantly found in gut metagenome sequences, hinting at a ligand that is highly relevant in this environment. This finding is a first step to determining the identity of this unknown ligand, and understanding how guanidine-IV-riboswitch-like structures can evolve to bind different ligands.

Structure-based insights into recognition and regulation of SAM-sensing riboswitches

Riboswitches are highly conserved RNA elements that located in the 5'-UTR of mRNAs, which undergo real-time structure conformational change to achieve the regulation of downstream gene expression by sensing their cognate ligands. S-adenosylmethionine (SAM) is a ubiquitous methyl donor for transmethylation reactions in all living organisms. SAM riboswitch is one of the most abundant riboswitches that bind to SAM with high affinity and selectivity, serving as regulatory modules in multiple metabolic pathways. To date, seven SAM-specific riboswitch classes that belong to four families, one SAM/SAH riboswitch and one SAH riboswitch have been identified. Each SAM riboswitch family has a well-organized tertiary core scaffold to support their unique ligand-specific binding pocket. In this review, we summarize the current research progress on the distribution, structure, ligand recognition and gene regulation mechanism of these SAM-related riboswitch families, and further discuss their evolutionary prospects and potential applications.

Bioinformatics and Genomic Analyses of the Suitability of Eight Riboswitches for Antibacterial Drug Targets

Antibiotic resistance (AR) is an acute problem that results in prolonged and debilitating illnesses. AR mortality worldwide is growing and causes a pressing need to research novel mechanisms of action and untested target molecules. This article presents in silico analyses of eight bacterial riboswitches for their suitability for antibacterial drug targets. Most bacterial riboswitches are located in the 5′-untranslated region of messenger RNAs, act as allosteric cis-acting gene control elements, and have not been found in humans before. Sensing metabolites, the riboswitches regulate the synthesis of vital cellular metabolites in various pathogenic bacteria. The analyses performed in this article represent a complete and informative genome-wide bioinformatics analysis of the adequacy of eight riboswitches as antibacterial drug targets in different pathogenic bacteria based on four criteria. Due to the ability of the riboswitch to control biosynthetic pathways and transport proteins of essential metabolites and the presence/absence of alternative biosynthetic pathways, we classified them into four groups based on their suitability for use as antibacterial drug targets guided by our in silico analyses. We concluded that some of them are promising targets for antibacterial drug discovery, such as the PreQ1, MoCo RNA, cyclic-di-GMP I, and cyclic-di-GMP II riboswitches.

The cotranscriptional folding landscape for two cyclic di-nucleotide-sensing riboswitches with highly homologous aptamer domains acting either as ON- or OFF-switches

Riboswitches are gene regulatory elements located in untranslated mRNA regions. They bind inducer molecules with high affinity and specificity. Cyclic-di-nucleotide-sensing riboswitches are major regulators of genes for the environment, membranes and motility (GEMM) of bacteria. Up to now, structural probing assays or crystal structures have provided insight into the interaction between cyclic-di-nucleotides and their corresponding riboswitches. ITC analysis, NMR analysis and computational modeling allowed us to gain a detailed understanding of the gene regulation mechanisms for the Cd1 (Clostridium difficile) and for the pilM (Geobacter metallireducens) riboswitches and their respective di-nucleotides c-di-GMP and c-GAMP. Binding capability showed a 25 nucleotide (nt) long window for pilM and a 61 nt window for Cd1. Within this window, binding affinities ranged from 35 μM to 0.25 μM spanning two orders of magnitude for Cd1 and pilM showing a strong dependence on competing riboswitch folds. Experimental results were incorporated into a Markov simulation to further our understanding of the transcriptional folding pathways of riboswitches. Our model showed the ability to predict riboswitch gene regulation and its dependence on transcription speed, pausing and ligand concentration.

Variants of the guanine riboswitch class exhibit altered ligand specificities for xanthine, guanine, or 2'-deoxyguanosine

The aptamer portions of previously reported riboswitch classes that sense guanine, adenine, or 2′-deoxyguanosine are formed by a highly similar three-stem junction with distinct nucleotide sequences in the regions joining the stems. The nucleotides in these joining regions form the major features of the selective ligand-binding pocket for each aptamer. Previously, we reported the existence of additional, rare variants of the predominant guanine-sensing riboswitch class that carry nucleotide differences in the ligand-binding pocket, suggesting that these RNAs have further diversified their structures and functions. Herein, we report the discovery and analysis of three naturally occurring variants of guanine riboswitches that are narrowly distributed across Firmicutes. These RNAs were identified using comparative sequence analysis methods, which also revealed that some of the gene associations for these variants are atypical for guanine riboswitches or their previously known natural variants. Binding assays demonstrate that the newfound variant riboswitch representatives recognize xanthine, guanine, or 2′-deoxyguanosine, with the guanine class exhibiting greater discrimination against related purines than the more common guanine riboswitch class reported previously. These three additional variant classes, together with the four previously discovered riboswitch classes that employ the same three-stem junction architecture, reveal how a simple structural framework can be diversified to expand the range of purine-based ligands sensed by RNA.

Combining Coarse-Grained Simulations and Single Molecule Analysis Reveals a Three-State Folding Model of the Guanidine-II Riboswitch

Riboswitch RNAs regulate gene expression by conformational changes induced by environmental conditions and specific ligand binding. The guanidine-II riboswitch is proposed to bind the small molecule guanidinium and to subsequently form a kissing loop interaction between the P1 and P2 hairpins. While an interaction was shown for isolated hairpins in crystallization and electron paramagnetic resonance experiments, an intrastrand kissing loop formation has not been demonstrated. Here, we report the first evidence of this interaction in cis in a ligand and Mg²⁺ dependent manner. Using single-molecule FRET spectroscopy and detailed structural information from coarse-grained simulations, we observe and characterize three interconvertible states representing an open and kissing loop conformation as well as a novel Mg²⁺ dependent state for the guanidine-II riboswitch from E. coli. The results further substantiate the proposed switching mechanism and provide detailed insight into the regulation mechanism for the guanidine-II riboswitch class. Combining single molecule experiments and coarse-grained simulations therefore provides a promising perspective in resolving the conformational changes induced by environmental conditions and to yield molecular insights into RNA regulation.

Identifying candidate structured RNAs in CRISPR operons

Noncoding RNAs with secondary structures play important roles in CRISPR-Cas systems. Many of these structures likely remain undiscovered. We used a large-scale comparative genomics approach to predict 156 novel candidate structured RNAs from 36,111 CRISPR-Cas systems. A number of these were found to overlap with coding genes, including palindromic candidates that overlapped with a variety of Cas genes in type I and III systems. Among these 156 candidates, we identified 46 new models of CRISPR direct repeats and 1 tracrRNA. This tracrRNA model occasionally overlapped with predicted cas9 coding regions, emphasizing the importance of expanding our search windows for novel structure RNAs in coding regions. We also demonstrated that the antirepeat sequence in this tracrRNA model can be used to accurately assign thousands of predicted CRISPR arrays to type II-C systems. This study highlights the importance of unbiased identification of candidate structured RNAs across CRISPR-Cas systems.

Cryo-EM advances in RNA structure determination

Cryo-electron microscopy (cryo-EM) has emerged as an unprecedented tool to resolve protein structures at atomic resolution. Structural insights of biological samples not accessible by conventional X-ray crystallography and NMR can be explored with cryo-EM because measurements are carried out under near-native crystal-free conditions, and large protein complexes with conformational and compositional heterogeneity are readily resolved. RNA has remained underexplored in cryo-EM, despite its essential role in various biological processes. This review highlights current challenges and recent progress in using cryo-EM single-particle analysis to determine protein-free RNA structures, enabled by improvement in sample preparation and integration of multiple structural and biochemical methods.

The application of riboswitch sequencing for human gut bacterial classification and identification

Bacteria, especially gut bacteria play important roles in human health and diseases. The classification of many bacterial genera by the 16S ribosomal RNA (rRNA) has failed due to its low inter-species resolution. Given the wide distribution of riboswitches in bacteria, they may help 16S rRNA differentiate closely related species. We found that among 28 groups of species that could not be distinguished by 16S rRNA, eight of them could be separated by the TPP riboswitch and other riboswitches. Moreover, the species in the 16S rRNA database and these riboswitch databases overlap, therefore, using riboswitch databases can help 16S rRNA better identify species. In addition, we used Klenow DNA polymerase and a pair of short primers to facilitate the library construction of TPP riboswitches for sequencing. The sequencing results showed that the TPP riboswitch could detect the major phyla similar to those detected by 16S rRNA. Therefore, the TPP riboswitch and other riboswitch classes could potentially be applied to gut bacteria classification.

Conserved Motifs and Domains in Members of Pospiviroidae

In 1985, Keese and Symons proposed a hypothesis on the sequence and secondary structure of viroids from the family Pospiviroidae: their secondary structure can be subdivided into five structural and functional domains and “viroids have evolved by rearrangement of domains between different viroids infecting the same cell and subsequent mutations within each domain”; this article is one of the most cited in the field of viroids. Employing the pairwise alignment method used by Keese and Symons and in addition to more recent methods, we tried to reproduce the original results and extent them to further members of Pospiviroidae which were unknown in 1985. Indeed, individual members of Pospiviroidae consist of a patchwork of sequence fragments from the family but the lengths of fragments do not point to consistent points of rearrangement, which is in conflict with the original hypothesis of fixed domain borders.

Characterization of Structure and Dynamics of the Guanidine-II Riboswitch from Escherichia coli by NMR Spectroscopy and Small-Angle X-ray Scattering (SAXS)

Riboswitches are regulatory RNA elements that undergo functionally important allosteric conformational switching upon binding of specific ligands. The here investigated guanidine‐II riboswitch binds the small cation, guanidinium, and forms a kissing loop‐loop interaction between its P1 and P2 hairpins. We investigated the structural changes to support previous studies regarding the binding mechanism. Using NMR spectroscopy, we confirmed the structure as observed in crystal structures and we characterized the kissing loop interaction upon addition of Mg²⁺ and ligand for the riboswitch aptamer from Escherichia coli. We further investigated closely related mutant constructs providing further insight into functional differences between the two (different) hairpins P1 and P2. Formation of intermolecular interactions were probed by small‐angle X‐ray scattering (SAXS) and NMR DOSY data. All data are consistent and show the formation of oligomeric states of the riboswitch induced by Mg²⁺ and ligand binding.

Small RNA Profiling in Mycobacterium Provides Insights Into Stress Adaptability

Mycobacteria encounter a number of environmental changes during infection and respond using different mechanisms. Small RNA (sRNA) is a post-transcriptionally regulatory system for gene functions and has been investigated in many other bacteria. This study used Mycobacterium tuberculosis and Mycobacterium bovis Bacillus Calmette-Guérin (BCG) infection models and sequenced whole bacterial RNAs before and after host cell infection. A comparison of differentially expressed sRNAs using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) and target prediction was carried out. Six pathogenically relevant stress conditions, growth rate, and morphology were used to screen and identify sRNAs. From these data, a subset of sRNAs was differentially expressed in multiple infection groups and stress conditions. Many were found associated with lipid metabolism. Among them, ncBCG427 was significantly downregulated when BCG entered into macrophages and was associated with increased biofilm formation. The reduction of virulence possibility depends on regulating lipid metabolism.

Small RNAs Asserting Big Roles in Mycobacteria

Tuberculosis (TB) is an infectious disease caused by Mycobacterium tuberculosis (Mtb), with 10.4 million new cases per year reported in the human population. Recent studies on the Mtb transcriptome have revealed the abundance of noncoding RNAs expressed at various phases of mycobacteria growth, in culture, in infected mammalian cells, and in patients. Among these noncoding RNAs are both small RNAs (sRNAs) between 50 and 350 nts in length and smaller RNAs (sncRNA) < 50 nts. In this review, we provide an up-to-date synopsis of the identification, designation, and function of these Mtb-encoded sRNAs and sncRNAs. The methodological advances including RNA sequencing strategies, small RNA antagonists, and locked nucleic acid sequence-specific RNA probes advancing the studies on these small RNA are described. Initial insights into the regulation of the small RNA expression and putative processing enzymes required for their synthesis and function are discussed. There are many open questions remaining about the biological and pathogenic roles of these small non-coding RNAs, and potential research directions needed to define the role of these mycobacterial noncoding RNAs are summarized.

Molecular Dynamics Simulations of the Aptamer Domain of Guanidinium Ion Binding Riboswitch ykkC-III: Structural Insights into the Discrimination of Cognate and Alternate Ligands

Guanidinium ion is a toxic cellular metabolite. The ykkC-III riboswitch, an mRNA stretch, regulates the gene expression by undergoing a conformational change in response to the binding of a free guanidinium ion and thereby plays a potentially important role in alleviating guanidinium toxicity in cells. An experimental crystal structure of the guanidinium-bound aptamer domain of the riboswitch from Thermobifida Fusca revealed the overall RNA architecture and mapped the specific noncovalent interactions that stabilize the ligand within the binding pocket aptamer. However, details of how the aptamer domain discriminates the cognate ligand from its closest structurally analogous physiological metabolites (arginine and urea), and how the binding of cognate ligand arrays information from the aptamer domain to the expression platform for regulating the gene expression, are not well understood. To fill this void, we perform a cumulative of 2 μs all-atom explicit-solvent molecular dynamics (MD) simulations on the full aptamer domain, augmented with quantum-chemical calculations on the ligand-binding pocket, to compare the structural and dynamical details of the guanidinium-bound state with the arginine or urea bound states, as well as the unbound (open) state. Analysis of the ligand-binding pocket reveals that due to unfavorable interactions with the binding-pocket residues, urea cannot bind the aptamer domain and thereby cannot alter the gene expression. Although interaction of the guanidyl moiety of arginine within the binding pocket is either comparable or stronger than the guanidinium ion, additional non-native hydrogen-bonding networks, as well as differences in the dynamical details of the arginine-bound state, explain why arginine cannot transmit the information from the aptamer domain to the expression platform. Based on our simulations, we propose a mechanism of how the aptamer domain communicates with the expression platform. Overall, our work provides interesting insights into the ligand recognition by a specific class of riboswitches and may hopefully inspire future studies to further understand the gene regulation by riboswitches.

Getting to the bottom of lncRNA mechanism: structure-function relationships

While long non-coding RNAs are known to play key roles in disease and development, relatively few structural studies have been performed for this important class of RNAs. Here, we review functional studies of long non-coding RNAs and expose the need for high-resolution 3-D structural studies, discussing the roles of long non-coding RNAs in the cell and how structure–function relationships might be used to elucidate further understanding. We then describe structural studies of other classes of RNAs using chemical probing, nuclear magnetic resonance, small-angle X-ray scattering, X-ray crystallography, and cryogenic electron microscopy (cryo-EM). Next, we review early structural studies of long non-coding RNAs to date and describe the way forward for the structural biology of long non-coding RNAs in terms of cryo-EM.

DIMPL: a bioinformatics pipeline for the discovery of structured noncoding RNA motifs in bacteria

SUMMARY

Recent efforts to identify novel bacterial structured noncoding RNA (ncRNA) motifs through searching long, GC-rich intergenic regions (IGRs) have revealed several new classes, including the recently validated HMP-PP riboswitch. The DIMPL (Discovery of Intergenic Motifs PipeLine) discovery pipeline described herein enables rapid extraction and selection of bacterial IGRs that are enriched for structured ncRNAs. Moreover, DIMPL automates the subsequent computational steps necessary for their functional identification.

AVAILABILITY AND IMPLEMENTATION

The DIMPL pipeline is freely available as a Docker image with an accompanying set of Jupyter notebooks. Full instructions for download and use are available at https://github.com/breakerlab/dimpl.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Siblings or doppelgängers? Deciphering the evolution of structured cis-regulatory RNAs beyond homology

Structured cis-regulatory RNAs have evolved across all domains of life, highlighting the utility and plasticity of RNA as a regulatory molecule. Homologous RNA sequences and structures often have similar functions, but homology may also be deceiving. The challenges that derive from trying to assign function to structure and vice versa are not trivial. Bacterial riboswitches, viral and eukaryotic IRESes, CITEs, and 3′ UTR elements employ an array of mechanisms to exert their effects. Bioinformatic searches coupled with biochemical and functional validation have elucidated some shared and many unique ways cis-regulators are employed in mRNA transcripts. As cis-regulatory RNAs are resolved in greater detail, it is increasingly apparent that shared homology can mask the full spectrum of mRNA cis-regulator functional diversity. Furthermore, similar functions may be obscured by lack of obvious sequence similarity. Thus looking beyond homology is crucial for furthering our understanding of RNA-based regulation.

Guanidine-II aptamer conformations and ligand binding modes through the lens of molecular simulation

Regulation of gene expression via riboswitches is a widespread mechanism in bacteria. Here, we investigate ligand binding of a member of the guanidine sensing riboswitch family, the guanidine-II riboswitch (Gd-II). It consists of two stem–loops forming a dimer upon ligand binding. Using extensive molecular dynamics simulations we have identified conformational states corresponding to ligand-bound and unbound states in a monomeric stem–loop of Gd-II and studied the selectivity of this binding. To characterize these states and ligand-dependent conformational changes we applied a combination of dimensionality reduction, clustering, and feature selection methods. In absence of a ligand, the shape of the binding pocket alternates between the conformation observed in presence of guanidinium and a collapsed conformation, which is associated with a deformation of the dimerization interface. Furthermore, the structural features responsible for the ability to discriminate against closely related analogs of guanidine are resolved. Based on these insights, we propose a mechanism that couples ligand binding to aptamer dimerization in the Gd-II system, demonstrating the value of computational methods in the field of nucleic acids research.

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.

Comprehensive discovery of novel structured noncoding RNAs in 26 bacterial genomes

Comparative sequence analysis methods are highly effective for uncovering novel classes of structured noncoding RNAs (ncRNAs) from bacterial genomic DNA sequence datasets. Previously, we developed a computational pipeline to more comprehensively identify structured ncRNA representatives from individual bacterial genomes. This search process exploits the fact that genomic regions serving as templates for the transcription of structured RNAs tend to be present in longer than average noncoding 'intergenic regions' (IGRs) that are enriched in G and C nucleotides compared to the remainder of the genome. In the present study, we apply this computational pipeline to identify structured ncRNA candidates from 26 diverse bacterial species. Numerous novel structured ncRNA motifs were discovered, including several riboswitch candidates, one whose ligand has been identified and others that have yet to be experimentally validated. Our findings support recent predictions that hundreds of novel ribo-switch classes and other ncRNAs remain undiscovered among the limited number of bacterial species whose genomes have been completely sequenced.

Identification of 11 candidate structured noncoding RNA motifs in humans by comparative genomics

Background

Only 1.5% of the human genome encodes proteins, while large part of the remaining encodes noncoding RNAs (ncRNA). Many ncRNAs form structures and perform many important functions. Accurately identifying structured ncRNAs in the human genome and discovering their biological functions remain a major challenge.

Results

Here, we have established a pipeline (CM-line) with the following features for analyzing the large genomes of humans and other animals. First, we selected species with larger genetic distances to facilitate the discovery of covariations and compatible mutations. Second, we used CMfinder, which can generate useful alignments even with low sequence conservation. Third, we removed repetitive sequences and known structured ncRNAs to reduce the workload of CMfinder. Fourth, we used Infernal to find more representatives and refine the structure. We reported 11 classes of structured ncRNA candidates with significant covariations in humans. Functional analysis showed that these ncRNAs may have variable functions. Some may regulate circadian clock genes through poly (A) signals (PAS); some may regulate the elongation factor (EEF1A) and the T-cell receptor signaling pathway by cooperating with RNA binding proteins.

Conclusions

By searching for important features of RNA structure from large genomes, the CM-line has revealed the existence of a variety of novel structured ncRNAs. Functional analysis suggests that some newly discovered ncRNA motifs may have biological functions. The pipeline we have established for the discovery of structured ncRNAs and the identification of their functions can also be applied to analyze other large genomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07474-9.

Widespread bacterial utilization of guanidine as nitrogen source

Guanidine is sensed by at least four different classes of riboswitches that are widespread in bacteria. However, only very few insights into physiological roles of guanidine exist. Genes predominantly regulated by guanidine riboswitches are Gdx transporters exporting the compound from the bacterial cell. In addition, urea/guanidine carboxylases and associated hydrolases and ABC transporters are often found combined in guanidine-inducible operons. We noted that the associated ABC transporters are configured to function as importers, challenging the current view that riboswitches solely control the detoxification of guanidine in bacteria. We demonstrate that the carboxylase pathway enables utilization of guanidine as sole nitrogen source. We isolated three enterobacteria (Raoultella terrigena, Klebsiella michiganensis, and Erwinia rhapontici) that utilize guanidine efficiently as N-source. Proteome analyses show that the expression of a carboxylase, associated hydrolases and transport genes is strongly induced by guanidine. Finding two urea/guanidine carboxylase enzymes in E. rhapontici, we demonstrate that the riboswitch-controlled carboxylase displays specificity toward guanidine, whereas the other enzyme prefers urea. We characterize the distribution of riboswitch-associated carboxylases and Gdx exporters in bacterial habitats by analyzing available metagenome data. The findings represent a paradigm shift from riboswitch-controlled detoxification of guanidine to the uptake and assimilation of this enigmatic nitrogen-rich compound.

A computational approach to identify CRISPR-Cas loci in the complete genomes of the lichen-associated Burkholderia sp. PAMC28687 and PAMC26561

The genus Burkholderia and its strains PAMC28687 and PAMC26561 are lichen-associated bacteria isolated from the Antarctic region. Our study is the first to provide the genome sequence of the Burkholderia sp. PAMC26561 strain. The genus Burkholderia includes bacteria that are pathogenic to plants, animals, and humans. Computational analysis of complete genomes of strains from the uncategorized Burkholderia group was performed using the NCBI databank and PATRIC (for genome sequence information) and CRISPRCasFinder (online and offline versions) software in order to predict the CRISPR loci and Cas genes. The RNAfold Webserver online software was used to predict RNA secondary structures. Our study showed that strain MSMB0852 (plasmid) possesses CRISPR-Cas system Class 2, and two lichen-associated strains, PAMC28687 (chromosome I) and PAMC26561 (chromosome I), possess CRISPR-Cas system Class 1. Additionally, only the two lichen-associated strains possess a variety of Cas genes.

Discovery and characterization of a fourth class of guanidine riboswitches

Riboswitches are RNAs that specifically sense a small molecule and regulate genes accordingly. The recent discovery of guanidine-binding riboswitches revealed the biological significance of this compound, and uncovered genes related to its biology. For example, certain sugE genes encode guanidine exporters and are activated by the riboswitches to reduce toxic levels of guanidine in the cell. In order to study guanidine biology and riboswitches, we applied a bioinformatics strategy for discovering additional guanidine riboswitches by searching for new candidate motifs associated with sugE genes. Based on in vitro and in vivo experiments, we determined that one of our six best candidates is a new structural class of guanidine riboswitches. The expression of a genetic reporter was induced 80-fold in response to addition of 5 mM guanidine in Staphylococcus aureus. This new class, called the guanidine-IV riboswitch, reveals additional guanidine-associated protein domains that are extremely rarely or never associated with previously established guanidine riboswitches. Among these protein domains are two transporter families that are structurally distinct from SugE, and could represent novel types of guanidine exporters. These results establish a new metabolite-binding RNA, further validate a bioinformatics method for finding riboswitches and suggest substrate specificities for as-yet uncharacterized transporter proteins.

Structural diversity and phylogenetic distribution of valyl tRNA-like structures in viruses

Viruses commonly use specifically folded RNA elements that interact with both host and viral proteins to perform functions important for diverse viral processes. Examples are found at the 3' termini of certain positive-sense ssRNA virus genomes where they partially mimic tRNAs, including being aminoacylated by host cell enzymes. Valine-accepting tRNA-like structures (TLS^Val) are an example that share some clear homology with canonical tRNAs but have several important structural differences. Although many examples of TLS^Val have been identified, we lacked a full understanding of their structural diversity and phylogenetic distribution. To address this, we undertook an in-depth bioinformatic and biochemical investigation of these RNAs, guided by recent high-resolution structures of a TLS^Val We cataloged many new examples in plant-infecting viruses but also in unrelated insect-specific viruses. Using biochemical and structural approaches, we verified the secondary structure of representative TLS^Val substrates and tested their ability to be valylated, confirming previous observations of structural heterogeneity within this class. In a few cases, large stem-loop structures are inserted within variable regions located in an area of the TLS distal to known host cell factor binding sites. In addition, we identified one virus whose TLS has switched its anticodon away from valine, causing a loss of valylation activity; the implications of this remain unclear. These results refine our understanding of the structural and functional mechanistic details of tRNA mimicry and how this may be used in viral infection.

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.

Do the P1 and P2 hairpins of the Guanidine-II riboswitch interact?

Riboswitches regulate genes by adopting different structures in responds to metabolite binding. The guanidine-II riboswitch is the smallest representative of the ykkC class with the mechanism of its function being centred on the idea that its two stem loops P1 and P2 form a kissing hairpin interaction upon binding of guanidinium (Gdm⁺). This mechanism is based on in-line probing experiments with the full-length riboswitch and crystal structures of the truncated stem loops P1 and P2. However, the crystal structures reveal only the formation of the homodimers P1 | P1 and P2 | P2 but not of the proposed heterodimer P1 | P2. Here, site-directed spin labeling (SDSL) in combination with Pulsed Electron–Electron Double Resonance (PELDOR or DEER) is used to study their structures in solution and how they change upon binding of Gdm⁺. It is found that both hairpins adopt different structures in solution and that binding of Gdm⁺ does indeed lead to the formation of the heterodimer but alongside the homodimers in a statistical 1:2:1 fashion. These results do thus support the proposed switching mechanism.

Advances in RNA 3D Structure Modeling Using Experimental Data

RNA is a unique bio-macromolecule that can both record genetic information and perform biological functions in a variety of molecular processes, including transcription, splicing, translation, and even regulating protein function. RNAs adopt specific three-dimensional conformations to enable their functions. Experimental determination of high-resolution RNA structures using x-ray crystallography is both laborious and demands expertise, thus, hindering our comprehension of RNA structural biology. The computational modeling of RNA structure was a milestone in the birth of bioinformatics. Although computational modeling has been greatly improved over the last decade showing many successful cases, the accuracy of such computational modeling is not only length-dependent but also varies according to the complexity of the structure. To increase credibility, various experimental data were integrated into computational modeling. In this review, we summarize the experiments that can be integrated into RNA structure modeling as well as the computational methods based on these experimental data. We also demonstrate how computational modeling can help the experimental determination of RNA structure. We highlight the recent advances in computational modeling which can offer reliable structure models using high-throughput experimental data.

Former orphan riboswitches reveal unexplored areas of bacterial metabolism, signaling, and gene control processes

Comparative sequence analyses have been used to discover numerous classes of structured noncoding RNAs, some of which are riboswitches that specifically recognize small-molecule or elemental ion ligands and influence expression of adjacent downstream genes. Determining the correct identity of the ligand for a riboswitch candidate typically is aided by an understanding of the genes under its regulatory control. Riboswitches whose ligands were straightforward to identify have largely been associated with well-characterized metabolic pathways, such as coenzyme or amino acid biosynthesis. Riboswitch candidates whose ligands resist identification, collectively known as orphan riboswitches, are often associated with genes coding for proteins of unknown function, or genes for various proteins with no established link to one another. The cognate ligands for 16 former orphan riboswitch motifs have been identified to date. The successful pursuit of the ligands for these classes has provided insight into areas of biology that are not yet fully explored, such as ion homeostasis, signaling networks, and other previously underappreciated biochemical or physiological processes. Herein we discuss the strategies and methods used to match ligands with orphan riboswitch classes, and overview the lessons learned to inform and motivate ongoing efforts to identify ligands for the many remaining candidates.

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.

Rational engineering of transcriptional riboswitches leads to enhanced metabolite levels in Bacillus subtilis

Many metabolic pathways in bacteria are regulated by metabolite sensing riboswitches that exert their control at the level of transcription employing a termination-antitermination mechanism. These riboswitches represent engineering targets to modulate expression of genes and operons relevant for the biotechnological production of commercially relevant compounds. We show that removal of the transcriptional riboswitches that control purine biosynthesis and riboflavin biosynthesis in Bacillus subtilis leads to auxotrophic strains. As an alternative, we report a rational approach for engineering transcriptional riboswitches independently from the availability of structural data. This approach consists in the identification and deletion of a key nucleotide sequence exclusively involved in transcription termination without affecting formation of other secondary and tertiary structures, which can be involved in other functions. To demonstrate the efficacy of our approach, we tested it with regard to deregulation of the purine and the riboflavin biosynthetic pathways in B. subtilis. Following validation of the engineered transcriptional riboswitches using specialized reporter strains, our approach was implemented into a B. subtilis wild-type strain employing CRISPR-Cas9 genome editing. The resulting purine and riboflavin production strains were characterized at the level of gene expression, metabolite synthesis and growth, and a substantial enhancement was measured at each level. Moreover, applying our approach to deregulate the purine pathway of an industrial riboflavin overproducing strain with impaired growth led to an increase in biomass by 53%, which resulted in an enhanced total production of riboflavin in the culture.

Nucleobase mutants of a bacterial preQ1-II riboswitch that uncouple metabolite sensing from gene regulation

Riboswitches are a class of nonprotein-coding RNAs that directly sense cellular metabolites to regulate gene expression. They are model systems for analyzing RNA-ligand interactions and are established targets for antibacterial agents. Many studies have analyzed the ligand-binding properties of riboswitches, but this work has outpaced our understanding of the underlying chemical pathways that govern riboswitch-controlled gene expression. To address this knowledge gap, we prepared 15 mutants of the preQ1-II riboswitch-a structurally and biochemically well-characterized HLout pseudoknot that recognizes the metabolite prequeuosine1 (preQ1). The mutants span the preQ1-binding pocket through the adjoining Shine-Dalgarno sequence (SDS) and include A-minor motifs, pseudoknot-insertion helix P4, U·A-U base triples, and canonical G-C pairs in the anti-SDS. As predicted-and confirmed by in vitro isothermal titration calorimetry measurements-specific mutations ablated preQ1 binding, but most aberrant binding effects were corrected by compensatory mutations. In contrast, functional analysis in live bacteria using a riboswitch-controlled GFPuv-reporter assay revealed that each mutant had a deleterious effect on gene regulation, even when compensatory changes were included. Our results indicate that effector binding can be uncoupled from gene regulation. We attribute loss of function to defects in a chemical interaction network that links effector binding to distal regions of the fold that support the gene-off RNA conformation. Our findings differentiate effector binding from biological function, which has ramifications for riboswitch characterization. Our results are considered in the context of synthetic ligands and drugs that bind tightly to riboswitches without eliciting a biological response.

Transposon mutagenesis in Mycobacterium kansasii links a small RNA gene to colony morphology and biofilm formation and identifies 9,885 intragenic insertions that do not compromise colony outgrowth

Mycobacterium kansasii (Mk) is a resilient opportunistic human pathogen that causes tuberculosis‐like chronic pulmonary disease and mortality stemming from comorbidities and treatment failure. The standard treatment of Mk infections requires costly, long‐term, multidrug courses with adverse side effects. The emergence of drug‐resistant isolates further complicates the already challenging drug therapy regimens and threatens to compromise the future control of Mk infections. Despite the increasingly recognized global burden of Mk infections, the biology of this opportunistic pathogen remains essentially unexplored. In particular, studies reporting gene function or generation of defined mutants are scarce. Moreover, no transposon (Tn) mutagenesis tool has been validated for use in Mk, a situation limiting the repertoire of genetic approaches available to accelerate the dissection of gene function and the generation of gene knockout mutants in this poorly characterized pathogen. In this study, we validated the functionality of a powerful Tn mutagenesis tool in Mk and used this tool in conjunction with a forward genetic screen to establish a previously unrecognized role of a conserved mycobacterial small RNA gene of unknown function in colony morphology features and biofilm formation. We also combined Tn mutagenesis with next‐generation sequencing to identify 12,071 Tn insertions that do not compromise viability in vitro. Finally, we demonstrated the susceptibility of the Galleria mellonella larva to Mk, setting the stage for further exploration of this simple and economical infection model system to the study of this pathogen.

Classification of Riboswitch Families Using Block Location-Based Feature Extraction (BLBFE) Method

Purpose: Riboswitches are special non-coding sequences usually located in mRNAs' un-translated regions and regulate gene expression and consequently cellular function. Furthermore, their interaction with antibiotics has been recently implicated. This raises more interest in development of bioinformatics tools for riboswitch studies. Herein, we describe the development and employment of novel block location-based feature extraction (BLBFE) method for classification of riboswitches. Methods: We have already developed and reported a sequential block finding (SBF) algorithm which, without operating alignment methods, identifies family specific sequential blocks for riboswitch families. Herein, we employed this algorithm for 7 riboswitch families including lysine, cobalamin, glycine, SAM-alpha, SAM-IV, cyclic-di-GMP-I and SAH. Then the study was extended toward implementation of BLBFE method for feature extraction. The outcome features were applied in various classifiers including linear discriminant analysis (LDA), probabilistic neural network (PNN), decision tree and k-nearest neighbors (KNN) classifiers for classification of the riboswitch families. The performance of the classifiers was investigated according to performance measures such as correct classification rate (CCR), accuracy, sensitivity, specificity and f-score. Results: As a result, average CCR for classification of riboswitches was 87.87%. Furthermore, application of BLBFE method in 4 classifiers displayed average accuracies of 93.98% to 96.1%, average sensitivities of 76.76% to 83.61%, average specificities of 96.53% to 97.69% and average f-scores of 74.9% to 81.91%. Conclusion: Our results approved that the proposed method of feature extraction; i.e. BLBFE method; can be successfully used for classification and discrimination of the riboswitch families with high CCR, accuracy, sensitivity, specificity and f-score values.

Cryo-EM structure of a 40 kDa SAM-IV riboswitch RNA at 3.7 Å resolution

Specimens below 50 kDa have generally been considered too small to be analyzed by single-particle cryo-electron microscopy (cryo-EM). The high flexibility of pure RNAs makes it difficult to obtain high-resolution structures by cryo-EM. In bacteria, riboswitches regulate sulfur metabolism through binding to the S-adenosylmethionine (SAM) ligand and offer compelling targets for new antibiotics. SAM-I, SAM-I/IV, and SAM-IV are the three most commonly found SAM riboswitches, but the structure of SAM-IV is still unknown. Here, we report the structures of apo and SAM-bound SAM-IV riboswitches (119-nt, ~40 kDa) to 3.7 Å and 4.1 Å resolution, respectively, using cryo-EM. The structures illustrate homologies in the ligand-binding core but distinct peripheral tertiary contacts in SAM-IV compared to SAM-I and SAM-I/IV. Our results demonstrate the feasibility of resolving small RNAs with enough detail to enable detection of their ligand-binding pockets and suggest that cryo-EM could play a role in structure-assisted drug design for RNA.

Mycobacterium tuberculosis 6C sRNA binds multiple mRNA targets via C-rich loops independent of RNA chaperones

Bacterial small regulatory RNAs (sRNAs) are the most abundant class of post-transcriptional regulators and have been well studied in Gram-negative bacteria. Little is known about the functions and mechanisms of sRNAs in high GC Gram-positive bacteria including Mycobacterium and Streptomyces. Here, we performed an in-depth study of 6C sRNA of Mycobacterium tuberculosis, which is conserved among high GC Gram-positive bacteria. Forty-seven genes were identified as possible direct targets of 6C sRNA and 15 of them were validated using an in vivo translational lacZ fusion system. We found that 6C sRNA plays a pleotropic role and regulates genes involved in various cellular processes, including DNA replication and protein secretion. Mapping the interactions of 6C sRNA with mRNA targets panD and dnaB revealed that the C-rich loops of 6C sRNA act as direct binding sites to mRNA targets. Unlike in Gram-negative bacteria where RNA binding proteins Hfq and ProQ are required, the interactions of 6C sRNA with mRNAs appear to be independent of RNA chaperones. Our findings suggest that the multiple G–C pairings between single stranded regions are sufficient to establish stable interactions between 6C sRNA and mRNA targets, providing a mechanism for sRNAs in high GC Gram-positive bacteria.

Regulation of Bacterial Gene Expression by Transcription Attenuation

A wide variety of mechanisms that control gene expression in bacteria are based on conditional transcription termination. Generally, in these mechanisms, a transcription terminator is located between a promoter and a downstream gene(s), and the efficiency of the terminator is controlled by a regulatory effector that can be a metabolite, protein, or RNA. The most common type of regulation involving conditional termination is transcription attenuation, in which the primary regulatory target is an essential element of a single terminator. The terminator can be either intrinsic or Rho dependent, with each presenting unique regulatory targets. Transcription attenuation mechanisms can be divided into five classes based primarily on the manner in which transcription termination is rendered conditional. This review summarizes each class of control mechanisms from a historical perspective, describes important examples in a physiological context and the current state of knowledge, highlights major advances, and discusses expectations of future discoveries.

Structural prediction of RNA switches using conditional base-pair probabilities

An RNA switch triggers biological functions by toggling between two conformations. RNA switches include bacterial riboswitches, where ligand binding can stabilize a bound structure. For RNAs with only one stable structure, structural prediction usually just requires a straightforward free energy minimization, but for an RNA switch, the prediction of a less stable alternative structure is often computationally costly and even problematic. The current sampling-clustering method predicts stable and alternative structures by partitioning structures sampled from the energy landscape into two clusters, but it is very time-consuming. Instead, we predict the alternative structure of an RNA switch from conditional probability calculations within the energy landscape. First, our method excludes base pairs related to the most stable structure in the energy landscape. Then, it detects stable stems (“seeds”) in the remaining landscape. Finally, it folds an alternative structure prediction around a seed. While having comparable riboswitch classification performance, the conditional-probability computations had fewer adjustable parameters, offered greater predictive flexibility, and were more than one thousand times faster than the sampling step alone in sampling-clustering predictions, the competing standard. Overall, the described approach helps traverse thermodynamically improbable energy landscapes to find biologically significant substructures and structures rapidly and effectively.

Two-Component Systems in Francisella Species

Bacteria alter gene expression in response to changes in their environment through various mechanisms that include signal transduction systems. These signal transduction systems use membrane histidine kinase with sensing domains to mediate phosphotransfer to DNA-binding proteins that alter the level of gene expression. Such regulators are called two-component systems (TCSs). TCSs integrate external signals and information from stress pathways, central metabolism and other global regulators, thus playing an important role as part of the overall regulatory network. This review will focus on the knowledge of TCSs in the Gram-negative bacterium, Francisella tularensis, a biothreat agent with a wide range of potential hosts and a significant ability to cause disease. While TCSs have been well-studied in several bacterial pathogens, they have not been well-studied in non-model organisms, such as F. tularensis and its subspecies, whose canonical TCS content surprisingly ranges from few to none. Additionally, of those TCS genes present, many are orphan components, including KdpDE, QseC, QseB/PmrA, and an unnamed two-component system (FTN_1452/FTN_1453). We discuss recent advances in this field related to the role of TCSs in Francisella physiology and pathogenesis and compare the TCS genes present in human virulent versus. environmental species and subspecies of Francisella.

Guanidine Riboswitch-Regulated Efflux Transporters Protect Bacteria against Ionic Liquid Toxicity

Plant cell walls contain a renewable, nearly limitless supply of sugar that could be used to support microbial production of commodity chemicals and biofuels. Imidazolium ionic liquid (IIL) solvents are among the best reagents for gaining access to the sugars in this otherwise recalcitrant biomass. However, the sugars from IIL-treated biomass are inevitably contaminated with residual IILs that inhibit growth in bacteria and yeast, blocking biochemical production by these organisms. IIL toxicity is, therefore, a critical roadblock in many industrial biosynthetic pathways. Although several IIL-tolerant (IIL^T) bacterial and yeast isolates have been identified in nature, few genetic mechanisms have been identified. In this study, we identified two IIL^T Bacillus isolates as well as a spontaneous IIL^T Escherichia coli lab strain that are tolerant to high levels of two widely used IILs. We demonstrate that all three IIL^T strains contain one or more pumps of the small multidrug resistance (SMR) family, and two of these strains contain mutations that affect an adjacent regulatory guanidine riboswitch. Furthermore, we show that the regulation of E. coli sugE by the guanidine II riboswitch can be exploited to promote IIL tolerance by the simple addition of guanidine to the medium. Our results demonstrate the critical role that transporter genes play in IIL tolerance in their native bacterial hosts. The study presented here is another step in engineering IIL tolerance into industrial strains toward overcoming this key gap in biofuels and industrial biochemical production processes.IMPORTANCE This study identifies bacteria that are tolerant to ionic liquid solvents used in the production of biofuels and industrial biochemicals. For industrial microbiology, it is essential to find less-harmful reagents and microbes that are resistant to their cytotoxic effects. We identified a family of small multidrug resistance efflux transporters, which are responsible for the tolerance of these strains. We also found that this resistance can be caused by mutations in the sequences of guanidine-specific riboswitches that regulate these efflux pumps. Extending this knowledge, we demonstrated that guanidine itself can promote ionic liquid tolerance. Our findings will inform genetic engineering strategies that improve conversion of cellulosic sugars into biofuels and biochemicals in processes where low concentrations of ionic liquids surpass bacterial tolerance.

The electrifying physiology of Geobacter bacteria, 30 years on

The family Geobacteraceae, with its only valid genus Geobacter, comprises deltaproteobacteria ubiquitous in soil, sediments, and subsurface environments where metal reduction is an active process. Research for almost three decades has provided novel insights into environmental processes and biogeochemical reactions not previously known to be carried out by microorganisms. At the heart of the environmental roles played by Geobacter bacteria is their ability to integrate redox pathways and regulatory checkpoints that maximize growth efficiency with electron donors derived from the decomposition of organic matter while respiring metal oxides, particularly the often abundant oxides of ferric iron. This metabolic specialization is complemented by versatile metabolic reactions, respiratory chains, and sensory networks that allow specific members to adaptively respond to environmental cues to integrate organic and inorganic contaminants in their oxidative and reductive metabolism, respectively. Thus, Geobacteraceae are important members of the microbial communities that degrade hydrocarbon contaminants under iron-reducing conditions and that contribute, directly or indirectly, to the reduction of radionuclides, toxic metals, and oxidized species of nitrogen. Their ability to produce conductive pili as nanowires for discharging respiratory electrons to solid-phase electron acceptors and radionuclides, or for wiring cells in current-harvesting biofilms highlights the unique physiological traits that make these organisms attractive biological platforms for bioremediation, bioenergy, and bioelectronics application. Here we review some of the most notable physiological features described in Geobacter species since the first model representatives were recovered in pure culture. We provide a historical account of the environmental research that has set the foundation for numerous physiological studies and the laboratory tools that had provided novel insights into the role of Geobacter in the functioning of microbial communities from pristine and contaminated environments. We pay particular attention to latest research, both basic and applied, that has served to expand the field into new directions and to advance interdisciplinary knowledge. The electrifying physiology of Geobacter, it seems, is alive and well 30 years on.

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.