The Subviral RNA Database: a toolbox for viroids, the hepatitis delta virus and satellite RNAs research

Background  

Viroids, satellite RNAs, satellites viruses and the human hepatitis delta virus form the 'brotherhood' of the smallest known infectious RNA agents, known as the subviral RNAs. For most of these species, it is generally accepted that characteristics such as cell movement, replication, host specificity and pathogenicity are encoded in their RNA sequences and their resulting RNA structures. Although many sequences are indexed in publicly available databases, these sequence annotation databases do not provide the advanced searches and data manipulation capability for identifying and characterizing subviral RNA motifs.     

Small RNA regulation of ovule development in the cotton plant, <Emphasis Type="Italic">G. hirsutum</Emphasis>L

Background  

The involvement of small RNAs in cotton fiber development is under explored. The objective of this work was to directly clone, annotate, and analyze small RNAs of developing ovules to reveal the candidate small interfering RNA/microRNAs involved in cotton ovule and fiber development.     

Sampled ensemble neutrality as a feature to classify potential structured RNAs

Background

Structured RNAs have many biological functions ranging from catalysis of chemical reactions to gene regulation. Yet, many homologous structured RNAs display most of their conservation at the secondary or tertiary structure level. As a result, strategies for structured RNA discovery rely heavily on identification of sequences sharing a common stable secondary structure. However, correctly distinguishing structured RNAs from surrounding genomic sequence remains challenging, especially during de novo discovery. RNA also has a long history as a computational model for evolution due to the direct link between genotype (sequence) and phenotype (structure). From these studies it is clear that evolved RNA structures, like protein structures, can be considered robust to point mutations. In this context, an RNA sequence is considered robust if its neutrality (extent to which single mutant neighbors maintain the same secondary structure) is greater than that expected for an artificial sequence with the same minimum free energy structure.

Results

In this work, we bring concepts from evolutionary biology to bear on the structured RNA de novo discovery process. We hypothesize that alignments corresponding to structured RNAs should consist of neutral sequences. We evaluate several measures of neutrality for their ability to distinguish between alignments of structured RNA sequences drawn from Rfam and various decoy alignments. We also introduce a new measure of RNA structural neutrality, the structure ensemble neutrality (SEN). SEN seeks to increase the biological relevance of existing neutrality measures in two ways. First, it uses information from an alignment of homologous sequences to identify a conserved biologically relevant structure for comparison. Second, it only counts base-pairs of the original structure that are absent in the comparison structure and does not penalize the formation of additional base-pairs.

Conclusion

We find that several measures of neutrality are effective at separating structured RNAs from decoy sequences, including both shuffled alignments and flanking genomic sequence. Furthermore, as an independent feature classifier to identify structured RNAs, SEN yields comparable performance to current approaches that consider a variety of features including stability and sequence identity. Finally, SEN outperforms other measures of neutrality at detecting mutational robustness in bacterial regulatory RNA structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-014-1203-8) contains supplementary material, which is available to authorized users.     

A method for aligning RNA secondary structures and its application to RNA motif detection

Background  

Alignment of RNA secondary structures is important in studying functional RNA motifs. In recent years, much progress has been made in RNA motif finding and structure alignment. However, existing tools either require a large number of prealigned structures or suffer from high time complexities. This makes it difficult for the tools to process RNAs whose prealigned structures are unavailable or process very large RNA structure databases.     

Ribonucleoprotein particles containing non-coding Y RNAs, Ro60, La and nucleolin are not required for Y RNA function in DNA replication

Background

Ro ribonucleoprotein particles (Ro RNPs) consist of a non-coding Y RNA bound by Ro60, La and possibly other proteins. The physiological function of Ro RNPs is controversial as divergent functions have been reported for its different constituents. We have recently shown that Y RNAs are essential for the initiation of mammalian chromosomal DNA replication, whereas Ro RNPs are implicated in RNA stability and RNA quality control. Therefore, we investigate here the functional consequences of RNP formation between Ro60, La and nucleolin proteins with hY RNAs for human chromosomal DNA replication.

Methodology/Principal Findings

We first immunoprecipitated Ro60, La and nucleolin together with associated hY RNAs from HeLa cytosolic cell extract, and analysed the protein and RNA compositions of these precipitated RNPs by Western blotting and quantitative RT-PCR. We found that Y RNAs exist in several RNP complexes. One RNP comprises Ro60, La and hY RNA, and a different RNP comprises nucleolin and hY RNA. In addition about 50% of the Y RNAs in the extract are present outside of these two RNPs. Next, we immunodepleted these RNP complexes from the cytosolic extract and tested the ability of the depleted extracts to reconstitute DNA replication in a human cell-free system. We found that depletion of these RNP complexes from the cytosolic extract does not inhibit DNA replication in vitro. Finally, we tested if an excess of recombinant pure Ro or La protein inhibits Y RNA-dependent DNA replication in this cell-free system. We found that Ro60 and La proteins do not inhibit DNA replication in vitro.

Conclusions/Significance

We conclude that RNPs containing hY RNAs and Ro60, La or nucleolin are not required for the function of hY RNAs in chromosomal DNA replication in a human cell-free system, which can be mediated by Y RNAs outside of these RNPs. These data suggest that Y RNAs can support different cellular functions depending on associated proteins.     

A memory-efficient dynamic programming algorithm for optimal alignment of a sequence to an RNA secondary structure

Background  

Covariance models (CMs) are probabilistic models of RNA secondary structure, analogous to profile hidden Markov models of linear sequence. The dynamic programming algorithm for aligning a CM to an RNA sequence of length N is O(N ³) in memory. This is only practical for small RNAs.     

A population-based statistical approach identifies parameters characteristic of human microRNA-mRNA interactions

Background  

MicroRNAs are ~17–24 nt. noncoding RNAs found in all eukaryotes that degrade messenger RNAs via RNA interference (if they bind in a perfect or near-perfect complementarity to the target mRNA), or arrest translation (if the binding is imperfect). Several microRNA targets have been identified in lower organisms, but only one mammalian microRNA target has yet been validated experimentally.     

Ultra-high throughput sequencing-based small RNA discovery and discrete statistical biomarker analysis in a collection of cervical tumours and matched controls

Background  

Ultra-high throughput sequencing technologies provide opportunities both for discovery of novel molecular species and for detailed comparisons of gene expression patterns. Small RNA populations are particularly well suited to this analysis, as many different small RNAs can be completely sequenced in a single instrument run.     

Northern blotting analysis of microRNAs,their precursors and RNA interference triggers

Background  

Numerous microRNAs (miRNAs) have heterogeneous ends resulting from imprecise cleavages by processing nucleases and from various non-templated nucleotide additions. The scale of miRNA end-heterogeneity is best shown by deep sequencing data revealing not only the major miRNA variants but also those that occur in only minute amounts and are unlikely to be of functional importance. All RNA interference (RNAi) technology reagents that are expressed and processed in cells are also exposed to the same machinery generating end-heterogeneity of the released short interfering RNAs (siRNAs) or miRNA mimetics.     

RNAstructure: software for RNA secondary structure prediction and analysis

Background  

To understand an RNA sequence's mechanism of action, the structure must be known. Furthermore, target RNA structure is an important consideration in the design of small interfering RNAs and antisense DNA oligonucleotides. RNA secondary structure prediction, using thermodynamics, can be used to develop hypotheses about the structure of an RNA sequence.     

Controlled expression of functional miR-122 with a ligand inducible expression system

Background  

To study the biological function of miRNAs, and to achieve sustained or conditional gene silencing with siRNAs, systems that allow controlled expression of these small RNAs are desirable. Methods for cell delivery of siRNAs include transient transfection of synthetic siRNAs and expression of siRNAs in the form of short hairpins using constitutive RNA polymerase III promoters. Systems employing constitutive RNA polymerase II promoters have been used to express miRNAs. However, for many experimental systems these methods do not offer sufficient control over expression.     

A novel representation of RNA secondary structure based on element-contact graphs

Background  

Depending on their specific structures, noncoding RNAs (ncRNAs) play important roles in many biological processes. Interest in developing new topological indices based on RNA graphs has been revived in recent years, as such indices can be used to compare, identify and classify RNAs. Although the topological indices presented before characterize the main topological features of RNA secondary structures, information on RNA structural details is ignored to some degree. Therefore, it is necessity to identify topological features with low degeneracy based on complete and fine-grained RNA graphical representations.     

A structural analysis of <Emphasis Type="Italic">in vitro</Emphasis> catalytic activities of hammerhead ribozymes

Background  

Ribozymes are small catalytic RNAs that possess the dual functions of sequence-specific RNA recognition and site-specific cleavage. Trans-cleaving ribozymes can inhibit translation of genes at the messenger RNA (mRNA) level in both eukaryotic and prokaryotic systems and are thus useful tools for studies of gene function. However, identification of target sites for efficient cleavage poses a challenge. Here, we have considered a number of structural and thermodynamic parameters that can affect the efficiency of target cleavage, in an attempt to identify rules for the selection of functional ribozymes.     

Hierarchical folding of multiple sequence alignments for the prediction of structures and RNA-RNA interactions

Background  

Many regulatory non-coding RNAs (ncRNAs) function through complementary binding with mRNAs or other ncRNAs, e.g., microRNAs, snoRNAs and bacterial sRNAs. Predicting these RNA interactions is essential for functional studies of putative ncRNAs or for the design of artificial RNAs. Many ncRNAs show clear signs of undergoing compensating base changes over evolutionary time. Here, we postulate that a non-negligible part of the existing RNA-RNA interactions contain preserved but covarying patterns of interactions.     

An image processing approach to computing distances between RNA secondary structures dot plots

Background  

Computing the distance between two RNA secondary structures can contribute in understanding the functional relationship between them. When used repeatedly, such a procedure may lead to finding a query RNA structure of interest in a database of structures. Several methods are available for computing distances between RNAs represented as strings or graphs, but none utilize the RNA representation with dot plots. Since dot plots are essentially digital images, there is a clear motivation to devise an algorithm for computing the distance between dot plots based on image processing methods.