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

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.

Genome-Wide Computational Prediction and Analysis of Noncoding RNAs in Oleidesulfovibrio alaskensis G20

Noncoding RNAs (ncRNAs) play key roles in the regulation of important pathways, including cellular growth, stress management, signaling, and biofilm formation. Sulfate-reducing bacteria (SRB) contribute to huge economic losses causing microbial-induced corrosion through biofilms on metal surfaces. To effectively combat the challenges posed by SRB, it is essential to understand their molecular mechanisms of biofilm formation. This study aimed to identify ncRNAs in the genome of a model SRB, Oleidesulfovibrio alaskensis G20 (OA G20). Three in silico approaches revealed genome-wide distribution of 37 ncRNAs excluding tRNAs in the OA G20. These ncRNAs belonged to 18 different Rfam families. This study identified riboswitches, sRNAs, RNP, and SRP. The analysis revealed that these ncRNAs could play key roles in the regulation of several pathways of biosynthesis and transport involved in biofilm formation by OA G20. Three sRNAs, Pseudomonas P10, Hammerhead type II, and sX4, which were found in OA G20, are rare and their roles have not been determined in SRB. These results suggest that applying various computational methods could enrich the results and lead to the discovery of additional novel ncRNAs, which could lead to understanding the "rules of life of OA G20" during biofilm formation.

Deciphering phylogenetic relationships of and delimiting species boundaries within the controversial ciliate genus Conchophthirus using an integrative morpho-evo approach

The phylum Ciliophora (ciliates) comprises about 2600 symbiotic and over 5500 free-living species. The inclusion of symbiotic ciliates in phylogenetic analyses often challenges traditional classification frameworks due to their morphological adaptions to the symbiotic lifestyle. Conchophthirus is such a controversial obligate endocommensal genus whose affinities to other symbiotic and free-living scuticociliates are still poorly understood. Using uni- and multivariate morphometrics as well as 2D-based molecular and phylogenetic analyses, we attempted to test for the monophyly of Conchophthirus, study the boundaries of Conchophthirus species isolated from various bivalves at mesoscale, and reveal the phylogenetic relationships of Conchophthirus to other scuticociliates. Multidimensional analyses of morphometric and cell geometric data generated the same homogenous clusters, as did phylogenetic analyses based on 144 new sequences of two mitochondrial and five nuclear molecular markers. Conchophthirus is not closely related to 'core' scuticociliates represented by the orders Pleuronematida and Philasterida, as assumed in the past using morphological data. Nuclear and mitochondrial markers consistently showed the free-living Dexiotricha and the mouthless endosymbiotic Haptophrya to be the nearest relatives of Conchophthirus. These three highly morphologically and ecologically dissimilar genera represent an orphan clade from the early radiation of scuticociliates in molecular phylogenies.

First molecular evidence of hybridization in endosymbiotic ciliates (Protista, Ciliophora)

Hybridization is an important evolutionary process that can fuel diversification via formation of hybrid species or can lead to fusion of previously separated lineages by forming highly diverse species complexes. We provide here the first molecular evidence of hybridization in wild populations of ciliates, a highly diverse group of free-living and symbiotic eukaryotic microbes. The impact of hybridization was studied on the model of Plagiotoma, an obligate endosymbiont of the digestive tube of earthworms, using split decomposition analyses and species networks, 2D modeling of the nuclear rRNA molecules and compensatory base change analyses as well as multidimensional morphometrics. Gene flow slowed down and eventually hampered the diversification of Lumbricus-dwelling plagiotomids, which collapsed into a single highly variable biological entity, the P. lumbrici complex. Disruption of the species boundaries was suggested also by the continuum of morphological variability in the phenotypic space. On the other hand, hybridization conspicuously increased diversity in the nuclear rDNA cistron and somewhat weakened the host structural specificity of the P. lumbrici complex, whose members colonize a variety of phylogenetically closely related anecic and epigeic earthworms. By contrast, another recorded species, P. aporrectodeae sp. n., showed no signs of introgression, no variability in the rDNA cistron, and very high host specificity. These contrasting eco-evolutionary patterns indicate that hybridization might decrease the alpha-diversity by dissolving species boundaries, weaken the structural host specificity by broadening ecological amplitudes, and increase the nuclear rDNA variability by overcoming concerted evolution within the P. lumbrici species complex.

Nuclear and Mitochondrial SSU rRNA Genes Reveal Hidden Diversity of Haptophrya Endosymbionts in Freshwater Planarians and Challenge Their Traditional Classification in Astomatia

Like many other aquatic animals, freshwater planarians have also become partners of symbiotic ciliates from the class Oligohymenophorea. In the present study, we explored the hidden diversity and addressed the questionable systematic position of mouthless obligatory gut endosymbionts of freshwater planarians, using the nuclear and mitochondrial SSU rRNA genes. Although all isolated ciliates morphologically corresponded to a single species, molecular analyses suggested the existence of three genetically distinct entities: Haptophrya planariarum, Haptophrya dugesiarum nov. spec., and Haptophrya schmidtearum nov. spec. The two former species share the same planarian host, which indicates a speciation model involving one duplication event without host switching. Such a diversification pattern was recognized also in astome ciliates inhabiting megascolecid and glossoscolecid earthworms. The present multi-gene phylogenies along with the secondary structure of the mitochondrial 16S rRNA molecule, however, challenge the traditional classification of Haptophrya within the subclass Astomatia. Haptophrya very likely evolved from an orphan scuticociliate lineage by the loss of oral apparatus and by the transformation of the thigmotactic field into an adhesive sucker. Since astomy evolved multiple times independently within the Oligohymenophorea, the loss of cell mouth cannot be used as a sole argument for the assignment of Haptophrya to the Astomatia anymore.

RiboDraw: semiautomated two-dimensional drawing of RNA tertiary structure diagrams

Publishing, discussing, envisioning, modeling, designing and experimentally determining RNA three-dimensional (3D) structures involve preparation of two-dimensional (2D) drawings that depict critical functional features of the subject molecules, such as noncanonical base pairs and protein contacts. Here, we describe RiboDraw, new software for crafting these drawings. We illustrate the features of RiboDraw by applying it to several RNAs, including the Escherichia coli tRNA-Phe, the P4-P6 domain of Tetrahymena ribozyme, a -1 ribosomal frameshift stimulation element from beet western yellows virus and the 5' untranslated region of SARS-CoV-2. We show secondary structure diagrams of the 23S and 16S subunits of the E. coli ribosome that reflect noncanonical base pairs, ribosomal proteins and structural motifs, and that convey the relative positions of these critical features in 3D space. This software is a MATLAB package freely available at https://github.com/DasLab/RiboDraw.

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 .

RNAcentral 2021: secondary structure integration, improved sequence search and new member databases

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences that provides a single access point to 44 RNA resources and >18 million ncRNA sequences from a wide range of organisms and RNA types. RNAcentral now also includes secondary (2D) structure information for >13 million sequences, making RNAcentral the world's largest RNA 2D structure database. The 2D diagrams are displayed using R2DT, a new 2D structure visualization method that uses consistent, reproducible and recognizable layouts for related RNAs. The sequence similarity search has been updated with a faster interface featuring facets for filtering search results by RNA type, organism, source database or any keyword. This sequence search tool is available as a reusable web component, and has been integrated into several RNAcentral member databases, including Rfam, miRBase and snoDB. To allow for a more fine-grained assignment of RNA types and subtypes, all RNAcentral sequences have been annotated with Sequence Ontology terms. The RNAcentral database continues to grow and provide a central data resource for the RNA community. RNAcentral is freely available at https://rnacentral.org.

T-psi-C: user friendly database of tRNA sequences and structures

tRNAs have been widely studied for their role as genetic code decoders in the ribosome during translation, but have recently received new attention due to the discovery of novel roles beyond decoding, often in connection with human diseases. Yet, existing tRNA databases have not been updated for more than a decade, so they do not contain this new functional information and have not kept pace with the rate of discovery in this field. Therefore, a regularly updated database that contains information about newly discovered characteristics of tRNA molecules and can be regularly updated is strongly needed. Here, we report the creation of the T-psi-C database (http://tpsic.igcz.poznan.pl), an up-to-date collection of tRNA sequences that contains data obtained from high-throughput tRNA sequencing, e.g. all isoacceptors and isodecoders for human HEK293 cells. This database also contains 3D tRNA structures obtained from Protein Data Bank and generated using homology modeling. The T-psi-C database can be continuously updated by any member of the scientific community, and contains its own application programming interface (API), which allows users to retrieve or upload data in JSON format. Altogether, T-psi-C is user-friendly, easy to develop and an up-to-date source of knowledge about tRNAs.

RNAcentral: a hub of information for non-coding RNA sequences

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences, collating information on ncRNA sequences of all types from a broad range of organisms. We have recently added a new genome mapping pipeline that identifies genomic locations for ncRNA sequences in 296 species. We have also added several new types of functional annotations, such as tRNA secondary structures, Gene Ontology annotations, and miRNA-target interactions. A new quality control mechanism based on Rfam family assignments identifies potential contamination, incomplete sequences, and more. The RNAcentral database has become a vital component of many workflows in the RNA community, serving as both the primary source of sequence data for academic and commercial groups, as well as a source of stable accessions for the annotation of genomic and functional features. These examples are facilitated by an improved RNAcentral web interface, which features an updated genome browser, a new sequence feature viewer, and improved text search functionality. RNAcentral is freely available at https://rnacentral.org.

rPredictorDB: a predictive database of individual secondary structures of RNAs and their formatted plots

Secondary data structure of RNA molecules provides insights into the identity and function of RNAs. With RNAs readily sequenced, the question of their structural characterization is increasingly important. However, RNA structure is difficult to acquire. Its experimental identification is extremely technically demanding, while computational prediction is not accurate enough, especially for large structures of long sequences. We address this difficult situation with rPredictorDB, a predictive database of RNA secondary structures that aims to form a middle ground between experimentally identified structures in PDB and predicted consensus secondary structures in Rfam. The database contains individual secondary structures predicted using a tool for template-based prediction of RNA secondary structure for the homologs of the RNA families with at least one homolog with experimentally solved structure. Experimentally identified structures are used as the structural templates and thus the prediction has higher reliability than de novo predictions in Rfam. The sequences are downloaded from public resources. So far rPredictorDB covers 7365 RNAs with their secondary structures. Plots of the secondary structures use the Traveler package for readable display of RNAs with long sequences and complex structures, such as ribosomal RNAs. The RNAs in the output of rPredictorDB are extensively annotated and can be viewed, browsed, searched and downloaded according to taxonomic, sequence and structure data. Additionally, structure of user-provided sequences can be predicted using the templates stored in rPredictorDB.

RNA2Drawer: geometrically strict drawing of nucleic acid structures with graphical structure editing and highlighting of complementary subsequences

RNA structure prediction programs remain imperfect and many substructures are still identified by manual exploration, which is most efficiently conducted within an RNA structure drawing program. However, most nucleic acid structure drawing programs have limited capability for structure modification (i.e., breaking and forming new bonds between bases), often requiring that the structure notation be textually edited. RNA2Drawer was developed to allow for graphical structure editing while maintaining the geometry of a drawing (e.g., ellipsoid loops, stems with evenly stacked base pairs) throughout structural changes and manual adjustments to the layout by the user. In addition, the program allows for annotations such as colouring and circling of bases and drawing of tertiary interactions (e.g., pseudoknots). RNA2Drawer can also draw commonly desired elements such as an optionally flattened outermost loop and assists structure editing by automatically highlighting complementary subsequences, which facilitates the discovery of potentially new and alternative pairings, particularly tertiary pairings over long-distances, which are biologically critical in the genomes of many RNA viruses and cannot be accurately predicted by current structure prediction programs. Additionally, RNA2Drawer outputs drawings either as PNG files, or as PPTX and SVG files, such that every object of a drawing (e.g., bases, bonds) is an individual PPTX or SVG object, allowing for further manipulation in Microsoft PowerPoint or a vector graphics editor such as Adobe Illustrator. PowerPoint is the standard for presentations and is often used to create figures for publications, and RNA2Drawer is the first program to export drawings as PPTX files.

RiboSketch: versatile visualization of multi-stranded RNA and DNA secondary structure

Summary

Creating clear, visually pleasing 2D depictions of RNA and DNA strands and their interactions is important to facilitate and communicate insights related to nucleic acid structure. Here we present RiboSketch, a secondary structure image production application that enables the visualization of multistranded structures via layout algorithms, comprehensive editing capabilities, and a multitude of simulation modes. These interactive features allow RiboSketch to create publication quality diagrams for structures with a wide range of composition, size and complexity. The program may be run in any web browser without the need for installation, or as a standalone Java application.

Availability and implementation

https://rnastructure.cancer.gov/ribosketch.

Pairwise visual comparison of small RNA secondary structures with base pair probabilities

BACKGROUND

Predicted RNA secondary structures are typically visualized using dot-plots for base pair binding probabilities and planar graphs for unique structures, such as the minimum free energy structure. These are however difficult to analyze simultaneously.

RESULTS

This work introduces a compact unified view of the most stable conformation of an RNA secondary structure and its base pair probabilities, which is called the Circular Secondary Structure Base Pairs Probabilities Plot (CS2BP2-Plot). Along with our design we provide access to a web server implementation of our solution that facilitates pairwise comparison of short RNA (and DNA) sequences up to 200 base pairs. The web server first calculates the minimum free energy secondary structure and the base pair probabilities for up to 10 RNA or DNA sequences using RNAfold and then provides a two panel comparative view that includes CS2BP2-Plots along with the traditional graph, planar and circular diagrams obtained with VARNA. The CS2BP2-Plots include highlighting of the nucleotide differences between two selected sequences using ClustalW local alignments. We also provide descriptive statistics, dot-bracket secondary structure representations and ClustalW local alignments for compared sequences.

CONCLUSIONS

Using circular diagrams and colour and weight-coded arcs, we demonstrate how a single image can replace the state-of-the-art dual representations (dot-plots and minimum free energy structures) for base-pair probabilities of RNA secondary structures while allowing efficient exploration and comparison of different RNA conformations via a web server front end. With that, we provide the community, especially the biologically oriented, with an intuitive tool for ncRNA visualization. Web-server: https://nrcmonsrv01.nrc.ca/cs2bp2plot.

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/.