Structure and mechanism of E. coli RNA 2′,3′-cyclic phosphodiesterase

2H (two-histidine) phosphoesterase enzymes are distributed widely in all domains of life and are implicated in diverse RNA and nucleotide transactions, including the transesterification and hydrolysis of cyclic phosphates. Here we report a biochemical and structural characterization of the Escherichia coli 2H protein YapD, which was identified originally as a reversible transesterifying “nuclease/ligase” at RNA 2′,5′-phosphodiesters. We find that YapD is an “end healing” cyclic phosphodiesterase (CPDase) enzyme that hydrolyzes an _HORNA>p substrate with a 2′,3′-cyclic phosphodiester to a _HORNAp product with a 2′-phosphomonoester terminus, without concomitant end joining. Thus we rename this enzyme ThpR (two-histidine 2′,3′-cyclic phosphodiesterase acting on RNA). The 2.0 Å crystal structure of ThpR in a product complex with 2′-AMP highlights the roles of extended histidine-containing motifs ⁴³HxTxxF⁴⁸ and ¹²⁵HxTxxR¹³⁰ in the CPDase reaction. His43-Nε makes a hydrogen bond with the ribose O3′ leaving group, thereby implicating His43 as a general acid catalyst. His125-Nε coordinates the O1P oxygen of the AMP 2′-phosphate (inferred from geometry to derive from the attacking water nucleophile), pointing to His125 as a general base catalyst. Arg130 makes bidentate contact with the AMP 2′-phosphate, suggesting a role in transition-state stabilization. Consistent with these inferences, changing His43, His125, or Arg130 to alanine effaced the CPDase activity of ThpR. Phe48 makes a π–π stack on the adenine nucleobase. Mutating Phe28 to alanine slowed the CPDase by an order of magnitude. The tertiary structure and extended active site motifs of ThpR are conserved in a subfamily of bacterial and archaeal 2H enzymes.     

Further Insights into the Roles of GTP and the C Terminus of the Hepatitis C Virus Polymerase in the Initiation of RNA Synthesis

The hepatitis C virus (HCV) NS5b protein is an RNA-dependent RNA polymerase essential for replication of the viral RNA genome. In vitro and presumably in vivo, NS5b initiates RNA synthesis by a de novo mechanism. Different structural elements of NS5b have been reported to participate in RNA synthesis, especially a so-called “β-flap” and a C-terminal segment (designated “linker”) that connects the catalytic core of NS5b to a transmembrane anchor. High concentrations of GTP have also been shown to stimulate de novo RNA synthesis by HCV NS5b. Here we describe a combined structural and functional analysis of genotype 1 HCV-NS5b of strains H77 (subtype 1a), for which no structure has been previously reported, and J4 (subtype 1b). Our results highlight the linker as directly involved in lifting the first boundary to processive RNA synthesis, the formation of the first dinucleotide primer. The transition from this first dinucleotide primer state to processive RNA synthesis requires removal of the linker and of the β-flap with which it is shown to strongly interact in crystal structures of HCV NS5b. We find that GTP specifically stimulates this transition irrespective of its incorporation in neosynthesized RNA.     

Captured Segment Exchange: A Strategy for Custom Engineering Large Genomic Regions in Drosophila melanogaster

Site-specific recombinases (SSRs) are valuable tools for manipulating genomes. In Drosophila, thousands of transgenic insertions carrying SSR recognition sites have been distributed throughout the genome by several large-scale projects. Here we describe a method with the potential to use these insertions to make custom alterations to the Drosophila genome in vivo. Specifically, by employing recombineering techniques and a dual recombinase-mediated cassette exchange strategy based on the phiC31 integrase and FLP recombinase, we show that a large genomic segment that lies between two SSR recognition-site insertions can be “captured” as a target cassette and exchanged for a sequence that was engineered in bacterial cells. We demonstrate this approach by targeting a 50-kb segment spanning the tsh gene, replacing the existing segment with corresponding recombineered sequences through simple and efficient manipulations. Given the high density of SSR recognition-site insertions in Drosophila, our method affords a straightforward and highly efficient approach to explore gene function in situ for a substantial portion of the Drosophila genome.     

Benefits and limitations of a new genome‐based PCR‐RFLP genotyping assay (GB‐RFLP): A SNP‐based detection method for identification of species in extremely young adaptive radiations

High‐throughput DNA sequencing technologies make it possible now to sequence entire genomes relatively easily. Complete genomic information obtained by whole‐genome resequencing (WGS) can aid in identifying and delineating species even if they are extremely young, cryptic, or morphologically difficult to discern and closely related. Yet, for taxonomic or conservation biology purposes, WGS can remain cost‐prohibitive, too time‐consuming, and often constitute a “data overkill.” Rapid and reliable identification of species (and populations) that is also cost‐effective is made possible by species‐specific markers that can be discovered by WGS. Based on WGS data, we designed a PCR restriction fragment length polymorphism (PCR‐RFLP) assay for 19 Neotropical Midas cichlid populations (Amphilophus cf. citrinellus), that includes all 13 described species of this species complex. Our work illustrates that identification of species and populations (i.e., fish from different lakes) can be greatly improved by designing genetic markers using available “high resolution” genomic information. Yet, our work also shows that even in the best‐case scenario, when whole‐genome resequencing information is available, unequivocal assignments remain challenging when species or populations diverged very recently, or gene flow persists. In summary, we provide a comprehensive workflow on how to design RFPL markers based on genome resequencing data, how to test and evaluate their reliability, and discuss the benefits and pitfalls of our approach.     

Rapid and Precise Engineering of the Caenorhabditis elegans Genome with Lethal Mutation Co-Conversion and Inactivation of NHEJ Repair

As in other organisms, CRISPR/Cas9 methods provide a powerful approach for genome editing in the nematode Caenorhabditis elegans. Oligonucleotides are excellent repair templates for introducing substitutions and short insertions, as they are cost effective, require no cloning, and appear in other organisms to target changes by homologous recombination at DNA double-strand breaks (DSBs). Here, I describe a methodology in C. elegans to efficiently knock in epitope tags in 8–9 days, using a temperature-sensitive lethal mutation in the pha-1 gene as a co-conversion marker. I demonstrate that 60mer oligos with 29 bp of homology drive efficient knock-in of point mutations, and that disabling nonhomologous end joining by RNAi inactivation of the cku-80 gene significantly improves knock-in efficiency. Homology arms of 35–80 bp are sufficient for efficient editing and DSBs up to 54 bp away from the insertion site produced knock-ins. These findings will likely be applicable for a range of genome editing approaches in C. elegans, which will improve editing efficiency and minimize screening efforts.     

Sliced microRNA targets and precise loop-first processing of MIR319 hairpins revealed by analysis of the Physcomitrella patens degradome

Expression profiling of the 5′ ends of uncapped mRNAs (“degradome” sequencing) can be used to empirically catalog microRNA (miRNA) targets, to probe patterns of miRNA hairpin processing, to examine mRNA decay, and to analyze accumulation of endogenous short interfering RNA (siRNA) precursors. We sequenced and analyzed the degradome of the moss Physcomitrella patens, an important model system for functional genomic analyses in plant evolution. A total of 52 target mRNAs of 27 different Physcomitrella miRNA families were identified. Many targets of both more conserved and less conserved miRNA families encoded putative regulatory proteins. Remnants of MIRNA hairpin processing also populated the degradome data and indicated an unusual “loop-first” mode of precise processing for the MIR319 gene family. Precise loop-first processing was confirmed for native Physcomitrella, rice, and Arabidopsis MIR319 hairpins, as well as an Arabidopsis artificial MIRNA (aMIRNA) based upon a MIR319 backbone. MIR319 is thus a conserved exception to the general rule of loop-last processing of MIRNA hairpins. Loop-first MIR319 processing may contribute to the high efficacy of a widely used MIR319-based strategy for aMIRNA production in plants.     

Circular RNAs are abundant,conserved, and associated with ALU repeats

Circular RNAs composed of exonic sequence have been described in a small number of genes. Thought to result from splicing errors, circular RNA species possess no known function. To delineate the universe of endogenous circular RNAs, we performed high-throughput sequencing (RNA-seq) of libraries prepared from ribosome-depleted RNA with or without digestion with the RNA exonuclease, RNase R. We identified >25,000 distinct RNA species in human fibroblasts that contained non-colinear exons (a “backsplice”) and were reproducibly enriched by exonuclease degradation of linear RNA. These RNAs were validated as circular RNA (ecircRNA), rather than linear RNA, and were more stable than associated linear mRNAs in vivo. In some cases, the abundance of circular molecules exceeded that of associated linear mRNA by >10-fold. By conservative estimate, we identified ecircRNAs from 14.4% of actively transcribed genes in human fibroblasts. Application of this method to murine testis RNA identified 69 ecircRNAs in precisely orthologous locations to human circular RNAs. Of note, paralogous kinases HIPK2 and HIPK3 produce abundant ecircRNA from their second exon in both humans and mice. Though HIPK3 circular RNAs contain an AUG translation start, it and other ecircRNAs were not bound to ribosomes. Circular RNAs could be degraded by siRNAs and, therefore, may act as competing endogenous RNAs. Bioinformatic analysis revealed shared features of circularized exons, including long bordering introns that contained complementary ALU repeats. These data show that ecircRNAs are abundant, stable, conserved and nonrandom products of RNA splicing that could be involved in control of gene expression.     

A long-range RNA–RNA interaction between the 5′ and 3′ ends of the HCV genome

The RNA genome of the hepatitis C virus (HCV) contains multiple conserved structural cis domains that direct protein synthesis, replication, and infectivity. The untranslatable regions (UTRs) play essential roles in the HCV cycle. Uncapped viral RNAs are translated via an internal ribosome entry site (IRES) located at the 5′ UTR, which acts as a scaffold for recruiting multiple protein factors. Replication of the viral genome is initiated at the 3′ UTR. Bioinformatics methods have identified other structural RNA elements thought to be involved in the HCV cycle. The 5BSL3.2 motif, which is embedded in a cruciform structure at the 3′ end of the NS5B coding sequence, contributes to the three-dimensional folding of the entire 3′ end of the genome. It is essential in the initiation of replication. This paper reports the identification of a novel, strand-specific, long-range RNA–RNA interaction between the 5′ and 3′ ends of the genome, which involves 5BSL3.2 and IRES motifs. Mutants harboring substitutions in the apical loop of domain IIId or in the internal loop of 5BSL3.2 disrupt the complex, indicating these regions are essential in initiating the kissing interaction. No complex was formed when the UTRs of the related foot and mouth disease virus were used in binding assays, suggesting this interaction is specific for HCV sequences. The present data firmly suggest the existence of a higher-order structure that may mediate a protein-independent circularization of the HCV genome. The 5′–3′ end bridge may have a role in viral translation modulation and in the switch from protein synthesis to RNA replication.     

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.     

Thermodynamic examination of 1- to 5-nt purine bulge loops in RNA and DNA constructs

Bulge loops are common features of RNA structures that are involved in the formation of RNA tertiary structures and are often sites for interactions with proteins and ions. Minimal thermodynamic data currently exist on the bulge size and sequence effects. Using thermal denaturation methods, thermodynamic properties of 1- to 5-nt adenine and guanine bulge loop constructs were examined in 10 mM MgCl₂ or 1 M KCl. The ΔG37∘ loop parameters for 1- to 5-nt purine bulge loops in RNA constructs were between 3.07 and 5.31 kcal/mol in 1 M KCl buffer. In 10 mM magnesium ions, the ΔΔG° values relative to 1 M KCl were 0.47–2.06 kcal/mol more favorable for the RNA bulge loops. The ΔG37∘ loop parameters for 1- to 5-nt purine bulge loops in DNA constructs were between 4.54 and 5.89 kcal/mol. Only 4- and 5-nt guanine constructs showed significant change in stability for the DNA constructs in magnesium ions. A linear correlation is seen between the size of the bulge loop and its stability. New prediction models are proposed for 1- to 5-nt purine bulge loops in RNA and DNA in 1 M KCl. We show that a significant stabilization is seen for small bulge loops in RNA in the presence of magnesium ions. A prediction model is also proposed for 1- to 5-nt purine bulge loop RNA constructs in 10 mM magnesium chloride.     

“Protein aggregates” contain RNA and DNA,entrapped by misfolded proteins but largely rescued by slowing translational elongation

All neurodegenerative diseases feature aggregates, which usually contain disease‐specific diagnostic proteins; non‐protein constituents, however, have rarely been explored. Aggregates from SY5Y‐APP_Sw neuroblastoma, a cell model of familial Alzheimer''s disease, were crosslinked and sequences of linked peptides identified. We constructed a normalized “contactome” comprising 11 subnetworks, centered on 24 high‐connectivity hubs. Remarkably, all 24 are nucleic acid‐binding proteins. This led us to isolate and sequence RNA and DNA from Alzheimer''s and control aggregates. RNA fragments were mapped to the human genome by RNA‐seq and DNA by ChIP‐seq. Nearly all aggregate RNA sequences mapped to specific genes, whereas DNA fragments were predominantly intergenic. These nucleic acid mappings are all significantly nonrandom, making an artifactual origin extremely unlikely. RNA (mostly cytoplasmic) exceeded DNA (chiefly nuclear) by twofold to fivefold. RNA fragments recovered from AD tissue were ~1.5‐to 2.5‐fold more abundant than those recovered from control tissue, similar to the increase in protein. Aggregate abundances of specific RNA sequences were strikingly differential between cultured SY5Y‐APP_Sw glioblastoma cells expressing APOE3 vs. APOE4, consistent with APOE4 competition for E‐box/CLEAR motifs. We identified many G‐quadruplex and viral sequences within RNA and DNA of aggregates, suggesting that sequestration of viral genomes may have driven the evolution of disordered nucleic acid‐binding proteins. After RNA‐interference knockdown of the translational‐procession factor EEF2 to suppress translation in SY5Y‐APP_Sw cells, the RNA content of aggregates declined by >90%, while reducing protein content by only 30% and altering DNA content by ≤10%. This implies that cotranslational misfolding of nascent proteins may ensnare polysomes into aggregates, accounting for most of their RNA content.