Changes in exon-intron structure during vertebrate evolution affect the splicing pattern of exons

Exon-intron architecture is one of the major features directing the splicing machinery to the short exons that are located within long flanking introns. However, the evolutionary dynamics of exon-intron architecture and its impact on splicing is largely unknown. Using a comparative genomic approach, we analyzed 17 vertebrate genomes and reconstructed the ancestral motifs of both 3' and 5' splice sites, as also the ancestral length of exons and introns. Our analyses suggest that vertebrate introns increased in length from the shortest ancestral introns to the longest primate introns. An evolutionary analysis of splice sites revealed that weak splice sites act as a restrictive force keeping introns short. In contrast, strong splice sites allow recognition of exons flanked by long introns. Reconstruction of the ancestral state suggests these phenomena were not prevalent in the vertebrate ancestor, but appeared during vertebrate evolution. By calculating evolutionary rate shifts in exons, we identified cis-acting regulatory sequences that became fixed during the transition from early vertebrates to mammals. Experimental validations performed on a selection of these hexamers confirmed their regulatory function. We additionally revealed many features of exons that can discriminate alternative from constitutive exons. These features were integrated into a machine-learning approach to predict whether an exon is alternative. Our algorithm obtains very high predictive power (AUC of 0.91), and using these predictions we have identified and successfully validated novel alternatively spliced exons. Overall, we provide novel insights regarding the evolutionary constraints acting upon exons and their recognition by the splicing machinery.     

Evolution of exon-intron structure and alternative splicing in fruit flies and malarial mosquito genomes

Comparative analysis of alternative splicing of orthologous genes from fruit flies (Drosophila melanogaster and Drosophila pseudoobscura) and mosquito (Anopheles gambiae) demonstrated that both in the fruit fly genes and in fruit fly-mosquito comparisons, constitutive exons and splicing sites are more conserved than alternative ones. While >97% of constitutive D. melanogaster exons are conserved in D. pseudoobscura, only approximately 80% of alternative exons are conserved. Similarly, 77% of constitutive fruit fly exons are conserved in the mosquito genes, compared with <50% of alternative exons. Internal alternatives are more conserved than terminal ones. Retained introns are the least conserved, alternative acceptor sites are slightly more conserved than donor sites, and mutually exclusive exons are almost as conserved as constitutive exons. Cassette and mutually exclusive exons experience almost no intron insertions. We also observed cases of interconversion of various elementary alternatives, e.g., transformation of cassette exons into alternative sites. These results agree with the observations made earlier in human-mouse comparisons and demonstrate that the phenomenon of relatively low conservation of alternatively spliced regions may be universal, as it has been observed in different taxonomic groups (mammals and insects) and at various evolutionary distances.     

A systematic analysis of intronic sequences downstream of 5' splice sites reveals a widespread role for U-rich motifs and TIA1/TIAL1 proteins in alternative splicing regulation

To identify human intronic sequences associated with 5' splice site recognition, we performed a systematic search for motifs enriched in introns downstream of both constitutive and alternative cassette exons. Significant enrichment was observed for U-rich motifs within 100 nucleotides downstream of 5' splice sites of both classes of exons, with the highest enrichment between positions +6 and +30. Exons adjacent to U-rich intronic motifs contain lower frequencies of exonic splicing enhancers and higher frequencies of exonic splicing silencers, compared with exons not followed by U-rich intronic motifs. These findings motivated us to explore the possibility of a widespread role for U-rich motifs in promoting exon inclusion. Since cytotoxic granule-associated RNA binding protein (TIA1) and TIA1-like 1 (TIAL1; also known as TIAR) were previously shown in vitro to bind to U-rich motifs downstream of 5' splice sites, and to facilitate 5' splice site recognition in vitro and in vivo, we investigated whether these factors function more generally in the regulation of splicing of exons followed by U-rich intronic motifs. Simultaneous knockdown of TIA1 and TIAL1 resulted in increased skipping of 36/41 (88%) of alternatively spliced exons associated with U-rich motifs, but did not affect 32/33 (97%) alternatively spliced exons that are not associated with U-rich motifs. The increase in exon skipping correlated with the proximity of the first U-rich motif and the overall "U-richness" of the adjacent intronic region. The majority of the alternative splicing events regulated by TIA1/TIAL1 are conserved in mouse, and the corresponding genes are associated with diverse cellular functions. Based on our results, we estimate that approximately 15% of alternative cassette exons are regulated by TIA1/TIAL1 via U-rich intronic elements.     

Prediction of Mutant mRNA Splice Isoforms by Information Theory‐Based Exon Definition

Mutations that affect mRNA splicing often produce multiple mRNA isoforms, resulting in complex molecular phenotypes. Definition of an exon and its inclusion in mature mRNA relies on joint recognition of both acceptor and donor splice sites. This study predicts cryptic and exon‐skipping isoforms in mRNA produced by splicing mutations from the combined information contents (R_i, which measures binding‐site strength, in bits) and distribution of the splice sites defining these exons. The total information content of an exon (R_i,total) is the sum of the R_i values of its acceptor and donor splice sites, adjusted for the self‐information of the distance separating these sites, that is, the gap surprisal. Differences between total information contents of an exon (ΔR_i,total) are predictive of the relative abundance of these exons in distinct processed mRNAs. Constraints on splice site and exon selection are used to eliminate nonconforming and poorly expressed isoforms. Molecular phenotypes are computed by the Automated Splice Site and Exon Definition Analysis ( http://splice.uwo.ca ) server. Predictions of splicing mutations were highly concordant (85.2%; n = 61) with published expression data. In silico exon definition analysis will contribute to streamlining assessment of abnormal and normal splice isoforms resulting from mutations.     

Splicing in action: assessing disease causing sequence changes

Variations in new splicing regulatory elements are difficult to identify exclusively by sequence inspection and may result in deleterious effects on precursor (pre) mRNA splicing. These mutations can result in either complete skipping of the exon, retention of the intron, or the introduction of a new splice site within an exon or intron. Sometimes mutations that do not disrupt or create a splice site activate pre-existing pseudo splice sites, consistent with the proposal that introns contain splicing inhibitory sequences. These variants can also affect the fine balance of isoforms produced by alternatively spliced exons and in consequence cause disease. Available genomic pathology data reveal that we are still partly ignorant of the basic mechanisms that underlie the pre-mRNA splicing process. The fact that human pathology can provide pointers to new modulatory elements of splicing should be exploited.     

Dichotomous splicing signals in exon flanks

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

Utilisation of a cryptic non-canonical donor splice site of the gene encoding PARAFIBROMIN is associated with familial isolated primary hyperparathyroidism

More than 99% of all splice sites conform to consensus sequences that usually include the invariant dinucleotides gt and ag at the 5' and 3' ends of the introns, respectively. We report on the utilisation of a non-consensus (non-canonical) donor splice site within exon 1 of the HRPT2 gene in familial isolated primary hyperparathyroidism (FIHP). HRPT2 mutations are more frequently associated with the hyperparathyroidism-jaw tumour syndrome (HPT-JT). Patients with FIHP were identified to have a donor splice site mutation, IVS1+1 g-->a, and the consequences of this for RNA processing were investigated. The mutant mRNA lacked 30 bp and DNA sequence analysis revealed this to result from utilisation of an alternative cryptic non-canonical donor splice site (gaatgt) in exon 1 together with the normally occurring acceptor splice site in intron 1. Translation of this mutant mRNA predicted the in-frame loss of 10 amino acids in the encoded protein, termed PARAFIBROMIN. Thus, these FIHP patients are utilising a ga-ag splice site pair, which until recently was considered to be incompatible with splicing but is now known to occur as a rare (<0.02%) normal splicing variant.     

Identification of two distinct intron elements involved in alternative splicing of beta-tropomyosin pre-mRNA

The rat beta-tropomyosin gene encodes two isoforms, termed skeletal muscle beta-tropomyosin and fibroblast last tropomyosim 1 (TM-1), via an alternative RNA processing mechanism. The gene contains 11 exons. Exons 1-5 and exons 8 and 9 are common to all mRNAs expressed from the gene. Exons 6 and 11 are used in fibroblasts, as well as smooth muscle, whereas exons 7 and 10 are used only in skeletal muscle. In the present studies we focused on the mutually exclusive internal alternative splice choice involving exon 6 (fibroblast-type splice) and exon 7 (skeletal muscle-type splice). We have identified two distinct elements in the intron, upstream of exon 7, involved in splice site selection. The first element is comprised of a polypyrimidine tract located 89-143 nucleotides upstream of the 3' splice site, which specifies the location of the lariat branchpoints used, 144-153 nucleotides upstream of exon 7. The 3' splice site AG dinucleotide has no role in selection of these branchpoints. The second element is comprised of intron sequences located between the polypyrimidine tract and the 3' splice site of exon 7. It contains an important determinant in alternative splice site selection, because deletion of these sequences results in the use of the skeletal muscle-specific exon in nonmuscle cells. We propose that the use of lariat branchpoints located far upstream from a 3' splice site may be a general feature of some alternatively excised introns, reflecting the presence of regulatory sequences located between the lariat branch site and the 3' splice site. The data also indicate that alternative splicing of the rat beta-tropomyosin gene is regulated by a somewhat different mechanism from that described for rat alpha-tropomyosin gene and the transformer-2 gene of Drosophila melanogaster.     

Ordered partitioning reveals extended splice-site consensus information

Using recently available cDNA and genomic data (Berkeley Drosophila Genome Project; http://www.fruitfly.org), we computed a large sample of 10,057 Drosophila splice sites. An information-theoretic analysis of the nucleotide sequences adjacent to these splice sites showed a strong correlation between the sizes of introns and exons and the levels of information, which is a measure of sequence conservation. The strong correlation permitted us to determine extensive consensus sequences at the donor and acceptor sites of longer introns. These sequences were further refined and extended by examining the information in regions around splice sites that only partially matched the consensus. The correlation between length and information provided the basis for determining alternative consensus arrangements associated with shorter introns, as well as general base-composition preferences that likely promote spliceosome function. We also observed a correlation between information near splice sites and the lengths of nonadjacent introns, indicating that there are long-range effects spanning multiple introns. The ordered partitioning approach used in this analysis may become increasingly useful as large genomic data sets become available.     

A point mutation in the human CD45 gene associated with defective splicing of exon A

CD45 is a receptor-type protein tyrosine phosphatase involved in the regulation of lymphocyte activation. Different CD45 isoforms are generated by alternative splicing of three variable exons (A, B and C). The pattern of CD45 splicing depends upon cell type and state of activation. CD45RA isoforms (containing exon A-encoded sequences) can usually be found on a subset of resting T cells, but not on activated T cells. We have recently described a variant pattern of CD45RA expression which is characterized by continuous expression of CD45RA molecules on activated and memory T cells. Here, we demonstrate that this phenotype is associated with heterozygosity for a point mutation at nucleotide position 77 of exon A, leading to a C → G transition. This mutation does not change the protein sequence of the CD45RA isoform. We conclude that position 77 is part of a motif necessary for splicing of exon A, which supports the hypothesis that sequences within exons have significant effects on alternative splicing. The mutation of this motif might prevent binding of a transacting splice factor. In the heterozygous state, this mutation is not associated with impaired T cell reactivity. Functional consequences of the homozygous state remain to be elucidated.     

Detecting tissue-specific alternative splicing and disease-associated aberrant splicing of the PTCH gene with exon junction microarrays

Mutations in the human ortholog of Drosophila patched (PTCH) have been identified in patients with autosomal dominant nevoid basal cell carcinoma syndrome (NBCCS), characterized by minor developmental anomalies and an increased incidence of cancers such as medulloblastoma and basal cell carcinoma. We identified many isoforms of PTCH mRNA involving exons 1-5, exon 10 and a novel exon, 12b, generated by alternative splicing (AS), most of which have not been deposited in GenBank nor discussed earlier. To monitor splicing events of the PTCH gene, we designed oligonucleotide arrays on which exon probes and exon-exon junction probes as well as a couple of intron probes for the PTCH gene were placed in duplicate. Probe intensities were normalized on the basis of the total expression of PTCH and probe sensitivity. Tissue-specific regulation of AS identified with the microarrays closely correlated with the results obtained by RT-PCR. Of note, the novel exon, exon 12b, was specifically expressed in the brain and heart, especially in the cerebellum. Additionally, using these microarrays, we were able to detect disease-associated aberrant splicings of the PTCH gene in two patients with NBCCS. In both cases, cryptic splice donor sites located either in an exon or in an intron were activated because of the partial disruption of the consensus sequence for the authentic splice donor sites due to point mutations. Taken together, oligonucleotide microarrays containing exon junction probes are demonstrated to be a powerful tool to investigate tissue-specific regulation of AS and aberrant splicing taking place in genetic disorders.