Identification of miR-128 Target mRNAs That Are Expressed in B Cells Using a Modified Dual Luciferase Vector

MicroRNAs (miRNAs) are 21-25 nucleotide long non-coding ribonucleic acids that modulate gene expression by degrading transcripts or inhibiting translation. The miRNA miR-128, originally thought to be brain-specific, was later also found in immune cells. To identify a valuable immune cell model system to modulate endogenous miR-128 amounts and to validate predicted miR-128 target mRNAs in B cells, we first investigated miR-128 expression using Northern blot analysis in several cell lines representing different stages of B cell development. The results showed that only primary brain cells showed significant levels of mature miR-128. To study the function of miR-128 in immune cells, we modified dual luciferase vectors to allow easy transfer of 3' UTR fragments with predicted miR-128 binding sites from widely used single to dual luciferase vectors. Comparison of in silico predicted miR-128-regulated mRNAs in single and dual luciferase constructs yielded similar results, validating the dual luciferase vector for miRNA target analysis. Furthermore, we confirmed miR-128-regulated mRNAs identified in silico and in vivo using the Ago HITS-CLIP technique and known to be expressed in B cells using the dual luciferase assay. In conclusion, this study provides new insights into the expression and function of miR-128 by validating novel target mRNAs expressed in B cells and identifying additional pathways likely controlled by this miRNA in the immune system.

Local changes in viral RNA sequence drive global changes in influenza nucleoprotein binding

The genome of influenza A viruses (IAV) consists of eight negative-sense RNA segments that are coated by viral nucleoprotein (NP). Until recently, it was assumed that NP binds viral genomic RNA (vRNA) uniformly along the entire segment. However, genome-wide studies have revised the original model in that NP instead binds preferentially to certain regions of vRNA, while others are depleted for NP binding. Even strains with high sequence similarity exhibit distinct NP-binding profiles. Thus, it remains unknown how NP-binding specificity to vRNA is established. Here we introduced nucleotide changes to vRNA to examine whether primary sequence can affect NP binding. Our findings demonstrate that NP binding is indeed susceptible to sequence alterations, as NP peaks can be lost or appear de novo at mutated sites. Unexpectedly, nucleotide changes not only affect NP binding locally at the site of mutation, but also impact NP binding at distal regions that have not been modified. Taken together, our results suggest that NP binding is not regulated by primary sequence alone, but that a network formed by multiple segments governs the deposition of NP on vRNA.

HITS-CLIP analysis of human ALKBH8 reveals interactions with fully processed substrate tRNAs and with specific noncoding RNAs

Transfer RNAs acquire a large plethora of chemical modifications. Among those, modifications of the anticodon loop play important roles in translational fidelity and tRNA stability. Four human wobble U-containing tRNAs obtain 5-methoxycarbonylmethyluridine (mcm5U34) or 5-methoxycarbonylmethyl-2-thiouridine (mcm5s2U34), which play a role in decoding. This mark involves a cascade of enzymatic activities. The last step is mediated by alkylation repair homolog 8 (ALKBH8). In this study, we performed a transcriptome-wide analysis of the repertoire of ALKBH8 RNA targets. Using a combination of HITS-CLIP and RIP-seq analyses, we uncover ALKBH8-bound RNAs. We show that ALKBH8 targets fully processed and CCA modified tRNAs. Our analyses uncovered the previously known set of wobble U-containing tRNAs. In addition, both our approaches revealed ALKBH8 binding to several other types of noncoding RNAs, in particular C/D box snoRNAs.

Next-generation sequencing: A new avenue to understand viral RNA-protein interactions

The genomes of RNA viruses present an astonishing source of both sequence and structural diversity. From intracellular viral RNA-host interfaces to interactions between the RNA genome and structural proteins in virus particles themselves, almost the entire viral lifecycle is accompanied by a myriad of RNA-protein interactions that are required to fulfill their replicative potential. It is therefore important to characterize such rich and dynamic collections of viral RNA-protein interactions to understand virus evolution and their adaptation to their hosts and environment. Recent advances in next-generation sequencing technologies have allowed the characterization of viral RNA-protein interactions, including both transient and conserved interactions, where molecular and structural approaches have fallen short. In this review, we will provide a methodological overview of the high-throughput techniques used to study viral RNA-protein interactions, their biochemical mechanisms, and how they evolved from classical methods as well as one another. We will discuss how different techniques have fueled virus research to characterize how viral RNA and proteins interact, both locally and on a global scale. Finally, we will present examples on how these techniques influence the studies of clinically important pathogens such as HIV-1 and SARS-CoV-2.

Solid-Support Directional (SSD) RNA-Seq as a Companion Method to CLIP-Seq

CLIP-Seq (Deep Sequencing after in vivo Crosslinking and Immunoprecipitation, HITS-CLIP) has emerged as a key method for the study of RNA-binding proteins (RBPs), as it can scrutinize the RNAs bound by an RBP in vivo, with minimum manipulation of biological samples. CLIP-Seq is best used to reveal changes of the RNA cargo of an RBP and differences on binding patterns of the bound RNAs in living cells in different genetic backgrounds or after experimental treatment, rather than simply identifying RNA species. It is therefore crucial that a reference of the steady state levels of the RNAs present in the samples used for the CLIP-Seq experiment is included in the bioinformatic analysis. A simple directional RNA-Seq method was developed that uses the same oligonucleotides and the same PCR amplification steps as our CLIP-Seq method, which therefore can be analyzed using the same bioinformatic pipeline as the CLIP-Seq data. This greatly simplifies and streamlines the analysis process, and at the same time reduces the chances of protocol-specific artifacts and biases interfering with data interpretation. Some considerations on ways to integrate CLIP-Seq and RNA-Seq analyses are also provided herein.

Star-PAP RNA Binding Landscape Reveals Novel Role of Star-PAP in mRNA Metabolism That Requires RBM10-RNA Association

Star-PAP is a non-canonical poly(A) polymerase that selects mRNA targets for polyadenylation. Yet, genome-wide direct Star-PAP targets or the mechanism of specific mRNA recognition is still vague. Here, we employ HITS-CLIP to map the cellular Star-PAP binding landscape and the mechanism of global Star-PAP mRNA association. We show a transcriptome-wide association of Star-PAP that is diminished on Star-PAP depletion. Consistent with its role in the 3'-UTR processing, we observed a high association of Star-PAP at the 3'-UTR region. Strikingly, there is an enrichment of Star-PAP at the coding region exons (CDS) in 42% of target mRNAs. We demonstrate that Star-PAP binding de-stabilises these mRNAs indicating a new role of Star-PAP in mRNA metabolism. Comparison with earlier microarray data reveals that while UTR-associated transcripts are down-regulated, CDS-associated mRNAs are largely up-regulated on Star-PAP depletion. Strikingly, the knockdown of a Star-PAP coregulator RBM10 resulted in a global loss of Star-PAP association on target mRNAs. Consistently, RBM10 depletion compromises 3'-end processing of a set of Star-PAP target mRNAs, while regulating stability/turnover of a different set of mRNAs. Our results establish a global profile of Star-PAP mRNA association and a novel role of Star-PAP in the mRNA metabolism that requires RBM10-mRNA association in the cell.

Functional microRNA targetome undergoes degeneration-induced shift in the retina

BACKGROUND

MicroRNA (miRNA) play a significant role in the pathogenesis of complex neurodegenerative diseases including age-related macular degeneration (AMD), acting as post-transcriptional gene suppressors through their association with argonaute 2 (AGO2) - a key member of the RNA Induced Silencing Complex (RISC). Identifying the retinal miRNA/mRNA interactions in health and disease will provide important insight into the key pathways miRNA regulate in disease pathogenesis and may lead to potential therapeutic targets to mediate retinal degeneration.

METHODS

To identify the active miRnome targetome interactions in the healthy and degenerating retina, AGO2 HITS-CLIP was performed using a rodent model of photoreceptor degeneration. Analysis of publicly available single-cell RNA sequencing (scRNAseq) data was performed to identify the cellular location of AGO2 and key members of the microRNA targetome in the retina. AGO2 findings were verified by in situ hybridization (RNA) and immunohistochemistry (protein).

RESULTS

Analysis revealed a similar miRnome between healthy and damaged retinas, however, a shift in the active targetome was observed with an enrichment of miRNA involvement in inflammatory pathways. This shift was further demonstrated by a change in the seed binding regions of miR-124-3p, the most abundant retinal AGO2-bound miRNA, and has known roles in regulating retinal inflammation. Additionally, photoreceptor cluster miR-183/96/182 were all among the most highly abundant miRNA bound to AGO2. Following damage, AGO2 expression was localized to the inner retinal layers and more in the OLM than in healthy retinas, indicating a locational miRNA response to retinal damage.

CONCLUSIONS

This study provides important insight into the alteration of miRNA regulatory activity that occurs as a response to retinal degeneration and explores the miRNA-mRNA targetome as a consequence of retinal degenerations. Further characterisation of these miRNA/mRNA interactions in the context of the degenerating retina may provide an important insight into the active role these miRNA may play in diseases such as AMD.

AGO HITS-CLIP reveals distinct miRNA regulation of white and brown adipose tissue identity

MicroRNAs (miRNAs) are short, noncoding RNAs that associate with Argonaute (AGO) to influence mRNA stability and translation, thereby regulating cellular determination and phenotype. While several individual miRNAs have been shown to control adipocyte function, including energy storage in white fat and energy dissipation in brown fat, a comprehensive analysis of miRNA activity in these tissues has not been performed. We used high-throughput sequencing of RNA isolated by cross-linking immunoprecipitation (HITS-CLIP) to comprehensively characterize the network of high-confidence, in vivo mRNA:miRNA interactions across white and brown fat, revealing >20,000 unique AGO binding sites. When coupled with miRNA and mRNA sequencing, we found an inverse correlation between depot-enriched miRNAs and their targets. To illustrate the functionality of our HITS-CLIP data set in identifying specific miRNA:mRNA interactions, we show that miR-29 is a novel regulator of leptin, an adipocyte-derived hormone that coordinates food intake and energy homeostasis. Two independent miR-29 binding sites in the leptin 3' UTR were validated using luciferase assays, and miR-29 gain and loss of function modulated leptin mRNA and protein secretion in primary adipocytes. This work represents the only experimentally generated miRNA targetome in adipose tissue and identifies multiple regulatory pathways that may specify the unique identities of white and brown fat.

Argonaute-CLIP delineates versatile, functional RNAi networks in Aedes aegypti, a major vector of human viruses

Argonaute (AGO) proteins bind small RNAs to silence complementary RNA transcripts, and they are central to RNA interference (RNAi). RNAi is critical for regulation of gene expression and antiviral defense in Aedes aegypti mosquitoes, which transmit Zika, chikungunya, dengue, and yellow fever viruses. In mosquitoes, AGO1 mediates miRNA interactions, while AGO2 mediates siRNA interactions. We applied AGO-crosslinking immunoprecipitation (AGO-CLIP) for both AGO1 and AGO2, and we developed a universal software package for CLIP analysis (CLIPflexR), identifying 230 small RNAs and 5,447 small RNA targets that comprise a comprehensive RNAi network map in mosquitoes. RNAi network maps predicted expression levels of small RNA targets in specific tissues. Additionally, this resource identified unexpected, context-dependent AGO2 target preferences, including endogenous viral elements and 3'UTRs. Finally, contrary to current thinking, mosquito AGO2 repressed imperfect targets. These findings expand our understanding of small RNA networks and have broad implications for the study of antiviral RNAi.

HITS-CLIP reveals sex-differential RNA binding and alterative splicing regulation of SRm160 in Drosophila

Serine/arginine (SR)-rich proteins are critical for the regulation of alternative splicing (AS), which generates multiple mRNA isoforms from one gene and provides protein diversity for cell differentiation and tissue development. Genetic evidence suggests that Drosophila genital-specific overexpression of SR-related nuclear matrix protein of 160 kDa (SRm160), an SR protein with a PWI RNA-binding motif, causes defective development only in male flies and results in abnormal male genital structures and abnormal testis. However, the molecular characterization of SRm160 is limited. Using the high-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP) method in two sex-specific embryonic cell lines, S2 from the male and Kc from the female, we first identified the genome-wide RNA-binding characteristics of SRm160, which preferred binding to the exonic tri-nucleotide repeats GCA and AAC. We then validated this binding through both in vitro gel-shift assay and in vivo splicing of minigenes and found that SRm160 level affects AS of many transcripts. Furthermore, we identified 492 differential binding sites (DBS) of SRm160 varying between the two sex-specific cell lines. Among these DBS-containing genes, splicing factors were highly enriched, including transformer, a key regulator in the sex determination cascade. Analyses of fly mutants demonstrated that the SRm160 level affects AS isoforms of transformer. These findings shed crucial light on SRm160's RNA-binding specificity and regulation of AS in Drosophila sex determination and development.

LIN28A binds to meiotic gene transcripts and modulates their translation in male germ cells

The RNA-binding protein LIN28A is required for maintaining tissue homeostasis, including in the reproductive system, but the underlying mechanisms on how LIN28A regulates germline progenitors remain unclear. Here, we dissected LIN28A-binding targets using high-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation (HITS-CLIP) in the mouse testes. LIN28A preferentially binds to mRNA coding sequence (CDS) or 3'UTR regions at sites enriched with GGAG(A) sequences. Further investigation of Lin28a-null mouse testes indicated that meiosis-associated mRNAs bound by LIN28A were differentially expressed. Next, ribosome profiling revealed that the mRNA levels of these targets were significantly reduced in the polysome fractions, and their protein expression levels decreased, in Lin28a-null mouse testes, even when meiotic arrest in the null mouse testes was not apparent. Collectively, these findings provide a set of LIN28A-regulated target mRNAs, and show that LIN28A binding might be a mechanism through which LIN28A acts to regulate undifferentiated spermatogonia fates and male fertility in mammals.

An RNA Switch of a Large Exon of Ninein Is Regulated by the Neural Stem Cell Specific-RNA Binding Protein, Qki5

A set of tissue-specific splicing factors are thought to govern alternative splicing events during neural progenitor cell (NPC)-to-neuron transition by regulating neuron-specific exons. Here, we propose one such factor, RNA-binding protein Quaking 5 (Qki5), which is specifically expressed in the early embryonic neural stem cells. We performed mRNA-SEQ (Sequence) analysis using mRNAs obtained by developing cerebral cortices in Qk (Quaking) conditional knockout (cKO) mice. As expected, we found a large number of alternative splicing changes between control and conditional knockouts relative to changes in transcript levels. DAVID (The Database for Annotation, Visualization and Integrated Discovery) and Metascape analyses suggested that the affected spliced genes are involved in axon development and microtubule-based processes. Among these, the mRNA coding for the Ninein protein is listed as one of Qki protein-dependent alternative splicing targets. Interestingly, this exon encodes a very long polypeptide (2121 nt), and has been previously defined as a dynamic RNA switch during the NPC-to-neuron transition. Additionally, we validated that the regulation of this large exon is consistent with the Qki5-dependent alternative exon inclusion mode suggested by our previous Qki5 HITS-CLIP (high throughput sequencing-cross linking immunoprecipitation) analysis. Taken together, these data suggest that Qki5 is an important factor for alternative splicing in the NPC-to-neuron transition.

Investigating miRNA-lncRNA Interactions: Computational Tools and Resources

In the last two decades noncoding RNAs have been the recipients of increasing scientific interest. In particular, miRNAs, short (~22 nts) noncoding transcripts, have been thoroughly investigated since their essential role in posttranscriptional gene expression regulation had been established in the early 2000s. With the advent and the advancements of high-throughput sequencing technologies in recent years, long noncoding RNAs have also started to emerge as important actors in cellular functions and processes. Such transcripts, on average longer than 200 nt, whose functions have yet to be fully characterized, have recently been identified as regulatory elements of the RNAi pathway, harboring several miRNA response elements, uncovering the phenomena of competing endogenous RNAs (ceRNAs), or "sponge RNAs." The present chapter aims to provide a brief update on the actual biomedical relevance of ceRNAs, together with a summary of resources, tools, and practical examples of their application to aid researchers in the discovery and further elucidation of lncRNA-miRNA interactions.

Contemporary Ribonomics Methods for Viral microRNA Target Analysis

Numerous cellular processes are regulated by microRNAs (miRNAs), both cellular and viral. Elucidating the targets of miRNAs has become an active area of research. An important method in this field is cross-linking and immunoprecipitation (CLIP), where cultured cells or tissues are UV-irradiated to cross-link protein and nucleic acid, the RNA binding protein of interest is immunoprecipitated, and the RNAs pulled down with the protein are isolated, reverse-transcribed, and analyzed by sequencing. CLIP using antibody against Argonaute (Ago), which binds to both miRNA and mRNA as they interact in RISC, has allowed researchers to uncover a large number of miRNA targets. Coupled with high-throughput sequencing, CLIP has been useful for revealing miRNA targetomes for the γ-herpesviruses Kaposi's sarcoma-associated herpesvirus (KSHV) and Epstein-Barr virus (EBV). Variants on the CLIP protocol are described, with the benefits and drawbacks of each. In particular, the most recent methods involving RNA⁻RNA ligation to join the miRNA and its RNA target have aided in target identification. Lastly, data supporting biologically meaningful interactions between miRNAs and long non-coding RNAs (lncRNAs) are reviewed. In summary, ribonomics-based miRNA targetome analysis has expanded our understanding of miRNA targeting and has provided a rich resource for EBV and KSHV research with respect to pathogenesis and tumorigenesis.

omniCLIP: probabilistic identification of protein-RNA interactions from CLIP-seq data

CLIP-seq methods allow the generation of genome-wide maps of RNA binding protein – RNA interaction sites. However, due to differences between different CLIP-seq assays, existing computational approaches to analyze the data can only be applied to a subset of assays. Here, we present a probabilistic model called omniCLIP that can detect regulatory elements in RNAs from data of all CLIP-seq assays. omniCLIP jointly models data across replicates and can integrate background information. Therefore, omniCLIP greatly simplifies the data analysis, increases the reliability of results and paves the way for integrative studies based on data from different assays.

Non-Uniform and Non-Random Binding of Nucleoprotein to Influenza A and B Viral RNA

The genomes of influenza A and B viruses have eight, single-stranded RNA segments that exist in the form of a viral ribonucleoprotein complex in association with nucleoprotein (NP) and an RNA-dependent RNA polymerase complex. We previously used high-throughput RNA sequencing coupled with crosslinking immunoprecipitation (HITS-CLIP) to examine where NP binds to the viral RNA (vRNA) and demonstrated for two H1N1 strains that NP binds vRNA in a non-uniform, non-random manner. In this study, we expand on those initial observations and describe the NP-vRNA binding profile for a seasonal H3N2 and influenza B virus. We show that, similar to H1N1 strains, NP binds vRNA in a non-uniform and non-random manner. Each viral gene segment has a unique NP binding profile with areas that are enriched for NP association as well as free of NP-binding. Interestingly, NP-vRNA binding profiles have some conservation between influenza A viruses, H1N1 and H3N2, but no correlation was observed between influenza A and B viruses. Our study demonstrates the conserved nature of non-uniform NP binding within influenza viruses. Mapping of the NP-bound vRNA segments provides information on the flexible NP regions that may be involved in facilitating assembly.

Identification and removal of sequencing artifacts produced by mispriming during reverse transcription in multiple RNA-seq technologies

The quality of RNA sequencing data relies on specific priming by the primer used for reverse transcription (RT-primer). Nonspecific annealing of the RT-primer to the RNA template can generate reads with incorrect cDNA ends and can cause misinterpretation of data (RT mispriming). This kind of artifact in RNA-seq based technologies is underappreciated and currently no adequate tools exist to computationally remove them from published data sets. We show that mispriming can occur with as little as two bases of complementarity at the 3' end of the primer followed by intermittent regions of complementarity. We also provide a computational pipeline that identifies cDNA reads produced from RT mispriming, allowing users to filter them out from any aligned data set. Using this analysis pipeline, we identify thousands of mispriming events in a dozen published data sets from diverse technologies including short RNA-seq, total/mRNA-seq, HITS-CLIP, and GRO-seq. We further show how RT mispriming can lead to misinterpretation of data. In addition to providing a solution to computationally remove RT-misprimed reads, we also propose an experimental solution to completely avoid RT-mispriming by performing RNA-seq using thermostable group II intron derived reverse transcriptase (TGIRT-seq).

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.

cCLIP-Seq: Retrieval of Chimeric Reads from HITS-CLIP (CLIP-Seq) Libraries

HITS-CLIP (High-Throughput Sequencing after in vivo Crosslinking and Immunoprecipitation, CLIP-Seq) libraries contain fragments of the RNA sequences bound in vivo by an RNA binding protein (RBP). Such fragments, especially if they represent RNA duplexes bound in vivo by the RBP, can occasionally be ligated together to form chimeric CLIP tags. Chimeric CLIP tags from Argonaute CLIP libraries can provide the exact base pairing profiles of small RNAs with their target RNA sequences, thus solving a critical problem in the field of post-transcriptional regulation. We recently reported an analysis of chimeric reads from the Drosophila Piwi protein Aubergine, which revealed a novel mechanism for mRNA entrapment within germ RNP granules. We term this novel approach chimeric CLIP (cCLIP) and present here the main steps that a researcher can take after the acquisition of the deep sequencing data, for the identification of candidate chimeric reads in Piwi CLIP libraries. Extending the scope beyond small-RNA binding proteins, we believe that cCLIP can be utilized to elucidate the in vivo functions of RNA-binding proteins in general, and especially those that modulate RNA secondary structures. We, therefore, also describe aspects of the generalized chimeric read identification problem, which can find use in the analysis of the CLIP libraries of any RNA-binding protein.

An RNA-binding protein, Qki5, regulates embryonic neural stem cells through pre-mRNA processing in cell adhesion signaling

Cell type-specific transcriptomes are enabled by the action of multiple regulators, which are frequently expressed within restricted tissue regions. In the present study, we identify one such regulator, Quaking 5 (Qki5), as an RNA-binding protein (RNABP) that is expressed in early embryonic neural stem cells and subsequently down-regulated during neurogenesis. mRNA sequencing analysis in neural stem cell culture indicates that Qki proteins play supporting roles in the neural stem cell transcriptome and various forms of mRNA processing that may result from regionally restricted expression and subcellular localization. Also, our in utero electroporation gain-of-function study suggests that the nuclear-type Qki isoform Qki5 supports the neural stem cell state. We next performed in vivo transcriptome-wide protein-RNA interaction mapping to search for direct targets of Qki5 and elucidate how Qki5 regulates neural stem cell function. Combined with our transcriptome analysis, this mapping analysis yielded a bona fide map of Qki5-RNA interaction at single-nucleotide resolution, the identification of 892 Qki5 direct target genes, and an accurate Qki5-dependent alternative splicing rule in the developing brain. Last, our target gene list provides the first compelling evidence that Qki5 is associated with specific biological events; namely, cell-cell adhesion. This prediction was confirmed by histological analysis of mice in which Qki proteins were genetically ablated, which revealed disruption of the apical surface of the lateral wall in the developing brain. These data collectively indicate that Qki5 regulates communication between neural stem cells by mediating numerous RNA processing events and suggest new links between splicing regulation and neural stem cell states.

Preferential binding of a stable G3BP ribonucleoprotein complex to intron-retaining transcripts in mouse brain and modulation of their expression in the cerebellum

Neuronal granules play an important role in the localization and transport of translationally silenced messenger ribonucleoproteins in neurons. Among the factors associated with these granules, the RNA-binding protein G3BP1 (stress-granules assembly factor) is involved in neuronal plasticity and is induced in Alzheimer's disease. We immunopurified a stable complex containing G3BP1 from mouse brain and performed high-throughput sequencing and cross-linking immunoprecipitation to identify the associated RNAs. The G3BP-complex contained the deubiquitinating protease USP10, CtBP1 and the RNA-binding proteins Caprin-1, G3BP2a and splicing factor proline and glutamine rich, or PSF. The G3BP-complex binds preferentially to transcripts that retain introns, and to non-coding sequences like 3'-untranslated region and long non-coding RNAs. Specific transcripts with retained introns appear to be enriched in the cerebellum compared to the rest of the brain and G3BP1 depletion decreased this intron retention in the cerebellum of G3BP1 knockout mice. Among the enriched transcripts, we found an overrepresentation of genes involved in synaptic transmission, especially glutamate-related neuronal transmission. Notably, G3BP1 seems to repress the expression of the mature Grm5 (metabotropic glutamate receptor 5) transcript, by promoting the retention of an intron in the immature transcript in the cerebellum. Our results suggest that G3BP is involved in a new functional mechanism to regulate non-coding RNAs including intron-retaining transcripts, and thus have broad implications for neuronal gene regulation, where intron retention is widespread.

NOVA2-mediated RNA regulation is required for axonal pathfinding during development

The neuron specific RNA-binding proteins NOVA1 and NOVA2 are highly homologous alternative splicing regulators. NOVA proteins regulate at least 700 alternative splicing events in vivo, yet relatively little is known about the biologic consequences of NOVA action and in particular about functional differences between NOVA1 and NOVA2. Transcriptome-wide searches for isoform-specific functions, using NOVA1 and NOVA2 specific HITS-CLIP and RNA-seq data from mouse cortex lacking either NOVA isoform, reveals that NOVA2 uniquely regulates alternative splicing events of a series of axon guidance related genes during cortical development. Corresponding axonal pathfinding defects were specific to NOVA2 deficiency: Nova2-/- but not Nova1-/- mice had agenesis of the corpus callosum, and axonal outgrowth defects specific to ventral motoneuron axons and efferent innervation of the cochlea. Thus we have discovered that NOVA2 uniquely regulates alternative splicing of a coordinate set of transcripts encoding key components in cortical, brainstem and spinal axon guidance/outgrowth pathways during neural differentiation, with severe functional consequences in vivo.

Improvements to the HITS-CLIP protocol eliminate widespread mispriming artifacts

Background

 High-throughput sequencing of RNA isolated by crosslinking immunoprecipitation (HITS-CLIP) allows for high resolution, genome-wide mapping of RNA-binding proteins. This methodology is frequently used to validate predicted targets of microRNA binding, as well as direct targets of other RNA-binding proteins. Hence, the accuracy and sensitivity of binding site identification is critical. 

Results

We found that substantial mispriming during reverse transcription results in the overrepresentation of sequences complementary to the primer used for reverse transcription. Up to 45 % of peaks in publicly available HITS-CLIP libraries are attributable to this mispriming artifact, and the majority of libraries have detectable levels of mispriming. We also found that standard techniques for validating microRNA-target interactions fail to differentiate between artifactual peaks and physiologically relevant peaks.

Conclusions

Here, we present a modification to the HITS-CLIP protocol that effectively eliminates this artifact and improves the sensitivity and complexity of resulting libraries.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2675-5) contains supplementary material, which is available to authorized users.

CLIPSeqTools--a novel bioinformatics CLIP-seq analysis suite

Immunoprecipitation of RNA binding proteins (RBPs) after in vivo crosslinking, coupled with sequencing of associated RNA footprints (HITS-CLIP, CLIP-seq), is a method of choice for the identification of RNA targets and binding sites for RBPs. Compared with RNA-seq, CLIP-seq analysis is widely diverse and depending on the RBPs that are analyzed, the approaches vary significantly, necessitating the development of flexible and efficient informatics tools. In this study, we present CLIPSeqTools, a novel, highly flexible computational suite that can perform analysis from raw sequencing data with minimal user input. It contains a wide array of tools to provide an in-depth view of CLIP-seq data sets. It supports extensive customization and promotes improvization, a critical virtue, since CLIP-seq analysis is rarely well defined a priori. To highlight CLIPSeqTools capabilities, we used the suite to analyze Ago-miRNA HITS-CLIP data sets that we prepared from human brains.

Experimental and Computational Considerations in the Study of RNA-Binding Protein-RNA Interactions

After an RNA is transcribed, it undergoes a variety of processing steps that can change the encoded protein sequence (through alternative splicing and RNA editing), regulate the stability of the RNA, and control subcellular localization, timing, and rate of translation. The recent explosion in genomics techniques has enabled transcriptome-wide profiling of RNA processing in an unbiased manner. However, it has also brought with it both experimental challenges in developing improved methods to probe distinct processing steps, as well as computational challenges in data storage, processing, and analysis tools to enable large-scale interpretation in the genomics era. In this chapter we review experimental techniques and challenges in profiling various aspects of RNA processing, as well as recent efforts to develop analyses integrating multiple data sources and techniques to infer RNA regulatory networks.

Analyzing MiRNA-LncRNA Interactions

Long noncoding RNAs (lncRNAs) are noncoding transcripts usually longer than 200 nts that have recently emerged as one of the largest and significantly diverse RNA families. The biological role and functions of lncRNAs are still mostly uncharacterized. Their target-mimetic, sponge/decoy function on microRNAs was recently uncovered. miRNAs are a class of noncoding RNA species (~22 nts) that play a central role in posttranscriptional regulation of protein coding genes by mRNA cleavage, direct translational repression and/or mRNA destabilization. LncRNAs can act as miRNA sponges, reducing their regulatory effect on mRNAs. This function introduces an extra layer of complexity in the miRNA-target interaction network. This chapter focuses on the study of miRNA-lncRNA interactions with either in silico or experimentally supported analyses. The proposed methodologies can be appropriately adapted in order to become the backbone of advanced multistep functional miRNA analyses.

A multiprotein occupancy map of the mRNP on the 3' end of histone mRNAs

The animal replication-dependent (RD) histone mRNAs are coordinately regulated with chromosome replication. The RD-histone mRNAs are the only known cellular mRNAs that are not polyadenylated. Instead, the mature transcripts end in a conserved stem-loop (SL) structure. This SL structure interacts with the stem-loop binding protein (SLBP), which is involved in all aspects of RD-histone mRNA metabolism. We used several genomic methods, including high-throughput sequencing of cross-linked immunoprecipitate (HITS-CLIP) to analyze the RNA-binding landscape of SLBP. SLBP was not bound to any RNAs other than histone mRNAs. We performed bioinformatic analyses of the HITS-CLIP data that included (i) clustering genes by sequencing read coverage using CVCA, (ii) mapping the bound RNA fragment termini, and (iii) mapping cross-linking induced mutation sites (CIMS) using CLIP-PyL software. These analyses allowed us to identify specific sites of molecular contact between SLBP and its RD-histone mRNA ligands. We performed in vitro crosslinking assays to refine the CIMS mapping and found that uracils one and three in the loop of the histone mRNA SL preferentially crosslink to SLBP, whereas uracil two in the loop preferentially crosslinks to a separate component, likely the 3'hExo. We also performed a secondary analysis of an iCLIP data set to map UPF1 occupancy across the RD-histone mRNAs and found that UPF1 is bound adjacent to the SLBP-binding site. Multiple proteins likely bind the 3' end of RD-histone mRNAs together with SLBP.

elegans

The evolutionarily conserved gene lin-28 encodes an RNA-binding protein and is an important regulator of the proper temporal succession of several developmental events in both invertebrates and vertebrates. At the cellular level, LIN-28 promotes stemness and proliferation, and inhibits differentiation, a feature best illustrated by its ability to induce pluripotency when ectopically expressed in human fibroblasts in combination with NANOG, OCT4, and SOX2. Mammalian LIN28 functions in part by regulating processing of the let-7 microRNA through a GGAG binding site in the pre-let-7's distal loop region. However, many human and animal let-7 precursors lack the GGAG binding motif. In order to dissect the molecular mechanisms underlying its biological functions in a living animal, we identified a map of LIN-28 interactions with the transcriptome by in vivo HITS-CLIP in Caenorhabditis elegans. LIN-28 binds a large pool of messenger RNAs, and a substantial fraction of the bona fide LIN-28 targets are involved in aspects of animal development. Furthermore, our data show that LIN-28 regulates the expression of the let-7 microRNA by binding its primary transcript in a previously unknown region, revealing a novel regulatory mechanism.

Computational Methods for CLIP-seq Data Processing

RNA-binding proteins (RBPs) are at the core of post-transcriptional regulation and thus of gene expression control at the RNA level. One of the principal challenges in the field of gene expression regulation is to understand RBPs mechanism of action. As a result of recent evolution of experimental techniques, it is now possible to obtain the RNA regions recognized by RBPs on a transcriptome-wide scale. In fact, CLIP-seq protocols use the joint action of CLIP, crosslinking immunoprecipitation, and high-throughput sequencing to recover the transcriptome-wide set of interaction regions for a particular protein. Nevertheless, computational methods are necessary to process CLIP-seq experimental data and are a key to advancement in the understanding of gene regulatory mechanisms. Considering the importance of computational methods in this area, we present a review of the current status of computational approaches used and proposed for CLIP-seq data.

NOVA-dependent regulation of cryptic NMD exons controls synaptic protein levels after seizure

The neuronal RNA binding protein NOVA regulates splicing, shuttles to the cytoplasm, and co-localizes with target transcripts in dendrites, suggesting links between splicing and local translation. Here we identified >200 transcripts showing NOVA-dependent changes in abundance, but, surprisingly, HITS-CLIP revealed NOVA binds these RNAs in introns rather than 3′ UTRs. This led us to discover NOVA-regulated splicing of cryptic exons within these introns. These exons triggered nonsense mediated decay (NMD), as UPF1 and protein synthesis were required for NOVA's effect on RNA levels. Their regulation was dynamic and physiologically relevant. The NMD exons were regulated by seizures, which also induced changes in Nova subcellular localization and mediated large changes in synaptic proteins, including proteins implicated in familial epilepsy. Moreover, Nova haploinsufficient mice had spontaneous epilepsy. The data reveal a hidden means of dynamic RNA regulation linking electrical activity to splicing and protein output, and of mediating homeostatic excitation/inhibition balance in neurons.

DOI:http://dx.doi.org/10.7554/eLife.00178.001

After the DNA in a gene has been transcribed into messenger RNA, portions of the mRNA called introns are removed, and the remaining stretches of mRNA, which are known as exons, are spliced together. Within eukaryotic cells, a process known as alternative splicing allows a single gene to encode for multiple protein variants by ensuring that some exons are included in the final, modified mRNA, while other exons are excluded. This modified mRNA is then translated into proteins.

Eukaryotic cells also contain proteins that bind to RNA to regulate alternative splicing. These RNA-binding proteins are often found in both the cytoplasm and nucleus of cells, and their involvement in splicing may be linked to other processes in the cell such as mRNA localization and translation. It has also become clear over the past two decades that certain types of RNA-binding proteins, including NOVA proteins, are only found in neurons, and that these proteins have been best characterized as alternative splicing regulators. Recent work has also suggested that they also have important roles in regulating neuronal activity and development, and that their actions in neuronal nuclei and cytoplasm might be coordinated.

Now Eom et al. use the predictive power of a high throughput sequencing and crosslinking method termed HITS-CLIP to show that NOVA proteins can indirectly regulate cytoplasmic mRNA levels by regulating the process of alternative splicing in the nucleus to produce ‘cryptic’ exons in the brains of mice. The presence of these exons in the mRNA leads to the production of premature termination codons in the cytoplasm. These codons trigger a process called nonsense-mediated decay that involves identifying mRNA transcripts that contain nonsense mutations, and then degrading them. These cryptic exons were seen in mice missing the NOVA proteins, where they are expressed in abnormally high levels; in normal mice, these exons have not been seen before, hence they were termed ‘cryptic’.

Eom et al. also show that these cryptic exons are physiologically relevant by inducing epileptic seizures in mice. Following the seizures, they find that the NOVA proteins up-regulate and down-regulate the levels of different cryptic exons, leading to changes in the levels of the proteins encoded by these mRNAs, including proteins that inhibit further seizures. Overall the results indicate that, by controlling the production of various proteins in neurons, these previously unknown cryptic exons have important roles in the workings of the brain.

DOI:http://dx.doi.org/10.7554/eLife.00178.002

'Cryptic' exons reveal some of their secrets

By regulating the inclusion of ‘cryptic’ exons in messenger RNA in nerve cells, NOVA proteins are able to influence the abundance of the corresponding proteins.

The Neuronal Splicing Factor Nova Co-Localizes with Target RNAs in the Dendrite

Nova proteins are neuron-specific RNA binding proteins targeted by autoantibodies in a disorder manifest by failure of motor inhibition, and they regulate splicing and alternative 3' processing. Nova regulates splicing of RNAs encoding synaptic proteins, including the inhibitory glycine receptor alpha2 subunit (GlyRalpha2), and binds to others, including the GIRK2 channel. We found that Nova harbors functional NES and NLS elements, shuttles between the nucleus and cytoplasm, and that 50% of the protein localizes to the soma-dendritic compartment. Immunofluoresence and EM analysis of spinal cord motor neurons demonstrated that Nova co-localizes beneath synaptic contacts in dendrites with the same RNA, GlyRalpha2, whose splicing it regulates in the nucleus. HITS-CLIP identified intronic and 3' UTR sites where Nova binds to GlyRalpha2 and GIRK2 transcripts in the brain. This led directly to the identification of a 3' UTR localization element that mediates Nova-dependent localization of GIRK2 in primary neurons. These data demonstrate that HITS-CLIP can identify functional RNA localization elements, and they suggest new links between the regulation of nuclear RNA processing and mRNA localization.