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.

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.

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.

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.

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.

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.

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.

The 46° Halo and Its Arcs

Ice crystals in the form of right hexagonal prisms have faces that form 90 degrees prisms. Light rays were traced through these prism faces by computer calculation, and the light patterns that would be produced in the sky for a particular distribution of crystal orientations were simulated. Crystals with random orientations produce a 46 degrees halo. Hexagonal plate crystals with nearly horizontal end faces produce circumzenithal and circumhorizontal arcs. Hexagonal column crystals with horizontal axes produce supralateral and infralateral arcs. Plate crystals spinning about a horizontal axis that is a face diagonal of the crystal produce a series of arcs touching the 46 degrees halo. Each of these effects was simulated for several elevations of the sun.

Structure of the steroid-binding site of human placental estradiol-17beta-dehydrogenase

We recently demonstrated that human placental estradiol-17beta-dehydrogenase possesses a histidyl residue in the catalytic region of the active site by affinity-labeling studies with 16alpha-bromoacetoxyestradiol-3-methyl ether. We now report the synthesis of 12beta-bromoacetoxy-4-estrene-3,17-dione and its use in affinity labeling of the enzyme. The steroid was synthesized by incubation of 4-estrene-3,17-dione with Colletotrichum gloesporioides. The product was recrystallized from ethanol and structure assured by IR and NMR spectroscopy. The steroid is a substrate, which indicates that it binds at the active site. When the enzyme is incubated with a 150-fold molar excess of 12beta-bromoacetoxy-4-estrene-3,17-dione in potassium phosphate buffer at pH 7.0, the enzyme is inactivated in a time-dependent, irreversible manner. Inactivation follows pseudo-first-order kinetics with a half life of 18 hours. Analysis of a hydrolysate of the enzyme after inactivation with 12beta-bromo[2'-3H]acetoxy-4-estrene-3,17-dione reveals tritiated 1-, 3-, and 1,3-dicarboxymethylhistidine. The affinity labeling of a histidyl enzyme residue by both 16alpha- and 12beta-bromoacetoxy steroids localizes that residue near the point of catalysis and suggests that it may participate in the catalytic event.

The diurnal pattern of 17 --oestradiol in pregnancy

Plasma unconjugated 17β-oestradiol was measured at 8 a.m., 4 p.m. and 10 p. m. in eleven normal females during the third trimester of pregnancy. The late evening values were significantly lower than those observed at 8 a. m. Marked differences in plasma levels were noted for some individuals at the three periods of sampling. A slight diurnal variation in plasma 17β-oestradiol appears compatible with the present concept of oestrogen biosynthesis in late pregnancy.