Search

RNA guanine quadruplexes (rG4) assume important roles in post-transcriptional regulations of gene expression, which are often modulated by rG4-binding proteins. Hence, understanding the biological functions of rG4s requires the identification and functional characterizations of rG4-recognition proteins. By employing a bioinformatic approach based on the analysis of overlap between peaks obtained from rG4-seq analysis and those detected in >230 eCLIP-seq datasets for RNA-binding proteins generated from the ENCODE project, we identified a large number of candidate rG4-binding proteins. We showed that one of these proteins, G3BP1, is able to bind directly to rG4 structures with high affinity and selectivity, where the binding entails its C-terminal RGG domain and is further enhanced by its RRM domain. Additionally, our seCLIP-Seq data revealed that pyridostatin, a small-molecule rG4 ligand, could displace G3BP1 from mRNA in cells, with the most pronounced effects being observed for the 3-untranslated regions (3-UTR) of mRNAs. Moreover, luciferase reporter assay results showed that G3BP1 positively regulates mRNA stability through its binding with rG4 structures. Together, we identified a number of candidate rG4-binding proteins and validated that G3BP1 can bind directly with rG4 structures and regulate the stabilities of mRNAs.

Diversity-generating retroelements (DGRs) vary protein sequences to the greatest extent known in the natural world. These elements are encoded by constituents of the human microbiome and the microbial 'dark matter'. Variation occurs through adenine-mutagenesis, in which genetic information in RNA is reverse transcribed faithfully to cDNA for all template bases but adenine. We investigated the determinants of adenine-mutagenesis in the prototypical Bordetella bacteriophage DGR through an in vitro system composed of the reverse transcriptase bRT, Avd protein, and a specific RNA. We found that the catalytic efficiency for correct incorporation during reverse transcription by the bRT-Avd complex was strikingly low for all template bases, with the lowest occurring for adenine. Misincorporation across a template adenine was only somewhat lower in efficiency than correct incorporation. We found that the C6, but not the N1 or C2, purine substituent was a key determinant of adenine-mutagenesis. bRT-Avd was insensitive to the C6 amine of adenine but recognized the C6 carbonyl of guanine. We also identified two bRT amino acids predicted to nonspecifically contact incoming dNTPs, R74 and I181, as promoters of adenine-mutagenesis. Our results suggest that the overall low catalytic efficiency of bRT-Avd is intimately tied to its ability to carry out adenine-mutagenesis.

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences that provides a single access point to 44 RNA resources and >18 million ncRNA sequences from a wide range of organisms and RNA types. RNAcentral now also includes secondary (2D) structure information for >13 million sequences, making RNAcentral the world's largest RNA 2D structure database. The 2D diagrams are displayed using R2DT, a new 2D structure visualization method that uses consistent, reproducible and recognizable layouts for related RNAs. The sequence similarity search has been updated with a faster interface featuring facets for filtering search results by RNA type, organism, source database or any keyword. This sequence search tool is available as a reusable web component, and has been integrated into several RNAcentral member databases, including Rfam, miRBase and snoDB. To allow for a more fine-grained assignment of RNA types and subtypes, all RNAcentral sequences have been annotated with Sequence Ontology terms. The RNAcentral database continues to grow and provide a central data resource for the RNA community. RNAcentral is freely available at https://rnacentral.org.

Viruses are integral components of all ecosystems and microbiomes on Earth. Through pervasive infections of their cellular hosts, viruses can reshape microbial community structure and drive global nutrient cycling. Over the past decade, viral sequences identified from genomes and metagenomes have provided an unprecedented view of viral genome diversity in nature. Since 2016, the IMG/VR database has provided access to the largest collection of viral sequences obtained from (meta)genomes. Here, we present the third version of IMG/VR, composed of 18 373 cultivated and 2 314 329 uncultivated viral genomes (UViGs), nearly tripling the total number of sequences compared to the previous version. These clustered into 935 362 viral Operational Taxonomic Units (vOTUs), including 188 930 with two or more members. UViGs in IMG/VR are now reported as single viral contigs, integrated proviruses or genome bins, and are annotated with a new standardized pipeline including genome quality estimation using CheckV, taxonomic classification reflecting the latest ICTV update, and expanded host taxonomy prediction. The new IMG/VR interface enables users to efficiently browse, search, and select UViGs based on genome features and/or sequence similarity. IMG/VR v3 is available at https://img.jgi.doe.gov/vr, and the underlying data are available to download at https://genome.jgi.doe.gov/portal/IMG_VR.

RNAcentral is a comprehensive database of non-coding RNA (ncRNA) sequences that provides a single access point to 44 RNA resources and >18 million ncRNA sequences from a wide range of organisms and RNA types. RNAcentral now also includes secondary (2D) structure information for >13 million sequences, making RNAcentral the world's largest RNA 2D structure database. The 2D diagrams are displayed using R2DT, a new 2D structure visualization method that uses consistent, reproducible and recognizable layouts for related RNAs. The sequence similarity search has been updated with a faster interface featuring facets for filtering search results by RNA type, organism, source database or any keyword. This sequence search tool is available as a reusable web component, and has been integrated into several RNAcentral member databases, including Rfam, miRBase and snoDB. To allow for a more fine-grained assignment of RNA types and subtypes, all RNAcentral sequences have been annotated with Sequence Ontology terms. The RNAcentral database continues to grow and provide a central data resource for the RNA community. RNAcentral is freely available at https://rnacentral.org.

DNA methyltransferases interact with their CpG target sites in the context of variable flanking sequences. We investigated DNA methylation by the human DNMT3B catalytic domain using substrate pools containing CpX target sites in randomized flanking context and identified combined effects of CpG recognition and flanking sequence interaction together with complex contact networks involved in balancing the interaction with different flanking sites. DNA methylation rates were more affected by flanking sequences at non-CpG than at CpG sites. We show that T775 has an essential dynamic role in the catalytic mechanism of DNMT3B. Moreover, we identify six amino acid residues in the DNA-binding interface of DNMT3B (N652, N656, N658, K777, N779, and R823), which are involved in the equalization of methylation rates of CpG sites in favored and disfavored sequence contexts by forming compensatory interactions to the flanking residues including a CpG specific contact to an A at the +1 flanking site. Non-CpG flanking preferences of DNMT3B are highly correlated with non-CpG methylation patterns in human cells. Comparison of the flanking sequence preferences of human and mouse DNMT3B revealed subtle differences suggesting a co-evolution of flanking sequence preferences and cellular DNMT targets.

Recognition of highly degenerate mammalian splice sites by the core spliceosomal machinery is regulated by several protein factors that predominantly bind exonic splicing motifs. These are postulated to be single-stranded in order to be functional, yet knowledge of secondary structural features that regulate the exposure of exonic splicing motifs across the transcriptome is not currently available. Using transcriptome-wide RNA structural information we show that retained introns in mouse are commonly flanked by a short (≲70 nucleotide), highly base-paired segment upstream and a predominantly single-stranded exonic segment downstream. Splicing assays with select pre-mRNA substrates demonstrate that loops immediately upstream of the introns contain pre-mRNA-specific splicing enhancers, the substitution or hybridization of which impedes splicing. Additionally, the exonic segments flanking the retained introns appeared to be more enriched in a previously identified set of hexameric exonic splicing enhancer (ESE) sequences compared to their spliced counterparts, suggesting that base-pairing in the exonic segments upstream of retained introns could be a means for occlusion of ESEs. The upstream exonic loops of the test substrate promoted recruitment of splicing factors and consequent pre-mRNA structural remodeling, leading up to assembly of the early spliceosome. These results suggest that disruption of exonic stem-loop structures immediately upstream (but not downstream) of the introns regulate alternative splicing events, likely through modulating accessibility of splicing factors.

Human ADAR3 is a catalytically inactive member of the Adenosine Deaminase Acting on RNA (ADAR) protein family, whose active members catalyze A-to-I RNA editing in metazoans. Until now, the reasons for the catalytic incapability of ADAR3 has not been defined and its biological function rarely explored. Yet, its exclusive expression in the brain and involvement in learning and memory suggest a central role in the nervous system. Here we describe the engineering of a catalytically active ADAR3 enzyme using a combination of computational design and functional screening. Five mutations (A389V, V485I, E527Q, Q549R and Q733D) engender RNA deaminase in human ADAR3. By way of its catalytic activity, the ADAR3 pentamutant was used to identify potential binding targets for wild type ADAR3 in a human glioblastoma cell line. Novel ADAR3 binding sites discovered in this manner include the 3'-UTRs of the mRNAs encoding early growth response 1 (EGR1) and dual specificity phosphatase 1 (DUSP1); both known to be activity-dependent immediate early genes that respond to stimuli in the brain. Further studies reveal that the wild type ADAR3 protein can regulate transcript levels for DUSP1 and EGR1, suggesting a novel role ADAR3 may play in brain function.

Non-homologous end-joining (NHEJ) plays an important role in double-strand break (DSB) repair of DNA. Recent studies have shown that the error patterns of NHEJ are strongly biased by sequence context, but these studies were based on relatively few templates. To investigate this more thoroughly, we systematically profiled ∼1.16 million independent mutational events resulting from CRISPR/Cas9-mediated cleavage and NHEJ-mediated DSB repair of 6872 synthetic target sequences, introduced into a human cell line via lentiviral infection. We find that: (i) insertions are dominated by 1 bp events templated by sequence immediately upstream of the cleavage site, (ii) deletions are predominantly associated with microhomology and (iii) targets exhibit variable but reproducible diversity with respect to the number and relative frequency of the mutational outcomes to which they give rise. From these data, we trained a model that uses local sequence context to predict the distribution of mutational outcomes. Exploiting the bias of NHEJ outcomes towards microhomology mediated events, we demonstrate the programming of deletion patterns by introducing microhomology to specific locations in the vicinity of the DSB site. We anticipate that our results will inform investigations of DSB repair mechanisms as well as the design of CRISPR/Cas9 experiments for diverse applications including genome-wide screens, gene therapy, lineage tracing and molecular recording.

Localized arrays of proteins cooperatively assemble onto chromosomes to control DNA activity in many contexts. Binding cooperativity is often mediated by specific protein-protein interactions, but cooperativity through DNA structure is becoming increasingly recognized as an additional mechanism. During the site-specific DNA recombination reaction that excises phage λ from the chromosome, the bacterial DNA architectural protein Fis recruits multiple λ-encoded Xis proteins to the attR recombination site. Here, we report X-ray crystal structures of DNA complexes containing Fis + Xis, which show little, if any, contacts between the two proteins. Comparisons with structures of DNA complexes containing only Fis or Xis, together with mutant protein and DNA binding studies, support a mechanism for cooperative protein binding solely by DNA allostery. Fis binding both molds the minor groove to potentiate insertion of the Xis β-hairpin wing motif and bends the DNA to facilitate Xis-DNA contacts within the major groove. The Fis-structured minor groove shape that is optimized for Xis binding requires a precisely positioned pyrimidine-purine base-pair step, whose location has been shown to modulate minor groove widths in Fis-bound complexes to different DNA targets.

The availability of genome-wide epigenomic datasets enables in-depth studies of epigenetic modifications and their relationships with chromatin structures and gene expression. Various alignment tools have been developed to align nucleotide or protein sequences in order to identify structurally similar regions. However, there are currently no alignment methods specifically designed for comparing multi-track epigenomic signals and detecting common patterns that may explain functional or evolutionary similarities. We propose a new local alignment algorithm, EpiAlign, designed to compare chromatin state sequences learned from multi-track epigenomic signals and to identify locally aligned chromatin regions. EpiAlign is a dynamic programming algorithm that novelly incorporates varying lengths and frequencies of chromatin states. We demonstrate the efficacy of EpiAlign through extensive simulations and studies on the real data from the NIH Roadmap Epigenomics project. EpiAlign is able to extract recurrent chromatin state patterns along a single epigenome, and many of these patterns carry cell-type-specific characteristics. EpiAlign can also detect common chromatin state patterns across multiple epigenomes, and it will serve as a useful tool to group and distinguish epigenomic samples based on genome-wide or local chromatin state patterns.

Transfer RNAs (tRNAs) are ubiquitous across the tree of life. Although tRNA structure is highly conserved, there is still significant variation in sequence features between clades, isotypes and even isodecoders. This variation not only impacts translation, but as shown by a variety of recent studies, nontranslation-associated functions are also sensitive to small changes in tRNA sequence. Despite the rapidly growing number of sequenced genomes, there is a lack of tools for both small- and large-scale comparative genomics analysis of tRNA sequence features. Here, we have integrated over 150 000 tRNAs spanning all domains of life into tRNAviz, a web application for exploring and visualizing tRNA sequence features. tRNAviz implements a framework for determining consensus sequence features and can generate sequence feature distributions by isotypes, clades and anticodons, among other tRNA properties such as score. All visualizations are interactive and exportable. The web server is publicly available at http://trna.ucsc.edu/tRNAviz/.

Deletions in the 16.6 kb mitochondrial genome have been implicated in numerous disorders that often display muscular and/or neurological symptoms due to the high-energy demands of these tissues. We describe a catalogue of 4489 putative mitochondrial DNA (mtDNA) deletions, including their frequency and relative read rate, using a combinatorial approach of mitochondria-targeted PCR, next-generation sequencing, bioinformatics, post-hoc filtering, annotation, and validation steps. Our bioinformatics pipeline uses MapSplice, an RNA-seq splice junction detection algorithm, to detect and quantify mtDNA deletion breakpoints rather than mRNA splices. Analyses of 93 samples from postmortem brain and blood found (i) the 4977 bp 'common deletion' was neither the most frequent deletion nor the most abundant; (ii) brain contained significantly more deletions than blood; (iii) many high frequency deletions were previously reported in MitoBreak, suggesting they are present at low levels in metabolically active tissues and are not exclusive to individuals with diagnosed mitochondrial pathologies; (iv) many individual deletions (and cumulative metrics) had significant and positive correlations with age and (v) the highest deletion burdens were observed in major depressive disorder brain, at levels greater than Kearns-Sayre Syndrome muscle. Collectively, these data suggest the Splice-Break pipeline can detect and quantify mtDNA deletions at a high level of resolution.

The piRNA pathway is an adaptive mechanism that maintains genome stability by repression of selfish genomic elements. In the male germline of Drosophila melanogaster repression of Stellate genes by piRNAs generated from Supressor of Stellate (Su(Ste)) locus is required for male fertility, but both Su(Ste) piRNAs and their targets are absent in other Drosophila species. We found that D. melanogaster genome contains multiple X-linked non-coding genomic repeats that have sequence similarity to the protein-coding host gene vasa. In the male germline, these vasa-related AT-chX repeats produce abundant piRNAs that are antisense to vasa; however, vasa mRNA escapes silencing due to imperfect complementarity to AT-chX piRNAs. Unexpectedly, we discovered AT-chX piRNAs target vasa of Drosophila mauritiana in the testes of interspecies hybrids. In the majority of hybrid flies, the testes were strongly reduced in size and germline content. A minority of hybrids maintained wild-type array of premeiotic germ cells in the testes, but in them harmful Stellate genes were derepressed due to the absence of Su(Ste) piRNAs, and meiotic failures were observed. Thus, the piRNA pathway contributes to reproductive isolation between D. melanogaster and closely related species, causing hybrid male sterility via misregulation of two different host protein factors.

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons' inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon-intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. 'intron definition') or across short exons (i.e. 'exon definition'), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

Telomerase adds telomeric repeats to chromosome ends by processive copying of a template within the telomerase RNA bound to telomerase reverse transcriptase. Telomerase RNAs have single-stranded regions that separate the template from a 5' stem and 3' pseudoknot, and mammals gained additional stem P2a.1 separating the template from the pseudoknot. Using human telomerase, we show that the length of template 3'-flanking single-stranded RNA is a determinant of repeat addition processivity whereas template 5'-flanking single-stranded RNA and P2a.1 are critical for activity but not processivity. In comparison, requirements for the template sequence itself are confounding: different substitutions of the same position have strikingly different consequences, from improved processivity and activity to complete inactivation. We discovered that some altered-template sequences stabilize an alternative RNA conformation that precludes the pseudoknot by base-pairing of one pseudoknot strand to the template 3' end. Using mutations to reduce over-stability of the alternative conformation, we restore high activity and processivity to otherwise inactive altered-template telomerase ribonucleoproteins. In cells, over-stabilization or destabilization of the alternative state severely inhibited biogenesis of active telomerase. Our findings delineate roles for human telomerase RNA template-flanking regions, establish a biologically relevant pseudoknot-alternative RNA conformation, and expand the repertoire of human telomerase repeat synthesis.

CRISPOR.org is a web tool for genome editing experiments with the CRISPR-Cas9 system. It finds guide RNAs in an input sequence and ranks them according to different scores that evaluate potential off-targets in the genome of interest and predict on-target activity. The list of genomes is continuously expanded, with more 150 genomes added in the last two years. CRISPOR tries to provide a comprehensive solution from selection, cloning and expression of guide RNA as well as providing primers needed for testing guide activity and potential off-targets. Recent developments include batch design for genome-wide CRISPR and saturation screens, creating custom oligonucleotides for guide cloning and the design of next generation sequencing primers to test for off-target mutations. CRISPOR is available from http://crispor.org, including the full source code of the website and a stand-alone, command-line version.

As molecular and cellular therapies advance in the clinic, the role of genetic regulation is becoming increasingly important for controlling therapeutic potency and safety. The emerging field of mammalian synthetic biology provides promising tools for the construction of regulatory platforms that can intervene with endogenous pathways and control cell behavior. Recent work has highlighted the development of synthetic biological systems that integrate sensing of molecular signals to regulated therapeutic function in various disease settings. However, the toxicity and limited dosing of currently available molecular inducers have largely inhibited translation to clinical settings. In this work, we developed synthetic microRNA-based genetic systems that are controlled by the pharmaceutical drug leucovorin, which is readily available and safe for prolonged administration in clinical settings. We designed microRNA switches to target endogenous cytokine receptor subunits (IL-2Rβ and γc) that mediate various signaling pathways in T cells. We demonstrate the function of these control systems by effectively regulating T cell proliferation with the drug input. Each control system produced unique functional responses, and combinatorial targeting of multiple receptor subunits exhibited greater repression of cell growth. This work highlights the potential use of drug-responsive genetic control systems to improve the management and safety of cellular therapeutics.

Gene therapy by engineering patient's own blood cells to confer HIV resistance can potentially lead to a functional cure for AIDS. Toward this goal, we have previously developed an anti-HIV lentivirus vector that deploys a combination of shRNA, ribozyme and RNA decoy. To further improve this therapeutic vector against viral escape, we sought an additional reagent to target HIV integrase. Here, we report the development of a new strategy for selection and expression of aptamer for gene therapy. We developed a SELEX protocol (multi-tag SELEX) for selecting RNA aptamers against proteins with low solubility or stability, such as integrase. More importantly, we expressed these aptamers in vivo by incorporating them in the terminal loop of shRNAs. This novel strategy allowed efficient expression of the shRNA-aptamer fusions that targeted RNAs and proteins simultaneously. Expressed shRNA-aptamer fusions targeting HIV integrase or reverse transcriptase inhibited HIV replication in cell cultures. Viral inhibition was further enhanced by combining an anti-integrase aptamer with an anti-HIV Tat-Rev shRNA. This construct exhibited efficacy comparable to that of integrase inhibitor Raltegravir. Our strategy for the selection and expression of RNA aptamers can potentially extend to other gene therapy applications.

Diversity-generating retroelements (DGRs) are novel genetic elements that use reverse transcription to generate vast numbers of sequence variants in specific target genes. Here, we present a detailed comparative bioinformatic analysis that depicts the landscape of DGR sequences in nature as represented by data in GenBank. Over 350 unique DGRs are identified, which together form a curated reference set of putatively functional DGRs. We classify target genes, variable repeats and DGR cassette architectures, and identify two new accessory genes. The great variability of target genes implies roles of DGRs in many undiscovered biological processes. There is much evidence for horizontal transfers of DGRs, and we identify lineages of DGRs that appear to have specialized properties. Because GenBank contains data from only 10% of described species, the compilation may not be wholly representative of DGRs present in nature. Indeed, many DGR subtypes are present only once in the set and DGRs of the candidate phylum radiation bacteria, and Diapherotrites, Parvarchaeota, Aenigmarchaeota, Nanoarchaeota, Nanohaloarchaea archaea, are exceptionally diverse in sequence, with little information available about functions of their target genes. Nonetheless, this study provides a detailed framework for classifying and studying DGRs as they are uncovered and studied in the future.

Translation rate per mRNA molecule correlates positively with mRNA abundance. As a result, protein levels do not scale linearly with mRNA levels, but instead scale with the abundance of mRNA raised to the power of an 'amplification exponent'. Here we show that to quantitate translational control, the translation rate must be decomposed into two components. One, TRmD, depends on the mRNA level and defines the amplification exponent. The other, TRmIND, is independent of mRNA amount and impacts the correlation coefficient between protein and mRNA levels. We show that in Saccharomyces cerevisiae TRmD represents ∼20% of the variance in translation and directs an amplification exponent of 1.20 with a 95% confidence interval [1.14, 1.26]. TRmIND constitutes the remaining ∼80% of the variance in translation and explains ∼5% of the variance in protein expression. We also find that TRmD and TRmIND are preferentially determined by different mRNA sequence features: TRmIND by the length of the open reading frame and TRmD both by a ∼60 nucleotide element that spans the initiating AUG and by codon and amino acid frequency. Our work provides more appropriate estimates of translational control and implies that TRmIND is under different evolutionary selective pressures than TRmD.

i-Motifs (iMs) are non-canonical, four-stranded secondary structures formed by stacking of hemi-protonated CH+·C base pairs in cytosine-rich DNA sequences, predominantly at pH < 7. The presence of iM structures in cells was a matter of debate until the recent development of iM-specific antibody, iMab, which was instrumental for several studies that suggested the existence of iMs in live cells and their putative biological roles. We assessed the interaction of iMab with cytosine-rich oligonucleotides by biolayer interferometry (BLI), pull-down assay and bulk-FRET experiments. Our results suggest that binding of iMab to DNA oligonucleotides is governed by the presence of runs of at least two consecutive cytosines and is generally increased in acidic conditions, irrespectively of the capacity of the sequence to adopt, or not, an iM structure. Moreover, the results of the bulk-FRET assay indicate that interaction with iMab results in unfolding of iM structures even in acidic conditions, similarly to what has been observed with hnRNP K, well-studied single-stranded DNA binding protein. Taken together, our results strongly suggest that iMab actually binds to blocks of 2-3 cytosines in single-stranded DNA, and call for more careful interpretation of results obtained with this antibody., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Vitamin B12 is an essential cofactor in all domains of life and B12-sensing riboswitches are some of the most widely distributed riboswitches. Mycobacterium tuberculosis, the causative agent of tuberculosis, harbours two B12-sensing riboswitches. One controls expression of metE, encoding a B12-independent methionine synthase, the other controls expression of ppe2 of uncertain function. Here, we analysed ligand sensing, secondary structure and gene expression control of the metE and ppe2 riboswitches. Our results provide the first evidence of B12 binding by these riboswitches and show that they exhibit different preferences for individual isoforms of B12, use distinct regulatory and structural elements and act as translational OFF switches. Based on our results, we propose that the ppe2 switch represents a new variant of Class IIb B12-sensing riboswitches. Moreover, we have identified short translated open reading frames (uORFs) upstream of metE and ppe2, which modulate the expression of their downstream genes. Translation of the metE uORF suppresses MetE expression, while translation of the ppe2 uORF is essential for PPE2 expression. Our findings reveal an unexpected regulatory interplay between B12-sensing riboswitches and the translational machinery, highlighting a new level of cis-regulatory complexity in M. tuberculosis. Attention to such mechanisms will be critical in designing next-level intervention strategies., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

We present CGeNArate, a new model for molecular dynamics simulations of very long segments of B-DNA in the context of biotechnological or chromatin studies. The developed method uses a coarse-grained Hamiltonian with trajectories that are back-mapped to the atomistic resolution level with extreme accuracy by means of Machine Learning Approaches. The method is sequence-dependent and reproduces very well not only local, but also global physical properties of DNA. The efficiency of the method allows us to recover with a reduced computational effort high-quality atomic-resolution ensembles of segments containing many kilobases of DNA, entering into the gene range or even the entire DNA of certain cellular organelles., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The computational design of synthetic DNA sequences with designer in vivo properties is gaining traction in the field of synthetic genomics. We propose here a computational method which combines a kinetic Monte Carlo framework with a deep mutational screening based on deep learning predictions. We apply our method to build regular nucleosome arrays with tailored nucleosomal repeat lengths (NRL) in yeast. Our design was validated in vivo by successfully engineering and integrating thousands of kilobases long tandem arrays of computationally optimized sequences which could accommodate NRLs much larger than the yeast natural NRL (namely 197 and 237 bp, compared to the natural NRL of ∼165 bp). RNA-seq results show that transcription of the arrays can occur but is not driven by the NRL. The computational method proposed here delineates the key sequence rules for nucleosome positioning in yeast and should be easily applicable to other sequence properties and other genomes., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Generative probabilistic models emerge as a new paradigm in data-driven, evolution-informed design of biomolecular sequences. This paper introduces a novel approach, called Edge Activation Direct Coupling Analysis (eaDCA), tailored to the characteristics of RNA sequences, with a strong emphasis on simplicity, efficiency, and interpretability. eaDCA explicitly constructs sparse coevolutionary models for RNA families, achieving performance levels comparable to more complex methods while utilizing a significantly lower number of parameters. Our approach demonstrates efficiency in generating artificial RNA sequences that closely resemble their natural counterparts in both statistical analyses and SHAPE-MaP experiments, and in predicting the effect of mutations. Notably, eaDCA provides a unique feature: estimating the number of potential functional sequences within a given RNA family. For example, in the case of cyclic di-AMP riboswitches (RF00379), our analysis suggests the existence of approximately 1039 functional nucleotide sequences. While huge compared to the known <4000 natural sequences, this number represents only a tiny fraction of the vast pool of nearly 1082 possible nucleotide sequences of the same length (136 nucleotides). These results underscore the promise of sparse and interpretable generative models, such as eaDCA, in enhancing our understanding of the expansive RNA sequence space., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Tandem donor splice sites (5'ss) are unique regions with at least two GU dinucleotides serving as splicing cleavage sites. The Δ3 tandem 5'ss are a specific subclass of 5'ss separated by 3 nucleotides which can affect protein function by inserting/deleting a single amino acid. One 5'ss is typically preferred, yet factors governing particular 5'ss choice are not fully understood. A highly conserved exon 21 of the STAT3 gene was chosen as a model to study Δ3 tandem 5'ss splicing mechanisms. Based on multiple lines of experimental evidence, endogenous U1 snRNA most likely binds only to the upstream 5'ss. However, the downstream 5'ss is used preferentially, and the splice site choice is not dependent on the exact U1 snRNA binding position. Downstream 5'ss usage was sensitive to exact nucleotide composition and dependent on the presence of downstream regulatory region. The downstream 5'ss usage could be best explained by two novel interactions with endogenous U6 snRNA. U6 snRNA enables the downstream 5'ss usage in STAT3 exon 21 by two mechanisms: (i) binding in a novel non-canonical register and (ii) establishing extended Watson-Crick base pairing with the downstream regulatory region. This study suggests that U6:5'ss interaction is more flexible than previously thought., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Characterizing the binding behaviors of RNA-binding proteins (RBPs) is important for understanding their functional roles in gene expression regulation. However, current high-throughput experimental methods for identifying RBP targets, such as CLIP-seq and RNAcompete, usually suffer from the false negative issue. Here, we develop a deep boosting based machine learning approach, called DeBooster, to accurately model the binding sequence preferences and identify the corresponding binding targets of RBPs from CLIP-seq data. Comprehensive validation tests have shown that DeBooster can outperform other state-of-the-art approaches in RBP target prediction. In addition, we have demonstrated that DeBooster may provide new insights into understanding the regulatory functions of RBPs, including the binding effects of the RNA helicase MOV10 on mRNA degradation, the potentially different ADAR1 binding behaviors related to its editing activity, as well as the antagonizing effect of RBP binding on miRNA repression. Moreover, DeBooster may provide an effective index to investigate the effect of pathogenic mutations in RBP binding sites, especially those related to splicing events. We expect that DeBooster will be widely applied to analyze large-scale CLIP-seq experimental data and can provide a practically useful tool for novel biological discoveries in understanding the regulatory mechanisms of RBPs. The source code of DeBooster can be downloaded from http://github.com/dongfanghong/deepboost.

Mutator-like transposable elements (MULEs) are widespread across fungal, plant and animal species. Despite their abundance and importance as genetic tools in plants, the transposition mechanism of the MULE superfamily was previously unknown. Discovery of the Muta1 element from Aedes aegypti and its successful transposition in yeast facilitated the characterization of key steps in Muta1 transposition. Here we show that purified transposase binds specifically to the Muta1 ends and catalyzes excision through double strand breaks (DSB) and the joining of newly excised transposon ends with target DNA. In the process, the DSB forms hairpin intermediates on the flanking DNA side. Analysis of transposase proteins containing site-directed mutations revealed the importance of the conserved DDE motif and a W residue. The transposition pathway resembles that of the V(D)J recombination reaction and the mechanism of hAT and Transib transposases including the importance of the conserved W residue in both MULEs and hATs. In addition, yeast transposition and in vitro assays demonstrated that the terminal motif and subterminal repeats of the Muta1 terminal inverted repeat also influence Muta1 transposition. Collectively, our data provides new insights to understand the evolutionary relationships between MULE, hAT and Transib elements and the V(D)J recombinase.

Realizing the full potential of genome editing requires the development of efficient and broadly applicable methods for delivering programmable nucleases and donor templates for homology-directed repair (HDR). The RNA-guided Cas9 endonuclease can be introduced into cells as a purified protein in complex with a single guide RNA (sgRNA). Such ribonucleoproteins (RNPs) can facilitate the high-fidelity introduction of single-base substitutions via HDR following co-delivery with a single-stranded DNA oligonucleotide. However, combining RNPs with transgene-containing donor templates for targeted gene addition has proven challenging, which in turn has limited the capabilities of the RNP-mediated genome editing toolbox. Here, we demonstrate that combining RNP delivery with naturally recombinogenic adeno-associated virus (AAV) donor vectors enables site-specific gene insertion by homology-directed genome editing. Compared to conventional plasmid-based expression vectors and donor templates, we show that combining RNP and AAV donor delivery increases the efficiency of gene addition by up to 12-fold, enabling the creation of lineage reporters that can be used to track the conversion of striatal neurons from human fibroblasts in real time. These results thus illustrate the potential for unifying nuclease protein delivery with AAV donor vectors for homology-directed genome editing.

The genomic loci bound by the glucocorticoid receptor (GR), a hormone-activated transcription factor, show little overlap between cell types. To study the role of chromatin and sequence in specifying where GR binds, we used Bayesian modeling within the universe of accessible chromatin. Taken together, our results uncovered that although GR preferentially binds accessible chromatin, its binding is biased against accessible chromatin located at promoter regions. This bias can only be explained partially by the presence of fewer GR recognition sequences, arguing for the existence of additional mechanisms that interfere with GR binding at promoters. Therefore, we tested the role of H3K9ac, the chromatin feature with the strongest negative association with GR binding, but found that this correlation does not reflect a causative link. Finally, we find a higher percentage of promoter-proximal GR binding for genes regulated by GR across cell types than for cell type-specific target genes. Given that GR almost exclusively binds accessible chromatin, we propose that cell type-specific regulation by GR preferentially occurs via distal enhancers, whose chromatin accessibility is typically cell type-specific, whereas ubiquitous target gene regulation is more likely to result from binding to promoter regions, which are often accessible regardless of cell type examined.

Telomerase maintains telomere length at the ends of linear chromosomes using an integral telomerase RNA (TER) and telomerase reverse transcriptase (TERT). An essential part of TER is the template/pseudoknot domain (t/PK) which includes the template, for adding telomeric repeats, template boundary element (TBE), and pseudoknot, enclosed in a circle by stem 1. The Tetrahymena telomerase holoenzyme catalytic core (p65-TER-TERT) was recently modeled in our 9 Å resolution cryo-electron microscopy map by fitting protein and TER domains, including a solution NMR structure of the Tetrahymena pseudoknot. Here, we describe in detail the structure and folding of the isolated pseudoknot, which forms a compact structure with major groove U•A-U and novel C•G-A+ base triples. Base substitutions that disrupt the base triples reduce telomerase activity in vitro NMR studies also reveal that the pseudoknot does not form in the context of full-length TER in the absence of TERT, due to formation of a competing structure that sequesters pseudoknot residues. The residues around the TBE remain unpaired, potentially providing access by TERT to this high affinity binding site during an early step in TERT-TER assembly. A model for the assembly pathway of the catalytic core is proposed.

Adenosine deamination is one of the most prevalent post-transcriptional modifications in mRNA. In humans, ADAR1 and ADAR2 catalyze this modification and their malfunction correlates with disease. Recently our laboratory reported crystal structures of the human ADAR2 deaminase domain bound to duplex RNA revealing a protein loop that binds the RNA on the 5' side of the modification site. This 5' binding loop appears to be one contributor to substrate specificity differences between ADAR family members. In this study, we endeavored to reveal detailed structure-activity relationships in this loop to advance our understanding of RNA recognition by ADAR2. To achieve this goal, we established a high-throughput mutagenesis approach which allows rapid screening of ADAR variants in single yeast cells and provides quantitative evaluation for enzymatic activity. Using this approach, we determined the importance of specific amino acids at 19 different positions in the ADAR2 5' binding loop and revealed six residues that provide essential structural elements supporting the fold of the loop and key RNA-binding functional groups. This work provided new insight into RNA recognition by ADAR2 and established a new tool for defining structure-function relationships in ADAR reactions.

The 7SK small nuclear ribonucleoprotein (snRNP) sequesters and inactivates the positive transcription elongation factor b (P-TEFb), an essential eukaryotic mRNA transcription factor. The human La-related protein group 7 (hLARP7) is a constitutive component of the 7SK snRNP and localizes to the 3' terminus of the 7SK long noncoding RNA. hLARP7, and in particular its C-terminal domain (CTD), is essential for 7SK RNA stability and assembly with P-TEFb. The hLARP7 N-terminal La module binds and protects the 3' end from degradation, but the structural and functional role of its CTD is unclear. We report the solution NMR structure of the hLARP7 CTD and show that this domain contains an xRRM, a class of atypical RRM first identified in the Tetrahymena thermophila telomerase LARP7 protein p65. The xRRM binds the 3' end of 7SK RNA at the top of stem-loop 4 (SL4) and interacts with both unpaired and base-paired nucleotides. This study confirms that the xRRM is general to the LARP7 family of proteins and defines the binding site for hLARP7 on the 7SK RNA, providing insight into function.

High-throughput genome sequencing continues to grow the need for rapid, accurate genome annotation and tRNA genes constitute the largest family of essential, ever-present non-coding RNA genes. Newly developed tRNAscan-SE 2.0 has advanced the state-of-the-art methodology in tRNA gene detection and functional prediction, captured by rich new content of the companion Genomic tRNA Database. Previously, web-server tRNA detection was isolated from knowledge of existing tRNAs and their annotation. In this update of the tRNAscan-SE On-line resource, we tie together improvements in tRNA classification with greatly enhanced biological context via dynamically generated links between web server search results, the most relevant genes in the GtRNAdb and interactive, rich genome context provided by UCSC genome browsers. The tRNAscan-SE On-line web server can be accessed at http://trna.ucsc.edu/tRNAscan-SE/.

A set of new data types emerged from functional genomic assays, including ChIP-seq, DNase-seq, FAIRE-seq and others. The results are typically stored as genome-wide intensities (WIG/bigWig files) or functional genomic regions (peak/BED files). These data types present new challenges to big data science. Here, we present GeNemo, a web-based search engine for functional genomic data. GeNemo searches user-input data against online functional genomic datasets, including the entire collection of ENCODE and mouse ENCODE datasets. Unlike text-based search engines, GeNemo's searches are based on pattern matching of functional genomic regions. This distinguishes GeNemo from text or DNA sequence searches. The user can input any complete or partial functional genomic dataset, for example, a binding intensity file (bigWig) or a peak file. GeNemo reports any genomic regions, ranging from hundred bases to hundred thousand bases, from any of the online ENCODE datasets that share similar functional (binding, modification, accessibility) patterns. This is enabled by a Markov Chain Monte Carlo-based maximization process, executed on up to 24 parallel computing threads. By clicking on a search result, the user can visually compare her/his data with the found datasets and navigate the identified genomic regions. GeNemo is available at www.genemo.org.

Alternative splicing (AS) is a robust generator of mammalian transcriptome complexity. Splice site specification is controlled by interactions of cis-acting determinants on a transcript with specific RNA binding proteins. These interactions are frequently localized to the intronic U-rich polypyrimidine tracts (PPT) located 5' to the majority of splice acceptor junctions. αCPs (also referred to as polyC-binding proteins (PCBPs) and hnRNPEs) comprise a subset of KH-domain proteins with high affinity and specificity for C-rich polypyrimidine motifs. Here, we demonstrate that αCPs promote the splicing of a defined subset of cassette exons via binding to a C-rich subset of polypyrimidine tracts located 5' to the αCP-enhanced exonic segments. This enhancement of splice acceptor activity is linked to interactions of αCPs with the U2 snRNP complex and may be mediated by cooperative interactions with the canonical polypyrimidine tract binding protein, U2AF65. Analysis of αCP-targeted exons predicts a substantial impact on fundamental cell functions. These findings lead us to conclude that the αCPs play a direct and global role in modulating the splicing activity and inclusion of an array of cassette exons, thus driving a novel pathway of splice site regulation within the mammalian transcriptome.

GC skew is a measure of the strand asymmetry in the distribution of guanines and cytosines. GC skew favors R-loops, a type of three stranded nucleic acid structures that form upon annealing of an RNA strand to one strand of DNA, creating a persistent RNA:DNA hybrid. Previous studies show that GC skew is prevalent at thousands of human CpG island (CGI) promoters and transcription termination regions, which correspond to hotspots of R-loop formation. Here, we investigated the conservation of GC skew patterns in 60 sequenced chordates genomes. We report that GC skew is a conserved sequence characteristic of the CGI promoter class in vertebrates. Furthermore, we reveal that promoter GC skew peaks at the exon 1/ intron1 junction and that it is highly correlated with gene age and CGI promoter strength. Our data also show that GC skew is predictive of unmethylated CGI promoters in a range of vertebrate species and that it imparts significant DNA hypomethylation for promoters with intermediate CpG densities. Finally, we observed that terminal GC skew is conserved for a subset of vertebrate genes that tend to be located significantly closer to their downstream neighbors, consistent with a role for R-loop formation in transcription termination.

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The bendability of genomic DNA, which measures the DNA looping rate, is crucial for numerous biological processes of DNA. Recently, an advanced high-throughput technique known as 'loop-seq' has made it possible to measure the inherent cyclizability of DNA fragments. However, quantifying the bendability of large-scale DNA is costly, laborious, and time-consuming. To close the gap between rapidly evolving large language models and expanding genomic sequence information, and to elucidate the DNA bendability's impact on critical regulatory sequence motifs such as super-enhancers in the human genome, we introduce an innovative computational model, named MIXBend, to forecast the DNA bendability utilizing both nucleotide sequences and physicochemical properties. In MIXBend, a pre-trained language model DNABERT and convolutional neural network with attention mechanism are utilized to construct both sequence- and physicochemical-based extractors for the sophisticated refinement of DNA sequence representations. These bimodal DNA representations are then fed to a k-mer sequence-physicochemistry matching module to minimize the semantic gap between each modality. Lastly, a self-attention fusion layer is employed for the prediction of DNA bendability. In conclusion, the experimental results validate MIXBend's superior performance relative to other state-of-the-art methods. Additionally, MIXBend reveals both novel and known motifs from the yeast. Moreover, MIXBend discovers significant bendability fluctuations within super-enhancer regions and transcription factors binding sites in the human genome., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Custom oligonucleotides (oligos) are widely used reagents in biomedical research. Some common applications of oligos include polymerase chain reaction (PCR), sequencing, hybridization, microarray, and library construction. The reliability of oligos in such applications depends on their purity and specificity. Here, we report that commercially available oligos are frequently contaminated with nonspecific sequences (i.e. other unrelated oligonucleotides). Most of the oligos that we designed to amplify clustered regularly interspersed palindromic repeats (CRISPR) guide sequences contained nonspecific CRISPR guides. These contaminants were detected in research-grade oligos procured from eight commercial oligo-suppliers located in three different geographic regions of the world. Deep sequencing of some of the oligos revealed a variety of contaminants. Given the wide range of applications of oligos, the impact of oligo cross-contamination varies greatly depending on the field and the experimental method. Incorporating appropriate control experiments in research design can help ensure that the quality of oligo reagents meets the intended purpose. This can also minimize risk depending on the purposes for which the oligos are used., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

During their maturation, ribosomal RNAs (rRNAs) are decorated by hundreds of chemical modifications that participate in proper folding of rRNA secondary structures and therefore in ribosomal function. Along with pseudouridine, methylation of the 2'-hydroxyl ribose moiety (Nm) is the most abundant modification of rRNAs. The majority of Nm modifications in eukaryotes are placed by Fibrillarin, a conserved methyltransferase belonging to a ribonucleoprotein complex guided by C/D box small nucleolar RNAs (C/D box snoRNAs). These modifications impact interactions between rRNAs, tRNAs and mRNAs, and some are known to fine tune translation rates and efficiency. In this study, we built the first comprehensive map of Nm sites in Drosophila melanogaster rRNAs using two complementary approaches (RiboMethSeq and Nanopore direct RNA sequencing) and identified their corresponding C/D box snoRNAs by whole-transcriptome sequencing. We de novo identified 61 Nm sites, from which 55 are supported by both sequencing methods, we validated the expression of 106 C/D box snoRNAs and we predicted new or alternative rRNA Nm targets for 31 of them. Comparison of methylation level upon different stresses show only slight but specific variations, indicating that this modification is relatively stable in D. melanogaster. This study paves the way to investigate the impact of snoRNA-mediated 2'-O-methylation on translation and proteostasis in a whole organism., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Tandem-repetitive DNA (where two or more DNA bases are repeated numerous times) can adopt non-canonical secondary structures. Many of these structures are implicated in important biological processes. Human Satellite III (HSat3) is enriched for tandem repeats of the sequence ATGGA and is located in pericentromeric heterochromatin in many human chromosomes. Here, we investigate the secondary structure of the four-repeat HSat3 sequence 5'-ATGGA ATGGA ATGGA ATGGA-3' using X-ray crystallography, NMR, and biophysical methods. Circular dichroism spectroscopy, thermal stability, native PAGE, and analytical ultracentrifugation indicate that this sequence folds into a monomolecular hairpin with non-canonical base pairing and B-DNA characteristics at concentrations below 0.9 mM. NMR studies at 0.05-0.5 mM indicate that the hairpin is likely folded-over into a compact structure with high dynamics. Crystallographic studies at 2.5 mM reveal an antiparallel self-complementary duplex with the same base pairing as in the hairpin, extended into an infinite polymer. The non-canonical base pairing includes a G-G intercalation sandwiched by sheared A-G base pairs, leading to a cross-strand four guanine stack, so called guanine zipper. The guanine zippers are spaced throughout the structure by A-T/T-A base pairs. Our findings lend further insight into recurring structural motifs associated with the HSat3 and their potential biological functions., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

i-Motifs (iMs), are secondary structures formed in cytosine-rich DNA sequences and are involved in multiple functions in the genome. Although putative iM forming sequences are widely distributed in the human genome, the folding status and strength of putative iMs vary dramatically. Much previous research on iM has focused on assessing the iM folding properties using biophysical experiments. However, there are no dedicated computational tools for predicting the folding status and strength of iM structures. Here, we introduce a machine learning pipeline, iM-Seeker, to predict both folding status and structural stability of DNA iMs. The programme iM-Seeker incorporates a Balanced Random Forest classifier trained on genome-wide iMab antibody-based CUT&Tag sequencing data to predict the folding status and an Extreme Gradient Boosting regressor to estimate the folding strength according to both literature biophysical data and our in-house biophysical experiments. iM-Seeker predicts DNA iM folding status with a classification accuracy of 81% and estimates the folding strength with coefficient of determination (R2) of 0.642 on the test set. Model interpretation confirms that the nucleotide composition of the C-rich sequence significantly affects iM stability, with a positive correlation with sequences containing cytosine and thymine and a negative correlation with guanine and adenine., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Gene knockout strategies, RNAi and rescue experiments are all employed to study mammalian gene function. However, the disadvantages of these approaches include: loss of function adaptation, reduced viability and gene overexpression that rarely matches endogenous levels. Here, we developed an endogenous gene knockdown/rescue strategy that combines RNAi selectivity with a highly efficient CRISPR directed recombinant Adeno-Associated Virus (rAAV) mediated gene targeting approach to introduce allele-specific mutations plus an allele-selective siRNA Sensitive (siSN) site that allows for studying gene mutations while maintaining endogenous expression and regulation of the gene of interest. CRISPR/Cas9 plus rAAV targeted gene-replacement and introduction of allele-specific RNAi sensitivity mutations in the CDK2 and CDK1 genes resulted in a >85% site-specific recombination of Neo-resistant clones versus ∼8% for rAAV alone. RNAi knockdown of wild type (WT) Cdk2 with siWT in heterozygotic knockin cells resulted in the mutant Cdk2 phenotype cell cycle arrest, whereas allele specific knockdown of mutant CDK2 with siSN resulted in a wild type phenotype. Together, these observations demonstrate the ability of CRISPR plus rAAV to efficiently recombine a genomic locus and tag it with a selective siRNA sequence that allows for allele-selective phenotypic assays of the gene of interest while it remains expressed and regulated under endogenous control mechanisms.

RNA research and therapy relies primarily on synthetic RNAs. We employed recombinant RNA technology toward large-scale production of pre-miRNA agents in bacteria, but found the majority of target RNAs were not or negligibly expressed. We thus developed a novel strategy to achieve consistent high-yield biosynthesis of chimeric RNAs carrying various small RNAs (e.g. miRNAs, siRNAs and RNA aptamers), which was based upon an optimal noncoding RNA scaffold (OnRS) derived from tRNA fusion pre-miR-34a (tRNA/mir-34a). Multi-milligrams of chimeric RNAs (e.g. OnRS/miR-124, OnRS/GFP-siRNA, OnRS/Neg (scrambled RNA) and OnRS/MGA (malachite green aptamer)) were readily obtained from 1 l bacterial culture. Deep sequencing analyses revealed that mature miR-124 and target GFP-siRNA were selectively released from chimeric RNAs in human cells. Consequently, OnRS/miR-124 was active in suppressing miR-124 target gene expression and controlling cellular processes, and OnRS/GFP-siRNA was effective in knocking down GFP mRNA levels and fluorescent intensity in ES-2/GFP cells and GFP-transgenic mice. Furthermore, the OnRS/MGA sensor offered a specific strong fluorescence upon binding MG, which was utilized as label-free substrate to accurately determine serum RNase activities in pancreatic cancer patients. These results demonstrate that OnRS-based bioengineering is a common, robust and versatile strategy to assemble various types of small RNAs for broad applications.

Clustered regularly interspaced short palindromic repeat (CRISPR) RNA-guided nucleases have gathered considerable excitement as a tool for genome engineering. However, questions remain about the specificity of target site recognition. Cleavage specificity is typically evaluated by low throughput assays (T7 endonuclease I assay, target amplification followed by high-throughput sequencing), which are limited to a subset of potential off-target sites. Here, we used ChIP-seq to examine genome-wide CRISPR binding specificity at gRNA-specific and gRNA-independent sites for two guide RNAs. RNA-guided Cas9 binding was highly specific to the target site while off-target binding occurred at much lower intensities. Cas9-bound regions were highly enriched in NGG sites, a sequence required for target site recognition by Streptococcus pyogenes Cas9. To determine the relationship between Cas9 binding and endonuclease activity, we applied targeted sequence capture, which allowed us to survey 1200 genomic loci simultaneously including potential off-target sites identified by ChIP-seq and by computational prediction. A high frequency of indels was observed at both target sites and one off-target site, while no cleavage activity could be detected at other ChIP-bound regions. Our results confirm the high-specificity of CRISPR endonucleases and demonstrate that sequence capture can be used as a high-throughput genome-wide approach to identify off-target activity.

DNA-binding motifs that are recognized by transcription factors (TFs) have been well studied; however, challenges remain in determining the in vivo architecture of TF-DNA complexes on a genome-scale. Here, we determined the in vivo architecture of Escherichia coli arginine repressor (ArgR)-DNA complexes using high-throughput sequencing of exonuclease-treated chromatin-immunoprecipitated DNA (ChIP-exo). The ChIP-exo has a unique peak-pair pattern indicating 5' and 3' ends of ArgR-binding region. We identified 62 ArgR-binding loci, which were classified into three groups, comprising single, double and triple peak-pairs. Each peak-pair has a unique 93 base pair (bp)-long (±2 bp) ArgR-binding sequence containing two ARG boxes (39 bp) and residual sequences. Moreover, the three ArgR-binding modes defined by the position of the two ARG boxes indicate that DNA bends centered between the pair of ARG boxes facilitate the non-specific contacts between ArgR subunits and the residual sequences. Additionally, our approach may also reveal other fundamental structural features of TF-DNA interactions that have implications for studying genome-scale transcriptional regulatory networks.

Sequence-specific nucleases such as TALEN and the CRISPR/Cas9 system have so far been used to disrupt, correct or insert transgenes at precise locations in mammalian genomes. We demonstrate efficient 'knock-in' targeted replacement of multi-kilobase genes in human induced pluripotent stem cells (iPSC). Using a model system replacing endogenous human genes with their mouse counterpart, we performed a comprehensive study of targeting vector design parameters for homologous recombination. A 2.7 kilobase (kb) homozygous gene replacement was achieved in up to 11% of iPSC without selection. The optimal homology arm length was around 2 kb, with homology length being especially critical on the arm not adjacent to the cut site. Homologous sequence inside the cut sites was detrimental to targeting efficiency, consistent with a synthesis-dependent strand annealing (SDSA) mechanism. Using two nuclease sites, we observed a high degree of gene excisions and inversions, which sometimes occurred more frequently than indel mutations. While homozygous deletions of 86 kb were achieved with up to 8% frequency, deletion frequencies were not solely a function of nuclease activity and deletion size. Our results analyzing the optimal parameters for targeting vector design will inform future gene targeting efforts involving multi-kilobase gene segments, particularly in human iPSC.

MicroRNAs (miRNAs) typically downregulate protein expression from target mRNAs through limited base-pairing interactions between the 5' 'seed' region of the miRNA and the mRNA 3' untranslated region (3'UTR). In contrast to this established mode of action, the liver-specific human miR-122 binds at two sites within the hepatitis C viral (HCV) 5'UTR, leading to increased production of infectious virions. We show here that two copies of miR-122 interact with the HCV 5'UTR at partially overlapping positions near the 5' end of the viral transcript to form a stable ternary complex. Both miR-122 binding sites involve extensive base pairing outside of the seed sequence; yet, they have substantially different interaction affinities. Structural probing reveals changes in the architecture of the HCV 5'UTR that occur on interaction with miR-122. In contrast to previous reports, however, results using both the recombinant cytoplasmic exonuclease Xrn1 and liver cell extracts show that miR-122-mediated protection of the HCV RNA from degradation does not correlate with stimulation of viral propagation in vivo. Thus, the miR-122:HCV ternary complex likely functions at other steps critical to the viral life cycle.