Improved discovery of RNA-binding protein binding sites in eCLIP data using DEWSeq

irCLIP-RNP and Re-CLIP reveal patterns of dynamic protein assemblies on RNA

RNA-binding proteins (RBPs) control varied processes, including RNA splicing, stability, transport and translation1-3. Dysfunctional RNA-RBP interactions contribute to the pathogenesis of human disease1,4,5; however, characterizing the nature and dynamics of multiprotein assemblies on RNA has been challenging. Here, to address this, non-isotopic ligation-based ultraviolet-light-induced cross-linking and immunoprecipitation6 was combined with mass spectrometry (irCLIP-RNP) to identify RNA-dependent associated proteins (RDAPs) co-bound to RNA with any RBP of interest. irCLIP-RNP defined landscapes of multimeric protein assemblies on RNA, revealing patterns of RBP-RNA associations, including cell-type-selective combinatorial relationships between RDAPs and primary RBPs. irCLIP-RNP also defined dynamic RDAP remodelling in response to epidermal growth factor (EGF), revealing that EGF-induced recruitment of UPF1 adjacent to HNRNPC promotes splicing surveillance of cell proliferation mRNAs. To identify the RNAs simultaneously co-bound by multiple studied RBPs, a sequential immunoprecipitation irCLIP (Re-CLIP) method was also developed. Re-CLIP confirmed binding relationships observed in irCLIP-RNP and identified HNRNPC and UPF1 RBP co-binding on RND3 and DDX3X mRNAs. irCLIP-RNP and Re-CLIP provide a framework to identify and characterize dynamic RNA-protein assemblies in living cells.

Imp/IGF2BP and Syp/SYNCRIP temporal RNA interactomes uncover combinatorial networks of regulators of Drosophila brain development

Temporal patterning of neural progenitors is an evolutionarily conserved mechanism generating neural diversity. In Drosophila, postembryonic neurogenesis requires the RNA binding proteins (RBPs) Imp/IGF2BP and Syp/SYNCRIP. However, how they coachieve their function is not well understood. Here, we elucidate the in vivo temporal RNA interactome landscapes of Imp and Syp during larval brain development. Imp and Syp bind a highly overlapping set of conserved mRNAs encoding proteins involved in neurodevelopment. We identify transcripts differentially occupied by Imp/Syp over time, featuring a network of known and previously unknown candidate temporal regulators that are post-transcriptionally regulated by Imp/Syp. Furthermore, the physical and coevolutionary relationships between Imp and Syp binding sites reveal a combinatorial, rather than competitive, mode of molecular interplay. Our study establishes an in vivo framework for dissecting the temporal coregulation of RBP networks as well as providing a resource for understanding neural fate specification.

RNA helicase D1PAS1 resolves R-loops and forms a complex for mouse pachytene piRNA biogenesis required for male fertility

During meiosis, RNA polymerase II transcribes pachytene piRNA precursors with unusually long and unspliced transcripts from discrete autosomal loci in the mouse genome. Despite the importance of piRNA for male fertility and a well-defined maturation process, the transcriptional machinery remains poorly understood. Here, we document that D1PAS1, an ATP-dependent RNA helicase, is critical for pachytene piRNA expression from multiple genomic loci and subsequent translocation into the cytoplasm to ensure mature piRNA biogenesis. Depletion of D1PAS1 in gene-edited mice results in the accumulation of R-loops in pachytene spermatocytes, leading to DNA-damage-induced apoptosis, disruption of piRNA biogenesis, spermatogenic arrest, and male infertility. Transcriptome, genome-wide R-loop profiling, and proteomic analyses document that D1PAS1 regulates pachytene piRNA transcript elongation and termination. D1PAS1 subsequently forms a complex with nuclear export components to ensure pachytene piRNA precursor translocation from the nucleus to the cytoplasm for processing into small non-coding RNAs. Thus, our study defines D1PAS1 as a specific transcription activator that promotes R-loop unwinding and is a critical factor in pachytene piRNA biogenesis.

APC mutations dysregulate alternative polyadenylation in cancer

BACKGROUND

Alternative polyadenylation (APA) affects most human genes and is recurrently dysregulated in all studied cancers. However, the mechanistic origins of this dysregulation are incompletely understood.

RESULTS

We describe an unbiased analysis of molecular regulators of poly(A) site selection across The Cancer Genome Atlas and identify that colorectal adenocarcinoma is an outlier relative to all other cancer subtypes. This distinction arises from the frequent presence of loss-of-function APC mutations in colorectal adenocarcinoma, which are strongly associated with long 3' UTR expression relative to tumors lacking APC mutations. APC knockout similarly dysregulates APA in human colon organoids. By mining previously published APC eCLIP data, we show that APC preferentially binds G- and C-rich motifs just upstream of proximal poly(A) sites. Lastly, we find that reduced APC expression is associated with APA dysregulation in tumor types lacking recurrent APC mutations.

CONCLUSIONS

As APC has been previously identified as an RNA-binding protein that preferentially binds 3' UTRs during mouse neurogenesis, our results suggest that APC promotes proximal poly(A) site use and that APC loss and altered expression contribute to pervasive APA dysregulation in cancers.

Cytosolic RNA binding of the mitochondrial TCA cycle enzyme malate dehydrogenase

Several enzymes of intermediary metabolism have been identified to bind RNA in cells, with potential consequences for the bound RNAs and/or the enzyme. In this study, we investigate the RNA-binding activity of the mitochondrial enzyme malate dehydrogenase 2 (MDH2), which functions in the tricarboxylic acid (TCA) cycle and the malate-aspartate shuttle. We confirmed in cellulo RNA binding of MDH2 using orthogonal biochemical assays and performed enhanced cross-linking and immunoprecipitation (eCLIP) to identify the cellular RNAs associated with endogenous MDH2. Surprisingly, MDH2 preferentially binds cytosolic over mitochondrial RNAs, although the latter are abundant in the milieu of the mature protein. Subcellular fractionation followed by RNA-binding assays revealed that MDH2-RNA interactions occur predominantly outside of mitochondria. We also found that a cytosolically retained N-terminal deletion mutant of MDH2 is competent to bind RNA, indicating that mitochondrial targeting is dispensable for MDH2-RNA interactions. MDH2 RNA binding increased when cellular NAD+ levels (MDH2's cofactor) were pharmacologically diminished, suggesting that the metabolic state of cells affects RNA binding. Taken together, our data implicate an as yet unidentified function of MDH2-binding RNA in the cytosol.

BestCRM: An Exhaustive Search for Optimal Cis-Regulatory Modules in Promoters Accelerated by the Multidimensional Hash Function

The concept of cis-regulatory modules located in gene promoters represents today's vision of the organization of gene transcriptional regulation. Such modules are a combination of two or more single, short DNA motifs. The bioinformatic identification of such modules belongs to so-called NP-hard problems with extreme computational complexity, and therefore, simplifications, assumptions, and heuristics are usually deployed to tackle the problem. In practice, this requires, first, many parameters to be set before the search, and second, it leads to the identification of locally optimal results. Here, a novel method is presented, aimed at identifying the cis-regulatory elements in gene promoters based on an exhaustive search of all the feasible modules' configurations. All required parameters are automatically estimated using positive and negative datasets. To be computationally efficient, the search is accelerated using a multidimensional hash function, allowing the search to complete in a few hours on a regular laptop (for example, a CPU Intel i7, 3.2 GH, 32 Gb RAM). Tests on an established benchmark and real data show better performance of BestCRM compared to the available methods according to several metrics like specificity, sensitivity, AUC, etc. A great practical advantage of the method is its minimum number of input parameters-apart from positive and negative promoters, only a desired level of module presence in promoters is required.