Simultaneous assessment of the effects of exonic mutations on RNA splicing and protein functions

To simultaneously assess the effects of exonic mutations on RNA splicing and protein functions, we report here an intron-inclusive cDNA (Intinc) expression system. As a test model, twenty-four mutations in exon 9 of the phenylalanine hydroxylase (PAH) gene were examined in an Intinc expression plasmid composed of the PAH cDNA with the exon 9 flanked by its authentic introns. When the PAH enzyme activities from the Intinc plasmid-transfected cells were compared to those of a standard cDNA expression system, five mutations resulted in significant relative differences in PAH activities attributed to altered exon 9-inclusive mRNA levels. Two of the mutations affected exon recognition probably through splice site modifications and the remaining three affected experimentally verified exon splicing enhancer (ESE) motifs. The Intinc expression system allows not only a better link between mutation genotype to disease phenotype but also contributes to further understanding of molecular mechanisms of deleterious effects of mutations.     

Missense mutations in cancer suppressor gene TP53 are colocalized with exonic splicing enhancers (ESEs)

Mutation databases can be viewed as footprints of functional organization of a gene and thus can be used to infer its functional organization. We studied the association of exonic splicing enhancers (ESEs) with missense mutations in the tumor suppressor gene TP53 using the International Agency for Research on Cancer (IARC) mutation database. The goals of the study were: (i) to verify the hypothesis that deleterious missense mutations are colocalized with ESEs; (ii) to identify potentially functional ESE sites in the open reading frame (ORF) of the TP53. If some sequence functions as a splicing enhancer, then nucleotide substitutions in the site will disturb splicing, abrogate p53 function, and cause an increased susceptibility to cancer. Therefore, among cancers showing p53 mutations, more missense mutations are expected within functional ESE sites as compared to non-functional ESE motifs. Using several statistical tests, we found that missense mutations in TP53 are strongly colocalized with ESEs, and that only a small fraction of ESE sites contributes to the association. There are usually one or two ESEs per exon showing a statistically significant association with missense mutations--so-called significant ESE sites. In many respects significant ESE sites are different from those that do not show association with missense mutations. We found that positions of significant ESE sites are codon-dependent--significant ESEs preferentially start from the first position of a codon, whereas non-significant ESEs show no position dependence. Significant ESEs showed a more limited set of sequences compared to non-significant ESEs. These findings suggest that there is a limited number of missense mutations that influence ESE sites and our analysis provides further insight into the types of sites that harbor exonic enhancer elements.     

Missense mutations in hMLH1 and hMSH2 are associated with exonic splicing enhancers

There is a critical need to understand why missense mutations are deleterious. The deleterious effects of missense mutations are commonly attributed to their impact on primary amino acid sequence and protein structure. However, several recent studies have shown that some missense mutations are deleterious because they disturb cis-acting splicing elements-so-called "exonic splicing enhancers" (ESEs). It is not clear whether the ESE-related deleterious effects of missense mutations are common. We have evaluated colocalization of pathogenic missense mutations (found in affected individuals) with high-score ESE motifs in the human mismatch-repair genes hMSH2 and hMLH1. We found that pathogenic missense mutations in the hMSH2 and hMLH1 genes are located in ESE sites significantly more frequently than expected. Pathogenic missense mutations also tended to decrease ESE scores, thus leading to a higher propensity for splicing defects. In contrast, nonpathogenic missense mutations (polymorphisms found in unaffected individuals) and nonsense mutations are distributed randomly in relation to ESE sites. Comparison of the observed and expected frequencies of missense mutations in ESE sites shows that pathogenic effects of >/=20% of mutations in hMSH2 result from disruption of ESE sites and disturbed splicing. Similarly, pathogenic effects of >/=16% of missense mutations in the hMLH1 gene are ESE related. The colocalization of pathogenic missense mutations with ESE sites strongly suggests that their pathogenic effects are splicing related.     

Exon identity established through differential antagonism between exonic splicing silencer-bound hnRNP A1 and enhancer-bound SR proteins.

SR proteins recognize exonic splicing enhancer (ESE) elements and promote exon use, whereas certain hnRNP proteins bind to exonic splicing silencer (ESS) elements and block exon recognition. We investigated how ESS3 in HIV-1 tat exon 3 blocks splicing promoted by one SR protein (SC35) but not another (SF2/ASF). hnRNP A1 mediates silencing by binding initially to a required high-affinity site in ESS3, which then promotes further hnRNP A1 association with the upstream region of the exon. Both SC35 and SF2/ASF recognize upstream ESE motifs, but only SF2/ASF prevents secondary hnRNP A1 binding, presumably by blocking its cooperative propagation along the exon. The differential antagonism between a negative and two positive regulators exemplifies how inclusion of an alternative exon can be modulated.     

RNA folding affects the recruitment of SR proteins by mouse and human polypurinic enhancer elements in the fibronectin EDA exon

In humans, inclusion or exclusion of the fibronectin EDA exon is mainly regulated by a polypurinic enhancer element (exonic splicing enhancer [ESE]) and a nearby silencer element (exonic splicing silencer [ESS]). While human and mouse ESEs behave identically, mutations introduced into the homologous mouse ESS sequence result either in no change in splicing efficiency or in complete exclusion of the exon. Here, we show that this apparently contradictory behavior cannot be simply accounted for by a localized sequence variation between the two species. Rather, the nucleotide differences as a whole determine several changes in the respective RNA secondary structures. By comparing how the two different structures respond to homologous deletions in their putative ESS sequences, we show that changes in splicing behavior can be accounted for by a differential ESE display in the two RNAs. This is confirmed by RNA-protein interaction analysis of levels of SR protein binding to each exon. The immunoprecipitation patterns show the presence of complex multi-SR protein-RNA interactions that are lost with secondary-structure variations after the introduction of ESE and ESS variations. Taken together, our results demonstrate that the sequence context, in addition to the primary sequence identity, can heavily contribute to the making of functional units capable of influencing pre-mRNA splicing.     

The Exon Splicing Silencer in Human Immunodeficiency Virus Type 1 Tat Exon 3 Is Bipartite and Acts Early in Spliceosome Assembly

Inefficient splicing of human immunodeficiency virus type 1 (HIV-1) RNA is necessary to preserve unspliced and singly spliced viral RNAs for transport to the cytoplasm by the Rev-dependent pathway. Signals within the HIV-1 genome that control the rate of splicing include weak 3′ splice sites, exon splicing enhancers (ESE), and exon splicing silencers (ESS). We have previously shown that an ESS present within tat exon 2 (ESS2) and a suboptimal 3′ splice site together act to inhibit splicing at the 3′ splice site flanking tat exon 2. This occurs at an early step in spliceosome assembly. Splicing at the 3′ splice site flanking tat exon 3 is regulated by a bipartite element composed of an ESE and an ESS (ESS3). Here we show that ESS3 is composed of two smaller elements (AGAUCC and UUAG) that can inhibit splicing independently. We also show that ESS3 is more active in the context of a heterologous suboptimal splice site than of an optimal 3′ splice site. ESS3 inhibits splicing by blocking the formation of a functional spliceosome at an early step, since A complexes are not detected in the presence of ESS3. Competitor RNAs containing either ESS2 or ESS3 relieve inhibition of splicing of substrates containing ESS3 or ESS2. This suggests that a common cellular factor(s) may be required for the inhibition of tat mRNA splicing mediated by ESS2 and ESS3.     

ESEfinder: A web resource to identify exonic splicing enhancers

Point mutations frequently cause genetic diseases by disrupting the correct pattern of pre-mRNA splicing. The effect of a point mutation within a coding sequence is traditionally attributed to the deduced change in the corresponding amino acid. However, some point mutations can have much more severe effects on the structure of the encoded protein, for example when they inactivate an exonic splicing enhancer (ESE), thereby resulting in exon skipping. ESEs also appear to be especially important in exons that normally undergo alternative splicing. Different classes of ESE consensus motifs have been described, but they are not always easily identified. ESEfinder (http://exon.cshl.edu/ESE/) is a web-based resource that facilitates rapid analysis of exon sequences to identify putative ESEs responsive to the human SR proteins SF2/ASF, SC35, SRp40 and SRp55, and to predict whether exonic mutations disrupt such elements.     

Determinants of exon 7 splicing in the spinal muscular atrophy genes, SMN1 and SMN2

Spinal muscular atrophy is a neurodegenerative disorder caused by the deletion or mutation of the survival-of-motor-neuron gene, SMN1. An SMN1 paralog, SMN2, differs by a C→T transition in exon 7 that causes substantial skipping of this exon, such that SMN2 expresses only low levels of functional protein. A better understanding of SMN splicing mechanisms should facilitate the development of drugs that increase survival motor neuron (SMN) protein levels by improving SMN2 exon 7 inclusion. In addition, exonic mutations that cause defective splicing give rise to many genetic diseases, and the SMN1/2 system is a useful paradigm for understanding exon-identity determinants and alternative-splicing mechanisms. Skipping of SMN2 exon 7 was previously attributed either to the loss of an SF2/ASF–dependent exonic splicing enhancer or to the creation of an hnRNP A/B–dependent exonic splicing silencer, as a result of the C→T transition. We report the extensive testing of the enhancer-loss and silencer-gain models by mutagenesis, RNA interference, overexpression, RNA splicing, and RNA-protein interaction experiments. Our results support the enhancer-loss model but also demonstrate that hnRNP A/B proteins antagonize SF2/ASF–dependent ESE activity and promote exon 7 skipping by a mechanism that is independent of the C→T transition and is, therefore, common to both SMN1 and SMN2. Our findings explain the basis of defective SMN2 splicing, illustrate the fine balance between positive and negative determinants of exon identity and alternative splicing, and underscore the importance of antagonistic splicing factors and exonic elements in a disease context.     

The Splicing Efficiency of Activating HRAS Mutations Can Determine Costello Syndrome Phenotype and Frequency in Cancer

Costello syndrome (CS) may be caused by activating mutations in codon 12/13 of the HRAS proto-oncogene. HRAS p.Gly12Val mutations have the highest transforming activity, are very frequent in cancers, but very rare in CS, where they are reported to cause a severe, early lethal, phenotype. We identified an unusual, new germline p.Gly12Val mutation, c.35_36GC>TG, in a 12-year-old boy with attenuated CS. Analysis of his HRAS cDNA showed high levels of exon 2 skipping. Using wild type and mutant HRAS minigenes, we confirmed that c.35_36GC>TG results in exon 2 skipping by simultaneously disrupting the function of a critical Exonic Splicing Enhancer (ESE) and creation of an Exonic Splicing Silencer (ESS). We show that this vulnerability of HRAS exon 2 is caused by a weak 3’ splice site, which makes exon 2 inclusion dependent on binding of splicing stimulatory proteins, like SRSF2, to the critical ESE. Because the majority of cancer- and CS- causing mutations are located here, they affect splicing differently. Therefore, our results also demonstrate that the phenotype in CS and somatic cancers is not only determined by the different transforming potentials of mutant HRAS proteins, but also by the efficiency of exon 2 inclusion resulting from the different HRAS mutations. Finally, we show that a splice switching oligonucleotide (SSO) that blocks access to the critical ESE causes exon 2 skipping and halts proliferation of cancer cells. This unravels a potential for development of new anti-cancer therapies based on SSO-mediated HRAS exon 2 skipping.     

An extended inhibitory context causes skipping of exon 7 of SMN2 in spinal muscular atrophy

SMN1 and SMN2 represent the two nearly identical copies of the survival of motor neuron gene in humans. The most frequent cause of spinal muscular atrophy (SMA) is loss of SMN1 accompanied by the inability of SMN2 to compensate due to an inhibitory mutation at position 6 in exon 7 (C6U) that causes exon 7 exclusion. How this single exonic nucleotide regulates exon 7 recognition has been of major interest. Based on score matrices and in vitro assays, abrogation of an exonic splicing enhancer (ESE) associated with SF2/ASF has been considered as the cause of exon 7 exclusion. However, a recent report supports the creation of an exonic splicing silencer (ESS) associated with hnRNP A1 as the determining factor for exon 7 exclusion. Here we show that C6U strengthens an inhibitory context that covers a larger sequence than the hnRNP A1 binding site. The inhibitory context can also be strengthened by the addition of a G residue at the first position of exon 7 in SMN1, promoting exon 7 skipping despite the presence of SF2/ASF binding site. Through in vivo selection and a series of mutations we demonstrate that the strengthening of the extended inhibitory context at the 5' end of exon 7 is exercised through overlapping sequence motifs that collaborate to regulate exon usage.     

Functional Analysis of Mutations in Exon 9 of NF1 Reveals the Presence of Several Elements Regulating Splicing

Neurofibromatosis type 1 (NF1) is one of the most common human hereditary disorders, predisposing individuals to the development of benign and malignant tumors in the nervous system, as well as other clinical manifestations. NF1 is caused by heterozygous mutations in the NF1 gene and around 25% of the pathogenic changes affect pre-mRNA splicing. Since the molecular mechanisms affected by these mutations are poorly understood, we have analyzed the splicing mutations identified in exon 9 of NF1, which is particularly prone to such changes, to better define the possible splicing regulatory elements. Using a minigene approach, we studied the effect of five splicing mutations in this exon described in patients. These highlighted three regulatory motifs within the exon. An in vivo splicing analysis of an extensive collection of changes generated in the minigene demonstrated that the CG motif at c.910-911 is critical for the recognition of exon 9. We also found that the GC motif at c.945-946 is involved in exon recognition through SRSF2 and that this motif is part of a Composite Exon Splicing Regulatory Element made up of physically overlapping enhancer and silencer elements. Finally, through an in vivo splicing analysis and in vitro binding assays, we demonstrated that the c.1007G>A mutation creates an Exonic Splicing Silencer element that binds the hnRNPA1 protein. The complexity of the splicing regulatory elements present in exon 9 is most likely responsible for the fact that mutations in this region represent 25% of all exonic changes that affect splicing in the NF1 gene.     

Optimization of a weak 3' splice site counteracts the function of a bovine papillomavirus type 1 exonic splicing suppressor in vitro and in vivo

Alternative splicing is a critical component of the early to late switch in papillomavirus gene expression. In bovine papillomavirus type 1 (BPV-1), a switch in 3' splice site utilization from an early 3' splice site at nucleotide (nt) 3225 to a late-specific 3' splice site at nt 3605 is essential for expression of the major capsid (L1) mRNA. Three viral splicing elements have recently been identified between the two alternative 3' splice sites and have been shown to play an important role in this regulation. A bipartite element lies approximately 30 nt downstream of the nt 3225 3' splice site and consists of an exonic splicing enhancer (ESE), SE1, followed immediately by a pyrimidine-rich exonic splicing suppressor (ESS). A second ESE (SE2) is located approximately 125 nt downstream of the ESS. We have previously demonstrated that the ESS inhibits use of the suboptimal nt 3225 3' splice site in vitro through binding of cellular splicing factors. However, these in vitro studies did not address the role of the ESS in the regulation of alternative splicing. In the present study, we have analyzed the role of the ESS in the alternative splicing of a BPV-1 late pre-mRNA in vivo. Mutation or deletion of just the ESS did not significantly change the normal splicing pattern where the nt 3225 3' splice site is already used predominantly. However, a pre-mRNA containing mutations in SE2 is spliced predominantly using the nt 3605 3' splice site. In this context, mutation of the ESS restored preferential use of the nt 3225 3' splice site, indicating that the ESS also functions as a splicing suppressor in vivo. Moreover, optimization of the suboptimal nt 3225 3' splice site counteracted the in vivo function of the ESS and led to preferential selection of the nt 3225 3' splice site even in pre-mRNAs with SE2 mutations. In vitro splicing assays also showed that the ESS is unable to suppress splicing of a pre-mRNA with an optimized nt 3225 3' splice site. These data confirm that the function of the ESS requires a suboptimal upstream 3' splice site. A surprising finding of our study is the observation that SE1 can stimulate both the first and the second steps of splicing.     

Depletion of TDP 43 overrides the need for exonic and intronic splicing enhancers in the human apoA-II gene

Exon 3 of the human apolipoprotein A-II (apoA-II) gene is efficiently included in the mRNA although its acceptor site is significantly weak because of a peculiar (GU)₁₆ tract instead of a canonical polypyrimidine tract within the intron 2/exon 3 junction. Our previous studies demonstrated that the SR proteins ASF/SF2 and SC35 bind specifically an exonic splicing enhancer (ESE) within exon 3 and promote exon 3 splicing. In the present study, we show that the ESE is necessary only in the proper context. In addition, we have characterized two novel sequences in the flanking introns that modulate apoA-II exon 3 splicing. There is a G-rich element in intron 2 that interacts with hnRNPH1 and inhibits exon 3 splicing. The second is a purine rich region in intron 3 that binds SRp40 and SRp55 and promotes exon 3 inclusion in mRNA. We have also found that the (GU) repeats in the apoA-II context bind the splicing factor TDP-43 and interfere with exon 3 definition. Significantly, blocking of TDP-43 expression by small interfering RNA overrides the need for all the other cis-acting elements making exon 3 inclusion constitutive even in the presence of disrupted exonic and intronic enhancers. Altogether, our results suggest that exonic and intronic enhancers have evolved to balance the negative effects of the two silencers located in intron 2 and hence rescue the constitutive exon 3 inclusion in apoA-II mRNA.     

hnRNP A1 controls HIV-1 mRNA splicing through cooperative binding to intron and exon splicing silencers in the context of a conserved secondary structure

The removal of the second intron in the HIV-1 rev/tat pre-mRNAs, which involves the joining of splice site SD4 to SA7, is inhibited by hnRNP A1 by a mechanism that requires the intronic splicing silencer (ISS) and the exon splicing silencer (ESS3). In this study, we have determined the RNA secondary structure and the hnRNP A1 binding sites within the 3' splice site region by phylogenetic comparison and chemical/enzymatic probing. A biochemical characterization of the RNA/protein complexes demonstrates that hnRNP A1 binds specifically to primarily three sites, the ISS, a novel UAG motif in the exon splicing enhancer (ESE) and the ESS3 element, which are all situated in experimentally supported stem loop structures. A mutational analysis of the ISS region revealed that the core hnRNP A1 binding site directly overlaps with a major branchpoint used in splicing to SA7, thereby providing a direct explanation for the inhibition of U2 snRNP association with the pre-mRNA by hnRNP A1. Binding of hnRNP A1 to the ISS core site is inhibited by RNA structure but strongly stimulated by the exonic silencer, ESS3. Moreover, the ISS also stimulate binding of hnRNP A1 to the exonic splicing regulators ESS3 and the ESE. Our results suggest a model where a network is formed between hnRNP A1 molecules situated at discrete sites in the intron and exon and that these interactions preclude the recognition of essential splicing signals including the branch point.     

Exonic Splicing Mutations Are More Prevalent than Currently Estimated and Can Be Predicted by Using In Silico Tools

The identification of a causal mutation is essential for molecular diagnosis and clinical management of many genetic disorders. However, even if next-generation exome sequencing has greatly improved the detection of nucleotide changes, the biological interpretation of most exonic variants remains challenging. Moreover, particular attention is typically given to protein-coding changes often neglecting the potential impact of exonic variants on RNA splicing. Here, we used the exon 10 of MLH1, a gene implicated in hereditary cancer, as a model system to assess the prevalence of RNA splicing mutations among all single-nucleotide variants identified in a given exon. We performed comprehensive minigene assays and analyzed patient’s RNA when available. Our study revealed a staggering number of splicing mutations in MLH1 exon 10 (77% of the 22 analyzed variants), including mutations directly affecting splice sites and, particularly, mutations altering potential splicing regulatory elements (ESRs). We then used this thoroughly characterized dataset, together with experimental data derived from previous studies on BRCA1, BRCA2, CFTR and NF1, to evaluate the predictive power of 3 in silico approaches recently described as promising tools for pinpointing ESR-mutations. Our results indicate that ΔtESRseq and ΔHZ_EI-based approaches not only discriminate which variants affect splicing, but also predict the direction and severity of the induced splicing defects. In contrast, the ΔΨ-based approach did not show a compelling predictive power. Our data indicates that exonic splicing mutations are more prevalent than currently appreciated and that they can now be predicted by using bioinformatics methods. These findings have implications for all genetically-caused diseases.