The RNA structure alignment ontology

DecoyFinder: Identification of Contaminants in Sets of Homologous RNA Sequences

Motivation

RNA structure is essential for the function of many non-coding RNAs. Using multiple homologous sequences, which share structure and function, secondary structure can be predicted with much higher accuracy than with a single sequence. It can be difficult, however, to establish a set of homologous sequences when their structure is not yet known. We developed a method to identify sequences in a set of putative homologs that are in fact non-homologs.

Results

Previously, we developed TurboFold to estimate conserved structure using multiple, unaligned RNA homologs. Here, we report that the positive predictive value of TurboFold is significantly reduced by the presence of contamination by non-homologous sequences, although the reduction is less than 1%. We developed a method called DecoyFinder, which applies machine learning trained with features determined by TurboFold, to detect sequences that are not homologous with the other sequences in the set. This method can identify approximately 45% of non-homologous sequences, at a rate of 5% misidentification of true homologous sequences.

Availability

DecoyFinder and TurboFold are incorporated in RNAstructure, which is provided for free and open source under the GPL V2 license. It can be downloaded at http://rna.urmc.rochester.edu/RNAstructure.html.

Plant Ribonomics: Proteins in Search of RNA Partners

Research into the regulation of gene expression underwent a shift from focusing on DNA-binding proteins as key transcriptional regulators to RNA-binding proteins (RBPs) that come into play once transcription has been initiated. RBPs orchestrate all RNA-processing steps in the cell. To obtain a global view of in vivo targets, the RNA complement associated with particular RBPs is determined via immunoprecipitation of the RBP and subsequent identification of bound RNAs via RNA-seq. Here, we describe technical advances in identifying RBP in vivo targets and their binding motifs. We provide an up-to-date view of targets of nucleocytoplasmic RBPs collected in arabidopsis. We also discuss current experimental limitations and provide an outlook on how the approaches may advance our understanding of post-transcriptional networks.

Aptamer Selection Technology and Recent Advances

Over the last decade, aptamers have begun to find their way from basic research to diverse commercial applications. The development of diagnostics is even more widespread than clinical applications because aptamers do not have to be extensively modified to enhance their in vivo stability and pharmacokinetics in diagnostic assays. The increasing attention has propelled the technical progress of the in vitro selection technology (SELEX) to enhance the efficiency of developing aptamers for commercially interesting targets. This review highlights recent progress in the technical steps of a SELEX experiment with a focus on high-throughput next-generation sequencing and bioinformatics. Achievements have been made in the optimization of aptamer libraries, separation schemes, amplification of the selected libraries and the identification of aptamer sequences from enriched libraries.

JACALIN-LECTIN LIKE1 Regulates the Nuclear Accumulation of GLYCINE-RICH RNA-BINDING PROTEIN7, Influencing the RNA Processing of FLOWERING LOCUS C Antisense Transcripts and Flowering Time in Arabidopsis

Lectins selectively recognize sugars or glycans for defense in living cells, but less is known about their roles in the development process and the functional network with other factors. Here, we show that Arabidopsis (Arabidopsis thaliana) JACALIN-LECTIN LIKE1 (AtJAC1) functions in flowering time control. Loss of function of AtJAC1 leads to precocious flowering, whereas overexpression of AtJAC1 causes delayed flowering. AtJAC1 influences flowering through regulation of the key flowering repressor gene FLOWERING LOCUS C (FLC). Genetic analysis revealed that AtJAC1's function is mostly dependent on GLYCINE-RICH RNA-BINDING PROTEIN7 (GRP7), an upstream regulator of FLC. Biochemical and cell biological data indicated that AtJAC1 interacted physically with GRP7 specifically in the cytoplasm. AtJAC1 influences the nucleocytoplasmic distribution of GRP7, with predominant nuclear localization of GRP7 when AtJAC1 function is lost but retention of GRP7 in the cytoplasm when AtJAC1 is overexpressed. A temporal inducible assay suggested that AtJAC1's regulation of flowering could be compromised by the nuclear accumulation of GRP7. In addition, GRP7 binds to the antisense precursor messenger RNA of FLC through a conserved RNA motif. Loss of GRP7 function leads to the elevation of total FLC antisense transcripts and reduced proximal-distal polyadenylation ratio, as well as histone methylation changes in the FLC gene body region and increased total functional sense FLC transcript. Attenuating the direct binding of GRP7 with competing artificial RNAs leads to changes of FLC antisense precursor messenger RNA processing and flowering transition. Taken together, our study indicates that AtJAC1 coordinates with GRP7 in shaping plant development through the regulation of RNA processing in Arabidopsis.

Ambivalent covariance models

BACKGROUND

Evolutionary variations let us define a set of similar nucleic acid sequences as a family if these different molecules execute a common function. Capturing their sequence variation by using e. g. position specific scoring matrices significantly improves sensitivity of detection tools. Members of a functional (non-coding) RNA family are affected by these variations not only on the sequence, but also on the structural level. For example, some transfer-RNAs exhibit a fifth helix in addition to the typical cloverleaf structure. Current covariance models - the unrivaled homology search approach for structured RNA - do not benefit from structural variation within a family, but rather penalize it. This leads to artificial subdivision of families and loss of information in the RFAM database.

RESULTS

We propose an extension to the fundamental architecture of covariance models to allow for several, compatible consensus structures. The resulting models are called ambivalent covariance models. Evaluation on several RFAM families shows that coalescence of structural variation within a family by using ambivalent consensus models is superior to subdividing the family into multiple classical covariance models.

CONCLUSION

A prototype and source code is available at http://bibiserv.cebitec.uni-bielefeld.de/acms.

RNA-binding protein immunoprecipitation from whole-cell extracts

RNA-based regulation is increasingly recognized as an important factor shaping the cellular transcriptome. RNA-binding proteins that interact with cis-regulatory motifs within pre-mRNAs determine the fate of their targets. Understanding posttranscriptional networks controlled by an RNA-binding protein requires the identification of its immediate in vivo targets. Here we describe RNA immunoprecipitation in Arabidopsis thaliana. Transgenic plants expressing an RNA-binding protein fused to green fluorescent protein are treated with formaldehyde to "trap" RNAs in complexes with their physiological protein partners. A whole-cell extract is subjected to immunoprecipitation with an antibody against the GFP tag. In parallel, a mock immunoprecipitation is performed using an unrelated antibody. Coprecipitated RNAs are eluted from the immunoprecipitate and identified via real-time PCR. Enrichment relative to immunoprecipitation from plants expressing GFP only and mock immunoprecipitation with an unrelated antibody indicates specific binding.

The art of editing RNA structural alignments

Manual editing of RNA structural alignments may be considered more art than science, since it still requires an expert biologist to take multiple levels of information into account and be slightly creative when constructing high-quality alignments. Even though the task is rather tedious, it is rewarded by great insight into the evolution of structure and function of your favorite RNA molecule. In this chapter I will review the methods and considerations that go into constructing RNA structural alignments at the secondary and tertiary structure level; introduce software, databases, and algorithms that have proven useful in semiautomating the work process; and suggest future directions towards full automatization.

R3D Align web server for global nucleotide to nucleotide alignments of RNA 3D structures

The R3D Align web server provides online access to ‘RNA 3D Align’ (R3D Align), a method for producing accurate nucleotide-level structural alignments of RNA 3D structures. The web server provides a streamlined and intuitive interface, input data validation and output that is more extensive and easier to read and interpret than related servers. The R3D Align web server offers a unique Gallery of Featured Alignments, providing immediate access to pre-computed alignments of large RNA 3D structures, including all ribosomal RNAs, as well as guidance on effective use of the server and interpretation of the output. By accessing the non-redundant lists of RNA 3D structures provided by the Bowling Green State University RNA group, R3D Align connects users to structure files in the same equivalence class and the best-modeled representative structure from each group. The R3D Align web server is freely accessible at http://rna.bgsu.edu/r3dalign/.

Nucleotides adjacent to the ligand-binding pocket are linked to activity tuning in the purine riboswitch

Direct sensing of intracellular metabolite concentrations by riboswitch RNAs provides an economical and rapid means to maintain metabolic homeostasis. Since many organisms employ the same class of riboswitch to control different genes or transcription units, it is likely that functional variation exists in riboswitches such that activity is tuned to meet cellular needs. Using a bioinformatic approach, we have identified a region of the purine riboswitch aptamer domain that displays conservation patterns linked to riboswitch activity. Aptamer domain compositions within this region can be divided into nine classes that display a spectrum of activities. Naturally occurring compositions in this region favor rapid association rate constants and slow dissociation rate constants for ligand binding. Using X-ray crystallography and chemical probing, we demonstrate that both the free and bound states are influenced by the composition of this region and that modest sequence alterations have a dramatic impact on activity. The introduction of non-natural compositions result in the inability to regulate gene expression in vivo, suggesting that aptamer domain activity is highly plastic and thus readily tunable to meet cellular needs.

An hnRNP-like RNA-binding protein affects alternative splicing by in vivo interaction with transcripts in Arabidopsis thaliana

Alternative splicing (AS) of pre-mRNAs is an important regulatory mechanism shaping the transcriptome. In plants, only few RNA-binding proteins are known to affect AS. Here, we show that the glycine-rich RNA-binding protein AtGRP7 influences AS in Arabidopsis thaliana. Using a high-resolution RT–PCR-based AS panel, we found significant changes in the ratios of AS isoforms for 59 of 288 analyzed AS events upon ectopic AtGRP7 expression. In particular, AtGRP7 affected the choice of alternative 5′ splice sites preferentially. About half of the events are also influenced by the paralog AtGRP8, indicating that AtGRP7 and AtGRP8 share a network of downstream targets. For 10 events, the AS patterns were altered in opposite directions in plants with elevated AtGRP7 level or lacking AtGRP7. Importantly, RNA immunoprecipitation from plant extracts showed that several transcripts are bound by AtGRP7 in vivo and indeed represent direct targets. Furthermore, the effect of AtGRP7 on these AS events was abrogated by mutation of a single arginine that is required for its RNA-binding activity. This indicates that AtGRP7 impacts AS of these transcripts via direct interaction. As several of the AS events are also controlled by other splicing regulators, our data begin to provide insights into an AS network in Arabidopsis.

RNASTAR: an RNA STructural Alignment Repository that provides insight into the evolution of natural and artificial RNAs

Automated RNA alignment algorithms often fail to recapture the essential conserved sites that are critical for function. To assist in the refinement of these algorithms, we manually curated a set of 148 alignments with a total of 9600 unique sequences, in which each alignment was backed by at least one crystal or NMR structure. These alignments included both naturally and artificially selected molecules. We used principles of isostericity to improve the alignments from an average of 83%-94% isosteric base pairs. We expect that this alignment collection will assist in a wide range of benchmarking efforts and provide new insight into evolutionary principles governing change in RNA structural motifs. The improved alignments have been contributed to the Rfam database.

rCAD: A Novel Database Schema for the Comparative Analysis of RNA

Beyond its direct involvement in protein synthesis with mRNA, tRNA, and rRNA, RNA is now being appreciated for its significance in the overall metabolism and regulation of the cell. Comparative analysis has been very effective in the identification and characterization of RNA molecules, including the accurate prediction of their secondary structure. We are developing an integrative scalable data management and analysis system, the RNA Comparative Analysis Database (rCAD), implemented with SQL Server to support RNA comparative analysis. The platformagnostic database schema of rCAD captures the essential relationships between the different dimensions of information for RNA comparative analysis datasets. The rCAD implementation enables a variety of comparative analysis manipulations with multiple integrated data dimensions for advanced RNA comparative analysis workflows. In this paper, we describe details of the rCAD schema design and illustrate its usefulness with two usage scenarios.

RNA-mediated silencing mechanisms in mammalian cells

Noncoding RNAs are a structural component of the nuclear scaffold and have been implicated in controlling gene expression. In mammals, long noncoding RNAs contribute to the regulation of imprinted gene expression, dosage compensation, development, and tumorigenesis. RNA is also a component of pericentric heterochromatin, and transcripts have been identified at the chromosomal telomeres. The functions of noncoding RNAs are likely diverse, and their underlying mechanisms are just beginning to be understood. Several noncoding RNAs interact with chromatin-modifying complexes and might have a role in targeting chromatin modifications to specific regions of the genome. This suggests a prominent function of RNA in establishing histone modification and DNA methylation patterns in development. Studies on model systems such as X inactivation, the regulation of the Hox clusters, and genomic imprinting have begun to shed light on the role of noncoding RNAs in chromosomal organization and regulation of gene expression. Well-studied examples of noncoding RNAs include Xist, Air, Kcnq1ot1, HOTAIR, and Tsix. Here, a concise review of noncoding RNA function in mammals is given, and the present understanding and future directions of the field are summarized.

Meeting report of the RNA Ontology Consortium January 8-9, 2011

This report summarizes the proceedings of the structure mapping working group meeting of the RNA Ontology Consortium (ROC), held in Kona, Hawaii on January 8-9, 2011. The ROC hosted this workshop to facilitate collaborations among those researchers formalizing concepts in RNA, those developing RNA-related software, and those performing genome annotation and standardization. The workshop included three software presentations, extended round-table discussions, and the constitution of two new working groups, the first to address the need for better software integration and the second to discuss standardization and benchmarking of existing RNA annotation pipelines. These working groups have subsequently pursued concrete implementation of actions suggested during the discussion. Further information about the ROC and its activities can be found at http://roc.bgsu.edu/.

Sharing and archiving nucleic acid structure mapping data

Nucleic acids are particularly amenable to structural characterization using chemical and enzymatic probes. Each individual structure mapping experiment reveals specific information about the structure and/or dynamics of the nucleic acid. Currently, there is no simple approach for making these data publically available in a standardized format. We therefore developed a standard for reporting the results of single nucleotide resolution nucleic acid structure mapping experiments, or SNRNASMs. We propose a schema for sharing nucleic acid chemical probing data that uses generic public servers for storing, retrieving, and searching the data. We have also developed a consistent nomenclature (ontology) within the Ontology of Biomedical Investigations (OBI), which provides unique identifiers (termed persistent URLs, or PURLs) for classifying the data. Links to standardized data sets shared using our proposed format along with a tutorial and links to templates can be found at http://snrnasm.bio.unc.edu.

Boulder ALignment Editor (ALE): a web-based RNA alignment tool

Summary: The explosion of interest in non-coding RNAs, together with improvements in RNA X-ray crystallography, has led to a rapid increase in RNA structures at atomic resolution from 847 in 2005 to 1900 in 2010. The success of whole-genome sequencing has led to an explosive growth of unaligned homologous sequences. Consequently, there is a compelling and urgent need for user-friendly tools for producing structure-informed RNA alignments. Most alignment software considers the primary sequence alone; some specialized alignment software can also include Watson–Crick base pairs, but none adequately addresses the needs introduced by the rapid influx of both sequence and structural data. Therefore, we have developed the Boulder ALignment Editor (ALE), which is a web-based RNA alignment editor, designed for editing and assessing alignments using structural information. Some features of BoulderALE include the annotation and evaluation of an alignment based on isostericity of Watson–Crick and non-Watson–Crick base pairs, along with the collapsing (horizontally and vertically) of the alignment, while maintaining the ability to edit the alignment.

Availability:http://www.microbio.me/boulderale.

Contact:jesse.stombaugh@colorado.edu

TurboFold: iterative probabilistic estimation of secondary structures for multiple RNA sequences

Background

The prediction of secondary structure, i.e. the set of canonical base pairs between nucleotides, is a first step in developing an understanding of the function of an RNA sequence. The most accurate computational methods predict conserved structures for a set of homologous RNA sequences. These methods usually suffer from high computational complexity. In this paper, TurboFold, a novel and efficient method for secondary structure prediction for multiple RNA sequences, is presented.

Results

TurboFold takes, as input, a set of homologous RNA sequences and outputs estimates of the base pairing probabilities for each sequence. The base pairing probabilities for a sequence are estimated by combining intrinsic information, derived from the sequence itself via the nearest neighbor thermodynamic model, with extrinsic information, derived from the other sequences in the input set. For a given sequence, the extrinsic information is computed by using pairwise-sequence-alignment-based probabilities for co-incidence with each of the other sequences, along with estimated base pairing probabilities, from the previous iteration, for the other sequences. The extrinsic information is introduced as free energy modifications for base pairing in a partition function computation based on the nearest neighbor thermodynamic model. This process yields updated estimates of base pairing probability. The updated base pairing probabilities in turn are used to recompute extrinsic information, resulting in the overall iterative estimation procedure that defines TurboFold.

TurboFold is benchmarked on a number of ncRNA datasets and compared against alternative secondary structure prediction methods. The iterative procedure in TurboFold is shown to improve estimates of base pairing probability with each iteration, though only small gains are obtained beyond three iterations. Secondary structures composed of base pairs with estimated probabilities higher than a significance threshold are shown to be more accurate for TurboFold than for alternative methods that estimate base pairing probabilities. TurboFold-MEA, which uses base pairing probabilities from TurboFold in a maximum expected accuracy algorithm for secondary structure prediction, has accuracy comparable to the best performing secondary structure prediction methods. The computational and memory requirements for TurboFold are modest and, in terms of sequence length and number of sequences, scale much more favorably than joint alignment and folding algorithms.

Conclusions

TurboFold is an iterative probabilistic method for predicting secondary structures for multiple RNA sequences that efficiently and accurately combines the information from the comparative analysis between sequences with the thermodynamic folding model. Unlike most other multi-sequence structure prediction methods, TurboFold does not enforce strict commonality of structures and is therefore useful for predicting structures for homologous sequences that have diverged significantly. TurboFold can be downloaded as part of the RNAstructure package at http://rna.urmc.rochester.edu.

Nuclear and cytoplasmic localization of neural stem cell microRNAs

Although generally regarded as functional in the cytoplasm, a number of microRNAs (miRNAs) have been found in the nucleus, possibly with a role in gene regulation. Here we report that, in fact, a substantial fraction of all human miRNAs are present in the nucleus of neural stem cells. Further, subsets of these miRNAs display consistently higher standardized rank in the nucleus than in the cytoplasm of these cells, as identified with an RT-qPCR technology and confirmed by microarray analysis. Likewise, other miRNAs display higher cytoplasmic standardized ranks. Three samples were partitioned into nuclear and cytoplasmic fractions in six assays for 373 miRNAs. From the 100 most highly expressed miRNAs, standard scores of nuclear and cytoplasmic concentrations were determined. Among those, 21 miRNAs had all three nuclear standard scores higher than all three cytoplasmic scores; likewise, 31 miRNAs had consistently higher cytoplasmic scores. Random concentrations would result in only five in each set. Remarkably, if one miRNA has a high standard score in a compartment, then other miRNAs having the same 5' seeds and certain similar 3' end patterns are also highly scored in the same way. That is, in addition to the seed sequence, 3' sequence similarity criteria identify families of mature miRNAs with consistently high nuclear or cytoplasmic expression.

Predicting and modeling RNA architecture

A general approach for modeling the architecture of large and structured RNA molecules is described. The method exploits the modularity and the hierarchical folding of RNA architecture that is viewed as the assembly of preformed double-stranded helices defined by Watson-Crick base pairs and RNA modules maintained by non-Watson-Crick base pairs. Despite the extensive molecular neutrality observed in RNA structures, specificity in RNA folding is achieved through global constraints like lengths of helices, coaxiality of helical stacks, and structures adopted at the junctions of helices. The Assemble integrated suite of computer tools allows for sequence and structure analysis as well as interactive modeling by homology or ab initio assembly with possibilities for fitting within electronic density maps. The local key role of non-Watson-Crick pairs guides RNA architecture formation and offers metrics for assessing the accuracy of three-dimensional models in a more useful way than usual root mean square deviation (RMSD) values.

Spotlight on post-transcriptional control in the circadian system

An endogenous timing mechanism, the circadian clock, causes rhythmic expression of a considerable fraction of the genome of most organisms to optimally align physiology and behavior with their environment. Circadian clocks are self-sustained oscillators primarily based on transcriptional feedback loops and post-translational modification of clock proteins. It is increasingly becoming clear that regulation at the RNA level strongly impacts the cellular circadian transcriptome and proteome as well as the oscillator mechanism itself. This review focuses on posttranscriptional events, discussing RNA-binding proteins that, by influencing the timing of pre-mRNA splicing, polyadenylation and RNA decay, shape rhythmic expression profiles. Furthermore, recent findings on the contribution of microRNAs to orchestrating circadian rhythms are summarized.

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.

Rfam: Wikipedia, clans and the "decimal" release

The Rfam database aims to catalogue non-coding RNAs through the use of sequence alignments and statistical profile models known as covariance models. In this contribution, we discuss the pros and cons of using the online encyclopedia, Wikipedia, as a source of community-derived annotation. We discuss the addition of groupings of related RNA families into clans and new developments to the website. Rfam is available on the Web at http://rfam.sanger.ac.uk.

Global transcript profiling of transgenic plants constitutively overexpressing the RNA-binding protein AtGRP7

BACKGROUND

The clock-controlled RNA-binding protein AtGRP7 influences circadian oscillations of its own transcript at the post-transcriptional level. To identify additional targets that are regulated by AtGRP7, transcript profiles of transgenic plants constitutively overexpressing AtGRP7 (AtGRP7-ox) and wild type plants were compared.

RESULTS

Approximately 1.4% of the transcripts represented on the Affymetrix ATH1 microarray showed changes in steady-state abundance upon AtGRP7 overexpression. One third of the differentially expressed genes are controlled by the circadian clock, and they show a distinct bias of their phase: The up-regulated genes preferentially peak around dawn, roughly opposite to the AtGRP7 peak abundance whereas the down-regulated genes preferentially peak at the end of the day. Further, transcripts responsive to abiotic and biotic stimuli were enriched among AtGRP7 targets. Transcripts encoding the pathogenesis-related PR1 and PR2 proteins were elevated in AtGRP7-ox plants but not in plants overexpressing AtGRP7 with a point mutation in the RNA-binding domain, indicating that the regulation involves RNA binding activity of AtGRP7. Gene set enrichment analysis uncovered components involved in ribosome function and RNA metabolism among groups of genes upregulated in AtGRP7-ox plants, consistent with its role in post-transcriptional regulation.

CONCLUSION

Apart from regulating a suite of circadian transcripts in a time-of-day dependent manner AtGRP7, both directly and indirectly, affects other transcripts including transcripts responsive to abiotic and biotic stimuli. This suggests a regulatory role of AtGRP7 in the output of the endogenous clock and a complex network of transcripts responsive to external stimuli downstream of the AtGRP7 autoregulatory circuit.

R3D Align: global pairwise alignment of RNA 3D structures using local superpositions

MOTIVATION

Comparing 3D structures of homologous RNA molecules yields information about sequence and structural variability. To compare large RNA 3D structures, accurate automatic comparison tools are needed. In this article, we introduce a new algorithm and web server to align large homologous RNA structures nucleotide by nucleotide using local superpositions that accommodate the flexibility of RNA molecules. Local alignments are merged to form a global alignment by employing a maximum clique algorithm on a specially defined graph that we call the 'local alignment' graph.

RESULTS

The algorithm is implemented in a program suite and web server called 'R3D Align'. The R3D Align alignment of homologous 3D structures of 5S, 16S and 23S rRNA was compared to a high-quality hand alignment. A full comparison of the 16S alignment with the other state-of-the-art methods is also provided. The R3D Align program suite includes new diagnostic tools for the structural evaluation of RNA alignments. The R3D Align alignments were compared to those produced by other programs and were found to be the most accurate, in comparison with a high quality hand-crafted alignment and in conjunction with a series of other diagnostics presented. The number of aligned base pairs as well as measures of geometric similarity are used to evaluate the accuracy of the alignments.

AVAILABILITY

R3D Align is freely available through a web server http://rna.bgsu.edu/R3DAlign. The MATLAB source code of the program suite is also freely available for download at that location.

Computational approaches to 3D modeling of RNA

Many exciting discoveries have recently revealed the versatility of RNA and its importance in a variety of functions within the cell. Since the structural features of RNA are of major importance to their biological function, there is much interest in predicting RNA structure, either in free form or in interaction with various ligands, including proteins, metabolites and other molecules. In recent years, an increasing number of researchers have developed novel RNA algorithms for predicting RNA secondary and tertiary structures. In this review, we describe current experimental and computational advances and discuss recent ideas that are transforming the traditional view of RNA folding. To evaluate the performance of the most recent RNA 3D folding algorithms, we provide a comparative study in order to test the performance of available 3D structure prediction algorithms for an RNA data set of 43 structures of various lengths and motifs. We find that the algorithms vary widely in terms of prediction quality across different RNA lengths and topologies; most predictions have very large root mean square deviations from the experimental structure. We conclude by outlining some suggestions for future RNA folding research.

Visualization of multiple alignments, phylogenies and gene family evolution

Software for visualizing sequence alignments and trees are essential tools for life scientists. In this review, we describe the major features and capabilities of a selection of stand-alone and web-based applications useful when investigating the function and evolution of a gene family. These range from simple viewers, to systems that provide sophisticated editing and analysis functions. We conclude with a discussion of the challenges that these tools now face due to the flood of next generation sequence data and the increasingly complex network of bioinformatics information sources.