The splicing regulatory element, UGCAUG, is phylogenetically and spatially conserved in introns that flank tissue-specific alternative exons

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.

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.

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.

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.

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.

Unannotated splicing regulatory elements in deep intron space

Deep intron space harbors a diverse array of splicing regulatory elements that cooperate with better-known exon-proximal elements to enforce proper tissue-specific and development-specific pre-mRNA processing. Many deep intron elements have been highly conserved through vertebrate evolution, yet remain poorly annotated in the human genome. Recursive splicing exons (RS-exons) and intraexons promote noncanonical, multistep resplicing pathways in long introns, involving transient intermediate structures that are greatly underrepresented in RNA-seq datasets. Decoy splice sites and decoy exons act at a distance to inhibit splicing catalysis at annotated splice sites, with functional consequences such as exon skipping and intron retention. RNA:RNA bridges can juxtapose distant sequences within or across introns to activate deep intron splicing enhancers and silencers, to loop out exons to be skipped, or to select one member of a mutually exclusive set of exons. Similarly, protein bridges mediated by interactions among transcript-bound RNA binding proteins (RBPs) can modulate splicing outcomes. Experimental disruption of deep intron elements serving any of these functions can abrogate normal splicing, strongly suggesting that natural mutations of deep intron elements can do likewise to cause human disease. Understanding noncanonical splicing pathways and discovering deep intron regulatory signals, many of which map hundreds to many thousands of nucleotides from annotated splice junctions, is of great academic interest for basic scientists studying alternative splicing mechanisms. Hopefully, this knowledge coupled with increased analysis of deep intron sequences will also have important medical applications, as better interpretation of deep intron mutations may reveal new disease mechanisms and suggest new therapies. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing.

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.

Using secondary structure to predict the effects of genetic variants on alternative splicing

Accurate interpretation of genomic variants that alter RNA splicing is critical to precision medicine. We present a computational framework, Prediction of variant Effect on Percent Spliced In (PEPSI), that predicts the splicing impact of coding and noncoding variants for the Fifth Critical Assessment of Genome Interpretation (CAGI5) "Vex-seq" challenge. PEPSI is a random forest regression model trained on multiple layers of features associated with sequence conservation and regulatory sequence elements. Compared to other splicing defect prediction tools from the literature, our framework integrates secondary structure information in predicting variants that disrupt splicing regulatory elements (SREs). We applied our model to classify splice-disrupting variants among 2,094 single-nucleotide polymorphisms from the Exome Aggregation Consortium using model-predicted changes in percent spliced in (ΔPSI) associated with tested variants. Benchmarking our model against widely used state-of-the-art tools, we demonstrate that PEPSI achieves comparable performance in terms of sensitivity and precision. Moreover, we also show that using secondary structure context can help resolve several cases where changes in the counts of SREs do not correspond with the directionality of ΔPSI measured for tested variants.

Prader-Willi locus Snord116 RNA processing requires an active endogenous allele and neuron-specific splicing by Rbfox3/NeuN

Prader-Willi syndrome (PWS), an imprinted neurodevelopmental disorder characterized by metabolic, sleep and neuropsychiatric features, is caused by the loss of paternal SNORD116, containing only non-coding RNAs (ncRNAs). The primary SNORD116 transcript is processed into small nucleolar RNAs (snoRNAs), which localize to nucleoli, and their spliced host gene 116HG, which is retained at its site of transcription. While functional complementation of the SNORD116 ncRNAs is a desirable goal for treating PWS, the mechanistic requirements of SNORD116 RNA processing are poorly understood. Here we developed and tested a novel transgenic mouse which ubiquitously expresses Snord116 on both a wild-type and a Snord116 paternal deletion (Snord116+/-) background. Interestingly, while the Snord116 transgene was ubiquitously expressed in multiple tissues, splicing of the transgene and production of snoRNAs was limited to brain tissues. Knockdown of Rbfox3, encoding neuron-specific splicing factor neuronal nuclei (NeuN) in Snord116+/--derived neurons, reduced splicing of the transgene in neurons. RNA fluorescence in situ hybridization for 116HG revealed a single significantly larger signal in transgenic mice, demonstrating colocalization of transgenic and endogenous 116HG RNAs. Similarly, significantly increased snoRNA levels were detected in transgenic neuronal nucleoli, indicating that transgenic Snord116 snoRNAs were effectively processed and localized. In contrast, neither transgenic 116HG nor snoRNAs were detectable in either non-neuronal tissues or Snord116+/- neurons. Together, these results demonstrate that exogenous expression and neuron-specific splicing of the Snord116 locus are insufficient to rescue the genetic deficiency of Snord116 paternal deletion. Elucidating the mechanisms regulating Snord116 processing and localization is essential to develop effective gene replacement therapies for PWS.

Rbfox family proteins make the homo- and hetero-oligomeric complexes

Rbfox family of proteins that consists of Rbfox1, Rbfox2, and Rbfox3 in mammals regulates alternative pre-mRNA splicing in various tissues via direct binding to their RNA binding element. Although many studies have indicated the splicing activity of each member of the Rbfox family, the interactions of Rbfox family proteins are largely unknown. Here, we have investigated interactions among Rbfox family proteins. Co-immunoprecipitation (Co-IP) and GST-pull down assays confirmed that Rbfox proteins form homo and hetero complexes. Moreover, in vivo crosslinking using disuccinimidyl suberate treatment indicated that the Rbfox proteins form a dimer which then assembles with other proteins to form a large multiprotein complex. Duolink in situ proximity ligation (PLA) assay revealed that neuron specific Rbfox3 protein interacts with other Rbfox family proteins. This study is the first to provide an evidence that Rbfox family proteins form homo- and hetero-oligomeric complexes in vivo.

Cell-Specific RNA Binding Protein Rbfox2 Regulates CaV2.2 mRNA Exon Composition and CaV2.2 Current Size

The majority of multiexon mammalian genes contain alternatively spliced exons that have unique expression patterns in different cell populations and that have important cell functions. The expression profiles of alternative exons are controlled by cell-specific splicing factors that can promote exon inclusion or exon skipping but with few exceptions we do not know which specific splicing factors control the expression of alternatively spliced exons of known biological function. Many ion channel genes undergo extensive alternative splicing including Cacna1b that encodes the voltage-gated CaV2.2 α1 subunit. Alternatively spliced exon 18a in Cacna1b RNA encodes 21 amino acids in the II-III loop of CaV2.2, and its expression differs across the nervous system and over development. Genome-wide, protein-RNA binding analyses coupled to high-throughput RNA sequencing show that RNA binding Fox (Rbfox) proteins associate with CaV2.2 (Cacna1b) pre-mRNAs. Here, we link Rbfox2 to suppression of e18a. We show increased e18a inclusion in CaV2.2 mRNAs: (1) after siRNA knockdown of Rbfox2 in a neuronal cell line and (2) in RNA from sympathetic neurons of adult compared to early postnatal mice. By immunoprecipitation of Rbfox2-RNA complexes followed by qPCR, we demonstrate reduced Rbfox2 binding upstream of e18a in RNA from sympathetic neurons of adult compared to early postnatal mice. CaV2.2 currents in cell lines and in sympathetic neurons expressing only e18a-CaV2.2 are larger compared to currents from those expressing only Δ18a-CaV2.2. We conclude that Rbfox2 represses e18a inclusion during pre-mRNA splicing of CaV2.2, limiting the size of CaV2.2 currents early in development in certain neuronal populations.

A saga of cancer epigenetics: linking epigenetics to alternative splicing

The discovery of an increasing number of alternative splicing events in the human genome highlighted that ∼94% of genes generate alternatively spliced transcripts that may produce different protein isoforms with diverse functions. It is now well known that several diseases are a direct and indirect consequence of aberrant splicing events in humans. In addition to the conventional mode of alternative splicing regulation by 'cis' RNA-binding sites and 'trans' RNA-binding proteins, recent literature provides enormous evidence for epigenetic regulation of alternative splicing. The epigenetic modifications may regulate alternative splicing by either influencing the transcription elongation rate of RNA polymerase II or by recruiting a specific splicing regulator via different chromatin adaptors. The epigenetic alterations and aberrant alternative splicing are known to be associated with various diseases individually, but this review discusses/highlights the latest literature on the role of epigenetic alterations in the regulation of alternative splicing and thereby cancer progression. This review also points out the need for further studies to understand the interplay between epigenetic modifications and aberrant alternative splicing in cancer progression.

Satellite-like cells contribute to pax7-dependent skeletal muscle repair in adult zebrafish

Satellite cells, also known as muscle stem cells, are responsible for skeletal muscle growth and repair in mammals. Pax7 and Pax3 transcription factors are established satellite cell markers required for muscle development and regeneration, and there is great interest in identifying additional factors that regulate satellite cell proliferation, differentiation, and/or skeletal muscle regeneration. Due to the powerful regenerative capacity of many zebrafish tissues, even in adults, we are exploring the regenerative potential of adult zebrafish skeletal muscle. Here, we show that adult zebrafish skeletal muscle contains cells similar to mammalian satellite cells. Adult zebrafish satellite-like cells have dense heterochromatin, express Pax7 and Pax3, proliferate in response to injury, and show peak myogenic responses 4-5 days post-injury (dpi). Furthermore, using a pax7a-driven GFP reporter, we present evidence implicating satellite-like cells as a possible source of new muscle. In lieu of central nucleation, which distinguishes regenerating myofibers in mammals, we describe several characteristics that robustly identify newly-forming myofibers from surrounding fibers in injured adult zebrafish muscle. These characteristics include partially overlapping expression in satellite-like cells and regenerating myofibers of two RNA-binding proteins Rbfox2 and Rbfoxl1, known to regulate embryonic muscle development and function. Finally, by analyzing pax7a; pax7b double mutant zebrafish, we show that Pax7 is required for adult skeletal muscle repair, as it is in the mouse.

Developmental regulation of RNA processing by Rbfox proteins

The Rbfox genes encode an ancient family of sequence-specific RNA binding proteins (RBPs) that are critical developmental regulators in multiple tissues including skeletal muscle, cardiac muscle, and brain. The hallmark of Rbfox proteins is a single high-affinity RRM domain, highly conserved from insects to humans, that binds preferentially to UGCAUG motifs at diverse regulatory sites in pre-mRNA introns, mRNA 3'UTRs, and pre-miRNAs hairpin structures. Versatile regulatory circuits operate on Rbfox pre-mRNA and mRNA to ensure proper expression of Rbfox1 protein isoforms, which then act on the broader transcriptome to regulate alternative splicing networks, mRNA stability and translation, and microRNA processing. Complex Rbfox expression is encoded in large genes encompassing multiple promoters and alternative splicing options that govern spatiotemporal expression of structurally distinct and tissue-specific protein isoforms with different classes of RNA targets. Nuclear Rbfox1 is a candidate master regulator that binds intronic UGCAUG elements to impact splicing efficiency of target alternative exons, many in transcripts for other splicing regulators. Tissue-specificity of Rbfox-mediated alternative splicing is executed by combinatorial regulation through the integrated activity of Rbfox proteins and synergistic or antagonistic splicing factors. Studies in animal models show that Rbfox1-related genes are critical for diverse developmental processes including germ cell differentiation and memory in Drosophila, neuronal migration and function in mouse brain, myoblast fusion and skeletal muscle function, and normal heart function. Finally, genetic and biochemical evidence suggest that aberrations in Rbfox-regulated circuitry are risk factors for multiple human disorders, especially neurodevelopmental disorders including epilepsy and autism, and cardiac hypertrophy. WIREs RNA 2017, 8:e1398. doi: 10.1002/wrna.1398 For further resources related to this article, please visit the WIREs website.

Rbfox2 function in RNA metabolism is impaired in hypoplastic left heart syndrome patient hearts

Hypoplastic left heart syndrome (HLHS) is a fatal congenital heart disease in which the left side of the heart is underdeveloped, impairing the systemic circulation. Underdeveloped left ventricle exerts biomechanical stress on the right ventricle that can progress into heart failure. Genome-wide transcriptome changes have been identified at early stages in the right ventricle (RV) of infants with HLHS, although the molecular mechanisms remain unknown. Here, we demonstrate that the RNA binding protein Rbfox2, which is mutated in HLHS patients, is a contributor to transcriptome changes in HLHS patient RVs. Our results indicate that majority of transcripts differentially expressed in HLHS patient hearts have validated Rbfox2 binding sites. We show that Rbfox2 regulates mRNA levels of targets with 3’UTR binding sites contributing to aberrant gene expression in HLHS patients. Strikingly, the Rbfox2 nonsense mutation identified in HLHS patients truncates the protein, impairs its subcellular distribution and adversely affects its function in RNA metabolism. Overall, our findings uncover a novel role for Rbfox2 in controlling transcriptome in HLHS.

Splicing predictions reliably classify different types of alternative splicing

Alternative splicing is a key player in the creation of complex mammalian transcriptomes and its misregulation is associated with many human diseases. Multiple mRNA isoforms are generated from most human genes, a process mediated by the interplay of various RNA signature elements and trans-acting factors that guide spliceosomal assembly and intron removal. Here, we introduce a splicing predictor that evaluates hundreds of RNA features simultaneously to successfully differentiate between exons that are constitutively spliced, exons that undergo alternative 5' or 3' splice-site selection, and alternative cassette-type exons. Surprisingly, the splicing predictor did not feature strong discriminatory contributions from binding sites for known splicing regulators. Rather, the ability of an exon to be involved in one or multiple types of alternative splicing is dictated by its immediate sequence context, mainly driven by the identity of the exon's splice sites, the conservation around them, and its exon/intron architecture. Thus, the splicing behavior of human exons can be reliably predicted based on basic RNA sequence elements.

HITS-CLIP and integrative modeling define the Rbfox splicing-regulatory network linked to brain development and autism

The RNA binding proteins Rbfox1/2/3 regulate alternative splicing in the nervous system, and disruption of Rbfox1 has been implicated in autism. However, comprehensive identification of functional Rbfox targets has been challenging. Here, we perform HITS-CLIP for all three Rbfox family members in order to globally map, at a single-nucleotide resolution, their in vivo RNA interaction sites in the mouse brain. We find that the two guanines in the Rbfox binding motif UGCAUG are critical for protein-RNA interactions and crosslinking. Using integrative modeling, these interaction sites, combined with additional datasets, define 1,059 direct Rbfox target alternative splicing events. Over half of the quantifiable targets show dynamic changes during brain development. Of particular interest are 111 events from 48 candidate autism-susceptibility genes, including syndromic autism genes Shank3, Cacna1c, and Tsc2. Alteration of Rbfox targets in some autistic brains is correlated with downregulation of all three Rbfox proteins, supporting the potential clinical relevance of the splicing-regulatory network.

Rbfox proteins regulate alternative mRNA splicing through evolutionarily conserved RNA bridges

Alternative splicing (AS) enables programmed diversity of gene expression across tissues and development. We show here that binding in distal intronic regions (>500 nucleotides (nt) from any exon) by Rbfox splicing factors important in development is extensive and is an active mode of splicing regulation. Similarly to exon-proximal sites, distal sites contain evolutionarily conserved GCATG sequences and are associated with AS activation and repression upon modulation of Rbfox abundance in human and mouse experimental systems. As a proof of principle, we validated the activity of two specific Rbfox enhancers in KIF21A and ENAH distal introns and showed that a conserved long-range RNA-RNA base-pairing interaction (an RNA bridge) is necessary for Rbfox-mediated exon inclusion in the ENAH gene. Thus we demonstrate a previously unknown RNA-mediated mechanism for AS control by distally bound RNA-binding proteins.

Control of neuronal voltage-gated calcium ion channels from RNA to protein

Voltage-gated calcium ion (CaV) channels convert neuronal activity into rapid intracellular calcium signals to trigger a myriad of cellular responses. Their involvement in major neurological and psychiatric diseases, and importance as therapeutic targets, has propelled interest in subcellular-specific mechanisms that align CaV channel activity to specific tasks. Here, we highlight recent studies that delineate mechanisms controlling the expression of CaV channels at the level of RNA and protein. We discuss the roles of RNA editing and alternative pre-mRNA splicing in generating CaV channel isoforms with activities specific to the demands of individual cells; the roles of ubiquitination and accessory proteins in regulating CaV channel expression; and the specific binding partners that contribute to both pre- and postsynaptic CaV channel function.

Variants affecting exon skipping contribute to complex traits

DNA variants that affect alternative splicing and the relative quantities of different gene transcripts have been shown to be risk alleles for some Mendelian diseases. However, for complex traits characterized by a low odds ratio for any single contributing variant, very few studies have investigated the contribution of splicing variants. The overarching goal of this study is to discover and characterize the role that variants affecting alternative splicing may play in the genetic etiology of complex traits, which include a significant number of the common human diseases. Specifically, we hypothesize that single nucleotide polymorphisms (SNPs) in splicing regulatory elements can be characterized in silico to identify variants affecting splicing, and that these variants may contribute to the etiology of complex diseases as well as the inter-individual variability in the ratios of alternative transcripts. We leverage high-throughput expression profiling to 1) experimentally validate our in silico predictions of skipped exons and 2) characterize the molecular role of intronic genetic variations in alternative splicing events in the context of complex human traits and diseases. We propose that intronic SNPs play a role as genetic regulators within splicing regulatory elements and show that their associated exon skipping events can affect protein domains and structure. We find that SNPs we would predict to affect exon skipping are enriched among the set of SNPs reported to be associated with complex human traits.

Alternative splicing is a common eukaryotic cellular mechanism that allows for the production of multiple proteins from one gene and occurs in 40%–90% of all human genes. Alternative splicing has been shown to be important for many critical biological processes, including development, evolution, and even psychological behavior. Additionally, alternative splicing has been associated with 15%–50% of human genetic diseases, including breast cancer; however, the precise mechanism by which genetic variations regulate this process remains to be fully elucidated. In this study, we develop an integrative approach that utilizes sequence-based analysis and genome-wide expression profiling to identify genetic variations that may affect alternative splicing. We also evaluate their enrichment among established disease-associated variations. Our study provides insights into the functionality of these variations and emphasizes their importance for complex human traits and diseases.

Mutually exclusive splicing regulates the Nav 1.6 sodium channel function through a combinatorial mechanism that involves three distinct splicing regulatory elements and their ligands

Mutually exclusive splicing is a form of alternative pre-mRNA processing that consists in the use of only one of a set of two or more exons. We have investigated the mechanisms involved in this process for exon 18 of the Na(v) 1.6 sodium channel transcript and its significance regarding gene-expression regulation. The 18N exon (neonatal form) has a stop codon in phase and although the mRNA can be detected by amplification methods, the truncated protein has not been observed. The switch from 18N to 18A (adult form) occurs only in a restricted set of neural tissues producing the functional channel while other tissues display the mRNA with the 18N exon also in adulthood. We demonstrate that the mRNA species carrying the stop codon is subjected to Nonsense-Mediated Decay, providing a control mechanism of channel expression. We also map a string of cis-elements within the mutually exclusive exons and in the flanking introns responsible for their strict tissue and temporal specificity. These elements bind a series of positive (RbFox-1, SRSF1, SRSF2) and negative (hnRNPA1, PTB, hnRNPA2/B1, hnRNPD-like JKTBP) splicing regulatory proteins. These splicing factors, with the exception of RbFox-1, are ubiquitous but their levels vary during development and differentiation, ensuing unique sets of tissue and temporal levels of splicing factors. The combinatorial nature of these elements is highlighted by the dominance of the elements that bind the ubiquitous factors over the tissue specific RbFox-1.

Rare inherited A2BP1 deletion in a proband with autism and developmental hemiparesis

Ataxin 2 binding protein 1 (A2BP1 aka FOX1, RBFOX1) is an RNA binding protein responsible for regulation of pre-mRNA splicing events in a number of critical developmental genes expressed in muscle, heart and neuronal cells [Shibata et al. (2000); Mamm Genome 12:595-601; Jin et al. (2003); EMBO J 22:905-912; Underwood et al. (2005); Mol Cell Biol 25:10005-10016]. Rare copy number abnormalities of A2BP1 have been previously associated with cognitive impairment, attention deficit disorder and autism [Martin et al. (2007); Am J Med Gen Part B 144B:869-876; Elia et al. (2010); Mol Psychiatry 15:637-646.]. Using a 1M Illumina SNP microarray, we identified a 1.3 kb deletion in A2BP1, which was subsequently validated by quantitative PCR. Here we present an in depth case study of an individual with autism and mild developmental hemiparesis in whom the deletion was detected. This study provides further support for the possible role of rare copy number variants in A2BP1 in the development of autism and associated motor asymmetries.

VERSE: a varying effect regression for splicing elements discovery

Identification of splicing regulatory elements (SREs) deserves special attention because these cis-acting short sequences are vital parts of splicing code. The fact that a variety of other biological signals cooperatively govern the splicing pattern indicates the necessity of developing novel tools to incorporate information from multiple sources to improve splicing factor binding sites prediction. Under this context, we proposed a Varying Effect Regression for Splicing Elements (VERSE) to discover intronic SREs in the proximity of exon junctions by integrating other biological features. As a result, 1562 intronic SREs were identified in 16 human tissues, many of which overlapped with experimentally verified binding motifs for several well-known splicing factors, including FOX-1, PTB, hnRNP A/B, hnRNP F/H, and so on. The discovered tissue, region, and conservation preferences of the putative motifs demonstrate that splice site selection is a complicated process that needs subtle and delicate regulation. VERSE may serve as a powerful tool to not only discover SREs by incorporating additional informative signals but also precisely quantify their varying contribution under different biological contexts.

Common variation at BARD1 results in the expression of an oncogenic isoform that influences neuroblastoma susceptibility and oncogenicity

The mechanisms underlying genetic susceptibility at loci discovered by genome-wide association study (GWAS) approaches in human cancer remain largely undefined. In this study, we characterized the high-risk neuroblastoma association at the BRCA1-related locus, BARD1, showing that disease-associated variations correlate with increased expression of the oncogenically activated isoform, BARD1β. In neuroblastoma cells, silencing of BARD1β showed genotype-specific cytotoxic effects, including decreased substrate-adherence, anchorage-independence, and foci growth. In established murine fibroblasts, overexpression of BARD1β was sufficient for neoplastic transformation. BARD1β stabilized the Aurora family of kinases in neuroblastoma cells, suggesting both a mechanism for the observed effect and a potential therapeutic strategy. Together, our findings identify BARD1β as an oncogenic driver of high-risk neuroblastoma tumorigenesis, and more generally, they illustrate how robust GWAS signals offer genomic landmarks to identify molecular mechanisms involved in both tumor initiation and malignant progression. The interaction of BARD1β with the Aurora family of kinases lends strong support to the ongoing work to develop Aurora kinase inhibitors for clinically aggressive neuroblastoma.

RBFOX2 promotes protein 4.1R exon 16 selection via U1 snRNP recruitment

The erythroid differentiation-specific splicing switch of protein 4.1R exon 16, which encodes a spectrin/actin-binding peptide critical for erythrocyte membrane stability, is modulated by the differentiation-induced splicing factor RBFOX2. We have now characterized the mechanism by which RBFOX2 regulates exon 16 splicing through the downstream intronic element UGCAUG. Exon 16 possesses a weak 5' splice site (GAG/GTTTGT), which when strengthened to a consensus sequence (GAG/GTAAGT) leads to near-total exon 16 inclusion. Impaired RBFOX2 binding reduces exon 16 inclusion in the context of the native weak 5' splice site, but not the engineered strong 5' splice site, implying that RBFOX2 achieves its effect by promoting utilization of the weak 5' splice site. We further demonstrate that RBFOX2 increases U1 snRNP recruitment to the weak 5' splice site through direct interaction between its C-terminal domain (CTD) and the zinc finger region of U1C and that the CTD is required for the effect of RBFOX2 on exon 16 splicing. Our data suggest a novel mechanism for exon 16 5' splice site activation in which the binding of RBFOX2 to downstream intronic splicing enhancers stabilizes the pre-mRNA-U1 snRNP complex through interactions with U1C.

Association pattern mining of intron retention events in human based on hybrid learning machine

Alternative splicing is a main component of protein diversity, and aberrant splicing is known to be one of the main causes of genetic disorders such as cancer. Many statistical and computational approaches have identified several major factors that determine the splicing event, such as exon/intron length, splice site strength, and density of splicing enhancers or silencers. These factors may be correlated with one another and thus result in a specific type of splicing, but there has not been a systematic approach to extracting comprehensible association patterns. Here, we attempted to understand the decision making process of the learning machine on intron retention event. We adopted a hybrid learning machine approach using a random forest and association rule mining algorithm to determine the governing factors of intron retention events and their combined effect on decision-making processes. By quantifying all candidate features into five category values, we enhanced the understandability of generated rules. The interesting features found by the random forest algorithm are that only the adenine- and thymine-based triplets such as ATA, TTA, and ATT, but not the known intronic splicing enhancer GGG triplet is shown the significant features. The rules generated by the association rule mining algorithm also show that constitutive introns are generally characterized by high adenine- and thymine-based triplet frequency (level 3 and above), 3' and 5' splice site scores, exonic splicing silencer scores, and intron length, whereas retained introns are characterized by low-level counterpart scores.

Networking in a global world: establishing functional connections between neural splicing regulators and their target transcripts

Recent genome-wide analyses have indicated that almost all primary transcripts from multi-exon human genes undergo alternative pre-mRNA splicing (AS). Given the prevalence of AS and its importance in expanding proteomic complexity, a major challenge that lies ahead is to determine the functional specificity of isoforms in a cellular context. A significant fraction of alternatively spliced transcripts are regulated in a tissue- or cell-type-specific manner, suggesting that these mRNA variants likely function in the generation of cellular diversity. Complementary to these observations, several tissue-specific splicing factors have been identified, and a number of methodological advances have enabled the identification of large repertoires of target transcripts regulated by these proteins. An emerging theme is that tissue-specific splicing factors regulate coherent sets of splice variants in genes known to function in related biological pathways. This review focuses on the recent progress in our understanding of neural-specific splicing factors and their regulatory networks and outlines existing and emerging strategies for uncovering important biological roles for the isoforms that comprise these networks.

Synaptic signaling and aberrant RNA splicing in autism spectrum disorders

Interactions between presynaptic and postsynaptic cellular adhesion molecules (CAMs) drive synapse maturation during development. These trans-synaptic interactions are regulated by alternative splicing of CAM RNAs, which ultimately determines neurotransmitter phenotype. The diverse assortment of RNAs produced by alternative splicing generates countless protein isoforms necessary for guiding specialized cell-to-cell connectivity. Failure to generate the appropriate synaptic adhesion proteins is associated with disrupted glutamatergic and gamma-aminobutyric acid signaling, resulting in loss of activity-dependent neuronal plasticity, and risk for developmental disorders, including autism. While the majority of genetic mutations currently linked to autism are rare variants that change the protein-coding sequence of synaptic candidate genes, regulatory polymorphisms affecting constitutive and alternative splicing have emerged as risk factors in numerous other diseases, accounting for an estimated 40–60% of general disease risk. Here, we review the relationship between aberrant RNA splicing of synapse-related genes and autism spectrum disorders.

Systems analysis of alternative splicing and its regulation

Alternative splicing (AS) has emerged as a key mechanism that accounts for gene expression diversity in metazoan organisms. Splicing is tightly regulated by a repertoire of RNA and protein factors and RNA sequence elements that function in a cooperative manner. Systems-level experimental and computational approaches have been instrumental in establishing comprehensive profiles of transcript variants generated by AS. In addition, systems biology approaches are starting to define how combinatorial splicing regulation shapes the complex splicing phenotypes observed in different tissue types and developmental stages and under different conditions. Here, we review recent progress in these areas.

Model-based detection of alternative splicing signals

Motivation: Transcripts from ∼95% of human multi-exon genes are subject to alternative splicing (AS). The growing interest in AS is propelled by its prominent contribution to transcriptome and proteome complexity and the role of aberrant AS in numerous diseases. Recent technological advances enable thousands of exons to be simultaneously profiled across diverse cell types and cellular conditions, but require accurate identification of condition-specific splicing changes. It is necessary to accurately identify such splicing changes to elucidate the underlying regulatory programs or link the splicing changes to specific diseases.

Results: We present a probabilistic model tailored for high-throughput AS data, where observed isoform levels are explained as combinations of condition-specific AS signals. According to our formulation, given an AS dataset our tasks are to detect common signals in the data and identify the exons relevant to each signal. Our model can incorporate prior knowledge about underlying AS signals, measurement quality and gene expression level effects. Using a large-scale multi-tissue AS dataset, we demonstrate the advantage of our method over standard alternative approaches. In addition, we describe newly found tissue-specific AS signals which were verified experimentally, and discuss associated regulatory features.

Contact:yoseph@psi.utoronto.ca; frey@psi.utoronto.ca

Supplementary information:Supplementary data are available at Bioinformatics online.

Drosophila ataxin 2-binding protein 1 marks an intermediate step in the molecular differentiation of female germline cysts

In the Drosophila ovary, extrinsic signaling from the niche and intrinsic translational control machinery regulate the balance between germline stem cell maintenance and the differentiation of their daughters. However, the molecules that promote the continued stepwise development of ovarian germ cells after their exit from the niche remain largely unknown. Here, we report that the early development of germline cysts depends on the Drosophila homolog of the human ataxin 2-binding protein 1 (A2BP1) gene. Drosophila A2BP1 protein expression is first observed in the cytoplasm of 4-, 8- and 16-cell cysts, bridging the expression of the early differentiation factor Bam with late markers such as Orb, Rbp9 and Bruno encoded by arrest. The expression of A2BP1 is lost in bam, sans-fille (snf) and mei-P26 mutants, but is still present in other mutants such as rbp9 and arrest. A2BP1 alleles of varying strength produce mutant phenotypes that include germline counting defects and cystic tumors. Phenotypic analysis reveals that strong A2BP1 alleles disrupt the transition from mitosis to meiosis. These mutant cells continue to express high levels of mitotic cyclins and fail to express markers of terminal differentiation. Biochemical analysis reveals that A2BP1 isoforms bind to each other and associate with Bruno, a known translational repressor protein. These data show that A2BP1 promotes the molecular differentiation of ovarian germline cysts.

A network of conserved co-occurring motifs for the regulation of alternative splicing

Cis-acting short sequence motifs play important roles in alternative splicing. It is now possible to identify such sequence motifs as conserved sequence patterns in genome sequence alignments. Here, we report the systematic search for motifs in the neighboring introns of alternatively spliced exons by using comparative analysis of mammalian genome alignments. We identified 11 conserved sequence motifs that might be involved in the regulation of alternative splicing. These motifs are not only significantly overrepresented near alternatively spliced exons, but they also co-occur with each other, thus, forming a network of cis-elements, likely to be the basis for context-dependent regulation. Based on this finding, we applied the motif co-occurrence to predict alternatively skipped exons. We verified exon skipping in 29 cases out of 118 predictions (25%) by EST and mRNA sequences in the databases. For the predictions not verified by the database sequences, we confirmed exon skipping in 10 additional cases by using both RT–PCR experiments and the publicly available RNA-Seq data. These results indicate that even more alternative splicing events will be found with the progress of large-scale and high-throughput analyses for various tissue samples and developmental stages.

Computational identification of tissue-specific alternative splicing elements in mouse genes from RNA-Seq

Tissue-specific alternative splicing is a key mechanism for generating tissue-specific proteomic diversity in eukaryotes. Splicing regulatory elements (SREs) in pre-mature messenger RNA play a very important role in regulating alternative splicing. In this article, we use mouse RNA-Seq data to determine a positive data set where SREs are over-represented and a reliable negative data set where the same SREs are most likely under-represented for a specific tissue and then employ a powerful discriminative approach to identify SREs. We identified 456 putative splicing enhancers or silencers, of which 221 were predicted to be tissue-specific. Most of our tissue-specific SREs are likely different from constitutive SREs, since only 18% of our exonic splicing enhancers (ESEs) are contained in constitutive RESCUE-ESEs. A relatively small portion (20%) of our SREs is included in tissue-specific SREs in human identified in two recent studies. In the analysis of position distribution of SREs, we found that a dozen of SREs were biased to a specific region. We also identified two very interesting SREs that can function as an enhancer in one tissue but a silencer in another tissue from the same intronic region. These findings provide insight into the mechanism of tissue-specific alternative splicing and give a set of valuable putative SREs for further experimental investigations.

Exon-level microarray analyses identify alternative splicing programs in breast cancer

Protein isoforms produced by alternative splicing (AS) of many genes have been implicated in several aspects of cancer genesis and progression. These observations motivated a genome-wide assessment of AS in breast cancer. We accomplished this by measuring exon level expression in 31 breast cancer and nonmalignant immortalized cell lines representing luminal, basal, and claudin-low breast cancer subtypes using Affymetrix Human Junction Arrays. We analyzed these data using a computational pipeline specifically designed to detect AS with a low false-positive rate. This identified 181 splice events representing 156 genes as candidates for AS. Reverse transcription-PCR validation of a subset of predicted AS events confirmed 90%. Approximately half of the AS events were associated with basal, luminal, or claudin-low breast cancer subtypes. Exons involved in claudin-low subtype-specific AS were significantly associated with the presence of evolutionarily conserved binding motifs for the tissue-specific Fox2 splicing factor. Small interfering RNA knockdown of Fox2 confirmed the involvement of this splicing factor in subtype-specific AS. The subtype-specific AS detected in this study likely reflects the splicing pattern in the breast cancer progenitor cells in which the tumor arose and suggests the utility of assays for Fox-mediated AS in cancer subtype definition and early detection. These data also suggest the possibility of reducing the toxicity of protein-targeted breast cancer treatments by targeting protein isoforms that are not present in limiting normal tissues.

Deciphering the splicing code

Alternative splicing has a crucial role in the generation of biological complexity, and its misregulation is often involved in human disease. Here we describe the assembly of a 'splicing code', which uses combinations of hundreds of RNA features to predict tissue-dependent changes in alternative splicing for thousands of exons. The code determines new classes of splicing patterns, identifies distinct regulatory programs in different tissues, and identifies mutation-verified regulatory sequences. Widespread regulatory strategies are revealed, including the use of unexpectedly large combinations of features, the establishment of low exon inclusion levels that are overcome by features in specific tissues, the appearance of features deeper into introns than previously appreciated, and the modulation of splice variant levels by transcript structure characteristics. The code detected a class of exons whose inclusion silences expression in adult tissues by activating nonsense-mediated messenger RNA decay, but whose exclusion promotes expression during embryogenesis. The code facilitates the discovery and detailed characterization of regulated alternative splicing events on a genome-wide scale.

Fox-2 protein regulates the alternative splicing of scleroderma-associated lysyl hydroxylase 2 messenger RNA

OBJECTIVE

Scleroderma (systemic sclerosis [SSc]) is a complex connective tissue disorder characterized by hardening and thickening of the skin. One hallmark of scleroderma is excessive accumulation of collagen accompanied by increased levels of pyridinoline collagen crosslinks derived from hydroxylysine residues in the collagen telopeptide domains. Lysyl hydroxylase 2 (LH2), an important alternatively spliced enzyme in collagen biosynthesis, acts as a collagen telopeptide hydroxylase. Changes in the pattern of LH2 alternative splicing, favoring increased inclusion of the alternatively spliced LH2 exon 13A, thereby increasing the levels of the long transcript of LH2 (LH2[long]), are linked to scleroderma disease. This study was undertaken to examine the role played by RNA binding protein Fox-2 in regulating exon 13A inclusion, which leads to the generation of scleroderma-associated LH2(long) messenger RNA (mRNA).

METHODS

Phylogenetic sequence analysis of introns flanking exon 13A was performed. A tetracycline-inducible system in T-Rex 293 cells was used to induce Fox-2 protein, and endogenous LH2(long) mRNA was determined by reverse transcriptase-polymerase chain reaction. An LH2 minigene was designed, validated, and used in Fox-2 overexpression and mutagenesis experiments. Knockdown of Fox-2 was performed in mouse embryonic fibroblasts and in fibroblasts from SSc patients.

RESULTS

Overexpression of Fox-2 enhanced the inclusion of exon 13A and increased the generation of LH2(long) mRNA, whereas knockdown of Fox-2 decreased LH2(long) transcripts. Mutational analysis of an LH2 minigene demonstrated that 2 of the 4 Fox binding motifs flanking LH2 exon 13A are required for inclusion of exon 13A. In early passage fibroblasts derived from patients with scleroderma, the knockdown of Fox-2 protein significantly decreased the endogenous levels of LH2(long) mRNA.

CONCLUSION

Our findings indicate that Fox-2 plays an integral role in the regulation of LH2 splicing. Knockdown of Fox-2 and other methods to decrease the levels of fibrosis-associated LH2(long) mRNA in primary scleroderma cells may suggest a novel approach to strategies directed against scleroderma.

A method of predicting changes in human gene splicing induced by genetic variants in context of cis-acting elements

Background

Polymorphic variants and mutations disrupting canonical splicing isoforms are among the leading causes of human hereditary disorders. While there is a substantial evidence of aberrant splicing causing Mendelian diseases, the implication of such events in multi-genic disorders is yet to be well understood. We have developed a new tool (SpliceScan II) for predicting the effects of genetic variants on splicing and cis-regulatory elements. The novel Bayesian non-canonical 5'GC splice site (SS) sensor used in our tool allows inference on non-canonical exons.

Results

Our tool performed favorably when compared with the existing methods in the context of genes linked to the Autism Spectrum Disorder (ASD). SpliceScan II was able to predict more aberrant splicing isoforms triggered by the mutations, as documented in DBASS5 and DBASS3 aberrant splicing databases, than other existing methods. Detrimental effects behind some of the polymorphic variations previously associated with Alzheimer's and breast cancer could be explained by changes in predicted splicing patterns.

Conclusions

We have developed SpliceScan II, an effective and sensitive tool for predicting the detrimental effects of genomic variants on splicing leading to Mendelian and complex hereditary disorders. The method could potentially be used to screen resequenced patient DNA to identify de novo mutations and polymorphic variants that could contribute to a genetic disorder.

Splicing, cis genetic variation and disease

Splicing is a post-transcriptional modification of RNA during which introns are removed and exons are joined. Most of the mammalian genes undergo constitutive and alternative splicing events. In addition to the strong signals of the splice sites, splicing is influenced at a distance by a range of trans factors that interact with cis regulatory elements and influence the spliceosome. The intention of the present mini-review is to give some insights into the complexity of this interaction and to introduce the consequences of some kinds of detrimental genetic variation on alternative splicing and disease.

Fox-1 family of RNA-binding proteins

The Fox-1 family of RNA-binding proteins are evolutionarily conserved regulators of tissue-specific alternative splicing in metazoans. The Fox-1 family specifically recognizes the (U)GCAUG stretch in regulated exons or in flanking introns, and either promotes or represses target exons. Recent unbiased bioinformatics analyses of alternatively spliced exons and comparison of various vertebrate genomes identified the (U)GCAUG stretch as a highly conserved and widely distributed element enriched in intronic regions surrounding exons with altered inclusion in muscle, heart, and brain, consistent with specific expression of Fox-1 and Fox-2 in these tissues. Global identification of Fox-2 target RNAs in living cells revealed that many of the Fox-2 target genes themselves encode splicing regulators. Further systematic elucidation of target genes of the Fox-1 family and other splicing regulators in various tissues will lead to a comprehensive understanding of splicing regulatory networks.

Computational prediction of splicing regulatory elements shared by Tetrapoda organisms

Background

Auxiliary splicing sequences play an important role in ensuring accurate and efficient splicing by promoting or repressing recognition of authentic splice sites. These cis-acting motifs have been termed splicing enhancers and silencers and are located both in introns and exons. They co-evolved into an intricate splicing code together with additional functional constraints, such as tissue-specific and alternative splicing patterns. We used orthologous exons extracted from the University of California Santa Cruz multiple genome alignments of human and 22 Tetrapoda organisms to predict candidate enhancers and silencers that have reproducible and statistically significant bias towards annotated exonic boundaries.

Results

A total of 2,546 Tetrapoda enhancers and silencers were clustered into 15 putative core motifs based on their Markov properties. Most of these elements have been identified previously, but 118 putative silencers and 260 enhancers (~15%) were novel. Examination of previously published experimental data for the presence of predicted elements showed that their mutations in 21/23 (91.3%) cases altered the splicing pattern as expected. Predicted intronic motifs flanking 3' and 5' splice sites had higher evolutionary conservation than other sequences within intronic flanks and the intronic enhancers were markedly differed between 3' and 5' intronic flanks.

Conclusion

Difference in intronic enhancers supporting 5' and 3' splice sites suggests an independent splicing commitment for neighboring exons. Increased evolutionary conservation for ISEs/ISSs within intronic flanks and effect of modulation of predicted elements on splicing suggest functional significance of found elements in splicing regulation. Most of the elements identified were shown to have direct implications in human splicing and therefore could be useful for building computational splicing models in biomedical research.

An inducible change in Fox-1/A2BP1 splicing modulates the alternative splicing of downstream neuronal target exons

Neuronal depolarization and CaM kinase IV signaling alter the splicing of multiple exons in transcripts for ion channels, neurotransmitter receptors, and other synaptic proteins. These splicing changes are mediated in part by special CaM kinase-responsive RNA elements, within or adjacent to exons that are repressed in the initial phase of chronic depolarization. The splicing of many neuronal transcripts is also regulated by members of the Fox (Feminizing gene on X) protein family, and these Fox targets are also often proteins affecting synaptic activity. We show that Fox-1/Ataxin 2-Binding Protein 1 (A2BP1), a protein implicated in a variety of neurological diseases, can counteract the effects of chronic depolarization on splicing. We find that exon 19 of Fox-1 is itself repressed by depolarization. Fox-1 transcripts missing exon 19 encode a nuclear isoform of Fox-1 that progressively replaces the cytoplasmic Fox-1 isoform as cells are maintained depolarizing media. The resulting increase in nuclear Fox-1 leads to the reactivation of many Fox-1 target exons, including exon 5 of the NMDA receptor 1, that were initially repressed by the high-KCl medium. These results reveal a novel mechanism for the slow modulation of splicing as cells adapt to chronic stimuli: The subcellular localization of a splicing regulator is controlled through its own alternative splicing.

Genome-wide prediction of cis-acting RNA elements regulating tissue-specific pre-mRNA alternative splicing

Background

Human genes undergo various patterns of pre-mRNA splicing across different tissues. Such variation is primarily regulated by trans-acting factors that bind on exonic and intronic cis-acting RNA elements (CAEs). Here we report a computational method to mechanistically identify cis-acting RNA elements that contribute to the tissue-specific alternative splicing pattern. This method is an extension of our previous model, SplicingModeler, which predicts the significant CAEs that contribute to the splicing differences between two tissues. In this study, we introduce tissue-specific functional levels estimation step, which allows evaluating regulatory functions of predicted CAEs that are involved in more than two tissues.

Results

Using a publicly available Affymetrix Genechip^® Human Exon Array dataset, our method identifies 652 cis-acting RNA elements (CAEs) across 11 human tissues. About one third of predicted CAEs can be mapped to the known RBP (RNA binding protein) binding sites or match with other predicted exonic splicing regulator databases. Interestingly, the vast majority of predicted CAEs are in intronic regulatory regions. A noticeable exception is that many exonic elements are found to regulate the alternative splicing between cerebellum and testes. Most identified elements are found to contribute to the alternative splicing between two tissues, while some are important in multiple tissues. This suggests that genome-wide alternative splicing patterns are regulated by a combination of tissue-specific cis-acting elements and "general elements" whose functional activities are important but differ across multiple tissues.

Conclusion

In this study, we present a model-based computational approach to identify potential cis-acting RNA elements by considering the exon splicing variation as the combinatorial effects of multiple cis-acting regulators. This methodology provides a novel evaluation on the functional levels of cis-acting RNA elements by estimating their tissue-specific functions on various tissues.

Cancer-associated regulation of alternative splicing

Alternative splicing of pre-mRNA increases the diversity of protein functions. Here we show that about half of all active alternative splicing events in ovarian and breast tissues are changed in tumors, and many seem to be regulated by a single factor; sequence analysis revealed binding sites for the RNA binding protein FOX2 downstream of one-third of the exons skipped in cancer. High-resolution analysis of FOX2 binding sites defined the precise positions relative to alternative exons at which the protein may function as either a silencer or an enhancer. Most of the identified targets were shifted in the same direction by FOX2 depletion in cell lines as they were in breast and ovarian cancer tissues. Notably, we found expression of FOX2 itself is downregulated in ovarian cancer and its splicing is altered in breast cancer samples. These results suggest that the decreased expression of FOX2 in cancer tissues modulates splicing and controls proliferation.

A computational approach for genome-wide mapping of splicing factor binding sites

A computational method is presented for genome-wide mapping of splicing factor binding sites that considers both the genomic environment and evolutionary conservation.

Alternative splicing is regulated by splicing factors that serve as positive or negative effectors, interacting with regulatory elements along exons and introns. Here we present a novel computational method for genome-wide mapping of splicing factor binding sites that considers both the genomic environment and the evolutionary conservation of the regulatory elements. The method was applied to study the regulation of different alternative splicing events, uncovering an interesting network of interactions among splicing factors.

Prediction of alternatively skipped exons and splicing enhancers from exon junction arrays

Background

Alternative splicing of exons in a pre-mRNA transcript is an important mechanism which contributes to protein diversity in human. Arrays for detecting alternative splicing are available using several different probe designs, including those based on exon-junctions. In this work, we introduce a new method for predicting alternatively skipped exons from exon-junction arrays. Predictions based on our method are compared against controls and their sequences are analyzed to identify motifs important for regulating alternative splicing.

Results

Our comparison of several alternative methods shows that an exon-skipping score based on neighboring junctions best discriminates between positive and negative controls. Sequence analysis of our predicted exons confirms the presence of known splicing regulatory sequences. In addition, we also derive a set of development-related alternatively spliced genes based on fetal versus adult tissue comparisons and find that our predictions are consistent with their functional annotations. Ab initio motif finding algorithms are applied to identify several motifs that may be relevant for splicing during development.

Conclusion

This work describes a new method for analyzing exon-junction arrays, identifies sequence motifs that are specific for alternative and constitutive splicing and suggests a role for several known splicing factors and their motifs in developmental regulation.

Expression of 24,426 human alternative splicing events and predicted cis regulation in 48 tissues and cell lines

Alternative pre–messenger RNA splicing impacts development, physiology, and disease, but its regulation in humans is not well understood, partially due to the limited scale to which the expression of specific splicing events has been measured. We generated the first genome-scale expression compendium of human alternative splicing events using custom whole-transcript microarrays monitoring expression of 24,426 alternative splicing events in 48 diverse human samples. Over 11,700 genes and 9,500 splicing events were differentially expressed, providing a rich resource for studying splicing regulation. An unbiased, systematic screen of 21,760 4-mer to 7-mer words for cis-regulatory motifs identified 143 RNA 'words' enriched near regulated cassette exons, including six clusters of motifs represented by UCUCU, UGCAUG, UGCU, UGUGU, UUUU, and AGGG, which map to trans-acting regulators PTB, Fox, Muscleblind, CELF/CUG-BP, TIA-1, and hnRNP F/H, respectively. Each cluster showed a distinct pattern of genomic location and tissue specificity. For example, UCUCU occurs 110 to 35 nucleotides preceding cassette exons upregulated in brain and striated muscle but depleted in other tissues. UCUCU and UGCAUG appear to have similar function but independent action, occurring 5' and 3', respectively, of 33% of the cassette exons upregulated in skeletal muscle but co-occurring for only 2%.