Natural circularly permuted group II introns in bacteria produce RNA circles

ARTEMIS: a method for topology-independent superposition of RNA 3D structures and structure-based sequence alignment

Non-coding RNAs play a major role in diverse processes in living cells with their sequence and spatial structure serving as the principal determinants of their function. Superposition of RNA 3D structures is the most accurate method for comparative analysis of RNA molecules and for inferring structure-based sequence alignments. Topology-independent superposition is particularly relevant, as evidenced by structurally similar RNAs with sequence permutations such as tRNA and Y RNA. To date, state-of-the-art methods for RNA 3D structure superposition rely on intricate heuristics, and the potential for topology-independent superposition has not been exhausted. Recently, we introduced the ARTEM method for unrestrained pairwise superposition of RNA 3D modules and now we developed it further to solve the global RNA 3D structure alignment problem. Our new tool ARTEMIS significantly outperforms state-of-the-art tools in both sequentially-ordered and topology-independent RNA 3D structure superposition. Using ARTEMIS we discovered a helical packing motif to be preserved within different backbone topology contexts across various non-coding RNAs, including multiple ribozymes and riboswitches. We anticipate that ARTEMIS will be essential for elucidating the landscape of RNA 3D folds and motifs featuring sequence permutations that thus far remained unexplored due to limitations in previous computational approaches.

Dynamic conformation: Marching toward circular RNA function and application

Circular RNA is a group of covalently closed, single-stranded transcripts with unique biogenesis, stability, and conformation that play distinct roles in modulating cellular functions and also possess a great potential for developing circular RNA-based therapies. Importantly, due to its circular conformation, circular RNA generates distinct intramolecular base pairing that is different from the linear transcript. In this perspective, we review how circular RNA conformation can affect its turnover and modes of action, as well as what factors can modulate circular RNA conformation. We also discuss how understanding circular RNA conformation can facilitate learning about their functions as well as the remaining technological issues to further address their conformation. These efforts will ultimately inform the design of circular RNA-based platforms for biomedical applications.

Temperature-Dependent tRNA Modifications in Bacillales

Transfer RNA (tRNA) modifications are essential for the temperature adaptation of thermophilic and psychrophilic organisms as they control the rigidity and flexibility of transcripts. To further understand how specific tRNA modifications are adjusted to maintain functionality in response to temperature fluctuations, we investigated whether tRNA modifications represent an adaptation of bacteria to different growth temperatures (minimal, optimal, and maximal), focusing on closely related psychrophilic (P. halocryophilus and E. sibiricum), mesophilic (B. subtilis), and thermophilic (G. stearothermophilus) Bacillales. Utilizing an RNA sequencing approach combined with chemical pre-treatment of tRNA samples, we systematically profiled dihydrouridine (D), 4-thiouridine (s4U), 7-methyl-guanosine (m7G), and pseudouridine (Ψ) modifications at single-nucleotide resolution. Despite their close relationship, each bacterium exhibited a unique tRNA modification profile. Our findings revealed increased tRNA modifications in the thermophilic bacterium at its optimal growth temperature, particularly showing elevated levels of s4U8 and Ψ55 modifications compared to non-thermophilic bacteria, indicating a temperature-dependent regulation that may contribute to thermotolerance. Furthermore, we observed higher levels of D modifications in psychrophilic and mesophilic bacteria, indicating an adaptive strategy for cold environments by enhancing local flexibility in tRNAs. Our method demonstrated high effectiveness in identifying tRNA modifications compared to an established tool, highlighting its potential for precise tRNA profiling studies.

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.

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.

A linear and circular dual-conformation noncoding RNA involved in oxidative stress tolerance in Bacillus altitudinis

Circular RNAs have been extensively studied in eukaryotes, but their presence and/or biological functionality in bacteria are unclear. Here, we show that a regulatory noncoding RNA (DucS) exists in both linear and circular conformation in Bacillus altitudinis. The linear forms promote B. altitudinis tolerance to H2O2 stress, partly through increased translation of a stress-responsive gene, htrA. The 3' end sequences of the linear forms are crucial for RNA circularization, and formation of circular forms can decrease the levels of the regulatory linear cognates. Bioinformatic analysis of available RNA-seq datasets from 30 bacterial species revealed multiple circular RNA candidates, distinct from DucS, for all the examined species. Experiments testing for the presence of selected circular RNA candidates in four species successfully validated 7 out of 9 candidates from B. altitudinis and 4 out of 5 candidates from Bacillus paralicheniformis; However, none of the candidates tested for Bacillus subtilis and Escherichia coli were detected. Our work identifies a dual-conformation regulatory RNA in B. altitutidinis, and indicates that circular RNAs exist in diverse bacteria. However, circularization of specific RNAs does not seem to be conserved across species, and the circularization mechanisms and biological functionality of the circular forms remain unclear.

The Fate and Functionality of Alien tRNA Fragments in Culturing Medium and Cells of Escherichia coli

Numerous observations have supported the idea that various types of noncoding RNAs, including tRNA fragments (tRFs), are involved in communications between the host and its microbial community. The possibility of using their signaling function has stimulated the study of secreted RNAs, potentially involved in the interspecies interaction of bacteria. This work aimed at identifying such RNAs and characterizing their maturation during transport. We applied an approach that allowed us to detect oligoribonucleotides secreted by Prevotella copri (Segatella copri) or Rhodospirillum rubrum inside Escherichia coli cells. Four tRFs imported by E. coli cells co-cultured with these bacteria were obtained via chemical synthesis, and all of them affected the growth of E. coli. Their successive modifications in the culture medium and recipient cells were studied by high-throughput cDNA sequencing. Instead of the expected accidental exonucleolysis, in the milieu, we observed nonrandom cleavage by endonucleases continued in recipient cells. We also found intramolecular rearrangements of synthetic oligonucleotides, which may be considered traces of intermediate RNA circular isomerization. Using custom software, we estimated the frequency of such events in transcriptomes and secretomes of E. coli and observed surprising reproducibility in positions of such rare events, assuming the functionality of ring isoforms or their permuted derivatives in bacteria.

CircularSTAR3D: a stack-based RNA 3D structural alignment tool for circular matching

The functions of non-coding RNAs usually depend on their 3D structures. Therefore, comparing RNA 3D structures is critical in analyzing their functions. We noticed an interesting phenomenon that two non-coding RNAs may share similar substructures when rotating their sequence order. To the best of our knowledge, no existing RNA 3D structural alignment tools can detect this type of matching. In this article, we defined the RNA 3D structure circular matching problem and developed a software tool named CircularSTAR3D to solve this problem. CircularSTAR3D first uses the conserved stacks (consecutive base pairs with similar 3D structures) in the input RNAs to identify the circular matched internal loops and multiloops. Then it performs a local extension iteratively to obtain the whole circular matched substructures. The computational experiments conducted on a non-redundant RNA structure dataset show that circular matching is ubiquitous. Furthermore, we demonstrated the utility of CircularSTAR3D by detecting the conserved substructures missed by regular alignment tools, including structural motifs and conserved structures between riboswitches and ribozymes from different classes. We anticipate CircularSTAR3D to be a valuable supplement to the existing RNA 3D structural analysis techniques.

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.