HITS-CLIP and integrative modeling define the Rbfox splicing-regulatory network linked to brain development and autism

Convergent Evidence That ZNF804A Is a Regulator of Pre-messenger RNA Processing and Gene Expression

Genome-wide association studies have linked common variation in ZNF804A with an increased risk of schizophrenia. However, little is known about the biology of ZNF804A and its role in schizophrenia. Here, we investigate the function of ZNF804A using a variety of complementary molecular techniques. We show that ZNF804A is a nuclear protein that interacts with neuronal RNA splicing factors and RNA-binding proteins including RBFOX1, which is also associated with schizophrenia, CELF3/4, components of the ubiquitin-proteasome system and the ZNF804A paralog, GPATCH8. GPATCH8 also interacts with splicing factors and is localized to nuclear speckles indicative of a role in pre-messenger RNA (mRNA) processing. Sequence analysis showed that GPATCH8 contains ultraconserved, alternatively spliced poison exons that are also regulated by RBFOX proteins. ZNF804A knockdown in SH-SY5Y cells resulted in robust changes in gene expression and pre-mRNA splicing converging on pathways associated with nervous system development, synaptic contact, and cell adhesion. We observed enrichment (P = 1.66 × 10-9) for differentially spliced genes in ZNF804A-depleted cells among genes that contain RBFOX-dependent alternatively spliced exons. Differentially spliced genes in ZNF804A-depleted cells were also enriched for genes harboring de novo loss of function mutations in autism spectrum disorder (P = 6.25 × 10-7, enrichment 2.16) and common variant alleles associated with schizophrenia (P = .014), bipolar disorder and schizophrenia (P = .003), and autism spectrum disorder (P = .005). These data suggest that ZNF804A and its paralogs may interact with neuronal-splicing factors and RNA-binding proteins to regulate the expression of a subset of synaptic and neurodevelopmental genes.

Dynamic methylome of internal mRNA N7-methylguanosine and its regulatory role in translation

Over 150 types of RNA modifications are identified in RNA molecules. Transcriptome profiling is one of the key steps in decoding the epitranscriptomic panorama of these chemical modifications and their potential functions. N⁷-methylguanosine (m⁷G) is one of the most abundant modifications present in tRNA, rRNA and mRNA 5'cap, and has critical roles in regulating RNA processing, metabolism and function. Besides its presence at the cap position in mRNAs, m⁷G is also identified in internal mRNA regions. However, its transcriptome-wide distribution and dynamic regulation within internal mRNA regions remain unknown. Here, we have established m⁷G individual-nucleotide-resolution cross-linking and immunoprecipitation with sequencing (m⁷G miCLIP-seq) to specifically detect internal mRNA m⁷G modification. Using this approach, we revealed that m⁷G is enriched at the 5'UTR region and AG-rich contexts, a feature that is well-conserved across different human/mouse cell lines and mouse tissues. Strikingly, the internal m⁷G modification is dynamically regulated under both H₂O₂ and heat shock treatments, with remarkable accumulations in the CDS and 3'UTR regions, and functions in promoting mRNA translation efficiency. Consistently, a PCNA 3'UTR minigene reporter harboring the native m⁷G modification site displays both enriched m⁷G modification and increased mRNA translation upon H₂O₂ treatment compared to the m⁷G site-mutated minigene reporter (G to A). Taken together, our findings unravel the dynamic profiles of internal mRNA m⁷G methylome and highlight m⁷G as a novel epitranscriptomic marker with regulatory roles in translation.

Disruptive variants of CSDE1 associate with autism and interfere with neuronal development and synaptic transmission

RNA binding proteins are key players in posttranscriptional regulation and have been implicated in neurodevelopmental and neuropsychiatric disorders. Here, we report a significant burden of heterozygous, likely gene-disrupting variants in CSDE1 (encoding a highly constrained RNA binding protein) among patients with autism and related neurodevelopmental disabilities. Analysis of 17 patients identifies common phenotypes including autism, intellectual disability, language and motor delay, seizures, macrocephaly, and variable ocular abnormalities. HITS-CLIP revealed that Csde1-binding targets are enriched in autism-associated gene sets, especially FMRP targets, and in neuronal development and synaptic plasticity-related pathways. Csde1 knockdown in primary mouse cortical neurons leads to an overgrowth of the neurites and abnormal dendritic spine morphology/synapse formation and impaired synaptic transmission, whereas mutant and knockdown experiments in Drosophila result in defects in synapse growth and synaptic transmission. Our study defines a new autism-related syndrome and highlights the functional role of CSDE1 in synapse development and synaptic transmission.

Profiling allele-specific gene expression in brains from individuals with autism spectrum disorder reveals preferential minor allele usage

One fundamental but understudied mechanism of gene regulation in disease is allele-specific expression (ASE), the preferential expression of one allele. We leveraged RNA-sequencing data from human brain to assess ASE in autism spectrum disorder (ASD). When ASE is observed in ASD, the allele with lower population frequency (minor allele) is preferentially more highly expressed than the major allele, opposite to the canonical pattern. Importantly, genes showing ASE in ASD are enriched in those downregulated in ASD postmortem brains and in genes harboring de novo mutations in ASD. Two regions, 14q32 and 15q11, containing all known orphan C/D box small nucleolar RNAs (snoRNAs), are particularly enriched in shifts to higher minor allele expression. We demonstrate that this allele shifting enhances snoRNA-targeted splicing changes in ASD-related target genes in idiopathic ASD and 15q11-q13 duplication syndrome. Together, these results implicate allelic imbalance and dysregulation of orphan C/D box snoRNAs in ASD pathogenesis.

Rbfox1 up-regulation impairs BDNF-dependent hippocampal LTP by dysregulating TrkB isoform expression levels

Brain-derived neurotrophic factor (BDNF) is a potent modulator of brain synaptic plasticity. Signaling defects caused by dysregulation of its Ntrk2 (TrkB) kinase (TrkB.FL) and truncated receptors (TrkB.T1) have been linked to the pathophysiology of several neurological and neurodegenerative disorders. We found that upregulation of Rbfox1, an RNA binding protein associated with intellectual disability, epilepsy and autism, increases selectively hippocampal TrkB.T1 isoform expression. Physiologically, increased Rbfox1 impairs BDNF-dependent LTP which can be rescued by genetically restoring TrkB.T1 levels. RNA-seq analysis of hippocampi with upregulation of Rbfox1 in conjunction with the specific increase of TrkB.T1 isoform expression also shows that the genes affected by Rbfox1 gain of function are surprisingly different from those influenced by Rbfox1 deletion. These findings not only identify TrkB as a major target of Rbfox1 pathophysiology but also suggest that gain or loss of function of Rbfox1 regulate different genetic landscapes.

A Large Panel of Isogenic APP and PSEN1 Mutant Human iPSC Neurons Reveals Shared Endosomal Abnormalities Mediated by APP β-CTFs, Not Aβ

Familial Alzheimer's disease (fAD) results from mutations in the amyloid precursor protein (APP) and presenilin (PSEN1 and PSEN2) genes. Here we leveraged recent advances in induced pluripotent stem cell (iPSC) and CRISPR/Cas9 genome editing technologies to generate a panel of isogenic knockin human iPSC lines carrying APP and/or PSEN1 mutations. Global transcriptomic and translatomic profiling revealed that fAD mutations have overlapping effects on the expression of AD-related and endocytosis-associated genes. Mutant neurons also increased Rab5+ early endosome size. APP and PSEN1 mutations had discordant effects on Aβ production but similar effects on APP β C-terminal fragments (β-CTFs), which accumulate in all mutant neurons. Importantly, endosomal dysfunction correlated with accumulation of β-CTFs, not Aβ, and could be rescued by pharmacological modulation of β-secretase (BACE). These data display the utility of our mutant iPSCs in studying AD-related phenotypes in a non-overexpression human-based system and support mounting evidence that β-CTF may be critical in AD pathogenesis.

Inherited and De Novo Genetic Risk for Autism Impacts Shared Networks

We performed a comprehensive assessment of rare inherited variation in autism spectrum disorder (ASD) by analyzing whole-genome sequences of 2,308 individuals from families with multiple affected children. We implicate 69 genes in ASD risk, including 24 passing genome-wide Bonferroni correction and 16 new ASD risk genes, most supported by rare inherited variants, a substantial extension of previous findings. Biological pathways enriched for genes harboring inherited variants represent cytoskeletal organization and ion transport, which are distinct from pathways implicated in previous studies. Nevertheless, the de novo and inherited genes contribute to a common protein-protein interaction network. We also identified structural variants (SVs) affecting non-coding regions, implicating recurrent deletions in the promoters of DLG2 and NR3C2. Loss of nr3c2 function in zebrafish disrupts sleep and social function, overlapping with human ASD-related phenotypes. These data support the utility of studying multiplex families in ASD and are available through the Hartwell Autism Research and Technology portal.

Tissue- and cell-specific expression of a splice variant in the II-III cytoplasmic loop of Cacna1b

Presynaptic Ca_V2.2 (N‐type) channels are fundamental for transmitter release across the nervous system. The gene encoding Ca_V2.2 channels, Cacna1b, contains alternatively spliced exons that result in functionally distinct splice variants (e18a, e24a, e31a, and 37a/37b). Alternative splicing of the cassette exon 18a generates two mRNA transcripts (+e18a‐Cacna1b and ∆e18a‐Cacna1b). In this study, using novel mouse genetic models and in situ hybridization (BaseScope™), we confirmed that +e18a‐Cacna1b splice variants are expressed in monoaminergic regions of the midbrain. We expanded these studies and identified +e18a‐Cacna1b mRNA in deep cerebellar cells and spinal cord motor neurons. Furthermore, we determined that +e18a‐Cacna1b is enriched in cholecystokinin‐expressing interneurons. Our results provide key information to understand cell‐specific functions of Ca_V2.2 channels.

AGO CLIP Reveals an Activated Network for Acute Regulation of Brain Glutamate Homeostasis in Ischemic Stroke

Post-transcriptional regulation by microRNAs (miRNAs) is essential for complex molecular responses to physiological insult and disease. Although many disease-associated miRNAs are known, their global targets and culminating network effects on pathophysiology remain poorly understood. We applied Argonaute (AGO) crosslinking immunoprecipitation (CLIP) to systematically elucidate altered miRNA-target interactions in brain following ischemia and reperfusion (I/R) injury. Among 1,190 interactions identified, the most prominent was the cumulative loss of target regulation by miR-29 family members. Integration of translational and time-course RNA profiles revealed a dynamic mode of miR-29 target de-regulation, led by acute translational activation and a later increase in RNA levels, allowing rapid proteomic changes to take effect. These functional regulatory events rely on canonical and non-canonical miR-29 binding and engage glutamate reuptake signals, such as glial glutamate transporter (GLT-1), to control local glutamate levels. These results uncover a miRNA target network that acts acutely to maintain brain homeostasis after ischemic stroke.

Genetic Associations between Voltage-Gated Calcium Channels and Psychiatric Disorders

Psychiatric disorders are mental, behavioral or emotional disorders. These conditions are prevalent, one in four adults suffer from any type of psychiatric disorders world-wide. It has always been observed that psychiatric disorders have a genetic component, however, new methods to sequence full genomes of large cohorts have identified with high precision genetic risk loci for these conditions. Psychiatric disorders include, but are not limited to, bipolar disorder, schizophrenia, autism spectrum disorder, anxiety disorders, major depressive disorder, and attention-deficit and hyperactivity disorder. Several risk loci for psychiatric disorders fall within genes that encode for voltage-gated calcium channels (CaVs). Calcium entering through CaVs is crucial for multiple neuronal processes. In this review, we will summarize recent findings that link CaVs and their auxiliary subunits to psychiatric disorders. First, we will provide a general overview of CaVs structure, classification, function, expression and pharmacology. Next, we will summarize tools to study risk loci associated with psychiatric disorders. We will examine functional studies of risk variations in CaV genes when available. Finally, we will review pharmacological evidence of the use of CaV modulators to treat psychiatric disorders. Our review will be of interest for those studying pathophysiological aspects of CaVs.

Rbfox1 Regulates Synaptic Transmission through the Inhibitory Neuron-Specific vSNARE Vamp1

Dysfunction of the neuronal RNA binding protein RBFOX1 has been linked to epilepsy and autism spectrum disorders. Rbfox1 loss in mice leads to neuronal hyper-excitability and seizures, but the physiological basis for this is unknown. We identify the vSNARE protein Vamp1 as a major Rbfox1 target. Vamp1 is strongly downregulated in Rbfox1 Nes-cKO mice due to loss of 3' UTR binding by RBFOX1. Cytoplasmic Rbfox1 stimulates Vamp1 expression in part by blocking microRNA-9. We find that Vamp1 is specifically expressed in inhibitory neurons, and that both Vamp1 knockdown and Rbfox1 loss lead to decreased inhibitory synaptic transmission and E/I imbalance. Re-expression of Vamp1 selectively within interneurons rescues the electrophysiological changes in the Rbfox1 cKO, indicating that Vamp1 loss is a major contributor to the Rbfox1 Nes-cKO phenotype. The regulation of interneuron-specific Vamp1 by Rbfox1 provides a paradigm for broadly expressed RNA-binding proteins performing specialized functions in defined neuronal subtypes.

Monitored eCLIP: high accuracy mapping of RNA-protein interactions

CLIP-seq methods provide transcriptome-wide snapshots of RNA-protein interactions in live cells. Reverse transcriptases stopping at cross-linked nucleotides sign for RNA-protein binding sites. Reading through cross-linked positions results in false binding site assignments. In the ‘monitored enhanced CLIP’ (meCLIP) method, a barcoded biotinylated linker is ligated at the 5′ end of cross-linked RNA fragments to purify RNA prior to the reverse transcription. cDNAs keeping the barcode sequence correspond to reverse transcription read-throughs. Read through occurs in unpredictable proportions, representing up to one fourth of total reads. Filtering out those reads strongly improves reliability and precision in protein binding site assignment.

DM-RPIs: Predicting ncRNA-protein interactions using stacked ensembling strategy

ncRNA-protein interactions (ncRPIs) play an important role in a number of cellular processes, such as post-transcriptional modification, transcriptional regulation, disease progression and development. Since experimental methods are expensive and time-consuming to identify the ncRPIs, we proposed a computational method, Deep Mining ncRNA-Protein Interactions (DM-RPIs), for identifying the ncRPIs. In order to descending dimension and excavating hidden information from k-mer frequency of RNA and protein sequences, using the Deep Stacking Auto-encoders Networks (DSANs) model refined the raw data. Three common machine learning algorithms, Support Vector Machine (SVM), Random Forest (RF), and Convolution Neural Network (CNN), were separately trained as individual predictors and then the three individual predictors were integrated together using stacked ensembling strategy. Based on the RPI2241 dataset, DM-RPI obtains an accuracy of 0.851, precision of 0.852, sensitivity of 0.873, specificity of 0.826, and MCC of 0.701, which is promising and pioneering for the prediction of ncRPIs.

CLIPick: a sensitive peak caller for expression-based deconvolution of HITS-CLIP signals

High-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP, also called CLIP-Seq) has been used to map global RNA–protein interactions. However, a critical caveat of HITS-CLIP results is that they contain non-linear background noise—different extent of non-specific interactions caused by individual transcript abundance—that has been inconsiderately normalized, resulting in sacrifice of sensitivity. To properly deconvolute RNA–protein interactions, we have implemented CLIPick, a flexible peak calling pipeline for analyzing HITS-CLIP data, which statistically determines the signal-to-noise ratio for each transcript based on the expression-dependent background simulation. Comprising of streamlined Python modules with an easy-to-use standalone graphical user interface, CLIPick robustly identifies significant peaks and quantitatively defines footprint regions within which RNA–protein interactions were occurred. CLIPick outperforms other peak callers in accuracy and sensitivity, selecting the largest number of peaks particularly in lowly expressed transcripts where such marginal signals are hard to discriminate. Specifically, the application of CLIPick to Argonaute (Ago) HITS-CLIP data were sensitive enough to uncover extended features of microRNA target sites, and these sites were experimentally validated. CLIPick enables to resolve critical interactions in a wide spectrum of transcript levels and extends the scope of HITS-CLIP analysis. CLIPick is available at: http://clip.korea.ac.kr/clipick/

Modeling RNA-Binding Protein Specificity In Vivo by Precisely Registering Protein-RNA Crosslink Sites

RNA-binding proteins (RBPs) regulate post-transcriptional gene expression by recognizing short and degenerate sequence motifs in their target transcripts, but precisely defining their binding specificity remains challenging. Crosslinking and immunoprecipitation (CLIP) allows for mapping of the exact protein-RNA crosslink sites, which frequently reside at specific positions in RBP motifs at single-nucleotide resolution. Here, we have developed a computational method, named mCross, to jointly model RBP binding specificity while precisely registering the crosslinking position in motif sites. We applied mCross to 112 RBPs using ENCODE eCLIP data and validated the reliability of the discovered motifs by genome-wide analysis of allelic binding sites. Our analyses revealed that the prototypical SR protein SRSF1 recognizes clusters of GGA half-sites in addition to its canonical GGAGGA motif. Therefore, SRSF1 regulates splicing of a much larger repertoire of transcripts than previously appreciated, including HNRNPD and HNRNPDL, which are involved in multivalent protein assemblies and phase separation.

A mechanism underlying position-specific regulation of alternative splicing

Many RNA-binding proteins including a master regulator of splicing in developing brain and muscle, polypyrimidine tract-binding protein 1 (PTBP1), can either activate or repress alternative exons depending on the pre-mRNA recruitment position. When bound upstream or within regulated exons PTBP1 tends to promote their skipping, whereas binding to downstream sites often stimulates inclusion. How this switch is orchestrated at the molecular level is poorly understood. Using bioinformatics and biochemical approaches we show that interaction of PTBP1 with downstream intronic sequences can activate natural cassette exons by promoting productive docking of the spliceosomal U1 snRNP to a suboptimal 5' splice site. Strikingly, introducing upstream PTBP1 sites to this circuitry leads to a potent splicing repression accompanied by the assembly of an exonic ribonucleoprotein complex with a tightly bound U1 but not U2 snRNP. Our data suggest a molecular mechanism underlying the transition between a better-known repressive function of PTBP1 and its role as a bona fide splicing activator. More generally, we argue that the functional outcome of individual RNA contacts made by an RNA-binding protein is subject to extensive context-specific modulation.

Discovering and Mapping the Modified Nucleotides That Comprise the Epitranscriptome of mRNA

An important mechanism of gene expression regulation is the regulated modification of nucleotides in messenger RNA (mRNA). These modified nucleotides affect mRNA translation, stability, splicing, and other processes. A cluster of nucleotide modifications is found adjacent to the mRNA cap structure and another set can be found internally within transcripts. The most prominent modifications are methylations of adenosine to form either N ⁶-methyladenosine (m⁶A), an internal modified nucleotide, or N ⁶,2'-O-dimethyladenosine (m⁶Am), which is found exclusively at the first templated nucleotide of certain mRNAs. In addition, other rare modified nucleotides have been identified and together these form the epitranscriptomic code of mRNA. In the case of some modified nucleotides, the presence, location, or abundance is a subject of debate. Here, we review the methods that enable the discovery of modified nucleotides and how these approaches can be used to map epitranscriptomic modifications in mRNA.

A Bayesian framework that integrates multi-omics data and gene networks predicts risk genes from schizophrenia GWAS data

Genome-wide association studies (GWAS) have identified >100 schizophrenia (SCZ)-associated loci, but using these findings to illuminate disease biology remains a challenge. Here, we present integrative RIsk Gene Selector (iRIGS), a Bayesian framework that integrates multi-omics data and gene networks to infer risk genes in GWAS loci. By applying iRIGS to SCZ GWAS data, we predicted a set of high-confidence risk genes (HRGs), most of which are not the nearest genes to the GWAS index variants. HRGs account for a significantly enriched heritability estimated by stratified LD-score regression. Moreover, HRGs are predominantly expressed in brain tissues, especially prenatally, and are enriched for targets of approved drugs, suggesting opportunities to reposition existing drugs for SCZ. Thus, iRIGS can leverage accumulating functional genomics and GWAS data to advance understanding of SCZ etiology and potential therapeutics.

Defining the Genetic, Genomic, Cellular, and Diagnostic Architectures of Psychiatric Disorders

Studies of the genetics of psychiatric disorders have become one of the most exciting and fast-moving areas in human genetics. A decade ago, there were few reproducible findings, and now there are hundreds. In this review, we focus on the findings that have illuminated the genetic architecture of psychiatric disorders and the challenges of using these findings to inform our understanding of pathophysiology. The evidence is now overwhelming that psychiatric disorders are "polygenic"-that many genetic loci contribute to risk. With the exception of a subset of those with ASD, few individuals with a psychiatric disorder have a single, deterministic genetic cause; rather, developing a psychiatric disorder is influenced by hundreds of different genetic variants, consistent with a polygenic model. As progressively larger studies have uncovered more about their genetic architecture, the need to elucidate additional architectures has become clear. Even if we were to have complete knowledge of the genetic architecture of a psychiatric disorder, full understanding requires deep knowledge of the functional genomic architecture-the implicated loci impact regulatory processes that influence gene expression and the functional coordination of genes that control biological processes. Following from this is cellular architecture: of all brain regions, cell types, and developmental stages, where and when are the functional architectures operative? Given that the genetic architectures of different psychiatric disorders often strongly overlap, we are challenged to re-evaluate and refine the diagnostic architectures of psychiatric disorders using fundamental genetic and neurobiological data.

Methods to study RNA-protein interactions

Noncoding RNA sequences, including long noncoding RNAs, small nucleolar RNAs, and untranslated mRNA regions, accomplish many of their diverse functions through direct interactions with RNA-binding proteins (RBPs). Recent efforts have identified hundreds of new RBPs that lack known RNA-binding domains, thus underscoring the complexity and diversity of RNA-protein complexes. Recent progress has expanded the number of methods for studying RNA-protein interactions in two general categories: approaches that characterize proteins bound to an RNA of interest (RNA-centric), and those that examine RNAs bound to a protein of interest (protein-centric). Each method has unique strengths and limitations, which makes it important to select optimal approaches for the biological question being addressed. Here we review methods for the study of RNA-protein interactions, with a focus on their suitability for specific applications.

Conserved functions of RNA-binding proteins in muscle

Animals require different types of muscle for survival, for example for circulation, motility, reproduction and digestion. Much emphasis in the muscle field has been placed on understanding how transcriptional regulation generates diverse types of muscle during development. Recent work indicates that alternative splicing and RNA regulation are as critical to muscle development, and altered function of RNA-binding proteins causes muscle disease. Although hundreds of genes predicted to bind RNA are expressed in muscles, many fewer have been functionally characterized. We present a cross-species view summarizing what is known about RNA-binding protein function in muscle, from worms and flies to zebrafish, mice and humans. In particular, we focus on alternative splicing regulated by the CELF, MBNL and RBFOX families of proteins. We discuss the systemic nature of diseases associated with loss of RNA-binding proteins in muscle, focusing on mis-regulation of CELF and MBNL in myotonic dystrophy. These examples illustrate the conservation of RNA-binding protein function and the marked utility of genetic model systems in understanding mechanisms of RNA regulation.

RBP-Maps enables robust generation of splicing regulatory maps

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.

Toward Identifying Subnetworks from FBF Binding Landscapes in Caenorhabditis Spermatogenic or Oogenic Germlines

Metazoan PUF (Pumilio and FBF) RNA-binding proteins regulate various biological processes, but a common theme across phylogeny is stem cell regulation. In Caenorhabditis elegans, FBF (fem-3 Binding Factor) maintains germline stem cells regardless of which gamete is made, but FBF also functions in the process of spermatogenesis. We have begun to "disentangle" these biological roles by asking which FBF targets are gamete-independent, as expected for stem cells, and which are gamete-specific. Specifically, we compared FBF iCLIP binding profiles in adults making sperm to those making oocytes. Normally, XX adults make oocytes. To generate XX adults making sperm, we used a fem-3(gf) mutant requiring growth at 25°; for comparison, wild-type oogenic hermaphrodites were also raised at 25°. Our FBF iCLIP data revealed FBF binding sites in 1522 RNAs from oogenic adults and 1704 RNAs from spermatogenic adults. More than half of these FBF targets were independent of germline gender. We next clustered RNAs by FBF-RNA complex frequencies and found four distinct blocks. Block I RNAs were enriched in spermatogenic germlines, and included validated target fog-3, while Block II and III RNAs were common to both genders, and Block IV RNAs were enriched in oogenic germlines. Block II (510 RNAs) included almost all validated FBF targets and was enriched for cell cycle regulators. Block III (21 RNAs) was enriched for RNA-binding proteins, including previously validated FBF targets gld-1 and htp-1 We suggest that Block I RNAs belong to the FBF network for spermatogenesis, and that Blocks II and III are associated with stem cell functions.

Differential NOVA2-Mediated Splicing in Excitatory and Inhibitory Neurons Regulates Cortical Development and Cerebellar Function

RNA-binding proteins (RBPs) regulate genetic diversity, but the degree to which they do so in individual cell types in vivo is unknown. We developed NOVA2 cTag-crosslinking and immunoprecipitation (CLIP) to generate functional RBP-RNA maps from different neuronal populations in the mouse brain. Combining cell type datasets from Nova2-cTag and Nova2 conditional knockout mice revealed differential NOVA2 regulatory actions on alternative splicing (AS) on the same transcripts expressed in different neurons. This includes functional differences in transcripts expressed in cortical and cerebellar excitatory versus inhibitory neurons, where we find NOVA2 is required for, respectively, development of laminar structure, motor coordination, and synapse formation. We also find that NOVA2-regulated AS is coupled to NOVA2 regulation of intron retention in hundreds of transcripts, which can sequester the trans-acting splicing factor PTBP2. In summary, cTag-CLIP complements single-cell RNA sequencing (RNA-seq) studies by providing a means for understanding RNA regulation of functional cell diversity.

Expression of RNA-binding protein Rbfox1l demarcates a restricted population of dorsal telencephalic neurons within the adult zebrafish brain

Rbfox RNA-binding proteins are expressed in the adult mammalian brain and are required for proper brain development and function. Studies in mice and humans have implicated Rbfox1/RBFOX1 in autism, neuronal excitation and epilepsy, and Rbfox2/RBFOX2 in cerebellar development. The zebrafish has emerged as a prominent model system for brain study, possessing neuroanatomical conservation with mammals and an extensive capacity for adult neurogenesis and plasticity. In this study, we characterize Rbfox1l and Rbfox2 expression in the adult zebrafish brain. While Rbfox2 is expressed broadly, Rbfox1l is expressed in restricted populations of neurons in the dorsal telencephalon and cerebellum. In the dorsal telencephalon, Rbfox1l is expressed in a specific population of neurons spanning Dm and Dc regions. In the cerebellum, Rbfox1l and Rbfox2 are expressed in the Purkinje cell layer, reminiscent of Rbfox1 and Rbfox2 expression in the mammalian cerebellum. Our findings motivate future studies of Rbfox function in the zebrafish brain.

Evaluation of Post-transcriptional Gene Regulation in Pancreatic Cancer Cells: Studying RNA Binding Proteins and Their mRNA Targets

Post-transcriptional regulation of gene expression through interaction between RNA binding proteins (RBPs) and target mRNAs have gained considerable interest over the last decade. Altered expression of RBPs as detected in pancreatic ductal adenocarcinoma (PDAC) cells alters mRNA processing, and in turn, the entire transcriptome and proteome. Thus, this gene regulatory mechanism can regulate important pro-oncogenic signaling pathways (e.g., TP53, WEE1, and c-MYC) in PDAC cells. Ribonucleoprotein immunoprecipitation assays (RNP-IP or RIP) are a modified immunoprecipitation method to study physical interactions between RBPs and their mRNA targets. As a first step to explore RBP interactomes and define novel therapeutic targets and dysregulated pathways in disease, RIPs are a sensitive and established molecular biology technique used to isolate and differentiate bound transcripts to RBPs in a variety of experimental conditions. This chapter describes an up-to-date, detailed protocol for performing this assay in mammalian cytoplasmic extracts (i.e., PDAC cells), and reviews current methods to validate target binding sites such as electrophoretic mobility shift assay (EMSA) and cross-linking immunoprecipitation polymerase chain reaction (CLIP-PCR).

Polygenic risk score for schizophrenia is not strongly associated with the expression of specific genes or gene sets

BACKGROUND

The polygenic risk score (PRS) is derived from single nucleotide polymorphisms (SNPs) including those that are genome-wide significant and also including a large number of others more weakly associated with schizophrenia. Such variants are widely dispersed, though concentrated near genes expressed in the brain, and it has been proposed that these SNP associations result from impacts on cell regulatory networks that ultimately affect the expression or function of a modest number of 'core' genes. A previous study showed association of some genome-wide association study-significant variants with expression of a number of genes, by examining pairwise correlations of gene expression with SNP genotypes.

METHODS

The present study used data downloaded from the CommonMind Consortium site, consisting of SNP genotypes and RNAseq expression data from the dorsolateral prefrontal cortex, to examine whether the expression of individual genes or sets of genes correlated with PRS in 207 controls and 209 schizophrenia cases.

RESULTS

Although the PRS was significantly associated with phenotype, the correlations with genes and gene sets followed distributions expected by chance. Thus, this analysis failed to show that the PRS captures a cumulative effect of multiple variants impacting the expression of a small number of genes and it failed to focus attention on a small number of genes of biological relevance.

CONCLUSION

The multiple SNP associations observed in schizophrenia may result from other mechanisms, including effects mediated indirectly through environmental risk factors.

Prader-Willi locus Snord116 RNA processing requires an active endogenous allele and neuron-specific splicing by Rbfox3/NeuN

Prader-Willi syndrome (PWS), an imprinted neurodevelopmental disorder characterized by metabolic, sleep and neuropsychiatric features, is caused by the loss of paternal SNORD116, containing only non-coding RNAs (ncRNAs). The primary SNORD116 transcript is processed into small nucleolar RNAs (snoRNAs), which localize to nucleoli, and their spliced host gene 116HG, which is retained at its site of transcription. While functional complementation of the SNORD116 ncRNAs is a desirable goal for treating PWS, the mechanistic requirements of SNORD116 RNA processing are poorly understood. Here we developed and tested a novel transgenic mouse which ubiquitously expresses Snord116 on both a wild-type and a Snord116 paternal deletion (Snord116+/-) background. Interestingly, while the Snord116 transgene was ubiquitously expressed in multiple tissues, splicing of the transgene and production of snoRNAs was limited to brain tissues. Knockdown of Rbfox3, encoding neuron-specific splicing factor neuronal nuclei (NeuN) in Snord116+/--derived neurons, reduced splicing of the transgene in neurons. RNA fluorescence in situ hybridization for 116HG revealed a single significantly larger signal in transgenic mice, demonstrating colocalization of transgenic and endogenous 116HG RNAs. Similarly, significantly increased snoRNA levels were detected in transgenic neuronal nucleoli, indicating that transgenic Snord116 snoRNAs were effectively processed and localized. In contrast, neither transgenic 116HG nor snoRNAs were detectable in either non-neuronal tissues or Snord116+/- neurons. Together, these results demonstrate that exogenous expression and neuron-specific splicing of the Snord116 locus are insufficient to rescue the genetic deficiency of Snord116 paternal deletion. Elucidating the mechanisms regulating Snord116 processing and localization is essential to develop effective gene replacement therapies for PWS.

Expanded Insights Into Mechanisms of Gene Expression and Disease Related Disruptions

Definitive molecular diagnoses in disorders apparently due to genetic or genomic defects are still lacking in a significant number of investigated cases, despite use of studies designed to discover defects in the protein coding regions of the genome. Increasingly studies are being designed to search for defects in the non-protein coding genome, and for alterations in gene expression. Here we review new insights into genomic elements involved in control of gene expression, including methods to analyze chromatin that is accessible for transcription factor binding, enhancers, chromatin looping, transcription, RNA binding proteins, and alternative splicing. We review new studies on levels of genome organization, including the occurrence of transcriptional domains and their boundary elements. Information is presented on specific malformation syndromes that arise due to structural genomic changes that impact the non-protein coding genome and sometimes impact specific transcriptional domains. We also review convergence of genome-wide association with studies of gene expression, discoveries related to expression quantitative trait loci and splicing quantitative trait loci and the relevance of these to specific complex common diseases. Aspects of epigenetic mechanisms and clinical applications of analyses of methylation signatures are also discussed.

Rbfox1 Mediates Cell-type-Specific Splicing in Cortical Interneurons

Cortical interneurons display a remarkable diversity in their morphology, physiological properties, and connectivity. Elucidating the molecular determinants underlying this heterogeneity is essential for understanding interneuron development and function. We discovered that alternative splicing differentially regulates the integration of somatostatin- and parvalbumin-expressing interneurons into nascent cortical circuits through the cell-type-specific tailoring of mRNAs. Specifically, we identified a role for the activity-dependent splicing regulator Rbfox1 in the development of cortical interneuron-subtype-specific efferent connectivity. Our work demonstrates that Rbfox1 mediates largely non-overlapping alternative splicing programs within two distinct but related classes of interneurons.

Genome-wide association study identifies RBFOX1 locus influencing brain glucose metabolism

BACKGROUND

Fluorodeoxyglucose f18 positron emission tomography (18F-FDG PET) is regarded as the only functional neuroimaging biomarker for degeneration which can be used to increase the certainty of Alzheimer's disease (AD) pathophysiological process in research settings or as an optional clinical tool where available. Although a decline in FDG metabolism was confirmed in some regions known to be associated with AD, there was little known about the genetic association of FDG metabolism in AD cohorts. In this study, we present the first genome-wide association study (GWAS) analysis of brain FDG metabolism.

METHODS

A total of 222 individuals were included from the Alzheimer's Disease Neuroimaging Initiative 1 (ADNI-1) cohort. All subjects were restricted to non-Hispanic Caucasians and met all quality control (QC) criteria. Associations of 18F-FDG with the genetic variants were assessed using PLINK 1.07 under the additive genetic model. Genome-wide associations were visualized using a software program R 3.2.3.

RESULTS

One significant SNP rs12444565 in RNA-binding Fox1 (RBFOX1) was found to have a strong association with 18F-FDG (P=6.06×10-8). Rs235141, rs79037, rs12526331 and rs12529764 were identified as four suggestive loci associated with 18F-FDG.

CONCLUSIONS

Our study results suggest that a genome-wide significant SNP (rs12444565) in the RBFOX1, and four suggestive loci (rs235141, rs79037, rs12526331 and rs12529764) are associated with 18F-FDG.

Coregulation of alternative splicing by hnRNPM and ESRP1 during EMT

The epithelial-mesenchymal transition (EMT) is a fundamental developmental process that is abnormally activated in cancer metastasis. Dynamic changes in alternative splicing occur during EMT. ESRP1 and hnRNPM are splicing regulators that promote an epithelial splicing program and a mesenchymal splicing program, respectively. The functional relationships between these splicing factors in the genome scale remain elusive. Comparing alternative splicing targets of hnRNPM and ESRP1 revealed that they coregulate a set of cassette exon events, with the majority showing discordant splicing regulation. Discordant splicing events regulated by hnRNPM show a positive correlation with splicing during EMT; however, concordant events do not, indicating the role of hnRNPM in regulating alternative splicing during EMT is more complex than previously understood. Motif enrichment analysis near hnRNPM-ESRP1 coregulated exons identifies guanine-uridine rich motifs downstream from hnRNPM-repressed and ESRP1-enhanced exons, supporting a general model of competitive binding to these cis-elements to antagonize alternative splicing. The set of coregulated exons are enriched in genes associated with cell migration and cytoskeletal reorganization, which are pathways associated with EMT. Splicing levels of coregulated exons are associated with breast cancer patient survival and correlate with gene sets involved in EMT and breast cancer subtyping. This study identifies complex modes of interaction between hnRNPM and ESRP1 in regulation of splicing in disease-relevant contexts.

A weighted burden test using logistic regression for integrated analysis of sequence variants, copy number variants and polygenic risk score

Previously described methods of analysis allow variants in a gene to be weighted more highly according to rarity and/or predicted function and then for the variant contributions to be summed into a gene-wise risk score, which can be compared between cases and controls using a t-test. However, this does not allow incorporating covariates into the analysis. Schizophrenia is an example of an illness where there is evidence that different kinds of genetic variation can contribute to risk, including common variants contributing to a polygenic risk score (PRS), very rare copy number variants (CNVs) and sequence variants. A logistic regression approach has been implemented to compare the gene-wise risk scores between cases and controls, while incorporating as covariates population principal components, the PRS and the presence of pathogenic CNVs and sequence variants. A likelihood ratio test is performed, comparing the likelihoods of logistic regression models with and without this score. The method was applied to an ethnically heterogeneous exome-sequenced sample of 6000 controls and 5000 schizophrenia cases. In the raw analysis, the test statistic is inflated but inclusion of principal components satisfactorily controls for this. In this dataset, the inclusion of the PRS and effect from CNVs and sequence variants had only small effects. The set of genes which are FMRP targets showed some evidence for enrichment of rare, functional variants among cases (p = 0.0005). This approach can be applied to any disease in which different kinds of genetic and non-genetic risk factors make contributions to risk.

Identifying the genetic risk factors for treatment response to lurasidone by genome-wide association study: A meta-analysis of samples from three independent clinical trials

A genome-wide association study (GWAS) of response of schizophrenia patients to the atypical antipsychotic drug, lurasidone, based on two double-blind registration trials, identified SNPs from four classes of genes as predictors of efficacy, but none were genome wide significant (GWS). After inclusion of data from a third lurasidone trial, meta-analysis identified a GWS marker and other findings consistent with our first study. The primary end-point was change in Total Positive and Negative Syndrome Scale (PANSS) between baseline and last observation carried forward. rs4736253, a genetic locus near KCNK9, encoding the K_2P9.1 potassium channel, with a role in cognition and neurodevelopment, was the top marker in patients of European ancestry (EUR) (n = 264), reaching GWS (p = 4.78 × 10^-8). rs10180106 (p = 4.92 × 10^-7), located at an intron region of CTNNA2, a SCZ risk gene important for dendritic spine stabilization, was one of other best response markers for EUR patients. SNPs at STXBP5L (rs511841, p = 2.63 × 10^-7) were the top markers for patients of African ancestry (n = 158). The association between PTPRD, NRG1, and MAGI1 previously reported to be related to response to lurasidone in the first two trials, showed a trend of significant association in the third trial. None of these genetic loci showed significant associations with clinical response in the corresponding placebo groups (n = 107 for EUR; n = 58 for AFR). This meta-analysis yielded the first GWAS-based GWS biomarker for lurasidone response and additional support for the conclusion that genes related to synaptic biology and/or risk for SCZ are the strongest predictors of response to lurasidone in schizophrenia patients.

Mapping transcriptome-wide protein-RNA interactions to elucidate RNA regulatory programs

BACKGROUND

Our understanding of post-transcriptional gene regulation has increased exponentially with the development of robust methods to define protein-RNA interactions across the transcriptome. In this review, we highlight the evolution and successful applications of crosslinking and immunoprecipitation (CLIP) methods to interrogate protein-RNA interactions in a transcriptome-wide manner.

RESULTS

Here, we survey the vast array of in vitro and in vivo approaches used to identify protein-RNA interactions, including but not limited to electrophoretic mobility shift assays, systematic evolution of ligands by exponential enrichment (SELEX), and RIP-seq. We particularly emphasize the advancement of CLIP technologies, and detail protocol improvements and computational tools used to analyze the output data. Importantly, we discuss how profiling protein-RNA interactions can delineate biological functions including splicing regulation, alternative polyadenylation, cytoplasmic decay substrates, and miRNA targets.

CONCLUSIONS

In summary, this review summarizes the benefits of characterizing RNA-protein networks to further understand the regulation of gene expression and disease pathogenesis. Our review comments on how future CLIP technologies can be adapted to address outstanding questions related to many aspects of RNA metabolism and further advance our understanding of RNA biology.

Genome-wide association study and identification of chromosomal enhancer maps in multiple brain regions related to autism spectrum disorder

Autism spectrum disorder (ASD) is a complex developmental disorder with strong genetic components involved. Recent studies have demonstrated the importance of non-coding regulatory variants for complex diseases. To explore the roles of chromosomal enhancer regions in the pathogenesis of ASD, we conducted an integrative analysis of genome-wide association study (GWAS) and brain region related enhancer-gene networks for ASD. The GWAS data of ASD were driven from a published study, involving 7,387 ASD cases and 8,567 controls. The enhancer-gene networks of eight brain regions were used here. The GWAS of ASD was first merged respectively with the enhancer datasets of eight brain regions. Pathway enrichment analysis was then performed to detect ASD associated pathways based on the enhancer-related single nucleotide polymorphism (SNPs) of each brain region. We detected multiple genes with brain region specific or common association signals, such as PGM3 (P value = 1.93 × 10^-5 ) and RWDD2A (P value = 1.93 × 10^-5 ) for hippocampus middle, and ENPP4 (all P values <0.05), and ENPP5 (all P values <0.05) for seven brain regions. By comparing the pathway enrichment analysis results of various brain regions, several cross brain regions pathways were detected for ASD, such as REACTOME_POTASSIUM_CHANNELS (all P values <0.05) for six brain regions and KEGG_CELL_ADHESION_MOLECULES_CAMS (all P values <0.05) for seven brain regions. In addition, several pathways were also identified for specific brain regions, such as REACTOME_CD28_DEPENDENT_PI3K_AKT_SIGNALING (P value = 4.00 × 10^-3 ) for angular gyrus, REACTOME_SIGNALING_BY_CONSTITUTIVELY_ACTIVE_EGFR (P value = 2.22 × 10^-3 ) for anterior caudate, and KEGG_PRION_DISEASES (P value = 1.00 × 10^-4 ) for germinal matrix. Our results provide novel clues for understanding the genetic basis of ASD. Autism Research 2019, 12: 26-32. © 2018 International Society for Autism Research, Wiley Periodicals, Inc. LAY SUMMARY: ASD is a complex developmental disorder with strong genetic components, but the pathogenesis of ASD is still unclear. Using the latest GWAS data and enhancer map, we explored the brain region related biological pathways associated with ASD. Our results provide novel clues for revealing the functional relevance of enhancer variants with ASD and understanding the genetic basis of ASD.

Cell type-specific CLIP reveals that NOVA regulates cytoskeleton interactions in motoneurons

Background

Alternative RNA processing plays an essential role in shaping cell identity and connectivity in the central nervous system. This is believed to involve differential regulation of RNA processing in various cell types. However, in vivo study of cell type-specific post-transcriptional regulation has been a challenge. Here, we describe a sensitive and stringent method combining genetics and CLIP (crosslinking and immunoprecipitation) to globally identify regulatory interactions between NOVA and RNA in the mouse spinal cord motoneurons.

Results

We developed a means of undertaking motoneuron-specific CLIP to explore motoneuron-specific protein–RNA interactions relative to studies of the whole spinal cord in mouse. This allowed us to pinpoint differential RNA regulation specific to motoneurons, revealing a major role for NOVA in regulating cytoskeleton interactions in motoneurons. In particular, NOVA specifically promotes the palmitoylated isoform of the cytoskeleton protein Septin 8 in motoneurons, which enhances dendritic arborization.

Conclusions

Our study demonstrates that cell type-specific RNA regulation is important for fine tuning motoneuron physiology and highlights the value of defining RNA processing regulation at single cell type resolution.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1493-2) contains supplementary material, which is available to authorized users.

The Future of Cross-Linking and Immunoprecipitation (CLIP)

To understand the assembly and functional outcomes of protein-RNA regulation, it is crucial to precisely identify the positions of such interactions. Cross-linking and immunoprecipitation (CLIP) serves this purpose by exploiting covalent protein-RNA cross-linking and RNA fragmentation, along with a series of stringent purification and quality control steps to prepare complementary DNA (cDNA) libraries for sequencing. Here we describe the core steps of CLIP, its primary variations, and the approaches to data analysis. We present the application of CLIP to studies of specific cell types in genetically engineered mice and discuss the mechanistic and physiologic insights that have already been gained from studies using CLIP. We conclude by discussing the future opportunities for CLIP, including studies of human postmortem tissues from disease patients and controls, RNA epigenetic modifications, and RNA structure. These and other applications of CLIP will continue to unravel fundamental gene regulatory mechanisms while providing important biologic and clinically relevant insights.

LIN28 Selectively Modulates a Subclass of Let-7 MicroRNAs

LIN28 is a bipartite RNA-binding protein that post-transcriptionally inhibits the biogenesis of let-7 microRNAs to regulate development and influence disease states. However, the mechanisms of let-7 suppression remain poorly understood because LIN28 recognition depends on coordinated targeting by both the zinc knuckle domain (ZKD), which binds a GGAG-like element in the precursor, and the cold shock domain (CSD), whose binding sites have not been systematically characterized. By leveraging single-nucleotide-resolution mapping of LIN28 binding sites in vivo, we determined that the CSD recognizes a (U)GAU motif. This motif partitions the let-7 microRNAs into two subclasses, precursors with both CSD and ZKD binding sites (CSD+) and precursors with ZKD but no CSD binding sites (CSD-). LIN28 in vivo recognition-and subsequent 3' uridylation and degradation-of CSD+ precursors is more efficient, leading to their stronger suppression in LIN28-activated cells and cancers. Thus, CSD binding sites amplify the regulatory effects of LIN28.

Precise temporal regulation of alternative splicing during neural development

Alternative splicing (AS) is one crucial step of gene expression that must be tightly regulated during neurodevelopment. However, the precise timing of developmental splicing switches and the underlying regulatory mechanisms are poorly understood. Here we systematically analyze the temporal regulation of AS in a large number of transcriptome profiles of developing mouse cortices, in vivo purified neuronal subtypes, and neurons differentiated in vitro. Our analysis reveals early-switch and late-switch exons in genes with distinct functions, and these switches accurately define neuronal maturation stages. Integrative modeling suggests that these switches are under direct and combinatorial regulation by distinct sets of neuronal RNA-binding proteins including Nova, Rbfox, Mbnl, and Ptbp. Surprisingly, various neuronal subtypes in the sensory systems lack Nova and/or Rbfox expression. These neurons retain the "immature" splicing program in early-switch exons, affecting numerous synaptic genes. These results provide new insights into the organization and regulation of the neurodevelopmental transcriptome.

Genetic identification of brain cell types underlying schizophrenia

With few exceptions, the marked advances in knowledge about the genetic basis of schizophrenia have not converged on findings that can be confidently used for precise experimental modeling. By applying knowledge of the cellular taxonomy of the brain from single-cell RNA sequencing, we evaluated whether the genomic loci implicated in schizophrenia map onto specific brain cell types. We found that the common-variant genomic results consistently mapped to pyramidal cells, medium spiny neurons (MSNs) and certain interneurons, but far less consistently to embryonic, progenitor or glial cells. These enrichments were due to sets of genes that were specifically expressed in each of these cell types. We also found that many of the diverse gene sets previously associated with schizophrenia (genes involved in synaptic function, those encoding mRNAs that interact with FMRP, antipsychotic targets, etc.) generally implicated the same brain cell types. Our results suggest a parsimonious explanation: the common-variant genetic results for schizophrenia point at a limited set of neurons, and the gene sets point to the same cells. The genetic risk associated with MSNs did not overlap with that of glutamatergic pyramidal cells and interneurons, suggesting that different cell types have biologically distinct roles in schizophrenia.

ZFP36 RNA-binding proteins restrain T cell activation and anti-viral immunity

Dynamic post-transcriptional control of RNA expression by RNA-binding proteins (RBPs) is critical during immune response. ZFP36 RBPs are prominent inflammatory regulators linked to autoimmunity and cancer, but functions in adaptive immunity are less clear. We used HITS-CLIP to define ZFP36 targets in mouse T cells, revealing unanticipated actions in regulating T-cell activation, proliferation, and effector functions. Transcriptome and ribosome profiling showed that ZFP36 represses mRNA target abundance and translation, notably through novel AU-rich sites in coding sequence. Functional studies revealed that ZFP36 regulates early T-cell activation kinetics cell autonomously, by attenuating activation marker expression, limiting T cell expansion, and promoting apoptosis. Strikingly, loss of ZFP36 in vivo accelerated T cell responses to acute viral infection and enhanced anti-viral immunity. These findings uncover a critical role for ZFP36 RBPs in restraining T cell expansion and effector functions, and suggest ZFP36 inhibition as a strategy to enhance immune-based therapies.

The immune system must quickly respond to anything that may cause disease – from cancerous cells to viruses. For instance, a type of white blood cell called a T cell patrols the body, looking for potential threats. If a T cell identifies such a threat, it “activates” and undergoes various changes so that it can help to eliminate the problem.

One way that T cells change is by switching on different genes to make specific proteins. The information in the genes is first used as a template to produce a molecule called a messenger RNA (mRNA), which is then translated to build proteins. So-called RNA-binding proteins help control events before, during and after the translation stage in the process. Previous studies have shown that one particular RNA-binding protein, called ZFP36, controls the translation of proteins that are important for how the immune system recognizes the body’s own tissue and deals with cancer cells. However, it was less clear if it also helped T cells to activate and defeat viruses.

Now, using cutting-edge technology, Moore et al. have identified thousands of new mRNAs controlled by ZFP36 in mice, many of which did indeed make proteins that help T cells activate and spread throughout the body. Further experiments showed that mice that lack ZFP36 in the T cells were much quicker at responding to viruses than other mice. This suggests that ZFP36 actually restrains T cells and slows down the body’s immune system.

Knowing more about how T cells work could lead to new treatments for diseases; it may, for example, allow scientists to engineer T cells to better attack cancer cells, However, other studies have shown that mice without ZFP36 often go on to develop autoimmune diseases, which result from the immune system attacking healthy cells by mistake. As such, it seems that there is a fine line between improving the body’s immune system and increasing the risk of autoimmune diseases, and that RNA-binding proteins play an important role in managing this delicate balance.

rbFOX1/MBNL1 competition for CCUG RNA repeats binding contributes to myotonic dystrophy type 1/type 2 differences

Myotonic dystrophy type 1 and type 2 (DM1, DM2) are caused by expansions of CTG and CCTG repeats, respectively. RNAs containing expanded CUG or CCUG repeats interfere with the metabolism of other RNAs through titration of the Muscleblind-like (MBNL) RNA binding proteins. DM2 follows a more favorable clinical course than DM1, suggesting that specific modifiers may modulate DM severity. Here, we report that the rbFOX1 RNA binding protein binds to expanded CCUG RNA repeats, but not to expanded CUG RNA repeats. Interestingly, rbFOX1 competes with MBNL1 for binding to CCUG expanded repeats and overexpression of rbFOX1 partly releases MBNL1 from sequestration within CCUG RNA foci in DM2 muscle cells. Furthermore, expression of rbFOX1 corrects alternative splicing alterations and rescues muscle atrophy, climbing and flying defects caused by expression of expanded CCUG repeats in a Drosophila model of DM2.

CLIP-related methodologies and their application to retrovirology

Virtually every step of HIV-1 replication and numerous cellular antiviral defense mechanisms are regulated by the binding of a viral or cellular RNA-binding protein (RBP) to distinct sequence or structural elements on HIV-1 RNAs. Until recently, these protein-RNA interactions were studied largely by in vitro binding assays complemented with genetics approaches. However, these methods are highly limited in the identification of the relevant targets of RBPs in physiologically relevant settings. Development of crosslinking-immunoprecipitation sequencing (CLIP) methodology has revolutionized the analysis of protein-nucleic acid complexes. CLIP combines immunoprecipitation of covalently crosslinked protein-RNA complexes with high-throughput sequencing, providing a global account of RNA sequences bound by a RBP of interest in cells (or virions) at near-nucleotide resolution. Numerous variants of the CLIP protocol have recently been developed, some with major improvements over the original. Herein, we briefly review these methodologies and give examples of how CLIP has been successfully applied to retrovirology research.

Control of CNS functions by RNA-binding proteins in neurological diseases

Purpose of Review

This review summarizes recent studies on the molecular mechanisms of RNA binding proteins (RBPs) that control neurological functions and pathogenesis in various neurodevelopmental and neurodegenerative diseases, including autism spectrum disorders, schizophrenia, Alzheimer's disease, amyotrophic lateral sclerosis, frontotemporal dementia, and spinocerebellar ataxia.

Recent Findings

RBPs are critical players in gene expression that regulate every step of posttranscriptional modifications. Recent genome-wide approaches revealed that many proteins associate with RNA, but do not contain any known RNA binding motifs. Additionally, many causal and risk genes of neurodevelopmental and neurodegenerative diseases are RBPs. Development of high-throughput sequencing methods has mapped out the fingerprints of RBPs on transcripts and provides unprecedented potential to discover new mechanisms of neurological diseases. Insights into how RBPs modulate neural development are important for designing effective therapies for numerous neurodevelopmental and neurodegenerative diseases.

Summary

RBPs have diverse mechanisms for modulating RNA processing and, thereby, controlling neurogenesis. Understanding the role of disease-associated RBPs in neurogenesis is vital for developing novel treatments for neurological diseases.

Distinct and Modular Organization of Protein Interacting Sites in Long Non-coding RNAs

Background: Long non-coding RNAs (lncRNAs), are being reported to be extensively involved in diverse regulatory roles and have exhibited numerous disease associations. LncRNAs modulate their function through interaction with other biomolecules in the cell including DNA, RNA, and proteins. The availability of genome-scale experimental datasets of RNA binding proteins (RBP) motivated us to understand the role of lncRNAs in terms of its interactions with these proteins. In the current report, we demonstrate a comprehensive study of interactions between RBP and lncRNAs at a transcriptome scale through extensive analysis of the crosslinking and immunoprecipitation (CLIP) experimental datasets available for 70 RNA binding proteins.

Results: Our analysis suggests that density of interaction sites for these proteins was significantly higher for specific sub-classes of lncRNAs when compared to protein-coding transcripts. We also observe a positional preference of these RBPs across lncRNA and protein coding transcripts in addition to a significant co-occurrence of RBPs having similar functions, suggesting a modular organization of these elements across lncRNAs.

Conclusion: The significant enrichment of RBP sites across some lncRNA classes is suggestive that these interactions might be important in understanding the functional role of lncRNA. We observed a significant enrichment of RBPs which are involved in functional roles such as silencing, splicing, mRNA processing, and transport, indicating the potential participation of lncRNAs in such processes.

Plant Ribonomics: Proteins in Search of RNA Partners

Research into the regulation of gene expression underwent a shift from focusing on DNA-binding proteins as key transcriptional regulators to RNA-binding proteins (RBPs) that come into play once transcription has been initiated. RBPs orchestrate all RNA-processing steps in the cell. To obtain a global view of in vivo targets, the RNA complement associated with particular RBPs is determined via immunoprecipitation of the RBP and subsequent identification of bound RNAs via RNA-seq. Here, we describe technical advances in identifying RBP in vivo targets and their binding motifs. We provide an up-to-date view of targets of nucleocytoplasmic RBPs collected in arabidopsis. We also discuss current experimental limitations and provide an outlook on how the approaches may advance our understanding of post-transcriptional networks.

Weighted Burden Analysis of Exome-Sequenced Case-Control Sample Implicates Synaptic Genes in Schizophrenia Aetiology

A previous study of exome-sequenced schizophrenia cases and controls reported an excess of singleton, gene-disruptive variants among cases, concentrated in particular gene sets. The dataset included a number of subjects with a substantial Finnish contribution to ancestry. We have reanalysed the same dataset after removal of these subjects and we have also included non-singleton variants of all types using a weighted burden test which assigns higher weights to variants predicted to have a greater effect on protein function. We investigated the same 31 gene sets as previously and also 1454 GO gene sets. The reduced dataset consisted of 4225 cases and 5834 controls. No individual variants or genes were significantly enriched in cases but 13 out of the 31 gene sets were significant after Bonferroni correction and the “FMRP targets” set produced a signed log p value (SLP) of 7.1. The gene within this set with the highest SLP, equal to 3.4, was FYN, which codes for a tyrosine kinase which phosphorylates glutamate metabotropic receptors and ionotropic NMDA receptors, thus modulating their trafficking, subcellular distribution and function. In the most recent GWAS of schizophrenia it was identified as a “prioritized candidate gene”. Two of the subunits of the NMDA receptor which are substrates of FYN are coded for by GRIN1 (SLP = 1.7) and GRIN2B (SLP = 2.1). Of note, for some sets there was a substantial enrichment of non-singleton variants. Of 1454 GO gene sets, three were significant after Bonferroni correction. Identifying specific genes and variants will depend on genotyping them in larger samples and/or demonstrating that they cosegregate with illness within pedigrees.

Kaposi's Sarcoma-Associated Herpesvirus K8 Is an RNA Binding Protein That Regulates Viral DNA Replication in Coordination with a Noncoding RNA

Kaposi's sarcoma-associated herpesvirus (KSHV) lytic replication and constant primary infection of fresh cells are crucial for viral tumorigenicity. The virus-encoded bZIP family protein K8 plays an important role in viral DNA replication in both viral reactivation and de novo infection. The mechanism underlying the functional role of K8 in the viral life cycle is elusive. Here, we report that K8 is an RNA binding protein that also associates with many other proteins, including other RNA binding proteins. Many protein-protein interactions involving K8 are mediated by RNA. Using a UV cross-linking and immunoprecipitation (CLIP) procedure combined with high-throughput sequencing, RNAs that are associated with K8 in BCBL-1 cells were identified, including both viral (PAN, T1.4, T0.7, etc.) and cellular (MALAT-1, MRP, 7SK, etc.) RNAs. An RNA binding motif in K8 was defined, and mutation of the motif abolished the ability of K8 to bind to many noncoding RNAs, as well as viral DNA replication during de novo infection, suggesting that the K8 functions in viral replication are carried out through RNA association. The functions of K8 and associated T1.4 RNA were investigated in detail, and the results showed that T1.4 mediates the binding of K8 to ori-Lyt DNA. The T1.4-K8 complex physically bound to KSHV ori-Lyt DNA and recruited other proteins and cofactors to assemble a replication complex. Depletion of T1.4 abolished DNA replication in primary infection. These findings provide mechanistic insights into the role of K8 in coordination with T1.4 RNA in regulating KSHV DNA replication during de novo infection.IMPORTANCE Genomewide analyses of the mammalian transcriptome revealed that a large proportion of sequence previously annotated as noncoding regions is actually transcribed and gives rise to stable RNAs. The emergence of a large number of noncoding RNAs suggests that functional RNA-protein complexes, e.g., ribosomes or spliceosomes, are not ancient relics of the last ribo-organism but would be well adapted to a regulatory role in biology. K8 has been puzzling because of its unique characteristics, such as multiple regulatory roles in gene expression and DNA replication without DNA binding capability. This study reveals the mechanism underlying its regulatory role by demonstrating that K8 is an RNA binding protein that binds to DNA and initiates DNA replication in coordination with a noncoding RNA. It is suggested that many K8 functions, if not all, are carried out through its associated RNAs.