RIP: an mRNA localization technique

MiR-574-5p activates human TLR8 to promote autoimmune signaling and lupus

Endosomal single-stranded RNA-sensing Toll-like receptor-7/8 (TLR7/8) plays a pivotal role in inflammation and immune responses and autoimmune diseases. However, the mechanisms underlying the initiation of the TLR7/8-mediated autoimmune signaling remain to be fully elucidated. Here, we demonstrate that miR-574-5p is aberrantly upregulated in tissues of lupus prone mice and in the plasma of lupus patients, with its expression levels correlating with the disease activity. miR-574-5p binds to and activates human hTLR8 or its murine ortholog mTlr7 to elicit a series of MyD88-dependent immune and inflammatory responses. These responses include the overproduction of cytokines and interferons, the activation of STAT1 signaling and B lymphocytes, and the production of autoantigens. In a transgenic mouse model, the induction of miR-574-5p overexpression is associated with increased secretion of antinuclear and anti-dsDNA antibodies, increased IgG and C3 deposit in the kidney, elevated expression of inflammatory genes in the spleen. In lupus-prone mice, lentivirus-mediated silencing of miR-574-5p significantly ameliorates major symptoms associated with lupus and lupus nephritis. Collectively, these results suggest that the miR-574-5p-hTLR8/mTlr7 signaling is an important axis of immune and inflammatory responses, contributing significantly to the development of lupus and lupus nephritis.

COPI coatomer subunit α-COP interacts with the RNA binding protein Nucleolin via a C-terminal dilysine motif

The COPI coatomer subunit α-COP has been shown to co-precipitate mRNA in multiple settings, but it was unclear whether the interaction with mRNA was direct or mediated by interaction with an adapter protein. The COPI complex often interacts with proteins via C-terminal dilysine domains. A search for candidate RNA binding proteins with C-terminal dilysine motifs yielded Nucleolin, which terminates in a KKxKxx sequence. This protein was an especially intriguing candidate as it has been identified as an interacting partner for Survival Motor Neuron protein (SMN). Loss of SMN causes the neurodegenerative disease Spinal Muscular Atrophy. We have previously shown that SMN and α-COP interact and co-migrate in axons, and that overexpression of α-COP reduced phenotypic severity in cell culture and animal models of SMA. We show here that in an mRNA independent manner, endogenous Nucleolin co-precipitates endogenous α-COP and ε-COP but not β-COP which may reflect an interaction with the so-called B-subcomplex rather a complete COPI heptamer. The ability of Nucleolin to bind to α-COP requires the presence of the C-terminal KKxKxx domain of Nucleolin. Furthermore, we have generated a point mutant in the WD40 domain of α-COP which eliminates its ability to co-precipitate Nucleolin but does not interfere with precipitation of partners mediated by non-KKxKxx motifs such as the kainate receptor subunit 2. We propose that via interaction between the C-terminal dilysine motif of Nucleolin and the WD40 domain of α-COP, Nucleolin acts an adaptor to allow α-COP to interact with a population of mRNA.

STAU2 binds a complex RNA cargo that changes temporally with production of diverse intermediate progenitor cells during mouse corticogenesis

STAU2 is a double-stranded RNA-binding protein enriched in the nervous system. During asymmetric divisions in the developing mouse cortex, STAU2 preferentially distributes into the intermediate progenitor cell (IPC), delivering RNA molecules that can impact IPC behavior. Corticogenesis occurs on a precise time schedule, raising the hypothesis that the cargo STAU2 delivers into IPCs changes over time. To test this, we combine RNA-immunoprecipitation with sequencing (RIP-seq) over four stages of mouse cortical development, generating a comprehensive cargo profile for STAU2. A subset of the cargo was 'stable', present at all stages, and involved in chromosome organization, macromolecule localization, translation and DNA repair. Another subset was 'dynamic', changing with cortical stage, and involved in neurogenesis, cell projection organization, neurite outgrowth, and included cortical layer markers. Notably, the dynamic STAU2 cargo included determinants of IPC versus neuronal fates and genes contributing to abnormal corticogenesis. Knockdown of one STAU2 target, Taf13, previously linked to microcephaly and impaired myelination, reduced oligodendrogenesis in vitro. We conclude that STAU2 contributes to the timing of corticogenesis by binding and delivering complex and temporally regulated RNA cargo into IPCs.

Episo: quantitative estimation of RNA 5-methylcytosine at isoform level by high-throughput sequencing of RNA treated with bisulfite

Motivation

RNA 5-methylcytosine (m⁵C) is a type of post-transcriptional modification that may be involved in numerous biological processes and tumorigenesis. RNA m⁵C can be profiled at single-nucleotide resolution by high-throughput sequencing of RNA treated with bisulfite (RNA-BisSeq). However, the exploration of transcriptome-wide profile and potential function of m⁵C in splicing remains to be elucidated due to lack of isoform level m⁵C quantification tool.

Results

We developed a computational package to quantify Epitranscriptomal RNA m⁵C at the transcript isoform level (named Episo). Episo consists of three tools: mapper, quant and Bisulfitefq, for mapping, quantifying and simulating RNA-BisSeq data, respectively. The high accuracy of Episo was validated using an improved m⁵C-specific methylated RNA immunoprecipitation (meRIP) protocol, as well as a set of in silico experiments. By applying Episo to public human and mouse RNA-BisSeq data, we found that the RNA m⁵C is not evenly distributed among the transcript isoforms, implying the m⁵C may subject to be regulated at isoform level.

Availability and implementation

Episo is released under the GNU GPLv3+ license. The resource code Episo is freely accessible from https://github.com/liujunfengtop/Episo (with Tophat/cufflink) and https://github.com/liujunfengtop/Episo/tree/master/Episo_Kallisto (with Kallisto).

Supplementary information

Supplementary data are available at Bioinformatics online.

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.

Engineering Structurally Interacting RNA (sxRNA)

RNA-based three-way junctions (3WJs) are naturally occurring structures found in many functional RNA molecules including rRNA, tRNA, snRNA and ribozymes. 3WJs are typically characterized as resulting from an RNA molecule folding back on itself in cis but could also form in trans when one RNA, for instance a microRNA binds to a second structured RNA, such as a mRNA. Trans-3WJs can influence the final shape of one or both of the RNA molecules and can thus provide a means for modulating the availability of regulatory motifs including potential protein or microRNA binding sites. Regulatory 3WJs generated in trans represent a newly identified regulatory category that we call structurally interacting RNA or sxRNA for convenience. Here we show that they can be rationally designed using familiar cis-3WJ examples as a guide. We demonstrate that an sxRNA "bait" sequence can be designed to interact with a specific microRNA "trigger" sequence, creating a regulatable RNA-binding protein motif that retains its functional activity. Further, we show that when placed downstream of a coding sequence, sxRNA can be used to switch "ON" translation of that sequence in the presence of the trigger microRNA and the amount of translation corresponded with the amount of microRNA present.

Stem Cell Ribonomics: RNA-Binding Proteins and Gene Networks in Stem Cell Differentiation

Stem cells are undifferentiated cells with the ability to self-renew and the potential to differentiate into all body cell types. Stem cells follow a developmental genetic program and are able to respond to alterations in the environment through various signaling pathways. The mechanisms that control these processes involve the activation of transcription followed by a series of post-transcriptional events. These post-transcriptional steps are mediated by the interaction of RNA-binding proteins (RBPs) with defined subpopulations of RNAs creating a regulatory gene network. Characterizing these RNA-protein networks is essential to understanding the regulatory mechanisms underlying the control of stem cell fate. Ribonomics is the combination of classical biochemical purification protocols with the high-throughput identification of transcripts applied to the functional characterization of RNA-protein complexes. Here, we describe the different approaches that can be used in a ribonomic approach and how they have contributed to understanding the function of several RBPs with central roles in stem cell biology.

Identification of common genetic variants controlling transcript isoform variation in human whole blood

An understanding of the genetic variation underlying transcript splicing is essential to dissect the molecular mechanisms of common disease. The available evidence from splicing quantitative trait locus (sQTL) studies has been limited to small samples. We performed genome-wide screening to identify SNPs that might control mRNA splicing in whole blood collected from 5,257 Framingham Heart Study participants. We identified 572,333 cis sQTLs involving 2,650 unique genes. Many sQTL-associated genes (40%) undergo alternative splicing. Using the National Human Genome Research Institute (NHGRI) genome-wide association study (GWAS) catalog, we determined that 528 unique sQTLs were significantly enriched for 8,845 SNPs associated with traits in previous GWAS. In particular, we found 395 (4.5%) GWAS SNPs with evidence of cis sQTLs but not gene-level cis expression quantitative trait loci (eQTLs), suggesting that sQTL analysis could provide additional insights into the functional mechanism underlying GWAS results. Our findings provide an informative sQTL resource for further characterizing the potential functional roles of SNPs that control transcript isoforms relevant to common diseases.

RIP-seq of BmAgo2-associated small RNAs reveal various types of small non-coding RNAs in the silkworm, Bombyx mori

Background

Small non-coding RNAs (ncRNAs) are important regulators of gene expression in eukaryotes. Previously, only microRNAs (miRNAs) and piRNAs have been identified in the silkworm, Bombyx mori. Furthermore, only ncRNAs (50-500nt) of intermediate size have been systematically identified in the silkworm.

Results

Here, we performed a systematic identification and analysis of small RNAs (18-50nt) associated with the Bombyx mori argonaute2 (BmAgo2) protein. Using RIP-seq, we identified various types of small ncRNAs associated with BmAGO2. These ncRNAs showed a multimodal length distribution, with three peaks at ~20nt, ~27nt and ~33nt, which included tRNA-, transposable element (TE)-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. The tRNA-derived fragments (tRFs) were found at an extremely high abundance and accounted for 69.90% of the BmAgo2-associated small RNAs. Northern blotting confirmed that many tRFs were expressed or up-regulated only in the BmNPV-infected cells, implying that the tRFs play a prominent role by binding to BmAgo2 during BmNPV infection. Additional evidence suggested that there are potential cleavage sites on the D, anti-codon and TψC loops of the tRNAs. TE-derived small RNAs and piRNAs also accounted for a significant proportion of the BmAgo2-associated small RNAs, suggesting that BmAgo2 could be involved in the maintenance of genome stability by suppressing the activities of transposons guided by these small RNAs. Finally, Northern blotting was also used to confirm the Bombyx 5.8 s rRNA-derived small RNAs, demonstrating that various novel small RNAs exist in the silkworm.

Conclusions

Using an RIP-seq method in combination with Northern blotting, we identified various types of small RNAs associated with the BmAgo2 protein, including tRNA-, TE-, rRNA-, snoRNA- and snRNA-derived small RNAs as well as miRNAs and piRNAs. Our findings provide new clues for future functional studies of the role of small RNAs in insect development and evolution.

An in-depth map of polyadenylation sites in cancer

We present a comprehensive map of over 1 million polyadenylation sites and quantify their usage in major cancers and tumor cell lines using direct RNA sequencing. We built the Expression and Polyadenylation Database to enable the visualization of the polyadenylation maps in various cancers and to facilitate the discovery of novel genes and gene isoforms that are potentially important to tumorigenesis. Analyses of polyadenylation sites indicate that a large fraction (∼30%) of mRNAs contain alternative polyadenylation sites in their 3' untranslated regions, independent of the cell type. The shortest 3' untranslated region isoforms are preferentially upregulated in cancer tissues, genome-wide. Candidate targets of alternative polyadenylation-mediated upregulation of short isoforms include POLR2K, and signaling cascades of cell-cell and cell-extracellular matrix contact, particularly involving regulators of Rho GTPases. Polyadenylation maps also helped to improve 3' untranslated region annotations and identify candidate regulatory marks such as sequence motifs, H3K36Me3 and Pabpc1 that are isoform dependent and occur in a position-specific manner. In summary, these results highlight the need to go beyond monitoring only the cumulative transcript levels for a gene, to separately analysing the expression of its RNA isoforms.

Two-timing zinc finger transcription factors liaising with RNA

Classical zinc fingers (ZFs) are one of the most common protein domains in higher eukaryotes and have been known for almost 30 years to act as sequence-specific DNA-binding domains. This knowledge has come, however, from the study of a small number of archetypal proteins, and a larger picture is beginning to emerge that ZF functions are far more diverse than originally suspected. Here, we review the evidence that a subset of ZF proteins live double lives, binding to both DNA and RNA targets and frequenting both the cytoplasm and the nucleus. This duality can create an important additional level of gene regulation that serves to connect transcriptional and post-transcriptional control.