Aberrant 3' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

Norepinephrine deficiency is caused by combined abnormal mRNA processing and defective protein trafficking of dopamine beta-hydroxylase

Human norepinephrine (NE) deficiency (or dopamine β-hydroxylase (DBH) deficiency) is a rare congenital disorder of primary autonomic failure, in which neurotransmitters NE and epinephrine are undetectable. Although potential pathogenic mutations, such as a common splice donor site mutation (IVS1+2T→C) and various missense mutations, in NE deficiency patients were identified, molecular mechanisms underlying this disease remain unknown. Here, we show that the IVS1+2T→C mutation results in a non-detectable level of DBH protein production and that all three missense mutations tested lead to the DBH protein being trapped in the endoplasmic reticulum (ER). Supporting the view that mutant DBH induces an ER stress response, exogenous expression of mutant DBH dramatically induced expression of BiP, a master ER chaperone. Furthermore, we found that a pharmacological chaperone, glycerol, significantly rescued defective trafficking of mutant DBH proteins. Taken together, we propose that NE deficiency is caused by the combined abnormal processing of DBH mRNA and defective protein trafficking and that this disease could be treated by a pharmacological chaperone(s).

Identification of a MIP mutation that activates a cryptic acceptor splice site in the 3' untranslated region

PURPOSE

To investigate the consequence of a major intrinsic protein MIP splice-site mutation (c.607-1G>A) in a four-generation Chinese pedigree afflicted with autosomal dominant congenital cataracts (ADCC).

METHODS

Both a mutated minigene with c.607-1G>A, and a wild-type minigene were constructed using the pTARGET mammalian expression vector. They were transiently transfected into HeLa cells and human lens epithelial cells, respectively. After 48 h incubation, RNA extraction and RT-PCR analysis were performed and PCR products were separated and confirmed by sequencing. Structural models of the wild-type and the mutant aquaporin 0 (AQP0) were generated and analyzed using SWISS-MODEL.

RESULTS

The G>A transition activated a cryptic acceptor splice site (c.965-966) in the 3' untranslated region (3' UTR), resulting in the absence of the coding region and most of the 3'UTR in exon 4 of the mature mRNA. Moreover, homology modeling of the mutant protein suggested that the sixth transmembrane helix and carboxyl terminus were replaced with the Leu-His-Ser tripeptide (AQP0-LHS).

CONCLUSIONS

The MIP splice-site mutation (c.607-1G>A) activates a cryptic acceptor splice site in the 3' UTR, which may result in substitution of the sixth transmembrane helix and carboxyl terminus for AQP0-LHS. To our knowledge, this is the first report of activation of a cryptic splice site in the 3' UTR in a human disease gene.

Genome-wide data-mining of candidate human splice translational efficiency polymorphisms (STEPs) and an online database

BACKGROUND

Variation in pre-mRNA splicing is common and in some cases caused by genetic variants in intronic splicing motifs. Recent studies into the insulin gene (INS) discovered a polymorphism in a 5' non-coding intron that influences the likelihood of intron retention in the final mRNA, extending the 5' untranslated region and maintaining protein quality. Retention was also associated with increased insulin levels, suggesting that such variants--splice translational efficiency polymorphisms (STEPs)--may relate to disease phenotypes through differential protein expression. We set out to explore the prevalence of STEPs in the human genome and validate this new category of protein quantitative trait loci (pQTL) using publicly available data.

METHODOLOGY/PRINCIPAL FINDINGS

Gene transcript and variant data were collected and mined for candidate STEPs in motif regions. Sequences from transcripts containing potential STEPs were analysed for evidence of splice site recognition and an effect in expressed sequence tags (ESTs). 16 publicly released genome-wide association data sets of common diseases were searched for association to candidate polymorphisms with HapMap frequency data. Our study found 3324 candidate STEPs lying in motif sequences of 5' non-coding introns and further mining revealed 170 with transcript evidence of intron retention. 21 potential STEPs had EST evidence of intron retention or exon extension, as well as population frequency data for comparison.

CONCLUSIONS/SIGNIFICANCE

Results suggest that the insulin STEP was not a unique example and that many STEPs may occur genome-wide with potentially causal effects in complex disease. An online database of STEPs is freely accessible at http://dbstep.genes.org.uk/.

DBASS3 and DBASS5: databases of aberrant 3'- and 5'-splice sites

DBASS3 and DBASS5 provide comprehensive repositories of new exon boundaries that were induced by pathogenic mutations in human disease genes. Aberrant 5'- and 3'-splice sites were activated either by mutations in the consensus sequences of natural exon-intron junctions (cryptic sites) or elsewhere ('de novo' sites). DBASS3 and DBASS5 currently contain approximately 900 records of cryptic and de novo 3'- and 5'-splice sites that were produced by over a thousand different mutations in approximately 360 genes. DBASS3 and DBASS5 data can be searched by disease phenotype, gene, mutation, location of aberrant splice sites in introns and exons and their distance from authentic counterparts, by bibliographic references and by the splice-site strength estimated with several prediction algorithms. The user can also retrieve reference sequences of both aberrant and authentic splice sites with the underlying mutation. These data will facilitate identification of introns or exons frequently involved in aberrant splicing, mutation analysis of human disease genes and study of germline or somatic mutations that impair RNA processing. Finally, this resource will be useful for fine-tuning splice-site prediction algorithms, better definition of auxiliary splicing signals and design of new reporter assays. DBASS3 and DBASS5 are freely available at http://www.dbass.org.uk/.

Targeted genome-wide enrichment of functional regions

Only a small fraction of large genomes such as that of the human contains the functional regions such as the exons, promoters, and polyA sites. A platform technique for selective enrichment of functional genomic regions will enable several next-generation sequencing applications that include the discovery of causal mutations for disease and drug response. Here, we describe a powerful platform technique, termed “functional genomic fingerprinting” (FGF), for the multiplexed genomewide isolation and analysis of targeted regions such as the exome, promoterome, or exon splice enhancers. The technique employs a fixed part of a uniquely designed Fixed-Randomized primer, while the randomized part contains all the possible sequence permutations. The Fixed-Randomized primers bind with full sequence complementarity at multiple sites where the fixed sequence (such as the splice signals) occurs within the genome, and multiplex amplify many regions bounded by the fixed sequences (e.g., exons). Notably, validation of this technique using cardiac myosin binding protein-C (MYBPC3) gene as an example strongly supports the application and efficacy of this method. Further, assisted by genomewide computational analyses of such sequences, the FGF technique may provide a unique platform for high-throughput sample production and analysis of targeted genomic regions by the next-generation sequencing techniques, with powerful applications in discovering disease and drug response genes.

Some novel intron positions in conserved Drosophila genes are caused by intron sliding or tandem duplication

BACKGROUND

Positions of spliceosomal introns are often conserved between remotely related genes. Introns that reside in non-conserved positions are either novel or remnants of frequent losses of introns in some evolutionary lineages. A recent gain of such introns is difficult to prove. However, introns verified as novel are needed to evaluate contemporary processes of intron gain.

RESULTS

We identified 25 unambiguous cases of novel intron positions in 31 Drosophila genes that exhibit near intron pairs (NIPs). Here, a NIP consists of an ancient and a novel intron position that are separated by less than 32 nt. Within a single gene, such closely-spaced introns are very unlikely to have coexisted. In most cases, therefore, the ancient intron position must have disappeared in favour of the novel one. A survey for NIPs among 12 Drosophila genomes identifies intron sliding (migration) as one of the more frequent causes of novel intron positions. Other novel introns seem to have been gained by regional tandem duplications of coding sequences containing a proto-splice site.

CONCLUSIONS

Recent intron gains sometimes appear to have arisen by duplication of exonic sequences and subsequent intronization of one of the copies. Intron migration and exon duplication together may account for a significant amount of novel intron positions in conserved coding sequences.

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.

Genomic features defining exonic variants that modulate splicing

A comparative analysis of SNPs and their exonic and intronic environments identifies the features predictive of splice affecting variants.

Background

Single point mutations at both synonymous and non-synonymous positions within exons can have severe effects on gene function through disruption of splicing. Predicting these mutations in silico purely from the genomic sequence is difficult due to an incomplete understanding of the multiple factors that may be responsible. In addition, little is known about which computational prediction approaches, such as those involving exonic splicing enhancers and exonic splicing silencers, are most informative.

Results

We assessed the features of single-nucleotide genomic variants verified to cause exon skipping and compared them to a large set of coding SNPs common in the human population, which are likely to have no effect on splicing. Our findings implicate a number of features important for their ability to discriminate splice-affecting variants, including the naturally occurring density of exonic splicing enhancers and exonic splicing silencers of the exon and intronic environment, extensive changes in the number of predicted exonic splicing enhancers and exonic splicing silencers, proximity to the splice junctions and evolutionary constraint of the region surrounding the variant. By extending this approach to additional datasets, we also identified relevant features of variants that cause increased exon inclusion and ectopic splice site activation.

Conclusions

We identified a number of features that have statistically significant representation among exonic variants that modulate splicing. These analyses highlight putative mechanisms responsible for splicing outcome and emphasize the role of features important for exon definition. We developed a web-tool, Skippy, to score coding variants for these relevant splice-modulating features.

Novel exon nucleotide substitution at the splice junction causes a neonatal Marfan syndrome

The fibrillin-1 gene (FBN1) mutations are associated with a broad spectrum of disorders including Marfan syndrome (MFS) and show great clinical heterogeneity. An underrepresentation for mutations leading to premature termination codon (PTC) in FBN1 exons 24-32 was found in neonatal or severe MFS but the underlying cause was unclear. This study thoroughly examined two FBN1 mutations on exons 24-32 region to illustrate the molecular mechanisms underlying these FBN1 mutations on MFS etiology. Two nucleotide substitutions, c.3208G> C, the last nucleotide of exon 26, and c.3209A>G, the first nucleotide of exon 27, affecting the same amino acid, p.D1070H and p.D1070G, respectively, gave very different phenotypes. We demonstrate that c.3208G>C generates two alternatively spliced transcripts, while c.3209A>G does not affect the splicing. We further demonstrate that the aberrantly spliced transcripts do not go through nonsense-mediated decay, but rather produce unstable, premature protein peptides that are degraded by endoplasmic reticulum associated degradation. The molecular mechanism outlined here defines a model for the pathogenesis of PTC-containing mutation within the exons 24-32 of FBN1 in MFS. Furthermore, our data suggest that PTC mutation within this region may lead to early lethality in neonatal MFS.

A method of predicting changes in human gene splicing induced by genetic variants in context of cis-acting elements

Background

Polymorphic variants and mutations disrupting canonical splicing isoforms are among the leading causes of human hereditary disorders. While there is a substantial evidence of aberrant splicing causing Mendelian diseases, the implication of such events in multi-genic disorders is yet to be well understood. We have developed a new tool (SpliceScan II) for predicting the effects of genetic variants on splicing and cis-regulatory elements. The novel Bayesian non-canonical 5'GC splice site (SS) sensor used in our tool allows inference on non-canonical exons.

Results

Our tool performed favorably when compared with the existing methods in the context of genes linked to the Autism Spectrum Disorder (ASD). SpliceScan II was able to predict more aberrant splicing isoforms triggered by the mutations, as documented in DBASS5 and DBASS3 aberrant splicing databases, than other existing methods. Detrimental effects behind some of the polymorphic variations previously associated with Alzheimer's and breast cancer could be explained by changes in predicted splicing patterns.

Conclusions

We have developed SpliceScan II, an effective and sensitive tool for predicting the detrimental effects of genomic variants on splicing leading to Mendelian and complex hereditary disorders. The method could potentially be used to screen resequenced patient DNA to identify de novo mutations and polymorphic variants that could contribute to a genetic disorder.

Dual effect of a single nucleotide polymorphism in the first intron of the porcine secreted phosphoprotein 1 gene: allele-specific binding of C/EBP beta and activation of aberrant splicing

Background

Secreted phosphoprotein 1 (SPP1 or Osteopontin, OPN) is a multifunctional matricellular glycoprotein involved in development and regeneration of skeletal muscle. Previously, we have demonstrated that porcine SPP1 shows breed-related differential mRNA expression during myogenesis. With the aim to identify putative contributing cis-regulatory DNA variation we resequenced the 5' upstream region of the gene in the respective breeds Pietrain and Duroc. We found two single nucleotide polymorphisms (SNP; [GenBank:M84121]: g.1804C>T and g.3836A>G). We focused our investigation on the SNP g.3836A>G, because in silico analysis and knowledge about the regulation of SPP1 suggested an effect of this SNP on a CCAAT/enhancer binding protein beta (C/EBPβ) responsive transcriptional enhancer.

Results

Using electrophoretic mobility shift assay we demonstrated that, similar to human SPP1, the 3' terminal end of the first intron of porcine SPP1 harbors a C/EBPβ binding site and showed that this binding site is negatively affected by the mutant G allele. Genotyping of 48 fetuses per breed revealed that the G allele segregated exclusively in Duroc fetuses with a frequency of 57 percent. Using real-time quantitative PCR we showed that, consistent with its negative effect on a transcriptional enhancer element, the G allele tends to decrease mRNA abundance of SPP1 in the fetal musculus longissimus dorsi (~1.3 fold; P ≥ 0.1).

Moreover, we showed that the SNP g.3836A>G leads to ubiquitous aberrant splicing of the first intron by generating a de novo and activating a cryptic splice acceptor site. Aberrantly spliced transcripts comprise about half of the SPP1 messages expressed by the G allele. Both aberrant splice variants differ from the native transcript by insertions in the leader sequences which do not change the reading frame of SPP1.

Conclusion

At the 3' terminal end of the first intron of the porcine SPP1 we identified a unique, dually functional SNP g.3836A>G. This SNP affects the function of the SPP1 gene at the DNA level by affecting a C/EBPβ binding site and at the RNA level by activating aberrant splicing of the first intron, and thus represents an interesting DNA-marker to study phenotypic effects of SPP1 DNA-variation.

Transposable elements in disease-associated cryptic exons

Transposable elements (TEs) make up a half of the human genome, but the extent of their contribution to cryptic exon activation that results in genetic disease is unknown. Here, a comprehensive survey of 78 mutation-induced cryptic exons previously identified in 51 disease genes revealed the presence of TEs in 40 cases (51%). Most TE-containing exons were derived from short interspersed nuclear elements (SINEs), with Alus and mammalian interspersed repeats (MIRs) covering >18 and >16% of the exonized sequences, respectively. The majority of SINE-derived cryptic exons had splice sites at the same positions of the Alu/MIR consensus as existing SINE exons and their inclusion in the mRNA was facilitated by phylogenetically conserved changes that improved both traditional and auxiliary splicing signals, thus marking intronic TEs amenable for pathogenic exonization. The overrepresentation of MIRs among TE exons is likely to result from their high average exon inclusion levels, which reflect their strong splice sites, a lack of splicing silencers and a high density of enhancers, particularly (G)AA(G) motifs. These elements were markedly depleted in antisense Alu exons, had the most prominent position on the exon-intron gradient scale and are proposed to promote exon definition through enhanced tertiary RNA interactions involving unpaired (di)adenosines. The identification of common mechanisms by which the most dynamic parts of the genome contribute both to new exon creation and genetic disease will facilitate detection of intronic mutations and the development of computational tools that predict TE hot-spots of cryptic exon activation.

Disease-causing mutations improving the branch site and polypyrimidine tract: pseudoexon activation of LINE-2 and antisense Alu lacking the poly(T)-tail

Cryptic exons or pseudoexons are typically activated by point mutations that create GT or AG dinucleotides of new 5' or 3' splice sites in introns, often in repetitive elements. Here we describe two cases of tetrahydrobiopterin deficiency caused by mutations improving the branch point sequence and polypyrimidine tracts of repeat-containing pseudoexons in the PTS gene. In the first case, we demonstrate a novel pathway of antisense Alu exonization, resulting from an intronic deletion that removed the poly(T)-tail of antisense AluSq. The deletion brought a favorable branch point sequence within proximity of the pseudoexon 3' splice site and removed an upstream AG dinucleotide required for the 3' splice site repression on normal alleles. New Alu exons can thus arise in the absence of poly(T)-tails that facilitated inclusion of most transposed elements in mRNAs by serving as polypyrimidine tracts, highlighting extraordinary flexibility of Alu repeats in shaping intron-exon structure. In the other case, a PTS pseudoexon was activated by an A>T substitution 9 nt upstream of its 3' splice site in a LINE-2 sequence, providing the first example of a disease-causing exonization of the most ancient interspersed repeat. These observations expand the spectrum of mutational mechanisms that introduce repetitive sequences in mature transcripts and illustrate the importance of intronic mutations in alternative splicing and phenotypic variability of hereditary disorders.

Ab initio prediction of mutation-induced cryptic splice-site activation and exon skipping

Mutations that affect splicing of precursor messenger RNAs play a major role in the development of hereditary diseases. Most splicing mutations have been found to eliminate GT or AG dinucleotides that define the 5' and 3' ends of introns, leading to exon skipping or cryptic splice-site activation. Although accurate description of the mis-spliced transcripts is critical for predicting phenotypic consequences of these alterations, their exact nature in affected individuals cannot often be determined experimentally. Using a comprehensive collection of exons that sustained cryptic splice-site activation or were skipped as a result of splice-site mutations, we have developed a multivariate logistic discrimination procedure that distinguishes the two aberrant splicing outcomes from DNA sequences. The new algorithm was validated using an independent sample of exons and implemented as a free online utility termed CRYP-SKIP (http://www.dbass.org.uk/cryp-skip/). The web application takes up one or more mutated alleles, each consisting of one exon and flanking intronic sequences, and provides a list of important predictor variables and their values, the overall probability of activating cryptic splice vs exon skipping, and the location and intrinsic strength of predicted cryptic splice sites in the input sequence. These results will facilitate phenotypic prediction of splicing mutations and provide further insights into splicing enhancer and silencer elements and their relative importance for splice-site selection in vivo.

Deleterious variants of FIG4, a phosphoinositide phosphatase, in patients with ALS

Mutations of the lipid phosphatase FIG4 that regulates PI(3,5)P(2) are responsible for the recessive peripheral-nerve disorder CMT4J. We now describe nonsynonymous variants of FIG4 in 2% (9/473) of patients with amyotrophic lateral sclerosis (ALS) and primary lateral sclerosis (PLS). Heterozygosity for a deleterious allele of FIG4 appears to be a risk factor for ALS and PLS, extending the list of known ALS genes and increasing the clinical spectrum of FIG4-related diseases.

Prediction and assessment of splicing alterations: implications for clinical testing

Sequence variants that may result in splicing alterations are a particular class of inherited variants for which consequences can be more readily assessed, using a combination of bioinformatic prediction methods and in vitro assays. There is also a general agreement that a variant would invariably be considered pathogenic on the basis of convincing evidence that it results in transcript(s) carrying a premature stop codon or an in-frame deletion disrupting known functional domain(s). This commentary discusses current practices used to assess the clinical significance of this class of variants, provides suggestions to improve assessment, and highlights the issues involved in routine assessment of potential splicing aberrations. We conclude that classification of sequence variants that may alter splicing is greatly enhanced by supporting in vitro analysis. Additional studies that assess large numbers of variants for induction of splicing aberrations and exon skipping are needed to define the contribution of splicing/exon skipping to cancer and disease. These studies will also provide the impetus for development of algorithms that better predict splicing patterns. To facilitate variant classification and development of more specific bioinformatic tools, we call for the deposition of all laboratory data from splicing analyses into national and international databases.

Genome-wide survey of allele-specific splicing in humans

Background

Accurate mRNA splicing depends on multiple regulatory signals encoded in the transcribed RNA sequence. Many examples of mutations within human splice regulatory regions that alter splicing qualitatively or quantitatively have been reported and allelic differences in mRNA splicing are likely to be a common and important source of phenotypic diversity at the molecular level, in addition to their contribution to genetic disease susceptibility. However, because the effect of a mutation on the efficiency of mRNA splicing is often difficult to predict, many mutations that cause disease through an effect on splicing are likely to remain undiscovered.

Results

We have combined a genome-wide scan for sequence polymorphisms likely to affect mRNA splicing with analysis of publicly available Expressed Sequence Tag (EST) and exon array data. The genome-wide scan uses published tools and identified 30,977 SNPs located within donor and acceptor splice sites, branch points and exonic splicing enhancer elements. For 1,185 candidate splicing polymorphisms the difference in splicing between alternative alleles was corroborated by publicly available exon array data from 166 lymphoblastoid cell lines. We developed a novel probabilistic method to infer allele-specific splicing from EST data. The method uses SNPs and alternative mRNA isoforms mapped to EST sequences and models both regulated alternative splicing as well as allele-specific splicing. We have also estimated heritability of splicing and report that a greater proportion of genes show evidence of splicing heritability than show heritability of overall gene expression level. Our results provide an extensive resource that can be used to assess the possible effect on splicing of human polymorphisms in putative splice-regulatory sites.

Conclusion

We report a set of genes showing evidence of allele-specific splicing from an integrated analysis of genomic polymorphisms, EST data and exon array data, including several examples for which there is experimental evidence of polymorphisms affecting splicing in the literature. We also present a set of novel allele-specific splicing candidates and discuss the strengths and weaknesses of alternative technologies for inferring the effect of sequence variants on mRNA splicing.

Patterns of missplicing caused by RB1 gene mutations in patients with retinoblastoma and association with phenotypic expression

We have analyzed RNA from retinoblastoma patients and unaffected carriers with various RB1 gene mutations to determine the patterns of missplicing and associations with phenotypic expression. Most sequence alterations in or in the neighborhood of conserved splice signals that we tested resulted in simple exon skipping (15 mutations) or intron inclusion (new acceptor AG-sites, four mutations) as expected. Two mutations resulted in skipping of a neighboring exon (exon 11), a complex pattern indicating competition for correct lariat formation. We observed no activation of a cryptic splice site but found that a recurrent missense mutation in exon 7 creates a new splice site (two families). RT-PCR analysis enabled us to confirm the presence and to characterize the transcriptional consequences of gross insertions and deletions in the RB1 gene in six patients, including two patients with mutational mosaicism. We also used RT-PCR analysis to search for unknown mutations in 15 patients and identified three oncogenic point mutations deep in introns. Two of these mutations are recurrent thus indicating that, despite the vast extent of the introns of the RB1 gene, few bases are effective targets for oncogenic mutations. When analyzing associations between phenotypic expression (16 families) and mutational consequences we observed no link to the presence or absence of a premature termination codon in the mutant transcript. However, the location of a mutation relative to the splice sequence has a strong and consistent influence on phenotypic expression.

Computational analysis of splicing errors and mutations in human transcripts

Background

Most retained introns found in human cDNAs generated by high-throughput sequencing projects seem to result from underspliced transcripts, and thus they capture intermediate steps of pre-mRNA splicing. On the other hand, mutations in splice sites cause exon skipping of the respective exon or activation of pre-existing cryptic sites. Both types of events reflect properties of the splicing mechanism.

Results

The retained introns were significantly shorter than constitutive ones, and skipped exons are shorter than exons with cryptic sites. Both donor and acceptor splice sites of retained introns were weaker than splice sites of constitutive introns. The authentic acceptor sites affected by mutations were significantly weaker in exons with activated cryptic sites than in skipped exons. The distance from a mutated splice site to the nearest equivalent site is significantly shorter in cases of activated cryptic sites compared to exon skipping events. The prevalence of retained introns within genes monotonically increased in the 5'-to-3' direction (more retained introns close to the 3'-end), consistent with the model of co-transcriptional splicing. The density of exonic splicing enhancers was higher, and the density of exonic splicing silencers lower in retained introns compared to constitutive ones and in exons with cryptic sites compared to skipped exons.

Conclusion

Thus the analysis of retained introns in human cDNA, exons skipped due to mutations in splice sites and exons with cryptic sites produced results consistent with the intron definition mechanism of splicing of short introns, co-transcriptional splicing, dependence of splicing efficiency on the splice site strength and the density of candidate exonic splicing enhancers and silencers. These results are consistent with other, recently published analyses.

Conserved and species-specific alternative splicing in mammalian genomes

Background

Alternative splicing has been shown to be one of the major evolutionary mechanisms for protein diversification and proteome expansion, since a considerable fraction of alternative splicing events appears to be species- or lineage-specific. However, most studies were restricted to the analysis of cassette exons in pairs of genomes and did not analyze functionality of the alternative variants.

Results

We analyzed conservation of human alternative splice sites and cassette exons in the mouse and dog genomes. Alternative exons, especially minor-isofom ones, were shown to be less conserved than constitutive exons. Frame-shifting alternatives in the protein-coding regions are less conserved than frame-preserving ones. Similarly, the conservation of alternative sites is highest for evenly used alternatives, and higher when the distance between the sites is divisible by three. The rate of alternative-exon and site loss in mouse is slightly higher than in dog, consistent with faster evolution of the former. The evolutionary dynamics of alternative sites was shown to be consistent with the model of random activation of cryptic sites.

Conclusion

Consistent with other studies, our results show that minor cassette exons are less conserved than major-alternative and constitutive exons. However, our study provides evidence that this is caused not only by exon birth, but also lineage-specific loss of alternative exons and sites, and it depends on exon functionality.

Global control of aberrant splice-site activation by auxiliary splicing sequences: evidence for a gradient in exon and intron definition

Auxiliary splicing signals play a major role in the regulation of constitutive and alternative pre-mRNA splicing, but their relative importance in selection of mutation-induced cryptic or de novo splice sites is poorly understood. Here, we show that exonic sequences between authentic and aberrant splice sites that were activated by splice-site mutations in human disease genes have lower frequencies of splicing enhancers and higher frequencies of splicing silencers than average exons. Conversely, sequences between authentic and intronic aberrant splice sites have more enhancers and less silencers than average introns. Exons that were skipped as a result of splice-site mutations were smaller, had lower SF2/ASF motif scores, a decreased availability of decoy splice sites and a higher density of silencers than exons in which splice-site mutation activated cryptic splice sites. These four variables were the strongest predictors of the two aberrant splicing events in a logistic regression model. Elimination or weakening of predicted silencers in two reporters consistently promoted use of intron-proximal splice sites if these elements were maintained at their original positions, with their modular combinations producing expected modification of splicing. Together, these results show the existence of a gradient in exon and intron definition at the level of pre-mRNA splicing and provide a basis for the development of computational tools that predict aberrant splicing outcomes.

Analbuminemia Produced by a Novel Splicing Mutation

Analbuminemia is a rare autosomal recessive disorder manifested by the absence or severe reduction of circulating human serum albumin in homozygous or compound heterozygous individuals. It is an allelic heterogeneous defect, caused by a variety of mutations within the albumin gene. The analbuminemic condition was diagnosed in a Turkish female infant on the basis of low albumin concentration (∼9.0 g/L). The albumin gene was screened by single-strand conformation polymorphism and heteroduplex analysis and submitted to direct sequencing. The proband was found to be homozygous for a T→C transition at nucleotide 13381, the 2nd base of intron 11. The effect of this previously unreported mutation, which inactivates the strongly conserved GT dinucleotide at the 5′ splice site consensus sequence of intron 11, was evaluated by examining the cDNA obtained by reverse transcription-PCR from the albumin mRNA extracted from the proband leukocytes. This analysis revealed that the mutation, named Bartin for the geographical origin of the patient’s family, results in the skipping of exon 11. The subsequent frameshift within exon 12 originates a premature stop codon located 5 codons downstream at position 411. The predicted translation product would consist of 410 amino acids. This novel extensive cDNA alteration is responsible for the analbuminemic trait.

Aberrant 5' splice sites in human disease genes: mutation pattern, nucleotide structure and comparison of computational tools that predict their utilization

Despite a growing number of splicing mutations found in hereditary diseases, utilization of aberrant splice sites and their effects on gene expression remain challenging to predict. We compiled sequences of 346 aberrant 5′splice sites (5′ss) that were activated by mutations in 166 human disease genes. Mutations within the 5′ss consensus accounted for 254 cryptic 5′ss and mutations elsewhere activated 92 de novo 5′ss. Point mutations leading to cryptic 5′ss activation were most common in the first intron nucleotide, followed by the fifth nucleotide. Substitutions at position +5 were exclusively G>A transitions, which was largely attributable to high mutability rates of C/G>T/A. However, the frequency of point mutations at position +5 was significantly higher than that observed in the Human Gene Mutation Database, suggesting that alterations of this position are particularly prone to aberrant splicing, possibly due to a requirement for sequential interactions with U1 and U6 snRNAs. Cryptic 5′ss were best predicted by computational algorithms that accommodate nucleotide dependencies and not by weight-matrix models. Discrimination of intronic 5′ss from their authentic counterparts was less effective than for exonic sites, as the former were intrinsically stronger than the latter. Computational prediction of exonic de novo 5′ss was poor, suggesting that their activation critically depends on exonic splicing enhancers or silencers. The authentic counterparts of aberrant 5′ss were significantly weaker than the average human 5′ss. The development of an online database of aberrant 5′ss will be useful for studying basic mechanisms of splice-site selection, identifying splicing mutations and optimizing splice-site prediction algorithms.

A deep intronic mutation in the RB1 gene leads to intronic sequence exonisation

Familial forms of retinoblastoma, an embryonic neoplasm of retinal origin, are caused by constitutional mutations of the RB1 gene. In this paper, we describe a family with retinoblastoma affecting two brothers with no previous family history of cancer. Complete RB1 mutational screening including point mutation and large rearrangement screening failed to demonstrate any mutation. The whole coding sequence was therefore investigated at the cDNA level, demonstrating a 103 bp intronic insertion between exons 23 and 24, leading to subsequent frameshift and premature termination of translation. This intronic exonisation was caused by a deep intronic mutation in intron 23 generating a cryptic 3' splice site. This is the first report of a deep intronic mutation in RB1 and is a proof of concept that some undetected RB1 mutations should be investigated at the cDNA level, particularly in hereditary forms of retinoblastoma.