Expanded palette of RNA base editors for comprehensive RBP-RNA interactome studies.

RNA binding proteins (RBPs) are key regulators of RNA processing and cellular function. Technologies to discover RNA targets of RBPs such as TRIBE (targets of RNA binding proteins identified by editing) and STAMP (surveying targets by APOBEC1 mediated profiling) utilize fusions of RNA base-editors (rBEs) to RBPs to circumvent the limitations of immunoprecipitation (CLIP)-based methods that require enzymatic digestion and large amounts of input material. To broaden the repertoire of rBEs suitable for editing-based RBP-RNA interaction studies, we have devised experimental and computational assays in a framework called PRINTER (protein-RNA interaction-based triaging of enzymes that edit RNA) to assess over thirty A-to-I and C-to-U rBEs, allowing us to identify rBEs that expand the characterization of binding patterns for both sequence-specific and broad-binding RBPs. We also propose specific rBEs suitable for dual-RBP applications. We show that the choice between single or multiple rBEs to fuse with a given RBP or pair of RBPs hinges on the editing biases of the rBEs and the binding preferences of the RBPs themselves. We believe our study streamlines and enhances the selection of rBEs for the next generation of RBP-RNA target discovery.

Expanded palette of RNA base editors for comprehensive RBP-RNA interactome studies.

RNA binding proteins (RBPs) are key regulators of RNA processing and cellular function. Technologies to discover RNA targets of RBPs such as TRIBE (targets of RNA binding proteins identified by editing) and STAMP (surveying targets by APOBEC1 mediated profiling) utilize fusions of RNA base-editors (rBEs) to RBPs to circumvent the limitations of immunoprecipitation (CLIP)-based methods that require enzymatic digestion and large amounts of input material. To broaden the repertoire of rBEs suitable for editing-based RBP-RNA interaction studies, we have devised experimental and computational assays in a framework called PRINTER (protein-RNA interaction-based triaging of enzymes that edit RNA) to assess over thirty A-to-I and C-to-U rBEs, allowing us to identify rBEs that expand the characterization of binding patterns for both sequence-specific and broad-binding RBPs. We also propose specific rBEs suitable for dual-RBP applications. We show that the choice between single or multiple rBEs to fuse with a given RBP or pair of RBPs hinges on the editing biases of the rBEs and the binding preferences of the RBPs themselves. We believe our study streamlines and enhances the selection of rBEs for the next generation of RBP-RNA target discovery.

A super-enhancer-regulated RNA-binding protein cascade drives pancreatic cancer.

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy in need of new therapeutic options. Using unbiased analyses of super-enhancers (SEs) as sentinels of core genes involved in cell-specific function, here we uncover a druggable SE-mediated RNA-binding protein (RBP) cascade that supports PDAC growth through enhanced mRNA translation. This cascade is driven by a SE associated with the RBP heterogeneous nuclear ribonucleoprotein F, which stabilizes protein arginine methyltransferase 1 (PRMT1) to, in turn, control the translational mediator ubiquitin-associated protein 2-like. All three of these genes and the regulatory SE are essential for PDAC growth and coordinately regulated by the Myc oncogene. In line with this, modulation of the RBP network by PRMT1 inhibition reveals a unique vulnerability in Myc-high PDAC patient organoids and markedly reduces tumor growth in male mice. Our study highlights a functional link between epigenetic regulation and mRNA translation and identifies components that comprise unexpected therapeutic targets for PDAC.

RANK ligand converts the NCoR/HDAC3 co-repressor to a PGC1ß- and RNA-dependent co-activator of osteoclast gene expression.

The nuclear receptor co-repressor (NCoR) complex mediates transcriptional repression dependent on histone deacetylation by histone deacetylase 3 (HDAC3) as a component of the complex. Unexpectedly, we found that signaling by the receptor activator of nuclear factor κB (RANK) converts the NCoR/HDAC3 co-repressor complex to a co-activator of AP-1 and NF-κB target genes that are required for mouse osteoclast differentiation. Accordingly, the dominant function of NCoR/HDAC3 complexes in response to RANK signaling is to activate, rather than repress, gene expression. Mechanistically, RANK signaling promotes RNA-dependent interaction of the transcriptional co-activator PGC1ß with the NCoR/HDAC3 complex, resulting in the activation of PGC1ß and inhibition of HDAC3 activity for acetylated histone H3. Non-coding RNAs Dancr and Rnu12, which are associated with altered human bone homeostasis, promote NCoR/HDAC3 complex assembly and are necessary for RANKL-induced osteoclast differentiation in vitro. These findings may be prototypic for signal-dependent functions of NCoR in other biological contexts.

Skipper analysis of eCLIP datasets enables sensitive detection of constrained translation factor binding sites.

Technology for crosslinking and immunoprecipitation (CLIP) followed by sequencing (CLIP-seq) has identified the transcriptomic targets of hundreds of RNA-binding proteins in cells. To increase the power of existing and future CLIP-seq datasets, we introduce Skipper, an end-to-end workflow that converts unprocessed reads into annotated binding sites using an improved statistical framework. Compared with existing methods, Skipper on average calls 210%-320% more transcriptomic binding sites and sometimes >1,000% more sites, providing deeper insight into post-transcriptional gene regulation. Skipper also calls binding to annotated repetitive elements and identifies bound elements for 99% of enhanced CLIP experiments. We perform nine translation factor enhanced CLIPs and apply Skipper to learn determinants of translation factor occupancy, including transcript region, sequence, and subcellular localization. Furthermore, we observe depletion of genetic variation in occupied sites and nominate transcripts subject to selective constraint because of translation factor occupancy. Skipper offers fast, easy, customizable, and state-of-the-art analysis of CLIP-seq data.

MECP2-related pathways are dysregulated in a cortical organoid model of myotonic dystrophy.

Myotonic dystrophy type 1 (DM1) is a multisystem, autosomal-dominant inherited disorder caused by CTG microsatellite repeat expansions (MREs) in the 3' untranslated region of the dystrophia myotonica-protein kinase (DMPK) gene. Despite its prominence as the most common adult-onset muscular dystrophy, patients with congenital to juvenile-onset forms of DM1 can present with debilitating neurocognitive symptoms along the autism spectrum, characteristic of possible in utero cortical defects. However, the molecular mechanism by which CTG MREs lead to these developmental central nervous system (CNS) manifestations is unknown. Here, we showed that CUG foci found early in the maturation of three-dimensional (3D) cortical organoids from DM1 patient-derived induced pluripotent stem cells (iPSCs) cause hyperphosphorylation of CUGBP Elav-like family member 2 (CELF2) protein. Integrative single-cell RNA sequencing and enhanced cross-linking and immunoprecipitation (eCLIP) analysis revealed that reduced CELF2 protein-RNA substrate interactions results in misregulation of genes critical for excitatory synaptic signaling in glutamatergic neurons, including key components of the methyl-CpG binding protein 2 (MECP2) pathway. Comparisons to MECP2(y/-) cortical organoids revealed convergent molecular and cellular defects such as glutamate toxicity and neuronal loss. Our findings provide evidence suggesting that early-onset DM1 might involve neurodevelopmental disorder-associated pathways and identify N-methyl-d-aspartic acid (NMDA) antagonists as potential treatment avenues for neuronal defects in DM1.

Crosstalk between CRISPR-Cas9 and the human transcriptome.

CRISPR-Cas9 expression independent of its cognate synthetic guide RNA (gRNA) causes widespread genomic DNA damage in human cells. To investigate whether Cas9 can interact with endogenous human RNA transcripts independent of its guide, we perform eCLIP (enhanced CLIP) of Cas9 in human cells and find that Cas9 reproducibly interacts with hundreds of endogenous human RNA transcripts. This association can be partially explained by a model built on gRNA secondary structure and sequence. Critically, transcriptome-wide Cas9 binding sites do not appear to correlate with published genome-wide Cas9 DNA binding or cut-site loci under gRNA co-expression. However, even under gRNA co-expression low-affinity Cas9-human RNA interactions (which we term CRISPR crosstalk) do correlate with published elevated transcriptome-wide RNA editing. Our findings do not support the hypothesis that human RNAs can broadly guide Cas9 to bind and cleave human genomic DNA, but they illustrate a cellular and RNA impact likely inherent to CRISPR-Cas systems.

Discovery and functional interrogation of SARS-CoV-2 protein-RNA interactions.

The COVID-19 pandemic is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The betacoronvirus has a positive sense RNA genome which encodes for several RNA binding proteins. Here, we use enhanced crosslinking and immunoprecipitation to investigate SARS-CoV-2 protein interactions with viral and host RNAs in authentic virus-infected cells. SARS-CoV-2 proteins, NSP8, NSP12, and nucleocapsid display distinct preferences to specific regions in the RNA viral genome, providing evidence for their shared and separate roles in replication, transcription, and viral packaging. SARS-CoV-2 proteins expressed in human lung epithelial cells bind to 4773 unique host coding RNAs. Nine SARS-CoV-2 proteins upregulate target gene expression, including NSP12 and ORF9c, whose RNA substrates are associated with pathways in protein N-linked glycosylation ER processing and mitochondrial processes. Furthermore, siRNA knockdown of host genes targeted by viral proteins in human lung organoid cells identify potential antiviral host targets across different SARS-CoV-2 variants. Conversely, NSP9 inhibits host gene expression by blocking mRNA export and dampens cytokine productions, including interleukin-1α/ß. Our viral protein-RNA interactome provides a catalog of potential therapeutic targets and offers insight into the etiology of COVID-19 as a safeguard against future pandemics.

Transcriptome-wide identification of RNA-binding protein binding sites using seCLIP-seq.

Discovery of interaction sites between RNA-binding proteins (RBPs) and their RNA targets plays a critical role in enabling our understanding of how these RBPs control RNA processing and regulation. Cross-linking and immunoprecipitation (CLIP) provides a generalizable, transcriptome-wide method by which RBP/RNA complexes are purified and sequenced to identify sites of intermolecular contact. By simplifying technical challenges in prior CLIP methods and incorporating the generation of and quantitative comparison against size-matched input controls, the single-end enhanced CLIP (seCLIP) protocol allows for the profiling of these interactions with high resolution, efficiency and scalability. Here, we present a step-by-step guide to the seCLIP method, detailing critical steps and offering insights regarding troubleshooting and expected results while carrying out the ~4-d protocol. Furthermore, we describe a comprehensive bioinformatics pipeline that offers users the tools necessary to process two replicate datasets and identify reproducible and significant peaks for an RBP of interest in ~2 d.

Discovery and functional interrogation of SARS-CoV-2 protein-RNA interactions.

The COVID-19 pandemic is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The betacoronvirus has a positive sense RNA genome which encodes for several RNA binding proteins. Here, we use enhanced crosslinking and immunoprecipitation to investigate SARS-CoV-2 protein interactions with viral and host RNAs in authentic virus-infected cells. SARS-CoV-2 proteins, NSP8, NSP12, and nucleocapsid display distinct preferences to specific regions in the RNA viral genome, providing evidence for their shared and separate roles in replication, transcription, and viral packaging. SARS-CoV-2 proteins expressed in human lung epithelial cells bind to 4773 unique host coding RNAs. Nine SARS-CoV-2 proteins upregulate target gene expression, including NSP12 and ORF9c, whose RNA substrates are associated with pathways in protein N-linked glycosylation ER processing and mitochondrial processes. Furthermore, siRNA knockdown of host genes targeted by viral proteins in human lung organoid cells identify potential antiviral host targets across different SARS-CoV-2 variants. Conversely, NSP9 inhibits host gene expression by blocking mRNA export and dampens cytokine productions, including interleukin-1α/ß. Our viral protein-RNA interactome provides a catalog of potential therapeutic targets and offers insight into the etiology of COVID-19 as a safeguard against future pandemics.

Context-dependent functional compensation between Ythdf m6A reader proteins.

The N6-methyladenosine (m6A) modification is the most prevalent post-transcriptional mRNA modification, regulating mRNA decay and splicing. It plays a major role during normal development, differentiation, and disease progression. The modification is regulated by a set of writer, eraser, and reader proteins. The YTH domain family of proteins consists of three homologous m6A-binding proteins, Ythdf1, Ythdf2, and Ythdf3, which were suggested to have different cellular functions. However, their sequence similarity and their tendency to bind the same targets suggest that they may have overlapping roles. We systematically knocked out (KO) the Mettl3 writer, each of the Ythdf readers, and the three readers together (triple-KO). We then estimated the effect in vivo in mouse gametogenesis, postnatal viability, and in vitro in mouse embryonic stem cells (mESCs). In gametogenesis, Mettl3-KO severity is increased as the deletion occurs earlier in the process, and Ythdf2 has a dominant role that cannot be compensated by Ythdf1 or Ythdf3, due to differences in readers' expression pattern across different cell types, both in quantity and in spatial location. Knocking out the three readers together and systematically testing viable offspring genotypes revealed a redundancy in the readers' role during early development that is Ythdf1/2/3 gene dosage-dependent. Finally, in mESCs there is compensation between the three Ythdf reader proteins, since the resistance to differentiate and the significant effect on mRNA decay occur only in the triple-KO cells and not in the single KOs. Thus, we suggest a new model for the Ythdf readers function, in which there is profound dosage-dependent redundancy when all three readers are equivalently coexpressed in the same cell types.

Principles of RNA processing from analysis of enhanced CLIP maps for 150 RNA binding proteins.

BACKGROUND: A critical step in uncovering rules of RNA processing is to study the in vivo regulatory networks of RNA binding proteins (RBPs). Crosslinking and immunoprecipitation (CLIP) methods enable mapping RBP targets transcriptome-wide, but methodological differences present challenges to large-scale analysis across datasets. The development of enhanced CLIP (eCLIP) enabled the mapping of targets for 150 RBPs in K562 and HepG2, creating a unique resource of RBP interactomes profiled with a standardized methodology in the same cell types. RESULTS: Our analysis of 223 eCLIP datasets reveals a range of binding modalities, including highly resolved positioning around splicing signals and mRNA untranslated regions that associate with distinct RBP functions. Quantification of enrichment for repetitive and abundant multicopy elements reveals 70% of RBPs have enrichment for non-mRNA element classes, enables identification of novel ribosomal RNA processing factors and sites, and suggests that association with retrotransposable elements reflects multiple RBP mechanisms of action. Analysis of spliceosomal RBPs indicates that eCLIP resolves AQR association after intronic lariat formation, enabling identification of branch points with single-nucleotide resolution, and provides genome-wide validation for a branch point-based scanning model for 3' splice site recognition. Finally, we show that eCLIP peak co-occurrences across RBPs enable the discovery of novel co-interacting RBPs. CONCLUSIONS: This work reveals novel insights into RNA biology by integrated analysis of eCLIP profiling of 150 RBPs with distinct functions. Further, our quantification of both mRNA and other element association will enable further research to identify novel roles of RBPs in regulating RNA processing.