Promoter architecture modulates CFTR exon 9 skipping

Using hybrid minigene experiments, we have investigated the role of the promoter architecture on the regulation of two alternative spliced exons, cystic fibrosis transmembrane regulator (CFTR) exon 9 and fibronectin extra domain-A (EDB). A specific alternative splicing pattern corresponded to each analyzed promoter. Promoter-dependent sensitivity to cotransfected regulatory splicing factor SF2/ASF was observed only for the CFTR exon 9, whereas that of the EDB was refractory to promoter-mediated regulation. Deletion in the CFTR minigene of the downstream intronic splicing silencer element binding SF2/ASF abolished the specific promoter-mediated response to this splicing factor. A systematic analysis of the regulatory cis-acting elements showed that in the presence of suboptimal splice sites or by deletion of exonic enhancer elements the promoter-dependent sensitivity to splicing factor-mediated inhibition was lost. However, the basal regulatory effect of each promoter was preserved. The complex relationships between the promoter-dependent sensitivity to SF2 modulated by the exon 9 definition suggest a kinetic model of promoter-dependent alternative splicing regulation that possibly involves differential RNA polymerase II elongation.     

U5 snRNA interacts with exon sequences at 5' and 3' splice sites.

U5 snRNA is an essential pre-mRNA splicing factor whose function remains enigmatic. Specific mutations in a conserved single-stranded loop sequence in yeast U5 snRNA can activate cleavage of G1----A mutant pre-mRNAs at aberrant 5' splice sites and facilitate processing of dead-end lariat intermediates to mRNA. Activation of aberrant 5' cleavage sites involves base pairing between U5 snRNA and nucleotides upstream of the cleavage site. Processing of dead-end lariat intermediates to mRNA correlates with base pairing between U5 and the first two bases in exon 2. The loop sequence in U5 snRNA may therefore by intimately involved in the transesterification reactions at 5' and 3' splice sites. This pattern of interactions is strikingly reminiscent of exon recognition events in group II self-splicing introns and is consistent with the notion that U5 snRNA may be related to a specific functional domain from a group II-like self-splicing ancestral intron.     

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.     

U1-like snRNAs lacking complementarity to canonical 5' splice sites

We have detected a surprising heterogeneity among human spliceosomal U1 small nuclear RNA (snRNA). Most interestingly, we have identified three U1 snRNA variants that lack complementarity to the canonical 5' splice site (5'SS) GU dinucleotide. Furthermore, we have observed heterogeneity among the identified variant U1 snRNA genes caused by single nucleotide polymorphism (SNP). The identified snRNAs were ubiquitously expressed in a variety of human tissues representing different stages of development and displayed features of functional spliceosomal snRNAs, i.e., trimethylated cap structures, association with Sm proteins and presence in nuclear RNA-protein complexes. The unanticipated heterogeneity among spliceosomal snRNAs could contribute to the complexity of vertebrates by expanding the coding capacity of their genomes.     

Nuclear factor TDP-43 and SR proteins promote in vitro and in vivo CFTR exon 9 skipping

Alternative splicing of human cystic fibrosis transmembrane conductance regulator (CFTR) exon 9 is regulated by a combination of cis-acting elements distributed through the exon and both flanking introns (IVS8 and IVS9). Several studies have identified in the IVS8 intron 3' splice site a regulatory element that is composed of a polymorphic (TG)m(T)n repeated sequence. At present, no cellular factors have been identified that recognize this element. We have identified TDP-43, a nuclear protein not previously described to bind RNA, as the factor binding specifically to the (TG)m sequence. Transient TDP-43 overexpression in Hep3B cells results in an increase in exon 9 skipping. This effect is more pronounced with concomitant overexpression of SR proteins. Antisense inhibition of endogenous TDP-43 expression results in increased inclusion of exon 9, providing a new therapeutic target to correct aberrant splicing of exon 9 in CF patients. The clinical and biological relevance of this finding in vivo is demonstrated by our characterization of a CF patient carrying a TG10T9(DeltaF508)/TG13T3(wt) genotype leading to a disease-causing high proportion of exon 9 skipping.     

Biased exon/intron distribution of cryptic and de novo 3' splice sites

We compiled sequences of previously published aberrant 3′ splice sites (3′ss) that were generated by mutations in human disease genes. Cryptic 3′ss, defined here as those resulting from a mutation of the 3′YAG consensus, were more frequent in exons than in introns. They clustered in ~20 nt region adjacent to authentic 3′ss, suggesting that their under-representation in introns is due to a depletion of AG dinucleotides in the polypyrimidine tract (PPT). In contrast, most aberrant 3′ss that were induced by mutations outside the 3′YAG consensus (designated ‘de novo’) were in introns. The activation of intronic de novo 3′ss was largely due to AG-creating mutations in the PPT. In contrast, exonic de novo 3′ss were more often induced by mutations improving the PPT, branchpoint sequence (BPS) or distant auxiliary signals, rather than by direct AG creation. The Shapiro–Senapathy matrix scores had a good prognostic value for cryptic, but not de novo 3′ss. Finally, AG-creating mutations in the PPT that produced aberrant 3′ss upstream of the predicted BPS in vivo shared a similar ‘BPS-new AG’ distance. Reduction of this distance and/or the strength of the new AG PPT in splicing reporter pre-mRNAs improved utilization of authentic 3′ss, suggesting that AG-creating mutations that are located closer to the BPS and are preceded by weaker PPT may result in less severe splicing defects.     

Intrinsic differences between authentic and cryptic 5' splice sites

Cryptic splice sites are used only when use of a natural splice site is disrupted by mutation. To determine the features that distinguish authentic from cryptic 5′ splice sites (5′ss), we systematically analyzed a set of 76 cryptic 5′ss derived from 46 human genes. These cryptic 5′ss have a similar frequency distribution in exons and introns, and are usually located close to the authentic 5′ss. Statistical analysis of the strengths of the 5′ss using the Shapiro and Senapathy matrix revealed that authentic 5′ss have significantly higher score values than cryptic 5′ss, which in turn have higher values than the mutant ones. β-Globin provides an interesting exception to this rule, so we chose it for detailed experimental analysis in vitro. We found that the sequences of the β-globin authentic and cryptic 5′ss, but not their surrounding context, determine the correct 5′ss choice, although their respective scores do not reflect this functional difference. Our analysis provides a statistical basis to explain the competitive advantage of authentic over cryptic 5′ss in most cases, and should facilitate the development of tools to reliably predict the effect of disease-associated 5′ss-disrupting mutations at the mRNA level.     

5' exon requirement for self-splicing of the Tetrahymena thermophila pre-ribosomal RNA and identification of a cryptic 5' splice site in the 3' exon

The intervening sequence (IVS) of the Tetrahymena thermophila ribosomal RNA precursor undergoes accurate self-splicing in vitro. The work presented here examines the requirement for Tetrahymena rRNA sequences in the 5' exon for the accuracy and efficiency of splicing. Three plasmids were constructed with nine, four and two nucleotides of the natural 5' exon sequence, followed by the IVS and 26 nucleotides of the Tetrahymena 3' exon. RNA was transcribed from these plasmids in vitro and tested for self-splicing activity. The efficiency of splicing, as measured by the production of ligated exons, is reduced as the natural 5' exon sequence is replaced with plasmid sequences. Accurate splicing persists even when only four nucleotides of the natural 5' exon sequence remain. When only two nucleotides of the natural exon remain, no ligated exons are observed. As the efficiency of the normal reaction diminishes, novel RNA species are produced in increasing amounts. The novel RNA species were examined and found to be products of aberrant reactions of the precursor RNA. Two of these aberrant reactions involve auto-addition of GTP to sites six nucleotides and 52 nucleotides downstream from the 3' splice site. The former site occurs just after the sequence GGU, and may indicate the existence of a GGU-binding site within the IVS RNA. The latter site follows the sequence CUCU, which is identical with the four nucleotides preceding the 5' splice site. This observation led to a model where where the CUCU sequence in the 3' exon acts as a cryptic 5' splice site. The model predicted the existence of a circular RNA containing the first 52 nucleotides of the 3' exon. A small circular RNA was isolated and partially sequenced and found to support the model. So, a cryptic 5' splice site can function even if it is located downstream from the 3' splice site. Precursor RNA labeled at its 5' end, presumably by a GTP exchange reaction mediated by the IVS, is also described.     

Identification of alternative 5'/3' splice sites based on the mechanism of splice site competition

Alternative splicing plays an important role in regulating gene expression. Currently, most efficient methods use expressed sequence tags or microarray analysis for large-scale detection of alternative splicing. However, it is difficult to detect all alternative splice events with them because of their inherent limitations. Previous computational methods for alternative splicing prediction could only predict particular kinds of alternative splice events. Thus, it would be highly desirable to predict alternative 5'/3' splice sites with various splicing levels using genomic sequences alone. Here, we introduce the competition mechanism of splice sites selection into alternative splice site prediction. This approach allows us to predict not only rarely used but also frequently used alternative splice sites. On a dataset extracted from the AltSplice database, our method correctly classified approximately 70% of the splice sites into alternative and constitutive, as well as approximately 80% of the locations of real competitors for alternative splice sites. It outperforms a method which only considers features extracted from the splice sites themselves. Furthermore, this approach can also predict the changes in activation level arising from mutations in flanking cryptic splice sites of a given splice site. Our approach might be useful for studying alternative splicing in both computational and molecular biology.     

Determinants of the inherent strength of human 5' splice sites

We previously showed that the authentic 5' splice site (5'ss) of the first exon in the human beta-globin gene is intrinsically stronger than a cryptic 5'ss located 16 nucleotides upstream. Here we examined by mutational analysis the contribution of individual 5'ss nucleotides to discrimination between these two 5'ss. Based on the in vitro splicing efficiencies of a panel of 26 wild-type and mutant substrates in two separate 5'ss competition assays, we established a hierarchy of 5'ss and grouped them into three functional subclasses: strong, intermediate, and weak. Competition between two 5'ss from different subclasses always resulted in selection of the 5'ss that belongs to the stronger subclass. Moreover, each subclass has different characteristic features. Strong and intermediate 5'ss can be distinguished by their predicted free energy of base-pairing to the U1 snRNA 5' terminus (DeltaG). Whereas the extent of splicing via the strong 5'ss correlates well with the DeltaG, this is not the case for competition between intermediate 5'ss. Weak 5'ss were used only when the competing authentic 5'ss was inactivated by mutation. These results indicate that extensive complementarity to U1 snRNA exerts a dominant effect for 5'ss selection, but in the case of competing 5'ss with similarly modest complementarity to U1, the role of other 5'ss features is more prominent. This study reveals the importance of additional submotifs present in certain 5'ss sequences, whose characterization will be critical for understanding 5'ss selection in human genes.     

Modulation of exon skipping by high-affinity hnRNP A1-binding sites and by intron elements that repress splice site utilization

The RNA-binding protein hnRNP A1 is a splicing regulator produced by exclusion of alternative exon 7B from the A1 pre-mRNA. Each intron flanking exon 7B contains a high-affinity A1-binding site. The A1-binding elements promote exon skipping in vivo, activate distal 5' splice site selection in vitro and improve the responsiveness of pre-mRNAs to increases in the concentration of A1. Whereas the glycine-rich C-terminal domain of A1 is not required for binding, it is essential to activate the distal 5' splice site. Because A1 complexes can interact simultaneously with two A1-binding sites, we propose that an interaction between bound A1 proteins facilitates the pairing of distant splice sites. Based on the distribution of putative A1-binding sites in various pre-mRNAs, an A1-mediated change in pre-mRNA conformation may help define the borders of mammalian introns. We also identify an intron element which represses the 3' splice site of exon 7B. The activity of this element is mediated by a factor distinct from A1. Our results suggest that exon 7B skipping results from the concerted action of several intron elements that modulate splice site recognition and pairing.     

Cooperation of 5' and 3' processing sites as well as intron and exon sequences in calcitonin exon recognition.

We have previously shown that the calcitonin (CT)-encoding exon 4 of the human calcitonin/calcitonin gene-related peptide I (CGRP-I) gene (CALC-I gene) is surrounded by suboptimal processing sites. At the 5' end of exon 4 a weak 3' splice site is present because of an unusual branch acceptor nucleotide (U) and a weak poly(A) site is present at the 3' end of exon 4. For CT-specific RNA processing two different exon enhancer elements, A and B, located within exon 4 are required. In this study we have investigated the cooperation of these elements in CT exon recognition and inclusion by transient transfection into 293 cells of CALC-I minigene constructs. Improvement of the strength of the 3' splice site in front of exon 4 by the branchpoint mutation U-->A reduces the requirement for the presence of exon enhancer elements within exon 4 for CT-specific RNA processing, irrespective of the length of exon 4. Replacement of the exon 4 poly(A) site with a 5' splice site does not result in CT exon recognition, unless also one or more exon enhancer elements and/or the branchpoint mutation U-->A in front of exon 4 are present. This indicates that terminal and internal exons are recognised in a similar fashion. The number of additional enhancing elements that are required for CT exon recognition depends on the strength of the 5' splice site. Deletion of a large part of intron 4 also leads to partial exon 4 skipping. All these different elements contribute to CT exon recognition and inclusion. The CT exon is recognised as a whole entity and the sum of the strengths of the different elements determines recognition as an exon. Curiously, in one of our constructs a 5' splice site at the end of exon 4 is either ignored by the splicing machinery of the cell or recognised as a splice donor or as a splice acceptor site.     

Analysis and recognition of 5' UTR intron splice sites in human pre-mRNA

Prediction of splice sites in non-coding regions of genes is one of the most challenging aspects of gene structure recognition. We perform a rigorous analysis of such splice sites embedded in human 5' untranslated regions (UTRs), and investigate correlations between this class of splice sites and other features found in the adjacent exons and introns. By restricting the training of neural network algorithms to 'pure' UTRs (not extending partially into protein coding regions), we for the first time investigate the predictive power of the splicing signal proper, in contrast to conventional splice site prediction, which typically relies on the change in sequence at the transition from protein coding to non-coding. By doing so, the algorithms were able to pick up subtler splicing signals that were otherwise masked by 'coding' noise, thus enhancing significantly the prediction of 5' UTR splice sites. For example, the non-coding splice site predicting networks pick up compositional and positional bias in the 3' ends of non-coding exons and 5' non-coding intron ends, where cytosine and guanine are over-represented. This compositional bias at the true UTR donor sites is also visible in the synaptic weights of the neural networks trained to identify UTR donor sites. Conventional splice site prediction methods perform poorly in UTRs because the reading frame pattern is absent. The NetUTR method presented here performs 2-3-fold better compared with NetGene2 and GenScan in 5' UTRs. We also tested the 5' UTR trained method on protein coding regions, and discovered, surprisingly, that it works quite well (although it cannot compete with NetGene2). This indicates that the local splicing pattern in UTRs and coding regions is largely the same. The NetUTR method is made publicly available at www.cbs.dtu.dk/services/NetUTR.     

Conserved signals around the 5' splice sites in eukaryotic nuclear precursor mRNAs: G-runs are frequent in the introns and C in the exons near both 5' and 3' splice sites

It is known that the GT doublet is well conserved at the 5' exon/intron splice junction and is frequently embedded in the AGGT quartet. Although only the underlined G is invariable, splicing and ligation are accurately executed. In this work we search for additional conserved potential signals which may aid in 5' splice site recognition. Extensive searches which are not limited to a preconceived consensus sequence are carried out. We investigate the distributions of the 256 quartets in a 1000 nucleotide span around the 5' splice sites in approximately 1700 eukaryotic nuclear precursor mRNAs. Several potential signals are noted. Of particular interest are quartets containing runs of G, e.g., G4, G3T, G3C, G3A and AG3 in the intron immediately downstream and some C-containing quartets in the exon upstream of the 5' splice site. In an analogous calculation, (A)GGG(A) has also been found to be frequent in the intron, 60 nucleotides upstream and (A)CCC(A) in the exon downstream of the 3' splice site. These results are consistent with the recent indications that exon sequences may play a role in efficient splicing. Some models are proposed.