Searching for splicing motifs

Evolution of splicing regulatory networks in Drosophila

The proteome expanding effects of alternative pre-mRNA splicing have had a profound impact on eukaryotic evolution. The events that create this diversity can be placed into four major classes: exon skipping, intron retention, alternative 5′ splice sites, and alternative 3′ splice sites. Although the regulatory mechanisms and evolutionary pressures among alternative splicing classes clearly differ, how these differences affect the evolution of splicing regulation remains poorly characterized. We used RNA-seq to investigate splicing differences in D. simulans, D. sechellia, and three strains of D. melanogaster. Regulation of exon skipping and tandem alternative 3′ splice sites (NAGNAGs) were more divergent than other splicing classes. Splicing regulation was most divergent in frame-preserving events and events in noncoding regions. We further determined the contributions of cis- and trans-acting changes in splicing regulatory networks by comparing allele-specific splicing in F1 interspecific hybrids, because differences in allele-specific splicing reflect changes in cis-regulatory element activity. We find that species-specific differences in intron retention and alternative splice site usage are primarily attributable to changes in cis-regulatory elements (median ∼80% cis), whereas species-specific exon skipping differences are driven by both cis- and trans-regulatory divergence (median ∼50% cis). These results help define the mechanisms and constraints that influence splicing regulatory evolution and show that networks regulating the four major classes of alternative splicing diverge through different genetic mechanisms. We propose a model in which differences in regulatory network architecture among classes of alternative splicing affect the evolution of splicing regulation.

Regulation of alternative pre-mRNA splicing

Alternative splicing plays a prevalent role in generating functionally diversified proteomes from genomes with a more limited repertoire of protein-coding genes. Alternative splicing is frequently regulated with cell type or developmental specificity and in response to signaling pathways, and its mis-regulation can lead to disease. Co-regulated programs of alternative splicing involve interplay between a host of cis-acting transcript features and trans-acting RNA-binding proteins. Here, we review the current state of understanding of the logic and mechanism of regulated alternative splicing and indicate how this understanding can be exploited to manipulate splicing for therapeutic purposes.

Conservation/Mutation in the splice sites of cytokine receptor genes of mouse and human

Conservation/mutation in the intronic initial and terminal hexanucleotides was studied in 26 orthologous cytokine receptor genes of Mouse and Human. Introns began and ended with the canonical dinucleotides GT and AG, respectively. Identical configurations were found in 57% of the 5′ hexanucleotides and 28% of the 3′ hexanucleotides. The actual conservation percentages of the individual variable nucleotides at each position in the hexanucleotides were determined, and the theoretical rates of conservation of groups of three nucleotides were calculated under the hypothesis of a mutual evolutionary independence of the neighboring nucleotides (random association). Analysis of the actual conservation of groups of variable nucleotides showed that, at 5′, GTGAGx was significantly more expressed and GTAAGx was significantly less expressed, as compared to the random association. At 3′, TTTxAG and xTGCAG were overexpressed as compared to a random association. Study of Mouse and Human transcript variants involving the splice sites showed that most variants were not inherited from the common ancestor but emerged during the process of speciation. In some variants the silencing of a terminal hexanucleotide determined skipping of the downstream exon; in other variants the constitutive splicing hexanucleotide was replaced by another potential, in-frame, splicing hexanucleotide, leading to alterations of exon lengths.

Functional analysis of deep intronic SNP rs13438494 in intron 24 of PCLO gene

The single nucleotide polymorphism (SNP) rs13438494 in intron 24 of PCLO was significantly associated with bipolar disorder in a meta-analysis of genome-wide association studies. In this study, we performed functional minigene analysis and bioinformatics prediction of splicing regulatory sequences to characterize the deep intronic SNP rs13438494. We constructed minigenes with A and C alleles containing exon 24, intron 24, and exon 25 of PCLO to assess the genetic effect of rs13438494 on splicing. We found that the C allele of rs13438494 reduces the splicing efficiency of the PCLO minigene. In addition, prediction analysis of enhancer/silencer motifs using the Human Splice Finder web tool indicated that rs13438494 induces the abrogation or creation of such binding sites. Our results indicate that rs13438494 alters splicing efficiency by creating or disrupting a splicing motif, which functions by binding of splicing regulatory proteins, and may ultimately result in bipolar disorder in affected people.

Complexity of the alternative splicing landscape in plants

Alternative splicing (AS) of precursor mRNAs (pre-mRNAs) from multiexon genes allows organisms to increase their coding potential and regulate gene expression through multiple mechanisms. Recent transcriptome-wide analysis of AS using RNA sequencing has revealed that AS is highly pervasive in plants. Pre-mRNAs from over 60% of intron-containing genes undergo AS to produce a vast repertoire of mRNA isoforms. The functions of most splice variants are unknown. However, emerging evidence indicates that splice variants increase the functional diversity of proteins. Furthermore, AS is coupled to transcript stability and translation through nonsense-mediated decay and microRNA-mediated gene regulation. Widespread changes in AS in response to developmental cues and stresses suggest a role for regulated splicing in plant development and stress responses. Here, we review recent progress in uncovering the extent and complexity of the AS landscape in plants, its regulation, and the roles of AS in gene regulation. The prevalence of AS in plants has raised many new questions that require additional studies. New tools based on recent technological advances are allowing genome-wide analysis of RNA elements in transcripts and of chromatin modifications that regulate AS. Application of these tools in plants will provide significant new insights into AS regulation and crosstalk between AS and other layers of gene regulation.

Functional analysis of a large set of BRCA2 exon 7 variants highlights the predictive value of hexamer scores in detecting alterations of exonic splicing regulatory elements

Exonic variants can alter pre-mRNA splicing either by changing splice sites or by modifying splicing regulatory elements. Often these effects are difficult to predict and are only detected by performing RNA analyses. Here, we analyzed, in a minigene assay, 26 variants identified in the exon 7 of BRCA2, a cancer predisposition gene. Our results revealed eight new exon skipping mutations in this exon: one directly altering the 5' splice site and seven affecting potential regulatory elements. This brings the number of splicing regulatory mutations detected in BRCA2 exon 7 to a total of 11, a remarkably high number considering the total number of variants reported in this exon (n = 36), all tested in our minigene assay. We then exploited this large set of splicing data to test the predictive value of splicing regulator hexamers' scores recently established by Ke et al. (). Comparisons of hexamer-based predictions with our experimental data revealed high sensitivity in detecting variants that increased exon skipping, an important feature for prescreening variants before RNA analysis. In conclusion, hexamer scores represent a promising tool for predicting the biological consequences of exonic variants and may have important applications for the interpretation of variants detected by high-throughput sequencing.

Intragenic DNA methylation in transcriptional regulation, normal differentiation and cancer

Ever since the discovery of DNA methylation at cytosine residues, the role of this so called fifth base has been extensively studied and debated. Until recently, the majority of DNA methylation studies focused on the analysis of CpG islands associated to promoter regions. However, with the upcoming possibilities to study DNA methylation in a genome-wide context, this epigenetic mark can now be studied in an unbiased manner. As a result, recent studies have shown that not only promoters but also intragenic and intergenic regions are widely modulated during physiological processes and disease. In particular, it is becoming increasingly clear that DNA methylation in the gene body is not just a passive witness of gene transcription but it seems to be actively involved in multiple gene regulation processes. In this review we discuss the potential role of intragenic DNA methylation in alternative promoter usage, regulation of short and long non-coding RNAs, alternative RNA processing, as well as enhancer activity. Furthermore, we summarize how the intragenic DNA methylome is modified both during normal cell differentiation and neoplastic transformation.

Prediction of clustered RNA-binding protein motif sites in the mammalian genome

Sequence-specific interactions of RNA-binding proteins (RBPs) with their target transcripts are essential for post-transcriptional gene expression regulation in mammals. However, accurate prediction of RBP motif sites has been difficult because many RBPs recognize short and degenerate sequences. Here we describe a hidden Markov model (HMM)-based algorithm mCarts to predict clustered functional RBP-binding sites by effectively integrating the number and spacing of individual motif sites, their accessibility in local RNA secondary structures and cross-species conservation. This algorithm learns and quantifies rules of these features, taking advantage of a large number of in vivo RBP-binding sites obtained from cross-linking and immunoprecipitation data. We applied this algorithm to study two representative RBP families, Nova and Mbnl, which regulate tissue-specific alternative splicing through interacting with clustered YCAY and YGCY elements, respectively, and predicted their binding sites in the mouse transcriptome. Despite the low information content in individual motif elements, our algorithm made specific predictions for successful experimental validation. Analysis of predicted sites also revealed cases of extensive and distal RBP-binding sites important for splicing regulation. This algorithm can be readily applied to other RBPs to infer their RNA-regulatory networks. The software is freely available at http://zhanglab.c2b2.columbia.edu/index.php/MCarts.

Combined computational-experimental analyses of CFTR exon strength uncover predictability of exon-skipping level

With the increased number of identified nucleotide sequence variations in genes, the current challenge is to classify them as disease causing or neutral. These variants of unknown clinical significance can alter multiple processes, from gene transcription to RNA splicing or protein function. Using an approach combining several in silico tools, we identified some exons presenting weaker splicing motifs than other exons in the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene. These exons exhibit higher rates of basal skipping than exons harboring no identifiable weak splicing signals using minigene assays. We then screened 19 described mutations in three different exons, and identified exon-skipping substitutions. These substitutions induced higher skipping levels in exons having one or more weak splicing motifs. Indeed, this level remained under 2% for exons with strong splicing motifs and could reach 40% for exons having at least one weak motif. Further analysis revealed a functional exon splicing enhancer within exon 3 that was associated with the SR protein SF2/ASF and whose disruption induced exon skipping. Exon skipping was confirmed in vivo in two nasal epithelial cell brushing samples. Our approach, which point out exons with some splicing signals weaknesses, will help spot splicing mutations of clinical relevance.

Alternative splicing: a pivotal step between eukaryotic transcription and translation

Alternative splicing was discovered simultaneously with splicing over three decades ago. Since then, an enormous body of evidence has demonstrated the prevalence of alternative splicing in multicellular eukaryotes, its key roles in determining tissue- and species-specific differentiation patterns, the multiple post- and co-transcriptional regulatory mechanisms that control it, and its causal role in hereditary disease and cancer. The emerging evidence places alternative splicing in a central position in the flow of eukaryotic genetic information, between transcription and translation, in that it can respond not only to various signalling pathways that target the splicing machinery but also to transcription factors and chromatin structure.

A biophysical model for identifying splicing regulatory elements and their interactions

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

Identification of an intronic splicing regulatory element involved in auto-regulation of alternative splicing of SCL33 pre-mRNA

In Arabidopsis, pre-mRNAs of serine/arginine-rich (SR) proteins undergo extensive alternative splicing (AS). However, little is known about the cis-elements and trans-acting proteins involved in regulating AS. Using a splicing reporter (GFP-intron-GFP), consisting of the GFP coding sequence interrupted by an alternatively spliced intron of SCL33, we investigated whether cis-elements within this intron are sufficient for AS, and which SR proteins are necessary for regulated AS. Expression of the splicing reporter in protoplasts faithfully produced all splice variants from the intron, suggesting that cis-elements required for AS reside within the intron. To determine which SR proteins are responsible for AS, the splicing pattern of the GFP-intron-GFP reporter was investigated in protoplasts of three single and three double mutants of SR genes. These analyses revealed that SCL33 and a closely related paralog, SCL30a, are functionally redundant in generating specific splice variants from this intron. Furthermore, SCL33 protein bound to a conserved sequence in this intron, indicating auto-regulation of AS. Mutations in four GAAG repeats within the conserved region impaired generation of the same splice variants that are affected in the scl33 scl30a double mutant. In conclusion, we have identified the first intronic cis-element involved in AS of a plant SR gene, and elucidated a mechanism for auto-regulation of AS of this intron.

Interactions of SR45, an SR-like protein, with spliceosomal proteins and an intronic sequence: insights into regulated splicing

SR45 is a serine/arginine-rich (SR)-like protein with two arginine/serine-rich (RS) domains. We have previously shown that SR45 regulates alternative splicing (AS) by differential selection of 5' and 3' splice sites. However, it is unknown how SR45 regulates AS. To gain mechanistic insights into the roles of SR45 in splicing, we screened a yeast two-hybrid library with SR45. This screening resulted in the isolation of two spliceosomal proteins, U1-70K and U2AF(35) b that are known to function in 5' and 3' splice site selection, respectively. This screen not only confirmed our prior observation that U1-70K and SR45 interact, but also helped to identify an additional interacting partner (U2AF(35) ). In vitro and in vivo analyses revealed an interaction of SR45 with both paralogs of U2AF(35) . Furthermore, we show that the RS1 and RS2 domains of SR45, and not the RNA recognition motif (RRM) domain, associate independently with both U2AF(35) proteins. Interaction studies among U2AF(35) paralogs and between U2AF(35) and U1-70K revealed that U2AF(35) can form homo- or heterodimers and that U2AF(35) proteins can associate with U1-70K. Using RNA probes from SR30 intron 10, whose splicing is altered in the sr45 mutant, we show that SR45 and U2AF(35) b bind to different parts of the intron, with a binding site for SR45 in the 5' region and two binding regions, each ending with a known 3' splice site, for U2AF(35) b. These results suggest that SR45 recruits U1snRNP and U2AF to 5' and 3' splice sites, respectively, by interacting with pre-mRNA, U1-70K and U2AF(35) and modulates AS.

An integrated regulatory network reveals pervasive cross-regulation among transcription and splicing factors

Traditionally the gene expression pathway has been regarded as being comprised of independent steps, from RNA transcription to protein translation. To date there is increasing evidence of coupling between the different processes of the pathway, specifically between transcription and splicing. To study the interplay between these processes we derived a transcription-splicing integrated network. The nodes of the network included experimentally verified human proteins belonging to three groups of regulators: transcription factors, splicing factors and kinases. The nodes were wired by instances of predicted transcriptional and alternative splicing regulation. Analysis of the network indicated a pervasive cross-regulation among the nodes; specifically, splicing factors are significantly more connected by alternative splicing regulatory edges relative to the two other subgroups, while transcription factors are more extensively controlled by transcriptional regulation. Furthermore, we found that splicing factors are the most regulated of the three regulatory groups and are subject to extensive combinatorial control by alternative splicing and transcriptional regulation. Consistent with the network results, our bioinformatics analyses showed that the subgroup of kinases have the highest density of predicted phosphorylation sites. Overall, our systematic study reveals that an organizing principle in the logic of integrated networks favor the regulation of regulatory proteins by the specific regulation they conduct. Based on these results, we propose a new regulatory paradigm postulating that gene expression regulation of the master regulators in the cell is predominantly achieved by cross-regulation.

Intronic features that determine the selection of the 3' splice site

Most eukaryotic primary transcripts include segments, or introns, that will be accurately removed during RNA biogenesis. This process, known as pre-messenger RNA splicing, is catalyzed by the spliceosome, accurately selecting a set of intronic marks from others apparently equivalent. This identification is critical, as incorrectly spliced RNAs can be toxic for the organism. One of these marks, the dinucleotide AG, signals the intronic 3' end, or 3' splice site (ss). In this review we will focus on those intronic features that have an impact on 3' ss selection. These include the location and type of neighboring sequences, and their distance to the 3' end. We will see that their interplay is needed to select the right intronic end, and that this can be modulated by additional intronic elements that contribute to alternative splicing, whereby diverse RNAs can be generated from identical precursors. This complexity, still poorly understood, is fundamental for the accuracy of gene expression. In addition, a clear knowledge of 3' ss selection is needed to fully decipher the coding potential of genomes.

VERSE: a varying effect regression for splicing elements discovery

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

Deciphering the plant splicing code: experimental and computational approaches for predicting alternative splicing and splicing regulatory elements

Extensive alternative splicing (AS) of precursor mRNAs (pre-mRNAs) in multicellular eukaryotes increases the protein-coding capacity of a genome and allows novel ways to regulate gene expression. In flowering plants, up to 48% of intron-containing genes exhibit AS. However, the full extent of AS in plants is not yet known, as only a few high-throughput RNA-Seq studies have been performed. As the cost of obtaining RNA-Seq reads continues to fall, it is anticipated that huge amounts of plant sequence data will accumulate and help in obtaining a more complete picture of AS in plants. Although it is not an onerous task to obtain hundreds of millions of reads using high-throughput sequencing technologies, computational tools to accurately predict and visualize AS are still being developed and refined. This review will discuss the tools to predict and visualize transcriptome-wide AS in plants using short-reads and highlight their limitations. Comparative studies of AS events between plants and animals have revealed that there are major differences in the most prevalent types of AS events, suggesting that plants and animals differ in the way they recognize exons and introns. Extensive studies have been performed in animals to identify cis-elements involved in regulating AS, especially in exon skipping. However, few such studies have been carried out in plants. Here, we review the current state of research on splicing regulatory elements (SREs) and briefly discuss emerging experimental and computational tools to identify cis-elements involved in regulation of AS in plants. The availability of curated alternative splice forms in plants makes it possible to use computational tools to predict SREs involved in AS regulation, which can then be verified experimentally. Such studies will permit identification of plant-specific features involved in AS regulation and contribute to deciphering the splicing code in plants.

Bioinformatics and mutations leading to exon skipping

Our knowledge about human genes and the consequences of mutations leading to human genetic diseases has drastically improved over the last few years. It has been recognized that many mutations are indeed pathogenic because they impact the mRNA rather than the protein itself. With our better understanding of the very complex mechanism of splicing, various bioinformatics tools have been developed. They are now frequently used not only to search for sequence motifs corresponding to splicing signals (splice sites, branch points, ESE, and ESS) but also to predict the impact of mutations on these signals. We now need to address the impact of mutations that affect the splicing process, as their consequences could vary from the activation of cryptic signals to the skipping of one or multiple exons. Despite the major developments of the bioinformatics field coupled to experimental data generated on splicing, it is today still not possible to efficiently predict the consequences of mutations impacting splicing signals, especially to predict if they will lead to exon skipping or to cryptic splice site activation.

Alternative Splicing of Fibroblast Growth Factor Receptor IgIII Loops in Cancer

Alternative splicing of the IgIII loop of fibroblast growth factor receptors (FGFRs) 1–3 produces b- and c-variants of the receptors with distinctly different biological impact based on their distinct ligand-binding spectrum. Tissue-specific expression of these splice variants regulates interactions in embryonic development, tissue maintenance and repair, and cancer. Alterations in FGFR2 splicing are involved in epithelial mesenchymal transition that produces invasive, metastatic features during tumor progression. Recent research has elucidated regulatory factors that determine the splice choice both on the level of exogenous signaling events and on the RNA-protein interaction level. Moreover, methodology has been developed that will enable the in depth analysis of splicing events during tumorigenesis and provide further insight on the role of FGFR 1–3 IIIb and IIIc in the pathophysiology of various malignancies. This paper aims to summarize expression patterns in various tumor types and outlines possibilities for further analysis and application.

Identification of cis-regulatory sequence variations in individual genome sequences

Functional contributions of cis-regulatory sequence variations to human genetic disease are numerous. For instance, disrupting variations in a HNF4A transcription factor binding site upstream of the Factor IX gene contributes causally to hemophilia B Leyden. Although clinical genome sequence analysis currently focuses on the identification of protein-altering variation, the impact of cis-regulatory mutations can be similarly strong. New technologies are now enabling genome sequencing beyond exomes, revealing variation across the non-coding 98% of the genome responsible for developmental and physiological patterns of gene activity. The capacity to identify causal regulatory mutations is improving, but predicting functional changes in regulatory DNA sequences remains a great challenge. Here we explore the existing methods and software for prediction of functional variation situated in the cis-regulatory sequences governing gene transcription and RNA processing.

Efficient internal exon recognition depends on near equal contributions from the 3' and 5' splice sites

Pre-mRNA splicing is carried out by the spliceosome, which identifies exons and removes intervening introns. In vertebrates, most splice sites are initially recognized by the spliceosome across the exon, because most exons are small and surrounded by large introns. This gene architecture predicts that efficient exon recognition depends largely on the strength of the flanking 3' and 5' splice sites. However, it is unknown if the 3' or the 5' splice site dominates the exon recognition process. Here, we test the 3' and 5' splice site contributions towards efficient exon recognition by systematically replacing the splice sites of an internal exon with sequences of different splice site strengths. We show that the presence of an optimal splice site does not guarantee exon inclusion and that the best predictor for exon recognition is the sum of both splice site scores. Using a genome-wide approach, we demonstrate that the combined 3' and 5' splice site strengths of internal exons provide a much more significant separator between constitutive and alternative exons than either the 3' or the 5' splice site strength alone.

Quantitative evaluation of all hexamers as exonic splicing elements

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

Design principles for bifunctional targeted oligonucleotide enhancers of splicing

Controlling the patterns of splicing of specific genes is an important goal in the development of new therapies. We have shown that the splicing of a refractory exon, SMN2 exon 7, could be increased in fibroblasts derived from patients with spinal muscular atrophy by using bifunctional targeted oligonucleotide enhancers of splicing (TOES) oligonucleotides that anneal to the exon and contain a ‘tail’ of enhancer sequences that recruit activating proteins. We show here that there are striking agreements between the effects of oligonucleotides on splicing in vitro and on both splicing and SMN2 protein expression in patient-derived fibroblasts, indicating that the effects on splicing are the major determinant of success. Increased exon inclusion depends on the number, sequence and chemistry of the motifs that bind the activator protein SRSF1, but it is not improved by increasing the strength of annealing to the target site. The optimal oligonucleotide increases protein levels in transfected fibroblasts by a mean value of 2.6-fold (maximum 4.6-fold), and after two rounds of transfection the effect lasted for a month. Oligonucleotides targeted to the upstream exon (exon 6 in SMN) are also effective. We conclude that TOES oligonucleotides are highly effective reagents for restoring the splicing of refractory exons and can act across long introns.

Context-dependent splicing regulation: exon definition, co-occurring motif pairs and tissue specificity

Splicing is a crucial process in gene expression in higher organisms because: 1) most vertebrate genes contain introns; and 2) alternative splicing is primarily responsible for increasing proteomic complexity and functional diversity. Intron definition, the coordination across an intron, is a mandatory step in the splicing process. However, exon definition, the coordination across an exon, is also thought to be required for the splicing of most vertebrate exons. Recent investigations of exon definition complexes provide insights into splicing dynamics. That splicing regulators act in a context-dependent mode is supported by a large collection of evidence. Splicing contexts generally can be classified as cis-element and trans-element contexts. A widespread cis-element context is defined by co-occurring motif pairs to which splicing regulatory factors bind to direct specific molecular interactions. Splicing regulation is also coordinated by trans-element contexts as exemplified by tissue specific splicing, where alternative exons can be coordinately regulated by a few splicing factors, the expression and/or activity of which are concertedly higher or lower in the corresponding tissues.

Analysis of the single-nucleotide polymorphism in the 5'UTR and part of intron I of the sheep MSTN gene

The myostatin (MSTN) gene region encompassing the 5'UTR and part of intron I was sequenced in animals of two herds of Latvian Darkhead sheep to extend data on the ovine MSTN gene polymorphism and to provide information useful for local breed conservation. Two and four polymorphic loci were revealed in the 5'UTR and intron I. Four and five local haplotypes were constructed, respectively. The genotyping data obtained and that previously reported for the same genomic region were combined in one dataset for the haplotype analysis. Recombination events were detected between loci (c.-40, c.-37) in the 5'UTR and (c.373+18, c.373+101) and (c.373+101, c.373+241) in intron I. Single-nucleotide polymorphisms at c.373+249 and c.373+323 appear to be involved in the strong linkage (p < 0.01). Linkage blocks (c.373+241, c.373+243) and (c.373+241, c.373+259) were revealed at nominal (p < 0.05) level of probability. Haplotype-specific patterns of the transcription factor binding sites predicted in silico were constructed to evaluate a putative functional significance of the particular alleles and haplotypes. A nucleotide at c.373+18 was shown to influence the pre-mRNA secondary structure. DNA curvature predicted in silico for allele c.373+101C was proven experimentally. A possible impact of the particular polymorphisms on the transcription and/or splicing efficiency is discussed.

Using bioinformatics to predict the functional impact of SNVs

MOTIVATION

The past decade has seen the introduction of fast and relatively inexpensive methods to detect genetic variation across the genome and exponential growth in the number of known single nucleotide variants (SNVs). There is increasing interest in bioinformatics approaches to identify variants that are functionally important from millions of candidate variants. Here, we describe the essential components of bioinformatics tools that predict functional SNVs.

RESULTS

Bioinformatics tools have great potential to identify functional SNVs, but the black box nature of many tools can be a pitfall for researchers. Understanding the underlying methods, assumptions and biases of these tools is essential to their intelligent application.

Epigenetics in alternative pre-mRNA splicing

Alternative splicing plays critical roles in differentiation, development, and disease and is a major source for protein diversity in higher eukaryotes. Analysis of alternative splicing regulation has traditionally focused on RNA sequence elements and their associated splicing factors, but recent provocative studies point to a key function of chromatin structure and histone modifications in alternative splicing regulation. These insights suggest that epigenetic regulation determines not only what parts of the genome are expressed but also how they are spliced.

Systems analysis of alternative splicing and its regulation

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

Developments in RNA splicing and disease

Pre-mRNA processing, including 5'-end capping, splicing, editing, and polyadenylation, consists of a series of orchestrated and primarily cotranscriptional steps that ensure both the high fidelity and extreme diversity characteristic of eukaryotic gene expression. Alternative splicing and editing allow relatively small genomes to encode vast proteomic arrays while alternative 3'-end formation enables variations in mRNA localization, translation, and stability. Of course, this mechanistic complexity comes at a high price. Mutations in the myriad of RNA sequence elements that regulate mRNA biogenesis, as well as the trans-acting factors that act upon these sequences, underlie a number of human diseases. In this review, we focus on one of these key RNA processing steps, splicing, to highlight recent studies that describe both conventional and novel pathogenic mechanisms that underlie muscle and neurological diseases.

Quantification of CFTR transcripts

Quantification and analysis of CFTR transcripts is of crucial importance not only for cystic fibrosis (CF) diagnosis and prognosis, but also in evaluating the efficiency of various therapeutic approaches to CF, including gene therapy. Reverse transcription (RT) followed by quantitative polymerase chain reaction (qPCR) is at present the most sensitive method for transcript abundance measurement. Classical RNA-based methods require significant expression levels in target samples for appropriate analysis, thus PCR-based methods have evolved towards reliable quantification. In this chapter we describe and discuss several protocols for the quantitative analysis of CFTR transcripts, including those variants that result from alternative splicing.

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

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

A recombination hotspot leads to sequence variability within a novel gene (AK005651) and contributes to type 1 diabetes susceptibility

More than 25 loci have been linked to type 1 diabetes (T1D) in the nonobese diabetic (NOD) mouse, but identification of the underlying genes remains challenging. We describe here the positional cloning of a T1D susceptibility locus, Idd11, located on mouse chromosome 4. Sequence analysis of a series of congenic NOD mouse strains over a critical 6.9-kb interval in these mice and in 25 inbred strains identified several haplotypes, including a unique NOD haplotype, associated with varying levels of T1D susceptibility. Haplotype diversity within this interval between congenic NOD mouse strains was due to a recombination hotspot that generated four crossover breakpoints, including one with a complex conversion tract. The Idd11 haplotype and recombination hotspot are located within a predicted gene of unknown function, which exhibits decreased expression in relevant tissues of NOD mice. Notably, it was the recombination hotspot that aided our mapping of Idd11 and confirms that recombination hotspots can create genetic variation affecting a common polygenic disease. This finding has implications for human genetic association studies, which may be affected by the approximately 33,000 estimated hotspots in the genome.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Subtle discrepancies of SF2/ASF ESE sequence motif among human tissues: A computational approach

The intron removal during the pre-mRNA splicing in higher eukaryotes requires the accurate identification of the two splice sites at the ends of the exons, or exon definition. However, the consensus sequences at the splice sites provide insufficient information to distinguish true splice sites from the large number of the false ones that populate the primary transcripts. Additional information is provided by cis-acting regulatory sequences that serve to enhance or repress splicing, and that may be exonic or intronic in nature: the splicing enhancers and the splicing silencers, respectively. In this study, we tested by computational and statistical approaches if the exonic splicing enhancer motif binding to the SF2/ASF SR protein is conserved among several groups of human genes. The results showed that the SF2/ASF ESE consensus was conserved between genes within the same chromosome, within different chromosomes and between different levels of muscular cells differentiation. However, this motif displays subtle variations within the consensus sequence between genes expressed in different tissues. These results can emphasize the presence of different translational isoforms of the SFRS1 gene encoding for the SF2/ASF, or different post-translational protein maturations in different tissues. This tissular discrepancy can also account for the alternative splicing of several genes between tissues.

An antisense microwalk reveals critical role of an intronic position linked to a unique long-distance interaction in pre-mRNA splicing

Here we report a novel finding of an antisense oligonucleotide (ASO) microwalk in which we examined the position-specific role of intronic residues downstream from the 5' splice site (5' ss) of SMN2 exon 7, skipping of which is associated with spinal muscular atrophy (SMA), a leading genetic cause of infant mortality. Our results revealed the inhibitory role of a cytosine residue at the 10th intronic position ((10)C), which is neither conserved nor associated with any known splicing motif. Significance of (10)C emerged from the splicing pattern of SMN2 exon 7 in presence of a 14-mer ASO (L14) that sequestered two adjacent hnRNP A1 motifs downstream from (10)C and yet promoted SMN2 exon 7 skipping. Another 14-mer ASO (F14) that sequestered both, (10)C and adjacent hnRNP A1 motifs, led to a strong stimulation of SMN2 exon 7 inclusion. The inhibitory role of (10)C was found to be tightly linked to its unpaired status and specific positioning immediately upstream of a RNA:RNA helix formed between the targeting ASO and its intronic target. Employing a heterologous context as well as changed contexts of SMN2 intron 7, we show that the inhibitory effect of unpaired (10)C is dependent upon a long-distance interaction involving downstream intronic sequences. Our report furnishes one of the rare examples in which an ASO-based approach could be applied to unravel the critical role of an intronic position that may not belong to a linear motif and yet play significant role through long-distance interactions.

SFmap: a web server for motif analysis and prediction of splicing factor binding sites

Alternative splicing (AS) is a post-transcriptional process considered to be responsible for the huge diversity of proteins in higher eukaryotes. AS events are regulated by different splicing factors (SFs) that bind to sequence elements on the RNA. SFmap is a web server for predicting putative SF binding sites in genomic data (http://sfmap.technion.ac.il). SFmap implements the COS(WR) algorithm, which computes similarity scores for a given regulatory motif based on information derived from its sequence environment and its evolutionary conservation. Input for SFmap is a human genomic sequence or a list of sequences in FASTA format that can either be uploaded from a file or pasted into a window. SFmap searches within a given sequence for significant hits of binding motifs that are either stored in our database or defined by the user. SFmap results are provided both as a text file and as a graphical web interface.

A comprehensive survey of human polymorphisms at conserved splice dinucleotides and its evolutionary relationship with alternative splicing

BACKGROUND

Alternative splicing (AS) is a key molecular process that endows biological functions with diversity and complexity. Generally, functional redundancy leads to the generation of new functions through relaxation of selective pressure in evolution, as exemplified by duplicated genes. It is also known that alternatively spliced exons (ASEs) are subject to relaxed selective pressure. Within consensus sequences at the splice junctions, the most conserved sites are dinucleotides at both ends of introns (splice dinucleotides). However, a small number of single nucleotide polymorphisms (SNPs) occur at splice dinucleotides. An intriguing question relating to the evolution of AS diversity is whether mutations at splice dinucleotides are maintained as polymorphisms and produce diversity in splice patterns within the human population. We therefore surveyed validated SNPs in the database dbSNP located at splice dinucleotides of all human genes that are defined by the H-Invitational Database.

RESULTS

We found 212 validated SNPs at splice dinucleotides (sdSNPs); these were confirmed to be consistent with the GT-AG rule at either allele. Moreover, 53 of them were observed to neighbor ASEs (AE dinucleotides). No significant differences were observed between sdSNPs at AE dinucleotides and those at constitutive exons (CE dinucleotides) in SNP properties including average heterozygosity, SNP density, ratio of predicted alleles consistent with the GT-AG rule, and scores of splice sites formed with the predicted allele. We also found that the proportion of non-conserved exons was higher for exons with sdSNPs than for other exons.

CONCLUSIONS

sdSNPs are found at CE dinucleotides in addition to those at AE dinucleotides, suggesting two possibilities. First, sdSNPs at CE dinucleotides may be robust against sdSNPs because of unknown mechanisms. Second, similar to sdSNPs at AE dinucleotides, those at CE dinucleotides cause differences in AS patterns because of the arbitrariness in the classification of exons into alternative and constitutive type that varies according to the dataset. Taking into account the absence of differences in sdSNP properties between those at AE and CE dinucleotides, the increased proportion of non-conserved exons found in exons flanked by sdSNPs suggests the hypothesis that sdSNPs are maintained at the splice dinucleotides of newly generated exons at which negative selection pressure is relaxed.

Novel sequence variations in LAMA2 andSGCG genes modulating cis-acting regulatory elements and RNA secondary structure

In this study, we detected new sequence variations in LAMA2 and SGCG genes in 5 ethnic populations, and analysed their effect on enhancer composition and mRNA structure. PCR amplification and DNA sequencing were performed and followed by bioinformatics analyses using ESEfinder as well as MFOLD software. We found 3 novel sequence variations in the LAMA2 (c.3174+22_23insAT and c.6085 +12delA) and SGCG (c. ^* 102A/C) genes. These variations were present in 210 tested healthy controls from Tunisian, Moroccan, Algerian, Lebanese and French populations suggesting that they represent novel polymorphisms within LAMA2 and SGCG genes sequences. ESEfinder showed that the c. ^* 102A/C substitution created a new exon splicing enhancer in the 3'UTR of SGCG genes, whereas the c.6085 +12delA deletion was situated in the base pairing region between LAMA2 mRNA and the U1snRNA spliceosomal components. The RNA structure analyses showed that both variations modulated RNA secondary structure. Our results are suggestive of correlations between mRNA folding and the recruitment of spliceosomal components mediating splicing, including SR proteins. The contribution of common sequence variations to mRNA structural and functional diversity will contribute to a better study of gene expression.

Alternative splicing: role of pseudoexons in human disease and potential therapeutic strategies

What makes a nucleotide sequence an exon (or an intron) is a question that still lacks a satisfactory answer. Indeed, most eukaryotic genes are full of sequences that look like perfect exons, but which are nonetheless ignored by the splicing machinery (hence the name 'pseudoexons'). The existence of these pseudoexons has been known since the earliest days of splicing research, but until recently the tendency has been to view them as an interesting, but rather rare, curiosity. In recent years, however, the importance of pseudoexons in regulating splicing processes has been steadily revalued. Even more importantly, clinically oriented screening studies that search for splicing mutations are beginning to uncover a situation where aberrant pseudoexon inclusion as a cause of human disease is more frequent than previously thought. Here we aim to provide a review of the mechanisms that lead to pseudoexon activation in human genes and how the various cis- and trans-acting cellular factors regulate their inclusion. Moreover, we list the potential therapeutic approaches that are being tested with the aim of inhibiting their inclusion in the final mRNA molecules.

Functional properties and evolutionary splicing constraints on a composite exonic regulatory element of splicing in CFTR exon 12

In general, splicing regulatory elements are defined as Enhancers or Silencers depending on their positive or negative effect upon exon inclusion. Often, these sequences are usually present separate from each other in exonic/intronic sequences. The Composite Exonic Splicing Regulatory Elements (CERES) represent an extreme physical overlap of enhancer/silencer activity. As a result, when CERES elements are mutated the consequences on the splicing process are difficult to predict. Here, we show that the functional activity of the CERES2 sequence in CFTR exon 12 is regulated by the binding, in very close proximity to each other, of several SR and hnRNP proteins. Moreover, our results show that practically the entire exon 12 sequence context participate in its definition. The consequences of this situation can be observed at the evolutionary level by comparing changes in conservation of different splicing elements in different species. In conclusion, our study highlights how it is increasingly difficult to define many exonic sequences by simply breaking them down in isolated enhancer/silencer or even neutral elements. The real picture is close to one of continuous competition between positive and negative factors where affinity for the target sequences and other dynamic factors decide the inclusion or exclusion of the exon.

Divergence of exonic splicing elements after gene duplication and the impact on gene structures

An analysis of human exonic splicing elements in duplicated genes reveals their important role in the generation of new gene structures.

Background

The origin of new genes and their contribution to functional novelty has been the subject of considerable interest. There has been much progress in understanding the mechanisms by which new genes originate. Here we examine a novel way that new gene structures could originate, namely through the evolution of new alternative splicing isoforms after gene duplication.

Results

We studied the divergence of exonic splicing enhancers and silencers after gene duplication and the contributions of such divergence to the generation of new splicing isoforms. We found that exonic splicing enhancers and exonic splicing silencers diverge especially fast shortly after gene duplication. About 10% and 5% of paralogous exons undergo significantly asymmetric evolution of exonic splicing enhancers and silencers, respectively. When compared to pre-duplication ancestors, we found that there is a significant overall loss of exonic splicing enhancers and the magnitude increases with duplication age. Detailed examination reveals net gains and losses of exonic splicing enhancers and silencers in different copies and paralog clusters after gene duplication. Furthermore, we found that exonic splicing enhancer and silencer changes are mainly caused by synonymous mutations, though nonsynonymous changes also contribute. Finally, we found that exonic splicing enhancer and silencer divergence results in exon splicing state transitions (from constitutive to alternative or vice versa), and that the proportion of paralogous exon pairs with different splicing states also increases over time, consistent with previous predictions.

Conclusions

Our results suggest that exonic splicing enhancer and silencer changes after gene duplication have important roles in alternative splicing divergence and that these changes contribute to the generation of new gene structures.

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.

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

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

A short antisense oligonucleotide masking a unique intronic motif prevents skipping of a critical exon in spinal muscular atrophy

Spinal muscular atrophy (SMA) is the leading genetic cause of infant mortality. Most SMA cases are associated with the low levels of SMN owing to deletion of Survival Motor Neuron 1 (SMN1). SMN2, a nearly identical copy of SMN1, fails to compensate for the loss of SMN1 due to predominant skipping of exon 7. Hence, correction of aberrant splicing of SMN2 exon 7 holds the potential for cure of SMA. Here we report an 8-mer antisense oligonucleotide (ASO) to have a profound stimulatory response on correction of aberrant splicing of SMN2 exon 7 by binding to a unique GC-rich sequence located within intron 7 of SMN2. We confirm that the splicing-switching ability of this short ASO comes with a high degree of specificity and reduced off-target effect compared to larger ASOs targeting the same sequence. We further demonstrate that a single low nanomolar dose of this 8-mer ASO substantially increases the levels of SMN and a host of factors including Gemin 2, Gemin 8, ZPR1, hnRNP Q and Tra2-beta1 known to be down-regulated in SMA. Our findings underscore the advantages and unmatched potential of very short ASOs in splicing modulation in vivo.

An analysis of the human chemokine CXC receptor 4 gene

In this article we analyze some of the structural characteristics of the coding section and the intron of the human chemokine CXC receptor 4 (a 7-transmembrane receptor) pre-mRNA. In the coding sequence the frequencies of the individual nucleotides do not depart significantly from 0.25, while in the intron the frequencies of the As and Gs are significantly lower and higher, respectively, than expected from a random distribution. Analysis of the pattern of association of nucleotides into triplets or couples shows that some triplets or couples occur with frequencies significantly higher or lower than expected when assuming a random association of nucleotides. In particular, in the intron combinations of the same nucleotide are over-represented. 7-or-more nucleotide repeats occur in both the coding section and the intron with frequencies which exceed the confidence limits for a random distribution. For the coding sequence this is possibly explained by the alternans of relatively similar hydrophobic-coding sections and relatively similar intervening intracellular and extracellular hydrophilic-coding sections. 7-or-more nucleotide repeats in reverse order and in reverse/complemented order occur in the intron, but not in the coding section, with frequencies which significantly exceed a random distribution. The numerous intronic repeats in reverse/complemented order may be of relevance for the secondary structure of the intron and might be one important element of the integrated splicing code.