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

Identification and validation of lncRNA mutation hotspot SNPs associated with myasthenia gravis susceptibility

Background

Myasthenia gravis (MG) is an autoimmune disorder caused by antibodies that target the postsynaptic muscle membrane. Recent evidence suggests that genetic variants and long noncoding RNAs (lncRNAs) play crucial roles in the pathogenesis of MG. The purpose of this study was to investigate the associations between lncRNA-related single-nucleotide polymorphisms (SNPs) and MG susceptibility in Chinese populations.

Methods

First, we identified lncRNA mutation hotspot regions based on the improved Kolmogorov‒Smirnov test and the cumulative hypergeometric distribution principle. Next, we further identified lncRNA mutation hotspot SNPs by calculating conservative scores. Finally, experiments were conducted to verify the associations between lncRNA mutation hotspot SNPs and MG susceptibility. A total of 82 patients with MG and 82 healthy controls were recruited for genotyping of lncRNA mutation hotspot SNPs using the SNaPshot technique. Quantitative real-time PCR was used to investigate lncRNA expression in 34 patients with MG and 37 healthy controls.

Results

In the multistep calculation, 14 candidate SNPs of 3 lncRNAs (AL031686.1, NONHSAT028539.2 and AC245014.3) in MG were identified as mutation hotspot SNPs. The genotyping results of the 14 SNPs in our study revealed no statistically significant differences in the frequencies of genotypes and alleles between patients with MG and controls. However, in the lncRNA AL031686.1, rs1000383 and rs6094353 were in perfect linkage disequilibrium (LD) and were associated with an increased risk of ocular MG. Additionally, rs6094347 was associated with an increased risk of ocular MG. Nevertheless, no SNP was found to be associated with factors such as sex, age, the presence or absence of thymoma, or the genetic model of MG. Further experiments revealed that NONHSAT028539.2 expression was upregulated in peripheral blood mononuclear cells (PBMCs) from patients with MG compared with those from healthy controls.

Conclusion

In our study, we did not find an association between the 14 mutation hotspot SNPs of lncRNAs and susceptibility to MG. However, we observed that the rs6094347 and rs1000383/rs6094353 polymorphisms in the lncRNA AL031686.1 were associated with the risk of ocular MG.

Homozygous splice-site variant in ENPP1 underlies generalized arterial calcification of infancy

ENPP1 (ectonucleotide pyrophosphatase/phosphodiesterase 1) plays a critical role by converting extracellular ATP to AMP, generating extracellular PPi, a potential inhibitor of calcification. Pathogenic variants in the ENPP1 cause generalized arterial calcification of infancy (GACI [OMIM 208000]). GACI, is an ultra-rare disease characterized by early-onset calcification of large and medium-sized arteries, leading to severe cardiovascular complications such as heart failure, pulmonary stenosis (PS), hypertension, and more. In this study, we report a novel homozygous splice-site pathogenic variant in ENPP1 (NM_006208, c.2230 + 5G > A; p.Asp701Asnfs*2) residing in C-terminal nuclease-like domain (NLD) of ENPP1 protein in a Pakistani family diagnosed with severe valvular PS and mild right ventricular hypertrophy (RVH). cDNA assays confirmed the skipping of exon 21, and the splice product underwent nonsense-mediated decay. Functional studies on fibroblasts from the patient demonstrated increased calcification and decreased enzymatic activity of ENPP1, recapitulating the hallmarks of GACI. By combining genetic analysis with the in vitro study, we substantiate that ENPP1:c.2230 + 5G > A variant is pathogenic, underscoring its role in the development of GACI.

Advanced variant classification framework reduces the false positive rate of predicted loss-of-function variants in population sequencing data

Predicted loss of function (pLoF) variants are often highly deleterious and play an important role in disease biology, but many pLoF variants may not result in loss of function (LoF). Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines' PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 genes associated with autosomal-recessive disease from the Genome Aggregation Database (gnomAD v.2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in a low proportion expressed across transcripts (pext) scored region, or the presence of cryptic in-frame splice rescues. Variants predicted to evade LoF or to be potential artifacts were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of pLoF variants predicted as likely not LoF/not LoF, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

Advanced variant classification framework reduces the false positive rate of predicted loss of function (pLoF) variants in population sequencing data

Predicted loss of function (pLoF) variants are highly deleterious and play an important role in disease biology, but many of these variants may not actually result in loss-of-function. Here we present a framework that advances interpretation of pLoF variants in research and clinical settings by considering three categories of LoF evasion: (1) predicted rescue by secondary sequence properties, (2) uncertain biological relevance, and (3) potential technical artifacts. We also provide recommendations on adjustments to ACMG/AMP guidelines's PVS1 criterion. Applying this framework to all high-confidence pLoF variants in 22 autosomal recessive disease-genes from the Genome Aggregation Database (gnomAD, v2.1.1) revealed predicted LoF evasion or potential artifacts in 27.3% (304/1,113) of variants. The major reasons were location in the last exon, in a homopolymer repeat, in low per-base expression (pext) score regions, or the presence of cryptic splice rescues. Variants predicted to be potential artifacts or to evade LoF were enriched for ClinVar benign variants. PVS1 was downgraded in 99.4% (162/163) of LoF evading variants assessed, with 17.2% (28/163) downgraded as a result of our framework, adding to previous guidelines. Variant pathogenicity was affected (mostly from likely pathogenic to VUS) in 20 (71.4%) of these 28 variants. This framework guides assessment of pLoF variants beyond standard annotation pipelines, and substantially reduces false positive rates, which is key to ensure accurate LoF variant prediction in both a research and clinical setting.

Combining genetic constraint with predictions of alternative splicing to prioritize deleterious splicing in rare disease studies

BACKGROUND

Despite numerous molecular and computational advances, roughly half of patients with a rare disease remain undiagnosed after exome or genome sequencing. A particularly challenging barrier to diagnosis is identifying variants that cause deleterious alternative splicing at intronic or exonic loci outside of canonical donor or acceptor splice sites.

RESULTS

Several existing tools predict the likelihood that a genetic variant causes alternative splicing. We sought to extend such methods by developing a new metric that aids in discerning whether a genetic variant leads to deleterious alternative splicing. Our metric combines genetic variation in the Genome Aggregate Database with alternative splicing predictions from SpliceAI to compare observed and expected levels of splice-altering genetic variation. We infer genic regions with significantly less splice-altering variation than expected to be constrained. The resulting model of regional splicing constraint captures differential splicing constraint across gene and exon categories, and the most constrained genic regions are enriched for pathogenic splice-altering variants. Building from this model, we developed ConSpliceML. This ensemble machine learning approach combines regional splicing constraint with multiple per-nucleotide alternative splicing scores to guide the prediction of deleterious splicing variants in protein-coding genes. ConSpliceML more accurately distinguishes deleterious and benign splicing variants than state-of-the-art splicing prediction methods, especially in "cryptic" splicing regions beyond canonical donor or acceptor splice sites.

CONCLUSION

Integrating a model of genetic constraint with annotations from existing alternative splicing tools allows ConSpliceML to prioritize potentially deleterious splice-altering variants in studies of rare human diseases.

Performance evaluation of differential splicing analysis methods and splicing analytics platform construction

A proportion of previously defined benign variants or variants of uncertain significance in humans, which are challenging to identify, may induce an abnormal splicing process. An increasing number of methods have been developed to predict splicing variants, but their performance has not been completely evaluated using independent benchmarks. Here, we manually sourced ∼50 000 positive/negative splicing variants from > 8000 studies and selected the independent splicing variants to evaluate the performance of prediction methods. These methods showed different performances in recognizing splicing variants in donor and acceptor regions, reminiscent of different weight coefficient applications to predict novel splicing variants. Of these methods, 66.67% exhibited higher specificities than sensitivities, suggesting that more moderate cut-off values are necessary to distinguish splicing variants. Moreover, the high correlation and consistent prediction ratio validated the feasibility of integration of the splicing prediction method in identifying splicing variants. We developed a splicing analytics platform called SPCards, which curates splicing variants from publications and predicts splicing scores of variants in genomes. SPCards also offers variant-level and gene-level annotation information, including allele frequency, non-synonymous prediction and comprehensive functional information. SPCards is suitable for high-throughput genetic identification of splicing variants, particularly those located in non-canonical splicing regions.

Donor Splice Site Variant in SLC9A6 Causes Christianson Syndrome in a Lithuanian Family: A Case Report

Background and Objectives: The pathogenic variants of SLC9A6 are a known cause of a rare, X-linked neurological disorder called Christianson syndrome (CS). The main characteristics of CS are developmental delay, intellectual disability, and neurological findings. This study investigated the genetic basis and explored the molecular changes that led to CS in two male siblings presenting with intellectual disability, epilepsy, behavioural problems, gastrointestinal dysfunction, poor height, and weight gain. Materials and Methods: Next-generation sequencing of a tetrad was applied to identify the DNA changes and Sanger sequencing of proband’s cDNA was used to evaluate the impact of a splice site variant on mRNA structure. Bioinformatical tools were used to investigate SLC9A6 protein structure changes. Results: Sequencing and bioinformatical analysis revealed a novel donor splice site variant (NC_000023.11(NM_001042537.1):c.899 + 1G > A) that leads to a frameshift and a premature stop codon. Protein structure modelling showed that the truncated protein is unlikely to form any functionally relevant SLC9A6 dimers. Conclusions: Molecular and bioinformatical analysis revealed the impact of a novel donor splice site variant in the SLC9A6 gene that leads to truncated and functionally disrupted protein causing the phenotype of CS in the affected individuals.

Genome-wide identification of exon extension/shrinkage events induced by splice-site-creating mutations

Mutations that affect phenotypes have been identified primarily as those that directly alter amino acid sequences or disrupt splice sites. However, some mutations not located in functionally important sites can also affect phenotypes, such as splice-site-creating mutations (SCMs). To investigate how frequent exon extension/shrinkage events induced by SCMs occur in normal individuals, we used personal genome sequencing data and transcriptome data of the corresponding individuals and identified 371 exon extension/shrinkage events in normal individuals. This number was about three times higher than the number of pseudo-exon activation events identified in the previous study. The average numbers of exon extension and exon shrinkage events in each sample were 3.3 and 11.2, respectively. We also evaluated the impact of exon extension/shrinkage events on the resulting transcripts and their protein products and found that 40.2% of the identified events may have possible functional impacts by either generating premature termination codons in transcripts or affecting protein domains. Our results indicated that a certain fraction of SCMs identified in this study can be pathogenic mutations by creating novel splice sites.

Background splicing as a predictor of aberrant splicing in genetic disease

Mutations of splice sites, auxiliary splicing elements and the splicing machinery cause a wide range of genetic disease. Here we report that many of the complex effects of splicing mutations can be predicted from background splicing information, with emphasis on BRCA1, BRCA2 and DMD. Background splicing arises from very low level splicing between rarely used background splice sites and from low-level exon skipping between intron splice sites. We show how this information can be downloaded from the Snaptron database of spliced RNA, which we then compared with databases of human splice site mutations. We report that inactivating mutations of intron splice sites typically caused the non-mutated partner splice site to splice to a known background splice site in over 90% of cases and to the strongest background splice site in the large majority of cases. Consequently, background splicing information can usefully predict the effects of splice site mutations, which include cryptic splice activation and single or multiple exon skipping. In addition, de novo splice sites and splice sites involved in pseudoexon formation, recursive splicing and aberrant splicing in cancer show a 90% match to background splice sites, so establishing that the enhancement of background splicing causes a wide range of splicing aberrations. We also discuss how background splicing information can identify cryptic splice sites that might be usefully targeted by antisense oligonucleotides (ASOs) and how it might indicate possible multiple exon skipping side effects of ASOs designed to induce single exon skipping.

What's Wrong in a Jump? Prediction and Validation of Splice Site Variants

Alternative splicing (AS) is a crucial process to enhance gene expression driving organism development. Interestingly, more than 95% of human genes undergo AS, producing multiple protein isoforms from the same transcript. Any alteration (e.g., nucleotide substitutions, insertions, and deletions) involving consensus splicing regulatory sequences in a specific gene may result in the production of aberrant and not properly working proteins. In this review, we introduce the key steps of splicing mechanism and describe all different types of genomic variants affecting this process (splicing variants in acceptor/donor sites or branch point or polypyrimidine tract, exonic, and deep intronic changes). Then, we provide an updated approach to improve splice variants detection. First, we review the main computational tools, including the recent Machine Learning-based algorithms, for the prediction of splice site variants, in order to characterize how a genomic variant interferes with splicing process. Next, we report the experimental methods to validate the predictive analyses are defined, distinguishing between methods testing RNA (transcriptomics analysis) or proteins (proteomics experiments). For both prediction and validation steps, benefits and weaknesses of each tool/procedure are accurately reported, as well as suggestions on which approaches are more suitable in diagnostic rather than in clinical research.

IntSplice2: Prediction of the Splicing Effects of Intronic Single-Nucleotide Variants Using LightGBM Modeling

Prediction of the effect of a single-nucleotide variant (SNV) in an intronic region on aberrant pre-mRNA splicing is challenging except for an SNV affecting the canonical GU/AG splice sites (ss). To predict pathogenicity of SNVs at intronic positions −50 (Int-50) to −3 (Int-3) close to the 3’ ss, we developed light gradient boosting machine (LightGBM)-based IntSplice2 models using pathogenic SNVs in the human gene mutation database (HGMD) and ClinVar and common SNVs in dbSNP with 0.01 ≤ minor allelic frequency (MAF) < 0.50. The LightGBM models were generated using features representing splicing cis-elements. The average recall/sensitivity and specificity of IntSplice2 by fivefold cross-validation (CV) of the training dataset were 0.764 and 0.884, respectively. The recall/sensitivity of IntSplice2 was lower than the average recall/sensitivity of 0.800 of IntSplice that we previously made with support vector machine (SVM) modeling for the same intronic positions. In contrast, the specificity of IntSplice2 was higher than the average specificity of 0.849 of IntSplice. For benchmarking (BM) of IntSplice2 with IntSplice, we made a test dataset that was not used to train IntSplice. After excluding the test dataset from the training dataset, we generated IntSplice2-BM and compared it with IntSplice using the test dataset. IntSplice2-BM was superior to IntSplice in all of the seven statistical measures of accuracy, precision, recall/sensitivity, specificity, F1 score, negative predictive value (NPV), and matthews correlation coefficient (MCC). We made the IntSplice2 web service at https://www.med.nagoya-u.ac.jp/neurogenetics/IntSplice2.

A novel mutation in ext2 caused hereditary multiple exostoses through reducing the synthesis of heparan sulfate

Hereditary multiple exostoses (HME) is a rare skeletal disorder characterized by the formation of multiple benign cartilage-capped tumors, usually in the metaphyseal region of the long bones. Over 70% of HME cases arise from monoallelic mutations in either of the two genes encoding the heparan sulfate (HS) synthesis enzymes, ext1 and ext2. To identify more HME-associated mutations, genomic DNA from members of five independent consanguineous families with HME was sequenced with whole exome sequencing (WES). A novel heterozygous splice site mutation (c.1173+2T>A) in ext2 was detected in all three affected members of family V. Further study showed that the novel mutation caused exon 7 of ext2 mRNA to be skipped during splicing and caused a frameshift after the codon for Arg360, which results in the appearance of new 43 codons, followed by a termination codon. Although the resulting truncated protein was still localized to the Golgi, similar to the full-length EXT2, its HS synthesis activity decreased by 40%. In this study, a novel splice site mutation in ext2 was identified and suggested to be a pathogenic mutation of HME, which may expand the genetic etiology spectrum of HME and may be helpful for clinical genetic counseling and prenatal diagnosis.

Exome sequencing implicates genetic disruption of prenatal neuro-gliogenesis in sporadic congenital hydrocephalus

Congenital hydrocephalus (CH), characterized by enlarged brain ventricles, is considered a disease of excessive cerebrospinal fluid (CSF) accumulation and thereby treated with neurosurgical CSF diversion with high morbidity and failure rates. The poor neurodevelopmental outcomes and persistence of ventriculomegaly in some post-surgical patients highlight our limited knowledge of disease mechanisms. Through whole-exome sequencing of 381 patients (232 trios) with sporadic, neurosurgically treated CH, we found that damaging de novo mutations account for >17% of cases, with five different genes exhibiting a significant de novo mutation burden. In all, rare, damaging mutations with large effect contributed to ~22% of sporadic CH cases. Multiple CH genes are key regulators of neural stem cell biology and converge in human transcriptional networks and cell types pertinent for fetal neuro-gliogenesis. These data implicate genetic disruption of early brain development, not impaired CSF dynamics, as the primary pathomechanism of a significant number of patients with sporadic CH.

Transposon clusters as substrates for aberrant splice-site activation

Transposed elements (TEs) have dramatically shaped evolution of the exon-intron structure and significantly contributed to morbidity, but how recent TE invasions into older TEs cooperate in generating new coding sequences is poorly understood. Employing an updated repository of new exon-intron boundaries induced by pathogenic mutations, termed DBASS, here we identify novel TE clusters that facilitated exon selection. To explore the extent to which such TE exons maintain RNA secondary structure of their progenitors, we carried out structural studies with a composite exon that was derived from a long terminal repeat (LTR78) and AluJ and was activated by a C > T mutation optimizing the 5ʹ splice site. Using a combination of SHAPE, DMS and enzymatic probing, we show that the disease-causing mutation disrupted a conserved AluJ stem that evolved from helix 3.3 (or 5b) of 7SL RNA, liberating a primordial GC 5ʹ splice site from the paired conformation for interactions with the spliceosome. The mutation also reduced flexibility of conserved residues in adjacent exon-derived loops of the central Alu hairpin, revealing a cross-talk between traditional and auxilliary splicing motifs that evolved from opposite termini of 7SL RNA and were approximated by Watson-Crick base-pairing already in organisms without spliceosomal introns. We also identify existing Alu exons activated by the same RNA rearrangement. Collectively, these results provide valuable TE exon models for studying formation and kinetics of pre-mRNA building blocks required for splice-site selection and will be useful for fine-tuning auxilliary splicing motifs and exon and intron size constraints that govern aberrant splice-site activation.

A homozygote variant in the tRNA splicing endonuclease subunit 54 causes pontocerebellar hypoplasia in a consanguineous Iranian family

Background

Homozygous loss‐of‐function mutations in TSEN54 (tRNA splicing endonuclease subunit 54; OMIM: 608755) cause different types of pontocerebellar hypoplasias (PCH) including PCH2, PCH4, and PCH5. The study aimed to determine the possible genetic factors contributing to PCH phenotypes in two affected male infants in an Iranian family.

Methods

We subjected two affected individuals in a consanguineous Iranian family. To systematically investigate the susceptible gene(s), whole‐exome sequencing was performed on the proband and a novel identified variant was confirmed by Sanger sequencing. We also analyzed 26 relatives in three generations using PCR‐restriction fragment length polymorphism (PCR‐RFLP) followed and confirmed by Sanger sequencing.

Results

Physical and medical examinations confirmed PCH in the patients. Besides, the proband showed bilateral moderate sensorineural hearing loss and structural heart defects as the novel phenotypes. The molecular findings also verified that two affected individuals were homozygote for the novel synonymous variant, NM_207346.2: c.1170G>A; p.(Val390Val), in TSEN54. PCR‐RFLP and Sanger sequencing elucidated that the parents and 16 relatives were heterozygote for the novel variant.

Conclusion

We identified a novel synonymous variant, c.1170G>A, in TSEN54 associated with PCH in an Iranian family. Based on this study, we strongly suggest using “TSENopathies” to show the overlapped phenotypes among different types of PCH resulted from TSEN causative mutations.

A novel splicing mutation in SLC9A6 in a boy with Christianson syndrome

A loss of function mutation in SLC9A6 (Xq26.3) is responsible for Christianson syndrome in males. We identified a novel splicing mutation (NM_006359.2:c.1141-8C>A) of SLC9A6 in a seven-year-old boy with microcephaly, severe developmental delay, and intractable epilepsy. Functional analysis found multiple aberrant transcripts, none of which maintained the canonical open reading frame. Computer prediction tools, however, failed to detect all of the aberrant transcripts.

A novel splice mutation induces exon skipping of the EXT1 gene in patients with hereditary multiple exostoses

The molecular mechanism of hereditary multiple exostoses (HME) remains ambiguous and a limited number of studies have investigated the pathogenic mechanism of mutations in patients with HME. In the present study, a novel heterozygous splice mutation (c.1284+2del) in exostosin glycosyltransferase 1 (EXT1) gene was identified in a three-generation family with HME. Bioinformatics and TA clone-sequencing indicated that the splice site mutation would result in exon 4 skipping. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) revealed that the expression levels of wild-type EXT1/EXT2 mRNA in patients with HME were significantly decreased, compared with normal control participants (P<0.05). Abnormal EXT1 transcript lacking exon 4 (EXT1-DEL) and full-length EXT1 mRNA (EXT1-FL) were overexpressed in 293-T cells and Cos-7 cells using lentivirus infection. RT-qPCR demonstrated that the expression level of EXT1-DEL was significantly increased, compared with EXT1-FL (17.032 vs. 6.309, respectively; P<0.05). The protein encoded by EXT1-DEL was detected by western blot analysis, and the level was increased, compared with EXT1-FL protein expression. Immunofluorescence indicated that the protein encoded by EXT1-DEL was located in the cytoplasm of Cos-7 cells, which was consistent with the localization of the EXT1-FL protein. In conclusion, the present study identified a novel splice mutation that causes exon 4 skipping during mRNA splicing and causes reduced expression of EXT1/EXT2. The mutation in EXT1-DEL generated a unique peptide that is located in the cytoplasm in vitro, and it expands the mutation spectrum and provides molecular genetic evidence for a novel pathogenic mechanism of HME.

Identification and characterization of 40 novel hydroxymethylbilane synthase mutations that cause acute intermittent porphyria

Acute intermittent porphyria (AIP), an autosomal dominant disorder due to the half-normal activity of hydroxymethylbilane synthase (HMBS), is characterized by acute neurovisceral attacks that are precipitated by factors that induce heme biosynthesis. Molecular diagnosis is the most sensitive and specific diagnostic test for AIP, and importantly, it permits the identification of asymptomatic family members for genetic counseling and avoidance of precipitating factors. Here, we report the identification of 40 novel HMBS mutations, including 11 missense, four nonsense, 16 small insertions or deletions, eight consensus splice site mutations, and a complex insertion-deletion mutation in unrelated individuals with AIP. Prokaryotic expression of the missense mutations demonstrated that all mutants had ≤5% of expressed wildtype activity, except for c.1039G>C (p.A347P), which had 51% residual HMBS activity but was markedly thermolabile. Of note, the mutation c.612G>T (p.Q204H) altered the last nucleotide of exon 10, which resulted in an alternative HMBS transcript with an in-frame nine base-pair deletion at the 3'-terminus of exon 10 (encoding protein Q204HΔ3). When expressed, Q204HΔ3 and an in-frame three base-pair deletion (c.639_641delTGC) had no detectable HMBS activity. Western blot analyses and mapping of the missense mutations on the human HMBS crystal structure revealed that mutations near the active site or at the dimerization interface resulted in stably expressed proteins, while most that altered surface residues resulted in unstable proteins, presumably due to improper protein folding. These studies identified novel pathogenic HMBS mutations and expanded the molecular heterogeneity of AIP.

Splicing mutations in human genetic disorders: examples, detection, and confirmation

Precise pre-mRNA splicing, essential for appropriate protein translation, depends on the presence of consensus "cis" sequences that define exon-intron boundaries and regulatory sequences recognized by splicing machinery. Point mutations at these consensus sequences can cause improper exon and intron recognition and may result in the formation of an aberrant transcript of the mutated gene. The splicing mutation may occur in both introns and exons and disrupt existing splice sites or splicing regulatory sequences (intronic and exonic splicing silencers and enhancers), create new ones, or activate the cryptic ones. Usually such mutations result in errors during the splicing process and may lead to improper intron removal and thus cause alterations of the open reading frame. Recent research has underlined the abundance and importance of splicing mutations in the etiology of inherited diseases. The application of modern techniques allowed to identify synonymous and nonsynonymous variants as well as deep intronic mutations that affected pre-mRNA splicing. The bioinformatic algorithms can be applied as a tool to assess the possible effect of the identified changes. However, it should be underlined that the results of such tests are only predictive, and the exact effect of the specific mutation should be verified in functional studies. This article summarizes the current knowledge about the "splicing mutations" and methods that help to identify such changes in clinical diagnosis.

Activation of a cryptic 5' splice site reverses the impact of pathogenic splice site mutations in the spinal muscular atrophy gene

Spinal muscular atrophy (SMA) is caused by deletions or mutations of the Survival Motor Neuron 1 (SMN1) gene coupled with predominant skipping of SMN2 exon 7. The only approved SMA treatment is an antisense oligonucleotide that targets the intronic splicing silencer N1 (ISS-N1), located downstream of the 5' splice site (5'ss) of exon 7. Here, we describe a novel approach to exon 7 splicing modulation through activation of a cryptic 5'ss (Cr1). We discovered the activation of Cr1 in transcripts derived from SMN1 that carries a pathogenic G-to-C mutation at the first position (G1C) of intron 7. We show that Cr1-activating engineered U1 snRNAs (eU1s) have the unique ability to reprogram pre-mRNA splicing and restore exon 7 inclusion in SMN1 carrying a broad spectrum of pathogenic mutations at both the 3'ss and 5'ss of the exon 7. Employing a splicing-coupled translation reporter, we demonstrate that mRNAs generated by an eU1-induced activation of Cr1 produce full-length SMN. Our findings underscore a wider role for U1 snRNP in splicing regulation and reveal a novel approach for the restoration of SMN exon 7 inclusion for a potential therapy of SMA.

Alternative Splicing in Genetic Diseases: Improved Diagnosis and Novel Treatment Options

Alternative splicing is an important mechanism to regulate gene expression and to expand the repertoire of gene products in order to accommodate an increase in complexity of multicellular organisms. It needs to be precisely regulated, which is achieved via RNA structure, splicing factors, transcriptional regulation, and chromatin. Changes in any of these factors can lead to disease. These may include the core spliceosome, splicing enhancer/repressor sequences and their interacting proteins, the speed of transcription by RNA polymerase II, and histone modifications. While the basic principle of splicing is well understood, it is still very difficult to predict splicing outcome, due to the multiple levels of regulation. Current molecular diagnostics mainly uses Sanger sequencing of exons, or next-generation sequencing of gene panels or the whole exome. Functional analysis of potential splicing variants is scarce, and intronic variants are often not considered. This likely results in underestimation of the percentage of splicing variants. Understanding how sequence variants may affect splicing is not only crucial for confirmation of diagnosis and for genetic counseling, but also for the development of novel treatment options. These include small molecules, transsplicing, antisense oligonucleotides, and gene therapy. Here we review the current state of molecular mechanisms of splicing regulation and how deregulation can lead to human disease, diagnostics to detect splicing variants, and novel treatment options based on splicing correction.

Mutations of Pre-mRNA Splicing Regulatory Elements: Are Predictions Moving Forward to Clinical Diagnostics?

For more than three decades, researchers have known that consensus splice sites alone are not sufficient regulatory elements to provide complex splicing regulation. Other regulators, so-called splicing regulatory elements (SREs) are needed. Most importantly, their sequence variants often underlie the development of various human disorders. However, due to their variable location and high degeneracy, these regulatory sequences are also very difficult to recognize and predict. Many different approaches aiming to identify SREs have been tried, often leading to the development of in silico prediction tools. While these tools were initially expected to be helpful to identify splicing-affecting mutations in genetic diagnostics, we are still quite far from meeting this goal. In fact, most of these tools are not able to accurately discern the SRE-affecting pathological variants from those not affecting splicing. Nonetheless, several recent evaluations have given appealing results (namely for EX-SKIP, ESRseq and Hexplorer predictors). In this review, we aim to summarize the history of the different approaches to SRE prediction, and provide additional validation of these tools based on patients' clinical data. Finally, we evaluate their usefulness for diagnostic settings and discuss the challenges that have yet to be met.

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Alternative splicing of U2AF1 reveals a shared repression mechanism for duplicated exons

The auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF) facilitates branch point (BP) recognition and formation of lariat introns. The gene for the 35-kD subunit of U2AF gives rise to two protein isoforms (termed U2AF35a and U2AF35b) that are encoded by alternatively spliced exons 3 and Ab, respectively. The splicing recognition sequences of exon 3 are less favorable than exon Ab, yet U2AF35a expression is higher than U2AF35b across tissues. We show that U2AF35b repression is facilitated by weak, closely spaced BPs next to a long polypyrimidine tract of exon Ab. Each BP lacked canonical uridines at position -2 relative to the BP adenines, with efficient U2 base-pairing interactions predicted only for shifted registers reminiscent of programmed ribosomal frameshifting. The BP cluster was compensated by interactions involving unpaired cytosines in an upstream, EvoFold-predicted stem loop (termed ESL) that binds FUBP1/2. Exon Ab inclusion correlated with predicted free energies of mutant ESLs, suggesting that the ESL operates as a conserved rheostat between long inverted repeats upstream of each exon. The isoform-specific U2AF35 expression was U2AF65-dependent, required interactions between the U2AF-homology motif (UHM) and the α6 helix of U2AF35, and was fine-tuned by exon Ab/3 variants. Finally, we identify tandem homologous exons regulated by U2AF and show that their preferential responses to U2AF65-related proteins and SRSF3 are associated with unpaired pre-mRNA segments upstream of U2AF-repressed 3′ss. These results provide new insights into tissue-specific subfunctionalization of duplicated exons in vertebrate evolution and expand the repertoire of exon repression mechanisms that control alternative splicing.

IntSplice: prediction of the splicing consequences of intronic single-nucleotide variations in the human genome

Precise spatiotemporal regulation of splicing is mediated by splicing cis-elements on pre-mRNA. Single-nucleotide variations (SNVs) affecting intronic cis-elements possibly compromise splicing, but no efficient tool has been available to identify them. Following an effect-size analysis of each intronic nucleotide on annotated alternative splicing, we extracted 105 parameters that could affect the strength of the splicing signals. However, we could not generate reliable support vector regression models to predict the percent-splice-in (PSI) scores for normal human tissues. Next, we generated support vector machine (SVM) models using 110 parameters to directly differentiate pathogenic SNVs in the Human Gene Mutation Database and normal SNVs in the dbSNP database, and we obtained models with a sensitivity of 0.800±0.041 (mean and s.d.) and a specificity of 0.849±0.021. Our IntSplice models were more discriminating than SVM models that we generated with Shapiro-Senapathy score and MaxEntScan::score3ss. We applied IntSplice to a naturally occurring and nine artificial intronic mutations in RAPSN causing congenital myasthenic syndrome. IntSplice correctly predicted the splicing consequences for nine of the ten mutants. We created a web service program, IntSplice (http://www.med.nagoya-u.ac.jp/neurogenetics/IntSplice) to predict splicing-affecting SNVs at intronic positions from -50 to -3.

Genetic Testing Confirmed the Early Diagnosis of X-Linked Hypophosphatemic Rickets in a 7-Month-Old Infant

Loss-of-function mutations in the phosphate regulating gene with homologies to endopeptidases on the X-chromosome (PHEX) have been causally associated with X-linked hypophosphatemic rickets (XLHR). The early diagnosis of XLHR in infants is challenging when it is based solely on clinical features and biochemical findings. We report a 7-month-old boy with a family history of hypophosphatemic rickets., who demonstrated early clinical evidence of rickets, although serial biochemical findings could not definitively confirm rickets. A sequencing assay targeting the PHEX gene was first performed on the mother’s DNA to screen for mutations in the 5′UTR, 22 coding exons, and the exon-intron junctions. Targeted mutation analysis and mRNA studies were subsequently performed on the boys’ DNA to investigate the pathogenicity of the identified mutation. Genetic screening of the PHEX gene revealed a novel mutation, c.1080-2A>C, at the splice acceptor site in intron 9. The detection of an aberrant mRNA transcript with skipped (loss of) exon 10 establishes its pathogenicity and confirms the diagnosis of XLHR in this infant. Genetic testing of the PHEX gene resulted in early diagnosis of XLHR, thus enabling initiation of therapy and prevention of progressive rachitic changes in the infant.

Analysis of HGD Gene Mutations in Patients with Alkaptonuria from the United Kingdom: Identification of Novel Mutations

Alkaptonuria (AKU) is a rare autosomal recessive disorder with incidence ranging from 1:100,000 to 1:250,000. The disorder is caused by a deficiency of the enzyme homogentisate 1,2-dioxygenase (HGD), which results from defects in the HGD gene. This enzyme converts homogentisic acid to maleylacetoacetate and has a major role in the catabolism of phenylalanine and tyrosine. To elucidate the mutation spectrum of the HGD gene in patients with alkaptonuria from 42 patients attending the National Alkaptonuria Centre, 14 exons of the HGD gene and the intron-exon boundaries were analysed by PCR-based sequencing. A total of 34 sequence variants was observed, confirming the genetic heterogeneity of AKU. Of these mutations, 26 were missense substitutions and four splice site mutations. There were two deletions and one duplication giving rise to frame shifts and one substitution abolishing the translation termination codon (no stop). Nine of the mutations were previously unreported novel variants. Using computational approaches based on the 3D structure, these novel mutations are predicted to affect the activity of the protein complex through destabilisation of the individual protomer structure or through disruption of protomer-protomer interactions.

An ENU mutagenesis screen identifies novel and known genes involved in epigenetic processes in the mouse

Background

We have used a sensitized ENU mutagenesis screen to produce mouse lines that carry mutations in genes required for epigenetic regulation. We call these lines Modifiers of murine metastable epialleles (Mommes).

Results

We report a basic molecular and phenotypic characterization for twenty of the Momme mouse lines, and in each case we also identify the causative mutation. Three of the lines carry a mutation in a novel epigenetic modifier, Rearranged L-myc fusion (Rlf), and one gene, Rap-interacting factor 1 (Rif1), has not previously been reported to be involved in transcriptional regulation in mammals. Many of the other lines are novel alleles of known epigenetic regulators. For two genes, Rlf and Widely-interspaced zinc finger (Wiz), we describe the first mouse mutants. All of the Momme mutants show some degree of homozygous embryonic lethality, emphasizing the importance of epigenetic processes. The penetrance of lethality is incomplete in a number of cases. Similarly, abnormalities in phenotype seen in the heterozygous individuals of some lines occur with incomplete penetrance.

Conclusions

Recent advances in sequencing enhance the power of sensitized mutagenesis screens to identify the function of previously uncharacterized factors and to discover additional functions for previously characterized proteins. The observation of incomplete penetrance of phenotypes in these inbred mutant mice, at various stages of development, is of interest. Overall, the Momme collection of mouse mutants provides a valuable resource for researchers across many disciplines.

Hypophosphatemic rickets caused by a novel splice donor site mutation and activation of two cryptic splice donor sites in the PHEX gene

X-linked hypophosphatemic rickets (XLH) is the most common inherited form of rickets. XLH is caused by inactivating mutations in the PHEX gene and is transmitted as an X-linked dominant disorder. We investigated PHEX mutation in a sporadic Turkish girl with hypophosphatemic rickets. The patient was 2 years of age with a complaint of inability to walk. She had bowing of legs and growth retardation. Laboratory data showed normal calcium, low phosphate with markedly elevated ALP, and low phosphate renal tubular reabsorption. She was treated with Calcitriol 0.5 mg/kg/day and oral phosphate supplement with good response. The entire coding region of PHEX gene was sequenced from patient's peripheral leukocyte DNA and a novel 13 bp deletion at the donor splice site of exon5 was found (c.663+12del). Instead of using the donor splice site of intron 4 to splice out exon 5 and intron 5, the spliceosome utilized two nearby cryptic donor splice sites (5' splice site) to splice out intron 4, resulting in two smaller transcripts. Both of them could not translate into functional proteins due to frameshift. Her parents did not carry the mutation, indicating that this is a de novo PHEX mutation likely resulting from mutagenesis of X chromosome in paternal germ cells. We conclude that c.663+12del is a novel mutation that can activate nearby cryptic 5' splice sites. The selection of cryptic 5' splice sites adds the complexity of cell's splicing mechanisms. The current study extends the database of PHEX mutation and cryptic 5' splice sites.

Genetic variation and RNA binding proteins: tools and techniques to detect functional polymorphisms

At its most fundamental level the goal of genetics is to connect genotype to phenotype. This question is asked at a basic level evaluating the role of genes and pathways in genetic model organism. Increasingly, this question is being asked in the clinic. Genomes of individuals and populations are being sequenced and compared. The challenge often comes at the stage of analysis. The variant positions are analyzed with the hope of understanding human disease. However after a genome or exome has been sequenced, the researcher is often deluged with hundreds of potentially relevant variations. Traditionally, amino-acid changing mutations were considered the tractable class of disease-causing mutations; however, mutations that disrupt noncoding elements are the subject of growing interest. These noncoding changes are a major avenue of disease (e.g., one in three hereditary disease alleles are predicted to affect splicing). Here, we review some current practices of medical genetics, the basic theory behind biochemical binding and functional assays, and then explore technical advances in how variations that alter RNA protein recognition events are detected and studied. These advances are advances in scale-high-throughput implementations of traditional biochemical assays that are feasible to perform in any molecular biology laboratory. This chapter utilizes a case study approach to illustrate some methods for analyzing polymorphisms. The first characterizes a functional intronic SNP that deletes a high affinity PTB site using traditional low-throughput biochemical and functional assays. From here we demonstrate the utility of high-throughput splicing and spliceosome assembly assays for screening large sets of SNPs and disease alleles for allelic differences in gene expression. Finally we perform three pilot drug screens with small molecules (G418, tetracycline, and valproic acid) that illustrate how compounds that rescue specific instances of differential pre-mRNA processing can be discovered.

Antisense mediated splicing modulation for inherited metabolic diseases: challenges for delivery

In the past few years, research in targeted mutation therapies has experienced significant advances, especially in the field of rare diseases. In particular, the efficacy of antisense therapy for suppression of normal, pathogenic, or cryptic splice sites has been demonstrated in cellular and animal models and has already reached the clinical trials phase for Duchenne muscular dystrophy. In different inherited metabolic diseases, splice switching oligonucleotides (SSOs) have been used with success in patients' cells to force pseudoexon skipping or to block cryptic splice sites, in both cases recovering normal transcript and protein and correcting the enzyme deficiency. However, future in vivo studies require individual approaches for delivery depending on the gene defect involved, given the different patterns of tissue and organ expression. Herein we review the state of the art of antisense therapy targeting RNA splicing in metabolic diseases, grouped according to their expression patterns-multisystemic, hepatic, or in central nervous system (CNS)-and summarize the recent progress achieved in the field of in vivo delivery of oligonucleotides to each organ or system. Successful body-wide distribution of SSOs and preferential distribution in the liver after systemic administration have been reported in murine models for different diseases, while for CNS limited data are available, although promising results with intratechal injections have been achieved.

Missing genetic risk in neural tube defects: can exome sequencing yield an insight?

BACKGROUND

Neural tube defects (NTD) have a strong genetic component, with up to 70% of variance in human prevalence determined by heritable factors. Although the identification of causal DNA variants by sequencing candidate genes from functionally relevant pathways and model organisms has provided some success, alternative approaches are demanded.

METHODS

Next generation sequencing platforms are facilitating the production of massive amounts of sequencing data, primarily from the protein coding regions of the genome, at a faster rate and cheaper cost than has previously been possible. These platforms are permitting the identification of variants (de novo, rare, and common) that are drivers of NYTD etiology, and the cost of the approach allows for the screening of increased numbers of affected and unaffected individuals from NTD families and in simplex cases.

CONCLUSION

The next generation sequencing platforms represent a powerful tool in the armory of the genetics researcher to identify the causal genetic basis of NTDs.

Alport syndrome caused by a COL4A5 deletion and exonization of an adjacent AluY

Mutation-induced activation of splice sites in intronic repetitive sequences has contributed significantly to the evolution of exon–intron structure and genetic disease. Such events have been associated with mutations within transposable elements, most frequently in mutation hot-spots of Alus. Here, we report a case of Alu exonization resulting from a 367-nt genomic COL4A5 deletion that did not encompass any recognizable transposed element, leading to the Alport syndrome. The deletion brought to proximity the 5′ splice site of COL4A5 exon 33 and a cryptic 3′ splice site in an antisense AluY copy in intron 32. The fusion exon was depleted of purines and purine-rich splicing enhancers, but had low levels of intramolecular secondary structure, was flanked by short introns and had strong 5′ and Alu-derived 3′ splice sites, apparently compensating poor composition and context of the new exon. This case demonstrates that Alu splice sites can be activated by outlying deletions, highlighting Alu versatility in shaping the exon–intron organization and expanding the spectrum of mutational mechanisms that introduce repetitive sequences in mRNAs.

A splice mutation and mRNA decay of EXT2 provoke hereditary multiple exostoses

BACKGROUND

Hereditary multiple exostoses (HME) is an autosomal dominant disease. The classical paradigm of mutation screening seeks to relate alterations in the exostosin glycosyltransferase genes, EXT1 and EXT2, which are responsible for over 70% of HME cases. However, the pathological significance of the majority of these mutations is often unclear.

METHODS

In a Chinese family with HME, EXT1 and EXT2 genes were screened by direct sequencing. The consequence of a detected mutant was predicted by in silico analysis and confirmed by mRNA analysis. The EXT1 and EXT2 mRNA and protein levels and the HS patterns in the HME patients were compared with those in healthy controls.

RESULTS

A heterozygous transition (c.743+1G>A) in the EXT2 gene, which co-segregated with the HME phenotype in this family, was identified. The G residue at position +1 in intron 4 of EXT2 was predicted to be a 5' donor splice site. The mRNA analysis revealed an alternative transcript with a cryptic splice site 5 bp downstream of the wild-type site, which harbored a premature stop codon. However, the predicted truncated protein was not detected by western blot analysis. Decay of the mutant mRNA was shown by clone sequencing and quantification analysis. The corresponding downregulation of the EXT2 mRNA will contribute to the abnormal EXT1/EXT2 ratio and HS pattern that were detected in the patients with HME.

CONCLUSION

The heterozygous mutation c.743+1G>A in the EXT2 gene causes HME as a result of abnormal splicing, mRNA decay, and the resulting haploinsufficiency of EXT2.

PNPLA6 mutations cause Boucher-Neuhauser and Gordon Holmes syndromes as part of a broad neurodegenerative spectrum

Boucher-Neuhäuser and Gordon Holmes syndromes are clinical syndromes defined by early-onset ataxia and hypogonadism plus chorioretinal dystrophy (Boucher-Neuhäuser syndrome) or brisk reflexes (Gordon Holmes syndrome). Here we uncover the genetic basis of these two syndromes, demonstrating that both clinically distinct entities are allelic for recessive mutations in the gene PNPLA6. In five of seven Boucher-Neuhäuser syndrome/Gordon Holmes syndrome families, we identified nine rare conserved and damaging mutations by applying whole exome sequencing. Further, by dissecting the complex clinical presentation of Boucher-Neuhäuser syndrome and Gordon Holmes syndrome into its neurological system components, we set out to analyse an additional 538 exomes from families with ataxia (with and without hypogonadism), pure and complex hereditary spastic paraplegia, and Charcot-Marie-Tooth disease type 2. We identified four additional PNPLA6 mutations in spastic ataxia and hereditary spastic paraplegia families, revealing that Boucher-Neuhäuser and Gordon Holmes syndromes in fact represent phenotypic clusters on a spectrum of neurodegenerative diseases caused by mutations in PNPLA6. Structural analysis indicates that the majority of mutations falls in the C-terminal phospholipid esterase domain and likely inhibits the catalytic activity of PNPLA6, which provides the precursor for biosynthesis of the neurotransmitter acetylcholine. Our findings show that PNPLA6 influences a manifold of neuronal systems, from the retina to the cerebellum, upper and lower motor neurons and the neuroendocrine system, with damage of this protein causing an extraordinarily broad continuous spectrum of associated neurodegenerative disease.

In silico tools for splicing defect prediction: a survey from the viewpoint of end users

RNA splicing is the process during which introns are excised and exons are spliced. The precise recognition of splicing signals is critical to this process and mutations affecting splicing comprise a considerable proportion of genetic disease etiology. Analysis of RNA samples from the patient is the most straightforward and reliable method to detect splicing defects. However, currently the technical limitation prohibits its use in routine clinical practice. In silico tools that predict potential consequences of splicing mutations may be useful in daily diagnostic activities. In this review, we provide medical geneticists with some basic insights into some of the most popular in silico tools for splicing defect prediction, from the viewpoint of end-users. Bioinformaticians in relevant areas who are working on huge datasets may also benefit from this review. Specifically, we focus on those tools whose primary goal is to predict the impact of mutations within the 5′ and 3′ splicing consensus regions: the algorithms used by different tools as well as their major advantages and disadvantages are briefly introduced; the formats of their input and output are summarized; and the interpretation, evaluation, and prospection are also discussed.

A novel splice-site mutation in angiotensin I-converting enzyme (ACE) gene, c.3691+1G>A (IVS25+1G>A), causes a dramatic increase in circulating ACE through deletion of the transmembrane anchor

BACKGROUND

Angiotensin-converting enzyme (ACE) (EC 4.15.1) metabolizes many biologically active peptides and plays a key role in blood pressure regulation and vascular remodeling. Elevated ACE levels are associated with different cardiovascular and respiratory diseases.

METHODS AND RESULTS

Two Belgian families with a 8-16-fold increase in blood ACE level were incidentally identified. A novel heterozygous splice site mutation of intron 25 - IVS25+1G>A (c.3691+1G>A) - cosegregating with elevated plasma ACE was identified in both pedigrees. Messenger RNA analysis revealed that the mutation led to the retention of intron 25 and Premature Termination Codon generation. Subjects harboring the mutation were mostly normotensive, had no left ventricular hypertrophy or cardiovascular disease. The levels of renin-angiotensin-aldosterone system components in the mutated cases and wild-type controls were similar, both at baseline and after 50 mg captopril. Compared with non-affected members, quantification of ACE surface expression and shedding using flow cytometry assay of dendritic cells derived from peripheral blood monocytes of affected members, demonstrated a 50% decrease and 3-fold increase, respectively. Together with a dramatic increase in circulating ACE levels, these findings argue in favor of deletion of transmembrane anchor, leading to direct secretion of ACE out of cells.

CONCLUSIONS

We describe a novel mutation of the ACE gene associated with a major familial elevation of circulating ACE, without evidence of activation of the renin-angiotensin system, target organ damage or cardiovascular complications. These data are consistent with the hypothesis that membrane-bound ACE, rather than circulating ACE, is responsible for Angiotensin II generation and its cardiovascular consequences.

Pick one, but be quick: 5' splice sites and the problems of too many choices

Splice site selection is fundamental to pre-mRNA splicing and the expansion of genomic coding potential. 5' Splice sites (5'ss) are the critical elements at the 5' end of introns and are extremely diverse, as thousands of different sequences act as bona fide 5'ss in the human transcriptome. Most 5'ss are recognized by base-pairing with the 5' end of the U1 small nuclear RNA (snRNA). Here we review the history of research on 5'ss selection, highlighting the difficulties of establishing how base-pairing strength determines splicing outcomes. We also discuss recent work demonstrating that U1 snRNA:5'ss helices can accommodate noncanonical registers such as bulged duplexes. In addition, we describe the mechanisms by which other snRNAs, regulatory proteins, splicing enhancers, and the relative positions of alternative 5'ss contribute to selection. Moreover, we discuss mechanisms by which the recognition of numerous candidate 5'ss might lead to selection of a single 5'ss and propose that protein complexes propagate along the exon, thereby changing its physical behavior so as to affect 5'ss selection.

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.

PITX2 and FOXC1 spectrum of mutations in ocular syndromes

Anterior segment dysgenesis (ASD) encompasses a broad spectrum of developmental conditions affecting anterior ocular structures and associated with an increased risk for glaucoma. Various systemic anomalies are often observed in ASD conditions such as Axenfeld-Rieger syndrome (ARS) and De Hauwere syndrome. We report DNA sequencing and copy number analysis of PITX2 and FOXC1 in 76 patients with syndromic or isolated ASD and related conditions. PITX2 mutations and deletions were found in 24 patients with dental and/or umbilical anomalies seen in all. Seven PITX2-mutant alleles were novel including c.708_730del, the most C-terminal mutation reported to date. A second case of deletion of the distant upstream but not coding region of PITX2 was identified, highlighting the importance of this recently discovered mechanism for ARS. FOXC1 deletions were observed in four cases, three of which demonstrated hearing and/or heart defects, including a patient with De Hauwere syndrome; no nucleotide mutations in FOXC1 were identified. Review of the literature identified several other patients with 6p25 deletions and features of De Hauwere syndrome. The 1.3-Mb deletion of 6p25 presented here defines the critical region for this phenotype and includes the FOXC1, FOXF2, and FOXQ1 genes. In summary, PITX2 or FOXC1 disruptions explained 63% of ARS and 6% of other ASD in our cohort; all affected patients demonstrated additional systemic defects with PITX2 mutations showing a strong association with dental and/or umbilical anomalies and FOXC1 with heart and hearing defects. FOXC1 deletion was also found to be associated with De Hauwere syndrome.

Identification of activated cryptic 5' splice sites using structure profiles and odds measure

The activation of cryptic 5′ splice sites (5′ SSs) is often related to human hereditary diseases. The DNA-based mutation screening strategies are commonly used to recognize the cryptic 5′ SSs, because features of the local DNA sequence can influence the choice of cryptic 5′ SSs. To improve the identification of the cryptic 5′ SSs, we developed a structure-based method, named SPO (structure profiles and odds measure), which combines two parameters, the structural feature derived from hydroxyl radical cleavage pattern and odds measure, to assess the likelihood of a cryptic 5′ SS activation in competing with its paired authentic 5′ SS. Compared to the current tools for identifying activated cryptic 5′ SSs, the SPO algorithm achieves higher prediction accuracy than the other methods, including MaxEnt, MDD, Markov model, weight matrix model, Shapiro and Senapathy matrix, R_i and ΔG. In addition, the predicted ΔSPO scores from the SPO algorithm exhibited a greater degree of correlation with the strength of cryptic 5′ SS activation than that measured from the other seven methods. In conclusion, the SPO algorithm provides an optimal identification of cryptic 5′ SSs, can be applied in designing mutagenesis experiments for various splicing events and may be helpful to investigate the relationship between structural variants and human hereditary diseases.

A novel intronic mutation in SHOX causes short stature by disrupting a splice acceptor site: direct demonstration of aberrant splicing by expression of a minigene in HEK-293T cells

SHOX, the short stature homeobox-containing gene, encodes a critical regulatory protein controlling long bone growth. We examined patients in one family, identified an intronic mutation, and expressed SHOX minigenes in HEK293T cells to characterize the effect on gene splicing. We identified a novel mutation at position -3 (c.-432-3C>A;g.6120C>A) of the intron 1 splice acceptor site; three short (height Z-score -2.4 to -1.7) children were heterozygous and the father (height Z-score -3.4) was homozygous. A wild-type minigene produced alternative transcripts; one utilized the normal splice site between intron 1 and exon 2, the other a cryptic splice site in exon 2. Mutant SHOX minigene generated only the smaller transcript. The exon 2 acceptor splice site is weak; an alternative transcript is normally produced using a downstream cryptic splice site. The c.-432-3C>A mutation causes further weakening, and the cryptic splice site is preferentially utilized, resulting in SHOX deficiency and short stature.

Functional analysis of splice site mutations in the human hairless (HR) gene using a minigene assay

BACKGROUND

Congenital atrichia is a rare autosomal recessive form of isolated alopecia which is caused by mutations in the human hairless (HR) gene. Patients are born with normal hair that is shed almost completely and irreversibly during the first weeks of life.

OBJECTIVES

To investigate the molecular genetic basis of congenital atrichia in two patients, and to analyse the functional consequences of one newly identified and all seven previously identified HR splice site mutations using a minigene assay.

METHODS

Molecular analysis of the HR gene was performed by direct DNA sequencing. To analyse the functional consequences of the splice site mutations, the respective sequences were cloned into a vector which allows directed splicing. After transfection of COS7 cells, isolation of RNA and cDNA synthesis, sequencing was performed to analyse the products.

RESULTS

Two novel mutations were identified: an insertion in exon 2 (c.485insT; p.C162LfsX17), and a splice site mutation (c.2847-1G>A). In vitro analysis revealed aberrant splicing for all eight of the investigated HR splice site mutations. Comparison with the results of two biocomputational programs (neural network splice server and CRYP-SKIP) and calculation of consensus values revealed that the predictions of these two programs were consistent in only five and two of the eight mutations, respectively.

CONCLUSIONS

This is the first report to analyse the consequences of HR splice site mutations using a cell-based in vitro assay. The results highlight the importance of performing splicing experiments to clarify the consequences of putative splice site mutations.

Compensatory signals associated with the activation of human GC 5' splice sites

GC 5′ splice sites (5′ss) are present in ∼1% of human introns, but factors promoting their efficient selection are poorly understood. Here, we describe a case of X-linked agammaglobulinemia resulting from a GC 5′ss activated by a mutation in BTK intron 3. This GC 5′ss was intrinsically weak, yet it was selected in >90% primary transcripts in the presence of a strong and intact natural GT counterpart. We show that efficient selection of this GC 5′ss required a high density of GAA/CAA-containing splicing enhancers in the exonized segment and was promoted by SR proteins 9G8, Tra2β and SC35. The GC 5′ss was efficiently inhibited by splice-switching oligonucleotides targeting either the GC 5′ss itself or the enhancer. Comprehensive analysis of natural GC-AG introns and previously reported pathogenic GC 5′ss showed that their efficient activation was facilitated by higher densities of splicing enhancers and lower densities of silencers than their GT 5′ss equivalents. Removal of the GC-AG introns was promoted to a minor extent by the splice-site strength of adjacent exons and inhibited by flanking Alu repeats, with the first downstream Alus located on average at a longer distance from the GC 5′ss than other transposable elements. These results provide new insights into the splicing code that governs selection of noncanonical splice sites.

Cryptic splice sites and split genes

We describe a new program called cryptic splice finder (CSF) that can reliably identify cryptic splice sites (css), so providing a useful tool to help investigate splicing mutations in genetic disease. We report that many css are not entirely dormant and are often already active at low levels in normal genes prior to their enhancement in genetic disease. We also report a fascinating correlation between the positions of css and introns, whereby css within the exons of one species frequently match the exact position of introns in equivalent genes from another species. These results strongly indicate that many introns were inserted into css during evolution and they also imply that the splicing information that lies outside some introns can be independently recognized by the splicing machinery and was in place prior to intron insertion. This indicates that non-intronic splicing information had a key role in shaping the split structure of eukaryote genes.

A novel synonymous substitution in the GCK gene causes aberrant splicing in an Italian patient with GCK-MODY phenotype

GCK gene analysis in an Italian MODY patient revealed a novel synonymous substitution in exon 4 (c.459T>G; p.Pro153Pro) resulting in an aberrant transcript lacking the last eight codons of the same exon. Our findings emphazise the importance of not underestimating synonymous variations when screening for disease-causing mutations.

Prediction of single-nucleotide substitutions that result in exon skipping: identification of a splicing silencer in BRCA1 exon 6

Missense, nonsense, and translationally silent mutations can inactivate genes by altering the inclusion of mutant exons in mRNA, but their overall frequency among disease-causing exonic substitutions is unknown. Here, we have tested missense and silent mutations deposited in the BRCA1 mutation databases of unclassified variants for their effects on exon inclusion. Analysis of 21 BRCA1 variants using minigene assays revealed a single exon-skipping mutation c.231G>T. Comprehensive mutagenesis of an adjacent 12-nt segment showed that this silent mutation resulted in a higher level of exon skipping than the 35 other single-nucleotide substitutions. Exon inclusion levels of mutant constructs correlated significantly with predicted splicing enhancers/silencers, prompting the development of two online utilities freely available at http://www.dbass.org.uk. EX-SKIP quickly estimates which allele is more susceptible to exon skipping, whereas HOT-SKIP examines all possible mutations at each exon position and identifies candidate exon-skipping positions/substitutions. We demonstrate that the distribution of exon-skipping and disease-associated substitutions previously identified in coding regions was biased toward top-ranking HOT-SKIP mutations. Finally, we show that proteins 9G8, SC35, SF2/ASF, Tra2, and hnRNP A1 were associated with significant alterations of BRCA1 exon 6 inclusion in the mRNA. Together, these results facilitate prediction of exonic substitutions that reduce exon inclusion in mature transcripts.

Overexpression of adapted U1snRNA in patients' cells to correct a 5' splice site mutation in propionic acidemia

Splicing defects account for 16% of the mutant alleles in the PCCA and PCCB genes, encoding both subunits of the propionyl-CoA carboxylase (PCC) enzyme, defective in propionic acidemia, one of the most frequent organic acidemias causing variable neurological impairment. Most of the splicing mutations identified affect the conserved 3' splice (3' ss) or 5' splice (5' ss) sites, the latter predictably through lowering the strength of base pairing with U1snRNA. Among the 5' ss mutations we have focused on the c.1209+3A>G (IVS13+3A>G) mutation in the PCCA gene, identified in four patients (three homozygous and one heterozygous) of common geographical origin and causing exon 13 skipping. To study the potential of splicing modulation to restore PCC function, we analyzed the effect of transient transfections in patients' cells with modified U1snRNA adapted to compensate the mutant change and other mismatches at different positions of the 5' ss. Using this strategy normal transcript could be efficiently recovered with the concomitant disappearance of the aberrant exon skipping transcript, as observed after standard RT-PCR and sequence analysis or using fluorescent primers and semiquantitative RT-PCR. Different efficiencies with up to 100% exon inclusion were observed depending on the transfection conditions and specifically on the adapted U1snRNA used, confirming previously reported dependencies between nucleotides at the 5' ss for its correct recognition by the spliceosome. The reversal of the splicing defect did not result in a significant increase in enzyme activity, suggesting other factors must be taken into account for the application of overexpression of splice factors such as U1 as therapeutic strategy for splice defects.

Ex vivo splicing assays of mutations at noncanonical positions of splice sites in USHER genes

Molecular diagnosis in Usher syndrome type 1 and 2 patients led to the identification of 21 sequence variations located in noncanonical positions of splice sites in MYO7A, CDH23, USH1C, and USH2A genes. To establish experimentally the splicing pattern of these substitutions, whose impact on splicing is not always predictable by available softwares, ex vivo splicing assays were performed. The branch-point mapping strategy was also used to investigate further a putative branch-point mutation in USH2A intron 43. Aberrant splicing was demonstrated for 16 of the 21 (76.2%) tested sequence variations. The mutations resulted more frequently in activation of a nearby cryptic splice site or use of a de novo splice site than exon skipping (37.5%). This study allowed the reclassification as splicing mutations of one silent (c.7872G>A (p.Glu2624Glu) in CDH23) and four missense mutations (c.2993G>A (p.Arg998Lys) in USH2A, c.592G>A (p.Ala198Thr), c.3503G>C [p.Arg1168Pro], c.5944G>A (p.Gly1982Arg) in MYO7A), whereas it provided clues about a role in structure/function in four other cases: c.802G>A (p.Gly268Arg), c.653T>A (p.Val218Glu) (USH2A), and c.397C>T (p.His133Tyr), c.3502C>T (p.Arg1168Trp) (MYO7A). Our data provide insights into the contribution of splicing mutations in Usher genes and illustrate the need to define accurately their splicing outcome for diagnostic purposes.