Rbfox proteins regulate alternative mRNA splicing through evolutionarily conserved RNA bridges

MicroRNA mechanisms instructing Purkinje cell specification

MicroRNAs (miRNAs) are critical for brain development; however, if, when, and how miRNAs drive neuronal subtype specification remains poorly understood. To address this, we engineered technologies with vastly improved spatiotemporal resolution that allow the dissection of cell-type-specific miRNA-target networks. Fast and reversible miRNA loss of function showed that miRNAs are necessary for Purkinje cell (PC) differentiation, which previously appeared to be miRNA independent, and identified distinct critical miRNA windows for dendritogenesis and climbing fiber synaptogenesis, structural features defining PC identity. Using new mouse models that enable miRNA-target network mapping in rare cell types, we uncovered PC-specific post-transcriptional programs. Manipulation of these programs revealed that the PC-enriched miR-206 and targets Shank3, Prag1, En2, and Vash1, which are uniquely repressed in PCs, are critical regulators of PC-specific dendritogenesis and synaptogenesis, with miR-206 knockdown and target overexpression partially phenocopying miRNA loss of function. Our results suggest that gene expression regulation by miRNAs, beyond transcription, is critical for neuronal subtype specification.

Peak analysis of cell-free RNA finds recurrently protected narrow regions with clinical potential

BACKGROUND

Cell-free RNAs (cfRNAs) can be detected in biofluids and have emerged as valuable disease biomarkers. Accurate identification of the fragmented cfRNA signals, especially those originating from pathological cells, is crucial for understanding their biological functions and clinical value. However, many challenges still need to be addressed for their application, including developing specific analysis methods and translating cfRNA fragments with biological support into clinical applications.

RESULTS

We present cfPeak, a novel method combining statistics and machine learning models to detect the fragmented cfRNA signals effectively. When test in real and artificial cfRNA sequencing (cfRNA-seq) data, cfPeak shows an improved performance compared with other applicable methods. We reveal that narrow cfRNA peaks preferentially overlap with protein binding sites, vesicle-sorting sites, structural sites, and novel small non-coding RNAs (sncRNAs). When applied in clinical cohorts, cfPeak identified cfRNA peaks in patients' plasma that enable cancer detection and are informative of cancer types and metastasis.

CONCLUSIONS

Our study fills the gap in the current small cfRNA-seq analysis at fragment-scale and builds a bridge to the scientific discovery in cfRNA fragmentomics. We demonstrate the significance of finding low abundant tissue-derived signals in small cfRNA and prove the feasibility for application in liquid biopsy.

RBFOX1 Regulates Calcium Signaling and Enhances SERCA2 Translation

RBFOX1 is an RNA-binding protein that regulates alternative splicing and RNA processing in the neurons, skeletal muscle, and heart. We intended to define the role of RBFOX1 in regulating calcium homeostasis to maintain normal cardiac function. We generated cardiomyocyte-specific Rbfox1 gene-deletion mice (cKO). The cardiomyocyte-specific deletion of RBFOX1 was confirmed by Western blotting and immunohistochemistry. The cKO mice showed mild hypertrophy and depressed cardiac function under homeostatic conditions, which did not deteriorate with age. Pressure overload by trans-aortic constriction (TAC) caused exaggerated cardiac hypertrophy and accelerated heart failure in cKO compared with wild-type mice. We performed Western blotting to assess the expression of important Ca2+-handling proteins, which showed alterations in the phosphorylation of PLN and CAMKII and decreased expression of SERCA2. We measured the Ca2+ dynamics and noted significantly delayed Ca2+ reuptake into the sarcoplasmic reticulum. Importantly, the decrease in SERCA2 expression was not due to reduced mRNA expression or altered splicing. To assess the possibility of the post-transcriptional regulation of SERCA2 expression by RBFOX1, we performed RNA immunoprecipitation (RIP), which showed the binding of RBFOX1 protein to Serca2 mRNA, which was confirmed in luciferase assays with the Serca2a 3'-untranslated region fused to luciferase. Finally, we performed a puromycin incorporation experiment, which showed that RBFOX1 enhances SERCA2 protein translation. Our results show that RBFOX1 plays a crucial role in regulating the expression of Ca2+-handling genes to maintain normal cardiac function. We show an important post-transcriptional role of RBFOX1 in regulating SERCA2 expression.

Unveiling the regulatory role of GRP7 in ABA signal-mediated mRNA translation efficiency regulation

Abscisic acid (ABA) is a crucial phytohormone involved in plant growth and stress responses. While the transcriptional regulation triggered by ABA is well-documented, its effects on translational regulation have been less studied. Through Ribo-seq and RNA-seq analyses, we find that ABA treatment not only influences gene expression at the mRNA level but also significantly impacts mRNA translation efficiency (TE) in Arabidopsis thaliana. ABA inhibits global mRNA translation via its core signaling pathway, which includes ABA receptors, protein phosphatase 2Cs (PP2Cs), and SNF1-related protein kinase 2 s (SnRK2s). Upon ABA treatment, Glycine-rich RNA-binding proteins 7 and 8 (GRP7&8) protein levels decrease due to both reduced mRNA level and decreased TE, which diminishes their association with polysomes and leads to a global decline in mRNA TE. The absence of GRP7&8 results in a global impairment of ABA-regulated translational changes, linking ABA signaling to GRP7-dependent modulation of mRNA translation. The regulation of GRP7 on TE relies significantly on its direct binding to target mRNAs. Moreover, mRNA translation efficiency under drought stress is partially dependent on the ABA-GRP7&8 pathways. Collectively, our study reveals GRP7's role downstream of SnRK2s in mediating translation regulation in ABA signaling, offering a model for ABA-triggered multi-route regulation of environmental adaptation.

Integrated multi-omic characterizations of the synapse reveal RNA processing factors and ubiquitin ligases associated with neurodevelopmental disorders

The molecular composition of the excitatory synapse is incompletely defined due to its dynamic nature across developmental stages and neuronal populations. To address this gap, we apply proteomic mass spectrometry to characterize the synapse in multiple biological models, including the fetal human brain and human induced pluripotent stem cell (hiPSC)-derived neurons. To prioritize the identified proteins, we develop an orthogonal multi-omic screen of genomic, transcriptomic, interactomic, and structural data. This data-driven framework identifies proteins with key molecular features intrinsic to the synapse, including characteristic patterns of biophysical interactions and cross-tissue expression. The multi-omic analysis captures synaptic proteins across developmental stages and experimental systems, including 493 synaptic candidates supported by proteomics. We further investigate three such proteins that are associated with neurodevelopmental disorders-Cullin 3 (CUL3), DEAD-box helicase 3 X-linked (DDX3X), and Y-box binding protein-1 (YBX1)-by mapping their networks of physically interacting synapse proteins or transcripts. Our study demonstrates the potential of an integrated multi-omic approach to more comprehensively resolve the synaptic architecture.

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We find frequent splicing order differences between alleles and uncover significant single-nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This includes SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also have a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, show the most allele-specific splicing orders, which frequently co-occur with allele-specific AS, APA, or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Clusters of deep intronic RbFox motifs embedded in large assembly of splicing regulators sequences regulate alternative splicing

The RbFox RNA binding proteins regulate alternative splicing of genes governing mammalian development and organ function. They bind to the RNA sequence (U)GCAUG with high affinity but also non-canonical secondary motifs in a concentration dependent manner. However, the hierarchical requirement of RbFox motifs, which are widespread in the genome, is still unclear. Here we show that deep intronic, tightly clustered RbFox1 motifs cooperate and are important regulators of alternative exons splicing. Bioinformatic analysis revealed that (U)GCAUG-clusters are widely present in both mouse and human genes and are embedded in sequences binding the large assembly of splicing regulators (LASR). Integrative data analysis from eCLIP and RNAseq experiments showed a global increase in RNA isoform modulation of genes with Rbfox1 eCLIP-peaks associated with these clusters. Experimentally, by employing recombineering mutagenesis in a bacterial artificial chromosome containing the NTrk2 mouse region subjected to alternative splicing we showed that tightly clustered (U)GCAUG motifs in the middle of 50 Kb introns are necessary for RbFox1 regulation of NTrk2 gene isoforms expression. Moreover, clustered (U)GCAUG-motifs promote the recruitment of RbFox1 proteins to form a Rbfox1/LASR complex required for splicing. These data suggest that clustered, distal intronic Rbfox-binding motifs embedded in LASR binding sequences are important determinants of RbFox1 function in the mammalian genome and provide a target for identification of pathogenic mutations.

Influenza virus antagonizes self sensing by RIG-I to enhance viral replication

Innate immune sensors must finely distinguish pathogens from the host to mount a response only during infection. RIG-I is cytoplasmic sensor that surveils for foreign RNAs. When activated, RIG-I triggers a broad antiviral response that is a major regulator of RNA virus infection. Here were show that RIG-I not only bound viral RNAs, but was activated by host RNAs to amplify the antiviral state. These were primarily non-coding RNAs transcribed by RNA polymerase III. They were benign under normal conditions but became immunogenic during influenza virus infection where they signaled via RIG-I to suppress viral replication. This same class of RNAs was bound by influenza virus nucleoprotein (NP), which normally functions to encapsidate the viral genome. NP interacted with RIG-I and antagonized sensing of self RNAs to counter innate immune responses. Overall, these results demonstrate that self sensing is strategically deployed by the cell to amplify the antiviral response and reveal a newly identified viral countermeasure that disrupts RIG-I activation by host RNAs.

The Rbfox1/LASR complex controls alternative pre-mRNA splicing by recognition of multipart RNA regulatory modules

The Rbfox proteins regulate alternative pre-mRNA splicing by binding to the RNA element GCAUG. In the nucleus, most of Rbfox is bound to the large assembly of splicing regulators (LASR), a complex of RNA-binding proteins that recognize additional RNA motifs. However, it remains unclear how the different subunits of the Rbfox/LASR complex act together to bind RNA and regulate splicing. We used a nuclease protection assay to map the transcriptome-wide footprints of Rbfox1/LASR on nascent cellular RNA. In addition to GCAUG, Rbfox1/LASR binds RNA motifs for LASR subunits hnRNPs M, H/F, and C and Matrin3. These elements are often arranged in tandem, forming multipart modules of RNA motifs. To distinguish contact sites of Rbfox1 from the LASR subunits, we analyzed a mutant Rbfox1(F125A) that has lost RNA binding but remains associated with LASR. Rbfox1(F125A)/LASR complexes no longer interact with GCAUG but retain binding to RNA elements for LASR. Splicing analyses reveal that in addition to activating exons through adjacent GCAUG elements, Rbfox can also stimulate exons near binding sites for LASR subunits. Minigene experiments demonstrate that these diverse elements produce a combined regulatory effect on a target exon. These findings illuminate how a complex of RNA-binding proteins can decode combinatorial splicing regulatory signals by recognizing groups of tandem RNA elements.

Mapping snoRNA-target RNA interactions in an RNA-binding protein-dependent manner with chimeric eCLIP

BACKGROUND

Small nucleolar RNAs (snoRNAs) are non-coding RNAs that function in ribosome and spliceosome biogenesis, primarily by guiding modifying enzymes to specific sites on ribosomal RNA (rRNA) and spliceosomal RNA (snRNA). However, many orphan snoRNAs remain uncharacterized, with unidentified or unvalidated targets, and studies on additional snoRNA-associated proteins are limited.

RESULTS

We adapted an enhanced chimeric eCLIP approach to comprehensively profile snoRNA-target RNA interactions using both core and accessory snoRNA-binding proteins as baits. Using core snoRNA-binding proteins, we confirmed most annotated snoRNA-rRNA and snoRNA-snRNA interactions in mouse and human cell lines and called novel, high-confidence interactions for orphan snoRNAs. While some of these interactions result in chemical modification, others may have modification-independent functions. We showed that snoRNA ribonucleoprotein complexes containing certain accessory proteins, like WDR43 and NOLC1, enriched for specific subsets of snoRNA-target RNA interactions with distinct roles in ribosome and spliceosome biogenesis. Notably, we discovered that SNORD89 guides 2'-O-methylation at two neighboring sites in U2 snRNA that fine-tune splice site recognition.

CONCLUSIONS

Chimeric eCLIP of snoRNA-associating proteins enables a comprehensive framework for studying snoRNA-target interactions in an RNA-binding protein-dependent manner, revealing novel interactions and regulatory roles in RNA biogenesis.

Aberrant splicing in Huntington's disease accompanies disrupted TDP-43 activity and altered m6A RNA modification

Huntington's disease (HD) is caused by a CAG repeat expansion in the HTT gene, leading to altered gene expression. However, the mechanisms leading to disrupted RNA processing in HD remain unclear. Here we identify TDP-43 and the N6-methyladenosine (m6A) writer protein METTL3 to be upstream regulators of exon skipping in multiple HD systems. Disrupted nuclear localization of TDP-43 and cytoplasmic accumulation of phosphorylated TDP-43 occurs in HD mouse and human brains, with TDP-43 also co-localizing with HTT nuclear aggregate-like bodies distinct from mutant HTT inclusions. The binding of TDP-43 onto RNAs encoding HD-associated differentially expressed and aberrantly spliced genes is decreased. Finally, m6A RNA modification is reduced on RNAs abnormally expressed in the striatum of HD R6/2 mouse brain, including at clustered sites adjacent to TDP-43 binding sites. Our evidence supports TDP-43 loss of function coupled with altered m6A modification as a mechanism underlying alternative splicing in HD.

ECD functions as a novel RNA-binding protein to regulate mRNA splicing

The human ecdysoneless protein (ECD) plays an essential role in the regulation of cell cycle and cell survival. ECD has been implicated in RNA splicing through its association with the protein components of splicing complex. Here, using electrophoretic mobility shift assay and mutational analysis, we demonstrate that ECD directly binds to RNA through its N-terminal region, specifically using amino acids 135-148. Using enhanced CLIP-seq analyses in human cells, we identified a large repertoire of mRNAs bound to ECD. RNA-seq analyses revealed that ECD depletion in cells leads to widespread RNA splicing aberrations associated with alterations in gene expression. Significantly, we demonstrate that ECD mediates mRNA splicing by directly binding to RNA sequences located near splicing sites. Mechanistically, we demonstrate that ECD directly binds to U5 small nuclear RNA (snRNA), and this interaction is critical for maintaining the expression of key protein components of U5 small nuclear protein (snRNP) complex. Notably, RNA binding defective mutant of ECD fails to rescue downregulated levels of U5 snRNP components or cell proliferation block induced by ECD knockout. Collectively, we provide compelling evidence that ECD regulates RNA splicing by directly associating with RNAs, and the RNA binding activity of ECD is essential for its function.

LARP6 regulates the mRNA translation of fibrogenic genes in liver fibrosis

Metabolic syndrome and excessive alcohol consumption result in liver injury and fibrosis, which is characterized by increased collagen production by activated Hepatic Stellate Cells (HSCs). LARP6, an RNA-binding protein, was shown to facilitate collagen production. However, LARP6 expression and functionality as a regulator of fibrosis development in a disease relevant model remains elusive. By using snRNA-sequencing, we show that LARP6 is upregulated mainly in HSCs of liver fibrosis patients. Moreover, LARP6 knockdown in human HSCs suppresses fibrogenic gene expression. By integrating eCLIP analysis and ribosome profiling in HSCs, we show that LARP6 interacts with mature mRNAs comprising over 300 genes, including RNA structural elements within COL1A1 , COL1A2 , and COL3A1 to regulate mRNA expression and translation. Furthermore, LARP6 knockdown in HSC attenuates fibrosis development in human liver spheroids. Altogether, our results suggest that targeting LARP6 in human HSCs may provide new strategies for anti-fibrotic therapy.

Highlights

LARP6 is upregulated in liver fibrosis, mainly in HSCs.LARP6 knockdown in human HSCs reduces liver fibrosis development.Of the hundreds of gene targets, LARP6 interacts most with collagen mRNAs.LARP6 regulates mRNA translation via interaction with 5'UTRs.

Transcriptome and translatome profiling of Col-0 and grp7grp8 under ABA treatment in Arabidopsis

Abscisic acid (ABA) is a crucial phytohormone that regulates plant growth and stress responses. While substantial knowledge exists about transcriptional regulation, the molecular mechanisms underlying ABA-triggered translational regulation remain unclear. Recent advances in deep sequencing of ribosome footprints (Ribo-seq) enable the mapping and quantification of mRNA translation efficiency. Additionally, RNA-binding proteins (RBPs) play essential roles in translational regulation by interacting with target RNA molecules, making the identification of binding sites via UV crosslinking and immunoprecipitation (CLIP) critical for understanding RBP function. Glycine-rich RNA-binding proteins (GRPs), a prominent class of RBPs in plants, are responsive to ABA. In this study, RNA-seq and Ribo-seq analyses were conducted on 3-day-old Col-0 and grp7grp8 seedlings of Arabidopsis thaliana, treated with either ABA or mock solutions. These analyses facilitated deep sequencing of total mRNA and mRNA fragments protected by translating ribosomes. Additionally, CLIP-seq analysis of pGRP7::GRP7-GFP grp7-1 identified RNA bound by GRP7. This multi-omics dataset allows for a comprehensive investigation of the plant's response to ABA from various perspectives, providing a significant resource for studying ABA-regulated mRNA translation efficiency.

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.

Dual modes of ZFC3H1 confer selectivity in nuclear RNA sorting

The export and degradation pathways compete to sort nuclear RNAs, yet the default pathway remains unclear. Sorting of mature RNAs to degradation, facilitated by the exosome co-factor poly(A) exosome targeting (PAXT), is particularly challenging for their resemblance to mRNAs intended for translation. Here, we unveil that ZFC3H1, a core PAXT component, is co-transcriptionally loaded onto the first exon/intron of RNA precursors (pre-RNAs). Interestingly, this initial loading does not lead to pre-RNA degradation, as ZFC3H1 adopts a "closed" conformation, effectively blocking exosome recruitment. As processing progresses, RNA fate can be reshaped. Longer RNAs with more exons are allowed for nuclear export. By contrast, short RNAs with fewer exons preferentially recruit transient PAXT components ZC3H3 and RBM26/27 to the 3' end, triggering ZFC3H1 "opening" and subsequent exosomal degradation. Together, the decoupled loading and activation of ZFC3H1 pre-configures RNA fate for decay while still allowing a switch to nuclear export, depending on mature RNA features.

Alternative splicing of PBRM1 mediates resistance to PD-1 blockade therapy in renal cancer

Alternative pre-mRNA splicing (AS) is a biological process that results in proteomic diversity. However, implications of AS alterations in cancer remain poorly understood. Herein, we performed a comprehensive AS analysis in cancer driver gene transcripts across fifteen cancer types and found global alterations in inclusion rates of the PBAF SWI/SNF chromatin remodeling complex subunit Polybromo 1 (PBRM1) exon 27 (E27) in most types of cancer tissues compared with those in normal tissues. Further analysis confirmed that PBRM1 E27 is excluded by the direct binding of RBFOX2 to intronic UGCAUG elements. In addition, the E27-included PBRM1 isoform upregulated PD-L1 expression via enhanced PBAF complex recruitment to the PD-L1 promoter. PBRM1 wild-type patients with clear cell renal cell carcinoma were resistant to PD-1 blockade therapy when they expressed low RBFOX2 mRNA levels. Overall, our study suggests targeting of RBFOX2-mediated AS of PBRM1 as a potential therapeutic strategy for immune checkpoint blockade.

The YTHDF Proteins Shape the Brain Gene Signatures of Alzheimer's Disease

The gene signatures of Alzheimer's Disease (AD) brains reflect an output of a complex interplay of genetic, epigenetic, epi-transcriptomic, and post-transcriptional regulations. To identify the most significant factor that shapes the AD brain signature, we developed a machine learning model (DEcode-tree) to integrate cellular and molecular factors explaining differential gene expression in AD. Our model indicates that YTHDF proteins, the canonical readers of N6-methyladenosine RNA modification (m6A), are the most influential predictors of the AD brain signature. We then show that protein modules containing YTHDFs are downregulated in human AD brains, and knocking out YTHDFs in iPSC-derived neural cells recapitulates the AD brain gene signature in vitro . Furthermore, eCLIP-seq analysis revealed that YTHDF proteins influence AD signatures through both m6A-dependent and independent pathways. These results indicate the central role of YTHDF proteins in shaping the gene signature of AD brains.

An imbalance between proliferation and differentiation underlies the development of microRNA-defective pineoblastoma

Mutations in the microRNA processing genes DICER1 and DROSHA drive several cancers that resemble embryonic progenitors. To understand how microRNAs regulate tumorigenesis, we ablated Drosha or Dicer1 in the developing pineal gland to emulate the pathogenesis of pineoblastoma, a brain tumor that resembles undifferentiated precursors of the pineal gland. Accordingly, these mice develop pineal tumors marked by loss of microRNAs, including the let-7/miR-98-5p family, and de-repression of microRNA target genes. Pineal tumors driven by loss of Drosha or Dicer1 mimic tumors driven by Rb1 loss, as they exhibit upregulation of S-phase genes and homeobox transcription factors that regulate pineal development. Blocking proliferation of these tumors facilitates expression of pinealocyte maturation markers, with a concomitant reduction in embryonic markers. Select embryonic markers remain elevated, however, as the microRNAs that normally repress these target genes remain absent. One such microRNA target gene is the oncofetal transcription factor Plagl2, which regulates expression of pro-growth genes, and inhibiting their signaling impairs tumor growth. Thus, we demonstrate that tumors driven by loss of microRNA processing may be therapeutically targeted by inhibiting downstream drivers of proliferation.

Identification of RNA structures and their roles in RNA functions

The development of high-throughput RNA structure profiling methods in the past decade has greatly facilitated our ability to map and characterize different aspects of RNA structures transcriptome-wide in cell populations, single cells and single molecules. The resulting high-resolution data have provided insights into the static and dynamic nature of RNA structures, revealing their complexity as they perform their respective functions in the cell. In this Review, we discuss recent technical advances in the determination of RNA structures, and the roles of RNA structures in RNA biogenesis and functions, including in transcription, processing, translation, degradation, localization and RNA structure-dependent condensates. We also discuss the current understanding of how RNA structures could guide drug design for treating genetic diseases and battling pathogenic viruses, and highlight existing challenges and future directions in RNA structure research.

Mapping snoRNA-target RNA interactions in an RNA binding protein-dependent manner with chimeric eCLIP

Small nucleolar RNAs (snoRNAs) are non-coding RNAs that function in ribosome and spliceosome biogenesis, primarily by guiding modifying enzymes to specific sites on ribosomal RNA (rRNA) and spliceosomal RNA (snRNA). However, many orphan snoRNAs remain uncharacterized, with unidentified or unvalidated targets, and studies on additional snoRNA-associated proteins are limited. We adapted an enhanced chimeric eCLIP approach to comprehensively profile snoRNA-target RNA interactions using both core and accessory snoRNA binding proteins as baits. Using core snoRNA binding proteins, we confirmed most annotated snoRNA-rRNA and snoRNA-snRNA interactions in mouse and human cell lines and called novel, high-confidence interactions for orphan snoRNAs. While some of these interactions result in chemical modification, others may have modification-independent functions. We then showed that snoRNA ribonucleoprotein complexes containing certain accessory proteins, like WDR43 and NOLC1, enriched for specific subsets of snoRNA-target RNA interactions with distinct roles in ribosome and spliceosome biogenesis. Notably, we discovered that SNORD89 guides 2'-O-methylation at two neighboring sites in U2 snRNA that are important for activating splicing, but also appear to ensure imperfect splicing for a subset of near-constitutive exons. Thus, chimeric eCLIP of snoRNA-associating proteins enables a comprehensive framework for studying snoRNA-target interactions in an RNA binding protein-dependent manner, revealing novel interactions and regulatory roles in RNA biogenesis.

The RNA-binding Selectivity of the RGG/RG Motifs of hnRNP U is Abolished by Elements Within the C-terminal Intrinsically Disordered Region

The abundant nuclear protein hnRNP U interacts with a broad array of RNAs along with DNA and protein to regulate nuclear chromatin architecture. The RNA-binding activity is achieved via a disordered ∼100 residue C-terminal RNA-binding domain (RBD) containing two distinct RGG/RG motifs. Although the RNA-binding capabilities of RGG/RG motifs have been widely reported, less is known about hnRNP U's RNA-binding selectivity. Furthermore, while it is well established that hnRNP U binds numerous nuclear RNAs, it remains unknown whether it selectively recognizes sequence or structural motifs in target RNAs. To address this question, we performed equilibrium binding assays using fluorescence anisotropy (FA) and electrophoretic mobility shift assays (EMSAs) to quantitatively assess the ability of human hnRNP U RBD to interact with segments of cellular RNAs identified from eCLIP data. These RNAs often, but not exclusively, contain poly-uridine or 5'-AGGGAG sequence motifs. Detailed binding analysis of several target RNAs reveal that the hnRNP U RBD binds RNA in a promiscuous manner with high affinity for a broad range of structured RNAs, but with little preference for any distinct sequence motif. In contrast, the isolated RGG/RG of hnRNP U motif exhibits a strong preference for G-quadruplexes, similar to that observed for other RGG motif bearing peptides. These data reveal that the hnRNP U RBD attenuates the RNA binding selectivity of its core RGG motifs to achieve an extensive RNA interactome. We propose that a critical role of RGG/RG motifs in RNA biology is to alter binding affinity or selectivity of adjacent RNA-binding domains.

Large-scale evaluation of the ability of RNA-binding proteins to activate exon inclusion

RNA-binding proteins (RBPs) modulate alternative splicing outcomes to determine isoform expression and cellular survival. To identify RBPs that directly drive alternative exon inclusion, we developed tethered function luciferase-based splicing reporters that provide rapid, scalable and robust readouts of exon inclusion changes and used these to evaluate 718 human RBPs. We performed enhanced cross-linking immunoprecipitation, RNA sequencing and affinity purification-mass spectrometry to investigate a subset of candidates with no prior association with splicing. Integrative analysis of these assays indicates surprising roles for TRNAU1AP, SCAF8 and RTCA in the modulation of hundreds of endogenous splicing events. We also leveraged our tethering assays and top candidates to identify potent and compact exon inclusion activation domains for splicing modulation applications. Using these identified domains, we engineered programmable fusion proteins that outperform current artificial splicing factors at manipulating inclusion of reporter and endogenous exons. This tethering approach characterizes the ability of RBPs to induce exon inclusion and yields new molecular parts for programmable splicing control.

Genetic regulation of nascent RNA maturation revealed by direct RNA nanopore sequencing

Quantitative trait loci analyses have revealed an important role for genetic variants in regulating alternative splicing (AS) and alternative cleavage and polyadenylation (APA) in humans. Yet, these studies are generally performed with mature mRNA, so they report on the outcome rather than the processes of RNA maturation and thus may overlook how variants directly modulate pre-mRNA processing. The order in which the many introns of a human gene are removed can substantially influence AS, while nascent RNA polyadenylation can affect RNA stability and decay. However, how splicing order and poly(A) tail length are regulated by genetic variation has never been explored. Here, we used direct RNA nanopore sequencing to investigate allele-specific pre-mRNA maturation in 12 human lymphoblastoid cell lines. We found frequent splicing order differences between alleles and uncovered significant single nucleotide polymorphism (SNP)-splicing order associations in 17 genes. This included SNPs located in or near splice sites as well as more distal intronic and exonic SNPs. Moreover, several genes showed allele-specific poly(A) tail lengths, many of which also had a skewed allelic abundance ratio. HLA class I transcripts, which encode proteins that play an essential role in antigen presentation, showed the most allele-specific splicing orders, which frequently co-occurred with allele-specific AS, APA or poly(A) tail length differences. Together, our results expose new layers of genetic regulation of pre-mRNA maturation and highlight the power of long-read RNA sequencing for allele-specific analyses.

Full-Length Transcriptome Construction and Systematic Characterization of Virulence Factor-Associated Isoforms in Vairimorpha (Nosema) Ceranae

Vairimorpha (Nosema) ceranae is a single-cellular fungus that obligately infects the midgut epithelial cells of adult honeybees, causing bee microsporidiosis and jeopardizing bee health and production. This work aims to construct the full-length transcriptome of V. ceranae and conduct a relevant investigation using PacBio single-molecule real-time (SMRT) sequencing technology. Following PacBio SMRT sequencing, 41,950 circular consensus (CCS) were generated, and 25,068 full-length non-chimeric (FLNC) reads were then detected. After polishing, 4387 high-quality, full-length transcripts were gained. There are 778, 2083, 1202, 1559, 1457, 1232, 1702, and 3896 full-length transcripts that could be annotated to COG, GO, KEGG, KOG, Pfam, Swiss-Prot, eggNOG, and Nr databases, respectively. Additionally, 11 alternative splicing (AS) events occurred in 6 genes were identified, including 1 alternative 5' splice-site and 10 intron retention. The structures of 225 annotated genes in the V. ceranae reference genome were optimized, of which 29 genes were extended at both 5' UTR and 3' UTR, while 90 and 106 genes were, respectively, extended at the 5' UTR as well as 3' UTR. Furthermore, a total of 29 high-confidence lncRNAs were obtained, including 12 sense-lncRNAs, 10 lincRNAs, and 7 antisense-lncRNAs. Taken together, the high-quality, full-length transcriptome of V. ceranae was constructed and annotated, the structures of annotated genes in the V. ceranae reference genome were improved, and abundant new genes, transcripts, and lncRNAs were discovered. Findings from this current work offer a valuable resource and a crucial foundation for molecular and omics research on V. ceranae.

The Rbfox1/LASR complex controls alternative pre-mRNA splicing by recognition of multi-part RNA regulatory modules

The Rbfox proteins regulate alternative pre-mRNA splicing by binding to the RNA element GCAUG. In the nucleus, most of Rbfox is bound to LASR, a complex of RNA-binding proteins that recognize additional RNA motifs. However, it remains unclear how the different subunits of the Rbfox/LASR complex act together to bind RNA and regulate splicing. We used a nuclease-protection assay to map the transcriptome-wide footprints of Rbfox1/LASR on nascent cellular RNA. In addition to GCAUG, Rbfox1/LASR binds RNA containing motifs for LASR subunits hnRNPs M, H/F, C, and Matrin3. These elements are often arranged in tandem, forming multi-part modules of RNA motifs. To distinguish contact sites of Rbfox1 from the LASR subunits, we analyzed a mutant Rbfox1(F125A) that has lost RNA binding but remains associated with LASR. Rbfox1(F125A)/LASR complexes no longer interact with GCAUG but retain binding to RNA elements for LASR. Splicing analyses reveal that in addition to activating exons through adjacent GCAUG elements, Rbfox can also stimulate exons near binding sites for LASR subunits. Mini-gene experiments demonstrate that these diverse elements produce a combined regulatory effect on a target exon. These findings illuminate how a complex of RNA-binding proteins can decode combinatorial splicing regulatory signals by recognizing groups of tandem RNA elements.

ePRINT: exonuclease assisted mapping of protein-RNA interactions

RNA-binding proteins (RBPs) regulate key aspects of RNA processing including alternative splicing, mRNA degradation and localization by physically binding RNA molecules. Current methods to map these interactions, such as CLIP, rely on purifying single proteins at a time. Our new method, ePRINT, maps RBP-RNA interaction networks on a global scale without purifying individual RBPs. ePRINT uses exoribonuclease XRN1 to precisely map the 5' end of the RBP binding site and uncovers direct and indirect targets of an RBP of interest. Importantly, ePRINT can also uncover RBPs that are differentially activated between cell fate transitions, including neural progenitor differentiation into neurons.

MAPP unravels frequent co-regulation of splicing and polyadenylation by RNA-binding proteins and their dysregulation in cancer

Maturation of eukaryotic pre-mRNAs via splicing and polyadenylation is modulated across cell types and conditions by a variety of RNA-binding proteins (RBPs). Although there exist over 1,500 RBPs in human cells, their binding motifs and functions still remain to be elucidated, especially in the complex environment of tissues and in the context of diseases. To overcome the lack of methods for the systematic and automated detection of sequence motif-guided pre-mRNA processing regulation from RNA sequencing (RNA-Seq) data we have developed MAPP (Motif Activity on Pre-mRNA Processing). Applying MAPP to RBP knock-down experiments reveals that many RBPs regulate both splicing and polyadenylation of nascent transcripts by acting on similar sequence motifs. MAPP not only infers these sequence motifs, but also unravels the position-dependent impact of the RBPs on pre-mRNA processing. Interestingly, all investigated RBPs that act on both splicing and 3' end processing exhibit a consistently repressive or activating effect on both processes, providing a first glimpse on the underlying mechanism. Applying MAPP to normal and malignant brain tissue samples unveils that the motifs bound by the PTBP1 and RBFOX RBPs coordinately drive the oncogenic splicing program active in glioblastomas demonstrating that MAPP paves the way for characterizing pre-mRNA processing regulators under physiological and pathological conditions.

A phage nucleus-associated RNA-binding protein is required for jumbo phage infection

Large-genome bacteriophages (jumbo phages) of the proposed family Chimalliviridae assemble a nucleus-like compartment bounded by a protein shell that protects the replicating phage genome from host-encoded restriction enzymes and DNA-targeting CRISPR-Cas nucleases. While the nuclear shell provides broad protection against host nucleases, it necessitates transport of mRNA out of the nucleus-like compartment for translation by host ribosomes, and transport of specific proteins into the nucleus-like compartment to support DNA replication and mRNA transcription. Here, we identify a conserved phage nuclear shell-associated protein that we term Chimallin C (ChmC), which adopts a nucleic acid-binding fold, binds RNA with high affinity in vitro, and binds phage mRNAs in infected cells. ChmC also forms phase-separated condensates with RNA in vitro. Targeted knockdown of ChmC using mRNA-targeting dCas13d results in accumulation of phage-encoded mRNAs in the phage nucleus, reduces phage protein production, and compromises virion assembly. Taken together, our data show that the conserved ChmC protein plays crucial roles in the viral life cycle, potentially by facilitating phage mRNA translocation through the nuclear shell to promote protein production and virion development.

CRISPR-dCas13d-based deep screening of proximal and distal splicing-regulatory elements

Pre-mRNA splicing, a key process in gene expression, can be therapeutically modulated using various drug modalities, including antisense oligonucleotides (ASOs). However, determining promising targets is hampered by the challenge of systematically mapping splicing-regulatory elements (SREs) in their native sequence context. Here, we use the catalytically inactive CRISPR-RfxCas13d RNA-targeting system (dCas13d/gRNA) as a programmable platform to bind SREs and modulate splicing by competing against endogenous splicing factors. SpliceRUSH, a high-throughput screening method, was developed to map SREs in any gene of interest using a lentivirus gRNA library that tiles the genetic region, including distal intronic sequences. When applied to SMN2, a therapeutic target for spinal muscular atrophy, SpliceRUSH robustly identifies not only known SREs but also a previously unknown distal intronic SRE, which can be targeted to alter exon 7 splicing using either dCas13d/gRNA or ASOs. This technology enables a deeper understanding of splicing regulation with applications for RNA-based drug discovery.

Transcriptome regulation by PARP13 in basal and antiviral states in human cells

The RNA-binding protein PARP13 is a primary factor in the innate antiviral response, which suppresses translation and drives decay of bound viral and host RNA. PARP13 interacts with many proteins encoded by interferon-stimulated genes (ISG) to activate antiviral pathways including co-translational addition of ISG15, or ISGylation. We performed enhanced crosslinking immunoprecipitation (eCLIP) and RNA-seq in human cells to investigate PARP13's role in transcriptome regulation for both basal and antiviral states. We find that the antiviral response shifts PARP13 target localization, but not its binding preferences, and that PARP13 supports the expression of ISGylation-related genes, including PARP13's cofactor, TRIM25. PARP13 associates with TRIM25 via RNA-protein interactions, and we elucidate a transcriptome-wide periodicity of PARP13 binding around TRIM25. Taken together, our study implicates PARP13 in creating and maintaining a cellular environment poised for an antiviral response through limiting PARP13 translation, regulating access to distinct mRNA pools, and elevating ISGylation machinery expression.

Plasma Circular RNAs as Biomarkers for Breast Cancer

Breast cancer (BC) is currently the most common neoplasm, the second leading cause of cancer death in women worldwide, and is a major health problem. The discovery of new biomarkers is crucial to improve our knowledge of breast cancer and strengthen our clinical approaches to diagnosis, prognosis, and follow-up. In recent decades, there has been increasing interest in circulating RNA (circRNA) as modulators of gene expression involved in tumor development and progression. The study of circulating circRNAs (ccircRNAs) in plasma may provide new non-invasive diagnostic, prognostic, and predictive biomarkers for BC. This review describes the latest findings on BC-associated ccircRNAs in plasma and their clinical utility. Several ccircRNAs in plasma have shown great potential as BC biomarkers, especially from a diagnostic point of view. Mechanistically, most of the reported BC-associated ccircRNAs are involved in the regulation of cell survival, proliferation, and invasion, mainly via MAPK/AKT signaling pathways. However, the study of circRNAs is a relatively new area of research, and a larger number of studies will be crucial to confirm their potential as plasma biomarkers and to understand their involvement in BC.

RNA binding proteins in cardiovascular development and disease

Congenital heart disease (CHD) is the most common birth defect affecting>1.35 million newborn babies worldwide. CHD can lead to prenatal, neonatal, postnatal lethality or life-long cardiac complications. RNA binding protein (RBP) mutations or variants are emerging as contributors to CHDs. RBPs are wizards of gene regulation and are major contributors to mRNA and protein landscape. However, not much is known about RBPs in the developing heart and their contributions to CHD. In this chapter, we will discuss our current knowledge about specific RBPs implicated in CHDs. We are in an exciting era to study RBPs using the currently available and highly successful RNA-based therapies and methodologies. Understanding how RBPs shape the developing heart will unveil their contributions to CHD. Identifying their target RNAs in the embryonic heart will ultimately lead to RNA-based treatments for congenital heart disease.

Pooled CRISPR Screening Identifies P-Bodies as Repressors of Cancer Epithelial-Mesenchymal Transition

Epithelial-mesenchymal transition (EMT) is a fundamental cellular process frequently hijacked by cancer cells to promote tumor progression, especially metastasis. EMT is orchestrated by a complex molecular network acting at different layers of gene regulation. In addition to transcriptional regulation, posttranscriptional mechanisms may also play a role in EMT. Here, we performed a pooled CRISPR screen analyzing the influence of 1,547 RNA-binding proteins on cell motility in colon cancer cells and identified multiple core components of P-bodies (PB) as negative modulators of cancer cell migration. Further experiments demonstrated that PB depletion by silencing DDX6 or EDC4 could activate hallmarks of EMT thereby enhancing cell migration in vitro as well as metastasis formation in vivo. Integrative multiomics analysis revealed that PBs could repress the translation of the EMT driver gene HMGA2, which contributed to PB-meditated regulation of EMT. This mechanism is conserved in other cancer types. Furthermore, endoplasmic reticulum stress was an intrinsic signal that induced PB disassembly and translational derepression of HMGA2. Taken together, this study has identified a function of PBs in the regulation of EMT in cancer.

SIGNIFICANCE

Systematic investigation of the influence of posttranscriptional regulation on cancer cell motility established a connection between P-body-mediated translational control and EMT, which could be therapeutically exploited to attenuate metastasis formation.

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

An alternative splicing signature defines the basal-like phenotype and predicts worse clinical outcome in pancreatic cancer

Pancreatic ductal adenocarcinoma (PDAC) is characterized by extremely poor prognosis. PDAC presents with molecularly distinct subtypes, with the basal-like one being associated with enhanced chemoresistance. Splicing dysregulation contributes to PDAC; however, its involvement in subtype specification remains elusive. Herein, we uncover a subtype-specific splicing signature associated with prognosis in PDAC and the splicing factor Quaking (QKI) as a determinant of the basal-like signature. Single-cell sequencing analyses highlight QKI as a marker of the basal-like phenotype. QKI represses splicing events associated with the classical subtype while promoting basal-like events associated with shorter survival. QKI favors a plastic, quasi-mesenchymal phenotype that supports migration and chemoresistance in PDAC organoids and cell lines, and its expression is elevated in high-grade primary tumors and metastatic lesions. These studies identify a splicing signature that defines PDAC subtypes and indicate that QKI promotes an undifferentiated, plastic phenotype, which renders PDAC cells chemoresistant and adaptable to environmental changes.

hnRNP A1 dysfunction alters RNA splicing and drives neurodegeneration in multiple sclerosis (MS)

Neurodegeneration is the primary driver of disease progression in multiple sclerosis (MS) resulting in permanent disability, creating an urgent need to discover its underlying mechanisms. Herein, we establish that dysfunction of the RNA binding protein heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) results in differential of binding to RNA targets causing alternative RNA splicing, which contributes to neurodegeneration in MS and its models. Using RNAseq of MS brains, we discovered differential expression and aberrant splicing of hnRNP A1 target RNAs involved in neuronal function and RNA homeostasis. We confirmed this in vivo in experimental autoimmune encephalomyelitis employing CLIPseq specific for hnRNP A1, where hnRNP A1 differentially binds and regulates RNA, including aberrantly spliced targets identified in human samples. Additionally, dysfunctional hnRNP A1 expression in neurons caused neurite loss and identical changes in splicing, corroborating hnRNP A1 dysfunction as a cause of neurodegeneration. Collectively, these data indicate hnRNP A1 dysfunction causes altered neuronal RNA splicing, resulting in neurodegeneration in MS.

Differential microRNA editing may drive target pathway switching in human temporal lobe epilepsy

MicroRNAs have emerged as important regulators of the gene expression landscape in temporal lobe epilepsy. The mechanisms that control microRNA levels and influence target choice remain, however, poorly understood. RNA editing is a post-transcriptional mechanism mediated by the adenosine acting on RNA (ADAR) family of proteins that introduces base modification that diversifies the gene expression landscape. RNA editing has been studied for the mRNA landscape but the extent to which microRNA editing occurs in human temporal lobe epilepsy is unknown. Here, we used small RNA-sequencing data to characterize the identity and extent of microRNA editing in human temporal lobe epilepsy brain samples. This detected low-to-high editing in over 40 of the identified microRNAs. Among microRNA exhibiting the highest editing was miR-376a-3p, which was edited in the seed region and this was predicted to significantly change the target pool. The edited form was expressed at lower levels in human temporal lobe epilepsy samples. We modelled the shift in editing levels of miR-376a-3p in human-induced pluripotent stem cell-derived neurons. Reducing levels of the edited form of miR-376a-3p using antisense oligonucleotides resulted in extensive gene expression changes, including upregulation of mitochondrial and metabolism-associated pathways. Together, these results show that differential editing of microRNAs may re-direct targeting and result in altered functions relevant to the pathophysiology of temporal lobe epilepsy and perhaps other disorders of neuronal hyperexcitability.

Comparative RNA Genomics

Over the last quarter of a century it has become clear that RNA is much more than just a boring intermediate in protein expression. Ancient RNAs still appear in the core information metabolism and comprise a surprisingly large component in bacterial gene regulation. A common theme with these types of mostly small RNAs is their reliance of conserved secondary structures. Large-scale sequencing projects, on the other hand, have profoundly changed our understanding of eukaryotic genomes. Pervasively transcribed, they give rise to a plethora of large and evolutionarily extremely flexible non-coding RNAs that exert a vastly diverse array of molecule functions. In this chapter we provide a-necessarily incomplete-overview of the current state of comparative analysis of non-coding RNAs, emphasizing computational approaches as a means to gain a global picture of the modern RNA world.

RBFOX2 deregulation promotes pancreatic cancer progression and metastasis through alternative splicing

RNA splicing is an important biological process associated with cancer initiation and progression. However, the contribution of alternative splicing to pancreatic cancer (PDAC) development is not well understood. Here, we identify an enrichment of RNA binding proteins (RBPs) involved in splicing regulation linked to PDAC progression from a forward genetic screen using Sleeping Beauty insertional mutagenesis in a mouse model of pancreatic cancer. We demonstrate downregulation of RBFOX2, an RBP of the FOX family, promotes pancreatic cancer progression and liver metastasis. Specifically, we show RBFOX2 regulates exon splicing events in transcripts encoding proteins involved in cytoskeletal remodeling programs. These exons are differentially spliced in PDAC patients, with enhanced exon skipping in the classical subtype for several RBFOX2 targets. RBFOX2 mediated splicing of ABI1, encoding the Abelson-interactor 1 adapter protein, controls the abundance and localization of ABI1 protein isoforms in pancreatic cancer cells and promotes the relocalization of ABI1 from the cytoplasm to the periphery of migrating cells. Using splice-switching antisense oligonucleotides (AONs) we demonstrate the ABI1 ∆Ex9 isoform enhances cell migration. Together, our data identify a role for RBFOX2 in promoting PDAC progression through alternative splicing regulation.

Master Transcription Factor Reprogramming Unleashes Selective Translation Promoting Castration Resistance and Immune Evasion in Lethal Prostate Cancer

Signaling rewiring allows tumors to survive therapy. Here we show that the decrease of the master regulator microphthalmia transcription factor (MITF) in lethal prostate cancer unleashes eukaryotic initiation factor 3B (eIF3B)-dependent translation reprogramming of key mRNAs conferring resistance to androgen deprivation therapy (ADT) and promoting immune evasion. Mechanistically, MITF represses through direct promoter binding eIF3B, which in turn regulates the translation of specific mRNAs. Genome-wide eIF3B enhanced cross-linking immunoprecipitation sequencing (eCLIP-seq) showed specialized binding to a UC-rich motif present in subsets of 5' untranslated regions. Indeed, translation of the androgen receptor and major histocompatibility complex I (MHC-I) through this motif is sensitive to eIF3B amount. Notably, pharmacologic targeting of eIF3B-dependent translation in preclinical models sensitizes prostate cancer to ADT and anti-PD-1 therapy. These findings uncover a hidden connection between transcriptional and translational rewiring promoting therapy-refractory lethal prostate cancer and provide a druggable mechanism that may transcend into effective combined therapeutic strategies.

SIGNIFICANCE

Our study shows that specialized eIF3B-dependent translation of specific mRNAs released upon downregulation of the master transcription factor MITF confers castration resistance and immune evasion in lethal prostate cancer. Pharmacologic targeting of this mechanism delays castration resistance and increases immune-checkpoint efficacy. This article is featured in Selected Articles from This Issue, p. 2489.

GCLiPP: global crosslinking and protein purification method for constructing high-resolution occupancy maps for RNA binding proteins

GCLiPP is a global RNA interactome capture method that detects RNA-binding protein (RBP) occupancy transcriptome-wide. GCLiPP maps RBP-occupied sites at a higher resolution than phase separation-based techniques. GCLiPP sequence tags correspond with known RBP binding sites and are enriched for sites detected by RBP-specific crosslinking immunoprecipitation (CLIP) for abundant cytosolic RBPs. Comparison of human Jurkat T cells and mouse primary T cells uncovers shared peaks of GCLiPP signal across homologous regions of human and mouse 3' UTRs, including a conserved mRNA-destabilizing cis-regulatory element. GCLiPP signal overlapping with immune-related SNPs uncovers stabilizing cis-regulatory regions in CD5, STAT6, and IKZF1.

RNA molecular recording with an engineered RNA deaminase

RNA deaminases are powerful tools for base editing and RNA molecular recording. However, the enzymes used in currently available RNA molecular recorders such as TRIBE, DART or STAMP have limitations due to RNA structure and sequence dependence. We designed a platform for directed evolution of RNA molecular recorders. We engineered an RNA A-to-I deaminase (an RNA adenosine base editor, rABE) that has high activity, low bias and low background. Using rABE, we present REMORA (RNA-encoded molecular recording in adenosines), wherein deamination by rABE writes a molecular record of RNA-protein interactions. By combining rABE with the C-to-U deaminase APOBEC1 and long-read RNA sequencing, we measured binding by two RNA-binding proteins on single messenger RNAs. Orthogonal RNA molecular recording of mammalian Pumilio proteins PUM1 and PUM2 shows that PUM1 competes with PUM2 for a subset of sites in cells. Furthermore, we identify transcript isoform-specific RNA-protein interactions driven by isoform changes distal to the binding site. The genetically encodable RNA deaminase rABE enables single-molecule identification of RNA-protein interactions with cell type specificity.

MicroRNA-mediated attenuation of branched-chain amino acid catabolism promotes ferroptosis in chronic kidney disease

Chronic kidney disease can develop from kidney injury incident to chemotherapy with cisplatin, which complicates the prognosis of cancer patients. MicroRNAs regulate gene expression by pairing with specific sets of messenger RNAs. Therefore, elucidating direct physical interactions between microRNAs and their target messenger RNAs can help decipher crucial biological processes associated with cisplatin-induced kidney injury. Through intermolecular ligation and transcriptome-wide sequencing, we here identify direct pairs of microRNAs and their target messenger RNAs in the kidney of male mice injured by cisplatin. We find that a group of cisplatin-induced microRNAs can target select messenger RNAs that affect the mitochondrial metabolic pathways in the injured kidney. Specifically, a cisplatin-induced microRNA, miR-429-3p, suppresses the pathway that catabolizes branched-chain amino acids in the proximal tubule, leading to cell death dependent on lipid peroxidation, called ferroptosis. Identification of miRNA-429-3p-mediated ferroptosis stimulation suggests therapeutic potential for modulating the branched-chain amino acid pathway in ameliorating cisplatin-induced kidney injury.

Modulation of insulin secretion by RBFOX2-mediated alternative splicing

Insulin secretion is a tightly regulated process that is vital for maintaining blood glucose homeostasis. Although the molecular components of insulin granule trafficking and secretion are well established, how they are regulated to rapidly fine-tune secretion in response to changing environmental conditions is not well characterized. Recent studies have determined that dysregulation of RNA-binding proteins (RBPs) and aberrant mRNA splicing occurs at the onset of diabetes. We demonstrate that the RBP, RBFOX2, is a critical regulator of insulin secretion through the alternative splicing of genes required for insulin granule docking and exocytosis. Conditional mutation of Rbfox2 in the mouse pancreas results in decreased insulin secretion and impaired blood glucose homeostasis. Consistent with defects in secretion, we observe reduced insulin granule docking and corresponding splicing defects in the SNARE complex components. These findings identify an additional mechanism for modulating insulin secretion in both healthy and dysfunctional pancreatic β cells.

Aberrant splicing in Huntington's disease via disrupted TDP-43 activity accompanied by altered m6A RNA modification

Huntington's disease (HD) is a neurodegenerative disorder caused by a CAG repeat expansion in the first exon of the HTT gene encoding huntingtin. Prior reports have established a correlation between CAG expanded HTT and altered gene expression. However, the mechanisms leading to disruption of RNA processing in HD remain unclear. Here, our analysis of the reported HTT protein interactome identifies interactions with known RNA-binding proteins (RBPs). Total, long-read sequencing and targeted RASL-seq of RNAs from cortex and striatum of the HD mouse model R6/2 reveals increased exon skipping which is confirmed in Q150 and Q175 knock-in mice and in HD human brain. We identify the RBP TDP-43 and the N6-methyladenosine (m6A) writer protein methyltransferase 3 (METTL3) to be upstream regulators of exon skipping in HD. Along with this novel mechanistic insight, we observe decreased nuclear localization of TDP-43 and cytoplasmic accumulation of phosphorylated TDP-43 in HD mice and human brain. In addition, TDP-43 co-localizes with HTT in human HD brain forming novel nuclear aggregate-like bodies distinct from mutant HTT inclusions or previously observed TDP-43 pathologies. Binding of TDP-43 onto RNAs encoding HD-associated differentially expressed and aberrantly spliced genes is decreased. Finally, m6A RNA modification is reduced on RNAs abnormally expressed in striatum from HD R6/2 mouse brain, including at clustered sites adjacent to TDP-43 binding sites. Our evidence supports TDP-43 loss of function coupled with altered m6A modification as a novel mechanism underlying alternative splicing/unannotated exon usage in HD and highlights the critical nature of TDP-43 function across multiple neurodegenerative diseases.

Role of stress granules in tumorigenesis and cancer therapy

Stress granules (SGs) are membrane-less organelles that cell forms via liquid-liquid phase separation (LLPS) under stress conditions such as oxidative stress, ER stress, heat shock and hypoxia. SG assembly is a stress-responsive mechanism by regulating gene expression and cellular signaling pathways. Cancer cells face various stress conditions in tumor microenvironment during tumorigenesis, while SGs contribute to hallmarks of cancer including proliferation, invasion, migration, avoiding apoptosis, metabolism reprogramming and immune evasion. Here, we review the connection between SGs and cancer development, the limitation of SGs on current cancer therapy and promising cancer therapeutic strategies targeting SGs in the future.

Splicing regulation of GFPT1 muscle-specific isoform and its roles in glucose metabolisms and neuromuscular junction

Glutamine:fructose-6-phosphate transaminase 1 (GFPT1) is the rate-limiting enzyme of the hexosamine biosynthetic pathway (HBP). A 54-bp exon 9 of GFPT1 is specifically included in skeletal and cardiac muscles to generate a long isoform of GFPT1 (GFPT1-L). We showed that SRSF1 and Rbfox1/2 cooperatively enhance, and hnRNP H/F suppresses, the inclusion of human GFPT1 exon 9 by modulating recruitment of U1 snRNP. Knockout (KO) of GFPT1-L in skeletal muscle markedly increased the amounts of GFPT1 and UDP-HexNAc, which subsequently suppressed the glycolytic pathway. Aged KO mice showed impaired insulin-mediated glucose uptake, as well as muscle weakness and fatigue likely due to abnormal formation and maintenance of the neuromuscular junction. Taken together, GFPT1-L is likely to be acquired in evolution in mammalian striated muscles to attenuate the HBP for efficient glycolytic energy production, insulin-mediated glucose uptake, and the formation and maintenance of the neuromuscular junction.

Hypoplastic Left Heart Syndrome: Signaling & Molecular Perspectives, and the Road Ahead

Hypoplastic left heart syndrome (HLHS) is a lethal congenital heart disease (CHD) affecting 8-25 per 100,000 neonates globally. Clinical interventions, primarily surgical, have improved the life expectancy of the affected subjects substantially over the years. However, the etiological basis of HLHS remains fundamentally unclear to this day. Based upon the existing paradigm of studies, HLHS exhibits a multifactorial mode of etiology mediated by a complicated course of genetic and signaling cascade. This review presents a detailed outline of the HLHS phenotype, the prenatal and postnatal risks, and the signaling and molecular mechanisms driving HLHS pathogenesis. The review discusses the potential limitations and future perspectives of studies that can be undertaken to address the existing scientific gap. Mechanistic studies to explain HLHS etiology will potentially elucidate novel druggable targets and empower the development of therapeutic regimens against HLHS in the future.

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.