RNA Secondary Structure Models Research Articles

Evaluation of the authenticity of a highly novel environmental sequence from boreal forest soil using ribosomal RNA secondary structure modeling

The number of sequences from both formally described taxa and uncultured environmental DNA deposited in the International Nucleotide Sequence Databases has increased substantially over the last two decades. Although the majority of these sequences represent authentic gene copies, there is evidence of DNA artifacts in these databases as well. These include lab artifacts, such as PCR chimeras, and biological artifacts such as pseudogenes or other paralogous sequences. Sequences that fall in basal positions in phylogenetic trees and appear distant from known sequences are particularly suspect. Phylogenetic analyses suggest that a novel sequence type (NS1) found in two boreal forest soil clone libraries belongs to the fungal kingdom but does not fall unambiguously within any known phylum. We have evaluated this sequence type using an array of secondary-structure analyses. To our knowledge, such analyses have never been used on environmental ribosomal sequences. Ribosomal secondary structure was modeled for four rRNA loci (ITS1, 5.8S, ITS2, 5′ LSU). These models were analyzed for the presence of conserved domains, conserved nucleotide motifs, and compensatory base changes. Minimal free energy (MFE) foldings and GC contents of sequences representing the major fungal clades, as well as NS1, were also compared. NS1 displays secondary rRNA structures consistent with other fungi and many, but not all, conserved nucleotide motifs found across eukaryotes. However, our analyses show that many other authentic sequences from basal fungi lack more of these conserved motifs than does NS1. Together our findings suggest that NS1 represents an authentic gene copy. The methods described here can be used on any rRNA-coding sequence, not just environmental fungal sequences. As new-generation sequencing methods that yield shorter sequences become more widely implemented, methods that evaluate sequence authenticity should also be more widely implemented. For fungi, the adjacent 5.8S and ITS2 loci should be prioritized. This region is not only suited to distinguishing between closely related species, but it is also more informative in terms of expected secondary structure.

Statistical Analysis of SHAPE-Directed RNA Secondary Structure Modeling

The ability to predict RNA secondary structure is fundamental for understanding and manipulating RNA function. The information obtained from selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) experiments greatly improves the accuracy of RNA secondary structure prediction. Recently, Das and colleagues [Kladwang, W., et al. (2011) Biochemistry 50, 8049-8056] proposed a "bootstrapping" approach for estimating the variance and helix-by-helix confidence levels of predicted secondary structures based on resampling (randomizing and summing) the measured SHAPE data. We show that the specific resampling approach described by Kladwang et al. introduces systematic errors and underestimates confidence in secondary structure prediction using SHAPE data. Instead, a leave-data-out jackknife approach better estimates the influence of a given experimental data set on SHAPE-directed secondary structure modeling. Even when 35% of the data were left out in the jackknife approach, the confidence levels of SHAPE-directed secondary structure prediction were significantly higher than those calculated by Das and colleagues using bootstrapping. Helix confidence levels were thus underestimated in the recent study, and the resampling approach implemented by Kladwang et al. is not an appropriate metric for evaluating SHAPE-directed secondary structure modeling.

Statistical Analysis of SHAPE-Directed RNA Secondary Structure Modeling

Statistical Analysis of SHAPE-Directed RNA Secondary Structure Modeling

Accurate SHAPE-directed RNA secondary structure modeling, including pseudoknots

Accurate SHAPE-directed RNA secondary structure modeling, including pseudoknots

Quantitative Dimethyl Sulfate Mapping for Automated RNA Secondary Structure Inference

For decades, dimethyl sulfate (DMS) mapping has informed manual modeling of RNA structure in vitro and in vivo. Here, we incorporate DMS data into automated secondary structure inference using an energy minimization framework developed for 2'-OH acylation (SHAPE) mapping. On six noncoding RNAs with crystallographic models, DMS-guided modeling achieves overall false negative and false discovery rates of 9.5% and 11.6%, respectively, comparable to or better than those of SHAPE-guided modeling, and bootstrapping provides straightforward confidence estimates. Integrating DMS-SHAPE data and including 1-cyclohexyl(2-morpholinoethyl) carbodiimide metho-p-toluene sulfonate (CMCT) reactivities provide small additional improvements. These results establish DMS mapping, an already routine technique, as a quantitative tool for unbiased RNA secondary structure modeling.

Database for bacterial group II introns

The Database for Bacterial Group II Introns (http://webapps2.ucalgary.ca/~groupii/index.html#) provides a catalogue of full-length, non-redundant group II introns present in bacterial DNA sequences in GenBank. The website is divided into three sections. The first section provides general information on group II intron properties, structures and classification. The second and main section lists information for individual introns, including insertion sites, DNA sequences, intron-encoded protein sequences and RNA secondary structure models. The final section provides tools for identification and analysis of intron sequences. These include a step-by-step guide to identify introns in genomic sequences, a local BLAST tool to identify closest intron relatives to a query sequence, and a boundary-finding tool that predicts 5′ and 3′ intron–exon junctions in an input DNA sequence. Finally, selected intron data can be downloaded in FASTA format. It is hoped that this database will be a useful resource not only to group II intron and RNA researchers, but also to microbiologists who encounter these unexpected introns in genomic sequences.

Analysis of the Free Energy in a Stochastic RNA Secondary Structure Model

There are two custom ways for predicting RNA secondary structures: minimizing the free energy of a conformation according to a thermodynamic model and maximizing the probability of a folding according to a stochastic model. In most cases, stochastic grammars are used for the latter alternative applying the maximum likelihood principle for determining a grammar's probabilities. In this paper, building on such a stochastic model, we will analyze the expected minimum free energy of an RNA molecule according to Turner's energy rules. Even if the parameters of our grammar are chosen with respect to structural properties of native molecules only (and therefore, independent of molecules' free energy), we prove formulae for the expected minimum free energy and the corresponding variance as functions of the molecule's size which perfectly fit the native behavior of free energies. This gives proof for a high quality of our stochastic model making it a handy tool for further investigations. In fact, the stochastic model for RNA secondary structures presented in this work has, for example, been used as the basis of a new algorithm for the (nonuniform) generation of random RNA secondary structures.

The highly variable mitochondrial small-subunit ribosomal RNA gene of Ophiostoma minus

Mitochondrial genomes in the true fungi are highly variable both in size and organization. Most of this size variation is due to the presence of introns and intron-encoded open reading frames (ORFs). The objectives for this work were to examine the mitochondrial small-subunit ribosomal RNA ( rns) gene of strains of Ophiostoma minus for the presence of introns and to characterize such introns and their encoded ORFs. DNA sequence analysis showed that among different strains of O. minus various rns gene exon/intron configurations can be observed. Based on comparative sequence analysis and RNA secondary structure modeling group I introns with LAGLIDADG ORFs were uncovered at positions mS569 and mS1224 and group II introns were present at positions mS379 and mS952. The mS379 group II intron encoded a fragmented reverse transcriptase (RT)-like ORF and the mS952 group II intron encoded a LAGLIDADG-type ORF. Examples of intron ORF degeneration due to frameshift mutations were observed. The mS379 group II intron is the first mitochondrial group II intron to have an ORF inserted within domain II, typically RT-like ORFs are inserted in domain IV. The evolutionary dynamics of the intron-encoded ORFs have also been examined.

Genome-Wide Double-Stranded RNA Sequencing Reveals the Functional Significance of Base-Paired RNAs in Arabidopsis

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA–producing substrates of RNA–DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA–DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA.

Correlation between recombination junctions and RNA secondary structure elements in poliovirus Sabin strains.

In order to test the hypothesis that RNA structural elements promote the distribution of certain types of recombination junctions in each one of the 2C and 3D poliovirus genomic regions (Sabin 3/Sabin 2 or Sabin 1 in 2C and Sabin 2/Sabin 1 or Sabin 3 in 3D), we searched in 2C and 3D regions of reference Sabin strains for high probability RNA structural elements that could promote recombination. Recombination junctions that were identified in clinical strains of this study, as well as in clinical strains of previous studies, were superimposed on RNA secondary structure models of 2C and 3D genomic regions. Furthermore, we created an in vitro model, based on double infection of cell-culture with two poliovirus strains, for the production and identification of recombinant Sabin strains in 2C and 3D regions. Our intention was to compare the results that refer to the correlation of recombination junctions and RNA secondary structures in 2C and 3D regions of clinical strains, with the respective results of the in vitro model. Most of the recombination junctions of the clinical strains were correlated with RNA secondary structure elements, which were identical between recombining Sabin strains, and also presented high predictive value. In consensus were, the respective results originated from the in vitro model. We propose that the distribution of specific types of recombination junctions in certain regions of Sabin strains is not fortuitous and is correlated with RNA secondary structure elements identical to both recombination partners. Furthermore, results of this study highlight an important role for the stem region of the RNA structure elements in promoting recombination.Electronic supplementary materialThe online version of this article (doi:10.1007/s11262-010-0512-5) contains supplementary material, which is available to authorized users.

Placement of attine ant-associated Pseudonocardia in a global Pseudonocardia phylogeny (Pseudonocardiaceae, Actinomycetales): a test of two symbiont-association models.

We reconstruct the phylogenetic relationships within the bacterial genus Pseudonocardia to evaluate two models explaining how and why Pseudonocardia bacteria colonize the microbial communities on the integument of fungus-gardening ant species (Attini, Formicidae). The traditional Coevolution-Codivergence model views the integument-colonizing Pseudonocardia as mutualistic microbes that are largely vertically transmitted between ant generations and that supply antibiotics that specifically suppress the garden pathogen Escovopsis. The more recent Acquisition model views Pseudonocardia as part of a larger integumental microbe community that frequently colonizes the ant integument from environmental sources (e.g., soil, plant material). Under this latter model, ant-associated Pseudonocardia may have diverse ecological roles on the ant integument (possibly ranging from pathogenic, to commensal, to mutualistic) and are not necessarily related to Escovopsis suppression. We test distinct predictions of these two models regarding the phylogenetic proximity of ant-associated and environmental Pseudonocardia. We amassed 16S-rRNA gene sequence information for 87 attine-associated and 238 environmental Pseudonocardia, aligned the sequences with the help of RNA secondary structure modeling, and reconstructed phylogenetic relationships using a maximum-likelihood approach. We present 16S-rRNA secondary structure models of representative Pseudonocardia species to improve sequence alignments and identify sequencing errors. Our phylogenetic analyses reveal close affinities and even identical sequence matches between environmental Pseudonocardia and ant-associated Pseudonocardia, as well as nesting of environmental Pseudonocardia in subgroups that were previously thought to be specialized to associate only with attine ants. The great majority of ant-associated Pseudonocardia are closely related to autotrophic Pseudonocardia and are placed in a large subgroup of Pseudonocardia that is known essentially only from cultured isolates (rather than cloned 16S sequences). The preponderance of the known ant-associated Pseudonocardia in this latter clade of culturable lineages may not necessarily reflect abundance of these Pseudonocardia types on the ants, but isolation biases when screening for Pseudonocardia (e.g., preferential isolation of autotrophic Pseudonocardia with minimum-nutrient media). The accumulated phylogenetic patterns and the possibility of isolation biases in previous work further erode support for the traditional Coevolution-Codivergence model and calls for continued revision of our understanding how and why Pseudonocardia colonize the microbial communities on the integument of fungus-gardening ant species.

A growth model for RNA secondary structures

A hierarchical model for the growth of planar arch structures for RNA secondary structuresis presented, and shown to be equivalent to a tree-growth model. Both models can besolved analytically, giving access to scaling functions for large molecules, andcorrections to scaling, checked by numerical simulations of up to 6500 bases. Theequivalence of both models should be helpful in understanding more general tree-growthprocesses.

The MC-Fold and MC-Sym pipeline infers RNA structure from sequence data

The classical RNA secondary structure model considers A.U and G.C Watson-Crick as well as G.U wobble base pairs. Here we substitute it for a new one, in which sets of nucleotide cyclic motifs define RNA structures. This model allows us to unify all base pairing energetic contributions in an effective scoring function to tackle the problem of RNA folding. We show how pipelining two computer algorithms based on nucleotide cyclic motifs, MC-Fold and MC-Sym, reproduces a series of experimentally determined RNA three-dimensional structures from the sequence. This demonstrates how crucial the consideration of all base-pairing interactions is in filling the gap between sequence and structure. We use the pipeline to define rules of precursor microRNA folding in double helices, despite the presence of a number of presumed mismatches and bulges, and to propose a new model of the human immunodeficiency virus-1 -1 frame-shifting element.

Predicting RNA secondary structures with pseudoknots by MCMC sampling

The most probable secondary structure of an RNA molecule, given the nucleotide sequence, can be computed efficiently if a stochastic context-free grammar (SCFG) is used as the prior distribution of the secondary structure. The structures of some RNA molecules contain so-called pseudoknots. Allowing all possible configurations of pseudoknots is not compatible with context-free grammar models and makes the search for an optimal secondary structure NP-complete. We suggest a probabilistic model for RNA secondary structures with pseudoknots and present a Markov-chain Monte-Carlo Method for sampling RNA structures according to their posterior distribution for a given sequence. We favor Bayesian sampling over optimization methods in this context, because it makes the uncertainty of RNA structure predictions assessable. We demonstrate the benefit of our method in examples with tmRNA and also with simulated data. McQFold, an implementation of our method, is freely available from http://www.cs.uni-frankfurt.de/~metzler/McQFold.

Biphasic Folding Kinetics of RNA Pseudoknots and Telomerase RNA Activity

Using a combined master equation and kinetic cluster approach, we investigate RNA pseudoknot folding and unfolding kinetics. The energetic parameters are computed from a recently developed Vfold model for RNA secondary structure and pseudoknot folding thermodynamics. The folding kinetics theory is based on the complete conformational ensemble, including all the native-like and non-native states. The predicted folding and unfolding pathways, activation barriers, Arrhenius plots, and rate-limiting steps lead to several findings. First, for the PK5 pseudoknot, a misfolded 5′ hairpin emerges as a stable kinetic trap in the folding process, and the detrapping from this misfolded state is the rate-limiting step for the overall folding process. The calculated rate constant and activation barrier agree well with the experimental data. Second, as an application of the model, we investigate the kinetic folding pathways for human telomerase RNA (hTR) pseudoknot. The predicted folding and unfolding pathways not only support the proposed role of conformational switch between hairpin and pseudoknot in hTR activity, but also reveal molecular mechanism for the conformational switch. Furthermore, for an experimentally studied hTR mutation, whose hairpin intermediate is destabilized, the model predicts a long-lived transient hairpin structure, and the switch between the transient hairpin intermediate and the native pseudoknot may be responsible for the observed hTR activity. Such finding would help resolve the apparent contradiction between the observed hTR activity and the absence of a stable hairpin.

RNA signals in the 3′ terminus of the genome of Equine arteritis virus are required for viral RNA synthesis

RNA virus genomes contain cis-acting sequences and structural elements involved in virus replication. Both full-length and subgenomic negative-strand RNA synthesis are initiated at the 3' terminus of the positive-strand genomic RNA of Equine arteritis virus (EAV). To investigate the molecular mechanism of EAV RNA synthesis, the RNA secondary structure of the 3'-proximal region of the genome was analysed by chemical and enzymic probing. Based on the RNA secondary structure model derived from this analysis, several deletions were engineered in a full-length cDNA copy of the viral genome. Two RNA domains were identified that are essential for virus replication and most likely play a key role in viral RNA synthesis. The first domain, located directly upstream of the 3' untranslated region (UTR) (nt 12610-12654 of the genome), is mainly single-stranded but contains one small stem-loop structure. The second domain is located within the 3' UTR (nt 12661-12690) and folds into a prominent stem-loop structure with a large loop region. The location of this stem-loop structure near the 3' terminus of the genome suggests that it may act as a recognition signal during the initiation of minus-strand RNA synthesis.

Erratum: Corrigendum: A miniature yeast telomerase RNA functions in vivo and reconstitutes activity in vitro

Nature Structural & Molecular Biology 12, 1072–1077 (2005); published online 20 November 2005; corrected after print 14 April 2006 In the supplementary information initially published online to accompany this article, the legend of Supplementary Figure 1 describes a 500-nt Mini-T RNA secondary structure model, but the model for the 436-nt Mini-T is shown.

Telomerase limits the extent of base pairing between template RNA and telomeric DNA

Telomerase is the ribonucleoprotein reverse transcriptase that adds telomeric DNA repeats to the ends of chromosomes. This involves annealing of the telomerase RNA template to the 3' end of the chromosome, reverse transcription of the RNA template by the telomerase reverse transcriptase polypeptide and translocation. Here, we overexpress and partially purify the catalytically active yeast telomerase core in its natural host and probe telomerase RNA base methylation accessibility with dimethyl sulphate in the presence and absence of a DNA substrate and after substrate elongation. The length of the RNA-DNA hybrid is kept constant at seven base pairs after primer binding and elongation. Thus, new base-pair formation at the 3' end of the substrate during elongation coincides with disruption of base-pair interactions at the other side of the template. Presumably, this circumvents the generation of an exceedingly high energy barrier for translocation and dissociation. Our analysis also corroborates recently proposed yeast telomerase RNA secondary structure models.

Expression of axin2 Is Regulated by the Alternative 5′-Untranslated Regions of Its mRNA

Axin2 is a negative regulator of Wnt/beta-catenin signaling with roles in early development and tumor suppression. We find that axin2 expression is regulated at both transcriptional and translational levels. The gene allows transcription of mRNAs with three alternative 5'-untranslated regions, and these are differentially expressed in various human cell types. These untranslated regions can differentially determine protein expression from messages by influencing mRNA stability and translational efficiency. We identify short upstream reading frames and structural motifs that are responsible for modulation of mRNA translational efficiencies. We show that the proportions of axin2 message expressing each 5'-untranslated region influence the amount of Axin2 protein expressed within cells. We discuss this complex regulation in the context of the function of Axin2 as a tumor suppressor.

Molecular evolution of some selected benthic foraminifera as inferred from sequences of the small subunit ribosomal DNA

Up to now the benthic foraminiferal classification is based mainly on morphological characters leading to controversies and systematic rearrangements. Previous molecular phylogenetic studies using nuclear encoded small subunit (SSU) rDNA as a marker were able to resolve some of the major divergence points. However, they partly contradicted the morphological findings and had difficulties to resolve convincingly the phylogenetic relationships, especially within the multilocular agglutinated and calcitic orders. To further enlarge our knowledge about these phylogenetic relationships, we collected specimens of deep-sea benthic taxa from the Atlantic, the Pacific, the Mediterranean, and the Sea of Marmara, and classified single individuals by a combination of morphological methods (scanning electron microscope (SEM) and light microscope) and molecular tools (sequencing of SSU rDNA). In addition, we developed RNA secondary structure models of the sequenced part of the SSU rDNA, which allowed to determine exactly homologous nucleotide positions and to optimize the alignments used for the phylogenetic analyses. To illustrate the improvement of the methodology, the analyses were carried out together with sequences used in several former studies. The analysis was performed in several steps to find always suitable outgroups for the consecutively following tree construction, starting with the conserved regions for the basal phylogeny and including an increasing amount of variable but alignable regions to improve the resolution of phylogenetically younger groups. Within the limited number of sequenced taxa, the results presented here are, in general, consistent with the fossil record and existing morphological classification. The astrorhizids (agglutinated tests) with the oldest fossil record form the most ancient foraminiferal group. The selected allogromid taxa with organic tests appear to have evolved parallel to a cluster comprising lituolids (agglutinated tests), the porcelaneous miliolids (high magnesium calcite), and multilocular agglutinated and calcitic foraminifers. Furthermore, the calcitic taxa (e.g., buliminids and rotaliids) and the genera Trochammina and Textularia with agglutinated tests appear to be descendants of the stratigraphically older agglutinated foraminifers. Differences to the morphological classification are mainly found at family and genus level within the calcitic foraminifera.