Selection and characterization of pre-mRNA splicing enhancers: identification of novel SR protein-specific enhancer sequences

SMN in spinal muscular atrophy and snRNP biogenesis

Ribonucleoprotein (RNP) complexes function in nearly every facet of cellular activity. The spliceosome is an essential RNP that accurately identifies introns and catalytically removes the intervening sequences, providing exquisite control of spatial, temporal, and developmental gene expressions. U-snRNPs are the building blocks for the spliceosome. A significant amount of insight into the molecular assembly of these essential particles has recently come from a seemingly unexpected area of research: neurodegeneration. Survival motor neuron (SMN) performs an essential role in the maturation of snRNPs, while the homozygous loss of SMN1 results in the development of spinal muscular atrophy (SMA), a devastating neurodegenerative disease. In this review, the function of SMN is examined within the context of snRNP biogenesis and evidence is examined which suggests that the SMN functional defects in snRNP biogenesis may account for the motor neuron pathology observed in SMA.

NSrp70 is a novel nuclear speckle-related protein that modulates alternative pre-mRNA splicing in vivo

Nuclear speckles are known to be the storage sites of mRNA splicing regulators. We report here the identification and characterization of a novel speckle protein, referred to as NSrp70, based on its subcellular localization and apparent molecular weight. This protein was first identified as CCDC55 by the National Institutes of Health Mammalian Gene Collection, although its function has not been assigned. NSrp70 was colocalized and physically interacted with SC35 and ASF/SF2 in speckles. NSrp70 has a putative RNA recognition motif, the RS-like region, and two coiled-coil domains, suggesting a role in RNA processing. Accordingly, using CD44, Tra2β1 and Fas constructs as splicing reporter minigenes, we found that NSrp70 modulated alternative splice site selection in vivo. The C-terminal 10 amino acids (531–540), including ⁵³⁶RD⁵³⁷, were identified as a novel nuclear localization signal, and the region spanning 290–471 amino acids was critical for speckle localization and binding to SC35 and ASF/SF2. The N-terminal region (107–161) was essential for the pre-mRNA splicing activity. Finally, we found that knockout of NSrp70 gene in mice led to a lack of progeny, including fetal embryos. Collectively, we demonstrate that NSrp70 is a novel splicing regulator and essentially required early stage of embryonic development.

Global analysis reveals SRp20- and SRp75-specific mRNPs in cycling and neural cells

Members of the SR protein family of RNA-binding proteins have numerous roles in mRNA metabolism, from transcription to translation. To understand how SR proteins coordinate gene regulation, comprehensive knowledge of endogenous mRNA targets is needed. Here we establish physiological expression of GFP-tagged SR proteins from stable transgenes. Using the GFP tag for immunopurification of mRNPs, mRNA targets of SRp20 and SRp75 were identified in cycling and neurally induced P19 cells. Genome-wide analysis showed that SRp20 and SRp75 associate with hundreds of distinct, functionally related groups of transcripts that change in response to neural differentiation. Knockdown of either SRp20 or SRp75 led to up- or downregulation of specific transcripts, including identified targets, and rescue by the GFP-tagged SR proteins proved their functionality. Thus, SR proteins contribute to the execution of gene-expression programs through their association with distinct endogenous mRNAs.

Stoichiometry of a regulatory splicing complex revealed by single-molecule analyses

Splicing is regulated by complex interactions of numerous RNA-binding proteins. The molecular mechanisms involved remain elusive, in large part because of ignorance regarding the numbers of proteins in regulatory complexes. Polypyrimidine tract-binding protein (PTB), which regulates tissue-specific splicing, represses exon 3 of alpha-tropomyosin through distant pyrimidine-rich tracts in the flanking introns. Current models for repression involve either PTB-mediated looping or the propagation of complexes between tracts. To test these models, we used single-molecule approaches to count the number of bound PTB molecules both by counting the number of bleaching steps of GFP molecules linked to PTB within complexes and by analysing their total emissions. Both approaches showed that five or six PTB molecules assemble. Given the domain structures, this suggests that the molecules occupy primarily multiple overlapping potential sites in the polypyrimidine tracts, excluding propagation models. As an alternative to direct looping, we propose that repression involves a multistep process in which PTB binding forms small local loops, creating a platform for recruitment of other proteins that bring these loops into close proximity.

HIV-1 mRNA 3' end processing is distinctively regulated by eIF3f, CDK11, and splice factor 9G8

A genetic screen previously identified the N-terminal 91 amino acids of the eukaryotic initiation factor 3 subunit f (N91-eIF3f) as a potent inhibitor of HIV-1 replication. Overexpression of N91-eIF3f or full-length eIF3f reduced the level of HIV-1 mRNAs in the infected cell. Here we show that N91-eIF3f and eIF3f act by specifically blocking the 3' end processing of the HIV-1 pre-mRNA both in vivo and in vitro. Furthermore, the results suggest that eIF3f mediates this restriction of HIV-1 expression through the previously unsuspected involvement of a set of factors that includes eIF3f, the SR protein 9G8, and the cyclin-dependent kinase 11 (CDK11). eIF3f affects HIV-1 3' end processing by modulating the sequence-specific recognition of the HIV-1 pre-mRNA by 9G8.

The SR protein family

The SR proteins are not only involved in pre-mRNA splicing but in mRNA export and the initiation of translation.

Summary

The processing of pre-mRNAs is a fundamental step required for the expression of most metazoan genes. Members of the family of serine/arginine (SR)-rich proteins are critical components of the machineries carrying out these essential processing events, highlighting their importance in maintaining efficient gene expression. SR proteins are characterized by their ability to interact simultaneously with RNA and other protein components via an RNA recognition motif (RRM) and through a domain rich in arginine and serine residues, the RS domain. Their functional roles in gene expression are surprisingly diverse, ranging from their classical involvement in constitutive and alternative pre-mRNA splicing to various post-splicing activities, including mRNA nuclear export, nonsense-mediated decay, and mRNA translation. These activities point up the importance of SR proteins during the regulation of mRNA metabolism.

Cryptic exon activation by disruption of exon splice enhancer: novel mechanism causing 3-methylcrotonyl-CoA carboxylase deficiency

3-Methylcrotonyl-CoA carboxylase (MCC) deficiency is an autosomal recessive disorder of leucine catabolism. MCC is a heteromeric mitochondrial enzyme composed of biotin-containing alpha (MCCA) and smaller beta (MCCB) subunits encoded by MCCA and MCCB, respectively. We report studies of the c.1054G-->A mutation in exon 11 of MCCB detected in the homozygous state in a patient with MCC deficiency. Sequence analysis of MCCB cDNA revealed two overlapping transcripts, one containing the normal 73 bp of exon 11 including the missense mutation c.1054G-->A (p.G352R), the other with exon 11 replaced by a 64-bp sequence from intron 10 (cryptic exon 10a) that maintains the reading frame and is flanked by acceptable splice consensus sites. In expression studies, we show that both transcripts lack detectable MCC activity. Western blot analysis showed slightly reduced levels of MCCB using the transcript containing the missense mutation, whereas no MCCB was detected with the transcript containing the cryptic exon 10a. Analysis of the region harboring the mutation revealed that the c.1054G-->A mutation is located in an exon splice enhancer sequence. Using MCCB minigene constructs to transfect MCCB-deficient fibroblasts, we demonstrate that the reduction in utilization of exon 11 associated with the c.1054G-->A mutation is due to alteration of this exon splice enhancer. Further, we show that optimization of the weak splice donor site of exon 11 corrects the splicing defect. To our knowledge, this is the first demonstration of a point mutation disrupting an exon splice enhancer that causes exon skipping along with utilization of a cryptic exon.

A deep intronic mutation in FGB creates a consensus exonic splicing enhancer motif that results in afibrinogenemia caused by aberrant mRNA splicing, which can be corrected in vitro with antisense oligonucleotide treatment

We previously described a novel homozygous point mutation (FGB c.115-600A>G) located deep within intron 1 of the fibrinogen beta gene (FGB), as a likely cause of afibrinogenemia. While this was the only mutation detected, its pathological mechanism was unclear. Here we show the mutation causes the inclusion of a 50-bp cryptic exon by creating a consensus heptad motif recognized by the spliceosome recruiting protein pre-mRNA splicing factor (SF2)/arginine/serine-rich alternative splicing factor (ASF) splicing factor 2/alternative splicing factor (SF2/ASF). Translation of the aberrant mRNA would result in truncation of the Bbeta chain, preventing fibrinogen synthesis. Selective introduction of a second mutation into the enhancer motif abolished the SF2/ASF binding motif and re-established normal pre-mRNA splicing. Subsequent introduction of antisense phosphorodiamidate morpholino oligonucleotides (PMOs) into transfected cells containing the mutant construct blocked the protein-RNA interaction and successfully restored normal splicing ( approximately 50% at 2 microM and approximately 90% at 10 microM). The molecular characterization of this case has revealed a unique disease mechanism, shown the importance of screening for deep intronic mutations, and provided evidence that antisense gene therapy is potentially practical for the treatment of diseases caused by this class of mutation.

The SR protein family of splicing factors: master regulators of gene expression

The SR protein family comprises a number of phylogenetically conserved and structurally related proteins with a characteristic domain rich in arginine and serine residues, known as the RS domain. They play significant roles in constitutive pre-mRNA splicing and are also important regulators of alternative splicing. In addition they participate in post-splicing activities, such as mRNA nuclear export, nonsense-mediated mRNA decay and mRNA translation. These wide-ranging roles of SR proteins highlight their importance as pivotal regulators of mRNA metabolism, and if these functions are disrupted, developmental defects or disease may result. Furthermore, animal models have shown a highly specific, non-redundant role for individual SR proteins in the regulation of developmental processes. Here, we will review the current literature to demonstrate how SR proteins are emerging as one of the master regulators of gene expression.

Ab initio prediction of mutation-induced cryptic splice-site activation and exon skipping

Mutations that affect splicing of precursor messenger RNAs play a major role in the development of hereditary diseases. Most splicing mutations have been found to eliminate GT or AG dinucleotides that define the 5' and 3' ends of introns, leading to exon skipping or cryptic splice-site activation. Although accurate description of the mis-spliced transcripts is critical for predicting phenotypic consequences of these alterations, their exact nature in affected individuals cannot often be determined experimentally. Using a comprehensive collection of exons that sustained cryptic splice-site activation or were skipped as a result of splice-site mutations, we have developed a multivariate logistic discrimination procedure that distinguishes the two aberrant splicing outcomes from DNA sequences. The new algorithm was validated using an independent sample of exons and implemented as a free online utility termed CRYP-SKIP (http://www.dbass.org.uk/cryp-skip/). The web application takes up one or more mutated alleles, each consisting of one exon and flanking intronic sequences, and provides a list of important predictor variables and their values, the overall probability of activating cryptic splice vs exon skipping, and the location and intrinsic strength of predicted cryptic splice sites in the input sequence. These results will facilitate phenotypic prediction of splicing mutations and provide further insights into splicing enhancer and silencer elements and their relative importance for splice-site selection in vivo.

SRp20 and CUG-BP1 modulate insulin receptor exon 11 alternative splicing

The insulin receptor (IR) exists as two isoforms, IR-A and IR-B, which result from alternative splicing of exon 11 in the primary transcript. This alternative splicing is cell specific, and the relative proportions of exon 11 isoforms also vary during development, aging, and different disease states. We have previously demonstrated that both intron 10 and exon 11 contain regulatory sequences that affect IR splicing both positively and negatively. In this study, we sought to define the precise sequence elements within exon 11 that control exon recognition and cellular factors that recognize these elements. Using minigenes carrying linker-scanning mutations within exon 11, we detected both exonic splicing enhancer and exonic splicing silencer elements. We identified binding of SRp20 and SF2/ASF to the exonic enhancers and CUG-BP1 to the exonic silencer by RNA affinity chromatography. Overexpression and knockdown studies with hepatoma and embryonic kidney cells demonstrated that SRp20 and SF2/ASF increase exon inclusion but that CUG-BP1 causes exon skipping. We found that CUG-BP1 also binds to an additional intronic splicing silencer, located at the 3' end of intron 10, to promote exon 11 skipping. Thus, we propose that SRp20, SF2/ASF, and CUG-BP1 act antagonistically to regulate IR alternative splicing in vivo and that the relative ratios of SRp20 and SF2/ASF to CUG-BP1 in different cells determine the degree of exon inclusion.

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

Background

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

Results

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

Conclusion

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

Automatic detection of exonic splicing enhancers (ESEs) using SVMs

Background

Exonic splicing enhancers (ESEs) activate nearby splice sites and promote the inclusion (vs. exclusion) of exons in which they reside, while being a binding site for SR proteins. To study the impact of ESEs on alternative splicing it would be useful to have a possibility to detect them in exons. Identifying SR protein-binding sites in human DNA sequences by machine learning techniques is a formidable task, since the exon sequences are also constrained by their functional role in coding for proteins.

Results

The choice of training examples needed for machine learning approaches is difficult since there are only few exact locations of human ESEs described in the literature which could be considered as positive examples. Additionally, it is unclear which sequences are suitable as negative examples. Therefore, we developed a motif-oriented data-extraction method that extracts exon sequences around experimentally or theoretically determined ESE patterns. Positive examples are restricted by heuristics based on known properties of ESEs, e.g. location in the vicinity of a splice site, whereas negative examples are taken in the same way from the middle of long exons. We show that a suitably chosen SVM using optimized sequence kernels (e.g., combined oligo kernel) can extract meaningful properties from these training examples. Once the classifier is trained, every potential ESE sequence can be passed to the SVM for verification. Using SVMs with the combined oligo kernel yields a high accuracy of about 90 percent and well interpretable parameters.

Conclusion

The motif-oriented data-extraction method seems to produce consistent training and test data leading to good classification rates and thus allows verification of potential ESE motifs. The best results were obtained using an SVM with the combined oligo kernel, while oligo kernels with oligomers of a certain length could be used to extract relevant features.

Genome-wide analyses of alternative splicing in plants: opportunities and challenges

Alternative splicing (AS) creates multiple mRNA transcripts from a single gene. While AS is known to contribute to gene regulation and proteome diversity in animals, the study of its importance in plants is in its early stages. However, recently available plant genome and transcript sequence data sets are enabling a global analysis of AS in many plant species. Results of genome analysis have revealed differences between animals and plants in the frequency of alternative splicing. The proportion of plant genes that have one or more alternative transcript isoforms is approximately 20%, indicating that AS in plants is not rare, although this rate is approximately one-third of that observed in human. The majority of plant AS events have not been functionally characterized, but evidence suggests that AS participates in important plant functions, including stress response, and may impact domestication and trait selection. The increasing availability of plant genome sequence data will enable larger comparative analyses that will identify functionally important plant AS events based on their evolutionary conservation, determine the influence of genome duplication on the evolution of AS, and discover plant-specific cis-elements that regulate AS. This review summarizes recent analyses of AS in plants, discusses the importance of further analysis, and suggests directions for future efforts.

Genetic and biochemical studies in Argentinean patients with variegate porphyria

BACKGROUND

A partial deficiency in Protoporphyrinogen oxidase (PPOX) produces the mixed disorder Variegate Porphyria (VP), the second acute porphyria more frequent in Argentina. Identification of patients with an overt VP is absolutely important because treatment depends on an accurate diagnosis but more critical is the identification of asymptomatic relatives to avoid acute attacks which may progress to death.

METHODS

We have studied at molecular level 18 new Argentinean patients biochemically diagnosed as VP. PPOX gene was amplified in one or in twelve PCR reactions. All coding exons, flanking intronic and promoter regions were manual or automatically sequenced. For RT-PCR studies RNA was retrotranscripted, amplified and sequenced. PPOX activity in those families carrying a new and uncharacterized mutation was performed.

RESULTS

All affected individuals harboured mutations in heterozygous state. Nine novel mutations and 3 already reported mutations were identified. Six of the novel mutations were single nucleotide substitutions, 2 were small deletions and one a small insertion. Three single nucleotide substitutions and the insertion were at exon-intron boundaries. Two of the single nucleotide substitutions, c.471G>A and c.807G>A and the insertion (c.388+3insT) were close to the splice donor sites in exons 5, 7 and intron 4 respectively. The other single nucleotide substitution was a transversion in the last base of intron 7, g.3912G>C (c.808-1G>C) so altering the consensus acceptor splice site. However, only in the first case the abnormal band showing the skipping of exon 5 was detected. The other single nucleotide substitutions were transversions: c.101A>T, c.995G>C and c.670 T>G that result in p.E34V, p.G332A and W224G aminoacid substitutions in exons 3, 10 and 7 respectively. Activity measurements indicate that these mutations reduced about 50% PPOX activity and also that they co-segregate with this reduced activity value. Two frameshift mutations, c.133delT and c.925delA, were detected in exons 3 and 9 respectively. The first leads to an early termination signal 22 codons downstream (p.S45fsX67) and the second leads to a stop codon 5 codons downstream (p.I309fsX314). One reported mutation was a missense mutation (p.G232R) and 2 were frameshift mutations: c.1082insC and 1043insT. The last mutation was detected in six new apparently unrelated Argentinean families.

CONCLUSION

Molecular analysis in available family members revealed 14 individuals who were silent carriers of VP. Molecular techniques represent the most accurate approach to identify unaffected carriers and to provide accurate genetic counselling for asymptomatic individuals. The initial screening includes the insertion search.

Searching for splicing motifs

Intron removal during pre-mRNA splicing in higher eukaryotes requires the accurate identification of the two splice sites at the ends of the exons, or exon definition. The sequences constituting the splice sites provide insufficient information to distinguish true splice sites from the greater number of false splice sites that populate transcripts. Additional information used for exon recognition resides in a large number of positively or negatively acting elements that lie both within exons and in the adjacent introns. The identification of such sequence motifs has progressed rapidly in recent years, such that extensive lists are now available for exonic splicing enhancers and exonic splicing silencers. These motifs have been identified both by empirical experiments and by computational predictions, the validity of the latter being confirmed by experimental verification. Molecular searches have been carried out either by the selection of sequences that bind to splicing factors, or enhance or silence splicing in vitro or in vivo. Computational methods have focused on sequences of 6 or 8 nucleotides that are over- or under-represented in exons, compared to introns or transcripts that do not undergo splicing. These various methods have sought to provide global definitions of motifs, yet the motifs are distinctive to the method used for identification and display little overlap. Astonishingly, at least three-quarters of a typical mRNA would be comprised of these motifs. A present challenge lies in understanding how the cell integrates this surfeit of information to generate what is usually a binary splicing decision.

The Alternative Splicing Mutation Database: a hub for investigations of alternative splicing using mutational evidence

Background

Some mutations in the internal regions of exons occur within splicing enhancers and silencers, influencing the pattern of alternative splicing in the corresponding genes. To understand how these sequence changes affect splicing, we created a database of these mutations.

Findings

The Alternative Splicing Mutation Database (ASMD) serves as a repository for all exonic mutations not associated with splicing junctions that measurably change the pattern of alternative splicing. In this initial published release (version 1.2), only human sequences are present, but the ASMD will grow to include other organisms, (see Availability and requirements section for the ASMD web address).

This relational database allows users to investigate connections between mutations and features of the surrounding sequences, including flanking sequences, RNA secondary structures and strengths of splice junctions. Splicing effects of the mutations are quantified by the relative presence of alternative mRNA isoforms with and without a given mutation. This measure is further categorized by the accuracy of the experimental methods employed. The database currently contains 170 mutations in 66 exons, yet these numbers increase regularly.

We developed an algorithm to derive a table of oligonucleotide Splicing Potential (SP) values from the ASMD dataset. We present the SP concept and tools in detail in our corresponding article.

Conclusion

The current data set demonstrates that mutations affecting splicing are located throughout exons and might be enriched within local RNA secondary structures. Exons from the ASMD have below average splicing junction strength scores, but the difference is small and is judged not to be significant.

Coevolutionary networks of splicing cis-regulatory elements

Accurate and efficient splicing of eukaryotic pre-mRNAs requires recognition by trans-acting factors of a complex array of cis-acting RNA elements. Here, we developed a generalized Bayesian network to model the coevolution of splicing cis elements in diverse eukaryotic taxa. Cross-exon but not cross-intron compensatory interactions between the 5' splice site (5'ss) and 3' splice site (3'ss) were observed in human/mouse, indicating that the exon is the primary evolutionary unit in mammals. Studied plants, fungi, and invertebrates exhibited exclusively cross-intron interactions, suggesting that intron definition drives evolution in these organisms. In mammals, 5'ss strength and the strength of several classes of exonic splicing silencers (ESSs) evolved in a correlated way, whereas specific exonic splicing enhancers (ESEs), including motifs associated with hTra2, SRp55, and SRp20, evolved in a compensatory manner relative to the 5'ss and 3'ss. Interactions between specific ESS or ESE motifs were not observed, suggesting that elements bound by different factors are not commonly interchangeable. Thus, the splicing elements defining exons coevolve in a way that preserves overall exon strength, allowing specific elements to substitute for loss or weakening of others.

Global control of aberrant splice-site activation by auxiliary splicing sequences: evidence for a gradient in exon and intron definition

Auxiliary splicing signals play a major role in the regulation of constitutive and alternative pre-mRNA splicing, but their relative importance in selection of mutation-induced cryptic or de novo splice sites is poorly understood. Here, we show that exonic sequences between authentic and aberrant splice sites that were activated by splice-site mutations in human disease genes have lower frequencies of splicing enhancers and higher frequencies of splicing silencers than average exons. Conversely, sequences between authentic and intronic aberrant splice sites have more enhancers and less silencers than average introns. Exons that were skipped as a result of splice-site mutations were smaller, had lower SF2/ASF motif scores, a decreased availability of decoy splice sites and a higher density of silencers than exons in which splice-site mutation activated cryptic splice sites. These four variables were the strongest predictors of the two aberrant splicing events in a logistic regression model. Elimination or weakening of predicted silencers in two reporters consistently promoted use of intron-proximal splice sites if these elements were maintained at their original positions, with their modular combinations producing expected modification of splicing. Together, these results show the existence of a gradient in exon and intron definition at the level of pre-mRNA splicing and provide a basis for the development of computational tools that predict aberrant splicing outcomes.

The impact of alternative splicing in vivo: mouse models show the way

Alternative splicing is widely believed to have a major impact on almost all biological processes since it increases proteome complexity and thereby controls protein function. Recently, gene targeting in mice has been used to create in vivo models to study the regulation and consequences of alternative splicing. The evidence accumulated so far argues for a nonredundant, highly specific role of individual splicing factors in mammalian development, and furthermore, demonstrates the importance of distinct protein isoforms in vivo. In this review, we will compare phenotypes of mouse models for alternative splicing to crystallize common themes and to put them into perspective with the available in vitro data.

Interplay between exonic splicing enhancers, mRNA processing, and mRNA surveillance in the dystrophic Mdx mouse

Background

Pre-mRNA splicing, the removal of introns from RNA, takes place within the spliceosome, a macromolecular complex composed of five small nuclear RNAs and a large number of associated proteins. Spliceosome assembly is modulated by the 5′ and 3′ splice site consensus sequences situated at the ends of each intron, as well as by exonic and intronic splicing enhancers/silencers recognized by SR and hnRNP proteins. Nonsense mutations introducing a premature termination codon (PTC) often result in the activation of cellular quality control systems that reduce mRNA levels or alter the mRNA splicing pattern. The mdx mouse, a commonly used genetic model for Duchenne muscular dystrophy (DMD), lacks dystrophin by virtue of a premature termination codon (PTC) in exon 23 that also severely reduces the level of dystrophin mRNA. However, the effect of the mutation on dystrophin RNA processing has not yet been described.

Methodology/Principal Finding

Using combinations of different biochemical and cellular assays, we found that the mdx mutation partially disrupts a multisite exonic splicing enhancer (ESE) that is recognized by a 40 kDa SR protein. In spite of the presence of an inefficient intron 22 3′ splice site containing the rare GAG triplet, the mdx mutation does not activate nonsense-associated altered splicing (NAS), but induces exclusively nonsense-mediated mRNA decay (NMD). Functional binding sites for SR proteins were also identified in exon 22 and 24, and in vitro experiments show that SR proteins can mediate direct association between exon 22, 23, and 24.

Conclusions/Significance

Our findings highlight the complex crosstalk between trans-acting factors, cis-elements and the RNA surveillance machinery occurring during dystrophin mRNA processing. Moreover, they suggest that dystrophin exon–exon interactions could play an important role in preventing mdx exon 23 skipping, as well as in facilitating the pairing of committed splice sites.

SR proteins and related factors in alternative splicing

SR proteins are a family of RNA binding proteins that contain a signature RS domain enriched with serine/arginine repeats. The RS domain is also found in many other proteins, which are collectively referred to as SR-related proteins. Several prototypical SR proteins are essential splicing factors, but the majority of RS domain-containing factors are characterized by their ability to alter splice site selection in vitro or in transfected cells. SR proteins and SR-related proteins are generally believed to modulate splice site selection via RNA recognition motif-mediated binding to exonic splicing enhancers and RS domain-mediated protein-protein and protein-RNA interactions during spliceosome assembly. However, the biological function of individual RS domain-containing splicing regulators is complex because of redundant as well as competitive functions, context-dependent effects and regulation by cotranscriptional and post-translational events. This chapter will focus on our current mechanistic understanding of alternative splicing regulation by SR proteins and SR-related proteins and will discuss some of the questions that remain to be addressed in future research.

Compensatory relationship between splice sites and exonic splicing signals depending on the length of vertebrate introns

BACKGROUND

The signals that determine the specificity and efficiency of splicing are multiple and complex, and are not fully understood. Among other factors, the relative contributions of different mechanisms appear to depend on intron size inasmuch as long introns might hinder the activity of the spliceosome through interference with the proper positioning of the intron-exon junctions. Indeed, it has been shown that the information content of splice sites positively correlates with intron length in the nematode, Drosophila, and fungi. We explored the connections between the length of vertebrate introns, the strength of splice sites, exonic splicing signals, and evolution of flanking exons.

RESULTS

A compensatory relationship is shown to exist between different types of signals, namely, the splice sites and the exonic splicing enhancers (ESEs). In the range of relatively short introns (approximately, < 1.5 kilobases in length), the enhancement of the splicing signals for longer introns was manifest in the increased concentration of ESEs. In contrast, for longer introns, this effect was not detectable, and instead, an increase in the strength of the donor and acceptor splice sites was observed. Conceivably, accumulation of A-rich ESE motifs beyond a certain limit is incompatible with functional constraints operating at the level of protein sequence evolution, which leads to compensation in the form of evolution of the splice sites themselves toward greater strength. In addition, however, a correlation between sequence conservation in the exon ends and intron length, particularly, in synonymous positions, was observed throughout the entire length range of introns. Thus, splicing signals other than the currently defined ESEs, i.e., potential new classes of ESEs, might exist in exon sequences, particularly, those that flank long introns.

CONCLUSION

Several weak but statistically significant correlations were observed between vertebrate intron length, splice site strength, and potential exonic splicing signals. Taken together, these findings attest to a compensatory relationship between splice sites and exonic splicing signals, depending on intron length.

Tau exon 6 is regulated by an intricate interplay of trans factors and cis elements, including multiple branch points

Tau is a microtubule-associated protein whose transcript undergoes complex regulated splicing in the mammalian nervous system. Exon 6 of the gene is an alternatively spliced cassette whose expression profile differs from that of the other tau regulated exons, implying the involvement of distinct regulatory factors. Previous work had established the existence and use of two additional 3' splice sites within exon 6 and the influence of splicing factors polypyrimidine binding protein (PTB) and U2AF on its splicing. The present work shows that exon 6 isoforms exist in distinct ratios in different compartments of the nervous system and that splicing of exon 6 is governed by multiple branch points, exonic cis elements and additional trans factors. Recent results show that tau exon 6 is specifically suppressed in the brains of people who suffer from myotonic dystrophy type 1. The understanding of how tau exon 6 splicing is regulated may give us insights into the disease.

Molecular basis of RNA recognition and TAP binding by the SR proteins SRp20 and 9G8

The sequence-specific RNA-binding proteins SRp20 and 9G8 are the smallest members of the serine- and arginine-rich (SR) protein family, well known for their role in splicing. They also play a role in mRNA export, in particular of histone mRNAs. We present the solution structures of the free 9G8 and SRp20 RNA recognition motifs (RRMs) and of SRp20 RRM in complex with the RNA sequence 5'CAUC3'. The SRp20-RNA structure reveals that although all 4 nt are contacted by the RRM, only the 5' cytosine is primarily recognized in a specific way. This might explain the numerous consensus sequences found by SELEX (systematic evolution of ligands by exponential enrichment) for the RRM of 9G8 and SRp20. Furthermore, we identify a short arginine-rich peptide adjacent to the SRp20 and 9G8 RRMs, which does not contact RNA but is necessary and sufficient for interaction with the export factor Tip-associated protein (TAP). Together, these results provide a molecular description for mRNA and TAP recognition by SRp20 and 9G8.

Inference of splicing regulatory activities by sequence neighborhood analysis

Sequence-specific recognition of nucleic-acid motifs is critical to many cellular processes. We have developed a new and general method called Neighborhood Inference (NI) that predicts sequences with activity in regulating a biochemical process based on the local density of known sites in sequence space. Applied to the problem of RNA splicing regulation, NI was used to predict hundreds of new exonic splicing enhancer (ESE) and silencer (ESS) hexanucleotides from known human ESEs and ESSs. These predictions were supported by cross-validation analysis, by analysis of published splicing regulatory activity data, by sequence-conservation analysis, and by measurement of the splicing regulatory activity of 24 novel predicted ESEs, ESSs, and neutral sequences using an in vivo splicing reporter assay. These results demonstrate the ability of NI to accurately predict splicing regulatory activity and show that the scope of exonic splicing regulatory elements is substantially larger than previously anticipated. Analysis of orthologous exons in four mammals showed that the NI score of ESEs, a measure of function, is much more highly conserved above background than ESE primary sequence. This observation indicates a high degree of selection for ESE activity in mammalian exons, with surprisingly frequent interchangeability between ESE sequences.

Gene expression involves a series of steps in which specific short DNA or RNA segments are recognized by nucleic acid–binding proteins. One step that is particularly prominent and complex in humans and other vertebrates is the removal of introns and the ligation of exons in the process of pre-mRNA splicing. To better understand the sequences in exons that regulate this process, the authors have developed a method termed Neighborhood Inference that predicts the splicing regulatory activity of RNA segments based on the known splicing enhancer or silencer activity of other segments that have closely neighboring sequences. This method is applied to predict hundreds of new exonic splicing regulatory elements, as well as splicing-neutral sequences. A number of these predictions were validated experimentally, indicating that the number of exonic splicing regulatory sequences is larger than previously suspected. Neighborhood Inference scoring is also used to show that selection on exonic splicing enhancers (ESEs) frequently allows conversion of one ESE sequence to another over evolutionary time periods, suggesting that ESEs are, to at least some degree, interchangeable in constitutively spliced exons. The methods described may also find application in the study of other biomolecular processes that involve sequence-specific nucleic acid–binding proteins.

A machine learning strategy to identify candidate binding sites in human protein-coding sequence

Background

The splicing of RNA transcripts is thought to be partly promoted and regulated by sequences embedded within exons. Known sequences include binding sites for SR proteins, which are thought to mediate interactions between splicing factors bound to the 5' and 3' splice sites. It would be useful to identify further candidate sequences, however identifying them computationally is hard since exon sequences are also constrained by their functional role in coding for proteins.

Results

This strategy identified a collection of motifs including several previously reported splice enhancer elements. Although only trained on coding exons, the model discriminates both coding and non-coding exons from intragenic sequence.

Conclusion

We have trained a computational model able to detect signals in coding exons which seem to be orthogonal to the sequences' primary function of coding for proteins. We believe that many of the motifs detected here represent binding sites for both previously unrecognized proteins which influence RNA splicing as well as other regulatory elements.

Alternative pre-mRNA splicing in the human system: unexpected role of repetitive sequences as regulatory elements

Alternative splicing is a process by which multiple messenger RNAs (mRNAs) are generated from a single pre-mRNA, resulting in functionally distinct protein products. This is accomplished by the differential recognition of splice sites in the pre-mRNA, often regulated in a tissue- or development-specific manner. Alternative splicing constitutes not only an important mechanism in controlling gene expression in humans, but also an essential source for increasing proteome diversity. In this review we summarize the underlying mechanistic principles, focussing on the cis-acting regulatory elements. In particular, the role of short sequence repeats, which are often polymorphic, in splicing regulation is discussed.

An increased specificity score matrix for the prediction of SF2/ASF-specific exonic splicing enhancers

Numerous disease-associated point mutations exert their effects by disrupting the activity of exonic splicing enhancers (ESEs). We previously derived position weight matrices to predict putative ESEs specific for four human SR proteins. The score matrices are part of ESEfinder, an online resource to identify ESEs in query sequences. We have now carried out a refined functional SELEX screen for motifs that can act as ESEs in response to the human SR protein SF2/ASF. The test BRCA1 exon under selection was internal, rather than the 3'-terminal IGHM exon used in our earlier studies. A naturally occurring heptameric ESE in BRCA1 exon 18 was replaced with two libraries of random sequences, one seven nucleotides in length, the other 14. Following three rounds of selection for in vitro splicing via internal exon inclusion, new consensus motifs and score matrices were derived. Many winner sequences were demonstrated to be functional ESEs in S100-extract-complementation assays with recombinant SF2/ASF. Motif-score threshold values were derived from both experimental and statistical analyses. Motif scores were shown to correlate with levels of exon inclusion, both in vitro and in vivo. Our results confirm and extend our earlier data, as many of the same motifs are recognized as ESEs by both the original and our new score matrix, despite the different context used for selection. Finally, we have derived an increased specificity score matrix that incorporates information from both of our SF2/ASF-specific matrices and that accurately predicts the exon-skipping phenotypes of deleterious point mutations.

The Tgif2 gene contains a retained intron within the coding sequence

Background

TGIF and TGIF2 are homeodomain proteins, which act as TGFβ specific Smad transcriptional corepressors. TGIF recruits general repressors including mSin3 and CtBP. The related TGIF2 protein functions in a similar manner, but does not bind CtBP. In addition to repressing TGFβ activated gene expression, TGIF and TGIF2 repress gene expression by binding directly to DNA. TGIF and TGIF2 share two major blocks of similarity, encompassing the homeodomain, and a conserved carboxyl terminal repression domain. Here we characterize two splice variants of the Tgif2 gene from mouse and demonstrate that the Tgif2 gene contains a retained intron.

Results

By PCR from mouse cDNA, we identified two alternate splice forms of the Tgif2 gene. One splice variant encodes the full length 237 amino acid Tgif2, whereas the shorter form results in the removal of 39 codons from the centre of the coding region. The generation of this alternate splice form occurs with the mouse RNA, but not the human, and both splice forms are present in all mouse tissues analyzed. Human and mouse Tgif2 coding sequences contain a retained intron, which in mouse Tgif2 is removed by splicing from around 25–50% of RNAs, as assessed by RT-PCR. This splicing event is dependent on sequences within the mouse Tgif2 coding sequence. Both splice forms of mouse Tgif2 encode proteins which are active transcriptional repressors, and can repress both TGFβ dependent and independent transcription. In addition, we show that human and mouse Tgif2 interact with the transcriptional corepressor mSin3.

Conclusion

These data demonstrate that the Tgif2 gene contains a retained intron, within the second coding exon. This retained intron is not removed from the human mRNA at a detectable level, but is spliced out in a significant proportion of mouse RNAs. This alternate splicing is dependent entirely on sequences within the mouse Tgif2 coding sequence, suggesting the presence of an exonic splicing enhancer. Both splice forms of mouse Tgif2 produce proteins which are functional transcriptional repressors.

Mutational screening of theCYP26A1 gene in patients with caudal regression syndrome

BACKGROUND

The retinoic acid (RA)-catabolizing enzyme Cyp26a1 plays an important role in protecting tailbud tissues from inappropriate exposure to RA. Cyp26a1-null animals exhibit caudal agenesis and spina bifida, imperforate anus, agenesis of the caudal portions of the digestive and urogenital tracts, and malformed lumbosacral skeletal elements. This phenotype closely resembles the most severe form of caudal agenesis in humans. In view of these findings, we investigated a potential involvement of the human CYP26A1 gene in the pathogenesis of caudal regression syndrome (CRS).

METHODS

Mutational screening of 49 CRS patients and 132 controls was performed using denaturing high-performance liquid chromatography and sequencing. Differences in the genotype and allele frequency of each SNP were evaluated by chi(2) analysis. The biological significance of the intronic variants was investigated by transfection assays of mutant constructs and by analysis of the splicing patterns with RT-PCR.

RESULTS

Mutational screening allowed us to identify 6 SNPs, 4 of which (447 C>G, 1134 G>A, IVS 1+10 G>C, and IVS 4+8 AG>GA) are new. In addition, we describe a novel 2-site haplotype consisting of the 2 intronic SNPs. Both single-locus and haplotype analyses revealed no association with increased risk for CRS. The consequences of the 2 intronic polymorphisms on the mRNA splicing process were also investigated. Moreover, using functional and computational methods we demonstrated that both of these intronic polymorphisms affect the intron splicing efficiency.

CONCLUSIONS

Our research did not provide evidence that CYP26A1 has implications for the pathogenesis of human CRS. However, the relationship between CRS risk and the CYP26A1 genotype requires further study with a larger number of genotyped subjects.

Systematic study of sequence motifs for RNA trans splicing in Trypanosoma brucei

mRNA maturation in Trypanosoma brucei depends upon trans splicing, and variations in trans-splicing efficiency could be an important step in controlling the levels of individual mRNAs. RNA splicing requires specific sequence elements, including conserved 5' splice sites, branch points, pyrimidine-rich regions [poly(Y) tracts], 3' splice sites (3'SS), and sometimes enhancer elements. To analyze sequence requirements for efficient trans splicing in the poly(Y) tract and around the 3'SS, we constructed a luciferase-beta-galactosidase double-reporter system. By testing approximately 90 sequences, we demonstrated that the optimum poly(Y) tract length is approximately 25 nucleotides. Interspersing a purely uridine-containing poly(Y) tract with cytidine resulted in increased trans-splicing efficiency, whereas purines led to a large decrease. The position of the poly(Y) tract relative to the 3'SS is important, and an AC dinucleotide at positions -3 and -4 can lead to a 20-fold decrease in trans splicing. However, efficient trans splicing can be restored by inserting a second AG dinucleotide downstream, which does not function as a splice site but may aid in recruitment of the splicing machinery. These findings should assist in the development of improved algorithms for computationally identifying a 3'SS and help to discriminate noncoding open reading frames from true genes in current efforts to annotate the T. brucei genome.

The origins of eukaryotic gene structure

Most of the phenotypic diversity that we perceive in the natural world is directly attributable to the peculiar structure of the eukaryotic gene, which harbors numerous embellishments relative to the situation in prokaryotes. The most profound changes include introns that must be spliced out of precursor mRNAs, transcribed but untranslated leader and trailer sequences (untranslated regions), modular regulatory elements that drive patterns of gene expression, and expansive intergenic regions that harbor additional diffuse control mechanisms. Explaining the origins of these features is difficult because they each impose an intrinsic disadvantage by increasing the genic mutation rate to defective alleles. To address these issues, a general hypothesis for the emergence of eukaryotic gene structure is provided here. Extensive information on absolute population sizes, recombination rates, and mutation rates strongly supports the view that eukaryotes have reduced genetic effective population sizes relative to prokaryotes, with especially extreme reductions being the rule in multicellular lineages. The resultant increase in the power of random genetic drift appears to be sufficient to overwhelm the weak mutational disadvantages associated with most novel aspects of the eukaryotic gene, supporting the idea that most such changes are simple outcomes of semi-neutral processes rather than direct products of natural selection. However, by establishing an essentially permanent change in the population-genetic environment permissive to the genome-wide repatterning of gene structure, the eukaryotic condition also promoted a reliable resource from which natural selection could secondarily build novel forms of organismal complexity. Under this hypothesis, arguments based on molecular, cellular, and/or physiological constraints are insufficient to explain the disparities in gene, genomic, and phenotypic complexity between prokaryotes and eukaryotes.

Cooperative actions of Tra2alpha with 9G8 and SRp30c in the RNA splicing of the gonadotropin-releasing hormone gene transcript

In earlier studies, we demonstrated that excision of the first intron (intron A) from the gonadotropin-releasing hormone (GnRH) transcript is highly cell type- and developmental stage-specific. The removal of GnRH intron A requires exonic splicing enhancers on exons 3 and 4 (ESE3 and ESE4, respectively). Tra2alpha,a serine/arginine-rich (SR)-like protein, specifically binds to ESE4, although it requires additional nuclear co-factors for efficient removal of this intron. In the present study, we demonstrate the cooperative action of multiple SR proteins in the regulation of GnRH pre-mRNA splicing. SRp30c specifically binds to both ESE3 and ESE4, whereas 9G8 binds to an element in exon 3 and strongly enhances the excision of GnRH intron A in the presence of minimal amount of other nuclear components. Interestingly, Tra2alpha can interact with either 9G8 or SRp30c, whereas no interaction between 9G8 and SRp30c is observed. Tra2alpha has an additive effect on the RNA binding of these proteins. Overexpression or knock-down of these three proteins in cultured cells further suggests their essential role in intron A excision activities, and their presence in GnRH neurons of the mouse preoptic area further strengthens this possibility. Together, these results indicate that interaction of Tra2alpha with 9G8 and SRp30c appears to be crucial for ESE-dependent GnRH pre-mRNA splicing, allowing efficient generation of mature mRNA in GnRH-producing cells.

Evidence for purifying selection against synonymous mutations in mammalian exonic splicing enhancers

Silent sites in mammals have classically been assumed to be free from selective pressures. Consequently, the synonymous substitution rate (Ks) is often used as a proxy for the mutation rate. Although accumulating evidence demonstrates that the assumption is not valid, the mechanism by which selection acts remain unclear. Recent work has revealed that the presence of exonic splicing enhancers (ESEs) in coding sequence might influence synonymous evolution. ESEs are predominantly located near intron-exon junctions, which may explain the reduced single-nucleotide polymorphism (SNP) density in these regions. Here we show that synonymous sites in putative ESEs evolve more slowly than the remaining exonic sequence. Differential mutabilities of ESEs do not appear to explain this difference. We observe that substitution frequency at fourfold synonymous sites decreases as one approaches the ends of exons, consistent with the existing SNP data. This gradient is at least in part explained by ESEs being more abundant near junctions. Between-gene variation in Ks is hence partly explained by the proportion of the gene that acts as an ESE. Given the relative abundance of ESEs and the reduced rates of synonymous divergence within them, we estimate that constraints on synonymous evolution within ESEs causes the true mutation rate to be underestimated by not more than approximately 8%. We also find that Ks outside of ESEs is much lower in alternatively spliced exons than in constitutive exons, implying that other causes of selection on synonymous mutations exist. Additionally, selection on ESEs appears to affect nonsynonymous sites and may explain why amino acid usage near intron-exon junctions is nonrandom.

Exon inclusion is dependent on predictable exonic splicing enhancers

We have previously formulated a list of approximately 2,000 RNA octamers as putative exonic splicing enhancers (PESEs) based on a statistical comparison of human exonic and nonexonic sequences (X. H. Zhang and L. A. Chasin, Genes Dev. 18:1241-1250, 2004). When inserted into a poorly spliced test exon, all eight tested octamers stimulated splicing, a result consistent with their identification as exonic splicing enhancers (ESEs). Here we present a much more stringent test of the validity of this list of PESEs. Twenty-two naturally occurring examples of nonoverlapping PESEs or PESE clusters were identified in six mammalian exons; five of the six exons tested are constitutively spliced. Each of the 22 individual PESEs or PESE clusters was disrupted by site-directed mutagenesis, usually by a single-base substitution. Eighteen of the 22 disruptions (82%) resulted in decreased splicing efficiency. In contrast, 24 control mutations had little or no effect on splicing. This high rate of success suggests that most PESEs function as ESEs in their natural context. Like most exons, these exons contain several PESEs. Since knocking out any one of several could produce a severalfold decrease in splicing efficiency, we conclude that there is little redundancy among ESEs in an exon and that they must work in concert to optimize splicing.

Distribution of SR protein exonic splicing enhancer motifs in human protein-coding genes

Exonic splicing enhancers (ESEs) are pre-mRNA cis-acting elements required for splice-site recognition. We previously developed a web-based program called ESEfinder that scores any sequence for the presence of ESE motifs recognized by the human SR proteins SF2/ASF, SRp40, SRp55 and SC35 (). Using ESEfinder, we have undertaken a large-scale analysis of ESE motif distribution in human protein-coding genes. Significantly higher frequencies of ESE motifs were observed in constitutive internal protein-coding exons, compared with both their flanking intronic regions and with pseudo exons. Statistical analysis of ESE motif frequency distributions revealed a complex relationship between splice-site strength and increased or decreased frequencies of particular SR protein motifs. Comparison of constitutively and alternatively spliced exons demonstrated slightly weaker splice-site scores, as well as significantly fewer ESE motifs, in the alternatively spliced group. Our results underline the importance of ESE-mediated SR protein function in the process of exon definition, in the context of both constitutive splicing and regulated alternative splicing.

A novel cryptic exon identified in the 3′ region of intron 2 of the human dystrophin gene

The dystrophin gene, which is mutated in Duchenne muscular dystrophy (DMD), is the largest known human gene and is characterized by the huge size of its introns. Intron 2, the second largest intron, is 170-kb long and has been shown to include a 140-bp cryptic exon (exon 2a) in its 5' region. The rest of this intron has no known function. In this study, we find that another cryptic exon, located in the 3' region of intron 2, is activated in a promoter- or tissue-specific manner. An unknown 98-bp insertion precisely between exons 2 and 3 was identified in one of the dystrophin mRNAs from lymphocytes of a DMD patient with a duplication of exon 2. This 98-bp sequence, located in the 3' region of intron 2, was found to possess a branch point, acceptor and donor splice-site consensus sequences, and an exonic splicing enhancer sequence, and thus is a novel exon, which we named "exon 2b." In lymphocytes, exon 2b incorporation was detected in the muscle-specific, promoter-driven transcript. Five of 20 normal human tissue mRNAs, including cardiac and skeletal muscle mRNAs, were confirmed to contain a fragment extending from exon 1 to exon 2b by reverse transcription PCR amplification, indicating that exon 2b is activated in a tissue-specific manner. This provides a clue to a novel cause of dystrophinopathy.

The switch in alternative splicing of cyclic AMP-response element modulator protein CREM{tau}2{alpha} (activator) to CREM{alpha} (repressor) in human myometrial cells is mediated by SRp40

The transcription factor cAMP-response element modulator (CREM) protein, plays a major role in cAMP-responsive gene regulation. Biological consequences resulting from the transcriptional stimuli of CREM are dictated by the expression of multiple protein isoforms generated by extensive alternative splicing of its precursor mRNA. We have previously shown that alternative splicing enables the expression of the CREM gene to be "switched" within the human myometrium during pregnancy from the production of CREMtau(2alpha), a potent transcriptional activator to the synthesis of CREMalpha, a transcriptional repressor. Furthermore we have recently reported that this change in the expression of CREM spliced variants is likely to have important ramifications on the regulation of downstream cAMP-response element-responsive target genes involved in uterine activity during gestation. We have investigated the splicing factors involved in controlling the expression of myometrial CREM splice variants. Data presented here from transient transfections indicate that the switch in the synthesis of CREMtau(2)alpha to CREMalpha that occurs during pregnancy is regulated primarily by an SR protein family member, SRp40. We also show that expression of this splicing factor is tightly regulated in the myometrium during pregnancy. SRp40 regulates the splicing of CREM via its interactions with multiple ESE motifs present in the alternatively exons of CREM. In vitro splicing and electrophoretic mobility shift assays were employed to confirm the functionality of the SRp40-binding ESEs, thus providing a mechanistic explanation of how SRp40 regulates the switch in splicing from production of CREMtau(2)alpha to CREMalpha.

Familial adenomatous polyposis: aberrant splicing due to missense or silent mutations in the APC gene

Familial adenomatous polyposis (FAP) is caused by germline mutations in the tumor suppressor gene APC. To date, the relevance of rare exonic single-base substitutions at nucleotide positions close to splice sites that are predicted to result in missense or silent (SNP) variants or substitutions in introns at splice-site positions that are not highly conserved has not been systematically examined in FAP patients. In 34 index patients, we identified 26 different heterozygous single-base substitutions at or close to the splice sites. We characterized five exonic mutations in exon 4 (c.423G>T), exon 14 (c.1956C>T, c.1957A>G, and c.1957A>C), and exon 15 (c.1959G>A) by transcript analysis and by splice-prediction programs (BDGP, SpliceSiteFinder, and ESEfinder). The splicing patterns of these variants were compared to those of 16 different substitutions at highly or less-conserved intronic splice-site positions, and to normal controls. In addition, we analyzed cosegregation of the variants with affected family members and examined the genotype-phenotype correlation. We could demonstrate that the four unclear variants in exon 4 and 14 that are predicted to result in missense or silent mutations in fact lead to complete exon skipping due to aberrant splicing; one possible explanation for this observed effect might be the disruption of exonic splicing enhancer (ESE) motifs. In contrast, the substitution at the first position of exon 15 seems to actually be a silent variant. We present the first systematic evaluation of different single-base substitutions in APC at or close to splice sites at transcript level. We show that the consequence of exonic mutations cannot be evaluated only by the predicted change in amino acid sequence but rather by the change at DNA level. The functional analysis of variants with unknown pathogenic effect plays an important role in increasing the mutation detection rate and achieving validation of predictive testing.

Understanding alternative splicing: towards a cellular code

In violation of the 'one gene, one polypeptide' rule, alternative splicing allows individual genes to produce multiple protein isoforms - thereby playing a central part in generating complex proteomes. Alternative splicing also has a largely hidden function in quantitative gene control, by targeting RNAs for nonsense-mediated decay. Traditional gene-by-gene investigations of alternative splicing mechanisms are now being complemented by global approaches. These promise to reveal details of the nature and operation of cellular codes that are constituted by combinations of regulatory elements in pre-mRNA substrates and by cellular complements of splicing regulators, which together determine regulated splicing pathways.

A duplication in the canine beta-galactosidase gene GLB1 causes exon skipping and GM1-gangliosidosis in Alaskan huskies

GM(1)-gangliosidosis is a lysosomal storage disease that is inherited as an autosomal recessive disorder, predominantly caused by structural defects in the beta-galactosidase gene (GLB1). The molecular cause of GM(1)-gangliosidosis in Alaskan huskies was investigated and a novel 19-bp duplication in exon 15 of the GLB1 gene was identified. The duplication comprised positions +1688-+1706 of the GLB1 cDNA. It partially disrupted a potential exon splicing enhancer (ESE), leading to exon skipping in a fraction of the transcripts. Thus, the mutation caused the expression of two different mRNAs from the mutant allele. One transcript contained the complete exon 15 with the 19-bp duplication, while the other transcript lacked exon 15. In the transcript containing exon 15 with the 19-bp duplication a premature termination codon (PTC) appeared, but due to its localization in the last exon of canine GLB1, nonsense-mediated RNA decay (NMD) did not occur. As a consequence of these molecular events two different truncated GLB1 proteins are predicted to be expressed from the mutant GLB1 allele. In heterozygous carrier animals the wild-type allele produces sufficient amounts of the active enzyme to prevent clinical signs of disease. In affected homozygous dogs no functional GLB1 is synthesized and G(M1)-gangliosidosis occurs.

Dichotomous splicing signals in exon flanks

Intronic elements flanking the splice-site consensus sequences are thought to play a role in pre-mRNA splicing. However, the generality of this role, the catalog of effective sequences, and the mechanisms involved are still lacking. Using molecular genetic tests, we first showed that the approximately 50-nt intronic flanking sequences of exons beyond the splice-site consensus are generally important for splicing. We then went on to characterize exon flank sequences on a genomic scale. The G+C content of flanks displayed a bimodal distribution reflecting an exaggeration of this base composition in flanks relative to the gene as a whole. We divided all exons into two classes according to their flank G+C content and used computational and statistical methods to define pentamers of high relative abundance and phylogenetic conservation in exon flanks. Upstream pentamers were often common to the two classes, whereas downstream pentamers were totally different. Upstream and downstream pentamers were often identical around low G+C exons, and in contrast, were often complementary around high G+C exons. In agreement with this complementarity, predicted base pairing was more frequent between the flanks of high G+C exons. Pseudo exons did not exhibit this behavior, but rather tended to form base pairs between flanks and exon bodies. We conclude that most exons require signals in their immediate flanks for efficient splicing. G+C content is a sequence feature correlated with many genetic and genomic attributes. We speculate that there may be different mechanisms for splice site recognition depending on G+C content.

A family based study implicates solute carrier family 1-member 3 (SLC1A3) gene in attention-deficit/hyperactivity disorder

BACKGROUND

The glutamatergic system, the major excitatory neurotransmitter system in the central nervous system (CNS) has been proposed as contributing a possible role in the etiology of attention deficit hyperactivity disorder (ADHD). This is based upon observations from animal, neuroimaging, neuroanatomical and neuropsychological studies. Genes related to glutamate function are therefore good functional candidates for this disorder. The SLC1A3 (Solute Carrier Family 1, member 3) gene encodes a glial glutamate transporter which maps to chromosome 5p12, a region of linkage that coincides in two published ADHD genome scans so far. SLC1A3 is thus both a functional and positional candidate gene for ADHD.

METHODS

We have undertaken detailed association analysis of SLC1A3 using a multi-stage approach for candidate gene analysis.

RESULTS

In a family-based sample (n = 299) we found a significant association between marker rs2269272 (p = .007) and ADHD. Two, two-marker haplotypes, rs2269272/rs3776581 (p = .016) and rs2269272/rs2032893 (p = .013) also yielded evidence of association.

CONCLUSIONS

The results of our study suggest that genetic variation in SLC1A3 may contribute to susceptibility to ADHD.

Serine/arginine-rich protein-dependent suppression of exon skipping by exonic splicing enhancers

The 5' and 3' splice sites within an intron can, in principle, be joined to those within any other intron during pre-mRNA splicing. However, exons are joined in a strict 5' to 3' linear order in constitutively spliced pre-mRNAs. Thus, specific mechanisms must exist to prevent the random joining of exons. Here we report that insertion of exon sequences into an intron can inhibit splicing to the downstream 3' splice site and that this inhibition is independent of intron size. The exon sequences required for splicing inhibition were found to be exonic enhancer elements, and their inhibitory activity requires the binding of serine/arginine-rich splicing factors. We conclude that exonic enhancers can act as barriers to prevent exon skipping and thereby may play a key role in ensuring the correct 5' to 3' linear order of exons in spliced mRNA.

Advances in the study of SR protein family

The name of SR proteins is derived from their typical RS domain that is rich in serine (Ser, S) and arginine (Arg, R). They are conserved in evolution. Up to now, 10 members of the SR protein family have been identified in humans. SR proteins contain one or two RNA binding motifs aside from the RS domain, and also possess special biochemical and immunological features. As to the functions of SR proteins, they facilitate the recruitment of the components of splicesome via protein-protein interaction to prompt the assembly of early splicesome; while in alternative splicing, tissue-specifically expressed SR protein along with the relative ratio of SR protein and heterogeneous nuclear ribonucleoprotein (hnRNP) is composed of two main regulative mechanisms to alternative splicing. Almost all of the biochemical functions are regulated by reversible phosphorylation.

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

Previous studies have identified UGCAUG as an intron splicing enhancer that is frequently located adjacent to tissue-specific alternative exons in the human genome. Here, we show that UGCAUG is phylogenetically and spatially conserved in introns that flank brain-enriched alternative exons from fish to man. Analysis of sequence from the mouse, rat, dog, chicken and pufferfish genomes revealed a strongly statistically significant association of UGCAUG with the proximal intron region downstream of brain-enriched alternative exons. The number, position and sequence context of intronic UGCAUG elements were highly conserved among mammals and in chicken, but more divergent in fish. Control datasets, including constitutive exons and non-tissue-specific alternative exons, exhibited a much lower incidence of closely linked UGCAUG elements. We propose that the high sequence specificity of the UGCAUG element, and its unique association with tissue-specific alternative exons, mark it as a critical component of splicing switch mechanism(s) designed to activate a limited repertoire of splicing events in cell type-specific patterns. We further speculate that highly conserved UGCAUG-binding protein(s) related to the recently described Fox-1 splicing factor play a critical role in mediating this specificity.