Defining the regulatory network of the tissue-specific splicing factors Fox-1 and Fox-2

Rbfox1 controls alternative splicing of focal adhesion genes in cardiac muscle cells

Alternative splicing is one of the major cellular processes that determine the tissue-specific expression of protein variants. However, it remains challenging to identify physiologically relevant and tissue-selective proteins that are generated by alternative splicing. Hence, we investigated the target spectrum of the splicing factor Rbfox1 in the cardiac muscle context in more detail. By using a combination of in silico target prediction and in-cell validation, we identified several focal adhesion proteins as alternative splicing targets of Rbfox1. We focused on the alternative splicing patterns of vinculin (metavinculin isoform) and paxillin (extended paxillin isoform) and identified both as potential Rbfox1 targets. Minigene analyses suggested that both isoforms are promoted by Rbfox1 due to binding in the introns. Focal adhesions play an important role in the cardiac muscle context, since they mainly influence cell shape, cytoskeletal organization, and cell-matrix association. Our data confirmed that depletion of Rbfox1 changed cardiomyoblast morphology, cytoskeletal organization, and multinuclearity after differentiation, which might be due to changes in alternative splicing of focal adhesion proteins. Hence, our results indicate that Rbfox1 promotes alternative splicing of focal adhesion genes in cardiac muscle cells, which might contribute to heart disease progression, where downregulation of Rbfox1 is frequently observed.

Reverse engineering neuron type-specific and type-orthogonal splicing-regulatory networks using single-cell transcriptomes

Cell type-specific alternative splicing (AS) enables differential gene isoform expression between diverse neuron types with distinct identities and functions. Current studies linking individual RNA-binding proteins (RBPs) to AS in a few neuron types underscore the need for holistic modeling. Here, we use network reverse engineering to derive a map of the neuron type-specific AS regulatory landscape from 133 mouse neocortical cell types defined by single-cell transcriptomes. This approach reliably inferred the regulons of 350 RBPs and their cell type-specific activities. Our analysis revealed driving factors delineating neuronal identities, among which we validated Elavl2 as a key RBP for MGE-specific splicing in GABAergic interneurons using an in vitro ESC differentiation system. We also identified a module of exons and candidate regulators specific for long- and short-projection neurons across multiple neuronal classes. This study provides a resource for elucidating splicing regulatory programs that drive neuronal molecular diversity, including those that do not align with gene expression-based classifications.

Positioning of pyrimidine motifs around cassette exons defines their PTB-dependent splicing in Arabidopsis

Alternative splicing (AS) is a complex process that generates transcript variants from a single pre-mRNA and is involved in numerous biological functions. Many RNA-binding proteins are known to regulate AS; however, little is known about the underlying mechanisms, especially outside the mammalian clade. Here, we show that polypyrimidine tract binding proteins (PTBs) from Arabidopsis thaliana regulate AS of cassette exons via pyrimidine (Py)-rich motifs close to the alternative splice sites. Mutational studies on three PTB-dependent cassette exon events revealed that only some of the Py motifs in this region are critical for AS. Moreover, in vitro binding of PTBs did not reflect a motif's impact on AS in vivo. Our mutational studies and bioinformatic investigation of all known PTB-regulated cassette exons from A. thaliana and human suggested that the binding position of PTBs relative to a cassette exon defines whether its inclusion or skipping is induced. Accordingly, exon skipping is associated with a higher frequency of Py stretches within the cassette exon, and in human also upstream of it, whereas exon inclusion is characterized by increased Py motif occurrence downstream of said exon. Enrichment of Py motifs downstream of PTB-activated 5' splice sites is also seen for PTB-dependent intron removal and alternative 5' splice site events from A. thaliana, suggesting this is a common step of exon definition. In conclusion, the position-dependent AS regulatory mechanism by PTB homologs has been conserved during the separate evolution of plants and mammals, while other critical features, in particular intron length, have considerably changed.

Alternative splicing in EMT and TGF-β signaling during cancer progression

Epithelial to mesenchymal transition (EMT) is a physiological process during development where epithelial cells transform to acquire mesenchymal characteristics, which allows them to migrate and colonize secondary tissues. Many cellular signaling pathways and master transcriptional factors exert a myriad of controls to fine tune this vital process to meet various developmental and physiological needs. Adding to the complexity of this network are post-transcriptional and post-translational regulations. Among them, alternative splicing has been shown to play important roles to drive EMT-associated phenotypic changes, including actin cytoskeleton remodeling, cell-cell junction changes, cell motility and invasiveness. In advanced cancers, transforming growth factor-β (TGF-β) is a major inducer of EMT and is associated with tumor cell metastasis, cancer stem cell self-renewal, and drug resistance. This review aims to provide an overview of recent discoveries regarding alternative splicing events and the involvement of splicing factors in the EMT and TGF-β signaling. It will emphasize the importance of various splicing factors involved in EMT and explore their regulatory mechanisms.

Aberrant splicing of CaV1.2 calcium channel induced by decreased Rbfox1 enhances arterial constriction during diabetic hyperglycemia

Diabetic hyperglycemia induces dysfunctions of arterial smooth muscle, leading to diabetic vascular complications. The CaV1.2 calcium channel is one primary pathway for Ca2+ influx, which initiates vasoconstriction. However, the long-term regulation mechanism(s) for vascular CaV1.2 functions under hyperglycemic condition remains unknown. Here, Sprague-Dawley rats fed with high-fat diet in combination with low dose streptozotocin and Goto-Kakizaki (GK) rats were used as diabetic models. Isolated mesenteric arteries (MAs) and vascular smooth muscle cells (VSMCs) from rat models were used to assess K+-induced arterial constriction and CaV1.2 channel functions using vascular myograph and whole-cell patch clamp, respectively. K+-induced vasoconstriction is persistently enhanced in the MAs from diabetic rats, and CaV1.2 alternative spliced exon 9* is increased, while exon 33 is decreased in rat diabetic arteries. Furthermore, CaV1.2 channels exhibit hyperpolarized current-voltage and activation curve in VSMCs from diabetic rats, which facilitates the channel function. Unexpectedly, the application of glycated serum (GS), mimicking advanced glycation end-products (AGEs), but not glucose, downregulates the expression of the splicing factor Rbfox1 in VSMCs. Moreover, GS application or Rbfox1 knockdown dynamically regulates alternative exons 9* and 33, leading to facilitated functions of CaV1.2 channels in VSMCs and MAs. Notably, GS increases K+-induced intracellular calcium concentration of VSMCs and the vasoconstriction of MAs. These results reveal that AGEs, not glucose, long-termly regulates CaV1.2 alternative splicing events by decreasing Rbfox1 expression, thereby enhancing channel functions and increasing vasoconstriction under diabetic hyperglycemia. This study identifies the specific molecular mechanism for enhanced vasoconstriction under hyperglycemia, providing a potential target for managing diabetic vascular complications.

Quantifying transcriptome diversity: a review

Following the central dogma of molecular biology, gene expression heterogeneity can aid in predicting and explaining the wide variety of protein products, functions and, ultimately, heterogeneity in phenotypes. There is currently overlapping terminology used to describe the types of diversity in gene expression profiles, and overlooking these nuances can misrepresent important biological information. Here, we describe transcriptome diversity as a measure of the heterogeneity in (1) the expression of all genes within a sample or a single gene across samples in a population (gene-level diversity) or (2) the isoform-specific expression of a given gene (isoform-level diversity). We first overview modulators and quantification of transcriptome diversity at the gene level. Then, we discuss the role alternative splicing plays in driving transcript isoform-level diversity and how it can be quantified. Additionally, we overview computational resources for calculating gene-level and isoform-level diversity for high-throughput sequencing data. Finally, we discuss future applications of transcriptome diversity. This review provides a comprehensive overview of how gene expression diversity arises, and how measuring it determines a more complete picture of heterogeneity across proteins, cells, tissues, organisms and species.

Mutation of two intronic nucleotides alters RNA structure and dynamics inhibiting MBNL1 and RBFOX1 regulated splicing of the Insulin Receptor

Alternative splicing (AS) of Exon 11 of the Insulin Receptor ( INSR ) is highly regulated and disrupted in several human disorders. To better understand INSR exon 11 AS regulation, splicing activity of an INSR exon 11 minigene reporter was measured across a gradient of the AS regulator muscleblind-like 1 protein (MBNL1). The RNA-binding protein Fox-1 (RBFOX1) was added to determine its impact on MBNL1-regulated splicing. The role of the RBFOX1 UGCAUG binding site within intron 11 was assessed across the MBNL1 gradient. Mutating the UGCAUG motif inhibited RBFOX1 regulation of exon 11 and had the unexpected effect of reducing MBNL1 regulation of this exon. Molecular dynamics simulations showed that exon 11 and the adjacent RNA adopts a dynamically stable conformation. Mutation of the RBFOX1 binding site altered RNA structure and dynamics, while a mutation that created an optimal MBNL1 binding site at the RBFOX1 site shifted the RNA back to wild type. An antisense oligonucleotide (ASO) was used to confirm the structure in this region of the pre-mRNA. This example of intronic mutations shifting pre-mRNA structure and dynamics to modulate splicing suggests RNA structure and dynamics should be taken into consideration for AS regulation and therapeutic interventions targeting pre-mRNA.

Dysregulation of alternative splicing in spinocerebellar ataxia type 1

Spinocerebellar ataxia type 1 is caused by an expansion of the polyglutamine tract in ATAXIN-1. Ataxin-1 is broadly expressed throughout the brain and is involved in regulating gene expression. However, it is not yet known if mutant ataxin-1 can impact the regulation of alternative splicing events. We performed RNA sequencing in mouse models of spinocerebellar ataxia type 1 and identified that mutant ataxin-1 expression abnormally leads to diverse splicing events in the mouse cerebellum of spinocerebellar ataxia type 1. We found that the diverse splicing events occurred in a predominantly cell autonomous manner. A majority of the transcripts with misregulated alternative splicing events were previously unknown, thus allowing us to identify overall new biological pathways that are distinctive to those affected by differential gene expression in spinocerebellar ataxia type 1. We also provide evidence that the splicing factor Rbfox1 mediates the effect of mutant ataxin-1 on misregulated alternative splicing and that genetic manipulation of Rbfox1 expression modifies neurodegenerative phenotypes in a Drosophila model of spinocerebellar ataxia type 1 in vivo. Together, this study provides novel molecular mechanistic insight into the pathogenesis of spinocerebellar ataxia type 1 and identifies potential therapeutic strategies for spinocerebellar ataxia type 1.

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.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Noncoding snoRNA host genes are a distinct subclass of long noncoding RNAs

Mammalian genomes are pervasively transcribed into different noncoding (nc)RNA classes, each one with its own hallmarks and exceptions. Some of them are nested into each other, such as host genes for small nucleolar RNAs (snoRNAs), which were long believed to simply act as molecular containers strictly facilitating snoRNA biogenesis. However, recent findings show that noncoding snoRNA host genes (ncSNHGs) display features different from those of 'regular' long ncRNAs (lncRNAs) and, more importantly, they can exert independent and unrelated functions to those of the encoded snoRNAs. Here, we review and summarize past and recent evidence that ncSNHGs form a defined subclass among the plethora of lncRNAs, and discuss future research that can further elucidate their biological relevance.

To RNA-binding and beyond: Emerging facets of the role of Rbfox proteins in development and disease

The RNA-binding Fox-1 homologue (Rbfox) proteins represent an ancient family of splicing factors, conserved through evolution. All members share an RNA recognition motif (RRM), and a particular affinity for the GCAUG signature in target RNA molecules. The role of Rbfox, as a splice factor, deciding the tissue-specific inclusion/exclusion of an exon, depending on its binding position on the flanking introns, is well known. Rbfox often acts in concert with other splicing factors, and forms splicing regulatory networks. Apart from this canonical role, recent studies show that Rbfox can also function as a transcription co-factor, and affects mRNA stability and translation. The repertoire of Rbfox targets is vast, including genes involved in the development of tissue lineages, such as neurogenesis, myogenesis, and erythropoeiesis, and molecular processes, including cytoskeletal dynamics, and calcium handling. A second layer of complexity is added by the fact that Rbfox expression itself is regulated by multiple mechanisms, and, in vertebrates, exhibits tissue-specific expression. The optimum dosage of Rbfox is critical, and its misexpression is etiological to various disease conditions. In this review, we discuss the contextual roles played by Rbfox as a tissue-specific regulator for the expression of many important genes with diverse functions, through the lens of the emerging data which highlights its involvement in many human diseases. Furthermore, we explore the mechanistic details provided by studies in model organisms, with emphasis on the work with Drosophila. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Turnover and Surveillance > Regulation of RNA Stability RNA Processing > Splicing Regulation/Alternative Splicing.

Alternative splicing outcomes across an RNA-binding protein concentration gradient

Alternative splicing (AS) is a dynamic RNA processing step that produces multiple RNA isoforms from a single pre-mRNA transcript and contributes to the complexity of the cellular transcriptome and proteome. This process is regulated through a network of cis-regulatory sequence elements and trans-acting factors, most-notably RNA binding proteins (RBPs). The muscleblind-like (MBNL) and RNA binding fox-1 homolog (RBFOX) are two well characterized families of RBPs that regulate fetal to adult AS transitions critical for proper muscle, heart, and central nervous system development. To better understand how the concentration of these RBPs influences AS transcriptome wide, we engineered a MBNL1 and RBFOX1 inducible HEK-293 cell line. Modest induction of exogenous RBFOX1 in this cell line modulated MBNL1-dependent AS outcomes in 3 skipped exon events, despite significant levels of endogenous RBFOX1 and RBFOX2. Due to background RBFOX levels, we conducted a focused analysis of dose-dependent MBNL1 skipped exon AS outcomes and generated transcriptome wide dose-response curves. Analysis of this data demonstrates that MBNL1-regulated exclusion events may require higher concentrations of MBNL1 protein to properly regulate AS outcomes compared to inclusion events and that multiple arrangements of YGCY motifs can produce similar splicing outcomes. These results suggest that rather than a simple relationship between the organization of RBP binding sites and a specific splicing outcome, that complex interaction networks govern both AS inclusion and exclusion events across a RBP gradient.

mRNA isoform balance in neuronal development and disease

Balanced mRNA isoform diversity and abundance are spatially and temporally regulated throughout cellular differentiation. The proportion of expressed isoforms contributes to cell type specification and determines key properties of the differentiated cells. Neurons are unique cell types with intricate developmental programs, characteristic cellular morphologies, and electrophysiological potential. Neuron-specific gene expression programs establish these distinctive cellular characteristics and drive diversity among neuronal subtypes. Genes with neuron-specific alternative processing are enriched in key neuronal functions, including synaptic proteins, adhesion molecules, and scaffold proteins. Despite the similarity of neuronal gene expression programs, each neuronal subclass can be distinguished by unique alternative mRNA processing events. Alternative processing of developmentally important transcripts alters coding and regulatory information, including interaction domains, transcript stability, subcellular localization, and targeting by RNA binding proteins. Fine-tuning of mRNA processing is essential for neuronal activity and maintenance. Thus, the focus of neuronal RNA biology research is to dissect the transcriptomic mechanisms that underlie neuronal homeostasis, and consequently, predispose neuronal subtypes to disease. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA in Disease and Development > RNA in Development.

The aptamer-based RNA-PROTAC

RNA-binding proteins (RBPs) dysfunction has been implicated in a number of diseases, and RBPs have traditionally been considered to be undruggable targets. Here, targeted degradation of RBPs is achieved based on the aptamer-based RNA-PROTAC, which consists of a genetically encoded RNA scaffold and a synthetic heterobifunctional molecule. The target RBPs can bind to their RNA consensus binding element (RCBE) on the RNA scaffold, while the small molecule can recruit E3 ubiquitin ligase to the RNA scaffold in a non-covalent manner, thereby inducing proximity-dependent ubiquitination and subsequent proteasome-mediated degradation of the target protein. Different RBPs targets, including LIN28A and RBFOX1, have been successfully degraded by simply replacing the RCBE module on the RNA scaffold. In addition, the simultaneous degradation of multiple target proteins has been realized by inserting more functional RNA oligonucleotides into the RNA scaffold.

A cell-type-specific alternative splicing regulator shapes synapse properties in a trans-synaptic manner

The specification of synaptic properties is fundamental for the function of neuronal circuits. "Terminal selector" transcription factors coordinate terminal gene batteries that specify cell-type-specific properties. Moreover, pan-neuronal splicing regulators have been implicated in directing neuronal differentiation. However, the cellular logic of how splicing regulators instruct specific synaptic properties remains poorly understood. Here, we combine genome-wide mapping of mRNA targets and cell-type-specific loss-of-function studies to uncover the contribution of the RNA-binding protein SLM2 to hippocampal synapse specification. Focusing on pyramidal cells and somatostatin (SST)-positive GABAergic interneurons, we find that SLM2 preferentially binds and regulates alternative splicing of transcripts encoding synaptic proteins. In the absence of SLM2, neuronal populations exhibit normal intrinsic properties, but there are non-cell-autonomous synaptic phenotypes and associated defects in a hippocampus-dependent memory task. Thus, alternative splicing provides a critical layer of gene regulation that instructs specification of neuronal connectivity in a trans-synaptic manner.

Liver RBFOX2 regulates cholesterol homeostasis via Scarb1 alternative splicing in mice

RNA alternative splicing (AS) expands the regulatory potential of eukaryotic genomes. The mechanisms regulating liver-specific AS profiles and their contribution to liver function are poorly understood. Here, we identify a key role for the splicing factor RNA-binding Fox protein 2 (RBFOX2) in maintaining cholesterol homeostasis in a lipogenic environment in the liver. Using enhanced individual-nucleotide-resolution ultra-violet cross-linking and immunoprecipitation, we identify physiologically relevant targets of RBFOX2 in mouse liver, including the scavenger receptor class B type I (Scarb1). RBFOX2 function is decreased in the liver in diet-induced obesity, causing a Scarb1 isoform switch and alteration of hepatocyte lipid homeostasis. Our findings demonstrate that specific AS programmes actively maintain liver physiology, and underlie the lipotoxic effects of obesogenic diets when dysregulated. Splice-switching oligonucleotides targeting this network alleviate obesity-induced inflammation in the liver and promote an anti-atherogenic lipoprotein profile in the blood, underscoring the potential of isoform-specific RNA therapeutics for treating metabolism-associated diseases.

A widespread length-dependent splicing dysregulation in cancer

Dysregulation of alternative splicing is a key molecular hallmark of cancer. However, the common features and underlying mechanisms remain unclear. Here, we report an intriguing length-dependent splicing regulation in cancers. By systematically analyzing the transcriptome of thousands of cancer patients, we found that short exons are more likely to be mis-spliced and preferentially excluded in cancers. Compared to other exons, cancer-associated short exons (CASEs) are more conserved and likely to encode in-frame low-complexity peptides, with functional enrichment in GTPase regulators and cell adhesion. We developed a CASE-based panel as reliable cancer stratification markers and strong predictors for survival, which is clinically useful because the detection of short exon splicing is practical. Mechanistically, mis-splicing of CASEs is regulated by elevated transcription and alteration of certain RNA binding proteins in cancers. Our findings uncover a common feature of cancer-specific splicing dysregulation with important clinical implications in cancer diagnosis and therapies.

Specialization of the photoreceptor transcriptome by Srrm3-dependent microexons is required for outer segment maintenance and vision

Retinal photoreceptors have a distinct transcriptomic profile compared to other neuronal subtypes, likely reflecting their unique cellular morphology and function in the detection of light stimuli by way of the ciliary outer segment. We discovered a layer of this molecular specialization by revealing that the vertebrate retina expresses the largest number of tissue-enriched microexons of all tissue types. A subset of these microexons is included exclusively in photoreceptor transcripts, particularly in genes involved in cilia biogenesis and vesicle-mediated transport. This microexon program is regulated by Srrm3, a paralog of the neural microexon regulator Srrm4. Despite the fact that both proteins positively regulate retina microexons in vitro, only Srrm3 is highly expressed in mature photoreceptors. Its deletion in zebrafish results in widespread down-regulation of microexon inclusion from early developmental stages, followed by other transcriptomic alterations, severe photoreceptor defects, and blindness. These results shed light on the transcriptomic specialization and functionality of photoreceptors, uncovering unique cell type-specific roles for Srrm3 and microexons with implications for retinal diseases.

RBFOX2-regulated TEAD1 alternative splicing plays a pivotal role in Hippo-YAP signaling

Alternative pre-mRNA splicing is key to proteome diversity; however, the biological roles of alternative splicing (AS) in signaling pathways remain elusive. Here, we focus on TEA domain transcription factor 1 (TEAD1), a YAP binding factor in the Hippo signaling pathway. Public database analyses showed that expression of YAP-TEAD target genes negatively correlated with the expression of a TEAD1 isoform lacking exon 6 (TEAD1ΔE6) but did not correlate with overall TEAD1 expression. We confirmed that the transcriptional activity and oncogenic properties of the full-length TEAD1 isoform were greater than those of TEAD1ΔE6, with the difference in transcription related to YAP interaction. Furthermore, we showed that RNA-binding Fox-1 homolog 2 (RBFOX2) promoted the inclusion of TEAD1 exon 6 via binding to the conserved GCAUG element in the downstream intron. These results suggest a regulatory mechanism of RBFOX2-mediated TEAD1 AS and provide insight into AS-specific modulation of signaling pathways.

The comprehensive and systematic identification of BLCA-specific SF-regulated, survival-related AS events

Bladder urothelial carcinoma (BLCA) is a complex disease with high morbidity and mortality. Changes in alternative splicing (AS) and splicing factor (SF) can affect gene expression, thus playing an essential role in tumorigenesis. This study downloaded 412 patients' clinical information and 433 samples of transcriptome profiling data from TCGA. And we collected 48 AS signatures from SpliceSeq. LASSO and Cox analyses were used for identifying survival-related AS events in BLCA. Finally, 1,645 OS-related AS events in 1,129 genes were validated by Kaplan-Meier (KM) survival analysis, ROC analysis, risk curve analysis, and independent prognostic analysis. Finally, our survey provides an AS-SF regulation network consisting of five SFs and 46 AS events. In the end, we profiled genes that AS occurred in pan-cancer and five SFs' expression in tumor and normal samples in BLCA. We selected CLIP-seq data for validation the interaction regulated by RBP. Our study paves the way for potential therapeutic targets of BLCA.

ASpediaFI: Functional Interaction Analysis of Alternative Splicing Events

Alternative splicing (AS) regulates biological processes governing phenotypes and diseases. Differential AS (DAS) gene test methods have been developed to investigate important exonic expression from high-throughput datasets. However, the DAS events extracted using statistical tests are insufficient to delineate relevant biological processes. In this study, we developed a novel application, Alternative Splicing Encyclopedia: Functional Interaction (ASpediaFI), to systemically identify DAS events and co-regulated genes and pathways. ASpediaFI establishes a heterogeneous interaction network of genes and their feature nodes (i.e., AS events and pathways) connected by co-expression or pathway gene set knowledge. Next, ASpediaFI explores the interaction network using the random walk with restart algorithm and interrogates the proximity from a query gene set. Finally, ASpediaFI extracts significant AS events, genes, and pathways. To evaluate the performance of our method, we simulated RNA sequencing (RNA-seq) datasets to consider various conditions of sequencing depth and sample size. The performance was compared with that of other methods. Additionally, we analyzed three public datasets of cancer patients or cell lines to evaluate how well ASpediaFI detects biologically relevant candidates. ASpediaFI exhibits strong performance in both simulated and public datasets. Our integrative approach reveals that DAS events that recognize a global co-expression network and relevant pathways determine the functional importance of spliced genes in the subnetwork. ASpediaFI is publicly available at https://bioconductor.org/packages/ASpediaFI.

Alternative Splicing of the Aryl Hydrocarbon Receptor Nuclear Translocator (ARNT) Is Regulated by RBFOX2 in Lymphoid Malignancies

Aberrant alternative splicing (AS) of pre-mRNAs promotes the development and proliferation of cancerous cells. Accordingly, we had previously observed higher levels of the aryl hydrocarbon receptor nuclear translocator (ARNT) spliced variant isoform 1 in human lymphoid malignancies compared to that in normal lymphoid cells, which is a consequence of increased inclusion of alternative exon 5. ARNT is a transcription factor that has been implicated in the survival of various cancers. Notably, we found that ARNT isoform 1 promoted the growth and survival of lymphoid malignancies, but the regulatory mechanism controlling ARNT AS is unclear. Here, we report cis- and trans-regulatory elements which are important for the inclusion of ARNT exon 5. Specifically, we identified recognition motifs for the RNA-binding protein RBFOX2, which are required for RBFOX2-mediated exon 5 inclusion. RBFOX2 upregulation was observed in lymphoid malignancies, correlating with the observed increase in ARNT exon 5 inclusion. Moreover, suppression of RBFOX2 significantly reduced ARNT isoform 1 levels and cell growth. These observations reveal RBFOX2 as a critical regulator of ARNT AS in lymphoid malignancies and suggest that blocking the ARNT-specific RBFOX2 motifs to decrease ARNT isoform 1 levels is a viable option for targeting the growth of lymphoid malignancies.

Autophagy regulation by RNA alternative splicing and implications in human diseases

Autophagy and RNA alternative splicing are two evolutionarily conserved processes involved in overlapping physiological and pathological processes. However, the extent of functional connection is not well defined. Here, we consider the role for alternative splicing and generation of autophagy-related gene isoforms in the regulation of autophagy in recent work. The impact of changes to the RNA alternative splicing machinery and production of alternative spliced isoforms on autophagy are reviewed with particular focus on disease relevance. The use of drugs targeting both alternative splicing and autophagy as well as the selective regulation of single autophagy-related protein isoforms, are considered as therapeutic strategies.

Both alternative splicing and autophagy are core cell biological processes, but where they intersect has received little attention. Here, the authors reflect on recent connections identified between these pathways and consider their impact on human disease.

Evolution Increases Primates Brain Complexity Extending RbFOX1 Splicing Activity to LSD1 Modulation

Recent branching (100 MYA) of the mammalian evolutionary tree has enhanced brain complexity and functions at the putative cost of increased emotional circuitry vulnerability. Thus, to better understand psychopathology, a burden for the modern society, novel approaches should exploit evolutionary aspects of psychiatric-relevant molecular pathways. A handful of genes is nowadays tightly associated to psychiatric disorders. Among them, neuronal-enriched RbFOX1 modifies the activity of synaptic regulators in response to neuronal activity, keeping excitability within healthy domains. We here dissect a higher primates-restricted interaction between RbFOX1 and the transcriptional corepressor Lysine Specific Demethylase 1 (LSD1/KDM1A). A single nucleotide variation (AA to AG) in LSD1 gene appeared in higher primates and humans, endowing RbFOX1 with the ability to promote the alternative usage of a novel 3' AG splice site, which extends LSD1 exon E9 in the upstream intron (E9-long). Exon E9-long regulates LSD1 levels by Nonsense-Mediated mRNA Decay. As reintroduction of the archaic LSD1 variant (AA) abolishes E9-long splicing, the novel 3' AG splice site is necessary for RbFOX1 to control LSD1 levels. LSD1 is a homeostatic immediate early genes (IEGs) regulator playing a relevant part in environmental stress-response. In primates and humans, inclusion of LSD1 as RbFOX1 target provides RbFOX1 with the additional ability to regulate the IEGs. These data, together with extensive RbFOX1 involvement in psychiatric disorders and its stress-dependent regulation in male mice, suggest the RbFOX1-LSD1-IEGs axis as an evolutionary recent psychiatric-relevant pathway. Notably, outside the nervous system, RbFOX2-dependent LSD1 modulation could be a candidate deregulated mechanism in cancer.SIGNIFICANCE STATEMENT To be better understood, anxiety and depression need large human genetics studies aimed at further resolving the often ambiguous, aberrant neuronal pathomechanisms that impact corticolimbic circuitry physiology. Several genetic associations of the alternative splicing regulator RbFOX1 with psychiatric conditions suggest homeostatic unbalance as a neuronal signature of psychopathology. Here we move a step forward, characterizing a disease-relevant higher primates-specific pathway by which RbFOX1 acquires the ability to regulate neuronal levels of Lysine Specific Demethylase 1, an epigenetic modulator of environmental stress response. Thus, two brain-enriched enzymes, independently shown to homeostatically protect neurons with a clear readout in terms of emotional behavior in lower mammals, establish in higher primates and humans a new functional cooperation enhancing the complexity of environmental adaptation and stress vulnerability.

Characterization of distinct circular RNA signatures in solid tumors

Background

Circular RNAs (circRNAs) are differentially expressed between normal and cancerous tissues, contributing to tumor initiation and progression. However, comprehensive landscape of dysregulated circRNAs across cancer types remains unclear.

Methods

In this study, we conducted Ribo-Zero transcriptome sequencing on tumor tissues and their adjacent normal samples including glioblastoma, esophageal squamous cell carcinoma, lung adenocarcinoma, thyroid cancer, colorectal cancer, gastric cancer and hepatocellular carcinoma. CIRCexplorer2 was employed to identify circRNAs and dysregulated circRNAs and genes were determined by DESeq2 package. The expression of hsa_circ_0072309 (circLIFR) was measured by reverse transcription and quantitative real-time PCR, and its effect on cell migration was examined by Transwell and wound healing assays. The role of circLIFR in tumor metastasis was evaluated via mouse models of tail-vein injection and spleen injection for lung and liver metastasis, respectively.

Results

Distinct circRNA expression signatures were identified among seven types of solid tumors, and the dysregulated circRNAs exhibited cancer-specific expression or shared common expression signatures across cancers. Bioinformatics analyses indicated that aberrant expression of host genes and/or RNA-binding proteins contributed to circRNA dysregulation in cancer. Finally, circLIFR was experimentally validated to be downregulated in six solid tumors and to significantly inhibit cell migration in vitro and tumor metastasis in vivo.

Conclusions

Our results provide a comprehensive landscape of differentially expressed circRNAs in solid tumors and highlight that circRNAs are extensively involved in cancer pathogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-022-01546-4.

The Splicing of the Mitochondrial Calcium Uniporter Genuine Activator MICU1 Is Driven by RBFOX2 Splicing Factor during Myogenic Differentiation

Alternative splicing, the process by which exons within a pre-mRNA transcript are differentially joined or skipped, is crucial in skeletal muscle since it is required both during myogenesis and in post-natal life to reprogram the transcripts of contractile proteins, metabolic enzymes, and transcription factors in functionally distinct muscle fiber types. The importance of such events is underlined by the numerosity of pathological conditions caused by alternative splicing aberrations. Importantly, many skeletal muscle Ca²⁺ homeostasis genes are also regulated by alternative splicing mechanisms, among which is the Mitochondrial Ca²⁺ Uniporter (MCU) genuine activator MICU1 which regulates MCU opening upon cell stimulation. We have previously shown that murine skeletal muscle MICU1 is subjected to alternative splicing, thereby generating a splice variant-which was named MICU1.1-that confers unique properties to the mitochondrial Ca²⁺ uptake and ensuring sufficient ATP production for muscle contraction. Here we extended the analysis of MICU1 alternative splicing to human tissues, finding two additional splicing variants that were characterized by their ability to regulate mitochondrial Ca²⁺ uptake. Furthermore, we found that MICU1 alternative splicing is induced during myogenesis by the splicing factor RBFOX2. These results highlight the complexity of the alternative splicing mechanisms in skeletal muscle and the regulation of mitochondrial Ca²⁺ among tissues.

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.

Regulation of Alternative Splicing during Epithelial-Mesenchymal Transition

Alternative splicing is an essential mechanism of gene regulation, giving rise to remarkable protein diversity in higher eukaryotes. Epithelial-mesenchymal transition (EMT) is a developmental process that plays an essential role in metazoan embryogenesis. Recent studies have revealed that alternative splicing serves as a fundamental layer of regulation that governs cells to undergo EMT. In this review, we summarize recent findings on the functional impact of alternative splicing in EMT and EMT-associated activities. We then discuss the regulatory mechanisms that control alternative splicing changes during EMT.

Rett syndrome linked to defects in forming the MeCP2/Rbfox/LASR complex in mouse models

Rett syndrome (RTT) is a severe neurological disorder and a leading cause of intellectual disability in young females. RTT is mainly caused by mutations found in the X-linked gene encoding methyl-CpG binding protein 2 (MeCP2). Despite extensive studies, the molecular mechanism underlying RTT pathogenesis is still poorly understood. Here, we report MeCP2 as a key subunit of a higher-order multiunit protein complex Rbfox/LASR. Defective MeCP2 in RTT mouse models disrupts the assembly of the MeCP2/Rbfox/LASR complex, leading to reduced binding of Rbfox proteins to target pre-mRNAs and aberrant splicing of Nrxns and Nlgn1 critical for synaptic plasticity. We further show that MeCP2 disease mutants display defective condensate properties and fail to promote phase-separated condensates with Rbfox proteins in vitro and in cultured cells. These data link an impaired function of MeCP2 with disease mutation in splicing control to its defective properties in mediating the higher-order assembly of the MeCP2/Rbfox/LASR complex.

MeCP2 mutations can cause Rett syndrome, a severe childhood neurological disorder. Here the authors show that MeCP2 mediates the higher-order assembly of a large splicing complex Rbfox/LASR, which is disrupted in the mouse models of Rett syndrome.

Transcriptome programs involved in the development and structure of the cerebellum

In the past two decades, mounting evidence has modified the classical view of the cerebellum as a brain region specifically involved in the modulation of motor functions. Indeed, clinical studies and engineered mouse models have highlighted cerebellar circuits implicated in cognitive functions and behavior. Furthermore, it is now clear that insults occurring in specific time windows of cerebellar development can affect cognitive performance later in life and are associated with neurological syndromes, such as Autism Spectrum Disorder. Despite its almost homogenous cytoarchitecture, how cerebellar circuits form and function is not completely elucidated yet. Notably, the apparently simple neuronal organization of the cerebellum, in which Purkinje cells represent the only output, hides an elevated functional diversity even within the same neuronal population. Such complexity is the result of the integration of intrinsic morphogenetic programs and extracellular cues from the surrounding environment, which impact on the regulation of the transcriptome of cerebellar neurons. In this review, we briefly summarize key features of the development and structure of the cerebellum before focusing on the pathways involved in the acquisition of the cerebellar neuron identity. We focus on gene expression and mRNA processing programs, including mRNA methylation, trafficking and splicing, that are set in motion during cerebellar development and participate to its physiology. These programs are likely to add new layers of complexity and versatility that are fundamental for the adaptability of cerebellar neurons.

The Genetic Control of Stoichiometry Underlying Autism

Autism is a common and complex neurologic disorder whose scientific underpinnings have begun to be established in the past decade. The essence of this breakthrough has been a focus on families, where genetic analyses are strongest, versus large-scale, case-control studies. Autism genetics has progressed in parallel with technology, from analyses of copy number variation to whole-exome sequencing (WES) and whole-genome sequencing (WGS). Gene mutations causing complete loss of function account for perhaps one-third of cases, largely detected through WES. This limitation has increased interest in understanding the regulatory variants of genes that contribute in more subtle ways to the disorder. Strategies combining biochemical analysis of gene regulation, WGS analysis of the noncoding genome, and machine learning have begun to succeed. The emerging picture is that careful control of the amounts of transcription, mRNA, and proteins made by key brain genes-stoichiometry-plays a critical role in defining the clinical features of autism.

RBFOX2/GOLIM4 Splicing Axis Activates Vesicular Transport Pathway to Promote Nasopharyngeal Carcinogenesis

20-30% of patients with nasopharyngeal carcinoma (NPC) develop distant metastasis or recurrence leading to poor survival, of which the underlying key molecular events have yet to be addressed. Here alternative splicing events in 85 NPC samples are profiled using transcriptome analysis and it is revealed that the long isoform of GOLIM4 (-L) with exon-7 is highly expressed in NPC and associated with poor prognosis. Lines of evidence demonstrate the pro-tumorigenic function of GOLIM4-L in NPC cells. It is further revealed that RBFOX2 binds to a GGAA motif in exon-7 and promotes its inclusion forming GOLIM4-L. RBFOX2 knockdown suppresses the tumorigenesis of NPC cells, phenocopying GOLIM4-L knockdown, which is significantly rescued by GOLIM4-L overexpression. High expression of RBFOX2 is correlated with the exon-7 inclusion of GOLIM4 in NPC biopsies and associated with worse prognosis. It is observed that RBFOX2 and GOLIM4 can influence vesicle-mediated transport through maintaining the organization of Golgi apparatus. Finally, it is revealed that RAB26 interacts with GOLIM4 and mediates its tumorigenic potentials in NPC cells. Taken together, the findings provide insights into how alternative splicing contributes to NPC development, by highlighting a functional link between GOLIM4-L and its splicing regulator RBFOX2 activating vesicle-mediated transport involving RAB26.

ERG transcription factors have a splicing regulatory function involving RBFOX2 that is altered in the EWS-FLI1 oncogenic fusion

ERG family proteins (ERG, FLI1 and FEV) are a subfamily of ETS transcription factors with key roles in physiology and development. In Ewing sarcoma, the oncogenic fusion protein EWS-FLI1 regulates both transcription and alternative splicing of pre-messenger RNAs. However, whether wild-type ERG family proteins might regulate splicing is unknown. Here, we show that wild-type ERG proteins associate with spliceosomal components, are found on nascent RNAs, and induce alternative splicing when recruited onto a reporter minigene. Transcriptomic analysis revealed that ERG and FLI1 regulate large numbers of alternative spliced exons (ASEs) enriched with RBFOX2 motifs and co-regulated by this splicing factor. ERG and FLI1 are associated with RBFOX2 via their conserved carboxy-terminal domain, which is present in EWS-FLI1. Accordingly, EWS-FLI1 is also associated with RBFOX2 and regulates ASEs enriched in RBFOX2 motifs. However, in contrast to wild-type ERG and FLI1, EWS-FLI1 often antagonizes RBFOX2 effects on exon inclusion. In particular, EWS-FLI1 reduces RBFOX2 binding to the ADD3 pre-mRNA, thus increasing its long isoform, which represses the mesenchymal phenotype of Ewing sarcoma cells. Our findings reveal a RBFOX2-mediated splicing regulatory function of wild-type ERG family proteins, that is altered in EWS-FLI1 and contributes to the Ewing sarcoma cell phenotype.

Rbfox2 mediates exon 11 inclusion in insulin receptor pre-mRNA splicing in hepatoma cells

Insulin receptor (IR) pre-mRNA undergoes alternative splicing that produces two isoforms, IR-A and IR-B. The ratio of IR-A to IR-B varies among tissues, which strongly suggests that IR mRNA alternative splicing is regulated in a tissue-specific manner. However, the precise molecular mechanism for IR alternative splicing remains to be elucidated, especially in liver. In this study, we have analyzed IR alternative splicing mechanism by preparing a mini-gene splicing reporter with rat genomic DNA. The splicing reporter that contains exon 11 and its flanking intronic sequences could reproduce alternative splicing pattern in rat hepatoma H4IIE cells. Introducing several deletions in introns of the reporter revealed that intron 11 contains the region near exon 11 essential to promote exon 11 inclusion. This region contains an UGCAUG sequence, a specific binding site for the Rbfox splicing regulator, and mutation in this sequence results in exon 11 skipping. Furthermore, RbFox2 knockdown in H4IIE cells enhanced exon 11 skipping of endogenous IR pre-mRNA. Lastly mutations in the SRSF3 binding site of exon11 together with the Rbfox2 binding site completely abolished exon 11 inclusion with a mini-gene reporter pre-mRNA. Our results indicate that RbFox2 and SRSF3 proteins mediate exon 11 inclusion in rat hepatoma cells.

RNA-PROTACs: Degraders of RNA-Binding Proteins

Defects in the functions of RNA binding proteins (RBPs) are at the origin of many diseases; however, targeting RBPs with conventional drugs has proven difficult. PROTACs are a new class of drugs that mediate selective degradation of a target protein through a cell's ubiquitination machinery. PROTACs comprise a moiety that binds the selected protein, conjugated to a ligand of an E3 ligase. Herein, we introduce RNA-PROTACs as a new concept in the targeting of RBPs. These chimeric structures employ small RNA mimics as targeting groups that dock the RNA-binding site of the RBP, whereupon a conjugated E3-recruiting peptide derived from the HIF-1α protein directs the RBP for proteasomal degradation. We performed a proof-of-concept demonstration with the degradation of two RBPs-a stem cell factor LIN28 and a splicing factor RBFOX1-and showed their use in cancer cell lines. The RNA-PROTAC approach opens the way to rapid, selective targeting of RBPs in a rational and general fashion.

Toxoplasma gondii Chinese I genotype Wh6 strain infection induces tau phosphorylation via activating GSK3β and causes hippocampal neuron apoptosis

Toxoplasma gondii (T. gondii) is a neurophilic and intracellular parasite that can affect plenty of vertebrate animals, including humans. Recent researches indicate that T. gondii infection is associated with neurodegenerative diseases such as Alzheimer's disease(AD). In addition, tau hyper-phosphorylation is a crucial event leading to the formation of nerve fiber tangles in AD. Despite the efforts to understand the interactions between T. gondii and AD, there are no clear results available so far. Here, we infected mice with the T. gondii of the Chinese 1 genotype Wh6 strain (TgCtwh6) for 60 days. Then we observed the formation of tissue cysts in the brain, the damage of neuron and the increased expression of phosphorylated tau (p-tau) in the hippocampal tissue of the mice. Similarly, we also found that p-tau, glycogen synthase kinase 3 beta (GSK3β), and phosphorylated GSK3β (p-GSK3β) were upregulated in vitro in TgCtwh6 challenged hippocampal neuron cell strain, HT22 cells. We noted a down-regulated expression of GSK3β,p-GSK3β, and p-tau in HT22 cells, which were pretreated with LiCl, an inhibitor of GSK3β. These data suggested that p-GSK3β may mediate tau phosphorylation after TgCtwh6 infection. Furthermore, TgCtwh6 infection also caused the increased expression of Bax and Caspase3, the decreased expression of Bcl-XL in HT22 cells, which had both early and late apoptosis. In all, our results indicated that TgCtwh6 infection not only led to phosphorylation of tau via activating GSK3β but also promoted hippocampal neuron apoptosis. Our research may partially reveal the mechanism with which TgCtwh6 induce neurofibrillary pathology.

Concentration-dependent splicing is enabled by Rbfox motifs of intermediate affinity

The Rbfox family of splicing factors regulate alternative splicing during animal development and in disease, impacting thousands of exons in the maturing brain, heart, and muscle. Rbfox proteins have long been known to bind to the RNA sequence GCAUG with high affinity, but just half of Rbfox binding sites contain a GCAUG motif in vivo. We incubated recombinant RBFOX2 with over 60,000 mouse and human transcriptomic sequences to reveal substantial binding to several moderate-affinity, non-GCAYG sites at a physiologically relevant range of RBFOX concentrations. We find that many of these “secondary motifs” bind Rbfox robustly in cells and that several together can exert regulation comparable to GCAUG in a trichromatic splicing reporter assay. Furthermore, secondary motifs regulate RNA splicing in neuronal development and in neuronal subtypes where cellular Rbfox concentrations are highest, enabling a second wave of splicing changes as Rbfox levels increase.

Cleft lip and cleft palate in Esrp1 knockout mice is associated with alterations in epithelial-mesenchymal crosstalk

Cleft lip is one of the most common human birth defects. However, there remain a limited number of mouse models of cleft lip that can be leveraged to characterize the genes and mechanisms that cause this disorder. Crosstalk between epithelial and mesenchymal cells underlies formation of the face and palate, but the basic molecular events mediating this crosstalk remain poorly understood. We previously demonstrated that mice lacking the epithelial-specific splicing factor Esrp1 have fully penetrant bilateral cleft lip and palate. In this study, we further investigated the mechanisms leading to cleft lip as well as cleft palate in both existing and new Esrp1 mutant mouse models. These studies included a detailed transcriptomic analysis of changes in ectoderm and mesenchyme in Esrp1-/- embryos during face formation. We identified altered expression of genes previously implicated in cleft lip and/or palate, including components of multiple signaling pathways. These findings provide the foundation for detailed investigations using Esrp1 mutant disease models to examine gene regulatory networks and pathways that are essential for normal face and palate development - the disruption of which leads to orofacial clefting in human patients.

An Intricate Connection between Alternative Splicing and Phenotypic Plasticity in Development and Cancer

During tumor progression, hypoxia, nutrient deprivation or changes in the extracellular environment (i.e., induced by anti-cancer drugs) elicit adaptive responses in cancer cells. Cellular plasticity increases the chance that tumor cells may survive in a challenging microenvironment, acquire new mechanisms of resistance to conventional drugs, and spread to distant sites. Re-activation of stem pathways appears as a significant cause of cellular plasticity because it promotes the acquisition of stem-like properties through a profound phenotypic reprogramming of cancer cells. In addition, it is a major contributor to tumor heterogeneity, depending on the coexistence of phenotypically distinct subpopulations in the same tumor bulk. Several cellular mechanisms may drive this fundamental change, in particular, high-throughput sequencing technologies revealed a key role for alternative splicing (AS). Effectively, AS is one of the most important pre-mRNA processes that increases the diversity of transcriptome and proteome in a tissue- and development-dependent manner. Moreover, defective AS has been associated with several human diseases. However, its role in cancer cell plasticity and tumor heterogeneity remains unclear. Therefore, unravelling the intricate relationship between AS and the maintenance of a stem-like phenotype may explain molecular mechanisms underlying cancer cell plasticity and improve cancer diagnosis and treatment.

Rbfox3/NeuN Regulates Alternative Splicing of Tau Exon 10

Alternative splicing of tau exon 10 generates tau isoforms with three or four microtubule-binding repeats, 3R-tau or 4R-tau, which are under developmental regulation. Dysregulation of tau exon 10 splicing is sufficient to cause neurodegenerative disorders. The RNA-binding Fox3 (Rbfox3), identified as NeuN, regulates RNA processing. However, whether Rbfox3/NeuN regulates tau exon 10 splicing is unknown. In the present study, we found that the developmental expression of 4R-tau coincided with the expression of Rbfox3 in rat brains. Rbfox3 enhanced tau exon 10 inclusion. Tau intron 10 contains UGCAUG, the conservative binding sequence of Rbfox3. Intron 10 of tau pre-mRNA was co-immunoprecipitated by Rbfox3/NeuN. Deletion mutants of the RNA recognition motif (RRM) or three RNA-binding sites of the RRM in Rbfox3/NeuN failed to enhance tau exon 10 inclusion. Rbfox3, specifically expressed in the fetal brain, did not affect tau exon 10 splicing. The level of Rbfox3/NeuN was reduced and was associated with the ratio of 4R-tau/3R-tau in the excitotoxic mouse brains induced by kainic acid. These findings suggest that Rbfox3/NeuN regulates the alternative splicing of tau exon 10 and that decreased Rbfox3/NeuN may lower the ratio of 4R-tau/3R-tau.

Modeling RNA-Binding Protein Specificity In Vivo by Precisely Registering Protein-RNA Crosslink Sites

RNA-binding proteins (RBPs) regulate post-transcriptional gene expression by recognizing short and degenerate sequence motifs in their target transcripts, but precisely defining their binding specificity remains challenging. Crosslinking and immunoprecipitation (CLIP) allows for mapping of the exact protein-RNA crosslink sites, which frequently reside at specific positions in RBP motifs at single-nucleotide resolution. Here, we have developed a computational method, named mCross, to jointly model RBP binding specificity while precisely registering the crosslinking position in motif sites. We applied mCross to 112 RBPs using ENCODE eCLIP data and validated the reliability of the discovered motifs by genome-wide analysis of allelic binding sites. Our analyses revealed that the prototypical SR protein SRSF1 recognizes clusters of GGA half-sites in addition to its canonical GGAGGA motif. Therefore, SRSF1 regulates splicing of a much larger repertoire of transcripts than previously appreciated, including HNRNPD and HNRNPDL, which are involved in multivalent protein assemblies and phase separation.

A slow transcription rate causes embryonic lethality and perturbs kinetic coupling of neuronal genes

The rate of RNA polymerase II (RNAPII) elongation has an important role in the control of alternative splicing (AS); however, the in vivo consequences of an altered elongation rate are unknown. Here, we generated mouse embryonic stem cells (ESCs) knocked in for a slow elongating form of RNAPII We show that a reduced transcriptional elongation rate results in early embryonic lethality in mice. Focusing on neuronal differentiation as a model, we observed that slow elongation impairs development of the neural lineage from ESCs, which is accompanied by changes in AS and in gene expression along this pathway. In particular, we found a crucial role for RNAPII elongation rate in transcription and splicing of long neuronal genes involved in synapse signaling. The impact of the kinetic coupling of RNAPII elongation rate with AS is greater in ESC-differentiated neurons than in pluripotent cells. Our results demonstrate the requirement for an appropriate transcriptional elongation rate to ensure proper gene expression and to regulate AS during development.

Bioinformatics challenges and perspectives when studying the effect of epigenetic modifications on alternative splicing

It is widely known that epigenetic modifications are important in regulating transcription, but several have also been reported in alternative splicing. The regulation of pre-mRNA splicing is important to explain proteomic diversity and the misregulation of splicing has been implicated in many diseases. Here, we give a brief overview of the role of epigenetics in alternative splicing and disease. We then discuss the bioinformatics methods that can be used to model interactions between epigenetic marks and regulators of splicing. These models can be used to identify alternative splicing and epigenetic changes across different phenotypes.

This article is part of a discussion meeting issue ‘Frontiers in epigenetic chemical biology’.

Specific inhibition of splicing factor activity by decoy RNA oligonucleotides

Alternative splicing, a fundamental step in gene expression, is deregulated in many diseases. Splicing factors (SFs), which regulate this process, are up- or down regulated or mutated in several diseases including cancer. To date, there are no inhibitors that directly inhibit the activity of SFs. We designed decoy oligonucleotides, composed of several repeats of a RNA motif, which is recognized by a single SF. Here we show that decoy oligonucleotides targeting splicing factors RBFOX1/2, SRSF1 and PTBP1, can specifically bind to their respective SFs and inhibit their splicing and biological activities both in vitro and in vivo. These decoy oligonucleotides present an approach to specifically downregulate SF activity in conditions where SFs are either up-regulated or hyperactive.

Alternative splicing, critical for gene expression, is deregulated in many diseases. Here the authors develop decoy oligonucleotides to specifically downregulate splicing factors activity.

An RNA Switch of a Large Exon of Ninein Is Regulated by the Neural Stem Cell Specific-RNA Binding Protein, Qki5

A set of tissue-specific splicing factors are thought to govern alternative splicing events during neural progenitor cell (NPC)-to-neuron transition by regulating neuron-specific exons. Here, we propose one such factor, RNA-binding protein Quaking 5 (Qki5), which is specifically expressed in the early embryonic neural stem cells. We performed mRNA-SEQ (Sequence) analysis using mRNAs obtained by developing cerebral cortices in Qk (Quaking) conditional knockout (cKO) mice. As expected, we found a large number of alternative splicing changes between control and conditional knockouts relative to changes in transcript levels. DAVID (The Database for Annotation, Visualization and Integrated Discovery) and Metascape analyses suggested that the affected spliced genes are involved in axon development and microtubule-based processes. Among these, the mRNA coding for the Ninein protein is listed as one of Qki protein-dependent alternative splicing targets. Interestingly, this exon encodes a very long polypeptide (2121 nt), and has been previously defined as a dynamic RNA switch during the NPC-to-neuron transition. Additionally, we validated that the regulation of this large exon is consistent with the Qki5-dependent alternative exon inclusion mode suggested by our previous Qki5 HITS-CLIP (high throughput sequencing-cross linking immunoprecipitation) analysis. Taken together, these data suggest that Qki5 is an important factor for alternative splicing in the NPC-to-neuron transition.

Subphenotype meta-analysis of testicular cancer genome-wide association study data suggests a role for RBFOX family genes in cryptorchidism susceptibility

STUDY QUESTION

Can subphenotype analysis of genome-wide association study (GWAS) data from subjects with testicular germ cell tumor (TGCT) provide insight into cryptorchidism (undescended testis, UDT) susceptibility?

SUMMARY ANSWER

Suggestive intragenic GWAS signals common to UDT, TGCT case-case and TGCT case-control analyses occur in genes encoding RBFOX RNA-binding proteins (RBPs) and their neurodevelopmental targets.

WHAT IS KNOWN ALREADY

UDT is a strong risk factor for TGCT, but while genetic risk factors for TGCT are well-known, genetic susceptibility to UDT is poorly understood and appears to be more complex.

STUDY DESIGN, SIZE, DURATION

We performed a secondary subphenotype analysis of existing GWAS data from the Testicular Cancer Consortium (TECAC) and compared these results with our previously published UDT GWAS data, and with data previously acquired from studies of the fetal rat gubernaculum.

PARTICIPANTS/MATERIALS, SETTING, METHODS

Studies from the National Cancer Institute (NCI), United Kingdom (UK) and University of Pennsylvania (Penn) that enrolled white subjects were the source of the TGCT GWAS data. We completed UDT subphenotype case-case (TGCT/UDT vs TGCT/non-UDT) and case-control (TGCT/UDT vs control), collectively referred to as 'TECAC' analyses, followed by a meta-analysis comprising 129 TGCT/UDT cases, 1771 TGCT/non-UDT cases, and 3967 unaffected controls. We reanalyzed our UDT GWAS results comprising 844 cases and 2718 controls by mapping suggestive UDT and TECAC signals (defined as P < 0.001) to genes using Ingenuity Pathway Analysis (IPA®). We compared associated pathways and enriched gene categories common to all analyses after Benjamini-Hochberg multiple testing correction, and analyzed transcript levels and protein expression using qRT-PCR and rat fetal gubernaculum confocal imaging, respectively.

MAIN RESULTS AND THE ROLE OF CHANCE

We found suggestive signals within 19 genes common to all three analyses, including RBFOX1 and RBFOX3, neurodevelopmental paralogs that encode RBPs targeting (U)GCATG-containing transcripts. Ten of the 19 genes participate in neurodevelopment and/or contribute to risk of neurodevelopmental disorders. Experimentally predicted RBFOX gene targets were strongly overrepresented among suggestive intragenic signals for the UDT (117 of 628 (19%), P = 3.5 × 10-24), TECAC case-case (129 of 711 (18%), P = 2.5 × 10-27) and TECAC case-control (117 of 679 (17%), P = 2 × 10-21) analyses, and a majority of the genes common to all three analyses (12 of 19 (63%), P = 3 × 10-9) are predicted RBFOX targets. Rbfox1, Rbfox2 and their encoded proteins are expressed in the rat fetal gubernaculum. Predicted RBFOX targets are also enriched among transcripts differentially regulated in the fetal gubernaculum during normal development (P = 3 × 10-31), in response to in vitro hormonal stimulation (P = 5 × 10-45) and in the cryptorchid LE/orl rat (P = 2 × 10-42).

LARGE SCALE DATA

GWAS data included in this study are available in the database of Genotypes and Phenotypes (dbGaP accession numbers phs000986.v1.p1 and phs001349.v1p1).

LIMITATIONS, REASONS FOR CAUTION

These GWAS data did not reach genome-wide significance for any individual analysis. UDT appears to have a complex etiology that also includes environmental factors, and such complexity may require much larger sample sizes than are currently available. The current methodology may also introduce bias that favors false discovery of larger genes.

WIDER IMPLICATIONS OF THE FINDINGS

Common suggestive intragenic GWAS signals suggest that RBFOX paralogs and other neurodevelopmental genes are potential UDT risk candidates, and potential TGCT susceptibility modifiers. Enrichment of predicted RBFOX targets among differentially expressed transcripts in the fetal gubernaculum additionally suggests a role for this RBP family in regulation of testicular descent. As RBFOX proteins regulate alternative splicing of Calca to generate calcitonin gene-related peptide, a protein linked to development and function of the gubernaculum, additional studies that address the role of these proteins in UDT are warranted.

STUDY FUNDING/COMPETING INTEREST(S)

The Eunice Kennedy Shriver National Institute for Child Health and Human Development (R01HD060769); National Center for Research Resources (P20RR20173), National Institute of General Medical Sciences (P20GM103464), Nemours Biomedical Research, the Testicular Cancer Consortium (U01CA164947), the Intramural Research Program of the NCI, a support services contract HHSN26120130003C with IMS, Inc., the Abramson Cancer Center at Penn, National Cancer Institute (CA114478), the Institute of Cancer Research, UK and the Wellcome Trust Case-Control Consortium (WTCCC) 2. None of the authors reports a conflict of interest.

TGF-β-induced alternative splicing of TAK1 promotes EMT and drug resistance

Transforming growth factor-β (TGF-β) is major inducer of epithelial to mesenchymal transition (EMT), which associates with cancer cell metastasis and resistance to chemotherapy and targeted drugs, through both transcriptional and non-transcriptional mechanisms. We previously reported that in cancer cells, heightened mitogenic signaling allows TGF-β-activated Smad3 to interact with poly(RC) binding protein 1 (PCBP1) and together they regulate many alternative splicing events that favors expression of protein isoforms essential for EMT, cytoskeletal rearrangement, and adherens junction signaling. Here, we show that the exclusion of TGF-β-activated kinase 1 (TAK1) variable exon 12 requires another RNA-binding protein, Fox-1 homolog 2 (Rbfox2), which binds intronic sequences in front of exon 12 independently of the Smad3-PCBP1 complex. Functionally, exon 12-excluded TAK1∆E12 and full length TAK1FL are distinct. The short isoform TAK1∆E12 is constitutively active and supports TGF-β-induced EMT and nuclear factor kappa B (NF-κB) signaling, whereas the full-length isoform TAK1FL promotes TGF-β-induced apoptosis. These observations offer a harmonious explanation for how a single TAK1 kinase can mediate the opposing responses of cell survival and apoptosis in response to TGF-β. They also reveal a propensity of the alternatively spliced TAK1 isoform TAK1∆E12 to cause drug resistance due to its activity in supporting EMT and NF-κB survival signaling.

Mapping of In Vivo RNA-Binding Sites by Ultraviolet (UV)-Cross-Linking Immunoprecipitation (CLIP)

RNA "CLIP" (cross-linking immunoprecipitation), the method by which RNA-protein complexes are covalently cross-linked and purified and the RNA sequenced, has attracted attention as a powerful means of developing genome-wide maps of direct, functional RNA-protein interaction sites. These maps have been used to identify points of regulation, and they hold promise for understanding the dynamics of RNA regulation in normal cell function and its dysregulation in disease.