Limited RNA Editing in Exons of Mouse Liver and Adipose

Several studies have investigated RNA–DNA differences (RDD), presumably due to RNA editing, with conflicting results. We report a rigorous analysis of RDD in exonic regions in mice, taking into account critical biases in RNA-Seq analysis. Using deep-sequenced F1 reciprocal inbred mice, we mapped 40 million RNA-Seq reads per liver sample and 180 million reads per adipose sample. We found 7300 apparent hepatic RDDs using a multiple-site mapping procedure, compared with 293 RDD found using a unique-site mapping procedure. After filtering for repeat sequence, splice junction proximity, undirectional strand, and extremity read bias, 63 RDD remained. In adipose tissue unique-site mapping identified 1667 RDD, and after applying the same four filters, 188 RDDs remained. In both tissues, the filtering procedure increased the proportion of canonical (A-to-I and C-to-U) editing events. The genomic DNA of 12 RDD sites among the potential 63 hepatic RDD was tested by Sanger sequencing, three of which proved to be due to unreferenced SNPs. We validated seven liver RDD with Sequenom technology, including two noncanonical, Gm5424 C-to-I(G) and Pisd I(G)-to-A RDD. Differences in diet, sex, or genetic background had very modest effects on RDD occurrence. Only a small number of apparent RDD sites overlapped between liver and adipose, indicating a high degree of tissue specificity. Our findings underscore the importance of properly filtering for bias in RNA-Seq investigations, including the necessity of confirming the DNA sequence to eliminate unreferenced SNPs. Based on our results, we conclude that RNA editing is likely limited to hundreds of events in exonic RNA in liver and adipose.     

An extra double-stranded RNA binding domain confers high activity to a squid RNA editing enzyme

RNA editing by adenosine deamination is particularly prevalent in the squid nervous system. We hypothesized that the squid editing enzyme might contain structural differences that help explain this phenomenon. As a first step, a squid adenosine deaminase that acts on RNA (sqADAR2a) cDNA and the gene that encodes it were cloned from the giant axon system. PCR and RNase protection assays showed that a splice variant of this clone (sqADAR2b) was also expressed in this tissue. Both versions are homologous to the vertebrate ADAR2 family. sqADAR2b encodes a conventional ADAR2 family member with an evolutionarily conserved deaminase domain and two double-stranded RNA binding domains (dsRBD). sqADAR2a differs from sqADAR2b by containing an optional exon that encodes an “extra” dsRBD. Both splice variants are expressed at comparable levels and are extensively edited, each in a unique pattern. Recombinant sqADAR2a and sqADAR2b, produced in Pichia pastoris, are both active on duplex RNA. Using a standard 48-h protein induction, both sqADAR2a and sqADAR2b exhibit promiscuous self-editing; however, this activity is particularly robust for sqADAR2a. By decreasing the induction time to 16 h, self-editing was mostly eliminated. We next tested the ability of sqADAR2a and sqADAR2b to edit two K⁺ channel mRNAs in vitro. Both substrates are known to be edited in squid. For each mRNA, sqADAR2a edited many more sites than sqADAR2b. These data suggest that the “extra” dsRBD confers high activity on sqADAR2a.     

Crystal Structures of the Human RNA Demethylase Alkbh5 Reveal Basis for Substrate Recognition

N⁶-Methylation of adenosine is the most ubiquitous and abundant modification of nucleoside in eukaryotic mRNA and long non-coding RNA. This modification plays an essential role in the regulation of mRNA translation and RNA metabolism. Recently, human AlkB homolog 5 (Alkbh5) and fat mass- and obesity-associated protein (FTO) were shown to erase this methyl modification on mRNA. Here, we report five high resolution crystal structures of the catalytic core of Alkbh5 in complex with different ligands. Compared with other AlkB proteins, Alkbh5 displays several unique structural features on top of the conserved double-stranded β-helix fold typical of this protein family. Among the unique features, a distinct “lid” region of Alkbh5 plays a vital role in substrate recognition and catalysis. An unexpected disulfide bond between Cys-230 and Cys-267 is crucial for the selective binding of Alkbh5 to single-stranded RNA/DNA by bringing a “flipping” motif toward the central β-helix fold. We generated a substrate binding model of Alkbh5 based on a demethylation activity assay of several structure-guided site-directed mutants. Crystallographic and biochemical studies using various analogs of α-ketoglutarate revealed that the active site cavity of Alkbh5 is much smaller than that of FTO and preferentially binds small molecule inhibitors. Taken together, our findings provide a structural basis for understanding the substrate recognition specificity of Alkbh5 and offer a foundation for selective drug design against AlkB members.     

RNA Catalyst as a Reporter for Screening Drugs against RNA Editing in Trypanosomes

Substantial progress has been made in determining the mechanism of mitochondrial RNA editing in trypanosomes. Similarly, considerable progress has been made in identifying the components of the editosome complex that catalyze RNA editing. However, it is still not clear how those proteins work together. Chemical compounds obtained from a high-throughput screen against the editosome may block or affect one or more steps in the editing cycle. Therefore, the identification of new chemical compounds will generate valuable molecular probes for dissecting the editosome function and assembly. In previous studies, in vitro editing assays were carried out using radio-labeled RNA. These assays are time consuming, inefficient and unsuitable for high-throughput purposes. Here, a homogenous fluorescence-based “mix and measure” hammerhead ribozyme in vitro reporter assay to monitor RNA editing, is presented. Only as a consequence of RNA editing of the hammerhead ribozyme a fluorescence resonance energy transfer (FRET) oligoribonucleotide substrate undergoes cleavage. This in turn results in separation of the fluorophore from the quencher thereby producing a signal. In contrast, when the editosome function is inhibited, the fluorescence signal will be quenched. This is a highly sensitive and simple assay that should be generally applicable to monitor in vitro RNA editing or high throughput screening of chemicals that can inhibit the editosome function.     

Correlation between seed longevity and RNA integrity in the embryos of rice seeds

The mature embryos of rice seeds contain translatable mRNAs required for the initial phase of germination. To clarify the relationship between seed longevity and RNA integrity in embryos, germinability and stability of embryonic RNAs were analyzed using the seeds of japonica rice cultivars subjected to controlled deterioration treatment (CDT) or long periods of storage. Degradation of RNA from embryos of a japonica rice cultivar “Nipponbare” was induced by CDT before the decline of the germination rate and we observed a positive relationship between seed germinability and integrity of embryonic RNAs. Moreover, this relationship was confirmed in the experiments using aged seeds from the “Nipponbare”, “Sasanishiki” and “Koshihikari” rice cultivars. In addition, the RNA integrity number (RIN) values, calculated using electrophoresis data and Agilent Bioanalyzer software, had a positive correlation with germinability (R²=0.75). Therefore, the stability of embryonic RNAs required for germination is involved in maintaining seed longevity over time and RIN values can serve as a quantitative indicator to evaluate germinability in rice.     

The RNA structure alignment ontology

Multiple sequence alignments are powerful tools for understanding the structures, functions, and evolutionary histories of linear biological macromolecules (DNA, RNA, and proteins), and for finding homologs in sequence databases. We address several ontological issues related to RNA sequence alignments that are informed by structure. Multiple sequence alignments are usually shown as two-dimensional (2D) matrices, with rows representing individual sequences, and columns identifying nucleotides from different sequences that correspond structurally, functionally, and/or evolutionarily. However, the requirement that sequences and structures correspond nucleotide-by-nucleotide is unrealistic and hinders representation of important biological relationships. High-throughput sequencing efforts are also rapidly making 2D alignments unmanageable because of vertical and horizontal expansion as more sequences are added. Solving the shortcomings of traditional RNA sequence alignments requires explicit annotation of the meaning of each relationship within the alignment. We introduce the notion of “correspondence,” which is an equivalence relation between RNA elements in sets of sequences as the basis of an RNA alignment ontology. The purpose of this ontology is twofold: first, to enable the development of new representations of RNA data and of software tools that resolve the expansion problems with current RNA sequence alignments, and second, to facilitate the integration of sequence data with secondary and three-dimensional structural information, as well as other experimental information, to create simultaneously more accurate and more exploitable RNA alignments.     

Domain-elongation NMR spectroscopy yields new insights into RNA dynamics and adaptive recognition

By simplifying the interpretation of nuclear magnetic resonance spin relaxation and residual dipolar couplings data, recent developments involving the elongation of RNA helices are providing new atomic insights into the dynamical properties that allow RNA structures to change functionally and adaptively. Domain elongation, in concert with spin relaxation measurements, has allowed the detailed characterization of a hierarchical network of local and collective motional modes occurring at nanosecond timescale that mirror the structural rearrangements that take place following adaptive recognition. The combination of domain elongation with residual dipolar coupling measurements has allowed the experimental three-dimensional visualization of very large amplitude rigid-body helix motions in HIV-1 transactivation response element (TAR) that trace out a highly choreographed trajectory in which the helices twist and bend in a correlated manner. The dynamic trajectory allows unbound TAR to sample many of its ligand bound conformations, indicating that adaptive recognition occurs by “conformational selection” rather than “induced fit.” These studies suggest that intrinsic flexibility plays essential roles directing RNA conformational changes along specific pathways.     

Six RNA Viruses and Forty-One Hosts: Viral Small RNAs and Modulation of Small RNA Repertoires in Vertebrate and Invertebrate Systems

We have used multiplexed high-throughput sequencing to characterize changes in small RNA populations that occur during viral infection in animal cells. Small RNA-based mechanisms such as RNA interference (RNAi) have been shown in plant and invertebrate systems to play a key role in host responses to viral infection. Although homologs of the key RNAi effector pathways are present in mammalian cells, and can launch an RNAi-mediated degradation of experimentally targeted mRNAs, any role for such responses in mammalian host-virus interactions remains to be characterized. Six different viruses were examined in 41 experimentally susceptible and resistant host systems. We identified virus-derived small RNAs (vsRNAs) from all six viruses, with total abundance varying from “vanishingly rare” (less than 0.1% of cellular small RNA) to highly abundant (comparable to abundant micro-RNAs “miRNAs”). In addition to the appearance of vsRNAs during infection, we saw a number of specific changes in host miRNA profiles. For several infection models investigated in more detail, the RNAi and Interferon pathways modulated the abundance of vsRNAs. We also found evidence for populations of vsRNAs that exist as duplexed siRNAs with zero to three nucleotide 3′ overhangs. Using populations of cells carrying a Hepatitis C replicon, we observed strand-selective loading of siRNAs onto Argonaute complexes. These experiments define vsRNAs as one possible component of the interplay between animal viruses and their hosts.