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.     

Splicing of phenylalanine hydroxylase (PAH) exon 11 is vulnerable: molecular pathology of mutations in PAH exon 11

In about 20-30% of phenylketonuria (PKU) patients, phenylalanine (Phe) levels can be controlled by cofactor 6R-tetrahydrobiopterin (BH(4)) administration. The phenylalanine hydroxylase (PAH) genotype has a predictive value concerning BH(4)-response and therefore a correct assessment of the mutation molecular pathology is important. Mutations that disturb the splicing of exons (e.g. interplay between splice site strength and regulatory sequences like exon splicing enhancers (ESEs)/exon splicing silencers (ESSs)) may cause different severity of PKU. In this study, we identified PAH exon 11 as a vulnerable exon and used patient derived lymphoblast cell lines and PAH minigenes to study the molecular defect that impacted pre-mRNA processing. We showed that the c.1144T>C and c.1066-3C>T mutations cause exon 11 skipping, while the c.1139C>T mutation is neutral or slightly beneficial. The c.1144T>C mutation resides in a putative splicing enhancer motif and binding by splicing factors SF2/ASF, SRp20 and SRp40 is disturbed. Additional mutations in potential splicing factor binding sites contributed to elucidate the pathogenesis of mutations in PAH exon 11. We suggest that PAH exon 11 is vulnerable due to a weak 3' splice site and that this makes exon 11 inclusion dependent on an ESE spanning position c.1144. Importantly, this implies that other mutations in exon 11 may affect splicing, since splicing is often determined by a fine balance between several positive and negative splicing regulatory elements distributed throughout the exon. Finally, we identified a pseudoexon in intron 11, which would have pathogenic consequences if activated by mutations or improved splicing conditions. Exonic mutations that disrupt splicing are unlikely to facilitate response to BH(4) and may lead to inconsistent genotype-phenotype correlations. Therefore, recognizing such mutations enhances our ability to predict the BH(4)-response.     

Functional analysis of cis-acting elements regulating the alternative splicing of human CFTR exon 9.

The rate of exon 9 exclusion from the cystic fibrosis transmembrane conductance regulator (CFTR) mRNA is associated with monosymptomatic forms of cystic fibrosis. Exon 9 alternative splicing is modulated by a polymorphic polythymidine tract within its 3' splice site. We have generated a minigene carrying human CFTR exon 9 with its flanking intronic sequences and set up an in vivo model to study the cis-acting DNA elements which modulate its splicing. Transfections into human cell lines showed that T5, but not T9 or T7 alleles, significantly increases the alternative splicing of exon 9. Moreover, we found that another polymorphic locus juxtaposed upstream of the T tract, and constituted by (TG)(n)repeats, can further modulate exon 9 skipping but only when activated by the T5 allele. Then, we extended our studies to the mouse CFTR exon 9 which does not show alternative splicing. Comparison of human and mouse introns 8 and 9 revealed a low homology between the two sequences and the absence of the human polymorphic loci within the mouse intron 3' splice site. We have tested a series of constructs where the whole human exon 9 with its flanking intronic sequences was replaced partially or completely by the murine counterpart. The transfections of these constructs in human and murine cell lines reveal that also sequences of the downstream intron 9 affect exon 9 definition and co-modulate, with the UG/U 3' splice site sequences, the extent of exon 9 skipping in CFTR mRNA.     

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.     

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.     

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.     

Porcine T-cell receptor beta-chain: a genomic sequence covering Dbeta1.1 to Cbeta2 gene segments and the diversity of cDNA expressed in piglets including novel alternative splicing products

Porcine TCRbeta-chain cDNA clones were isolated from thymic and peripheral blood lymphocytes of piglets. Using these nucleotide sequences, a genomic 18kbp sequence stretch covering Dbeta1 to Cbeta2 gene segments was identified, which revealed that the porcine TCRbeta-chain locus consists of two sets of Dbeta-Jbeta-Cbeta gene groups with each set having a Dbeta gene segment, seven Jbeta gene segments and a down stream Cbeta gene segment composed of four exons. This structure is consistent with other known mammalian TCRbeta-chain loci. With this genomic information, TCRbeta-chain clones from cDNA libraries were analyzed. Sixteen Vbeta gene segments were obtained accompanied by either Dbeta1 or Dbeta2 and by one of the nine Jbeta gene segments. Five different Cbeta cDNA sequences were obtained including four types of Cbeta1 sequences and one type of Cbeta2 sequence. The differences among the Cbeta1 sequences are either allelic polymorphisms or two splice variants, one being a product of exon1 splicing to exon3 (exon2 skipping), and another being an alternative splicing using a splice acceptor site newly discovered inside Cbeta1 exon4. The latter splice acceptor site was also found in human, mouse and horse all giving short cytoplasmic domain with Phe at their C-terminal ends. Other splicing products included trans-splicing of Jbeta2 to Cbeta1, non-functional splicing of two Jbeta gene segments in tandem and a part of Jbeta2.7-Cbeta2 intron to Cbeta2 exon1. Numerous examples of splice variants may suggest the involvement of splicing in generating TCRbeta-chain functional diversity.     

Functional significance of a deep intronic mutation in the ATM gene and evidence for an alternative exon 28a

Screening for ATM mutations is usually performed using genomic DNA as a template for PCR amplification across exonic regions, with the consequence that deep intronic sequences are not analyzed. Here we report a novel pseudoexon-retaining deep intronic mutation (IVS28-159A>G; g.75117A>G based on GenBank U82828.1) in a patient with ataxia-telangiectasia (A-T), as well as the identification of a previously unrecognized alternative exon in the ATM gene (exon 28a) expressed in lymphoblastoid cell lines (LCL) derived from normal individuals. cDNA analysis using the A-T patient's LCL showed the retention of two aberrant intronic segments of 112 and 190 nt between exons 28 and 29. Minigenes were constructed to determine the functional significance of two genomic changes in the region of aberrant splicing: IVS28-193C>T (g.75083C>T) and IVS28-159A>G, revealing that: 1) the first is a polymorphism; 2) IVS28-159A>G weakens the 5' splice site of the alternative exon 28a and activates a cryptic 5' splice site (ss) 83 nt downstream; and 3) wild-type constructs also retain a 29-nt segment (exon 28a) as part of both the 112- and 190-nt segments. Maximum entropy estimates of ss strengths corroborate the cDNA and minigene findings. Such mutations may prove relevant in planning therapy that targets specific splicing aberrations.     

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.     

Signals and their transduction pathways regulating alternative splicing: a new dimension of the human genome

The human genome project demonstrated that alternative splicing of genes is more the rule than the exception. Missplicing events are an important cause and indication of human disease. Changing alternative splicing patterns in response to an external stimulus seems to be a physiological process performed by many cells. Organisms regulate alternative splice site selection by changing the concentration and activity of splicing regulatory proteins. This is achieved by de novo protein synthesis, by regulation of the intracellular localization and by phosphorylation.     

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. ( 2011 ). 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.     

DHPLC-based mutation analysis of ENG and ALK-1 genes in HHT Italian population

Hereditary haemorrhagic telangiectasia (HHT or Rendu-Osler-Weber syndrome) is an autosomal dominant disorder characterized by localized angiodysplasia due to mutations in endoglin, ALK-1 gene, and a still unidentified locus. The lack of highly recurrent mutations, locus heterogeneity, and the presence of mutations in almost all coding exons of the two genes makes the screening for mutations time-consuming and costly. In the present study, we developed a DHPLC-based protocol for mutation detection in ALK1 and ENG genes through retrospective analysis of known sequence variants, 20 causative mutations and 11 polymorphisms, and a prospective analysis on 47 probands with unknown mutation. Overall DHPLC analysis identified the causative mutation in 61 out 66 DNA samples (92.4%). We found 31 different mutations in the ALK1 gene, of which 15 are novel, and 20, of which 12 are novel, in the ENG gene, thus providing for the first time the mutational spectrum in a cohort of Italian HHT patients. In addition, we characterized the splicing pattern of ALK1 gene in lymphoblastoid cells, both in normal controls and in two individuals carrying a mutation in the non-invariant -3 position of the acceptor splice site upstream exon 6 (c.626-3C>G). Functional essay demonstrated the existence, also in normal individuals, of a small proportion of ALK1 alternative splicing, due to exon 5 skipping, and the presence of further aberrant splicing isoforms in the individuals carrying the c.626-3C>G mutation.