Search

Type III CRISPR immune systems bind viral or plasmid RNA transcripts and activate Csm3/Cmr4 and Cas10 nucleases to uniquely cleave both invader RNA and DNA, respectively. Additionally, type III effector complexes generate cyclic oligoadenylate (cOA) signaling molecules to activate trans-acting, auxiliary Csm6/Csx1 ribonucleases, previously proposed to be non-specific in their in vivo RNA cleavage preference. Despite extensive in vitro studies, the nuclease requirements of type III systems in their native contexts remain poorly understood. Here we systematically investigated the in vivo roles for immunity of each of the three Streptococcus thermophilus (Sth) type III-A Cas nucleases and cOA signaling by challenging nuclease defective mutant strains with plasmid and phage infections. Our results reveal that RNA cleavage by Csm6 is both sufficient and essential for maintaining wild-type levels of immunity. Importantly, Csm6 RNase activity leads to immunity against even high levels of phage challenge without causing host cell dormancy or death. Transcriptomic analyses during phage infection indicated Csm6-mediated and crRNA-directed preferential cleavage of phage transcripts. Our findings highlight the critical role of Csm6 RNase activity in type III immunity and demonstrate specificity for invader RNA transcripts by Csm6 to ensure host cell survival upon phage infection., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

RNA interactome studies have revealed that hundreds of zinc-finger proteins (ZFPs) are candidate RNA-binding proteins (RBPs), yet their RNA substrates and functional significance remain largely uncharacterized. Here, we present a systematic multi-omics analysis of the DNA- and RNA-binding targets and regulatory roles of more than 100 ZFPs representing 37 zinc-finger families. We show that multiple ZFPs are previously unknown regulators of RNA splicing, alternative polyadenylation, stability, or translation. The examined ZFPs show widespread sequence-specific RNA binding and preferentially bind proximal to transcription start sites. Additionally, several ZFPs associate with their targets at both the DNA and RNA levels. We highlight ZNF277, a C2H2 ZFP that binds thousands of RNA targets and acts as a multi-functional RBP. We also show that ZNF473 is a DNA/RNA-associated protein that regulates the expression and splicing of cell cycle genes. Our results reveal diverse roles for ZFPs in transcriptional and post-transcriptional gene regulation., Competing Interests: Declaration of interests G.W.Y. is a scientific advisory board (SAB) member of Jumpcode Genomics and a co-founder, member of the board of directors, on the SAB, equity holder, and paid consultant for Eclipse BioInnovations. G.W.Y. is a Distinguished Visiting Professor at the National University of Singapore. G.W.Y.’s interests have been reviewed and approved by UCSD in accordance with its conflict-of-interest policies. B.R.G. is a co-founder and SAB member for RNAConnect., (Copyright © 2024 University of California San Diego. Published by Elsevier Inc. All rights reserved.)

The prevalence of obesity is increasing in older adults and contributes to age-related decline. Caloric restriction (CR) alleviates obesity phenotypes and delays the onset of age-related changes. However, how late in life organisms benefit from switching from a high-(H) to a low-calorie (L) diet is unclear. We transferred male flies from a H to a L (HL) diet or vice versa (LH) at different times during life. Both shifts immediately change fly rate of aging even when applied late in life. HL shift rapidly reduces fly mortality rate to briefly lower rate than in flies on a constant L diet, and extends lifespan. Transcriptomic analysis uncovers that flies aged on H diet have acquired increased stress response, which may have temporal advantage over flies aged on L diet and leads to rapid decrease in mortality rate after HL switch. Conversely, a LH shift increases mortality rate, which is temporarily higher than in flies aged on a H diet, and shortens lifespan. Unexpectedly, more abundant transcriptomic changes accompanied LH shift, including increase in ribosome biogenesis, stress response and growth. These changes reflect protection from sudden release of ROS, energy storage, and use of energy to growth, which all likely contribute to higher mortality rate. As the beneficial effects of CR on physiology and lifespan are conserved across many organisms, our study provides framework to study underlying mechanisms of CR interventions that counteract the detrimental effects of H diets and reduce rate of aging even when initiated later in life., Competing Interests: Competing interests statement:The authors declare no competing interest.

CRISPR-Cas systems provide heritable immunity against viruses and other mobile genetic elements by incorporating fragments of invader DNA into the host CRISPR array as spacers. Integration of new spacers is localized to the 5' end of the array, and in certain Gram-negative Bacteria this polarized localization is accomplished by the integration host factor. For most other Bacteria and Archaea, the mechanism for 5' end localization is unknown. Here we show that archaeal histones play a key role in directing integration of CRISPR spacers. In Pyrococcus furiosus, deletion of either histone A or B impairs integration. In vitro, purified histones are sufficient to direct integration to the 5' end of the CRISPR array. Archaeal histone tetramers and bacterial integration host factor induce similar U-turn bends in bound DNA. These findings indicate a co-evolution of CRISPR arrays with chromosomal DNA binding proteins and a widespread role for binding and bending of DNA to facilitate accurate spacer integration., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

Prokaryotes are under constant pressure from phage infection and thus have evolved multiple means of defense or evasion. While CRISPR-Cas constitutes a robust immune system and appears to be the predominant means of survival for Streptococcus thermophilus when facing lytic phage infection, other forms of phage resistance coexist in this species. Here, we show that S. thermophilus strains with deleted CRISPR-Cas loci can still give rise to phage-resistant clones following lytic phage challenge. Notably, non-CRISPR phage-resistant survivors had multiple mutations which would truncate or recode a membrane-anchored host protease, FtsH. Phage adsorption was dramatically reduced in FtsH mutants, implicating this protein in phage attachment. Phages were isolated which could bypass FtsH-based resistance through mutations predicted to alter tape measure protein translation. Together, these results identify key components in phage propagation that are subject to mutation in the molecular arms race between phage and host cell. IMPORTANCE Streptococcus thermophilus is an important organism for production of cultured dairy foods, but it is susceptible to lytic phages which can lead to failed products. Consequently, mechanisms for phage resistance are an active area of research. One such mechanism is CRISPR-Cas, and S. thermophilus is a model organism for the study of this form of adaptive immunity. Here, we expand on known mechanisms with our finding that spontaneous mutations in ftsH , a gene encoding a membrane-anchored protease, protected against phage infection by disrupting phage adsorption. In turn, mutations in phage tail protein genes allowed phages to overcome ftsH -based resistance. Our results identified components in phage propagation that are subject to mutation in the molecular arms race between phage and host., Competing Interests: The authors declare no conflict of interest.

RNA binding proteins (RBPs) are central regulators of gene expression implicated in all facets of RNA metabolism. As such, they play key roles in cellular physiology and disease etiology. Since different steps of post-transcriptional gene expression tend to occur in specific regions of the cell, including nuclear or cytoplasmic locations, defining the subcellular distribution properties of RBPs is an important step in assessing their potential functions. Here, we present the RBP Image Database, a resource that details the subcellular localization features of 301 RBPs in the human HepG2 and HeLa cell lines, based on the results of systematic immuno-fluorescence studies conducted using a highly validated collection of RBP antibodies and a panel of 12 markers for specific organelles and subcellular structures. The unique features of the RBP Image Database include: (i) hosting of comprehensive representative images for each RBP-marker pair, with ∼250,000 microscopy images; (ii) a manually curated controlled vocabulary of annotation terms detailing the localization features of each factor; and (iii) a user-friendly interface allowing the rapid querying of the data by target or annotation. The RBP Image Database is freely available at https://rnabiology.ircm.qc.ca/RBPImage/., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Pyrococcus furiosus is a hyperthermophilic archaeon with three effector CRISPR complexes (types I-A, I-B, and III-B) that each employ crRNAs derived from seven CRISPR arrays. Here, we investigate the CRISPR adaptation response to a newly discovered and self-transmissible plasmid, pT33.3. Transconjugant strains of Pyrococcus furiosus exhibited dramatically elevated levels of new spacer integration at CRISPR loci relative to the strain harboring a commonly employed, laboratory-constructed plasmid. High-throughput sequence analysis demonstrated that the vast majority of the newly acquired spacers were preferentially selected from DNA surrounding a particular region of the pT33.3 plasmid and exhibited a bi-directional pattern of strand bias that is a hallmark of primed adaptation by type I systems. We observed that one of the CRISPR arrays of our Pyrococcus furiosus laboratory strain encodes a spacer that closely matches the region of the conjugative plasmid that is targeted for adaptation. The hyper-adaptation phenotype was found to strictly depend both on the presence of this single matching spacer as well as the I-B effector complex, known to mediate primed adaptation. Our results indicate that Pyrococcus furiosus naturally encountered this conjugative plasmid or a related mobile genetic element in the past and responds to reinfection with robust primed adaptation., (© 2022. The Author(s), under exclusive licence to Springer Japan KK, part of Springer Nature.)

CRISPR-Cas systems are functionally diverse prokaryotic antiviral defense systems, which encompass six distinct types (I-VI) that each encode different effector Cas nucleases with distinct nucleic acid cleavage specificities. By harnessing the unique attributes of the various CRISPR-Cas systems, a range of innovative CRISPR-based DNA and RNA targeting tools and technologies have been developed. Here, we exploit the ability of type III-A CRISPR-Cas systems to carry out RNA-guided and sequence-specific target RNA cleavage for establishment of research tools for post-transcriptional control of gene expression. Type III-A systems from three bacterial species ( L. lactis , S. epidermidis , and S. thermophilus ) were each expressed on a single plasmid in E. coli , and the efficiency and specificity of gene knockdown was assessed by northern blot and transcriptomic analysis. We show that engineered type III-A modules can be programmed using tailored CRISPR RNAs to efficiently knock down gene expression of both coding and noncoding RNAs in vivo. Moreover, simultaneous degradation of multiple cellular mRNA transcripts can be directed by utilizing a CRISPR array expressing corresponding gene-targeting crRNAs. Our results demonstrate the utility of distinct type III-A modules to serve as specific and effective gene knockdown platforms in heterologous cells. This transcriptome engineering technology has the potential to be further refined and exploited for key applications including gene discovery and gene pathway analyses in additional prokaryotic and perhaps eukaryotic cells and organisms., (© 2022 Woodside et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Although reverse-transcriptase (RT) enzymes are critical reagents for research and biotechnology, their mechanical properties are not well understood. In particular, we know little about their relative speed and response to structural obstacles in the template. Commercial retroviral RTs stop at many positions along mixed sequence templates, resulting in truncated cDNA products that complicate downstream analysis. By contrast, group II intron-encoded RTs appear to copy long RNAs with high processivity and minimal stops. However, their speed, consistency and pausing behavior have not been explored. Here, we analyze RT velocity as the enzyme moves through heterogeneous sequences and structures that are embedded within a long noncoding RNA transcript. We observe that heterogeneities in the template are highly disruptive to primer extension by retroviral RTs. However, sequence composition and template structure have negligible effects on behavior of group II intron RTs, such as MarathonRT (MRT). Indeed, MRT copies long RNAs in a single pass, and displays synchronized primer extension at a constant speed of 25 nt/sec. In addition, it passes through stable RNA structural motifs without perturbation of velocity. Taken together, the results demonstrate that consistent, robust translocative behavior is a hallmark of group II intron-encoded RTs, some of which operate at high velocity., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Type III CRISPR-Cas systems have a unique mode of interference, involving crRNA-guided recognition of nascent RNA and leading to DNA and RNA degradation. How type III systems acquire new CRISPR spacers is currently not well understood. Here, we characterize CRISPR spacer uptake by a type III-A system within its native host, Streptococcus thermophilus. Adaptation by the type II-A system in the same host provided a basis for comparison. Cas1 and Cas2 proteins were critical for type III adaptation but deletion of genes responsible for crRNA biogenesis or interference did not detectably change spacer uptake patterns, except those related to host counter-selection. Unlike the type II-A system, type III spacers are acquired in a PAM- and orientation-independent manner. Interestingly, certain regions of plasmids and the host genome were particularly well-sampled during type III-A, but not type II-A, spacer uptake. These regions included the single-stranded origins of rolling-circle replicating plasmids, rRNA and tRNA encoding gene clusters, promoter regions of expressed genes and 5' UTR regions involved in transcription attenuation. These features share the potential to form DNA secondary structures, suggesting a preferred substrate for type III adaptation. Lastly, the type III-A system adapted to and protected host cells from lytic phage infection., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The RNA isoform repertoire is regulated by splicing factor (SF) expression, and alterations in SF levels are associated with disease. SFs contain ultraconserved poison exon (PE) sequences that exhibit greater identity across species than nearby coding exons, but their physiological role and molecular regulation is incompletely understood. We show that PEs in serine-arginine-rich (SR) proteins, a family of 14 essential SFs, are differentially spliced during induced pluripotent stem cell (iPSC) differentiation and in tumors versus normal tissues. We uncover an extensive cross-regulatory network of SR proteins controlling their expression via alternative splicing coupled to nonsense-mediated decay. We define sequences that regulate PE inclusion and protein expression of the oncogenic SF TRA2β using an RNA-targeting CRISPR screen. We demonstrate location dependency of RS domain activity on regulation of TRA2β-PE using CRISPR artificial SFs. Finally, we develop splice-switching antisense oligonucleotides to reverse the increased skipping of TRA2β-PE detected in breast tumors, altering breast cancer cell viability, proliferation, and migration., Competing Interests: Declaration of Interests A patent application concerning this work is in preparation., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE 1 and Roadmap Epigenomics 2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

The Encylopedia of DNA Elements (ENCODE) Project launched in 2003 with the long-term goal of developing a comprehensive map of functional elements in the human genome. These included genes, biochemical regions associated with gene regulation (for example, transcription factor binding sites, open chromatin, and histone marks) and transcript isoforms. The marks serve as sites for candidate cis-regulatory elements (cCREs) that may serve functional roles in regulating gene expression 1 . The project has been extended to model organisms, particularly the mouse. In the third phase of ENCODE, nearly a million and more than 300,000 cCRE annotations have been generated for human and mouse, respectively, and these have provided a valuable resource for the scientific community.

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome 1 . Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

CRISPR-Cas adaptive immune systems are used by prokaryotes to defend against invaders like viruses and other mobile genetic elements. Immune memories are stored in the form of 'spacers' which are short DNA sequences that are captured from invaders and added to the CRISPR array during a process called 'adaptation'. Spacers are transcribed and the resulting CRISPR (cr)RNAs assemble with different Cas proteins to form effector complexes that recognize matching nucleic acid and destroy it ('interference'). Adaptation can be 'naïve', i.e. independent of any existing spacer matches, or it can be 'primed', i.e. spurred by the crRNA-mediated detection of a complete or partial match to an invader sequence. Here we show that primed adaptation occurs in Pyrococcus furiosus. Although P. furiosus has three distinct CRISPR-Cas interference systems (I-B, I-A and III-B), only the I-B system and Cas3 were necessary for priming. Cas4, which is important for selection and processing of new spacers in naïve adaptation, was also essential for priming. Loss of either the I-B effector proteins or Cas3 reduced naïve adaptation. However, when Cas3 and all crRNP genes were deleted, uptake of correctly processed spacers was observed, indicating that none of these interference proteins are necessary for naïve adaptation., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Reverse transcriptase (RT) enzymes are indispensable tools for interrogating diverse aspects of RNA metabolism and transcriptome composition. Due to the growing interest in sequence and structural complexity of long RNA molecules, processive RT enzymes are now required for preserving linkage and information content in mixed populations of transcripts, and the low-processivity RT enzymes that are commercially available cannot meet this need. MarathonRT is encoded within a eubacterial group II intron, and it has been shown to efficiently copy highly structured long RNA molecules in a single pass. In this work, we systematically characterize MarathonRT as a tool enzyme and optimize its performance in a variety of applications that include single-cycle reverse transcription of long RNAs, dimethyl sulfate mutational profiling (DMS-MaP), selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP), using ultra-long amplicons and the detection of natural RNA base modifications. By diversifying MarathonRT reaction protocols, we provide an upgraded suite of tools for cutting-edge RNA research and clinical application., Competing Interests: Declaration of Competing Interest A patent application on MarathonRT has been filed by Yale University., (Copyright © 2020 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Background: A critical step in uncovering rules of RNA processing is to study the in vivo regulatory networks of RNA binding proteins (RBPs). Crosslinking and immunoprecipitation (CLIP) methods enable mapping RBP targets transcriptome-wide, but methodological differences present challenges to large-scale analysis across datasets. The development of enhanced CLIP (eCLIP) enabled the mapping of targets for 150 RBPs in K562 and HepG2, creating a unique resource of RBP interactomes profiled with a standardized methodology in the same cell types., Results: Our analysis of 223 eCLIP datasets reveals a range of binding modalities, including highly resolved positioning around splicing signals and mRNA untranslated regions that associate with distinct RBP functions. Quantification of enrichment for repetitive and abundant multicopy elements reveals 70% of RBPs have enrichment for non-mRNA element classes, enables identification of novel ribosomal RNA processing factors and sites, and suggests that association with retrotransposable elements reflects multiple RBP mechanisms of action. Analysis of spliceosomal RBPs indicates that eCLIP resolves AQR association after intronic lariat formation, enabling identification of branch points with single-nucleotide resolution, and provides genome-wide validation for a branch point-based scanning model for 3' splice site recognition. Finally, we show that eCLIP peak co-occurrences across RBPs enable the discovery of novel co-interacting RBPs., Conclusions: This work reveals novel insights into RNA biology by integrated analysis of eCLIP profiling of 150 RBPs with distinct functions. Further, our quantification of both mRNA and other element association will enable further research to identify novel roles of RBPs in regulating RNA processing.

We report here the complete genome sequence of Streptococcus ratti strain JH145. Streptococcus ratti is a cariogenic species of mutans streptococcus that has been isolated from rat and human teeth. The strain JH145, derived from strain BHT-2, is interesting for oral health because it does not produce cariogenic lactic acid but shows robust biofilm production., (Copyright © 2020 Garrett et al.)

CRISPR-Cas systems provide heritable immunity against viruses by capturing short invader DNA sequences, termed spacers, and incorporating them into the CRISPR loci of the prokaryotic host genome. Here, we investigate DNA elements that control accurate spacer uptake in the type II-A CRISPR locus of Streptococcus thermophilus. We determined that purified Cas1 and Cas2 proteins catalyze spacer integration with high specificity for CRISPR repeat junctions. We show that 10 bp of the CRISPR leader sequence is critical for stimulating polarized integration preferentially at the repeat proximal to the leader. Spacer integration proceeds through a two-step transesterification reaction where the 3' hydroxyl groups of the spacer target both repeat borders on opposite strands. The leader-proximal end of the repeat is preferentially targeted for the first site of integration through recognition of sequences spanning the leader-repeat junction. Subsequently, second-site integration at the leader-distal end of the repeat is specified by multiple determinants including a length-defining mechanism relying on a repeat element proximal to the second site of integration. Our results highlight the intrinsic ability of type II Cas1/Cas2 proteins to coordinate directional and site-specific spacer integration into the CRISPR locus to ensure precise duplication of the repeat required for CRISPR immunity., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Acquiring foreign spacer DNA into the CRISPR locus is an essential primary step of the CRISPR-Cas pathway in prokaryotes for developing host immunity to mobile genetic elements. Here, we investigate spacer integration in vitro using proteins from Pyrococcus furiosus and demonstrate that Cas1 and Cas2 are sufficient to accurately integrate spacers into a minimal CRISPR locus. Using high-throughput sequencing, we identified high frequency spacer integration occurring at the same CRISPR repeat border sites utilized in vivo, as well as at several non-CRISPR plasmid sequences which share features with repeats. Analysis of non-CRISPR integration sites revealed that Cas1 and Cas2 are directed to catalyze full-site spacer integration at specific DNA stretches where guanines and/or cytosines are 30 base pairs apart and the intervening sequence harbors several positionally conserved bases. Moreover, assaying a series of CRISPR repeat mutations, followed by sequencing of the integration products, revealed that the specificity of integration is primarily directed by sequences at the leader-repeat junction as well as an adenine-rich sequence block in the mid-repeat. Together, our results indicate that P. furiosus Cas1 and Cas2 recognize multiple sequence features distributed over a 30 base pair DNA region for accurate spacer integration at the CRISPR repeat., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Increasing evidence suggests that transcriptional control and chromatin activities at large involve regulatory RNAs, which likely enlist specific RNA-binding proteins (RBPs). Although multiple RBPs have been implicated in transcription control, it has remained unclear how extensively RBPs directly act on chromatin. We embarked on a large-scale RBP ChIP-seq analysis, revealing widespread RBP presence in active chromatin regions in the human genome. Like transcription factors (TFs), RBPs also show strong preference for hotspots in the genome, particularly gene promoters, where their association is frequently linked to transcriptional output. Unsupervised clustering reveals extensive co-association between TFs and RBPs, as exemplified by YY1, a known RNA-dependent TF, and RBM25, an RBP involved in splicing regulation. Remarkably, RBM25 depletion attenuates all YY1-dependent activities, including chromatin binding, DNA looping, and transcription. We propose that various RBPs may enhance network interaction through harnessing regulatory RNAs to control transcription., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Allele-specific protein-RNA binding is an essential aspect that may reveal functional genetic variants (GVs) mediating post-transcriptional regulation. Recently, genome-wide detection of in vivo binding of RNA-binding proteins is greatly facilitated by the enhanced crosslinking and immunoprecipitation (eCLIP) method. We developed a new computational approach, called BEAPR, to identify allele-specific binding (ASB) events in eCLIP-Seq data. BEAPR takes into account crosslinking-induced sequence propensity and variations between replicated experiments. Using simulated and actual data, we show that BEAPR largely outperforms often-used count analysis methods. Importantly, BEAPR overcomes the inherent overdispersion problem of these methods. Complemented by experimental validations, we demonstrate that the application of BEAPR to ENCODE eCLIP-Seq data of 154 proteins helps to predict functional GVs that alter splicing or mRNA abundance. Moreover, many GVs with ASB patterns have known disease relevance. Overall, BEAPR is an effective method that helps to address the outstanding challenge of functional interpretation of GVs.

Alternative splicing of pre-messenger RNA transcripts enables the generation of multiple protein isoforms from the same gene locus, providing a major source of protein diversity in mammalian genomes. RNA binding proteins (RBPs) bind to RNA to control splice site choice and define which exons are included in the resulting mature RNA transcript. However, depending on where the RBPs bind relative to splice sites, they can activate or repress splice site usage. To explore this position-specific regulation, in vivo binding sites identified by methods such as cross-linking and immunoprecipitation (CLIP) are integrated with alternative splicing events identified by RNA-seq or microarray. Merging these data sets enables the generation of a "splicing map," where CLIP signal relative to a merged meta-exon provides a simple summary of the position-specific effect of binding on splicing regulation. Here, we provide RBP-Maps, a software tool to simplify generation of these maps and enable researchers to rapidly query regulatory patterns of an RBP of interest. Further, we discuss various alternative approaches to generate such splicing maps, focusing on how decisions in construction (such as the use of peak versus read density, or whole-reads versus only single-nucleotide candidate crosslink positions) can affect the interpretation of these maps using example eCLIP data from the 150 RBPs profiled by the ENCODE consortium., (© 2019 Yee et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Adenosine-to-inosine (A-to-I) editing, mediated by the ADAR enzymes, diversifies the transcriptome by altering RNA sequences. Recent studies reported global changes in RNA editing in disease and development. Such widespread editing variations necessitate an improved understanding of the regulatory mechanisms of RNA editing. Here, we study the roles of >200 RNA-binding proteins (RBPs) in mediating RNA editing in two human cell lines. Using RNA-sequencing and global protein-RNA binding data, we identify a number of RBPs as key regulators of A-to-I editing. These RBPs, such as TDP-43, DROSHA, NF45/90 and Ro60, mediate editing through various mechanisms including regulation of ADAR1 expression, interaction with ADAR1, and binding to Alu elements. We highlight that editing regulation by Ro60 is consistent with the global up-regulation of RNA editing in systemic lupus erythematosus. Additionally, most key editing regulators act in a cell type-specific manner. Together, our work provides insights for the regulatory mechanisms of RNA editing., Competing Interests: G.W.Y. is a co-founder of Locana and Eclipse Bioinnovations and member of the scientific advisory boards of Locana, Eclipse Bioinnovations and Aquinnah Pharmaceuticals. E.V.N is a co-founder and member of the scientific advisory board of Eclipse BioInnovations. The terms of these arrangements have been reviewed and approved by the University of California, San Diego in accordance with its conflict of interest policies. The remaining authors declare no competing interests.

Translation and decay of eukaryotic mRNAs is controlled by shortening of the poly(A) tail and release of the poly(A)-binding protein Pab1/PABP. The Ccr4-Not complex contains two exonucleases-Ccr4 and Caf1/Pop2-that mediate mRNA deadenylation. Here, using a fully reconstituted biochemical system with proteins from the fission yeast Schizosaccharomyces pombe, we show that Pab1 interacts with Ccr4-Not, stimulates deadenylation, and differentiates the roles of the nuclease enzymes. Surprisingly, Pab1 release relies on Ccr4 activity. In agreement with this, in vivo experiments in budding yeast show that Ccr4 is a general deadenylase that acts on all mRNAs. In contrast, Caf1 only trims poly(A) not bound by Pab1. As a consequence, Caf1 is a specialized deadenylase required for the selective deadenylation of transcripts with lower rates of translation elongation and reduced Pab1 occupancy. These findings reveal a coupling between the rates of translation and deadenylation that is dependent on Pab1 and Ccr4-Not., (Copyright © 2018 MRC Laboratory of Molecular Biology. Published by Elsevier Inc. All rights reserved.)

RNA binding proteins (RBPs) orchestrate the production, processing, and function of mRNAs. Here, we present the affinity landscapes of 78 human RBPs using an unbiased assay that determines the sequence, structure, and context preferences of these proteins in vitro by deep sequencing of bound RNAs. These data enable construction of "RNA maps" of RBP activity without requiring crosslinking-based assays. We found an unexpectedly low diversity of RNA motifs, implying frequent convergence of binding specificity toward a relatively small set of RNA motifs, many with low compositional complexity. Offsetting this trend, however, we observed extensive preferences for contextual features distinct from short linear RNA motifs, including spaced "bipartite" motifs, biased flanking nucleotide composition, and bias away from or toward RNA structure. Our results emphasize the importance of contextual features in RNA recognition, which likely enable targeting of distinct subsets of transcripts by different RBPs that recognize the same linear motif., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

To achieve adaptive and heritable immunity against viruses and other mobile genetic elements, CRISPR-Cas systems must capture and store short DNA fragments (spacers) from these foreign elements into host genomic CRISPR arrays. This process is catalyzed by conserved Cas1/Cas2 integration complexes, but the specific roles of another highly conserved protein linked to spacer acquisition, the Cas4 nuclease, are just now emerging. Here, we show that two Cas4 nucleases (Cas4-1 and Cas4-2) play critical roles in CRISPR spacer acquisition in Pyrococcus furiosus. The nuclease activities of both Cas4 proteins are required to process protospacers to the correct size. Cas4-1 specifies the upstream PAM (protospacer adjacent motif), while Cas4-2 specifies the conserved downstream motif. Both Cas4 proteins ensure CRISPR spacer integration in a defined orientation leading to CRISPR immunity. Collectively, these findings provide in vivo evidence for critical roles of Cas4 nucleases in protospacer generation and functional spacer integration at CRISPR arrays., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Understanding the functional impact of genomic variants is a major goal of modern genetics and personalized medicine. Although many synonymous and non-coding variants act through altering the efficiency of pre-mRNA splicing, it is difficult to predict how these variants impact pre-mRNA splicing. Here, we describe a massively parallel approach we use to test the impact on pre-mRNA splicing of 2059 human genetic variants spanning 110 alternative exons. This method, called variant exon sequencing (Vex-seq), yields data that reinforce known mechanisms of pre-mRNA splicing, identifies variants that impact pre-mRNA splicing, and will be useful for increasing our understanding of genome function.

We report here the complete genome sequence of Streptococcus thermophilus DGCC 7710. S. thermophilus is widely used in industrial dairy production., (Copyright © 2018 Hatmaker et al.)

To acquire CRISPR-Cas immunity against invasive mobile genetic elements, prokaryotes must first integrate fragments of foreign DNA into their genomic CRISPR arrays for use in future invader silencing. Here, we found that the hyperthermophilic archaeaon, Pyrococcus furiosus, actively incorporates DNA fragments (spacers) from both plasmid (foreign) and host genome (self) sequences into its seven CRISPR loci. The majority of new spacers were derived from DNA immediately downstream from a 5'-CCN-3' protospacer adjacent motif (PAM) that is critical for invader targeting. Interestingly, spacers were preferentially acquired from genome or plasmid regions corresponding to active transposons, CRISPR loci, ribosomal RNA genes, rolling circle origins of replication, and areas where plasmids recombined with the host chromosome. A common feature of the highly sampled spacers is that they arise from DNA regions expected to undergo DNA nicking and/or double-strand breaks. Taken together with recent results from bacterial systems, our findings indicate that free DNA termini and PAMs are conserved features important for CRISPR spacer uptake in diverse prokaryotes and CRISPR-Cas systems. Moreover, lethal self-targeting by CRISPR systems may contribute to host genome stability by eliminating cells undergoing active transposon mobility or chromosomal uptake of autonomously replicating foreign mobile genetic elements., (© The Author(s) 2017. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The generation of inhibitory interneuron progenitors from human embryonic stem cells (ESCs) is of great interest due to their potential use in transplantation therapies designed to treat central nervous system disorders. The medial ganglionic eminence (MGE) is a transient embryonic structure in the ventral telencephalon that is a major source of cortical GABAergic inhibitory interneuron progenitors. These progenitors migrate tangentially to sites in the cortex and differentiate into a variety of interneuron subtypes, forming local synaptic connections with excitatory projection neurons to modulate activity of the cortical circuitry. The homeobox domain-containing transcription factor NKX2.1 is highly expressed in the MGE and pre-optic area of the ventral subpallium and is essential for specifying cortical interneuron fate. Using a combination of growth factor agonists and antagonists to specify ventral telencephalic fates, we previously optimized a protocol for the efficient generation of NKX2.1-positive MGE-like neural progenitors from human ESCs. To establish their identity, we now characterize the transcriptome of these MGE-like neural progenitors using RNA sequencing and demonstrate the capacity of these cells to differentiate into inhibitory interneurons in vitro using a neuron-astrocyte co-culture system. These data provide information on the potential origin of interneurons in the human brain.

RNA structure is intimately related to function, yet methods to identify base-paired RNA strands in a transcriptome-wide manner in cells have remained elusive. One recent paper in Cell and two in Molecular Cell describe related methods to identify RNA sequences that interact in living cells, setting the stage for breakthroughs in our understanding of RNA structure and function., (Copyright © 2016. Published by Elsevier Inc.)

Obtaining RNA-seq measurements involves a complex data analytical process with a large number of competing algorithms as options. There is much debate about which of these methods provides the best approach. Unfortunately, it is currently difficult to evaluate their performance due in part to a lack of sensitive assessment metrics. We present a series of statistical summaries and plots to evaluate the performance in terms of specificity and sensitivity, available as a R/Bioconductor package ( http://bioconductor.org/packages/rnaseqcomp ). Using two independent datasets, we assessed seven competing pipelines. Performance was generally poor, with two methods clearly underperforming and RSEM slightly outperforming the rest.

Transcriptome-wide maps of RNA binding protein (RBP)-RNA interactions by immunoprecipitation (IP)-based methods such as RNA IP (RIP) and crosslinking and IP (CLIP) are key starting points for evaluating the molecular roles of the thousands of human RBPs. A significant bottleneck to the application of these methods in diverse cell lines, tissues, and developmental stages is the availability of validated IP-quality antibodies. Using IP followed by immunoblot assays, we have developed a validated repository of 438 commercially available antibodies that interrogate 365 unique RBPs. In parallel, 362 short-hairpin RNA (shRNA) constructs against 276 unique RBPs were also used to confirm specificity of these antibodies. These antibodies can characterize subcellular RBP localization. With the burgeoning interest in the roles of RBPs in cancer, neurobiology, and development, these resources are invaluable to the broad scientific community. Detailed information about these resources is publicly available at the ENCODE portal (https://www.encodeproject.org/)., (Copyright © 2016 Elsevier Inc. All rights reserved.)

Spatial restriction of mRNA to distinct subcellular locations enables local regulation and synthesis of proteins. However, the organizing principles of mRNA localization remain poorly understood. Here we analyzed subcellular transcriptomes of neural projections and soma of primary mouse cortical neurons and two neuronal cell lines and found that alternative last exons (ALEs) often confer isoform-specific localization. Surprisingly, gene-distal ALE isoforms were four times more often localized to neurites than gene-proximal isoforms. Localized isoforms were induced during neuronal differentiation and enriched for motifs associated with muscleblind-like (Mbnl) family RNA-binding proteins. Depletion of Mbnl1 and/or Mbnl2 reduced localization of hundreds of transcripts, implicating Mbnls in localization of mRNAs to neurites. We provide evidence supporting a model in which the linkage between genomic position of ALEs and subcellular localization enables coordinated induction of localization-competent mRNA isoforms through a post-transcriptional regulatory program that is induced during differentiation and reversed in cellular reprogramming and cancer., (Copyright © 2016 Elsevier Inc. All rights reserved.)

In eukaryotic cells, RNAs exist as ribonucleoprotein particles (RNPs). Despite the importance of these complexes in many biological processes, including splicing, polyadenylation, stability, transportation, localization, and translation, their compositions are largely unknown. We affinity-purified 20 distinct RNA-binding proteins (RBPs) from cultured Drosophila melanogaster cells under native conditions and identified both the RNA and protein compositions of these RNP complexes. We identified "high occupancy target" (HOT) RNAs that interact with the majority of the RBPs we surveyed. HOT RNAs encode components of the nonsense-mediated decay and splicing machinery, as well as RNA-binding and translation initiation proteins. The RNP complexes contain proteins and mRNAs involved in RNA binding and post-transcriptional regulation. Genes with the capacity to produce hundreds of mRNA isoforms, ultracomplex genes, interact extensively with heterogeneous nuclear ribonuclear proteins (hnRNPs). Our data are consistent with a model in which subsets of RNPs include mRNA and protein products from the same gene, indicating the widespread existence of auto-regulatory RNPs. From the simultaneous acquisition and integrative analysis of protein and RNA constituents of RNPs, we identify extensive cross-regulatory and hierarchical interactions in post-transcriptional control., (© 2015 Stoiber et al.; Published by Cold Spring Harbor Laboratory Press.)

Alternative splicing is regulated by RNA binding proteins (RBPs) that recognize pre-mRNA sequence elements and activate or repress adjacent exons. Here, we used RNA interference and RNA-seq to identify splicing events regulated by 56 Drosophila proteins, some previously unknown to regulate splicing. Nearly all proteins affected alternative first exons, suggesting that RBPs play important roles in first exon choice. Half of the splicing events were regulated by multiple proteins, demonstrating extensive combinatorial regulation. We observed that SR and hnRNP proteins tend to act coordinately with each other, not antagonistically. We also identified a cross-regulatory network where splicing regulators affected the splicing of pre-mRNAs encoding other splicing regulators. This large-scale study substantially enhances our understanding of recent models of splicing regulation and provides a resource of thousands of exons that are regulated by 56 diverse RBPs., (© 2015 Brooks et al.; Published by Cold Spring Harbor Laboratory Press.)

In species with a heterogametic sex, population genetics theory predicts that DNA sequences on the X chromosome can evolve faster than comparable sequences on autosomes. Both neutral and nonneutral evolutionary processes can generate this pattern. Complex traits like gene expression are not predicted to have accelerated evolution by these theories, yet a "faster-X" pattern of gene expression divergence has recently been reported for both Drosophila and mammals. Here, we test the hypothesis that accelerated adaptive evolution of cis-regulatory sequences on the X chromosome is responsible for this pattern by comparing the relative contributions of cis- and trans-regulatory changes to patterns of faster-X expression divergence observed between strains and species of Drosophila with a range of divergence times. We find support for this hypothesis, especially among male-biased genes, when comparing different species. However, we also find evidence that trans-regulatory differences contribute to a faster-X pattern of expression divergence both within and between species. This contribution is surprising because trans-acting regulators of X-linked genes are generally assumed to be randomly distributed throughout the genome. We found, however, that X-linked transcription factors appear to preferentially regulate expression of X-linked genes, providing a potential mechanistic explanation for this result. The contribution of trans-regulatory variation to faster-X expression divergence was larger within than between species, suggesting that it is more likely to result from neutral processes than positive selection. These data show how accelerated evolution of both coding and noncoding sequences on the X chromosome can lead to accelerated expression divergence on the X chromosome relative to autosomes., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Short-read high-throughput RNA sequencing, though powerful, is limited in its ability to directly measure exon connectivity in mRNAs that contain multiple alternative exons located farther apart than the maximum read length. Here, we use the Oxford Nanopore MinION sequencer to identify 7,899 'full-length' isoforms expressed from four Drosophila genes, Dscam1, MRP, Mhc, and Rdl. These results demonstrate that nanopore sequencing can be used to deconvolute individual isoforms and that it has the potential to be a powerful method for comprehensive transcriptome characterization.

CRISPR-Cas immune systems function to defend prokaryotes against potentially harmful mobile genetic elements including viruses and plasmids. The multiple CRISPR-Cas systems (Types I, II, and III) each target destruction of foreign nucleic acids via structurally and functionally diverse effector complexes (crRNPs). CRISPR-Cas effector complexes are comprised of CRISPR RNAs (crRNAs) that contain sequences homologous to the invading nucleic acids and Cas proteins specific to each immune system type. We have previously characterized a crRNP in Pyrococcus furiosus (Pfu) that contains Cmr (Type III-B) Cas proteins associated with one of two size classes of crRNAs and cleaves complementary target RNAs. Here, we have isolated and characterized two additional native Pfu crRNPs containing either Csa (Type I-A) or Cst (Type I-G) Cas proteins and distinct profiles of associated crRNAs. For each complex, the Cas proteins were identified by mass spectrometry and immunoblotting and the crRNAs by RNA sequencing and Northern blot analysis. The crRNAs associated with both the Csa and Cst complexes originate from all seven Pfu CRISPR loci and contain identical 5' ends (8-nt repeat-derived 5' tag sequences) but heterogeneous 3' ends (containing variable amounts of downstream repeat sequences). These crRNA forms are distinct from Cmr-associated crRNAs, indicating different 3' end processing pathways following primary cleavage of common pre-crRNAs. Like other previously characterized Type I CRISPR-Cas effector complexes, we predict that the newly identified Pfu Csa and Cst crRNPs each function to target invading DNA, adding an additional layer of protection beyond that afforded by the previously characterized RNA targeting Cmr complex., (© 2015 Majumdar et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Recursive splicing is a process in which large introns are removed in multiple steps by re-splicing at ratchet points--5' splice sites recreated after splicing. Recursive splicing was first identified in the Drosophila Ultrabithorax (Ubx) gene and only three additional Drosophila genes have since been experimentally shown to undergo recursive splicing. Here we identify 197 zero nucleotide exon ratchet points in 130 introns of 115 Drosophila genes from total RNA sequencing data generated from developmental time points, dissected tissues and cultured cells. The sequential nature of recursive splicing was confirmed by identification of lariat introns generated by splicing to and from the ratchet points. We also show that recursive splicing is a constitutive process, that depletion of U2AF inhibits recursive splicing, and that the sequence and function of ratchet points are evolutionarily conserved in Drosophila. Finally, we identify four recursively spliced human genes, one of which is also recursively spliced in Drosophila. Together, these results indicate that recursive splicing is commonly used in Drosophila, occurs in humans, and provides insight into the mechanisms by which some large introns are removed.

Many RNAs, including pre-mRNAs and long non-coding RNAs, can be thousands of nucleotides long and undergo complex post-transcriptional processing. Multiple sites of alternative splicing within a single gene exponentially increase the number of possible spliced isoforms, with most human genes currently estimated to express at least ten. To understand the mechanisms underlying these complex isoform expression patterns, methods are needed that faithfully maintain long-range exon connectivity information in individual RNA molecules. In this study, we describe SeqZip, a methodology that uses RNA-templated DNA-DNA ligation to retain and compress connectivity between distant sequences within single RNA molecules. Using this assay, we test proposed coordination between distant sites of alternative exon utilization in mouse Fn1, and we characterize the extraordinary exon diversity of Drosophila melanogaster Dscam1.