A comprehensive computational characterization of conserved mammalian intronic sequences reveals conserved motifs associated with constitutive and alternative splicing

Exon-intron boundary inhibits m6A deposition, enabling m6A distribution hallmark, longer mRNA half-life and flexible protein coding

Regional bias of N6-methyladenosine (m6A) mRNA modification avoiding splice site region, calls for an open hypothesis whether exon-intron boundary could affect m6A deposition. By deep learning modeling, we find that exon-intron boundary represses a proportion (12% to 34%) of m6A deposition at adjacent exons (~100 nt to splice site). Experiments validate that m6A signal increases once the host gene does not undergo pre-mRNA splicing to produce the same mRNA. Inhibited m6A sites have higher m6A enhancers and lower m6A silencers locally and show high heterogeneity at different exons genome-widely, with only a small proportion (12% to 15%) of exons showing strong inhibition, enabling more stable mRNAs and flexible protein coding. m6A is majorly responsible for why mRNAs with more exons be more stable. Exon junction complex (EJC) only partially contributes to this exon-intron boundary m6A inhibition in some short internal exons, highlighting additional factors yet to be identified.

Polymorphic differences within HLA-C alleles contribute to alternatively spliced transcripts lacking exon 5

The HLA genes are amongst the most polymorphic in the human genome. Alternative splicing could add an extra layer of complexity, but has not been studied extensively. Here, we applied an RNA based approach to study the influence of allele polymorphism on alternative splicing of HLA‐C in peripheral blood. RNA was isolated from these peripheral cells, converted into cDNA and amplified specifically for 12 common HLA‐C allele groups. Through subsequent sequencing of HLA‐C, we observed alternative splicing variants of HLA‐C*04 and *16 that resulted in exon 5 skipping and were co‐expressed with the mature transcript. Investigation of intron 4 sequences of HLA‐C*04 and *16 compared with other HLA‐C alleles demonstrated no effect on predicted splice sites and branch point. To further investigate if the unique polymorphic positions in exon 5 of HLA‐C*04 or *16 may facilitate alternative splicing by acting on splicing regulatory elements (SRE), in‐silico splicing analysis was performed. While the HLA‐C*04 specific SNP in exon 5 had no effect on predicted exonic SRE, the HLA‐C*16 specific exon 5 SNP did alter exonic SRE. Our findings provide experimental and theoretical support for the concept that polymorphisms within the HLA‐C alleles influence the alternative splicing of HLA‐C.

GC-content biases in protein-coding genes act as an "mRNA identity" feature for nuclear export

It has long been observed that human protein-coding genes have a particular distribution of GC-content: the 5' end of these genes has high GC-content while the 3' end has low GC-content. In 2012, it was proposed that this pattern of GC-content could act as an mRNA identity feature that would lead to it being better recognized by the cellular machinery to promote its nuclear export. In contrast, junk RNA, which largely lacks this feature, would be retained in the nucleus and targeted for decay. Now two recent papers have provided evidence that GC-content does promote the nuclear export of many mRNAs in human cells.

At Least Two Genes Encode Many Variants of Irak3 in Rainbow Trout, but Neither the Full-Length Factor Nor Its Variants Interfere Directly With the TLR-Mediated Stimulation of Inflammation

The interleukin-1-receptor-associated kinase 3 (IRAK3) is known in mammals as a negative feedback regulator of NF-κB-mediated innate-immune mechanisms. Our RNA-seq experiments revealed a prototypic 1920-nt sequence encoding irak3 from rainbow trout (Oncorhynchus mykiss), as well as 20 variants that vary in length and nucleotide composition. Based on the DNA-sequence information from two closely related irak3 genes from rainbow trout and an irak3-sequence fragment from Atlantic salmon retrieved from public databases, we elucidated the underlying genetic causes for this striking irak3 diversity. Infecting rainbow trout with a lethal dose of Aeromonas salmonicida enhanced the expression of all variants in the liver, head kidney, and peripheral blood leucocytes. We analyzed the functional impact of the full-length factor and selected structural variants by overexpressing them in mammalian HEK-293 cells. The full-length factor enhanced the basal activity of NF-κB, but did not dampen the TLR2-signaling-induced levels of NF-κB activation. Increasing the basal NF-κB-activity through Irak3 apparently does not involve its C-terminal domain. However, more severely truncated factors had only a minor impact on the activity of NF-κB. The TLR2-mediated stimulation did not alter the spatial distribution of Irak3 inside the cells. In salmonid CHSE-214 cells, we observed that the Irak3-splice variant that prominently expresses the C-terminal domain significantly quenched the stimulation-dependent production of interleukin-1β and interleukin-8, but not the production of other immune regulators. We conclude that the different gene and splice variants of Irak3 from trout play distinct roles in the activation of immune-regulatory mechanisms.

Activation of cryptic splice sites in three patients with chronic granulomatous disease

Background

Chronic granulomatous disease (CGD) is a primary immune deficiency caused by mutations in the genes encoding the structural components of the phagocyte NADPH oxidase. As a result, the patients cannot generate sufficient amounts of reactive oxygen species required for killing pathogenic microorganisms.

Methods

We analyzed NADPH oxidase activity and component expression in neutrophils, performed genomic DNA and cDNA analysis, and used mRNA splicing prediction tools to evaluate the impact of mutations.

Results

In two patients with CGD, we had previously found mutations that cause aberrant pre‐mRNA splicing. In one patient an exonic mutation in a cryptic donor splice site caused the deletion of the 3' part of exon 6 from the mRNA of CYBB. This patient suffers from X‐linked CGD. The second patient, with autosomal CGD, has a mutation in the donor splice site of intron 1 of CYBA that activates a cryptic donor splice site downstream in intron 1, causing the insertion of intronic sequences in the mRNA. The third patient, recently analyzed, also with autosomal CGD, has a mutation in intron 4 of CYBA, 15 bp from the acceptor splice site. This mutation weakens a branch site and activates a cryptic acceptor splice site, causing the insertion of 14 intronic nucleotides into the mRNA.

Conclusion

We found three different mutations, one exonic, one in a donor splice site and one intronic, that all caused missplicing of pre‐mRNA. We analyzed these mutations with four different splice prediction programs and found that predictions of splice site strength, splice enhancer and splice silencer protein binding and branch site strength are all essential for correct prediction of pre‐mRNA splicing.

A Multiplexed Assay for Exon Recognition Reveals that an Unappreciated Fraction of Rare Genetic Variants Cause Large-Effect Splicing Disruptions

Mutations that lead to splicing defects can have severe consequences on gene function and cause disease. Here, we explore how human genetic variation affects exon recognition by developing a multiplexed functional assay of splicing using Sort-seq (MFASS). We assayed 27,733 variants in the Exome Aggregation Consortium (ExAC) within or adjacent to 2,198 human exons in the MFASS minigene reporter and found that 3.8% (1,050) of variants, most of which are extremely rare, led to large-effect splice-disrupting variants (SDVs). Importantly, we find that 83% of SDVs are located outside of canonical splice sites, are distributed evenly across distinct exonic and intronic regions, and are difficult to predict a priori. Our results indicate extant, rare genetic variants can have large functional effects on splicing at appreciable rates, even outside the context of disease, and MFASS enables their empirical assessment at scale.

Identification of intronic-splice site mutations in GATA4 gene in Indian patients with congenital heart disease

Congenital Heart Disease (CHD) is the most common birth defect among congenital anomalies that arise before birth. GATA4 transcription factor plays an important role in foetal heart development. Mutational analysis of GATA4 gene in CHD patients revealed five known heterozygous mutations (p.T355S, p.S377G, p.V380M, p.P394T and p.D425N) identified in exons 5 and 6 regions and fifteen intronic variants in the non-coding regions (g.76885T>C/Y,g.76937G>S, g.78343G>R, g.83073T>Y, g.83271C>A/M, g.83318G>K, g.83415G>R, g.83502A>C/M, g.84991G>R, g.85294C>Y, g.85342C>T/Y, g.86268A>R, g.87409G>A/R, g.87725T>Y, g.87813A>T/W). In silico analysis of these intronic variants identified two potential branch point mutations (g.83271C>A/M, g.86268A>R) and predicted effects of these on intronic splice sites as enhancer and silencer motifs. This study attempts to correlate the pattern of intronic variants of GATA4 gene which might provide new insights to unravel the possible molecular etiology of CHD.

Mutation in an exonic splicing enhancer site causing chronic granulomatous disease

In a male patient suffering from X-linked chronic granulomatous disease (CGD) we found a c.389G>T mutation in exon 5 of the CYBB gene. We have analyzed why 95% of the transcripts of this gene lacked exon 5, leading to a frameshift and premature termination codon. The mutation was located in a region comprising three putative exonic splicing enhancer binding sites, for SRSF1, SRFS2 and SRFS6, according to the ESEfinder Tool (http://rulai.cshl.edu/cgi-bin/tools/ESE3/esefinder.cgi). With the Analyser Splice Tool we calculated the probability of skipping of exon 5 in CYBB mRNA, and by means of Sroogle the number of putative binding motifs for splicing enhancer and splicing silencer proteins (http://astlab.tau.ac.il/index.php). These analyses clarify why this exon was skipped in the majority of the mRNA. The normally spliced transcript contains an amino acid change p.Arg130Leu. This poorly expressed transcript gives rise to a protein with low expression but presumably normal activity, leading to a respiratory burst activity in the patient's neutrophils of about 15% of normal.

Mutations of Pre-mRNA Splicing Regulatory Elements: Are Predictions Moving Forward to Clinical Diagnostics?

For more than three decades, researchers have known that consensus splice sites alone are not sufficient regulatory elements to provide complex splicing regulation. Other regulators, so-called splicing regulatory elements (SREs) are needed. Most importantly, their sequence variants often underlie the development of various human disorders. However, due to their variable location and high degeneracy, these regulatory sequences are also very difficult to recognize and predict. Many different approaches aiming to identify SREs have been tried, often leading to the development of in silico prediction tools. While these tools were initially expected to be helpful to identify splicing-affecting mutations in genetic diagnostics, we are still quite far from meeting this goal. In fact, most of these tools are not able to accurately discern the SRE-affecting pathological variants from those not affecting splicing. Nonetheless, several recent evaluations have given appealing results (namely for EX-SKIP, ESRseq and Hexplorer predictors). In this review, we aim to summarize the history of the different approaches to SRE prediction, and provide additional validation of these tools based on patients' clinical data. Finally, we evaluate their usefulness for diagnostic settings and discuss the challenges that have yet to be met.

Does conservation account for splicing patterns?

BACKGROUND

Alternative mRNA splicing is critical to proteomic diversity and tissue and species differentiation. Exclusion of cassette exons, also called exon skipping, is the most common type of alternative splicing in mammals.

RESULTS

We present a computational model that predicts absolute (though not tissue-differential) percent-spliced-in of cassette exons more accurately than previous models, despite not using any 'hand-crafted' biological features such as motif counts. We achieve nearly identical performance using only the conservation score (mammalian phastCons) of each splice junction normalized by average conservation over 100 bp of the corresponding flanking intron, demonstrating that conservation is an unexpectedly powerful indicator of alternative splicing patterns. Using this method, we provide evidence that intronic splicing regulation occurs predominantly within 100 bp of the alternative splice sites and that conserved elements in this region are, as expected, functioning as splicing regulators. We show that among conserved cassette exons, increased conservation of flanking introns is associated with reduced inclusion. We also propose a new definition of intronic splicing regulatory elements (ISREs) that is independent of conservation, and show that most ISREs do not match known binding sites or splicing factors despite being predictive of percent-spliced-in.

CONCLUSIONS

These findings suggest that one mechanism for the evolutionary transition from constitutive to alternative splicing is the emergence of cis-acting splicing inhibitors. The association of our ISREs with differences in splicing suggests the existence of novel RNA-binding proteins and/or novel splicing roles for known RNA-binding proteins.

Cytokine-related and sodium channel polymorphism as candidate predisposing factors for childhood encephalopathy FIRES/AERRPS

Febrile infection-related epilepsy syndrome (FIRES), or acute encephalitis with refractory, repetitive partial seizures (AERRPS), is an epileptic encephalopathy beginning with fever-mediated seizures. The etiology remains unclear. To elucidate the genetic background of FIRES/AERRPS (hereafter FIRES), we recruited 19 Japanese patients, genotyped polymorphisms of the IL1B, IL6, IL10, TNFA, IL1RN, SCN1A and SCN2A genes, and compared their frequency between the patients and controls. For IL1RN, the frequency of a variable number of tandem repeat (VNTR) allele, RN2, was significantly higher in the patients than in controls (p=0.0067), and A allele at rs4251981 in 5' upstream of IL1RN with borderline significance (p=0.015). Haplotype containing RN2 was associated with an increased risk of FIRES (OR 3.88, 95%CI 1.40-10.8, p=0.0057). For SCN1A, no polymorphisms showed a significant association, whereas a missense mutation, R1575C, was found in two patients. For SCN2A, the minor allele frequency of G allele at rs1864885 was higher in patients with borderline significance (p=0.011). We demonstrated the association of IL1RN haplotype containing RN2 with FIRES, and showed a possible association of IL1RN rs4251981 G>A and SCN2A rs1864885 A>G, in Japanese patients. These preliminary findings suggest the involvement of multiple genetic factors in FIRES, which needs to be confirmed by future studies in a larger number of FIRES cases.

Systematical identification of splicing regulatory cis-elements and cognate trans-factors

The majority of human genes undergo alternative splicing to generate multiple isoforms with distinct functions. This process is generally controlled by cis-acting splicing regulatory elements (SREs) that recruit trans-acting factors to promote or inhibit the use of nearby splice sites. The growing interest in understanding the regulatory rules of splicing necessitates the systematic identification of these SREs and their cognate protein factors using experimental and computational approaches. Here we describe a strategy to identify and analyze both cis-acting SREs and trans-acting splicing factors. This strategy involves a cell-based screen to identify SREs from a random sequences library and a modified RNA affinity purification approach to unbiasedly identify the splicing factors. These methods can be adopted to identify splicing enhancers or silencers in both exons and introns, and can be extended to different cultured cells. The resulting SREs and splicing factors can be further analyzed with a series of computational and experimental approaches. This approach will help us to collect a molecular part-list for splicing regulation, providing a rich data source that enables a better understanding of the "splicing code".

Quaking and PTB control overlapping splicing regulatory networks during muscle cell differentiation

Alternative splicing contributes to muscle development, but a complete set of muscle-splicing factors and their combinatorial interactions are unknown. Previous work identified ACUAA ("STAR" motif) as an enriched intron sequence near muscle-specific alternative exons such as Capzb exon 9. Mass spectrometry of myoblast proteins selected by the Capzb exon 9 intron via RNA affinity chromatography identifies Quaking (QK), a protein known to regulate mRNA function through ACUAA motifs in 3' UTRs. We find that QK promotes inclusion of Capzb exon 9 in opposition to repression by polypyrimidine tract-binding protein (PTB). QK depletion alters inclusion of 406 cassette exons whose adjacent intron sequences are also enriched in ACUAA motifs. During differentiation of myoblasts to myotubes, QK levels increase two- to threefold, suggesting a mechanism for QK-responsive exon regulation. Combined analysis of the PTB- and QK-splicing regulatory networks during myogenesis suggests that 39% of regulated exons are under the control of one or both of these splicing factors. This work provides the first evidence that QK is a global regulator of splicing during muscle development in vertebrates and shows how overlapping splicing regulatory networks contribute to gene expression programs during differentiation.

Function of alternative splicing

Almost all polymerase II transcripts undergo alternative pre-mRNA splicing. Here, we review the functions of alternative splicing events that have been experimentally determined. The overall function of alternative splicing is to increase the diversity of mRNAs expressed from the genome. Alternative splicing changes proteins encoded by mRNAs, which has profound functional effects. Experimental analysis of these protein isoforms showed that alternative splicing regulates binding between proteins, between proteins and nucleic acids as well as between proteins and membranes. Alternative splicing regulates the localization of proteins, their enzymatic properties and their interaction with ligands. In most cases, changes caused by individual splicing isoforms are small. However, cells typically coordinate numerous changes in 'splicing programs', which can have strong effects on cell proliferation, cell survival and properties of the nervous system. Due to its widespread usage and molecular versatility, alternative splicing emerges as a central element in gene regulation that interferes with almost every biological function analyzed.

A biophysical model for identifying splicing regulatory elements and their interactions

Alternative splicing (AS) of precursor mRNA (pre-mRNA) is a crucial step in the expression of most eukaryotic genes. Splicing factors (SFs) play an important role in AS regulation by binding to the cis-regulatory elements on the pre-mRNA. Although many splicing factors (SFs) and their binding sites have been identified, their combinatorial regulatory effects remain to be elucidated. In this paper, we derive a biophysical model for AS regulation that integrates combinatorial signals of cis-acting splicing regulatory elements (SREs) and their interactions. We also develop a systematic framework for model inference. Applying the biophysical model to a human RNA-Seq data set, we demonstrate that our model can explain 49.1%–66.5% variance of the data, which is comparable to the best result achieved by biophysical models for transcription. In total, we identified 119 SRE pairs between different regions of cassette exons that may regulate exon or intron definition in splicing, and 77 SRE pairs from the same region that may arise from a long motif or two different SREs bound by different SFs. Particularly, putative binding sites of polypyrimidine tract-binding protein (PTB), heterogeneous nuclear ribonucleoprotein (hnRNP) F/H and E/K are identified as interacting SRE pairs, and have been shown to be consistent with the interaction models proposed in previous experimental results. These results show that our biophysical model and inference method provide a means of quantitative modeling of splicing regulation and is a useful tool for identifying SREs and their interactions. The software package for model inference is available under an open source license.

An exon splice enhancer primes IGF2:IGF2R binding site structure and function evolution

Placental development and genomic imprinting coevolved with parental conflict over resource distribution to mammalian offspring. The imprinted genes IGF2 and IGF2R code for the growth promoter insulin-like growth factor 2 (IGF2) and its inhibitor, mannose 6-phosphate (M6P)/IGF2 receptor (IGF2R), respectively. M6P/IGF2R of birds and fish do not recognize IGF2. In monotremes, which lack imprinting, IGF2 specifically bound M6P/IGF2R via a hydrophobic CD loop. We show that the DNA coding the CD loop in monotremes functions as an exon splice enhancer (ESE) and that structural evolution of binding site loops (AB, HI, FG) improved therian IGF2 affinity. We propose that ESE evolution led to the fortuitous acquisition of IGF2 binding by M6P/IGF2R that drew IGF2R into parental conflict; subsequent imprinting may then have accelerated affinity maturation.

Intronic splicing enhancers, cognate splicing factors and context-dependent regulation rules

Most human genes produce multiple splicing isoforms with distinct functions. To systematically understand splicing regulation, we conducted an unbiased screen and identified >100 intronic splicing enhancers (ISEs) that were clustered by sequence similarity into six groups. All ISEs functioned in another cell type and heterologous introns, and their distribution and conservation patterns in different pre-mRNA regions are similar to exonic splicing silencers. Consistently all ISEs inhibited use of splice sites from exonic locations. The putative trans-factors of each ISE group were identified and validated. Five distinct ISE motifs were recognized by hnRNP H and F whose C-terminal domains were sufficient to render context-dependent activities of ISEs. The sixth group was controlled by factors that either activate or suppress splicing. This work provided a comprehensive picture of general ISE activities and provided new models of how a single element can function oppositely depending on its locations and binding factors.

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.

Identification and experimental validation of splicing regulatory elements in Drosophila melanogaster reveals functionally conserved splicing enhancers in metazoans

RNA sequence elements involved in the regulation of pre-mRNA splicing have previously been identified in vertebrate genomes by computational methods. Here, we apply such approaches to predict splicing regulatory elements in Drosophila melanogaster and compare them with elements previously found in the human, mouse, and pufferfish genomes. We identified 99 putative exonic splicing enhancers (ESEs) and 231 putative intronic splicing enhancers (ISEs) enriched near weak 5' and 3' splice sites of constitutively spliced introns, distinguishing between those found near short and long introns. We found that a significant proportion (58%) of fly enhancer sequences were previously reported in at least one of the vertebrates. Furthermore, 20% of putative fly ESEs were previously identified as ESEs in human, mouse, and pufferfish; while only two fly ISEs, CTCTCT and TTATAA, were identified as ISEs in all three vertebrate species. Several putative enhancer sequences are similar to characterized binding-site motifs for Drosophila and mammalian splicing regulators. To provide additional evidence for the function of putative ISEs, we separately identified 298 intronic hexamers significantly enriched within sequences phylogenetically conserved among 15 insect species. We found that 73 putative ISEs were among those enriched in conserved regions of the D. melanogaster genome. The functions of nine enhancer sequences were verified in a heterologous splicing reporter, demonstrating that these sequences are sufficient to enhance splicing in vivo. Taken together, these data identify a set of predicted positive-acting splicing regulatory motifs in the Drosophila genome and reveal regulatory sequences that are present in distant metazoan genomes.

Quantitative evaluation of all hexamers as exonic splicing elements

We describe a comprehensive quantitative measure of the splicing impact of a complete set of RNA 6-mer sequences by deep sequencing successfully spliced transcripts. All 4096 6-mers were substituted at five positions within two different internal exons in a 3-exon minigene, and millions of successfully spliced transcripts were sequenced after transfection of human cells. The results allowed the assignment of a relative splicing strength score to each mutant molecule. The effect of 6-mers on splicing often depended on their location; much of this context effect could be ascribed to the creation of different overlapping sequences at each site. Taking these overlaps into account, the splicing effect of each 6-mer could be quantified, and 6-mers could be designated as enhancers (ESEseqs) and silencers (ESSseqs), with an ESRseq score indicating their strength. Some 6-mers exhibited positional bias relative to the two splice sites. The distribution and conservation of these ESRseqs in and around human exons supported their classification. Predicted RNA secondary structure effects were also seen: Effective enhancers, silencers and 3' splice sites tend to be single stranded, and effective 5' splice sites tend to be double stranded. 6-mers that may form positive or negative synergy with another were also identified. Chromatin structure may also influence the splicing enhancement observed, as a good correspondence was found between splicing performance and the predicted nucleosome occupancy scores of 6-mers. This approach may prove of general use in defining nucleic acid regulatory motifs, substitute for functional SELEX in most cases, and provide insights about splicing mechanisms.

The population genetics of chronic kidney disease: insights from the MYH9-APOL1 locus

Many rare kidney disorders exhibit a monogenic, Mendelian pattern of inheritance. Population-based genetic studies have identified many genetic variants associated with an increased risk of developing common kidney diseases. Strongly associated variants have potential clinical uses as predictive markers and may advance our understanding of disease pathogenesis. These principles are elegantly illustrated by a region within chromosome 22q12 that has a strong association with common forms of kidney disease. Researchers had identified DNA sequence variants in this locus that were highly associated with an increased prevalence of common chronic kidney diseases in people of African ancestry. Initial research concentrated on MYH9 as the most likely candidate gene; however, population-based whole-genome analysis enabled two independent research teams to discover more strongly associated mutations in the neighboring APOL1 gene. The powerful evolutionary selection pressure of an infectious pathogen in West Africa favored the spread of APOL1 variants that protect against a lethal form of African sleeping sickness but are highly associated with an increased risk of kidney disease. We describe the data sources, process of discovery, and reasons for initial misidentification of the candidate gene, as well as the lessons that can be learned for future population genetics research.

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.

TIA1 prevents skipping of a critical exon associated with spinal muscular atrophy

Prevention of skipping of exon 7 during pre-mRNA splicing of Survival Motor Neuron 2 (SMN2) holds the promise for cure of spinal muscular atrophy (SMA), a leading genetic cause of infant mortality. Here, we report T-cell-restricted intracellular antigen 1 (TIA1) and TIA1-related (TIAR) proteins as intron-associated positive regulators of SMN2 exon 7 splicing. We show that TIA1/TIAR stimulate exon recognition in an entirely novel context in which intronic U-rich motifs are separated from the 5' splice site by overlapping inhibitory elements. TIA1 and TIAR are modular proteins with three N-terminal RNA recognition motifs (RRMs) and a C-terminal glutamine-rich (Q-rich) domain. Our results reveal that any one RRM in combination with a Q domain is necessary and sufficient for TIA1-associated regulation of SMN2 exon 7 splicing in vivo. We also show that increased expression of TIA1 counteracts the inhibitory effect of polypyrimidine tract binding protein, a ubiquitously expressed factor recently implicated in regulation of SMN exon 7 splicing. Our findings expand the scope of TIA1/TIAR in genome-wide regulation of alternative splicing under normal and pathological conditions.

Context-dependent robustness to 5' splice site polymorphisms in human populations

There has been growing evidence for extensive diversity of alternative splicing in human populations. Genetic variants within the 5' splice site can cause splicing differences among human individuals and constitute an important class of human disease mutations. In this study, we explored whether natural variations of splicing could reveal important signals of 5' splice site recognition. In seven lymphoblastoid cell lines of Asian, European and African ancestry, we identified 1174 single nucleotide polymorphisms (SNPs) within the consensus 5' splice site. We selected 129 SNPs predicted to significantly alter the splice site activity, and quantitatively examined their splicing impact in the seven individuals. Surprisingly, outside of the essential GT dinucleotide position, only ∼14% of the tested SNPs altered splicing. Bioinformatic and minigene analyses identified signals that could modify the impact of 5' splice site polymorphisms, most notably a strong 3' splice site and the presence of intronic motifs downstream of the 5' splice site. Strikingly, we found that the poly-G run, a known intronic splicing enhancer, was the most significantly enriched motif downstream of exons unaffected by 5' splice site SNPs. In TRIM62, the upstream 3' splice site and downstream intronic poly-G runs functioned redundantly to protect an exon from its 5' splice site polymorphism. Collectively, our study reveals widespread context-dependent robustness to 5' splice site polymorphisms in human transcriptomes. Consequently, certain exons are more susceptible to 5' splice site mutations. Additionally, our work demonstrates that genetic diversity of alternative splicing can provide significant insights into the splicing code of mammalian cells.

Fox-3 and PSF interact to activate neural cell-specific alternative splicing

Fox-1 family (Fox) proteins, which consist of Fox-1 (A2BP1), Fox-2 (Rbm9) and Fox-3 (NeuN) in mammals, bind to the RNA element UGCAUG and regulate alternative pre-mRNA splicing. However the mechanisms for Fox-regulated splicing are largely unknown. We analyzed the expression pattern of the three Fox proteins as well as neural cell-specific alternative splicing of a cassette exon N30 of nonmuscle myosin heavy chain (NMHC) II-B in the mouse central nervous system. Histological and biochemical analyses following fluorescence-activated cell sorting demonstrate a positive correlation of N30 inclusion and Fox-3 expression. Further, we identified polypyrimidine tract binding protein-associated splicing factor (PSF) as an interacting protein with Fox-3 by affinity-chromatography. In cultured cells, enhancement of N30 inclusion by Fox-3 depends on the presence of PSF. PSF enhances N30 inclusion in a UGCAUG-dependent manner, although it does not bind directly to this element. Fox-3 is recruited to the UGCAUG element downstream of N30 in the endogenous NMHC II-B transcript in a PSF-dependent manner. This study is the first to identify PSF as a coactivator of Fox proteins and provides evidence that the Fox-3 and PSF interaction is an integral part of the mechanism by which Fox proteins regulate activation of alternative exons via a downstream intronic enhancer.

iCLIP predicts the dual splicing effects of TIA-RNA interactions

The regulation of alternative splicing involves interactions between RNA-binding proteins and pre-mRNA positions close to the splice sites. T-cell intracellular antigen 1 (TIA1) and TIA1-like 1 (TIAL1) locally enhance exon inclusion by recruiting U1 snRNP to 5' splice sites. However, effects of TIA proteins on splicing of distal exons have not yet been explored. We used UV-crosslinking and immunoprecipitation (iCLIP) to find that TIA1 and TIAL1 bind at the same positions on human RNAs. Binding downstream of 5' splice sites was used to predict the effects of TIA proteins in enhancing inclusion of proximal exons and silencing inclusion of distal exons. The predictions were validated in an unbiased manner using splice-junction microarrays, RT-PCR, and minigene constructs, which showed that TIA proteins maintain splicing fidelity and regulate alternative splicing by binding exclusively downstream of 5' splice sites. Surprisingly, TIA binding at 5' splice sites silenced distal cassette and variable-length exons without binding in proximity to the regulated alternative 3' splice sites. Using transcriptome-wide high-resolution mapping of TIA-RNA interactions we evaluated the distal splicing effects of TIA proteins. These data are consistent with a model where TIA proteins shorten the time available for definition of an alternative exon by enhancing recognition of the preceding 5' splice site. Thus, our findings indicate that changes in splicing kinetics could mediate the distal regulation of alternative splicing.

Intronic motif pairs cooperate across exons to promote pre-mRNA splicing

BACKGROUND

A very early step in splice site recognition is exon definition, a process that is as yet poorly understood. Communication between the two ends of an exon is thought to be required for this step. We report genome-wide evidence for exons being defined through the combinatorial activity of motifs located in flanking intronic regions.

RESULTS

Strongly co-occurring motifs were found to specifically reside in four intronic regions surrounding a large number of human exons. These paired motifs occur around constitutive and alternative exons but not pseudo exons. Most co-occurring motifs are limited to intronic regions within 100 nucleotides of the exon. They are preferentially associated with weaker exons. Their pairing is conserved in evolution and they exhibit a lower frequency of single nucleotide polymorphism when paired. Paired motifs display specificity with respect to distance from the exon borders and in constitutive versus alternative splicing. Many resemble binding sites for heterogeneous nuclear ribonucleoproteins. Specific pairs are associated with tissue-specific genes, the higher expression of which coincides with that of the pertinent RNA binding proteins. Tested pairs acted synergistically to enhance exon inclusion, and this enhancement was found to be exon-specific.

CONCLUSIONS

The exon-flanking sequence pairs identified here by genomic analysis promote exon inclusion and may play a role in the exon definition step in pre-mRNA splicing. We propose a model in which multiple concerted interactions are required between exonic sequences and flanking intronic sequences to effect exon definition.

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.

Functional selection and systematic analysis of intronic splicing elements identify active sequence motifs and associated splicing factors

Despite the critical role of pre-mRNA splicing in generating proteomic diversity and regulating gene expression, the sequence composition and function of intronic splicing regulatory elements (ISREs) have not been well elucidated. Here, we employed a high-throughput in vivo Screening PLatform for Intronic Control Elements (SPLICE) to identify 125 unique ISRE sequences from a random nucleotide library in human cells. Bioinformatic analyses reveal consensus motifs that resemble splicing regulatory elements and binding sites for characterized splicing factors and that are enriched in the introns of naturally occurring spliced genes, supporting their biological relevance. In vivo characterization, including an RNAi silencing study, demonstrate that ISRE sequences can exhibit combinatorial regulatory activity and that multiple trans-acting factors are involved in the regulatory effect of a single ISRE. Our work provides an initial examination into the sequence characteristics and function of ISREs, providing an important contribution to the splicing code.

Relationship between gene compactness and base composition in rice and human genome

In human, highly expressed genes contain shorter and fewer introns and these have been attributed to selection for economy in transcription and translation. On the other hand, in plants, it has been shown that highly expressed genes tend to be longer than lowly expressed genes. Here, in this study, we analyzed compositional influence on genome organization in both rice and human. We demonstrated that, in GC rich rice genes, highly expressed genes are less compact than lowly expressed genes. In GC-poor class, there is no difference in gene compactness between highly and lowly expressed genes. However, the scenario is different for human as there is no influence of GC composition on gene compactness due to their expression levels. We also reported that, highly expressed rice GC-rich pre-mRNA tend to form less stable secondary structure than that of lowly expressed genes. However, on removing intronic sequences, highly expressed mRNA form a stable secondary structure as compared to lowly expressed GC-rich genes. We suggest that in GC-rich rice genes long introns are under selection for enhancing transcriptional efficiency by modulating pre-mRNA secondary structural stability. Thus evolutionary mechanisms behind genome organization are different between these two genomes (human and rice).

The pivotal roles of TIA proteins in 5' splice-site selection of alu exons and across evolution

More than 5% of alternatively spliced internal exons in the human genome are derived from Alu elements in a process termed exonization. Alus are comprised of two homologous arms separated by an internal polypyrimidine tract (PPT). In most exonizations, splice sites are selected from within the same arm. We hypothesized that the internal PPT may prevent selection of a splice site further downstream. Here, we demonstrate that this PPT enhanced the selection of an upstream 5′ splice site (5′ss), even in the presence of a stronger 5′ss downstream. Deletion of this PPT shifted selection to the stronger downstream 5′ss. This enhancing effect depended on the strength of the downstream 5′ss, on the efficiency of base-pairing to U1 snRNA, and on the length of the PPT. This effect of the PPT was mediated by the binding of TIA proteins and was dependent on the distance between the PPT and the upstream 5′ss. A wide-scale evolutionary analysis of introns across 22 eukaryotes revealed an enrichment in PPTs within ∼20 nt downstream of the 5′ss. For most metazoans, the strength of the 5′ss inversely correlated with the presence of a downstream PPT, indicative of the functional role of the PPT. Finally, we found that the proteins that mediate this effect, TIA and U1C, and in particular their functional domains, are highly conserved across evolution. Overall, these findings expand our understanding of the role of TIA1/TIAR proteins in enhancing recognition of exons, in general, and Alu exons, in particular.

Human genes are composed of functional regions, termed exons, separated by non-functional regions, termed introns. Intronic sequences may gradually accumulate mutations and subsequently become recognized by the splicing machinery as exons, a process termed exonization. Alu elements are prone to undergo exonization: more than 5% of alternatively spliced internal exons in the human genome originate from Alu elements. A typical Alu element is ∼300 nucleotides long, consisting of two arms separated by a polypyrimdine tract (PPT). Interestingly, in most cases, exonization occurs almost exclusively within either the right arm or the left, not both. Here we found that the PPT between the two arms serves as a binding site for TIA proteins and prevents the exon selection process from expanding into downstream regions. To obtain a wider overview of TIA function, we performed a cross-evolutionary analysis within 22 eukaryotes of this protein and of U1C, a protein known to interact with it, and found that functional regions of both these proteins were highly conserved. These findings highlight the pivotal role of TIA proteins in 5′ splice-site selection of Alu exons and exon recognition in general.

Regulation of vascular endothelial growth factor (VEGF) splicing from pro-angiogenic to anti-angiogenic isoforms: a novel therapeutic strategy for angiogenesis

Vascular endothelial growth factor (VEGF) is produced either as a pro-angiogenic or anti-angiogenic protein depending upon splice site choice in the terminal, eighth exon. Proximal splice site selection (PSS) in exon 8 generates pro-angiogenic isoforms such as VEGF₁₆₅, and distal splice site selection (DSS) results in anti-angiogenic isoforms such as VEGF₁₆₅b. Cellular decisions on splice site selection depend upon the activity of RNA-binding splice factors, such as ASF/SF2, which have previously been shown to regulate VEGF splice site choice. To determine the mechanism by which the pro-angiogenic splice site choice is mediated, we investigated the effect of inhibition of ASF/SF2 phosphorylation by SR protein kinases (SRPK1/2) on splice site choice in epithelial cells and in in vivo angiogenesis models. Epithelial cells treated with insulin-like growth factor-1 (IGF-1) increased PSS and produced more VEGF₁₆₅ and less VEGF₁₆₅b. This down-regulation of DSS and increased PSS was blocked by protein kinase C inhibition and SRPK1/2 inhibition. IGF-1 treatment resulted in nuclear localization of ASF/SF2, which was blocked by SPRK1/2 inhibition. Pull-down assay and RNA immunoprecipitation using VEGF mRNA sequences identified an 11-nucleotide sequence required for ASF/SF2 binding. Injection of an SRPK1/2 inhibitor reduced angiogenesis in a mouse model of retinal neovascularization, suggesting that regulation of alternative splicing could be a potential therapeutic strategy in angiogenic pathologies.

Chromatin organization marks exon-intron structure

An increasing body of evidence indicates that transcription and splicing are coupled, and it is accepted that chromatin organization regulates transcription. Little is known about the cross-talk between chromatin structure and exon-intron architecture. By analysis of genome-wide nucleosome-positioning data sets from humans, flies and worms, we found that exons show increased nucleosome-occupancy levels with respect to introns, a finding that we link to differential GC content and nucleosome-disfavoring elements between exons and introns. Analysis of genome-wide chromatin immunoprecipitation data in humans and mice revealed four specific post-translational histone modifications enriched in exons. Our findings indicate that previously described enrichment of H3K36me3 modifications in exons reflects a more fundamental phenomenon, namely increased nucleosome occupancy along exons. Our results suggest an RNA polymerase II-mediated cross-talk between chromatin structure and exon-intron architecture, implying that exon selection may be modulated by chromatin structure.

HRP-2, the Caenorhabditis elegans homolog of mammalian heterogeneous nuclear ribonucleoproteins Q and R, is an alternative splicing factor that binds to UCUAUC splicing regulatory elements

Alternative splicing is regulated by cis sequences in the pre-mRNA that serve as binding sites for trans-acting alternative splicing factors. In a previous study, we used bioinformatics and molecular biology to identify and confirm that the intronic hexamer sequence UCUAUC is a nematode alternative splicing regulatory element. In this study, we used RNA affinity chromatography to identify trans factors that bind to this sequence. HRP-2, the Caenorhabditis elegans homolog of human heterogeneous nuclear ribonucleoproteins Q and R, binds to UCUAUC in the context of unc-52 intronic regulatory sequences as well as to RNAs containing tandem repeats of this sequence. The three Us in the hexamer are the most important determinants of this binding specificity. We demonstrate, using RNA interference, that HRP-2 regulates the alternative splicing of two genes, unc-52 and lin-10, both of which have cassette exons flanked by an intronic UCUAUC motif. We propose that HRP-2 is a protein responsible for regulating alternative splicing through binding interactions with the UCUAUC sequence.

Evolution of alternative splicing regulation: changes in predicted exonic splicing regulators are not associated with changes in alternative splicing levels in primates

Alternative splicing is tightly regulated in a spatio-temporal and quantitative manner. This regulation is achieved by a complex interplay between spliceosomal (trans) factors that bind to different sequence (cis) elements. cis-elements reside in both introns and exons and may either enhance or silence splicing. Differential combinations of cis-elements allows for a huge diversity of overall splicing signals, together comprising a complex ‘splicing code’. Many cis-elements have been identified, and their effects on exon inclusion levels demonstrated in reporter systems. However, the impact of interspecific differences in these elements on the evolution of alternative splicing levels has not yet been investigated at genomic level. Here we study the effect of interspecific differences in predicted exonic splicing regulators (ESRs) on exon inclusion levels in human and chimpanzee. For this purpose, we compiled and studied comprehensive datasets of predicted ESRs, identified by several computational and experimental approaches, as well as microarray data for changes in alternative splicing levels between human and chimpanzee. Surprisingly, we found no association between changes in predicted ESRs and changes in alternative splicing levels. This observation holds across different ESR exon positions, exon lengths, and 5′ splice site strengths. We suggest that this lack of association is mainly due to the great importance of context for ESR functionality: many ESR-like motifs in primates may have little or no effect on splicing, and thus interspecific changes at short-time scales may primarily occur in these effectively neutral ESRs. These results underscore the difficulties of using current computational ESR prediction algorithms to identify truly functionally important motifs, and provide a cautionary tale for studies of the effect of SNPs on splicing in human disease.

SROOGLE: webserver for integrative, user-friendly visualization of splicing signals

Exons are typically only 140 nt in length and are surrounded by intronic oceans that are thousands of nucleotides long. Four core splicing signals, aided by splicing-regulatory sequences (SRSs), direct the splicing machinery to the exon/intron junctions. Many different algorithms have been developed to identify and score the four splicing signals and thousands of putative SRSs have been identified, both computationally and experimentally. Here we describe SROOGLE, a webserver that makes splicing signal sequence and scoring data available to the biologist in an integrated, visual, easily interpretable, and user-friendly format. SROOGLE's input consists of the sequence of an exon and flanking introns. The graphic browser output displays the four core splicing signals with scores based on nine different algorithms and highlights sequences belonging to 13 different groups of SRSs. The interface also offers the ability to examine the effect of point mutations at any given position, as well a range of additional metrics and statistical measures regarding each potential signal. SROOGLE is available at http://sroogle.tau.ac.il, and may also be downloaded as a desktop version.

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.

Improved identification of conserved cassette exons using Bayesian networks

Background

Alternative splicing is a major contributor to the diversity of eukaryotic transcriptomes and proteomes. Currently, large scale detection of alternative splicing using expressed sequence tags (ESTs) or microarrays does not capture all alternative splicing events. Moreover, for many species genomic data is being produced at a far greater rate than corresponding transcript data, hence in silico methods of predicting alternative splicing have to be improved.

Results

Here, we show that the use of Bayesian networks (BNs) allows accurate prediction of evolutionary conserved exon skipping events. At a stringent false positive rate of 0.5%, our BN achieves an improved true positive rate of 61%, compared to a previously reported 50% on the same dataset using support vector machines (SVMs). Incorporating several novel discriminative features such as intronic splicing regulatory elements leads to the improvement. Features related to mRNA secondary structure increase the prediction performance, corroborating previous findings that secondary structures are important for exon recognition. Random labelling tests rule out overfitting. Cross-validation on another dataset confirms the increased performance. When using the same dataset and the same set of features, the BN matches the performance of an SVM in earlier literature. Remarkably, we could show that about half of the exons which are labelled constitutive but receive a high probability of being alternative by the BN, are in fact alternative exons according to the latest EST data. Finally, we predict exon skipping without using conservation-based features, and achieve a true positive rate of 29% at a false positive rate of 0.5%.

Conclusion

BNs can be used to achieve accurate identification of alternative exons and provide clues about possible dependencies between relevant features. The near-identical performance of the BN and SVM when using the same features shows that good classification depends more on features than on the choice of classifier. Conservation based features continue to be the most informative, and hence distinguishing alternative exons from constitutive ones without using conservation based features remains a challenging problem.

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%.

Posttranscriptional regulation of gene networks by GU-rich elements and CELF proteins

GU-rich elements found in pre-mRNA and mRNA transcripts play diverse roles in the control of gene expression by regulating mRNA stability, translation and pre-mRNA processing. Regulatory GU-rich elements are highly conserved throughout evolution, and play major roles in development in diverse species from worms to mammals. The conservation of the GU-rich element allowed it to be identified as a sequence that was enriched in the 3' UTR of human transcripts that exhibited rapid mRNA decay. This element functions, at least in part, as a molecular target for members of the CELF family of RNA-binding proteins, which recruit other components of the cellular posttranscriptional gene regulatory machinery to the transcript. Depending on the context, binding to GU-rich sequences by CELF proteins direct a variety of posttranscriptional regulatory events, including deadenylation, mRNA decay, translation or pre-mRNA processing. Thus, GU-rich elements and CELF proteins serve multiple functions in gene expression regulation and define an important evolutionarily conserved posttranscriptional regulatory network.

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

The precise regulation of many alternative splicing (AS) events by specific splicing factors is essential to determine tissue types and developmental stages. However, the molecular basis of tissue-specific AS regulation and the properties of splicing regulatory networks (SRNs) are poorly understood. Here we comprehensively predict the targets of the brain- and muscle-specific splicing factor Fox-1 (A2BP1) and its paralog Fox-2 (RBM9) and systematically define the corresponding SRNs genome-wide. Fox-1/2 are conserved from worm to human, and specifically recognize the RNA element UGCAUG. We integrate Fox-1/2-binding specificity with phylogenetic conservation, splicing microarray data, and additional computational and experimental characterization. We predict thousands of Fox-1/2 targets with conserved binding sites, at a false discovery rate (FDR) of approximately 24%, including many validated experimentally, suggesting a surprisingly extensive SRN. The preferred position of the binding sites differs according to AS pattern, and determines either activation or repression of exon recognition by Fox-1/2. Many predicted targets are important for neuromuscular functions, and have been implicated in several genetic diseases. We also identified instances of binding site creation or loss in different vertebrate lineages and human populations, which likely reflect fine-tuning of gene expression regulation during evolution.

RNA binding specificity of Drosophila muscleblind

Members of the muscleblind family of RNA binding proteins found in Drosophila and mammals are key players in both the human disease myotonic dystrophy and the regulation of alternative splicing. Recently, the mammalian muscleblind-like protein, MBNL1, has been shown to have interesting RNA binding properties with both endogenous and disease-related RNA targets. Here we report the characterization of RNA binding properties of the Drosophila muscleblind protein Mbl. Mutagenesis of double-stranded CUG repeats demonstrated that Mbl requires pyrimidine-pyrimidine mismatches for binding and that the identity and location of the C-G and G-C base pairs within the repeats are essential for Mbl binding. Systematic evolution of ligands by exponential enrichment (SELEX) was used to identify RNA sequences that bind Mbl with much higher affinity than CUG repeats. The RNA sequences identified by SELEX are structured and contain a five-nucleotide consensus sequence of 5'-AGUCU-3'. RNase footprinting of one of the SELEX RNA sequences with Mbl showed that Mbl binds both double-stranded and single-stranded regions of the RNA. Three guanosines show the strongest footprint in the presence of Mbl; mutation of any of these three guanosines eliminates Mbl binding. It was also found that Mbl specifically bound a human MBNL1 RNA target, demonstrating the conservation of the muscleblind proteins in recognizing RNA targets. Our results reveal that Mbl recognizes complex RNA secondary structures.

Identification of motifs that function in the splicing of non-canonical introns

The enrichment of specific intronic splicing enhancers upstream of weak PY tracts suggests a novel mechanism for intron recognition that compensates for a weakened canonical pre-mRNA splicing motif.

Background

While the current model of pre-mRNA splicing is based on the recognition of four canonical intronic motifs (5' splice site, branchpoint sequence, polypyrimidine (PY) tract and 3' splice site), it is becoming increasingly clear that splicing is regulated by both canonical and non-canonical splicing signals located in the RNA sequence of introns and exons that act to recruit the spliceosome and associated splicing factors. The diversity of human intronic sequences suggests the existence of novel recognition pathways for non-canonical introns. This study addresses the recognition and splicing of human introns that lack a canonical PY tract. The PY tract is a uridine-rich region at the 3' end of introns that acts as a binding site for U2AF65, a key factor in splicing machinery recruitment.

Results

Human introns were classified computationally into low- and high-scoring PY tracts by scoring the likely U2AF65 binding site strength. Biochemical studies confirmed that low-scoring PY tracts are weak U2AF65 binding sites while high-scoring PY tracts are strong U2AF65 binding sites. A large population of human introns contains weak PY tracts. Computational analysis revealed many families of motifs, including C-rich and G-rich motifs, that are enriched upstream of weak PY tracts. In vivo splicing studies show that C-rich and G-rich motifs function as intronic splicing enhancers in a combinatorial manner to compensate for weak PY tracts.

Conclusion

The enrichment of specific intronic splicing enhancers upstream of weak PY tracts suggests that a novel mechanism for intron recognition exists, which compensates for a weakened canonical pre-mRNA splicing motif.

A systematic analysis of intronic sequences downstream of 5' splice sites reveals a widespread role for U-rich motifs and TIA1/TIAL1 proteins in alternative splicing regulation

To identify human intronic sequences associated with 5' splice site recognition, we performed a systematic search for motifs enriched in introns downstream of both constitutive and alternative cassette exons. Significant enrichment was observed for U-rich motifs within 100 nucleotides downstream of 5' splice sites of both classes of exons, with the highest enrichment between positions +6 and +30. Exons adjacent to U-rich intronic motifs contain lower frequencies of exonic splicing enhancers and higher frequencies of exonic splicing silencers, compared with exons not followed by U-rich intronic motifs. These findings motivated us to explore the possibility of a widespread role for U-rich motifs in promoting exon inclusion. Since cytotoxic granule-associated RNA binding protein (TIA1) and TIA1-like 1 (TIAL1; also known as TIAR) were previously shown in vitro to bind to U-rich motifs downstream of 5' splice sites, and to facilitate 5' splice site recognition in vitro and in vivo, we investigated whether these factors function more generally in the regulation of splicing of exons followed by U-rich intronic motifs. Simultaneous knockdown of TIA1 and TIAL1 resulted in increased skipping of 36/41 (88%) of alternatively spliced exons associated with U-rich motifs, but did not affect 32/33 (97%) alternatively spliced exons that are not associated with U-rich motifs. The increase in exon skipping correlated with the proximity of the first U-rich motif and the overall "U-richness" of the adjacent intronic region. The majority of the alternative splicing events regulated by TIA1/TIAL1 are conserved in mouse, and the corresponding genes are associated with diverse cellular functions. Based on our results, we estimate that approximately 15% of alternative cassette exons are regulated by TIA1/TIAL1 via U-rich intronic elements.

Splicing regulation: from a parts list of regulatory elements to an integrated splicing code

Alternative splicing of pre-mRNAs is a major contributor to both proteomic diversity and control of gene expression levels. Splicing is tightly regulated in different tissues and developmental stages, and its disruption can lead to a wide range of human diseases. An important long-term goal in the splicing field is to determine a set of rules or "code" for splicing that will enable prediction of the splicing pattern of any primary transcript from its sequence. Outside of the core splice site motifs, the bulk of the information required for splicing is thought to be contained in exonic and intronic cis-regulatory elements that function by recruitment of sequence-specific RNA-binding protein factors that either activate or repress the use of adjacent splice sites. Here, we summarize the current state of knowledge of splicing cis-regulatory elements and their context-dependent effects on splicing, emphasizing recent global/genome-wide studies and open questions.

Multifactorial interplay controls the splicing profile of Alu-derived exons

Exonization of Alu elements creates primate-specific genomic diversity. Here we combine bioinformatic and experimental methodologies to reconstruct the molecular changes leading to exon selection. Our analyses revealed an intricate network involved in Alu exonization. A typical Alu element contains multiple sites with the potential to serve as 5' splice sites (5'ss). First, we demonstrated the role of 5'ss strength in controlling exonization events. Second, we found that a cryptic 5'ss enhances the selection of a more upstream site and demonstrate that this is mediated by binding of U1 snRNA to the cryptic splice site, challenging the traditional role attributed to U1 snRNA of binding the 5'ss only. Third, we used a simple algorithm to identify specific sequences that determine splice site selection within specific Alu exons. Finally, by inserting identical exons within different sequences, we demonstrated the importance of flanking genomic sequences in determining whether an Alu exon will undergo exonization. Overall, our results demonstrate the complex interplay between at least four interacting layers that affect Alu exonization. These results shed light on the mechanism through which Alu elements enrich the primate transcriptome and allow a better understanding of the exonization process in general.

Alternative splicing of Alu exons--two arms are better than one

Alus, primate-specific retroelements, are the most abundant repetitive elements in the human genome. They are composed of two related but distinct monomers, left and right arms. Intronic Alu elements may acquire mutations that generate functional splice sites, a process called exonization. Most exonizations occur in right arms of antisense Alu elements, and are alternatively spliced. Here we show that without the left arm, exonization of the right arm shifts from alternative to constitutive splicing. This eliminates the evolutionary conserved isoform and may thus be selected against. We further show that insertion of the left arm downstream of a constitutively spliced non-Alu exon shifts splicing from constitutive to alternative. Although the two arms are highly similar, the left arm is characterized by weaker splicing signals and lower exonic splicing regulatory (ESR) densities. Mutations that improve these potential splice signals activate exonization and shift splicing from the right to the left arm. Collaboration between two or more putative splice signals renders the intronic left arm with a pseudo-exon function. Thus, the dimeric form of the Alu element fortuitously provides it with an evolutionary advantage, allowing enrichment of the primate transcriptome without compromising its original repertoire.