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.     

The separation between the 5′-3′ ends in long RNA molecules is short and nearly constant

RNA molecules play different roles in coding, decoding and gene expression regulation. Such roles are often associated to the RNA secondary or tertiary structures. The folding dynamics lead to multiple secondary structures of long RNA molecules, since an RNA molecule might fold into multiple distinct native states. Despite an ensemble of different structures, it has been theoretically proposed that the separation between the 5′ and 3′ ends of long single-stranded RNA molecules (ssRNA) remains constant, independent of their base content and length. Here, we present the first experimental measurements of the end-to-end separation in long ssRNA molecules. To determine this separation, we use single molecule Fluorescence Resonance Energy Transfer of fluorescently end-labeled ssRNA molecules ranging from 500 to 5500 nucleotides in length, obtained from two viruses and a fungus. We found that the end-to-end separation is indeed short, within 5–9 nm. It is remarkable that the separation of the ends of all RNA molecules studied remains small and similar, despite the origin, length and differences in their secondary structure. This implies that the ssRNA molecules are ‘effectively circularized’ something that might be a general feature of RNAs, and could result in fine-tuning for translation and gene expression regulation.     

Defining topological equivalences in macromolecules.

We describe a completely automated and objective method for defining topological equivalents in macromolecules. The method is based on well established techniques for identifying topologically and topographically equivalent atoms in small molecules and has been used in structural alignment of proteins and RNA molecules, and to extract fragments of molecules (protein secondary structures and RNA and DNA double helices) from structural databases consistent with some specified template structure.     

NTP-entry routes in multi-subunit RNA polymerases

The recent elucidation of crystal structures for multi-subunit RNA polymerases immediately revealed a mystery: how do nucleotide triphosphate (NTP) substrates reach an active site that is buried deep within the enzyme? The prevailing view is that NTPs enter through an approximately 20A-long secondary channel between the active site and the enzyme surface. Recently, an alternative view has been advocated; namely, NTPs enter the active site pre-bound to the DNA template from the downstream DNA portion of the main channel of the enzyme.     

RNA channelling by the archaeal exosome

Exosomes are complexes containing 3' --> 5' exoribonucleases that have important roles in processing, decay and quality control of various RNA molecules. Archaeal exosomes consist of a hexameric core of three active RNase PH subunits (ribosomal RNA processing factor (Rrp)41) and three inactive RNase PH subunits (Rrp42). A trimeric ring of subunits with putative RNA-binding domains (Rrp4/cep1 synthetic lethality (Csl)4) is positioned on top of the hexamer on the opposite side to the RNA degrading sites. Here, we present the 1.6 A resolution crystal structure of the nine-subunit exosome of Sulfolobus solfataricus and the 2.3 A structure of this complex bound to an RNA substrate designed to be partly trimmed rather than completely degraded. The RNA binds both at the active site on one side of the molecule and on the opposite side in the narrowest constriction of the central channel. Multiple substrate-binding sites and the entrapment of the substrate in the central channel provide a rationale for the processive degradation of extended RNAs and the stalling of structured RNAs.     

Structural characterization of naturally occurring RNA single mismatches

RNA is known to be involved in several cellular processes; however, it is only active when it is folded into its correct 3D conformation. The folding, bending and twisting of an RNA molecule is dependent upon the multitude of canonical and non-canonical secondary structure motifs. These motifs contribute to the structural complexity of RNA but also serve important integral biological functions, such as serving as recognition and binding sites for other biomolecules or small ligands. One of the most prevalent types of RNA secondary structure motifs are single mismatches, which occur when two canonical pairs are separated by a single non-canonical pair. To determine sequence–structure relationships and to identify structural patterns, we have systematically located, annotated and compared all available occurrences of the 30 most frequently occurring single mismatch-nearest neighbor sequence combinations found in experimentally determined 3D structures of RNA-containing molecules deposited into the Protein Data Bank. Hydrogen bonding, stacking and interaction of nucleotide edges for the mismatched and nearest neighbor base pairs are described and compared, allowing for the identification of several structural patterns. Such a database and comparison will allow researchers to gain insight into the structural features of unstudied sequences and to quickly look-up studied sequences.     

Statistics of RNA secondary structures

A statistical reference for RNA secondary structures with minimum free energies is computed by folding large ensembles of random RNA sequences. Four nucleotide alphabets are used: two binary alphabets, AU and GC, the biophysical AUGC and the synthetic GCXK alphabet. RNA secondary structures are made of structural elements, such as stacks, loops, joints, and free ends. Statistical properties of these elements are computed for small RNA molecules of chain lengths up to 100. The results of RNA structure statistics depend strongly on the particular alphabet chosen. The statistical reference is compared with the data derived from natural RNA molecules with similar base frequencies. Secondary structures are represented as trees. Tree editing provides a quantitative measure for the distance d_t, between two structures. We compute a structure density surface as the conditional probability of two structures having distance t given that their sequences have distance h. This surface indicates that the vast majority of possible minimum free energy secondary structures occur within a fairly small neighborhood of any typical (random) sequence. Correlation lengths for secondary structures in their tree representations are computed from probability densities. They are appropriate measures for the complexity of the sequence-structure relation. The correlation length also provides a quantitative estimate for the mean sensitivity of structures to point mutations. © 1993 John Wiley & Sons, Inc.     

The guide RNA database (3.0).

The RNA editing process within the mitochondria of kinetoplastid organisms is controlled by small, trans -acting RNA molecules referred to as guide RNAs. The guide RNA database is a compilation of published guide RNA sequences, currently containing 254 entries from 11 different organisms. Additional information includes RNA secondary and tertiary structure models, information on the gene localisation, literature citations and other relevant facts. The database can be accessed through the World Wide Web (WWW) at http://www.biochem.mpg.de/ goeringe/     

In vivo analysis of plant RNA structure: soybean 18S ribosomal and ribulose-1,5-bisphosphate carboxylase small subunit RNAs

A method to investigate the structure of RNA molecules within intact plant tissues has been developed. The RNA structures are analyzed using dimethyl sulfate (DMS), which modifies substituents of adenine and cytosine residues within single-stranded regions of RNA molecules. Reactive sites are identified by primer extension analysis. Using this procedure, an analysis of the secondary structure of the cytoplasmic 18S ribosomal RNA in soybean seedling leaves has been completed. DMS modification data are in good agreement with the phylogenetic structure predicted for soybean 18S rRNA. However, there are a few notable exceptions where residues thought to be involved in double-stranded regions in all 18S rRNAs are strongly modified in soybean leaf samples. These data taken together with the phylogenetic structure suggest that alternate structures may exist in vivo.The further applicability of this technique is demonstrated by comparing the modification pattern obtained in vivo to that obtained in vitro for a particular mRNA molecule encoding the small subunit of ribulose-1,5-bisphosphate carboxylase. The results obtained are compared to a predicted minimum energy secondary structure. The data indicate that the conformation of RNA molecules within the cell may not be reflected in a structural analysis of purified mRNA molecules.     

RNA STRAND: The RNA Secondary Structure and Statistical Analysis Database

Background  

The ability to access, search and analyse secondary structures of a large set of known RNA molecules is very important for deriving improved RNA energy models, for evaluating computational predictions of RNA secondary structures and for a better understanding of RNA folding. Currently there is no database that can easily provide these capabilities for almost all RNA molecules with known secondary structures.     

The guide RNA database.

Guide RNAs (gRNAs) are small, metabolically stable RNA molecules which perform a pivotal, template-like function during the RNA editing process in kinetoplastid protozoa. The gRNA database currently contains 250 guide RNA sequences as well as secondary and tertiary structure models and other relevant information. The database is made available as a hypertext document accessible via the World Wide Web (WWW) at the URL: http://www.biochem.mpg.de/ goeringe/     

Parallel chemical genetic and genome-wide RNAi screens identify cytokinesis inhibitors and targets

Cytokinesis involves temporally and spatially coordinated action of the cell cycle and cytoskeletal and membrane systems to achieve separation of daughter cells. To dissect cytokinesis mechanisms it would be useful to have a complete catalog of the proteins involved, and small molecule tools for specifically inhibiting them with tight temporal control. Finding active small molecules by cell-based screening entails the difficult step of identifying their targets. We performed parallel chemical genetic and genome-wide RNA interference screens in Drosophila cells, identifying 50 small molecule inhibitors of cytokinesis and 214 genes important for cytokinesis, including a new protein in the Aurora B pathway (Borr). By comparing small molecule and RNAi phenotypes, we identified a small molecule that inhibits the Aurora B kinase pathway. Our protein list provides a starting point for systematic dissection of cytokinesis, a direction that will be greatly facilitated by also having diverse small molecule inhibitors, which we have identified. Dissection of the Aurora B pathway, where we found a new gene and a specific small molecule inhibitor, should benefit particularly. Our study shows that parallel RNA interference and small molecule screening is a generally useful approach to identifying active small molecules and their target pathways.