The phenylalanine hydroxylase c.30C>G synonymous variation (p.G10G) creates a common exonic splicing silencer

All exons are not created equal-exon vulnerability determines the effect of exonic mutations on splicing

It is now widely accepted that aberrant splicing of constitutive exons is often caused by mutations affecting cis-acting splicing regulatory elements (SREs), but there is a misconception that all exons have an equal dependency on SREs and thus a similar vulnerability to aberrant splicing. We demonstrate that some exons are more likely to be affected by exonic splicing mutations (ESMs) due to an inherent vulnerability, which is context dependent and influenced by the strength of exon definition. We have developed VulExMap, a tool which is based on empirical data that can designate whether a constitutive exon is vulnerable. Using VulExMap, we find that only 25% of all exons can be categorized as vulnerable, whereas two-thirds of 359 previously reported ESMs in 75 disease genes are located in vulnerable exons. Because VulExMap analysis is based on empirical data on splicing of exons in their endogenous context, it includes all features important in determining the vulnerability. We believe that VulExMap will be an important tool when assessing the effect of exonic mutations by pinpointing whether they are located in exons vulnerable to ESMs.

Splice-Switching Antisense Oligonucleotides Correct Phenylalanine Hydroxylase Exon 11 Skipping Defects and Rescue Enzyme Activity in Phenylketonuria

The PAH gene encodes the hepatic enzyme phenylalanine hydroxylase (PAH), and its deficiency, known as phenylketonuria (PKU), leads to neurotoxic high levels of phenylalanine. PAH exon 11 is weakly defined, and several missense and intronic variants identified in patients affect the splicing process. Recently, we identified a novel intron 11 splicing regulatory element where U1snRNP binds, participating in exon 11 definition. In this work, we describe the implementation of an antisense strategy targeting intron 11 sequences to correct the effect of PAH mis-splicing variants. We used an in vitro assay with minigenes and identified splice-switching antisense oligonucleotides (SSOs) that correct the exon skipping defect of PAH variants c.1199+17G>A, c.1199+20G>C, c.1144T>C, and c.1066-3C>T. To examine the functional rescue induced by the SSOs, we generated a hepatoma cell model with variant c.1199+17G>A using CRISPR/Cas9. The edited cell line reproduces the exon 11 skipping pattern observed from minigenes, leading to reduced PAH protein levels and activity. SSO transfection results in an increase in exon 11 inclusion and corrects PAH deficiency. Our results provide proof of concept of the potential therapeutic use of a single SSO for different exonic and intronic splicing variants causing PAH exon 11 skipping in PKU.

In silico methods for predicting functional synonymous variants

Single nucleotide variants (SNVs) contribute to human genomic diversity. Synonymous SNVs are previously considered to be "silent," but mounting evidence has revealed that these variants can cause RNA and protein changes and are implicated in over 85 human diseases and cancers. Recent improvements in computational platforms have led to the development of numerous machine-learning tools, which can be used to advance synonymous SNV research. In this review, we discuss tools that should be used to investigate synonymous variants. We provide supportive examples from seminal studies that demonstrate how these tools have driven new discoveries of functional synonymous SNVs.

Code inside the codon: The role of synonymous mutations in regulating splicing machinery and its impact on disease

In eukaryotes, precise pre-mRNA processing, including alternative splicing, is essential to carry out the intricate protein translation process. Both point mutations (that alter the translated protein sequence) and synonymous mutations (that do not alter the translated protein sequence) are capable of affecting the splicing process. Synonymous mutations are known to affect gene expression via altering mRNA stability, mRNA secondary structure, splicing processes, and translational kinetics. In higher eukaryotes, precise splicing is regulated by three weakly conserved cis-elements, 5' and 3' splice sites and the branch site. Many other cis-acting elements (exonic/intronic splicing enhancers and silencers) and trans-acting splicing factors (serine and arginine-rich proteins and heterogeneous nuclear ribonucleoproteins) have also been found to enhance or suppress the splicing process. The appearance of synonymous mutations in cis-acting elements can alter the splicing process by changing the binding pattern of splicing factors to exonic splicing enhancers or silencer motifs. This results in exon skipping, intron retention, and various other forms of alternative splicing, eventually leading to the emergence of a wide range of diseases. The focus of this review is to elucidate the role of synonymous mutations and their impact on abnormal splicing mechanisms. Further, this study highlights the function of synonymous mutation in mediating abnormal splicing in cancer and development of X-linked, and autosomal inherited diseases.

Phenylalanine hydroxylase genotype-phenotype associations in the United States: A single center study

Phenylketonuria (PKU) is an autosomal recessive inborn error of metabolism caused by pathogenic variants in the phenylalanine hydroxylase gene (PAH). The correlation between genotype and phenotype can be complex and sometimes variable but often very useful for categorizing and predicting dietary tolerance and potential outcome. We reviewed medical records for 367 patients diagnosed with PKU or persistent mild hyperphenylalaninemia (MHP) between 1950 and 2015 who had PAH genotyping. In 351 we had the full PAH genotype as well as phenotypic characteristics such as phenylalanine (Phe) concentrations (at newborn screening, confirmation, and highest known), and dietary Phe tolerance. On 716 mutant chromosomes, including 14 in genotypes with only one identified variant, we identified 114 different pathogenic variants. The most frequent, p.R408W, was present in 15.4% of the alleles; other frequent variants were c.1315 + 1G > A (6.1%), p.I65T (5.7%), and p.R261Q (5.7%). Three variants, c.142 T > G (p.L48 V), c.615G > C (p.E205D), and c.1342_1345delCTCC, were novel. We used the phenotypic parameters of variants paired with null alleles (functional hemizygotes) to assign the variants as classic PKU, moderate PKU, mild PKU, MHP-gray zone, or MHP. We also included the phenotype association(s) for all of the full genotypes. In 103 patients, we also could assign sapropterin dihydrochloride responsiveness, which is a synthetic form of the tetrahydrobiopterin (BH4) PAH cofactor. This compilation from a single metabolic center provides further information on PAH variants in the United States and the correlations between genotype and phenotype.

Intronic PAH gene mutations cause a splicing defect by a novel mechanism involving U1snRNP binding downstream of the 5' splice site

Phenylketonuria (PKU), one of the most common inherited diseases of amino acid metabolism, is caused by mutations in the phenylalanine hydroxylase (PAH) gene. Recently, PAH exon 11 was identified as a vulnerable exon due to a weak 3’ splice site, with different exonic mutations affecting exon 11 splicing through disruption of exonic splicing regulatory elements. In this study, we report a novel intron 11 regulatory element, which is involved in exon 11 splicing, as revealed by the investigated pathogenic effect of variants c.1199+17G>A and c.1199+20G>C, identified in PKU patients. Both mutations cause exon 11 skipping in a minigene system. RNA binding assays indicate that binding of U1snRNP70 to this intronic region is disrupted, concomitant with a slightly increased binding of inhibitors hnRNPA1/2. We have investigated the effect of deletions and point mutations, as well as overexpression of adapted U1snRNA to show that this splicing regulatory motif is important for regulation of correct splicing at the natural 5’ splice site. The results indicate that U1snRNP binding downstream of the natural 5’ splice site determines efficient exon 11 splicing, thus providing a basis for development of therapeutic strategies to correct PAH exon 11 splicing mutations. In this work, we expand the functional effects of non-canonical intronic U1 snRNP binding by showing that it may enhance exon definition and that, consequently, intronic mutations may cause exon skipping by a novel mechanism, where they disrupt stimulatory U1 snRNP binding close to the 5’ splice site. Notably, our results provide further understanding of the reported therapeutic effect of exon specific U1 snRNA for splicing mutations in disease.

Splicing defects constitute a major cause of human disease. Mutations affecting conserved splicing sequences at exon-intron junctions are easily recognized as possibly pathogenic, whereas variants in exonic or intronic regions are difficult to classify without functional evidence provided by transcript analysis or in vitro analysis using minigenes. In this work, we sought out to study the pathogenicity of two novel intronic PAH variants identified in phenylketonuria patients. Both mutations resulted in exon skipping in minigenes. We demonstrate that U1snRNP70 binds to the intronic region and that this binding is abolished in the mutant sequences. Correction of the splicing defect was achieved using modified U1 snRNA perfectly complementary to each of the mutant sequences. The results extend the repertoire of natural U1 snRNP cellular functions by including its role as splicing enhancer via binding downstream of the natural 5’ splice site. In addition, our results correlate with the described therapeutic effect of modified U1snRNP for splicing mutations in different genes, thus having a significant impact in the development of specific therapies for splicing defects.

Melanoma cases demonstrate increased carrier frequency of phenylketonuria/hyperphenylalanemia mutations

Identifying novel melanoma genetic risk factors informs screening and prevention efforts. Mutations in the phenylalanine hydroxylase gene (the causative gene in phenylketonuria) lead to reduced pigmentation in untreated phenylketonuria patients, and reduced pigmentation is associated with greater melanoma risk. Therefore, we sought to characterize the relationship between phenylketonuria carrier status and melanoma risk. Using National Newborn Screening Reports, we determined the United States phenylketonuria/hyperphenylalanemia carrier frequency in Caucasians to be 1.76%. We examined three publically available melanoma datasets for germline mutations in the phenylalanine hydroxylase gene associated with classic phenylketonuria and/or hyperphenylalanemia. Mutations were identified in 29/814 melanoma patients, with a carrier frequency of 3.56%. There was a twofold enrichment (p-value = 3.4 × 10^-5 ) compared to the Caucasian frequency of hyperphenylalanemia/phenylketonuria carriers. These data demonstrate a novel association between phenylalanine hydroxylase carrier status and melanoma risk. Further, functional investigation is warranted to determine the link between phenylalanine hydroxylase mutations and melanomagenesis.

The prevalent deep intronic c. 639+919 G>A GLA mutation causes pseudoexon activation and Fabry disease by abolishing the binding of hnRNPA1 and hnRNP A2/B1 to a splicing silencer

Fabry disease is an X-linked recessive inborn disorder of the glycosphingolipid metabolism, caused by total or partial deficiency of the lysosomal α-galactosidase A enzyme due to mutations in the GLA gene. The prevalent c.639+919 G>A mutation in GLA leads to pathogenic insertion of a 57bp pseudoexon sequence from intron 4, which is responsible for the cardiac variant phenotype. In this study we investigate the splicing regulatory mechanism leading to GLA pseudoexon activation. Splicing analysis of GLA minigenes revealed that pseudoexon activation is influenced by cell-type. We demonstrate that the wild-type sequence harbors an hnRNP A1 and hnRNP A2/B1-binding exonic splicing silencer (ESS) overlapping the 5'splice site (5'ss) that prevents pseudoexon inclusion. The c.639+919 G>A mutation disrupts this ESS allowing U1 snRNP recognition of the 5'ss. We show that the wild-type GLA 5'ss motif with the ESS is also able to inhibit inclusion of an unrelated pseudoexon in the FGB gene, and that also in the FGB context inactivation of the ESS by the c.639+919 G>A mutation causes pseudoexon activation, underscoring the universal nature of the ESS. Finally, we demonstrate that splice switching oligonucleotide (SSO) mediated blocking of the pseudoexon 3'ss and 5'ss effectively restores normal GLA splicing. This indicates that SSO based splicing correction may be a therapeutic alternative in the treatment of Fabry disease.

Exploring Splicing-Switching Molecules For Seckel Syndrome Therapy

The c.2101A>G synonymous change (p.G674G) in the gene for ATR, a key player in the DNA-damage response, has been the first identified genetic cause of Seckel Syndrome (SS), an orphan disease characterized by growth and mental retardation. This mutation mainly causes exon 9 skipping, through an ill-defined mechanism. Through ATR minigene expression studies, we demonstrated that the detrimental effect of this mutation (6±1% of correct transcripts only) depends on the poor exon 9 definition (47±4% in the ATR^wt context), because the change was ineffective when the weak 5' or the 3' splice sites (ss) were strengthened (scores from 0.54 to 1) by mutagenesis. Interestingly, the exonic c.2101A nucleotide is conserved across species, and the SS-causing mutation is predicted to concurrently strengthen a Splicing Silencer (ESS) and weaken a Splicing Enhancer (ESE). Consistently, the artificial c.2101A>C change, predicted to weaken the ESE only, moderately impaired exon inclusion (28±7% of correct transcripts). The observation that an antisense oligonucleotide (AON^ATR) targeting the c.2101A position recovers exon inclusion in the mutated context supports a major role of the underlying ESS. A U1snRNA variant (U1^ATR) designed to perfectly base-pair the weak 5'ss, rescued exon inclusion (63±3%) in the ATR^SS-allele. Most importantly, upon lentivirus-mediated delivery, the U1^ATR partially rescued ATR mRNA splicing (from ~19% to ~54%) and protein (from negligible to ~6%) in embryonic fibroblasts derived from humanized ATR^SS mice. Altogether these data elucidate the molecular mechanisms of the ATR c.2101A>G mutation and identify two potential complementary RNA-based therapies for Seckel syndrome.

Global identification of hnRNP A1 binding sites for SSO-based splicing modulation

Background

Many pathogenic genetic variants have been shown to disrupt mRNA splicing. Besides splice mutations in the well-conserved splice sites, mutations in splicing regulatory elements (SREs) may deregulate splicing and cause disease. A promising therapeutic approach is to compensate for this deregulation by blocking other SREs with splice-switching oligonucleotides (SSOs). However, the location and sequence of most SREs are not well known.

Results

Here, we used individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP) to establish an in vivo binding map for the key splicing regulatory factor hnRNP A1 and to generate an hnRNP A1 consensus binding motif. We find that hnRNP A1 binding in proximal introns may be important for repressing exons. We show that inclusion of the alternative cassette exon 3 in SKA2 can be significantly increased by SSO-based treatment which blocks an iCLIP-identified hnRNP A1 binding site immediately downstream of the 5’ splice site. Because pseudoexons are well suited as models for constitutive exons which have been inactivated by pathogenic mutations in SREs, we used a pseudoexon in MTRR as a model and showed that an iCLIP-identified hnRNP A1 binding site downstream of the 5′ splice site can be blocked by SSOs to activate the exon.

Conclusions

The hnRNP A1 binding map can be used to identify potential targets for SSO-based therapy. Moreover, together with the hnRNP A1 consensus binding motif, the binding map may be used to predict whether disease-associated mutations and SNPs affect hnRNP A1 binding and eventually mRNA splicing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0279-9) contains supplementary material, which is available to authorized users.

The Splicing Efficiency of Activating HRAS Mutations Can Determine Costello Syndrome Phenotype and Frequency in Cancer

Costello syndrome (CS) may be caused by activating mutations in codon 12/13 of the HRAS proto-oncogene. HRAS p.Gly12Val mutations have the highest transforming activity, are very frequent in cancers, but very rare in CS, where they are reported to cause a severe, early lethal, phenotype. We identified an unusual, new germline p.Gly12Val mutation, c.35_36GC>TG, in a 12-year-old boy with attenuated CS. Analysis of his HRAS cDNA showed high levels of exon 2 skipping. Using wild type and mutant HRAS minigenes, we confirmed that c.35_36GC>TG results in exon 2 skipping by simultaneously disrupting the function of a critical Exonic Splicing Enhancer (ESE) and creation of an Exonic Splicing Silencer (ESS). We show that this vulnerability of HRAS exon 2 is caused by a weak 3' splice site, which makes exon 2 inclusion dependent on binding of splicing stimulatory proteins, like SRSF2, to the critical ESE. Because the majority of cancer- and CS- causing mutations are located here, they affect splicing differently. Therefore, our results also demonstrate that the phenotype in CS and somatic cancers is not only determined by the different transforming potentials of mutant HRAS proteins, but also by the efficiency of exon 2 inclusion resulting from the different HRAS mutations. Finally, we show that a splice switching oligonucleotide (SSO) that blocks access to the critical ESE causes exon 2 skipping and halts proliferation of cancer cells. This unravels a potential for development of new anti-cancer therapies based on SSO-mediated HRAS exon 2 skipping.

A smart tailor-made G-clip reporter for sensitive detection of G-triplet-containing sequences

Taking advantage of the intrinsic characteristics of G-triplet-containing sequences, a pioneering tailor-made clip-like reporter containing three-fourths of a G-quadruplex is established. The reporter can clip the G triplet in the target sequence through a recognition process to form a complete G-quadruplex structure.

Splice-shifting oligonucleotide (SSO) mediated blocking of an exonic splicing enhancer (ESE) created by the prevalent c.903+469T>C MTRR mutation corrects splicing and restores enzyme activity in patient cells

The prevalent c.903+469T>C mutation in MTRR causes the cblE type of homocystinuria by strengthening an SRSF1 binding site in an ESE leading to activation of a pseudoexon. We hypothesized that other splicing regulatory elements (SREs) are also critical for MTRR pseudoexon inclusion. We demonstrate that the MTRR pseudoexon is on the verge of being recognized and is therefore vulnerable to several point mutations that disrupt a fine-tuned balance between the different SREs. Normally, pseudoexon inclusion is suppressed by a hnRNP A1 binding exonic splicing silencer (ESS). When the c.903+469T>C mutation is present two ESEs abrogate the activity of the ESS and promote pseudoexon inclusion.

Blocking the 3′splice site or the ESEs by SSOs is effective in restoring normal splicing of minigenes and endogenous MTRR transcripts in patient cells. By employing an SSO complementary to both ESEs, we were able to rescue MTRR enzymatic activity in patient cells to approximately 50% of that in controls.

We show that several point mutations, individually, can activate a pseudoexon, illustrating that this mechanism can occur more frequently than previously expected. Moreover, we demonstrate that SSO blocking of critical ESEs is a promising strategy to treat the increasing number of activated pseudoexons.

Decoding mechanisms by which silent codon changes influence protein biogenesis and function

SCOPE

Synonymous codon usage has been a focus of investigation since the discovery of the genetic code and its redundancy. The occurrences of synonymous codons vary between species and within genes of the same genome, known as codon usage bias. Today, bioinformatics and experimental data allow us to compose a global view of the mechanisms by which the redundancy of the genetic code contributes to the complexity of biological systems from affecting survival in prokaryotes, to fine tuning the structure and function of proteins in higher eukaryotes. Studies analyzing the consequences of synonymous codon changes in different organisms have revealed that they impact nucleic acid stability, protein levels, structure and function without altering amino acid sequence. As such, synonymous mutations inevitably contribute to the pathogenesis of complex human diseases. Yet, fundamental questions remain unresolved regarding the impact of silent mutations in human disorders. In the present review we describe developments in this area concentrating on mechanisms by which synonymous mutations may affect protein function and human health.

PURPOSE

This synopsis illustrates the significance of synonymous mutations in disease pathogenesis. We review the different steps of gene expression affected by silent mutations, and assess the benefits and possible harmful effects of codon optimization applied in the development of therapeutic biologics.

PHYSIOLOGICAL AND MEDICAL RELEVANCE

Understanding mechanisms by which synonymous mutations contribute to complex diseases such as cancer, neurodegeneration and genetic disorders, including the limitations of codon-optimized biologics, provides insight concerning interpretation of silent variants and future molecular therapies.

Exposing synonymous mutations

Synonymous codon changes, which do not alter protein sequence, were previously thought to have no functional consequence. Although this concept has been overturned in recent years, there is no unique mechanism by which these changes exert biological effects. A large repertoire of both experimental and bioinformatic methods has been developed to understand the effects of synonymous variants. Results from this body of work have provided global insights into how biological systems exploit the degeneracy of the genetic code to control gene expression, protein folding efficiency, and the coordinated expression of functionally related gene families. Although it is now clear that synonymous variants are important in a variety of contexts, from human disease to the safety and efficacy of therapeutic proteins, there is no clear consensus on the approaches to identify and validate these changes. Here, we review the diverse methods to understand the effects of synonymous mutations.

Splice, insertion-deletion and nonsense mutations that perturb the phenylalanine hydroxylase transcript cause phenylketonuria in India

Phenylketonuria (PKU) is an autosomal recessive metabolic disorder caused by mutational inactivation of the phenylalanine hydroxylase (PAH) gene. Missense mutations are the most common PAH mutation type detected in PKU patients worldwide. We performed PAH mutation analysis in 27 suspected Indian PKU families (including 7 from our previous study) followed by structure and function analysis of specific missense and splice/insertion-deletion/nonsense mutations, respectively. Of the 27 families, disease-causing mutations were detected in 25. A total of 20 different mutations were identified of which 7 "unique" mutations accounted for 13 of 25 mutation positive families. The unique mutations detected exclusively in Indian PKU patients included three recurrent mutations detected in three families each. The 20 mutations included only 5 missense mutations in addition to 5 splice, 4 each nonsense and insertion-deletion mutations, a silent variant in coding region and a 3'UTR mutation. One deletion and two nonsense mutations were characterized to confirm significant reduction in mutant transcript levels possibly through activation of nonsense mediated decay. All missense mutations affected conserved amino acid residues and sequence and structure analysis suggested significant perturbations in the enzyme activity of respective mutant proteins. This is probably the first report of identification of a significantly low proportion of missense PAH mutations from PKU families and together with the presence of a high proportion of splice, insertion-deletion, and nonsense mutations, points to a unique PAH mutation profile in Indian PKU patients.

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

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

SMN2 exon 7 splicing is inhibited by binding of hnRNP A1 to a common ESS motif that spans the 3' splice site

Spinal Muscular Atrophy is caused by homozygous loss of SMN1 with phenotypic modulation by SMN2. SMN2 expresses only limited amounts of full-length transcript due to skipping of exon 7 caused by disruption of an SF2/ASF binding ESE. Additionally, hnRNP A1 has been reported to inhibit inclusion of SMN2 exon 7. We previously reported high similarity between the sequence spanning the 3' ss of SMN1 and SMN2 exon 7 and an hnRNP A1 binding ESS, which regulates MCAD exon 5 splicing. We show here that this 3' ss motif indeed functions as a crucial hnRNP A1 binding ESS, which inhibits inclusion of SMN1/2 exon 7 and is antagonized by the SMN1 ESE, but not by the inactive SMN2 sequence. Pull-down experiments revealed a specific interaction between hnRNP A1 and the 3' ss AG-dinucleotide, which could be disrupted by mutations shown to improve splicing in reporter minigenes. Genomic analyses revealed that in the human genome, 3' ss matching the SMN1/2 ESS motif region are much less abundant than 3' ss with a disrupted ESS motif. This indicates that this ESS may be a general splicing inhibitory motif, which binds hnRNP A1 and inhibits exon inclusion by binding to 3' ss harboring this ESS motif.