Large Noncoding RNAs in Bacteria

Evidence that ribosomal protein bS21 is a component of the OLE ribonucleoprotein complex

OLE RNAs represent a large and highly structured noncoding RNA (ncRNA) class that is mostly found in Gram-positive extremophiles and/or anaerobes of the Bacillota phylum. These ~600-nucleotide RNAs are among the most structurally complex and well-conserved large ncRNAs whose precise biochemical functions remain to be established. In Halalkalibacterium halodurans, OLE RNA is involved in the adaptation to various unfavourable growth conditions, including exposure to cold (≤20°C), ethanol (≥3% [v/v]), excess Mg2+ (≥4 mM), and non-glucose carbon/energy sources. OLE forms a ribonucleoprotein (RNP) complex with the OLE-associated proteins A, B and C, which are known to be essential for OLE RNP complex function in this species. Bacteria lacking OLE RNA (Δole) or a functional OLE RNP complex exhibit growth defects under the stresses listed above. Here, we demonstrate that ribosomal protein bS21 is a natural component of the OLE RNP complex and we map its precise RNA binding site. The presence of bS21 results in a conformational change in OLE RNA resembling a k-turn substructure previously reported to be relevant to the function of the OLE RNP complex. Mutational disruption of the bS21 protein or its OLE RNA binding site results in growth inhibition under cold and ethanol stress to the same extent as the deletion of the gene for OLE RNA. These findings are consistent with the hypothesis that bS21 is a biologically relevant component of the OLE RNP complex under a subset of stresses managed by the OLE RNP complex.

Transcriptional stochasticity reveals multiple mechanisms of long non-coding RNA regulation at the Xist-Tsix locus

Long noncoding RNAs (LncRNAs) are increasingly recognized as being involved in human physiology and diseases, but there is a lack of mechanistic understanding for the majority of lncRNAs. We comparatively test proposed mechanisms of antisense lncRNA regulation at the X-chromosome Inactivation (XCI) locus. We find that due to stochasticity in transcription, different mechanisms based on the act of transcription regulate Xist and Tsix at different levels of nascent transcription. At medium levels, RNA polymerases transcribe Xist and Tsix on each strand at the same transcription site and deposit significant amounts of the histone mark H3K36me3, which inhibits Xist. At high levels of nascent transcription, many RNA polymerases transcribe Xist or Tsix resulting in transcriptional interference. Therefore, lncRNA expression variability is not just a quirk of transcription but an important aspect of regulation that allows multiple mechanisms to be employed by the same gene locus within the same cell population.

Molecular process control for industrial biotechnology

The development of sustainable and economically competitive biotechnological processes is a central challenge of modern industrial biotechnology. Conventional strategies such as macroscopic and molecular bioprocess design are often insufficient to exploit their full potential. To circumvent this, molecular process control provides the missing link to further consolidate various optimization strategies to achieve multilayered process design. This review highlights the molecular mechanisms that can be exploited for molecular process control. These can either be endogenous or specifically implemented into the organism, and comprise regulatory mechanisms at the transcriptional, translational, and system levels. In addition to serving as a design tool to enhance existing bioprocesses, molecular process control is the gateway to biotechnological advances that will extend the boundaries of future process design.

Cryo-EM reveals mechanisms of natural RNA multivalency

Homo-oligomerization of biological macromolecules leads to functional assemblies that are critical to understanding various cellular processes. However, RNA quaternary structures have rarely been reported. Comparative genomics analysis has identified RNA families containing hundreds of sequences that adopt conserved secondary structures and likely fold into complex three-dimensional structures. In this study, we used cryo-electron microscopy (cryo-EM) to determine structures from four RNA families, including ARRPOF and OLE forming dimers and ROOL and GOLLD forming hexameric, octameric, and dodecameric nanostructures, at 2.6- to 4.6-angstrom resolutions. These homo-oligomeric assemblies reveal a plethora of structural motifs that contribute to RNA multivalency, including kissing-loop, palindromic base-pairing, A-stacking, metal ion coordination, pseudoknot, and minor-groove interactions. These results provide the molecular basis of intermolecular interactions driving RNA multivalency with potential functional relevance.

The role of noncoding RNAs in bacterial immunity

The evolutionary arms race between bacteria and phages has driven the development of diverse anti-phage defense mechanisms. Recent studies have identified noncoding RNAs (ncRNAs) as key players in bacteria-phage conflicts, including CRISPR-Cas, toxin-antitoxin (TA), and reverse transcriptase (RT)-based defenses; however, our understanding of their roles in immunity is still emerging. In this review, we explore the multifaceted roles of ncRNAs in bacterial immunity, offering insights into their contributions to defense and anti-defense mechanisms, their influence on immune regulatory networks, and potential biotechnological applications. Finally, we highlight key outstanding questions in the field to spark future research directions.

Cryo-EM Structure of raiA ncRNA From Clostridium Reveals a New RNA 3D Fold

Advancements in genome-wide sequence analysis have led to the discovery of numerous novel bacterial non-coding RNAs (ncRNAs). These ncRNAs have been categorized into various RNA families and classes based on their size, structure, function, and evolutionary relationships. One such ncRNA family, raiA, is notably abundant in the bacterial phyla Firmicutes and Actinobacteria and is remarkably well-conserved across many Gram-positive bacteria. In this study, we integrated cryo-electron microscopy single-particle analysis with computational modeling and biochemical techniques to elucidate the structural characteristics of raiA from Clostridium sp. CAG 138. Our findings reveal the globular 3D fold of raiA, providing valuable structural insights. This analysis paves the way for future investigations into the functional properties of raiA, potentially uncovering new regulatory mechanisms in bacterial ncRNAs.

The Escherichia coli ribB riboswitch senses flavin mononucleotide within a defined transcriptional window

Riboswitches are metabolite-binding RNA regulators that modulate gene expression at the levels of transcription and translation. One of the hallmarks of riboswitch regulation is that they undergo structural changes upon metabolite binding. While a lot of effort has been put to characterize how the metabolite is recognized by the riboswitch, there is still relatively little information regarding how ligand sensing is performed within a transcriptional context. Here, we study the ligand-dependent cotranscriptional folding of the FMN-sensing ribB riboswitch of Escherichia coli Using RNase H assays to study nascent ribB riboswitch transcripts, DNA probes targeting the P1 and sequestering stems indicate that FMN binding leads to the protection of these regions from RNase H cleavage, consistent with the riboswitch inhibiting translation initiation when bound to FMN. Our results show that ligand sensing is strongly affected by the position of elongating RNA polymerase, which is defining an FMN-binding transcriptional window that is bordered in its 3' extremity by a transcriptional pause site. Also, using successively overlapping DNA probes targeting a subdomain of the riboswitch, our data suggest the presence of a previously unsuspected helical region involving the 3' strand of the P1 stem. Our results show that this helical region is conserved across bacterial species, thus suggesting that this predicted structure, the anti*-P1 stem, is involved in the FMN-free conformation of the ribB riboswitch. Overall, our study further demonstrates that intricate folding strategies may be used by riboswitches to perform metabolite sensing during the transcriptional process.

Exploring the hidden hot world of long non-coding RNAs in thermophilic fungus using a robust computational pipeline

Long noncoding RNAs (lncRNAs) are versatile RNA molecules recently identified as key regulators of gene expression in response to environmental stress. Our primary focus in this study was to develop a robust computational pipeline for identifying structurally identical lncRNAs across replicates from publicly available bulk RNA-seq datasets. In order to demonstrate the effectiveness of the pipeline, we utilized the transcriptome of the thermophilic fungus Thermothelomyces thermophilus and assessed the expression pattern of lncRNAs in conjunction with Heat Shock Proteins (HSP), a well-known protein family critical for the cell's response to high temperatures. Our findings demonstrate that the identification of structurally identical transcripts among replicates in this thermophilic fungus ensures the reliability and accuracy of RNA studies, contributing to the validity of biological interpretations. Furthermore, the majority of lncRNAs exhibited a distinct expression pattern compared to HSPs. Our study contributes to advancing the understanding of the biological mechanisms comprising lncRNAs in thermophilic fungi.

Bioinformatic prediction of proteins relevant to functions of the bacterial OLE ribonucleoprotein complex

OLE (ornate, large, extremophilic) RNAs are members of a noncoding RNA class present in many Gram-positive, extremophilic bacteria. The large size, complex structure, and extensive sequence conservation of OLE RNAs are characteristics consistent with the hypothesis that they likely function as ribozymes. The OLE RNA representative from Halalkalibacterium halodurans is known to localize to the phospholipid membrane and requires at least three essential protein partners: OapA, OapB, and OapC. However, the precise biochemical functions of this unusual ribonucleoprotein (RNP) complex remain unknown. Genetic disruption of OLE RNA or its partners revealed that the complex is beneficial under diverse stress conditions. To search for additional links between OLE RNA and other cellular components, we used phylogenetic profiling to identify proteins that are either correlated or anticorrelated with the presence of OLE RNA in various bacterial species. This analysis revealed strong correlations between the essential protein-binding partners of OLE RNA and organisms that carry the ole gene. Similarly, proteins involved in sporulation are correlated, suggesting a potential role for the OLE RNP complex in spore formation. Intriguingly, the Mg2+ transporter MpfA is strongly anticorrelated with OLE RNA. Evidence indicates that MpfA is structurally related to OapA and therefore MpfA may serve as a functional replacement for some contributions otherwise performed by the OLE RNP complex in species that lack this device. Indeed, OLE RNAs might represent an ancient RNA class that enabled primitive organisms to sense and respond to major cellular stresses.IMPORTANCEOLE (ornate, large, extremophilic) RNAs were first reported nearly 20 years ago, and they represent one of the largest and most intricately folded noncoding RNA classes whose biochemical function remains to be established. Other RNAs with similar size, structural complexity, and extent of sequence conservation have proven to catalyze chemical transformations. Therefore, we speculate that OLE RNAs likewise operate as ribozymes and that they might catalyze a fundamental reaction that has persisted since the RNA World era-a time before the emergence of proteins in evolution. To seek additional clues regarding the function of OLE RNA, we undertook a computational effort to identify potential protein components of the OLE ribonucleoprotein (RNP) complex or other proteins that have functional links to this device. This analysis revealed known protein partners and several additional proteins that might be physically or functionally linked to the OLE RNP complex. Finally, we identified a Mg2+ transporter protein, MpfA, that strongly anticorrelates with the OLE RNP complex. This latter result suggests that MpfA might perform at least some functions that are like those carried out by the OLE RNP complex.

The discovery of novel noncoding RNAs in 50 bacterial genomes

Structured noncoding RNAs (ncRNAs) contribute to many important cellular processes involving chemical catalysis, molecular recognition and gene regulation. Few ncRNA classes are broadly distributed among organisms from all three domains of life, but the list of rarer classes that exhibit surprisingly diverse functions is growing. We previously developed a computational pipeline that enables the near-comprehensive identification of structured ncRNAs expressed from individual bacterial genomes. The regions between protein coding genes are first sorted based on length and the fraction of guanosine and cytidine nucleotides. Long, GC-rich intergenic regions are then examined for sequence and structural similarity to other bacterial genomes. Herein, we describe the implementation of this pipeline on 50 bacterial genomes from varied phyla. More than 4700 candidate intergenic regions with the desired characteristics were identified, which yielded 44 novel riboswitch candidates and numerous other putative ncRNA motifs. Although experimental validation studies have yet to be conducted, this rate of riboswitch candidate discovery is consistent with predictions that many hundreds of novel riboswitch classes remain to be discovered among the bacterial species whose genomes have already been sequenced. Thus, many thousands of additional novel ncRNA classes likely remain to be discovered in the bacterial domain of life.

Genetic disruption of the bacterial raiA motif noncoding RNA causes defects in sporulation and aggregation

Several structured noncoding RNAs in bacteria are essential contributors to fundamental cellular processes. Thus, discoveries of additional ncRNA classes provide opportunities to uncover and explore biochemical mechanisms relevant to other major and potentially ancient processes. A candidate structured ncRNA named the "raiA motif" has been found via bioinformatic analyses in over 2,500 bacterial species. The gene coding for the RNA typically resides between the raiA and comFC genes of many species of Bacillota and Actinomycetota. Structural probing of the raiA motif RNA from the Gram-positive anaerobe Clostridium acetobutylicum confirms key features of its sophisticated secondary structure model. Expression analysis of raiA motif RNA reveals that the RNA is constitutively produced but reaches peak abundance during the transition from exponential growth to stationary phase. The raiA motif RNA becomes the fourth most abundant RNA in C. acetobutylicum, excluding ribosomal RNAs and transfer RNAs. Genetic disruption of the raiA motif RNA causes cells to exhibit substantially decreased spore formation and diminished ability to aggregate. Restoration of normal cellular function in this knock-out strain is achieved by expression of a raiA motif gene from a plasmid. These results demonstrate that raiA motif RNAs normally participate in major cell differentiation processes by operating as a trans-acting factor.

Disruption of the bacterial OLE RNP complex impairs growth on alternative carbon sources

Ornate, large, extremophilic (OLE) RNAs comprise a class of large noncoding RNAs in bacteria whose members form a membrane-associated ribonucleoprotein (RNP) complex. This complex facilitates cellular adaptation to diverse stresses such as exposure to cold, short-chain alcohols, and elevated Mg2+ concentrations. Here, we report additional phenotypes exhibited by Halalkalibacterium halodurans (formerly called Bacillus halodurans) strains lacking functional OLE RNP complexes. Genetic disruption of the complex causes restricted growth compared to wild-type cells when cultured in minimal media (MM) wherein glucose is replaced with alternative carbon/energy sources. Genetic suppressor selections conducted in glutamate MM yielded isolates that carry mutations in or near genes relevant to Mn2+ homeostasis (ykoY and mntB), phosphate homeostasis (phoR), and putative multidrug resistance (bmrCD). These functional links between OLE RNA, carbon/energy management, and other fundamental processes including protein secretion are consistent with the hypothesis that the OLE RNP complex is a major contributor to cellular adaptation to unfavorable growth conditions.

Regulation of bacterial gene expression by non-coding RNA: It is all about time!

Commensal and pathogenic bacteria continuously evolve to survive in diverse ecological niches by efficiently coordinating gene expression levels in their ever-changing environments. Regulation through the RNA transcript itself offers a faster and more cost-effective way to adapt than protein-based mechanisms and can be leveraged for diagnostic or antimicrobial purposes. However, RNA can fold into numerous intricate, not always functional structures that both expand and obscure the plethora of roles that regulatory RNAs serve within the cell. Here, we review the current knowledge of bacterial non-coding RNAs in relation to their folding pathways and interactions. We posit that co-transcriptional folding of these transcripts ultimately dictates their downstream functions. Elucidating the spatiotemporal folding of non-coding RNAs during transcription therefore provides invaluable insights into bacterial pathogeneses and predictive disease diagnostics. Finally, we discuss the implications of co-transcriptional folding andapplications of RNAs for therapeutics and drug targets.

Evidence that OLE RNA is a component of a major stress-responsive ribonucleoprotein particle in extremophilic bacteria

OLE RNA is a ~600-nucleotide noncoding RNA present in many Gram-positive bacteria that thrive mostly in extreme environments, including elevated temperature, salt, and pH conditions. The precise biochemical functions of this highly conserved RNA remain unknown, but it forms a ribonucleoprotein (RNP) complex that localizes to cell membranes. Genetic disruption of the RNA or its essential protein partners causes reduced cell growth under various stress conditions. These phenotypes include sensitivity to short-chain alcohols, cold intolerance, reduced growth on sub-optimal carbon sources, and intolerance of even modest concentrations of Mg2+ . Thus, many bacterial species appear to employ OLE RNA as a component of an intricate RNP apparatus to monitor fundamental cellular processes and make physiological and metabolic adaptations. Herein we hypothesize that the OLE RNP complex is functionally equivalent to the eukaryotic TOR complexes, which integrate signals from various diverse pathways to coordinate processes central to cell growth, replication, and survival.

Transcriptional Landscapes of Herelleviridae Bacteriophages and Staphylococcus aureus during Phage Infection: An Overview

The issue of antibiotic resistance in healthcare worldwide has led to a pressing need to explore and develop alternative approaches to combat infectious diseases. Among these methods, phage therapy has emerged as a potential solution to tackle this growing challenge. Virulent phages of the Herelleviridae family, known for their ability to cause lysis of Staphylococcus aureus, a clinically significant pathogen frequently associated with multidrug resistance, have proven to be one of the most effective viruses utilized in phage therapy. In order to utilize phages for therapeutic purposes effectively, a thorough investigation into their physiology and mechanisms of action on infected cells is essential. The use of omics technologies, particularly total RNA sequencing, is a promising approach for analyzing the interaction between phages and their hosts, allowing for the assessment of both the behavior of the phage during infection and the cell's response. This review aims to provide a comprehensive overview of the physiology of the Herelleviridae family, utilizing existing analyses of their total phage transcriptomes. Additionally, it sheds light on the changes that occur in the metabolism of S. aureus when infected with virulent bacteriophages, contributing to a deeper understanding of the phage-host interaction.

Cotranscriptional Assembly and Native Purification of Large RNA-RNA Complexes for Structural Analyses

Recent technological developments such as cryogenic electron microscopy (Cryo-EM) and X-ray free electron lasers (XFEL) have significantly expanded the available toolkit to visualize large, complex noncoding RNAs and their complexes. Consequently, the quality of the RNA sample, as measured by its chemical monodispersity and conformational homogeneity, has become the bottleneck that frequently precludes effective structural analyses. Here we describe a general RNA sample preparation protocol that combines cotranscriptional RNA folding and RNA-RNA complex assembly, followed by native purification of stoichiometric complexes. We illustrate and discuss the utility of this versatile method in overcoming RNA misfolding and enabling the structural and mechanistic elucidations of the T-box riboswitch-tRNA complexes.

Compact Cas9d and HEARO enzymes for genome editing discovered from uncultivated microbes

Programmable, RNA-guided nucleases are diverse enzymes that have been repurposed for biotechnological applications. However, to further expand the therapeutic application of these tools there is a need for targetable systems that are small enough to be delivered efficiently. Here, we mined an extensive genome-resolved metagenomics database and identified families of uncharacterized RNA-guided, compact nucleases (between 450 and 1,050 aa). We report that Cas9d, a new CRISPR type II subtype, contains Zinc-finger motifs and high arginine content, features that we also found in nucleases related to HEARO effectors. These enzymes exhibit diverse biochemical characteristics and are broadly targetable. We show that natural Cas9d enzymes are capable of genome editing in mammalian cells with >90% efficiency, and further engineered nickase variants into the smallest base editors active in E. coli and human cells. Their small size, broad targeting potential, and translatability suggest that Cas9d and HEARO systems will enable a variety of genome editing applications.

Ornate, large, extremophilic (OLE) RNA forms a kink turn necessary for OapC protein recognition and RNA function

Ornate, large, extremophilic (OLE) RNAs represent a class of noncoding RNAs prevalent in Gram-positive, extremophilic/anaerobic bacterial species. OLE RNAs (∼600 nt), whose precise biochemical functions remain mysterious, form an intricate secondary structure interspersed with regions of highly conserved nucleotides. In the alkali-halophilic bacterium Bacillus halodurans, OLE RNA is a component of a ribonucleoprotein (RNP) complex involving at least two proteins named OapA and OapB, but additional components may exist that could point to functional roles for the RNA. Disruption of the genes for either OLE RNA, OapA, or OapB result in the inability of cells to overcome cold, alcohol, or Mg²⁺ stresses. In the current study, we used in vivo crosslinking followed by OLE RNA isolation to identify the protein YbxF as a potential additional partner in the OLE RNP complex. Notably, a mutation in the gene for this same protein was also reported to be present in a strain wherein the complex is nonfunctional. The B. halodurans YbxF (herein renamed OapC) is homologous to a bacterial protein earlier demonstrated to bind kink turn (k-turn) RNA structural motifs. In vitro RNA-protein binding assays reveal that OLE RNA forms a previously unrecognized k-turn that serves as the natural binding site for YbxF/OapC. Moreover, B. halodurans cells carrying OLE RNAs with disruptive mutations in the k-turn exhibit phenotypes identical to cells lacking functional OLE RNP complexes. These findings reveal that the YbxF/OapC protein of B. halodurans is important for the formation of a functional OLE RNP complex.

Prokaryotic ncRNAs: Master regulators of gene expression

ncRNA plays a very pivotal role in various biological activities ranging from gene regulation to controlling important developmental networks. It is imperative to note that this small molecule is not only present in all three domains of cellular life, but is an important modulator of gene regulation too in all these domains. In this review, we discussed various aspects of ncRNA biology, especially their role in bacteria. The last two decades of scientific research have proved that this molecule plays an important role in the modulation of various regulatory pathways in bacteria including the adaptive immune system and gene regulation. It is also very surprising to note that this small molecule is also employed in various processes related to the pathogenicity of virulent microorganisms.

Global Transcriptomic Analysis of Bacteriophage-Host Interactions between a Kayvirus Therapeutic Phage and Staphylococcus aureus

Kayviruses are polyvalent broad host range staphylococcal phages with a potential to combat staphylococcal infections. However, the implementation of rational phage therapy in medicine requires a thorough understanding of the interactions between bacteriophages and pathogens at omics level. To evaluate the effect of a phage used in therapy on its host bacterium, we performed differential transcriptomic analysis by RNA-Seq from bacteriophage K of genus Kayvirus infecting two Staphylococcus aureus strains, prophage-less strain SH1000 and quadruple lysogenic strain Newman. The temporal transcriptional profile of phage K was comparable in both strains except for a few loci encoding hypothetical proteins. Stranded sequencing revealed transcription of phage noncoding RNAs that may play a role in the regulation of phage and host gene expression. The transcriptional response of S. aureus to phage K infection resembles a general stress response with differential expression of genes involved in a DNA damage response. The host transcriptional changes involved upregulation of nucleotide, amino acid and energy synthesis and transporter genes and downregulation of host transcription factors. The interaction of phage K with variable genetic elements of the host showed slight upregulation of gene expression of prophage integrases and antirepressors. The virulence genes involved in adhesion and immune evasion were only marginally affected, making phage K suitable for therapy. IMPORTANCE Bacterium Staphylococcus aureus is a common human and veterinary pathogen that causes mild to life-threatening infections. As strains of S. aureus are becoming increasingly resistant to multiple antibiotics, the need to search for new therapeutics is urgent. A promising alternative to antibiotic treatment of staphylococcal infections is a phage therapy using lytic phages from the genus Kayvirus. Here, we present a comprehensive view on the phage-bacterium interactions on transcriptomic level that improves the knowledge of molecular mechanisms underlying the Kayvirus lytic action. The results will ensure safer usage of the phage therapeutics and may also serve as a basis for the development of new antibacterial strategies.

Genomic insights into the diversity of non-coding RNAs in Bacillus cereus sensu lato

Bacillus cereus sensu lato is a group of bacteria of medical and agricultural importance in different ecological niches and with controversial taxonomic relationships. Studying the composition of non-coding RNAs (ncRNAs) in several bacterial groups has been an important tool for identifying genetic information and better understanding genetic regulation towards environment adaptation. However, to date, no comparative genomics study of ncRNA has been performed in this group. Thus, this study aimed to identify and characterize the set of ncRNAs from 132 strains of Bacillus cereus, Bacillus thuringiensis and Bacillus anthracis to obtain an overview of the diversity and distribution of these genetic elements in these species. We observed that the number of ncRNAs differs in the chromosomes of the three species, but not in the plasmids, when species or phylogenetic clusters were compared. The prevailing functional/structural category was Cis-reg and the most frequent class was Riboswitch. However, in plasmids, the class Group II intron was the most frequent. Also, nine ncRNAs were selected for validation in the strain B. thuringiensis 407 by RT-PCR, which allowed to identify the expression of the ncRNAs. The wide distribution and diversity of ncRNAs in the B. cereus group, and more intensely in B. thuringiensis, may help improve the abilities of these species to adapt to various environmental changes. Further studies should address the expression of these genetic elements in different conditions.

Study of the diversity of 16S-23S rDNA internal transcribed spacer (ITS) typing of Escherichia coli strains isolated from various biotopes in Tunisia

We investigated the 16S-23S rRNA intergenic spacer region (ISR)-PCR and the phylogenetic PCR analyses of 150 Escherichia coli isolates as tools to explore their diversity, according to their sampling origins, and their relative dominance in these sampling sources. These genetic markers are used to explore phylogenetic and genetic relationships of these 150 E. coli isolates recovered from different environmental sources (water, food, animal, human and vegetables). These isolates are tested for their biochemical pattern and later genotyped through the 16S-23S rRNA intergenic spacer PCR amplification and their polymorphism investigation of PCR-amplified 16S-23S rDNA ITS. The main results of the pattern band profile revealed one to four DNA fragments. Distributing 150 E. coli isolates according to their ITS and using RS-PCR, revealed four genotypes and four subtypes. The DNA fragment size ranged from 450 to 550 bp. DNA band patterns analysis revealed considerable genetic diversity in interspecies. Thus, the 450 and 550 bp sizes of the common bands in all E. coli isolates are highly diversified. Genotype I appeared as the most frequent with 77.3% (116 isolates), genotype II with 12% (18 isolates); genotype III with 9.7% (14 isolates), and the IV rarely occurred with 4% (2 isolates). Distributing the E. coli phylogroups showed 84 isolates (56%) of group A, 35 isolates (23.3%) of group B1, 28 isolates (18.7%) of group B2 and only three isolates (2%) of group D.

Natural circularly permuted group II introns in bacteria produce RNA circles

Group II self-splicing introns are large structured RNAs that remove themselves from transcripts while simultaneously sealing the resulting gaps. Some representatives can subsequently reverse splice into DNA, accounting for their pervasive distribution in bacteria. The catalytically active tertiary structure of each group II intron is assembled from six domains that are arranged in a conserved order. Here, we report structural isomers of group II introns, called CP group II ribozymes, wherein the characteristic order of domains has been altered. Domains five and six, which normally reside at the 3' end of group II introns, instead occupy the 5' end to form circularly permuted variants. These unusual group II intron derivatives are catalytically active and generate large linear branched and small circular RNAs, reaction products that are markedly different from those generated by canonical group II introns. The biological role of CP group II ribozymes is currently unknown.

A workflow to identify novel proteins based on the direct mapping of peptide-spectrum-matches to genomic locations

Background

Small Proteins have received increasing attention in recent years. They have in particular been implicated as signals contributing to the coordination of bacterial communities. In genome annotations they are often missing or hidden among large numbers of hypothetical proteins because genome annotation pipelines often exclude short open reading frames or over-predict hypothetical proteins based on simple models. The validation of novel proteins, and in particular of small proteins (sProteins), therefore requires additional evidence. Proteogenomics is considered the gold standard for this purpose. It extends beyond established annotations and includes all possible open reading frames (ORFs) as potential sources of peptides, thus allowing the discovery of novel, unannotated proteins. Typically this results in large numbers of putative novel small proteins fraught with large fractions of false-positive predictions.

Results

We observe that number and quality of the peptide-spectrum matches (PSMs) that map to a candidate ORF can be highly informative for the purpose of distinguishing proteins from spurious ORF annotations. We report here on a workflow that aggregates PSM quality information and local context into simple descriptors and reliably separates likely proteins from the large pool of false-positive, i.e., most likely untranslated ORFs. We investigated the artificial gut microbiome model SIHUMIx, comprising eight different species, for which we validate 5114 proteins that have previously been annotated only as hypothetical ORFs. In addition, we identified 37 non-annotated protein candidates for which we found evidence at the proteomic and transcriptomic level. Half (19) of these candidates have close functional homologs in other species. Another 12 candidates have homologs designated as hypothetical proteins in other species. The remaining six candidates are short (< 100 AA) and are most likely bona fide novel proteins.

Conclusions

The aggregation of PSM quality information for predicted ORFs provides a robust and efficient method to identify novel proteins in proteomics data. The workflow is in particular capable of identifying small proteins and frameshift variants. Since PSMs are explicitly mapped to genomic locations, it furthermore facilitates the integration of transcriptomics data and other sources of genome-level information.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04159-8.

DRAGoM: Classification and Quantification of Noncoding RNA in Metagenomic Data

Noncoding RNAs (ncRNAs) play important regulatory and functional roles in microorganisms, such as regulation of gene expression, signaling, protein synthesis, and RNA processing. Hence, their classification and quantification are central tasks toward the understanding of the function of the microbial community. However, the majority of the current metagenomic sequencing technologies generate short reads, which may contain only a partial secondary structure that complicates ncRNA homology detection. Meanwhile, de novo assembly of the metagenomic sequencing data remains challenging for complex communities. To tackle these challenges, we developed a novel algorithm called DRAGoM (Detection of RNA using Assembly Graph from Metagenomic data). DRAGoM first constructs a hybrid graph by merging an assembly string graph and an assembly de Bruijn graph. Then, it classifies paths in the hybrid graph and their constituent readsinto differentncRNA families based on both sequence and structural homology. Our benchmark experiments show that DRAGoMcan improve the performance and robustness over traditional approaches on the classification and quantification of a wide class of ncRNA families.

The World of Stable Ribonucleoproteins and Its Mapping With Grad-Seq and Related Approaches

Macromolecular complexes of proteins and RNAs are essential building blocks of cells. These stable supramolecular particles can be viewed as minimal biochemical units whose structural organization, i.e., the way the RNA and the protein interact with each other, is directly linked to their biological function. Whether those are dynamic regulatory ribonucleoproteins (RNPs) or integrated molecular machines involved in gene expression, the comprehensive knowledge of these units is critical to our understanding of key molecular mechanisms and cell physiology phenomena. Such is the goal of diverse complexomic approaches and in particular of the recently developed gradient profiling by sequencing (Grad-seq). By separating cellular protein and RNA complexes on a density gradient and quantifying their distributions genome-wide by mass spectrometry and deep sequencing, Grad-seq charts global landscapes of native macromolecular assemblies. In this review, we propose a function-based ontology of stable RNPs and discuss how Grad-seq and related approaches transformed our perspective of bacterial and eukaryotic ribonucleoproteins by guiding the discovery of new RNA-binding proteins and unusual classes of noncoding RNAs. We highlight some methodological aspects and developments that permit to further boost the power of this technique and to look for exciting new biology in understudied and challenging biological models.

A cell-based ribozyme reporter system employing a chromosomally-integrated 5' exonuclease gene

Background

Bioinformatic genome surveys indicate that self-cleaving ribonucleic acids (ribozymes) appear to be widespread among all domains of life, although the functions of only a small number have been validated by biochemical methods. Alternatively, cell-based reporter gene assays can be used to validate ribozyme function. However, reporter activity can be confounded by phenomena unrelated to ribozyme-mediated cleavage of RNA.

Results

We established a ribozyme reporter system in Escherichia coli in which a significant reduction of reporter activity is manifest when an active ribozyme sequence is fused to the reporter gene and the expression of a foreign Bacillus subtilis RNaseJ1 5′ exonuclease is induced from a chromosomally-integrated gene in the same cell.

Conclusions

The reporter system could be useful for validating ribozyme function in candidate sequences identified from bioinformatics.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12860-021-00357-7.

Structure of a bacterial OapB protein with its OLE RNA target gives insights into the architecture of the OLE ribonucleoprotein complex

Bacterial noncoding RNAs (ncRNAs) play key roles in many biological processes including gene regulation, RNA processing and modification, and protein synthesis and translocation. OLE RNAs, found in many Gram-positive species, are one of the largest highly structured ncRNA classes whose biochemical functions remain unknown. In Bacillus halodurans, OLE RNAs interact with at least two proteins, OapA and OapB, which are required to assemble a functional OLE ribonucleoprotein (RNP) complex contributing to cellular responses to certain environmental stresses. We established X-ray structural models that reveal the sequence elements and tertiary structural features of OLE RNA that are critical for its specific recognition by OapB, which will aid future exploration of the biological and biochemical functions of the unusual OLE RNP complex.

The OLE (ornate, large, and extremophilic) RNA class is one of the most complex and well-conserved bacterial noncoding RNAs known to exist. This RNA is known to be important for bacterial responses to stress caused by short-chain alcohols, cold, and elevated Mg²⁺ concentrations. These biological functions have been shown to require the formation of a ribonucleoprotein (RNP) complex including at least two protein partners: OLE-associated protein A (OapA) and OLE-associated protein B (OapB). OapB directly binds OLE RNA with high-affinity and specificity and is believed to assist in assembling the functional OLE RNP complex. To provide the atomic details of OapB–OLE RNA interaction and to potentially reveal previously uncharacterized protein–RNA interfaces, we determined the structure of OapB from Bacillus halodurans alone and in complex with an OLE RNA fragment at resolutions of 1.0 Å and 2.0 Å, respectively. The structure of OapB exhibits a K-shaped overall architecture wherein its conserved KOW motif and additional unique structural elements of OapB form a bipartite RNA-binding surface that docks to the P13 hairpin and P12.2 helix of OLE RNA. These high-resolution structures elucidate the molecular contacts used by OapB to form a stable RNP complex and explain the high conservation of sequences and structural features at the OapB–OLE RNA-binding interface. These findings provide insight into the role of OapB in the assembly and biological function of OLE RNP complex and can guide the exploration of additional possible OLE RNA-binding interactions present in OapB.

BSGatlas: a unified Bacillus subtilis genome and transcriptome annotation atlas with enhanced information access

A large part of our current understanding of gene regulation in Gram-positive bacteria is based on Bacillus subtilis, as it is one of the most well studied bacterial model systems. The rapid growth in data concerning its molecular and genomic biology is distributed across multiple annotation resources. Consequently, the interpretation of data from further B. subtilis experiments becomes increasingly challenging in both low- and large-scale analyses. Additionally, B. subtilis annotation of structured RNA and non-coding RNA (ncRNA), as well as the operon structure, is still lagging behind the annotation of the coding sequences. To address these challenges, we created the B. subtilis genome atlas, BSGatlas, which integrates and unifies multiple existing annotation resources. Compared to any of the individual resources, the BSGatlas contains twice as many ncRNAs, while improving the positional annotation for 70 % of the ncRNAs. Furthermore, we combined known transcription start and termination sites with lists of known co-transcribed gene sets to create a comprehensive transcript map. The combination with transcription start/termination site annotations resulted in 717 new sets of co-transcribed genes and 5335 untranslated regions (UTRs). In comparison to existing resources, the number of 5' and 3' UTRs increased nearly fivefold, and the number of internal UTRs doubled. The transcript map is organized in 2266 operons, which provides transcriptional annotation for 92 % of all genes in the genome compared to the at most 82 % by previous resources. We predicted an off-target-aware genome-wide library of CRISPR-Cas9 guide RNAs, which we also linked to polycistronic operons. We provide the BSGatlas in multiple forms: as a website (https://rth.dk/resources/bsgatlas/), an annotation hub for display in the UCSC genome browser, supplementary tables and standardized GFF3 format, which can be used in large scale -omics studies. By complementing existing resources, the BSGatlas supports analyses of the B. subtilis genome and its molecular biology with respect to not only non-coding genes but also genome-wide transcriptional relationships of all genes.

A Family of Viral Satellites Manipulates Invading Virus Gene Expression and Can Affect Cholera Toxin Mobilization

Many viruses possess temporally unfolding gene expression patterns aimed at subverting host defenses, commandeering host metabolism, and ultimately producing a large number of progeny virions. High-throughput omics tools, such as RNA sequencing (RNA-seq), have dramatically enhanced the resolution of expression patterns during infection. Less studied have been viral satellites, mobile genomes that parasitize viruses. By performing RNA-seq on infection time courses, we have obtained the first time-resolved transcriptomes for bacteriophage satellites during lytic infection. Specifically, we have acquired transcriptomes for the lytic Vibrio cholerae phage ICP1 and all five known variants of ICP1's parasite, the phage inducible chromosomal island-like elements (PLEs). PLEs rely on ICP1 for both DNA replication and mobilization and abolish production of ICP1 progeny in infected cells. We investigated PLEs' impact on ICP1 gene expression and found that PLEs did not broadly restrict or reduce ICP1 gene expression. A major exception occurred in ICP1's capsid morphogenesis operon, which was downregulated by each of the PLE variants. Surprisingly, PLEs were also found to alter the gene expression of CTXΦ, the integrative phage that encodes cholera toxin and is necessary for virulence of toxigenic V. cholerae One PLE, PLE1, upregulated CTXΦ genes involved in replication and integration and boosted CTXΦ mobility following induction of the SOS response.IMPORTANCE Viral satellites are found in all domains of life and can have profound fitness effects on both the viruses they parasitize and the cells they reside in. In this study, we have acquired the first RNA sequencing (RNA-seq) transcriptomes of viral satellites outside plants, as well as the transcriptome of the phage ICP1, a predominant predator of pandemic Vibrio cholerae Capsid downregulation, previously observed in an unrelated phage satellite, is conserved among phage inducible chromosomal island-like elements (PLEs), suggesting that viral satellites are under strong selective pressure to reduce the capsid expression of their larger host viruses. Despite conserved manipulation of capsid expression, PLEs exhibit divergent effects on CTXΦ transcription and mobility. Our results demonstrate that PLEs can influence both their hosts' resistance to phage and the mobility of virulence-encoding elements, suggesting that PLEs can play a substantial role in shaping Vibrio cholerae evolution.

Exosomal Long Non-coding RNAs: Emerging Players in the Tumor Microenvironment

Recent advances in exosome biology have uncovered a significant role of exosomes in cancer and make them a determining factor in intercellular communication. Exosomes are types of extracellular vesicles that are involved in the communication between cells by exchanging various signaling molecules between the surrounding cells. Among various signaling molecules, long non-coding RNAs (lncRNAs), a type of non-coding RNA having a size of more than 200 nt in length and lacking protein-coding potential, have emerged as crucial regulators of intercellular communication. Tumor-derived exosomes containing various lncRNAs, known as exosomal lncRNAs, reprogram the microenvironment by regulating numerous cellular functions, including the regulation of gene transcription that favors cancer growth and progression, thus significantly determining the biological effects of exosomes. In addition, deregulated expression of lncRNAs is found in various human cancers and serves as a diagnostic biomarker to predict cancer type. The present review discusses the role of exosomal lncRNAs in the crosstalk between tumor cells and the surrounding cells of the microenvironment. Furthermore, we also discuss the involvement of exosomal lncRNAs within the tumor microenvironment in favoring tumor growth, metabolic reprogramming of tumor cells, and tumor-supportive autophagy. Therefore, lncRNAs can be used as a therapeutic target in the treatment of various human cancers.

Transcriptome-wide expression profiling of Sporothrix schenckii yeast and mycelial forms and the establishment of the Sporothrix Genome DataBase

Sporothrix schenckii is a dimorphic fungus existing as mould in the environment and as yeast in the host. The morphological shift between mycelial/yeast phases is crucial for its virulence, but the transcriptional networks implicated in dimorphic transition are still not fully understood. Here, we report the global transcriptomic differences occurring between mould and yeast phases of S. schenckii, including changes in gene expression profiles associated with these distinct cellular phenotypes. Moreover, we also propose a new genome annotation, which reveals a more complex transcriptional architecture than previously assumed. Using RNA-seq, we identified a total of 17 307 genes, of which 11 217 were classified as protein-encoding genes, whereas 6090 were designated as non-coding RNAs (ncRNAs). Approximately ~71 % of all annotated genes were found to overlap and the different-strand overlapping type was the most common. Gene expression analysis revealed that 8795 genes were differentially regulated among yeast and mould forms. Differential gene expression was also observed for antisense ncRNAs overlapping neighbouring protein-encoding genes. The release of transcriptome-wide data and the establishment of the Sporothrix Genome DataBase (http://sporothrixgenomedatabase.unime.it) represent an important milestone for Sporothrix research, because they provide a strong basis for future studies on the molecular pathways involved in numerous biological processes.

Bacillus halodurans OapB forms a high-affinity complex with the P13 region of the noncoding RNA OLE

Noncoding RNAs (ncRNAs) longer than 200 nucleotides are rare in bacteria, likely because bacterial genomes are under strong evolutionary pressures to maintain a small genome size. Of the long ncRNAs unique to bacteria, the OLE (ornate, large, extremophilic) RNA class is among the largest and most structurally complex. OLE RNAs form a ribonucleoprotein (RNP) complex by partnering with at least two proteins, OapA and OapB, that directly bind OLE RNA. The biochemical functions of the OLE RNP complex remain unknown, but are required for proper adaptation to certain environmental stresses, such as cold temperatures, short chain alcohols, and high magnesium concentrations. In the current study, we used electrophoretic mobility shift assays to examine the binding of OLE RNA fragments by OapB and found that OapB recognizes a small subregion of OLE RNA, including stem P13, with a dissociation constant (K_D ) of ∼700 pm Analyses with mutated RNA constructs, and the application of in vitro selection, revealed that strong binding of OLE RNA by OapB requires a stem containing a precisely located single-nucleotide bulge and a GNRA tetraloop. Although the vast majority of bacteria with the ole gene also have the oapB gene, there are many whose genomes contain oapB but lack ole, suggesting that OapB has other RNA partners in some species that might exhibit similar structural features.

The rare lncRNA GOLLD is widespread and structurally conserved among Mycobacterium tRNA arrays

Noncoding RNA (ncRNA) genes produce transcripts involved in a wide range of functions, including catalytic and regulatory functions. Besides, some transcripts have highly complex structures that may impact their activities. Among the largest bacterial ncRNAs, there is the rare GOLLD RNA, which is associated with tRNA genes and supposed to be chromosome- and phage-encoded in specialized groups of bacteria, including those from Lactobacillales and Actinomycetales orders. The only GOLLD structure was inferred from a variety of sequences, including many marine metagenomes. To explore GOLLD RNA in bacterial genomes, we mined the GOLLD gene in thousands of Mycobacterium and virus genomes using Infernal software. We identified this gene in 350 mycobacteria, including megaplasmids, and 39 bacteriophages, mainly in the genomic context of tRNA arrays. Mycobacterium GOLLD genes presented a high diversity and were distributed in three phylogenetic groups: (i) Mycobacterium exclusive; (ii) Mycobacterium and mycobacteriophages; and (iii) mycobacteriophage exclusive. We also determined the GOLLD secondary structure of each group using R2 R software based on GOLLD alignments generated by Infernal software. All GOLLD groups displayed a 3' half conserved structure, including utter E-loops pseudoknots substructures, also shared by non-Mycobacterium GOLLD while the 5' half motif was different among the groups. Here, we showed that the lncRNA GOLLD is widespread in Mycobacterium within tRNA arrays and corroborated the previously predicted GOLLD secondary structure.

Disruption of the OLE ribonucleoprotein complex causes magnesium toxicity in Bacillus halodurans

OLE RNAs represent an unusual class of bacterial noncoding RNAs common in Gram-positive anaerobes. The OLE RNA of the alkaliphile Bacillus halodurans is highly expressed and naturally interacts with at least two RNA-binding proteins called OapA and OapB. The phenotypes of the corresponding knockouts include growth inhibition when exposed to ethanol or other short-chain alcohols or when incubated at modestly reduced temperatures (e.g. 20°C). Intriguingly, the OapA 'PM1' mutant, which carries two amino acid changes to a highly conserved region, yields a dominant-negative phenotype that causes more severe growth defects under these same stress conditions. Herein, we report that the PM1 strain also exhibits extreme sensitivity to elevated Mg²⁺ concentrations, beginning as low as 2 mM. Suppressor mutants predominantly map to genes for aconitate hydratase and isocitrate dehydrogenase, which are expected to alter cellular citrate concentrations. Citrate reduces the severity of the Mg²⁺ toxicity phenotype, but neither the genomic mutations nor the addition of citrate to the medium overcomes ethanol toxicity or temperature sensitivity. These findings reveal that OLE RNA and its protein partners are involved in biochemical responses under several stress conditions, wherein the unusual sensitivity to Mg²⁺ can be independently suppressed by specific genomic mutations.

A conserved RNA structural motif for organizing topology within picornaviral internal ribosome entry sites

Picornaviral IRES elements are essential for initiating the cap-independent viral translation. However, three-dimensional structures of these elements remain elusive. Here, we report a 2.84-Å resolution crystal structure of hepatitis A virus IRES domain V (dV) in complex with a synthetic antibody fragment-a crystallization chaperone. The RNA adopts a three-way junction structure, topologically organized by an adenine-rich stem-loop motif. Despite no obvious sequence homology, the dV architecture shows a striking similarity to a circularly permuted form of encephalomyocarditis virus J-K domain, suggesting a conserved strategy for organizing the domain architecture. Recurrence of the motif led us to use homology modeling tools to compute a 3-dimensional structure of the corresponding domain of foot-and-mouth disease virus, revealing an analogous domain organizing motif. The topological conservation observed among these IRESs and other viral domains implicates a structured three-way junction as an architectural scaffold to pre-organize helical domains for recruiting the translation initiation machinery.

Application of Transcriptomics to Enhance Early Diagnostics of Mycobacterial Infections, with an Emphasis on Mycobacterium avium ssp. paratuberculosis

Mycobacteria cause a wide variety of disease in human and animals. Species that infect ruminants include M. bovis and M. avium ssp. paratuberculosis (MAP). MAP is the causative agent of Johne’s disease in ruminants, which is a chronic granulomatous enteric infection that leads to severe economic losses worldwide. Characteristic of MAP infection is the long, latent phase in which intermittent shedding can take place, while diagnostic tests are unable to reliably detect an infection in this stage. This leads to unnoticed dissemination within herds and the presence of many undetected, silent carriers, which makes the eradication of Johne’s disease difficult. To improve the control of MAP infection, research is aimed at improving early diagnosis. Transcriptomic approaches can be applied to characterize host-pathogen interactions during infection, and to develop novel biomarkers using transcriptional profiles. Studies have focused on the identification of specific RNAs that are expressed in different infection stages, which will assist in the development and clinical implementation of early diagnostic tests.

The Putative Smallest Introns in the Arabidopsis Genome

Most eukaryotic genes contain introns, which are noncoding sequences that are removed during premRNA processing. Introns are usually preserved across evolutionary time. However, the sizes of introns vary greatly. In Arabidopsis, some introns are longer than 10 kilo base pairs (bp) and others are predicted to be shorter than 10 bp. To identify the shortest intron in the genome, we analyzed the predicted introns in annotated version 10 of the Arabidopsis thaliana genome and found 103 predicted introns that are 30 bp or shorter, which make up only 0.08% of all introns in the genome. However, our own bioinformatics and experimental analyses found no evidence for the existence of these predicted introns. The predicted introns of 30–39 bp, 40–49 bp, and 50–59 bp in length are also rare and constitute only 0.07%, 0.2%, and 0.28% of all introns in the genome, respectively. An analysis of 30 predicted introns 31–59 bp long verified two in this range, both of which were 59 bp long. Thus, this study suggests that there is a limit to how small introns in A. thaliana can be, which is useful for the understanding of the evolution and processing of small introns in plants in general.

A second RNA-binding protein is essential for ethanol tolerance provided by the bacterial OLE ribonucleoprotein complex

OLE (ornate, large, extremophilic) RNAs comprise a class of structured noncoding RNAs (ncRNAs) found in many extremophilic bacteria species. OLE RNAs constitute one of the longest and most widespread bacterial ncRNA classes whose major biochemical function remains unknown. In the Gram-positive alkaliphile Bacillus halodurans, OLE RNA is abundant, and localizes to the cell membrane by association with the transmembrane OLE-associated protein called OapA (formerly OAP). These characteristics, along with the well-conserved sequence and structural features of OLE RNAs, suggest that the OLE ribonucleoprotein (RNP) complex performs important biological functions. B. halodurans strains lacking OLE RNA (∆ole) or OapA (∆oapA) are less tolerant of cold (20 °C) and short-chain alcohols (e.g., ethanol). Here, we describe the effects of a mutant OapA (called PM1) that more strongly inhibits growth under cold or ethanol stress compared with strains lacking the oapA gene, even when wild-type OapA is present. This dominant-negative effect of PM1 is reversed by mutations that render OLE RNA nonfunctional. This finding demonstrates that the deleterious PM1 phenotype requires an intact RNP complex, and suggests that the complex has one or more additional undiscovered components. A genetic screen uncovered PM1 phenotype suppressor mutations in the ybzG gene, which codes for a putative RNA-binding protein of unknown biological function. We observe that YbzG protein (also called OapB) selectively binds OLE RNA in vitro, whereas a mutant version of the protein is not observed to bind OLE RNA. Thus, YbzG/OapB is an important component of the functional OLE RNP complex in B. halodurans.