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.     

Seemingly neutral polymorphic variants may confer immunity to splicing-inactivating mutations: a synonymous SNP in exon 5 of MCAD protects from deleterious mutations in a flanking exonic splicing enhancer

The idea that point mutations in exons may affect splicing is intriguing and adds an additional layer of complexity when evaluating their possible effects. Even in the best-studied examples, the molecular mechanisms are not fully understood. Here, we use patient cells, model minigenes, and in vitro assays to show that a missense mutation in exon 5 of the medium-chain acyl-CoA dehydrogenase (MCAD) gene primarily causes exon skipping by inactivating a crucial exonic splicing enhancer (ESE), thus leading to loss of a functional protein and to MCAD deficiency. This ESE functions by antagonizing a juxtaposed exonic splicing silencer (ESS) and is necessary to define a suboptimal 3′ splice site. Remarkably, a synonymous polymorphic variation in MCAD exon 5 inactivates the ESS, and, although this has no effect on splicing by itself, it makes splicing immune to deleterious mutations in the ESE. Furthermore, the region of MCAD exon 5 that harbors these elements is nearly identical to the exon 7 region of the survival of motor neuron (SMN) genes that contains the deleterious silent mutation in SMN2, indicating a very similar and finely tuned interplay between regulatory elements in these two genes. Our findings illustrate a mechanism for dramatic context-dependent effects of single-nucleotide polymorphisms on gene-expression regulation and show that it is essential that potential deleterious effects of mutations on splicing be evaluated in the context of the relevant haplotype.     

Full-exon pyrosequencing screening of BRCA germline mutations in Mexican women with inherited breast and ovarian cancer

Hereditary breast cancer comprises 10% of all breast cancers. The most prevalent genes causing this pathology are BRCA1 and BRCA2 (breast cancer early onset 1 and 2), which also predispose to other cancers. Despite the outstanding relevance of genetic screening of BRCA deleterious variants in patients with a history of familial cancer, this practice is not common in Latin American public institutions. In this work we assessed mutations in the entire exonic and splice-site regions of BRCA in 39 patients with breast and ovarian cancer and with familial history of breast cancer or with clinical features suggestive for BRCA mutations by massive parallel pyrosequencing. First we evaluated the method with controls and found 41-485 reads per sequence in BRCA pathogenic mutations. Negative controls did not show deleterious variants, confirming the suitability of the approach. In patients diagnosed with cancer we found 4 novel deleterious mutations (c.2805_2808delAGAT and c.3124_3133delAGCAATATTA in BRCA1; c.2639_2640delTG and c.5114_5117delTAAA in BRCA2). The prevalence of BRCA mutations in these patients was 10.2%. Moreover, we discovered 16 variants with unknown clinical significance (11 in exons and 5 in introns); 4 were predicted as possibly pathogenic by in silico analyses, and 3 have not been described previously. This study illustrates how massive pyrosequencing technology can be applied to screen for BRCA mutations in the whole exonic and splice regions in patients with suspected BRCA-related cancers. This is the first effort to analyse the mutational status of BRCA genes on a Mexican-mestizo population by means of pyrosequencing.     

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.     

Linking C5 deficiency to an exonic splicing enhancer mutation

As an important component of the innate immune system, complement provides the initial response to prevent infections by pathogenic microorganisms. Patients with dysfunction of C5 display a propensity for severe recurrent infections. In this study, we present a patient with C5 deficiency demonstrated by immunochemical and functional analyses. Direct sequencing of all C5 exons displayed no mutation of obvious functional significance, except for an A to G transition in exon 10 predicting an exchange from lysine to arginine. This sequence alteration was present in only one allele of family members with a reduced serum C5 concentration and in both alleles of the patient with almost complete C5 deficiency, suggesting that this alteration may be producing the phenotype. Recent findings indicate that distinct nucleotide sequences, termed exonic splicing enhancers (ESEs), influence the splicing process. cDNA from all family members harboring the mutated allele showed skipping of exon 10, which resulted in a premature STOP codon, explaining the lack of C5 in the propositus. Sequence analysis of the mutated region revealed the substitution to be located within an ESE, as predicted by the RESCUE-ESE program. The altered ESE sequence is located close to the 5' splicing site and also lowers the predicted strength of the splice site itself. This apparently inconsequential sequence alteration represents a noncanonical splicing mutation altering an ESE. Our finding sheds a new light on the role of putative silent/conservative mutations in disease-associated genes.     

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.     

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.     

Binding of DAZAP1 and hnRNPA1/A2 to an exonic splicing silencer in a natural BRCA1 exon 18 mutant

A disease-causing G-to-T transversion at position +6 of BRCA1 exon 18 induces exclusion of the exon from the mRNA and, as has been suggested by in silico analysis, disrupts an ASF/SF2-dependent splicing enhancer. We show here using a pulldown assay with an internal standard that wild-type (WT) and mutant T6 sequences displayed similar ASF/SF2 binding efficiencies, which were significantly lower than that of a typical exonic splicing enhancer derived from the extra domain A exon of fibronectin. Overexpression or small interfering RNA (siRNA)-mediated depletion of ASF/SF2 did not affect the splicing of a WT BRCA1 minigene but resulted in an increase and decrease of T6 exon 18 inclusion, respectively. Furthermore, extensive mutation analysis using hybrid minigenes indicated that the T6 mutant creates a sequence with a prevalently inhibitory function. Indeed, RNA-protein interaction and siRNA experiments showed that the skipping of T6 BRCA1 exon 18 is due to the creation of a splicing factor-dependent silencer. This sequence specifically binds to the known repressor protein hnRNPA1/A2 and to DAZAP1, the involvement of which in splicing inhibition we have demonstrated. Our results indicate that the binding of the splicing factors hnRNPA1/A2 and DAZAP1 is the primary determinant of T6 BRCA1 exon 18 exclusion.     

A Retroelement Modifies Pre-mRNA Splicing: THE MURINE Glrbspa ALLELE IS A SPLICING SIGNAL POLYMORPHISM AMPLIFIED BY LONG INTERSPERSED NUCLEAR ELEMENT INSERTION

The glycine receptor-deficient mutant mouse spastic carries a full-length long interspersed nuclear element (LINE1) retrotransposon in intron 6 of the glycine receptor β subunit gene, Glrb(spa). The mutation arose in the C57BL/6J strain and is associated with skipping of exon 6 or a combination of the exons 5 and 6, thus resulting in a translational frameshift within the coding regions of the GlyR β subunit. The effect of the Glrb(spa) LINE1 insertion on pre-mRNA splicing was studied using a minigene approach. Sequence comparison as well as motif prediction and mutational analysis revealed that in addition to the LINE1 insertion the inactivation of an exonic splicing enhancer (ESE) within exon 6 is required for skipping of exon 6. Reconstitution of the ESE by substitution of a single residue was sufficient to prevent exon skipping. In addition to the ESE, two regions within the 5' and 3' UTR of the LINE1 were shown to be critical determinants for exon skipping, indicating that LINE1 acts as efficient modifier of subtle endogenous splicing phenotypes. Thus, the spastic allele of the murine glycine receptor β subunit gene is a two-hit mutation, where the hypomorphic alteration in an ESE is amplified by the insertion of a LINE1 element in the adjacent intron. Conversely, the LINE1 effect on splicing may be modulated by individual polymorphisms, depending on the insertional environment within the host genome.     

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.     

Two percent of men with early-onset prostate cancer harbor germline mutations in the BRCA2 gene

Studies of families with breast cancer have indicated that male carriers of BRCA2 mutations are at increased risk of prostate cancer, particularly at an early age. To evaluate the contribution of BRCA2 mutations to early-onset prostate cancer, we screened the complete coding sequence of BRCA2 for germline mutations, in 263 men with diagnoses of prostate cancer who were 相似文献   

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.     

Contribution of BRCA1 germline mutation in patients with sporadic breast cancer

Hereditary artifacts in BRCA1 gene have a significant contributory role in familial cases of breast cancer. However, its germline mutational penetrance in sporadic breast cancer cases with respect to Pakistani population has not yet been very well defined. This study was designed to assess the contributory role of germline mutations of this gene in sporadic cases of breast cancer. 150 cases of unilateral breast cancer patients, with no prior family history of breast cancer and no other disorders or diseases in general with age range 35–75 yrs, were included in this study.Mutational analysis for hot spots on Exon 2, 3 and 13 of BRCA1 was done by using Single Strand Conformational Polymorphism (SSCP). Sequence analysis revealed five variants (missense) and one novel splice site mutation at exon 13. No germline mutation was observed on the remaining exons with respect sporadic breast cancer cases in Pakistani population. A vast majority of breast cancer cases are sporadic; the present study may be helpful for designing a better genetic screening tool for germline BRCA mutations in sporadic breast cancer patients of Pakistani population. Further studies involving a screening of entire coding region of BRCA1 is required to explore the merits of genetic diagnosis and counseling in breast cancer patients.     

Interactions among SR proteins, an exonic splicing enhancer, and a lentivirus Rev protein regulate alternative splicing.

We examine here the roles of cellular splicing factors and virus regulatory proteins in coordinately regulating alternative splicing of the tat/rev mRNA of equine infectious anemia virus (EIAV). This bicistronic mRNA contains four exons; exons 1 and 2 encode Tat, and exons 3 and 4 encode Rev. In the absence of Rev expression, the four-exon mRNA is synthesized exclusively, but when Rev is expressed, exon 3 is skipped to produce an mRNA that contains only exons 1, 2, and 4. We identify a purine-rich exonic splicing enhancer (ESE) in exon 3 that promotes exon inclusion. Similar to other cellular ESEs that have been identified by other laboratories, the EIAV ESE interacted specifically with SR proteins, a group of serine/arginine-rich splicing factors that function in constitutive and alternative mRNA splicing. Substitution of purines with pyrimidines in the ESE resulted in a switch from exon inclusion to exon skipping in vivo and abolished binding of SR proteins in vitro. Exon skipping was also induced by expression of EIAV Rev. We show that Rev binds to exon 3 RNA in vitro, and while the precise determinants have not been mapped, Rev function in vivo and RNA binding in vitro indicate that the RNA element necessary for Rev responsiveness overlaps or is adjacent to the ESE. We suggest that EIAV Rev promotes exon skipping by interfering with SR protein interactions with RNA or with other splicing factors.     

Screening for genomic rearrangements in families with breast and ovarian cancer identifies BRCA1 mutations previously missed by conformation-sensitive gel electrophoresis or sequencing

The frequency of genomic rearrangements in BRCA1 was assessed in 42 American families with breast and ovarian cancer who were seeking genetic testing and who were subsequently found to be negative for BRCA1 and BRCA2 coding-region mutations. An affected individual from each family was tested by PCR for the exon 13 duplication (Puget et al. 1999a) and by Southern blot analysis for novel genomic rearrangements. The exon 13 duplication was detected in one family, and four families had other genomic rearrangements. A total of 5 (11. 9%) of the 42 families with breast/ovarian cancer who did not have BRCA1 and BRCA2 coding-region mutations had mutations in BRCA1 that were missed by conformation-sensitive gel electrophoresis or sequencing. Four of five families with BRCA1 genomic rearrangements included at least one individual with both breast and ovarian cancer; therefore, 4 (30.8%) of 13 families with a case of multiple primary breast and ovarian cancer had a genomic rearrangement in BRCA1. Families with genomic rearrangements had prior probabilities of having a BRCA1 mutation, ranging from 33% to 97% (mean 70%) (Couch et al. 1997). In contrast, in families without rearrangements, prior probabilities of having a BRCA1 mutation ranged from 7% to 92% (mean 37%). Thus, the prior probability of detecting a BRCA1 mutation may be a useful predictor when considering the use of Southern blot analysis for families with breast/ovarian cancer who do not have detectable coding-region mutations.     

BRCA1 mutations in a population-based study of breast cancer in Stockholm County

The mutation frequency of BRCA1 and BRCA2 in women with breast cancer varies according to family history, age at diagnosis and ethnicity. The contribution of BRCA1 and BRCA2 mutations in breast cancer populations, unselected for age and family history, has been examined in several studies reporting mutation frequencies between 1% and 12% by screening methods, population sizes, and to what extent the gene/s were screened differed in the studies. We wanted to clarify the proportion of breast cancer attributable to mutations in BRCA1 in an unselected breast cancer population from the Stockholm region. All incident breast cancer patients treated surgically in a 19-month period were eligible for the study and 70% (489/696) participated. Exon 11 of BRCA1 was screened for mutations using the protein truncation test, and the mutation frequency was estimated from that. In previous studies on high-risk families from Stockholm, more than 70% of the mutations were detected in exon 11. Two mutations were found, both in patients with a family history or their own medical history of ovarian cancer, giving a mutation frequency in exon 11 of 0.4% and an estimated BRCA1 mutation frequency of <1%. Mutations in BRCA1 in unselected breast cancer cases in our region are rare and likely to be found only in high-risk families. Our BRCA1 prevalence is the lowest of all studies on unselected breast cancer patients, probably reflecting the comparatively low rates detected also in high-risk breast cancer families from the region.     

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.