Evolution of secondary structure in the family of 7SL-like RNAs

Primate and rodent genomes are populated with hundreds of thousands copies of Alu and B1 elements dispersed by retroposition, i.e., by genomic reintegration of their reverse transcribed RNAs. These, as well as primate BC200 and rodent 4.5S RNAs, are ancestrally related to the terminal portions of 7SL RNA sequence. The secondary structure of 7SL RNA (an integral component of the signal recognition particle) is conserved from prokaryotes to distant eukaryotic species. Yet only in primates and rodents did this molecule give rise to retroposing Alu and B1 RNAs and to apparently functional BC200 and 4.5S RNAs. To understand this transition and the underlying molecular events, we examined, by comparative analysis, the evolution of RNA structure in this family of molecules derived from 7SL RNA.RNA sequences of different simian (mostly human) and prosimian Alu subfamilies as well as rodent B1 repeats were derived from their genomic consensus sequences taken from the literature and our unpublished results (prosimian and New World Monkey). RNA secondary structures were determined by enzymatic studies (new data on 4.5S RNA are presented) and/or energy minimization analyses followed by phylogenetic comparison. Although, with the exception of 4.5S RNA, all 7SL-derived RNA species maintain the cruciform structure of their progenitor, the details of 7SL RNA folding domains are modified to a different extent in various RNA groups. Novel motifs found in retropositionally active RNAs are conserved among Alu and B1 subfamilies in different genomes. In RNAs that do not proliferate by retroposition these motifs are modified further. This indicates structural adaptation of 7SL-like RNA molecules to novel functions, presumably mediated by specific interactions with proteins; these functions were either useful for the host or served the selfish propagation of RNA templates within the host genome.Abbreviations FAM fossil Alu element - FLAM free left Alu monomer - FRAM free right Alu monomer - L-Alu left Alu subunit - R-Alu right Alu subunit Correspondence to: D. LabudaDedicated to Dr. Robert Cedergren on the occasion of his 25th anniversary at the University of Montreal     

Expression of B7-H1 and B7-DC on the airway epithelium is enhanced by double-stranded RNA

Viral infection in the airway provokes various immune responses, including Th1 and Th2 responses, which are partly initiated by double-stranded RNA (dsRNA), a viral product for its replication. B7-H1 (PD-L1) and B7-DC (PD-L2) are B7-family molecules that bind to programmed death-1 (PD-1) on lymphocytes and are implicated in peripheral tolerance. We investigated the effect of dsRNA on the expression of B7-H1 and B7-DC on airway epithelial cell lines. B7-H1 and B7-DC were constitutively expressed on the cells, and their expression was profoundly upregulated by stimulation with an analog of viral dsRNA, polyinosinic-polycytidylic acid. B7-H1 and B7-DC were also upregulated by stimulation with IFN-gamma, IL-13, and the supernatant from T cell clones. A relatively high concentration of dexamethasone (1 microM) was required to suppress the upregulation of B7-H1 or B7-DC. These results suggest that epithelial B7-H1 and B7-DC play a role in virus-associated immune responses in the airways.     

Successive waves of fixation of B1 variants in rodent lineage history

Summary A new method of analyzing phylogenetic relations among members of sequence family (Quentin 1988) discriminates at least six possible B1 subfamilies in the mouse genome. Several additional and independent observations suggest that these grouping have evolutionary significance, and that successive waves of fixation of new variants occur during rodent lineage history. We have reason to believe that, in a genome, the founder sequences of different families of retroposons are in competition with regard to the amplification/fixation process.     

正在崛起的RNA纳米技术

RNA纳米技术得益于纽约大学西曼(Nadrian C.Seeman)教授开创的DNA纳米技术,RNA是由腺嘌呤(A)、尿嘧啶(U)、鸟嘌呤(G)和胞嘧啶(C)构成的一种核糖核酸高分子,与DNA的Watson-Crick碱基配对(A-T,G-C)的双螺旋链的结构不完全一样,RNA的二级结构里经常出现一些非传统的碱基配对如环环相互作用,这些非传统配对促使RNA分子折叠成刚性结构。本文综述了正在崛起的RNA纳米技术,列举了一些著名的实验,如郭培宣(Peixuan Guo)等从自然界的phi29噬菌体中发现的pRNA纳米马达是由六个小RNA分子构成的六环结构,Jaeger等发展了RNA构造术(RNA-tectonics),根据已知的RNA分子的碱基和非传统配对,他们设计利用小RNA分子构造二聚体、一维线性多聚体、和二维网状的七巧板迷宫(jigsaw puzzle)等图案,用tRNA分子或设计用几条RNA分子来构建多面体如立方体和八面体等立体结构等。RNA纳米技术正在崛起,它将在医学、生物技术、合成生物学和纳米技术领域扮演重要的角色。     

Evolution of B2 repeats: the muroid explosion

B2 repeats are a group of short interspersed elements (SINEs) specific for rodent genomes. Copy numbers were determined for different rodent genera. All the Muroid (rat, mouse, deer mouse, hamster, gerbil) rodent genomes analyzed exhibited 80,000–100,000 copies per haploid genome, whereas the squirrel genome contains only 2,500 copies, and fewer than 100 (if any) copies were observed for the Hystricognath rodents (guinea pig and nutria). These findings demonstrate that there was an explosion of amplification of B2 elements within muroid rodents. The similar copy number of B2 elements within the different muroid species could be explained by formation of a high proportion of the B2 elements prior to the divergence of the different muroid species. However, the 3-end of the B2 sequence is unique between murid and cricetid rodents suggesting that the majority of elements amplified after the divergence of these species. Also consistent with recent amplification of these elements in parallel within the muroid genomes is the finding that within mouse and rat there are distinct subfamilies of B2 repeats. The pattern of consistent parallel amplification of B2 elements in muroid species contrasts with the sporadic nature of ID repeat amplification in the same genomes. The consensus of the young mouse subfamily of elements corresponds to the B2 RNA that is preferentially transcribed in embryonic, tumor, and normal liver cells. The subfamily is young based on both its low divergence from the subfamily consensus sequence and the finding that the most recent B2 element insertions in the mouse genome are members of this subfamily.     

ExpaRNA-P: simultaneous exact pattern matching and folding of RNAs

Background

Identifying sequence-structure motifs common to two RNAs can speed up the comparison of structural RNAs substantially. The core algorithm of the existent approach ExpaRNA solves this problem for a priori known input structures. However, such structures are rarely known; moreover, predicting them computationally is no rescue, since single sequence structure prediction is highly unreliable.

Results

The novel algorithm ExpaRNA-P computes exactly matching sequence-structure motifs in entire Boltzmann-distributed structure ensembles of two RNAs; thereby we match and fold RNAs simultaneously, analogous to the well-known “simultaneous alignment and folding” of RNAs. While this implies much higher flexibility compared to ExpaRNA, ExpaRNA-P has the same very low complexity (quadratic in time and space), which is enabled by its novel structure ensemble-based sparsification. Furthermore, we devise a generalized chaining algorithm to compute compatible subsets of ExpaRNA-P’s sequence-structure motifs. Resulting in the very fast RNA alignment approach ExpLoc-P, we utilize the best chain as anchor constraints for the sequence-structure alignment tool LocARNA. ExpLoc-P is benchmarked in several variants and versus state-of-the-art approaches. In particular, we formally introduce and evaluate strict and relaxed variants of the problem; the latter makes the approach sensitive to compensatory mutations. Across a benchmark set of typical non-coding RNAs, ExpLoc-P has similar accuracy to LocARNA but is four times faster (in both variants), while it achieves a speed-up over 30-fold for the longest benchmark sequences (≈400nt). Finally, different ExpLoc-P variants enable tailoring of the method to specific application scenarios. ExpaRNA-P and ExpLoc-P are distributed as part of the LocARNA package. The source code is freely available at http://www.bioinf.uni-freiburg.de/Software/ExpaRNA-P.

Conclusions

ExpaRNA-P’s novel ensemble-based sparsification reduces its complexity to quadratic time and space. Thereby, ExpaRNA-P significantly speeds up sequence-structure alignment while maintaining the alignment quality. Different ExpaRNA-P variants support a wide range of applications.

Electronic supplementary material

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

The p53 target protein Wig-1 binds hnRNP A2/B1 and RNA Helicase A via RNA

The p53-induced Wig-1 gene encodes a double stranded RNA-binding zinc finger protein. We generated Saos-2 osteosarcoma cells expressing tetracycline-inducible Flag-tagged human Wig-1. Induction of Wig-1 expression by doxycycline inhibited cell growth in a long-term assay but did not cause any changes in cell cycle distribution nor increased fraction of apoptotic cells. Using co-immunoprecipitation and mass spectrometry, we identified two Wig-1-binding proteins, hnRNP A2/B1 and RNA Helicase A, both of which are involved in RNA processing. The binding was dependent on the presence of RNA. Our results establish a link between the p53 tumor suppressor and RNA processing via hnRNPA2/B1 and RNA Helicase A.     

RNA2D3D: A program for Generating,Viewing, and Comparing 3-Dimensional Models of RNA

Abstract

Using primary and secondary structure information of an RNA molecule, the program RNA2D3D automatically and rapidly produces a first-order approximation of a 3-dimensional conformation consistent with this information. Applicable to structures of arbitrary branching complexity and pseudoknot content, it features efficient interactive graphical editing for the removal of any overlaps introduced by the initial generating procedure and for making conformational changes favorable to targeted features and subsequent refinement. With emphasis on fast exploration of alternative 3D conformations, one may interactively add or delete base-pairs, adjacent stems can be coaxially stacked or unstacked, single strands can be shaped to accommodate special constraints, and arbitrary subsets can be defined and manipulated as rigid bodies. Compaction, whereby base stacking within stems is optimally extended into connecting single strands, is also available as a means of strategically making the structures more compact and revealing folding motifs. Subsequent refinement of the first-order approximation, of modifications, and for the imposing of tertiary constraints is assisted with standard energy refinement techniques. Previously determined coordinates for any part of the molecule are readily incorporated, and any part of the modeled structure can be output as a PDB or XYZ file. Illustrative applications in the areas of ribozymes, viral kissing loops, viral internal ribosome entry sites, and nanobiology are presented.     

人B7—1和B7—2cDNA的克隆及鉴定

目的：为探索性构建全新型重组人B7－PE40绿脓杆菌外毒素融合蛋白以长期诱导免疫耐受,本研究从急性B淋巴细胞白血病细胞株Raji中隆N－末端分别缺失34和16个氨基酸的人B7－1和B7－2基因胞外区,并构建含此基因的重组质粒。方法：根据B7－1和B7－2基因序列设计合成可增B7－1和B7－2cDNA的特异性引物,用RT－PCR的方法从Raji细胞总RNA中扩增B7－1和B7－2cDNA,并克隆至pGEM－T载体中,经酶切鉴定后再进行序列分析。结果和结论：从Raji细胞中扩增出预期 的624和675bp的B7－1和B7－2cDNA,将其克隆至pGEM-T载体中,分别经EcoRI/HindⅢ和BamHI/SphI双酶切电泳和序列分析确证,为进一步构建人B7－PE40外毒素融合蛋白奠定了基础。     

Improving the prediction of RNA secondary structure by detecting and assessing conserved stems

The prediction of RNA secondary structure can be facilitated by incorporating with comparative analysis of homologous sequences. However, most of existing comparative methods are vulnerable to alignment errors and thus are of low accuracy in practical application. Here we improve the prediction of RNA secondary structure by detecting and assessing conserved stems shared by all sequences in the alignment. Our method can be summarized by: 1) we detect possible stems in single RNA sequence using the so-called position matrix with which some possibly paired positions can be uncovered; 2) we detect conserved stems across multiple RNA sequences by multiplying the position matrices; 3) we assess the conserved stems using the Signal-to-Noise; 4) we compute the optimized secondary structure by incorporating the so-called reliable conserved stems with predictions by RNAalifold program. We tested our method on data sets of RNA alignments with known secondary structures. The accuracy, measured as sensitivity and specificity, of our method is greater than predictions by RNAalifold.     

Fragile X mental retardation protein interactions with the microtubule associated protein 1B RNA

Fragile X mental retardation syndrome, the most common form of inherited mental retardation, is caused by the absence of the fragile X mental retardation protein (FMRP). FMRP has been shown to use its arginine-glycine-glycine (RGG) box to bind to a subset of RNA targets that form a G quadruplex structure. We performed a detailed analysis of the interactions between the FMRP RGG box and the microtubule associated protein 1B (MAP1B) mRNA, a relevant in vivo FMRP target. We show that MAP1B RNA forms an intramolecular G quadruplex structure, which is bound with high affinity and specificity by the FMRP RGG box. We determined that hydrophobic interactions are important in the FMRP RGG box-MAP1B RNA association, with minor contributions from electrostatic interactions. Our findings that at low protein:RNA ratios the RNA G quadruplex structure is slightly stabilized, whereas at high ratios is unfolded, suggest a mechanism by which the FMRP concentration variation in response to a neurotransmitter stimulation event could act as a regulatory switch for the protein function, from translation repressor to translation activator.     

Insight into the intrinsic flexibility of the SL1 stem-loop from genomic RNA of HIV-1 as probed by molecular dynamics simulation

The SL1 stem-loop is the dimerization initiation site for linking the two copies of the RNA forming the HIV-1 genome. The 26 nucleotides stem contains a defect consisting on a highly conserved G-rich 1-3 asymmetrical internal loop, which is a major site for nucleocapsid protein binding. Several NMR attempts were undertaken to determine the internal loop structure in the SL1 monomer. However, the RNA constructs used in the different studies were largely mutated, in particular with replacement of the nine nucleotides apical loop by a tetraloop, and divergent results were obtained ranging from a rigid structure to a particularly large flexibility. To investigate the reasons for such discrepancies, we used molecular dynamics simulation of the SL1 monomer to probe the effect of mutations on the conformational stability of the internal loop and of the whole stem. It is found that in the wild-type sequence, the internal loop displays conformational variability originating mainly from the nine nucleotides apical loop flexibility that causes large conformational fluctuations (without changing the average structure) in the 7 bp duplex linking the apical and internal loops. The large amplitude atomic motions in the duplex are transmitted to the internal loop in which they induce conformational changes characterized by a labile hydrogen bond network such as G5 successively H-bonded to A29 and G30. On the contrary, with a four nucleotides apical loop, conformational fluctuations in the duplex are reduced by a factor of 2 and are not sufficiently energizing for promoting changes in the internal loop structure at the time scale of the simulations.     

RNase J1 endonuclease activity as a probe of RNA secondary structure

Reliable determination of RNA secondary structure depends on both computer algorithms and experimental probing of nucleotides in single- or double-stranded conformation. Here we describe the exploitation of the endonucleolytic activity of the Bacillus subtilis enzyme RNase J1 as a probe of RNA structure. RNase J1 cleaves in single-stranded regions and, in vitro at least, the enzyme has relatively relaxed nucleotide specificity. We confirmed the feasibility of the approach on an RNA of known structure, B. subtilis tRNA^Thr. We then used RNase J1 to solve the secondary structure of the 5′ end of the hbs mRNA. Finally, we showed that RNase J1 can also be used in footprinting experiments by probing the interaction between the 30S ribosomal subunit and the Shine–Dalgarno element of the hbs mRNA.     

Translational standby sites: how ribosomes may deal with the rapid folding kinetics of mRNA

We have previously shown that stable base-pairing at a translational initiation site in Escherichia coli can inhibit translation by competing with the binding of ribosomes. When the base-pairing is not too strong, this competition is won by the ribosomes, resulting in efficient translation from a structured ribosome binding site (RBS). We now re-examine these results in the light of RNA folding kinetics and find that the window during which a folded RBS is open is generally much too short to recruit a 30S ribosomal subunit from the cytoplasm. We argue that to achieve efficient expression, a 30S subunit must already be in contact with the mRNA while this is still folded, to shift into place as soon as the structure opens. Single-stranded regions flanking the structure may constitute a standby site, to which the 30S subunit can attach non-specifically. We propose a steady-state kinetic model for the early steps of translational initiation and use this to examine various quantitative aspects of standby binding. The kinetic model provides an explanation of why the earlier equilibrium competition model predicted implausibly high 30S-mRNA affinities. Because all RNA is structured to some degree, standby binding is probably a general feature of translational initiation.