Mathias Munschauer

Transcriptome-wide identification of RNA-binding protein and microRNA target sites by PAR-CLIP

RNA transcripts are subject to posttranscriptional gene regulation involving hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) expressed in a cell-type dependent fashion. We developed a cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs. The crosslinked sites are revealed by thymidine to cytidine transitions in the cDNAs prepared from immunopurified RNPs of 4-thiouridine-treated cells. We determined the binding sites and regulatory consequences for several intensely studied RBPs and miRNPs, including PUM2, QKI, IGF2BP1-3, AGO/EIF2C1-4 and TNRC6A-C. Our study revealed that these factors bind thousands of sites containing defined sequence motifs and have distinct preferences for exonic versus intronic or coding versus untranslated transcript regions. The precise mapping of binding sites across the transcriptome will be critical to the interpretation of the rapidly emerging data on genetic variation between individuals and how these variations contribute to complex genetic diseases.

Circular RNAs are a large class of animal RNAs with regulatory potency

Circular RNAs (circRNAs) in animals are an enigmatic class of RNA with unknown function. To explore circRNAs systematically, we sequenced and computationally analysed human, mouse and nematode RNA. We detected thousands of well-expressed, stable circRNAs, often showing tissue/developmental-stage-specific expression. Sequence analysis indicated important regulatory functions for circRNAs. We found that a human circRNA, antisense to the cerebellar degeneration-related protein 1 transcript (CDR1as), is densely bound by microRNA (miRNA) effector complexes and harbours 63 conserved binding sites for the ancient miRNA miR-7. Further analyses indicated that CDR1as functions to bind miR-7 in neuronal tissues. Human CDR1as expression in zebrafish impaired midbrain development, similar to knocking down miR-7, suggesting that CDR1as is a miRNA antagonist with a miRNA-binding capacity ten times higher than any other known transcript. Together, our data provide evidence that circRNAs form a large class of post-transcriptional regulators. Numerous circRNAs form by head-to-tail splicing of exons, suggesting previously unrecognized regulatory potential of coding sequences.

FMRP targets distinct mRNA sequence elements to regulate protein expression

Fragile X syndrome (FXS) is a multi-organ disease that leads to mental retardation, macro-orchidism in males and premature ovarian insufficiency in female carriers. FXS is also a prominent monogenic disease associated with autism spectrum disorders (ASDs). FXS is typically caused by the loss of fragile X mental retardation 1 (FMR1) expression, which codes for the RNA-binding protein FMRP. Here we report the discovery of distinct RNA-recognition elements that correspond to the two independent RNA-binding domains of FMRP, in addition to the binding sites within the messenger RNA targets for wild-type and I304N mutant FMRP isoforms and the FMRP paralogues FXR1P and FXR2P (also known as FXR1 and FXR2). RNA-recognition-element frequency, ratio and distribution determine target mRNA association with FMRP. Among highly enriched targets, we identify many genes involved in ASD and show that FMRP affects their protein levels in human cell culture, mouse ovaries and human brain. Notably, we discovered that these targets are also dysregulated in Fmr1−/− mouse ovaries showing signs of premature follicular overdevelopment. These results indicate that FMRP targets share signalling pathways across different cellular contexts. As the importance of signalling pathways in both FXS and ASD is becoming increasingly apparent, our results provide a ranked list of genes as basis for the pursuit of new therapeutic targets for these neurological disorders.

PAR-CliP-a method to identify transcriptome-wide the binding sites of RNA binding proteins

RNA transcripts are subjected to post-transcriptional gene regulation by interacting with hundreds of RNA-binding proteins (RBPs) and microRNA-containing ribonucleoprotein complexes (miRNPs) that are often expressed in a cell-type dependently. To understand how the interplay of these RNA-binding factors affects the regulation of individual transcripts, high resolution maps of in vivo protein-RNA interactions are necessary1. A combination of genetic, biochemical and computational approaches are typically applied to identify RNA-RBP or RNA-RNP interactions. Microarray profiling of RNAs associated with immunopurified RBPs (RIP-Chip)2 defines targets at a transcriptome level, but its application is limited to the characterization of kinetically stable interactions and only in rare cases3,4 allows to identify the RBP recognition element (RRE) within the long target RNA. More direct RBP target site information is obtained by combining in vivo UV crosslinking5,6 with immunoprecipitation7-9 followed by the isolation of crosslinked RNA segments and cDNA sequencing (CLIP)10. CLIP was used to identify targets of a number of RBPs11-17. However, CLIP is limited by the low efficiency of UV 254 nm RNA-protein crosslinking, and the location of the crosslink is not readily identifiable within the sequenced crosslinked fragments, making it difficult to separate UV-crosslinked target RNA segments from background non-crosslinked RNA fragments also present in the sample. We developed a powerful cell-based crosslinking approach to determine at high resolution and transcriptome-wide the binding sites of cellular RBPs and miRNPs that we term PAR-CliP (Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation) (see Fig. 1A for an outline of the method). The method relies on the incorporation of photoreactive ribonucleoside analogs, such as 4-thiouridine (4-SU) and 6-thioguanosine (6-SG) into nascent RNA transcripts by living cells. Irradiation of the cells by UV light of 365 nm induces efficient crosslinking of photoreactive nucleoside-labeled cellular RNAs to interacting RBPs. Immunoprecipitation of the RBP of interest is followed by isolation of the crosslinked and coimmunoprecipitated RNA. The isolated RNA is converted into a cDNA library and deep sequenced using Solexa technology. One characteristic feature of cDNA libraries prepared by PAR-CliP is that the precise position of crosslinking can be identified by mutations residing in the sequenced cDNA. When using 4-SU, crosslinked sequences thymidine to cytidine transition, whereas using 6-SG results in guanosine to adenosine mutations. The presence of the mutations in crosslinked sequences makes it possible to separate them from the background of sequences derived from abundant cellular RNAs. Application of the method to a number of diverse RNA binding proteins was reported in Hafner et al.18

Identification of LIN28B-bound mRNAs reveals features of target recognition and regulation

The conserved human LIN28 RNA-binding proteins function in development, maintenance of pluripotency and oncogenesis. We used PAR-CLIP and a newly developed variant of this method, iDo-PAR-CLIP, to identify LIN28B targets as well as sites bound by the individual RNA-binding domains of LIN28B in the human transcriptome at nucleotide resolution. The position of target binding sites reflected the known structural relative orientation of individual LIN28B-binding domains, validating iDo-PAR-CLIP. Our data suggest that LIN28B directly interacts with most expressed mRNAs and members of the let-7 microRNA family. The Lin28-binding motif detected in pre-let-7 was enriched in mRNA sequences bound by LIN28B. Upon LIN28B knockdown, cell proliferation and the cell cycle were strongly impaired. Quantitative shotgun proteomics of LIN28B depleted cells revealed significant reduction of protein synthesis from its RNA targets. Computational analyses provided evidence that the strength of protein synthesis reduction correlated with the location of LIN28B binding sites within target transcripts.

Differential protein occupancy profiling of the mRNA transcriptome

BackgroundRNA-binding proteins (RBPs) mediate mRNA biogenesis, translation and decay. We recently developed an approach to profile transcriptome-wide RBP contacts on polyadenylated transcripts by next-generation sequencing. A comparison of such profiles from different biological conditions has the power to unravel dynamic changes in protein-contacted cis-regulatory mRNA regions without a priori knowledge of the regulatory protein component.ResultsWe compared protein occupancy profiles of polyadenylated transcripts in MCF7 and HEK293 cells. Briefly, we developed a bioinformatics workflow to identify differential crosslinking sites in cDNA reads of 4-thiouridine crosslinked polyadenylated RNA samples. We identified 30,000 differential crosslinking sites between MCF7 and HEK293 cells at an estimated false discovery rate of 10%. 73% of all reported differential protein-RNA contact sites cannot be explained by local changes in exon usage as indicated by complementary RNA-seq data. The majority of differentially crosslinked positions are located in 3′ UTRs, show distinct secondary-structure characteristics and overlap with binding sites of known RBPs, such as ELAVL1. Importantly, mRNA transcripts with the most significant occupancy changes show elongated mRNA half-lives in MCF7 cells.ConclusionsWe present a global comparison of protein occupancy profiles from different cell types, and provide evidence for altered mRNA metabolism as a result of differential protein-RNA contacts. Additionally, we introduce POPPI, a bioinformatics workflow for the analysis of protein occupancy profiling experiments. Our work demonstrates the value of protein occupancy profiling for assessing cis-regulatory RNA sequence space and its dynamics in growth, development and disease.

Comprehensive Protein Interactome Analysis of a Key RNA Helicase: Detection of Novel Stress Granule Proteins

DDX6 (p54/RCK) is a human RNA helicase with central roles in mRNA decay and translation repression. To help our understanding of how DDX6 performs these multiple functions, we conducted the first unbiased, large-scale study to map the DDX6-centric protein-protein interactome using immunoprecipitation and mass spectrometry. Using DDX6 as bait, we identify a high-confidence and high-quality set of protein interaction partners which are enriched for functions in RNA metabolism and ribosomal proteins. The screen is highly specific, maximizing the number of true positives, as demonstrated by the validation of 81% (47/58) of the RNA-independent interactors through known functions and interactions. Importantly, we minimize the number of indirect interaction partners through use of a nuclease-based digestion to eliminate RNA. We describe eleven new interactors, including proteins involved in splicing which is an as-yet unknown role for DDX6. We validated and characterized in more detail the interaction of DDX6 with Nuclear fragile X mental retardation-interacting protein 2 (NUFIP2) and with two previously uncharacterized proteins, FAM195A and FAM195B (here referred to as granulin-1 and granulin-2, or GRAN1 and GRAN2). We show that NUFIP2, GRAN1, and GRAN2 are not P-body components, but re-localize to stress granules upon exposure to stress, suggesting a function in translation repression in the cellular stress response. Using a complementary analysis that resolved DDX6’s multiple complex memberships, we further validated these interaction partners and the presence of splicing factors. As DDX6 also interacts with the E3 SUMO ligase TIF1β, we tested for and observed a significant enrichment of sumoylation amongst DDX6’s interaction partners. Our results represent the most comprehensive screen for direct interaction partners of a key regulator of RNA life cycle and localization, highlighting new stress granule components and possible DDX6 functions—many of which are likely conserved across eukaryotes.

High-resolution profiling of protein occupancy on polyadenylated RNA transcripts

A key prerequisite to understand how gene regulatory processes are controlled by the interplay of RNA-binding proteins and ribonucleoprotein complexes with RNAs is the generation of comprehensive high-resolution maps of protein-RNA interactions. Recent advances in next-generation sequencing technology accelerated the development of various crosslinking and immunoprecipitation (CLIP) approaches to broadly identify RNA regions contacted by RNA-binding proteins. However these methods only consider single RNA-binding proteins and their contact sites, irrespective of the overall cis-regulatory sequence space contacted by other RNA interacting factors. Here we describe the application of protein occupancy profiling, a novel approach that globally displays the RNA contact sites of the poly(A)+ RNA-bound proteome. Protein occupancy profiling enables the generation of transcriptome-wide maps of protein-RNA interactions on polyadenylated transcripts and narrows the sequence search space for transcript regions involved in cis-regulation of gene expression in response to internal or external stimuli, altered cellular programs or disease.

New insights into the cellular temporal response to proteostatic stress

Maintaining a healthy proteome involves all layers of gene expression regulation. By quantifying temporal changes of the transcriptome, translatome, proteome, and RNA-protein interactome in cervical cancer cells, we systematically characterize the molecular landscape in response to proteostatic challenges. We identify shared and specific responses to misfolded proteins and to oxidative stress, two conditions that are tightly linked. We reveal new aspects of the unfolded protein response, including many genes that escape global translation shutdown. A subset of these genes supports rerouting of energy production in the mitochondria. We also find that many genes change at multiple levels, in either the same or opposing directions, and at different time points. We highlight a variety of putative regulatory pathways, including the stress-dependent alternative splicing of aminoacyl-tRNA synthetases, and protein-RNA binding within the 3’ untranslated region of molecular chaperones. These results illustrate the potential of this information-rich resource.

Quantifying multi-layered expression regulation in response to stress of the endoplasmic reticulum

The mammalian response to endoplasmic reticulum (ER) stress dynamically affects all layers of gene expression regulation. We quantified transcript and protein abundance along with footprints of ribosomes and non-ribosomal proteins for thousands of genes in cervical cancer cells responding to treatment with tunicamycin or hydrogen peroxide over an eight hour time course. We identify shared and stress-specific significant regulatory events at the transcriptional and post-transcriptional level and at different phases of the experiment. ER stress regulators increase transcription and translation at different times supporting an adaptive response. ER stress also induces translation of genes from serine biosynthesis and one-carbon metabolism indicating a shift in energy production. Discordant regulation of DNA repair genes suggests transcriptional priming in which delayed translation fine-tunes the early change in the transcriptome. Finally, case studies on stress-dependent alternative splicing and protein-mRNA binding demonstrate the ability of this resource to generate hypotheses for new regulatory mechanisms.