Intronic alternative splicing regulators identified by comparative genomics in nematodes

Many alternative splicing events are regulated by pentameric and hexameric intronic sequences that serve as binding sites for splicing regulatory factors. We hypothesized that intronic elements that regulate alternative splicing are under selective pressure for evolutionary conservation. Using a Wobble Aware Bulk Aligner genomic alignment of Caenorhabditis elegans and Caenorhabditis briggsae, we identified 147 alternatively spliced cassette exons that exhibit short regions of high nucleotide conservation in the introns flanking the alternative exon. In vivo experiments on the alternatively spliced let-2 gene confirm that these conserved regions can be important for alternative splicing regulation. Conserved intronic element sequences were collected into a dataset and the occurrence of each pentamer and hexamer motif was counted. We compared the frequency of pentamers and hexamers in the conserved intronic elements to a dataset of all C. elegans intron sequences in order to identify short intronic motifs that are more likely to be associated with alternative splicing. High-scoring motifs were examined for upstream or downstream preferences in introns surrounding alternative exons. Many of the high- scoring nematode pentamer and hexamer motifs correspond to known mammalian splicing regulatory sequences, such as (T)GCATG, indicating that the mechanism of alternative splicing regulation is well conserved in metazoans. A comparison of the analysis of the conserved intronic elements, and analysis of the entire introns flanking these same exons, reveals that focusing on intronic conservation can increase the sensitivity of detecting putative splicing regulatory motifs. This approach also identified novel sequences whose role in splicing is under investigation and has allowed us to take a step forward in defining a catalog of splicing regulatory elements for an organism. In vivo experiments confirm that one novel high-scoring sequence from our analysis, (T)CTATC, is important for alternative splicing regulation of the unc-52 gene.     

Genes encoding human killer-cell Ig-like receptors with D1 and D2 extracellular domains all contain untranslated pseudoexons encoding a third Ig-like domain

Human killer-cell immunoglobulin-like receptors (KIR) show three types of organization of their extracellular domains: D0-D1-D2 in KIR3D, D1-D2 in the majority of KIR2D, and D0-D2 in KIR2DL4 and the novel KIR2DL5. The gene for a KIR2DL3 variant, which has a D1-D2 structure, has been shown previously to have a nonexpressed region (pseudoexon 3) that is paralogous to the exon encoding the D0 domain of other KIR. This pseudoexon is not expressed because it is skipped during splicing of pre-mRNA. In this study, we demonstrate that all eight genes encoding human KIR with D1-D2 configuration (KIR2DL1-KIR2DL3, KIR2DS1-KIR2DS5) have similarly untranslated pseudoexons. Whereas the pseudoexons of four of these KIR genes bear nonsense mutations and/or altered splicing sites, the pseudoexons in the other four KIR genes have no major structural abnormalities, indicating that other mechanisms are responsible for inactivation of their exons 3. A comparison of the sequences on pseudoexons 3 with the paralogous expressed exons suggests that an exonic splicing enhancer may be necessary for the expression of exon 3 in KIR genes.     

Dual Role of G-runs and hnRNP F in the Regulation of a Mutation-Activated Pseudoexon in the Fibrinogen Gamma-Chain Transcript

Most pathological pseudoexon inclusion events originate from single activating mutations, suggesting that many intronic sequences are on the verge of becoming exons. However, the precise mechanisms controlling pseudoexon definition are still largely unexplored. Here, we investigated the cis-acting elements and trans-acting regulatory factors contributing to the regulation of a previously described fibrinogen gamma-chain (FGG) pseudoexon, which is activated by a deep-intronic mutation (IVS6-320A>T). This pseudoexon contains several G-run elements, which may be bound by heterogeneous nuclear ribonucleoproteins (hnRNPs) F and H. To explore the effect of these proteins on FGG pseudoexon inclusion, both silencing and overexpression experiments were performed in eukaryotic cells. While hnRNP H did not significantly affect pseudoexon splicing, hnRNP F promoted pseudoexon inclusion, indicating that these two proteins have only partially redundant functions. To verify the binding of hnRNP F and the possible involvement of other trans-acting splicing modulators, pulldown experiments were performed on the region of the pseudoexon characterized by both a G-run and enrichment for exonic splicing enhancers. This 25-bp-long region strongly binds hnRNP F/H and weakly interacts with Serine/Arginine-rich protein 40, which however was demonstrated to be dispensable for FGG pseudoexon inclusion in overexpression experiments. Deletion analysis, besides confirming the splicing-promoting role of the G-run within this 25-bp region, demonstrated that two additional hnRNP F binding sites might instead function as silencer elements. Taken together, our results indicate a major role of hnRNP F in regulating FGG pseudoexon inclusion, and strengthen the notion that G-runs may function either as splicing enhancers or silencers of the same exon.     

The intronic splicing code: multiple factors involved in ATM pseudoexon definition

Abundance of pseudo splice sites in introns can potentially give rise to innumerable pseudoexons, outnumbering the real ones. Nonetheless, these are efficiently ignored by the splicing machinery, a process yet to be understood completely. Although numerous 5′ splice site‐like sequences functioning as splicing silencers have been found to be enriched in predicted human pseudoexons, the lack of active pseudoexons pose a fundamental challenge to how these U1snRNP‐binding sites function in splicing inhibition. Here, we address this issue by focusing on a previously described pathological ATM pseudoexon whose inhibition is mediated by U1snRNP binding at intronic splicing processing element (ISPE), composed of a consensus donor splice site. Spliceosomal complex assembly demonstrates inefficient A complex formation when ISPE is intact, implying U1snRNP‐mediated unproductive U2snRNP recruitment. Furthermore, interaction of SF2/ASF with its motif seems to be dependent on RNA structure and U1snRNP interaction. Our results suggest a complex combinatorial interplay of RNA structure and trans‐acting factors in determining the splicing outcome and contribute to understanding the intronic splicing code for the ATM pseudoexon.     

Pre-spliceosomal binding of U1 small nuclear ribonucleoprotein (RNP) and heterogenous nuclear RNP E1 is associated with suppression of a growth hormone receptor pseudoexon

Pseudoexons occur frequently in the human genome. This paper characterizes a pseudoexon in the GH receptor gene. Inappropriate activation of this pseudoexon causes Laron syndrome. Using in vitro splicing assays, pseudoexon silencing was shown to require a combination of a weak 5' pseudosplice-site and splicing silencing elements within the pseudoexon. Immunoprecipitation experiments showed that specific binding of heterogenous nuclear ribonucleoprotein E1 (hnRNP E1) and U1 small nuclear ribonucleoprotein (snRNP) in the pre-spliceosomal complex was associated with silencing of pseudoexon splicing. The possible role of hnRNP E1 was further supported by RNA interference experiments in cultured cells. Immunoprecipitation experiments with three other pseudoexons suggested that pre-spliceosomal binding of U1 snRNP is a potential general mechanism of suppression of pseudoexons.     

RNA structure is a key regulatory element in pathological ATM and CFTR pseudoexon inclusion events

Genomic variations deep in the intronic regions of pre-mRNA molecules are increasingly reported to affect splicing events. However, there is no general explanation why apparently similar variations may have either no effect on splicing or cause significant splicing alterations. In this work we have examined the structural architecture of pseudoexons previously described in ATM and CFTR patients. The ATM case derives from the deletion of a repressor element and is characterized by an aberrant 5′ss selection despite the presence of better alternatives. The CFTR pseudoexon instead derives from the creation of a new 5′ss that is used while a nearby pre-existing donor-like sequence is never selected. Our results indicate that RNA structure is a major splicing regulatory factor in both cases. Furthermore, manipulation of the original RNA structures can lead to pseudoexon inclusion following the exposure of unused 5′ss already present in their wild-type intronic sequences and prevented to be recognized because of their location in RNA stem structures. Our data show that intrinsic structural features of introns must be taken into account to understand the mechanism of pseudoexon activation in genetic diseases. Our observations may help to improve diagnostics prediction programmes and eventual therapeutic targeting.     

Combinatorial control of a neuron-specific exon.

The mouse c-src gene contains a short neuron-specific exon, N1. N1 exon splicing is partly controlled by an intronic splicing enhancer sequence that activates splicing of a heterologous reporter exon in both neural and nonneural cells. Here we attempt to dissect all of the regulatory elements controlling the N1 exon and examine how these multiple elements work in combination. We show that the 3' splice site sequence upstream of exon N1 represses the activation of splicing by the downstream intronic enhancer. This repression is stronger in nonneural cells and these two regulatory sequences combine to make a reporter exon highly cell-type specific. Substitution of the 3' splice site of this test exon with sites from other exons indicates that activation by the enhancer is very dependent on the nature of the upstream 3' splice site. In addition, we identify a previously uncharacterized purine-rich sequence within exon N1 that cooperates with the downstream intronic enhancer to increase exon inclusion. Finally, different regulatory elements were tested in multiple cell lines of both neuronal and nonneuronal origin. The individual splicing regulatory sequences from the src gene vary widely in their activity between different cell lines. These results demonstrate how a simple cassette exon is controlled by a variety of regulatory elements that only in combination will produce the correct tissue specificity of 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 (http://rulai.cshl.edu/tools/ESE/). 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.     

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.     

Systematic identification and analysis of exonic splicing silencers

Exonic splicing silencers (ESSs) are cis-regulatory elements that inhibit the use of adjacent splice sites, often contributing to alternative splicing (AS). To systematically identify ESSs, an in vivo splicing reporter system was developed to screen a library of random decanucleotides. The screen yielded 141 ESS decamers, 133 of which were unique. The silencer activity of over a dozen of these sequences was also confirmed in a heterologous exon/intron context and in a second cell type. Of the unique ESS decamers, most could be clustered into groups to yield seven putative ESS motifs, some resembling known motifs bound by hnRNPs H and A1. Potential roles of ESSs in constitutive splicing were explored using an algorithm, ExonScan, which simulates splicing based on known or putative splicing-related motifs. ExonScan and related bioinformatic analyses suggest that these ESS motifs play important roles in suppression of pseudoexons, in splice site definition, and in AS.     

Exon creation and establishment in human genes

Background

A large proportion of species-specific exons are alternatively spliced. In primates, Alu elements play a crucial role in the process of exon creation but many new exons have appeared through other mechanisms. Despite many recent studies, it is still unclear which are the splicing regulatory requirements for de novo exonization and how splicing regulation changes throughout an exon's lifespan.

Results

Using comparative genomics, we have defined sets of exons with different evolutionary ages. Younger exons have weaker splice-sites and lower absolute values for the relative abundance of putative splicing regulators between exonic and adjacent intronic regions, indicating a less consolidated splicing regulation. This relative abundance is shown to increase with exon age, leading to higher exon inclusion. We show that this local difference in the density of regulators might be of biological significance, as it outperforms other measures in real exon versus pseudo-exon classification. We apply this new measure to the specific case of the exonization of anti-sense Alu elements and show that they are characterized by a general lack of exonic splicing silencers.

Conclusions

Our results suggest that specific sequence environments are required for exonization and that these can change with time. We propose a model of exon creation and establishment in human genes, in which splicing decisions depend on the relative local abundance of regulatory motifs. Using this model, we provide further explanation as to why Alu elements serve as a major substrate for exon creation in primates. Finally, we discuss the benefits of integrating such information in gene prediction.     

A novel intronic cis element, ISE/ISS-3, regulates rat fibroblast growth factor receptor 2 splicing through activation of an upstream exon and repression of a downstream exon containing a noncanonical branch point sequence

Mutually exclusive splicing of fibroblast growth factor receptor 2 (FGFR2) exons IIIb and IIIc yields two receptor isoforms, FGFR2-IIIb and -IIIc, with distinctly different ligand binding properties. Several RNA cis elements in the intron (intron 8) separating these exons have been described that are required for splicing regulation. Using a heterologous splicing reporter, we have identified a new regulatory element in this intron that confers cell-type-specific inclusion of an unrelated exon that mirrors its ability to promote cell-type-specific inclusion of exon IIIb. This element promoted inclusion of exon IIIb while at the same time silencing exon IIIc inclusion in cells expressing FGFR2-IIIb; hence, we have termed this element ISE/ISS-3 (for "intronic splicing enhancer-intronic splicing silencer 3"). Silencing of exon IIIc splicing by ISE/ISS-3 was shown to require a branch point sequence (BPS) using G as the primary branch nucleotide. Replacing a consensus BPS with A as the primary branch nucleotide resulted in constitutive splicing of exon IIIc. Our results suggest that the branch point sequence constitutes an important component that can contribute to the efficiency of exon definition of alternatively spliced cassette exons. Noncanonical branch points may thus facilitate cell-type-specific silencing of regulated exons by flanking cis elements.     

Characterization of the intronic splicing silencers flanking FGFR2 exon IIIb

The cell type-specific alternative splicing of FGFR2 pre-mRNA results in the mutually exclusive use of exons IIIb and IIIc, which leads to critically important differences in receptor function. The choice of exon IIIc in mesenchymal cells involves activation of this exon and repression of exon IIIb. This repression is mediated by the function of upstream and downstream intronic splicing silencers (UISS and DISS). Here we present a detailed characterization of the determinants of silencing function within UISS and DISS. We used a systematic mutational analysis, introducing deletions and substitutions to define discrete elements within these two silencers of exon IIIb. We show that UISS requires polypyrimidine tract-binding protein (PTB)-binding sites, which define the UISS1 sub-element, and an eight nucleotide sequence 5'-GCAGCACC-3' (UISS2) that is also required. Even though UISS2 does not bind PTB, the full UISS can be replaced with a synthetic silencer designed to provide optimal PTB binding. DISS is composed of a 5'-conserved sub-element (5'-CE) and two regions that contain multiple PTB sites and are functionally redundant (DISS1 and DISS2). DISS1 and DISS2 are separated by the activator sequence IAS2, and together these opposing elements form the intronic control element. Deletion of DISS in the FGFR2 exon IIIb context resulted in the near full inclusion of exon IIIb, and insertion of this silencer downstream of a heterologous exon with a weak 5' splice site was capable of repressing exon inclusion. Extensive deletion analysis demonstrated that the majority of silencing activity could be mapped to the conserved octamer CUCGGUGC within the 5'CE. Replacement of 5'CE and DISS1 with PTB-binding elements failed to restore repression of exon IIIb. We tested the importance of the relative position of the silencers and of the subelements within each silencer. Whereas UISS1, UISS2, DISS1, and DISS2 appear somewhat malleable, the 5'CE is rigid in terms of relative position and redundancy. Our data defined elements of function within the ISSs flanking exon IIIb and suggested that silencing of this exon is mediated by multiple trans-acting factors.