Mudskipper detects combinatorial RNA binding protein interactions in multiplexed CLIP data.

The uncovering of protein-RNA interactions enables a deeper understanding of RNA processing. Recent multiplexed crosslinking and immunoprecipitation (CLIP) technologies such as antibody-barcoded eCLIP (ABC) dramatically increase the throughput of mapping RNA binding protein (RBP) binding sites. However, multiplex CLIP datasets are multivariate, and each RBP suffers non-uniform signal-to-noise ratio. To address this, we developed Mudskipper, a versatile computational suite comprising two components: a Dirichlet multinomial mixture model to account for the multivariate nature of ABC datasets and a softmasking approach that identifies and removes non-specific protein-RNA interactions in RBPs with low signal-to-noise ratio. Mudskipper demonstrates superior precision and recall over existing tools on multiplex datasets and supports analysis of repetitive elements and small non-coding RNAs. Our findings unravel splicing outcomes and variant-associated disruptions, enabling higher-throughput investigations into diseases and regulation mediated by RBPs.

Decoding protein-RNA interactions using CLIP-based methodologies.

Protein-RNA interactions are central to all RNA processing events, with pivotal roles in the regulation of gene expression and cellular functions. Dysregulation of these interactions has been increasingly linked to the pathogenesis of human diseases. High-throughput approaches to identify RNA-binding proteins and their binding sites on RNA - in particular, ultraviolet crosslinking followed by immunoprecipitation (CLIP) - have helped to map the RNA interactome, yielding transcriptome-wide protein-RNA atlases that have contributed to key mechanistic insights into gene expression and gene-regulatory networks. Here, we review these recent advances, explore the effects of cellular context on RNA binding, and discuss how these insights are shaping our understanding of cellular biology. We also review the potential therapeutic applications arising from new knowledge of protein-RNA interactions.

Protocol to process crosslinking and immunoprecipitation data into annotated binding sites.

Here, we present a protocol for using Skipper, a pipeline designed to process crosslinking and immunoprecipitation (CLIP) data into annotated binding sites. We describe steps for partitioning annotated transcript regions and fitting data to a beta-binomial model to call windows of enriched binding. From raw CLIP data, we detail how users can map reproducible RNA-binding sites to call enriched windows and perform downstream analysis. This protocol supports optional customizations for different use cases. For complete details on the use and execution of this protocol, please refer to Boyle et al.1.

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.

HydRA: Deep-learning models for predicting RNA-binding capacity from protein interaction association context and protein sequence.

RNA-binding proteins (RBPs) control RNA metabolism to orchestrate gene expression and, when dysfunctional, underlie human diseases. Proteome-wide discovery efforts predict thousands of RBP candidates, many of which lack canonical RNA-binding domains (RBDs). Here, we present a hybrid ensemble RBP classifier (HydRA), which leverages information from both intermolecular protein interactions and internal protein sequence patterns to predict RNA-binding capacity with unparalleled specificity and sensitivity using support vector machines (SVMs), convolutional neural networks (CNNs), and Transformer-based protein language models. Occlusion mapping by HydRA robustly detects known RBDs and predicts hundreds of uncharacterized RNA-binding associated domains. Enhanced CLIP (eCLIP) for HydRA-predicted RBP candidates reveals transcriptome-wide RNA targets and confirms RNA-binding activity for HydRA-predicted RNA-binding associated domains. HydRA accelerates construction of a comprehensive RBP catalog and expands the diversity of RNA-binding associated domains.

Large-scale map of RNA binding protein interactomes across the mRNA life-cycle.

Messenger RNAs (mRNAs) interact with RNA-binding proteins (RBPs) in diverse ribonucleoprotein complexes (RNPs) during distinct life-cycle stages for their processing and maturation. While substantial attention has focused on understanding RNA regulation by assigning proteins, particularly RBPs, to specific RNA substrates, there has been considerably less exploration leveraging protein-protein interaction (PPI) methodologies to identify and study the role of proteins in mRNA life-cycle stages. To address this gap, we generated an RNA-aware RBP-centric PPI map across the mRNA life-cycle by immunopurification (IP-MS) of ~100 endogenous RBPs across the life-cycle in the presence or absence of RNase, augmented by size exclusion chromatography (SEC-MS). Aside from confirming 8,700 known and discovering 20,359 novel interactions between 1125 proteins, we determined that 73% of our IP interactions are regulated by the presence of RNA. Our PPI data enables us to link proteins to life-cycle stage functions, highlighting that nearly half of the proteins participate in at least two distinct stages. We show that one of the most highly interconnected proteins, ERH, engages in multiple RNA processes, including via interactions with nuclear speckles and the mRNA export machinery. We also demonstrate that the spliceosomal protein SNRNP200 participates in distinct stress granule-associated RNPs and occupies different RNA target regions in the cytoplasm during stress. Our comprehensive RBP-focused PPI network is a novel resource for identifying multi-stage RBPs and exploring RBP complexes in RNA maturation.

Proteomic discovery of chemical probes that perturb protein complexes in human cells.

Most human proteins lack chemical probes, and several large-scale and generalizable small-molecule binding assays have been introduced to address this problem. How compounds discovered in such "binding-first" assays affect protein function, nonetheless, often remains unclear. Here, we describe a "function-first" proteomic strategy that uses size exclusion chromatography (SEC) to assess the global impact of electrophilic compounds on protein complexes in human cells. Integrating the SEC data with cysteine-directed activity-based protein profiling identifies changes in protein-protein interactions that are caused by site-specific liganding events, including the stereoselective engagement of cysteines in PSME1 and SF3B1 that disrupt the PA28 proteasome regulatory complex and stabilize a dynamic state of the spliceosome, respectively. Our findings thus show how multidimensional proteomic analysis of focused libraries of electrophilic compounds can expedite the discovery of chemical probes with site-specific functional effects on protein complexes in human cells.

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.

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.

RNA Binding Proteins as Pioneer Determinants of Infection: Protective, Proviral, or Both?

As the first intracellular host factors that directly interact with the genomes of RNA viruses, RNA binding proteins (RBPs) have a profound impact on the outcome of an infection. Recent discoveries brought about by new methodologies have led to an unprecedented ability to peer into the earliest events between viral RNA and the RBPs that act upon them. These discoveries have sparked a re-evaluation of current paradigms surrounding RBPs and post-transcriptional gene regulation. Here, we highlight questions that have bloomed from the implementation of these novel approaches. Canonical RBPs can impact the fates of both cellular and viral RNA during infection, sometimes in conflicting ways. Noncanonical RBPs, some of which were first characterized via interactions with viral RNA, may encompass physiological roles beyond viral pathogenesis. We discuss how these RBPs might discriminate between an RNA of either cellular or viral origin and thus exert either pro- or antiviral effects-which is a particular challenge as viruses contain mechanisms to mimic molecular features of cellular RNA.

ELAVL1 primarily couples mRNA stability with the 3' UTRs of interferon-stimulated genes.

Upon pathogen detection, the innate immune system triggers signaling events leading to upregulation of pro-inflammatory and anti-microbial mRNA transcripts. RNA-binding proteins (RBPs) interact with these critical mRNAs and regulate their fates at the post-transcriptional level. One such RBP is ELAVL1. Although significant progress has been made in understanding how embryonic lethal vision-like protein 1 (ELAVL1) regulates mRNAs, its target repertoire and binding distribution within an immunological context remain poorly understood. We overlap four high-throughput approaches to define its context-dependent targets and determine its regulatory impact during immune activation. ELAVL1 transitions from binding overwhelmingly intronic sites to 3' UTR sites upon immune stimulation of cells, binding previously and newly expressed mRNAs. We find that ELAVL1 mediates the RNA stability of genes that regulate pathways essential to pathogen sensing and cytokine production. Our findings reveal the importance of examining RBP regulatory impact under dynamic transcriptomic events to understand their post-transcriptional regulatory roles within specific biological circuitries.

Viral crosslinking and solid-phase purification enables discovery of ribonucleoprotein complexes on incoming RNA virus genomes.

The initial interactions between incoming, pre-replicated virion RNA and host protein factors are important in infection and immunity. Yet currently there are no methods to study these crucial events. We established VIR-CLASP (VIRal Cross-Linking And Solid-phase Purification) to identify the primary viral RNA-host protein interactions. First, host cells are infected with 4-thiouridine (4SU)-labeled RNA viruses and irradiated with 365 nm light to crosslink 4SU-labeled viral genomes and interacting proteins from host or virus. The crosslinked RNA binding proteins (RBPs) are purified by solid-phase reversible immobilization (SPRI) beads with protein-denaturing buffers, and then identified by proteomics. With VIR-CLASP, only the incoming virion RNAs are labeled with 4SU, so crosslinking events specifically occur between proteins and pre-replicated virion RNA. Since solid-phase purification under protein-denaturing conditions, rather than sequence-specific nucleic acid purification, is used to pull-down total RNA and crosslinked RBPs, this method facilitates investigation of potentially all RNA viruses, regardless of RNA sequence. Preparation of 4SU-labeled virus takes ∼7 days and VIR-CLASP takes 1 day.

Discovery of Widespread Host Protein Interactions with the Pre-replicated Genome of CHIKV Using VIR-CLASP.

The primary interactions between incoming viral RNA genomes and host proteins are crucial to infection and immunity. Until now, the ability to study these events was lacking. We developed viral cross-linking and solid-phase purification (VIR-CLASP) to characterize the earliest interactions between viral RNA and cellular proteins. We investigated the infection of human cells using Chikungunya virus (CHIKV) and influenza A virus and identified hundreds of direct RNA-protein interactions. Here, we explore the biological impact of three protein classes that bind CHIKV RNA within minutes of infection. We find CHIKV RNA binds and hijacks the lipid-modifying enzyme fatty acid synthase (FASN) for pro-viral activity. We show that CHIKV genomes are N6-methyladenosine modified, and YTHDF1 binds and suppresses CHIKV replication. Finally, we find that the innate immune DNA sensor IFI16 associates with CHIKV RNA, reducing viral replication and maturation. Our findings have direct applicability to the investigation of potentially all RNA viruses.

Publisher Correction: Small molecule inhibition of cGAS reduces interferon expression in primary macrophages from autoimmune mice.

The previously published version of this Article contained errors in Fig. 6. In panel h the units of the x axis were incorrectly given as mM and should have been given as µM. Also, the IC50s for RU.365, RU.332 and RU.521 within panel h were incorrectly given as mM and should have been given as µM. These errors have been corrected in both the PDF and HTML versions of the Article.

Small molecule inhibition of cGAS reduces interferon expression in primary macrophages from autoimmune mice.

Cyclic GMP-AMP synthase is essential for innate immunity against infection and cellular damage, serving as a sensor of DNA from pathogens or mislocalized self-DNA. Upon binding double-stranded DNA, cyclic GMP-AMP synthase synthesizes a cyclic dinucleotide that initiates an inflammatory cellular response. Mouse studies that recapitulate causative mutations in the autoimmune disease Aicardi-Goutières syndrome demonstrate that ablating the cyclic GMP-AMP synthase gene abolishes the deleterious phenotype. Here, we report the discovery of a class of cyclic GMP-AMP synthase inhibitors identified by a high-throughput screen. These compounds possess defined structure-activity relationships and we present crystal structures of cyclic GMP-AMP synthase, double-stranded DNA, and inhibitors within the enzymatic active site. We find that a chemically improved member, RU.521, is active and selective in cellular assays of cyclic GMP-AMP synthase-mediated signaling and reduces constitutive expression of interferon in macrophages from a mouse model of Aicardi-Goutières syndrome. RU.521 will be useful toward understanding the biological roles of cyclic GMP-AMP synthase and can serve as a molecular scaffold for development of future autoimmune therapies.Upon DNA binding cyclic GMP-AMP synthase (cGAS) produces a cyclic dinucleotide, which leads to the upregulation of inflammatory genes. Here the authors develop small molecule cGAS inhibitors, functionally characterize them and present the inhibitor and DNA bound cGAS crystal structures, which will facilitate drug development.

Parsing the Interferon Transcriptional Network and Its Disease Associations.

Type 1 interferon (IFN) is a key mediator of organismal responses to pathogens, eliciting prototypical "interferon signature genes" that encode antiviral and inflammatory mediators. For a global view of IFN signatures and regulatory pathways, we performed gene expression and chromatin analyses of the IFN-induced response across a range of immunocyte lineages. These distinguished ISGs by cell-type specificity, kinetics, and sensitivity to tonic IFN and revealed underlying changes in chromatin configuration. We combined 1,398 human and mouse datasets to computationally infer ISG modules and their regulators, validated by genetic analysis in both species. Some ISGs are controlled by Stat1/2 and Irf9 and the ISRE DNA motif, but others appeared dependent on non-canonical factors. This regulatory framework helped to interpret JAK1 blockade pharmacology, different clusters being affected under tonic or IFN-stimulated conditions, and the IFN signatures previously associated with human diseases, revealing unrecognized subtleties in disease footprints, as affected by human ancestry.

Gene expression during the generation and activation of mouse neutrophils: implication of novel functional and regulatory pathways.

As part of the Immunological Genome Project (ImmGen), gene expression was determined in unstimulated (circulating) mouse neutrophils and three populations of neutrophils activated in vivo, with comparison among these populations and to other leukocytes. Activation conditions included serum-transfer arthritis (mediated by immune complexes), thioglycollate-induced peritonitis, and uric acid-induced peritonitis. Neutrophils expressed fewer genes than any other leukocyte population studied in ImmGen, and down-regulation of genes related to translation was particularly striking. However, genes with expression relatively specific to neutrophils were also identified, particularly three genes of unknown function: Stfa2l1, Mrgpr2a and Mrgpr2b. Comparison of genes up-regulated in activated neutrophils led to several novel findings: increased expression of genes related to synthesis and use of glutathione and of genes related to uptake and metabolism of modified lipoproteins, particularly in neutrophils elicited by thioglycollate; increased expression of genes for transcription factors in the Nr4a family, only in neutrophils elicited by serum-transfer arthritis; and increased expression of genes important in synthesis of prostaglandins and response to leukotrienes, particularly in neutrophils elicited by uric acid. Up-regulation of genes related to apoptosis, response to microbial products, NFkB family members and their regulators, and MHC class II expression was also seen, in agreement with previous studies. A regulatory model developed from the ImmGen data was used to infer regulatory genes involved in the changes in gene expression during neutrophil activation. Among 64, mostly novel, regulatory genes predicted to influence these changes in gene expression, Irf5 was shown to be important for optimal secretion of IL-10, IP-10, MIP-1α, MIP-1ß, and TNF-α by mouse neutrophils in vitro after stimulation through TLR9. This data-set and its analysis using the ImmGen regulatory model provide a basis for additional hypothesis-based research on the importance of changes in gene expression in neutrophils in different conditions.

The transcriptional landscape of αß T cell differentiation.

The differentiation of αßT cells from thymic precursors is a complex process essential for adaptive immunity. Here we exploited the breadth of expression data sets from the Immunological Genome Project to analyze how the differentiation of thymic precursors gives rise to mature T cell transcriptomes. We found that early T cell commitment was driven by unexpectedly gradual changes. In contrast, transit through the CD4(+)CD8(+) stage involved a global shutdown of housekeeping genes that is rare among cells of the immune system and correlated tightly with expression of the transcription factor c-Myc. Selection driven by major histocompatibility complex (MHC) molecules promoted a large-scale transcriptional reactivation. We identified distinct signatures that marked cells destined for positive selection versus apoptotic deletion. Differences in the expression of unexpectedly few genes accompanied commitment to the CD4(+) or CD8(+) lineage, a similarity that carried through to peripheral T cells and their activation, demonstrated by mass cytometry phosphoproteomics. The transcripts newly identified as encoding candidate mediators of key transitions help define the 'known unknowns' of thymocyte differentiation.

Identification of transcriptional regulators in the mouse immune system.

The differentiation of hematopoietic stem cells into cells of the immune system has been studied extensively in mammals, but the transcriptional circuitry that controls it is still only partially understood. Here, the Immunological Genome Project gene-expression profiles across mouse immune lineages allowed us to systematically analyze these circuits. To analyze this data set we developed Ontogenet, an algorithm for reconstructing lineage-specific regulation from gene-expression profiles across lineages. Using Ontogenet, we found differentiation stage-specific regulators of mouse hematopoiesis and identified many known hematopoietic regulators and 175 previously unknown candidate regulators, as well as their target genes and the cell types in which they act. Among the previously unknown regulators, we emphasize the role of ETV5 in the differentiation of γδ T cells. As the transcriptional programs of human and mouse cells are highly conserved, it is likely that many lessons learned from the mouse model apply to humans.

Conservation and divergence in the transcriptional programs of the human and mouse immune systems.

Much of the knowledge about cell differentiation and function in the immune system has come from studies in mice, but the relevance to human immunology, diseases, and therapy has been challenged, perhaps more from anecdotal than comprehensive evidence. To this end, we compare two large compendia of transcriptional profiles of human and mouse immune cell types. Global transcription profiles are conserved between corresponding cell lineages. The expression patterns of most orthologous genes are conserved, particularly for lineage-specific genes. However, several hundred genes show clearly divergent expression across the examined cell lineages, and among them, 169 genes did so even with highly stringent criteria. Finally, regulatory mechanisms--reflected by regulators' differential expression or enriched cis-elements--are conserved between the species but to a lower degree, suggesting that distinct regulation may underlie some of the conserved transcriptional responses.