Drawing and editing the secondary structure(s) of RNA

A 36-base hairpin within lncRNA DRAIC, which is modulated by alternative splicing, interacts with the IKKα coiled-coil domain and inhibits NF-κB and tumor cell phenotypes

A tumor-suppressive long noncoding RNA (lncRNA) DRAIC (down-regulated RNA in cancers) inhibits NF-κB activity and physically interacts with IKKα, a kinase component of the IKK complex, in several cancer types. Here we explore the precise molecular mechanism involved in this interaction and suppression. Using SHAPE-MaP, we identified a 36-nucleotide hairpin (A+B) within DRAIC that is necessary and sufficient for its anti-oncogenic function. RNA immunoprecipitation (RIP) and Electrophoretic mobility shift assays (EMSA) confirmed this hairpin physically interacts with the coiled coil domain of IKKα. A+B RNA has a high binding affinity (KD ∼1-7 nM) to the coiled-coil domain of IKKα. The binding of A+B disrupts the dimerization of NEMO and IKKα coiled-coil domains, a critical step for IKK action. Consistent with this, A+B inhibits the phosphorylation of the NF-κB inhibitor IκBα and suppresses NF-κB activity. Publicly available tumor RNAseq data revealed that alternative splicing modulates the presence of this critical hairpin: the inclusion of exon 4a (encoding one side of the A+B hairpin) in lung tumors correlates with reduced NF-κB activity. By demonstrating that the A+B hairpin is both necessary and sufficient to inhibit IKK and oncogenic phenotypes, this study underscores the centrality of IKKα interaction and NF-κB inhibition in DRAIC-mediated cancer suppression and indicates that the activity of this lncRNA is regulated by alternative splicing. This study also reveals the first example of a short RNA disrupting coiled-coil dimerization, offering a new approach to disrupt such dimerization in cancer biology.

R2DT: a comprehensive platform for visualizing RNA secondary structure

RNA secondary (2D) structure visualization is an essential tool for understanding RNA function. R2DT is a software package designed to visualize RNA 2D structures in consistent, recognizable, and reproducible layouts. The latest release, R2DT 2.0, introduces multiple significant features, including the ability to display position-specific information, such as single nucleotide polymorphisms or SHAPE reactivities. It also offers a new template-free mode allowing visualization of RNAs without pre-existing templates, alongside a constrained folding mode and support for animated visualizations. Users can interactively modify R2DT diagrams, either manually or using natural language prompts, to generate new templates or create publication-quality images. Additionally, R2DT features faster performance, an expanded template library, and a growing collection of compatible tools and utilities. Already integrated into multiple biological databases, R2DT has evolved into a comprehensive platform for RNA 2D visualization, accessible at https://r2dt.bio.

A 36-base hairpin within lncRNA DRAIC , modulated by alternative splicing, interacts with the IKKα coiled-coil domain and inhibits NF-κB and tumor cell phenotypes

The 1.7 kb DRAIC long noncoding RNA inhibits tumor growth, inhibits cancer cell invasion, migration, colony formation and interacts with IKK (IκB kinase) subunits, inhibiting the phosphorylation and degradation of the NF-κB inhibitor, IκB, to suppress the activation of NF-κB. Whether these activities are all linked is unclear. We used SHAPE-MaP to obtain the secondary structure of the lncRNA to perform structure-functions studies which identified the minimal region of DRAIC necessary for repressing NF-κB. A 36-nucleotide hairpin (A+B) within DRAIC inhibits NF-κB, inhibits IκB phosphorylation and binds specifically with the IKKα coiled- coil domain with a very high affinity: K D of ∼1-5 nM. This interaction weakens the dimerization of of the coiled coil domains of two IKK subunits, a dimerization that is indispensable for IKK activity. A+B is sufficient and necessary to inhibit the oncogenic phenotypes in multiple cancer cell-lines, demonstrating that interaction with IKK and inhibition of NF-κB is key for cancer suppression by DRAIC. Presence of this critical hairpin is modulated by alternative splicing the extends exon 4 to exon 4a of DRAIC and the expression of exon 4a in lung tumors is associated with low NF-κB activity. This is also the first demonstration that a short RNA can disrupt coiled- coil dimerization.

R2DT: A COMPREHENSIVE PLATFORM FOR VISUALISING RNA SECONDARY STRUCTURE

RNA secondary (2D) structure visualisation is an essential tool for understanding RNA function. R2DT is a software package designed to visualise RNA 2D structures in consistent, recognisable, and reproducible layouts. The latest release, R2DT 2.0, introduces multiple significant features, including the ability to display position-specific information, such as single nucleotide polymorphisms (SNPs) or SHAPE reactivities. It also offers a new template-free mode allowing visualisation of RNAs without pre-existing templates, alongside a constrained folding mode and support for animated visualisations. Users can interactively modify R2DT diagrams, either manually or using natural language prompts, to generate new templates or create publication-quality images. Additionally, R2DT features faster performance, an expanded template library, and a growing collection of compatible tools and utilities. Already integrated into multiple biological databases, R2DT has evolved into a comprehensive platform for RNA 2D visualisation, accessible at https://r2dt.bio.

RiboVision2: A Web Server for Advanced Visualization of Ribosomal RNAs

RiboVision2 is a web server designed to visualize phylogenetic, structural, and evolutionary properties of ribosomal RNAs simultaneously at the levels of primary, secondary, and three-dimensional structure and in the context of full ribosomal complexes. RiboVision2 instantly computes and displays a broad variety of data; it has no login requirements, is open-source, free for all users, and available at https://ribovision2.chemistry.gatech.edu.

RNAcanvas: interactive drawing and exploration of nucleic acid structures

Two-dimensional drawing of nucleic acid structures, particularly RNA structures, is fundamental to the communication of nucleic acids research. However, manually drawing structures is laborious and infeasible for structures thousands of nucleotides long. RNAcanvas automatically arranges residues into strictly shaped stems and loops while providing robust interactive editing features, including click-and-drag layout adjustment. Drawn elements are highly customizable in a point-and-click manner, including colours, fonts, size and shading, flexible numbering, and outlining of bases. Tertiary interactions can be drawn as draggable, curved lines. Leontis-Westhof notation for depicting non-canonical base-pairs is fully supported, as well as text labels for structural features (e.g. hairpins). RNAcanvas also has many unique features and performance optimizations for large structures that cannot be correctly predicted and require manual refinement based on the researcher's own analyses and expertise. To this end, RNAcanvas has point-and-click structure editing with real-time highlighting of complementary sequences and motif search functionality, novel features that greatly aid in the identification of putative long-range tertiary interactions, de novo analysis of local structures, and phylogenetic comparisons. For ease in producing publication quality figures, drawings can be exported in both SVG and PowerPoint formats. URL: https://rnacanvas.app.

Predicting the Structure of a Viroid : Structure, Structure Distribution, Consensus Structure, and Structure Drawing

Viroids are small non-coding RNAs that require a special sequence and structure to be replicated and transported by the host machinery. Many of these features can be predicted and later experimentally verified. Here, we will present workflows to predict viroid structures and draw the predicted structures in a pleasing and descriptive way using recently developed software.

R2DT is a framework for predicting and visualising RNA secondary structure using templates

Non-coding RNAs (ncRNA) are essential for all life, and their functions often depend on their secondary (2D) and tertiary structure. Despite the abundance of software for the visualisation of ncRNAs, few automatically generate consistent and recognisable 2D layouts, which makes it challenging for users to construct, compare and analyse structures. Here, we present R2DT, a method for predicting and visualising a wide range of RNA structures in standardised layouts. R2DT is based on a library of 3,647 templates representing the majority of known structured RNAs. R2DT has been applied to ncRNA sequences from the RNAcentral database and produced >13 million diagrams, creating the world's largest RNA 2D structure dataset. The software is amenable to community expansion, and is freely available at https://github.com/rnacentral/R2DT and a web server is found at https://rnacentral.org/r2dt .

Evolution of the genetic code in the mitochondria of Labyrinthulea (Stramenopiles)

Mitochondrial translation often exhibits departures from the standard genetic code, but the full spectrum of these changes has certainly not yet been described and the molecular mechanisms behind the changes in codon meaning are rarely studied. Here we report a detailed analysis of the mitochondrial genetic code in the stramenopile group Labyrinthulea (Labyrinthulomycetes) and their relatives. In the genus Aplanochytrium, UAG is not a termination codon but encodes tyrosine, in contrast to the unaffected meaning of the UAA codon. This change is evolutionarily independent of the reassignment of both UAG and UAA as tyrosine codons recently reported from two uncultivated labyrinthuleans (S2 and S4), which we show are not thraustochytrids as proposed before, but represent the clade LAB14 previously recognised in environmental 18S rRNA gene surveys. We provide rigorous evidence that the UUA codon in the mitochondria of all labyrinthuleans serves as a termination codon instead of encoding leucine, and propose that a sense-to-stop reassignment has also affected the AGG and AGA codons in the LAB14 clade. The distribution of the different forms of sense-to-stop and stop-to-sense reassignments correlates with specific modifications of the mitochondrial release factor mtRF2a in different subsets of labyrinthuleans, and with the unprecedented loss of mtRF1a in Aplanochytrium and perhaps also in the LAB14 clade, pointing towards a possible mechanistic basis of the code changes observed. Curiously, we show that labyrinthulean mitochondria also exhibit a sense-to-sense codon reassignment, manifested as AUA encoding methionine instead of isoleucine. Furthermore, we show that this change evolved independently in the uncultivated stramenopile lineage MAST8b, together with the reassignment of the AGR codons from arginine to serine. Altogether, our study has uncovered novel variants of the mitochondrial genetic code and previously unknown modifications of the mitochondrial translation machinery, further enriching our understanding of the rules governing the evolution of one of the central molecular process in the cell.

Evolution and Unprecedented Variants of the Mitochondrial Genetic Code in a Lineage of Green Algae

Mitochondria of diverse eukaryotes have evolved various departures from the standard genetic code, but the breadth of possible modifications and their phylogenetic distribution are known only incompletely. Furthermore, it is possible that some codon reassignments in previously sequenced mitogenomes have been missed, resulting in inaccurate protein sequences in databases. Here we show, considering the distribution of codons at conserved amino acid positions in mitogenome-encoded proteins, that mitochondria of the green algal order Sphaeropleales exhibit a diversity of codon reassignments, including previously missed ones and some that are unprecedented in any translation system examined so far, necessitating redefinition of existing translation tables and creating at least seven new ones. We resolve a previous controversy concerning the meaning the UAG codon in Hydrodictyaceae, which beyond any doubt encodes alanine. We further demonstrate that AGG, sometimes together with AGA, encodes alanine instead of arginine in diverse sphaeroplealeans. Further newly detected changes include Arg-to-Met reassignment of the AGG codon and Arg-to-Leu reassignment of the CGG codon in particular species. Analysis of tRNAs specified by sphaeroplealean mitogenomes provides direct support for and molecular underpinning of the proposed reassignments. Furthermore, we point to unique mutations in the mitochondrial release factor mtRF1a that correlate with changes in the use of termination codons in Sphaeropleales, including the two independent stop-to-sense UAG reassignments, the reintroduction of UGA in some Scenedesmaceae, and the sense-to-stop reassignment of UCA widespread in the group. Codon disappearance seems to be the main drive of the dynamic evolution of the mitochondrial genetic code in Sphaeropleales.

jViz.RNA 4.0-Visualizing pseudoknots and RNA editing employing compressed tree graphs

Previously, we have introduced an improved version of jViz.RNA which enabled faster and more stable RNA visualization by employing compressed tree graphs. However, the new RNA representation and visualization method required a sophisticated mechanism of pseudoknot visualization. In this work, we present our novel pseudoknot classification and implementation of pseudoknot visualization in the context of the new RNA graph model. We then compare our approach with other RNA visualization software, and demonstrate jViz.RNA 4.0's benefits compared to other software. Additionally, we introduce interactive editing functionality into jViz.RNA and demonstrate its benefits in exploring and building RNA structures. The results presented highlight the new high degree of utility jViz.RNA 4.0 now offers. Users are now able to visualize pseudoknotted RNA, manipulate the resulting automatic layouts to suit their individual needs, and change both positioning and connectivity of the RNA molecules examined. Care was taken to limit overlap between structural elements, particularly in the case of pseudoknots to ensure an intuitive and informative layout of the final RNA structure. Availability: The software is freely available at: https://jviz.cs.sfu.ca/.

RNApuzzler: efficient outerplanar drawing of RNA-secondary structures

Motivation

RNA secondary structure is a useful representation for studying the function of RNA, which captures most of the free energy of RNA folding. Using empirically determined energy parameters, secondary structures of nucleic acids can be efficiently computed by recursive algorithms. Several software packages supporting this task are readily available. As RNA secondary structures are outerplanar graphs, they can be drawn without intersection in the plane. Interpretation by the practitioner is eased when these drawings conform to a series of additional constraints beyond outerplanarity. These constraints are the reason why RNA drawing is difficult. Many RNA drawing algorithms therefore do not always produce intersection-free (outerplanar) drawings.

Results

To remedy this shortcoming we propose here the RNApuzzler algorithm which is guaranteed to produce intersection-free drawings. It is based on a drawing algorithm respecting constraints based on nucleotide distances (RNAturtle). We investigate relaxations of these constraints allowing for intersection-free drawings. Based on these relaxations, we implemented a fully automated, simple, and robust algorithm that produces aesthetic drawings adhering to previously established guidelines. We tested our algorithm using the RFAM database and found that we can compute intersection-free drawings of all RNAs therein efficiently.

Availability and implementation

The software can be accessed freely at: https://github.com/dwiegreffe/RNApuzzler.

Supplementary information

Supplementary data are available at Bioinformatics online.

TRAVeLer: a tool for template-based RNA secondary structure visualization

BACKGROUND

Visualization of RNA secondary structures is a complex task, and, especially in the case of large RNA structures where the expected layout is largely habitual, the existing visualization tools often fail to produce suitable visualizations. This led us to the idea to use existing layouts as templates for the visualization of new RNAs similarly to how templates are used in homology-based structure prediction.

RESULTS

This article introduces Traveler, a software tool enabling visualization of a target RNA secondary structure using an existing layout of a sufficiently similar RNA structure as a template. Traveler is based on an algorithm which converts the target and template structures into corresponding tree representations and utilizes tree edit distance coupled with layout modification operations to transform the template layout into the target one. Traveler thus accepts a pair of secondary structures and a template layout and outputs a layout for the target structure.

CONCLUSIONS

Traveler is a command-line open source tool able to quickly generate layouts for even the largest RNA structures in the presence of a sufficiently similar layout. It is available at http://github.com/davidhoksza/traveler .

Amalga-like virus infecting Antonospora locustae, a microsporidian pathogen of grasshoppers, plus related viruses associated with other arthropods

A previously reported Expressed Sequence Tag (EST) library from spores of microsporidian Antonospora locustae includes a number of clones with sequence similarities to plant amalgaviruses. Reexamining the sequence accessions from that library, we found additional such clones, contributing to a 3247-nt contig that approximates the length of an amalga-like virus genome. Using A. locustae spores stored from that previous study, and new ones obtained from the same source, we newly visualized the putative dsRNA genome of this virus and obtained amplicons yielding a 3387-nt complete genome sequence. Phylogenetic analyses suggested it as prototype strain of a new genus in family Amalgaviridae. The genome contains two partially overlapping long ORFs, with downstream ORF2 in the +1 frame relative to ORF1 and a proposed motif for +1 ribosomal frameshifting in the region of overlap. Subsequent database searches using the predicted fusion protein sequence of this new amalga-like virus identified related sequences in the transcriptome of a basal hexapod, the springtail species Tetrodontophora bielanensis. We speculate that this second new amalga-like virus (contig length, 3475 nt) likely also derived from a microsporidian, or related organism, which was associated with the springtail specimens at the time of sampling for transcriptome analysis. Other findings of interest include evidence that the ORF1 translation products of these two new amalga-like viruses contain a central region of predicted α-helical coiled coil, as recently reported for plant amalgaviruses, and transcriptome-based evidence for another new amalga-like virus in the transcriptome of another basal hexapod, the two-pronged bristletail species Campodea augens.

Revisiting the structure/function relationships of H/ACA(-like) RNAs: a unified model for Euryarchaea and Crenarchaea

A structural and functional classification of H/ACA and H/ACA-like motifs is obtained from the analysis of the H/ACA guide RNAs which have been identified previously in the genomes of Euryarchaea (Pyrococcus) and Crenarchaea (Pyrobaculum). A unified structure/function model is proposed based on the common structural determinants shared by H/ACA and H/ACA-like motifs in both Euryarchaea and Crenarchaea. Using a computational approach, structural and energetic rules for the guide:target RNA-RNA interactions are derived from structural and functional data on the H/ACA RNP particles. H/ACA(-like) motifs found in Pyrococcus are evaluated through the classification and their biological relevance is discussed. Extra-ribosomal targets found in both Pyrococcus and Pyrobaculum might support the hypothesis of a gene regulation mediated by H/ACA(-like) guide RNAs in archaea.