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.

Antisense Oligonucleotide Rescue of Deep-Intronic Variants Activating Pseudoexons in the 6-Pyruvoyl-Tetrahydropterin Synthase Gene.

We report two new 6-pyruvoyl-tetrahydropterin synthase splicing variants identified through genomic sequencing and transcript analysis in a patient with tetrahydrobiopterin deficiency, presenting with hyperphenylalaninemia and monoamine neurotransmitter deficiency. Variant c.243 + 3A>G causes exon 4 skipping. The deep-intronic c.164-672C>T variant creates a potential 5' splice site that leads to the inclusion of four overlapping pseudoexons, corresponding to exonizations of an antisense short interspersed nuclear element AluSq repeat sequence. Two of the identified pseudoexons have been reported previously, activated by different deep-intronic variants, and were also detected at residual levels in control cells. Interestingly, the predominant pseudoexon is nearly identical to a disease causing activated pseudoexon in the F8 gene, with the same 3' and 5' splice sites. Splice switching antisense oligonucleotides (SSOs) were designed to hybridize with splice sites and/or predicted binding sites for regulatory splice factors. Different SSOs corrected the aberrant pseudoexon inclusion, both in minigenes and in fibroblasts from patients carrying the new variant c.164-672C>T or the previously described c.164-716A>T. With SSO treatment PTPS protein was recovered, illustrating the therapeutic potential of the approach, for patients with different pseudoexon activating variants in the region. In addition, the natural presence of pseudoexons in the wild type context suggests the possibility of applying the antisense strategy in patients with hypomorphic PTS variants with the purpose of upregulating their expression to increase overall protein and activity.

Essential role of CK2α for the interaction and stability of replication fork factors during DNA synthesis and activation of the S-phase checkpoint.

The ataxia telangiectasia mutated and Rad3-related (ATR)-CHK1 pathway is the major signalling cascade activated in response to DNA replication stress. This pathway is associated with the core of the DNA replication machinery comprising CDC45, the replicative MCM2-7 hexamer, GINS (altogether forming the CMG complex), primase-polymerase (POLε, -α, and -δ) complex, and additional fork protection factors such as AND-1, CLASPIN (CLSPN), and TIMELESS/TIPIN. In this study, we report that functional protein kinase CK2α is critical for preserving replisome integrity and for mounting S-phase checkpoint signalling. We find that CDC45, CLSPN and MCM7 are novel CK2α interacting partners and these interactions are particularly important for maintenance of stable MCM7-CDC45, ATRIP-ATR-MCM7, and ATR-CLSPN protein complexes. Consistently, cells depleted of CK2α and treated with hydroxyurea display compromised replisome integrity, reduced chromatin binding of checkpoint mediator CLSPN, attenuated ATR-mediated S-phase checkpoint and delayed recovery of stalled forks. In further support of this, differential gene expression analysis by RNA-sequencing revealed that down-regulation of CK2α accompanies global shutdown of genes that are implicated in the S-phase checkpoint. These findings add to our understanding of the molecular mechanisms involved in DNA replication by showing that the protein kinase CK2α is essential for maintaining the stability of the replisome machinery and for optimizing ATR-CHK1 signalling activation upon replication stress.

Pseudoexon activation in disease by non-splice site deep intronic sequence variation - wild type pseudoexons constitute high-risk sites in the human genome.

Accuracy of pre-messenger RNA (pre-mRNA) splicing is crucial for normal gene expression. Complex regulation supports the spliceosomal distinction between authentic exons and the many seemingly functional splice sites delimiting pseudoexons. Pseudoexons are nonfunctional intronic sequences that can be activated for aberrant inclusion in mRNA, which may cause disease. Pseudoexon activation is very challenging to predict, in particular when activation occurs by sequence variants that alter the splicing regulatory environment without directly affecting splice sites. As pseudoexon inclusion often evades detection due to activation of nonsense-mediated mRNA decay, and because conventional diagnostic procedures miss deep intronic sequence variation, pseudoexon activation is a heavily underreported disease mechanism. Pseudoexon characteristics have mainly been studied based on in silico predicted sequences. Moreover, because recognition of sequence variants that create or strengthen splice sites is possible by comparison with well-established consensus sequences, this type of pseudoexon activation is by far the most frequently reported. Here we review all known human disease-associated pseudoexons that carry functional splice sites and are activated by deep intronic sequence variants located outside splice site sequences. We delineate common characteristics that make this type of wild type pseudoexons distinct high-risk sites in the human genome.

Vulnerable exons, like ACADM exon 5, are highly dependent on maintaining a correct balance between splicing enhancers and silencers.

It is now widely accepted that aberrant splicing of constitutive exons is often caused by mutations affecting cis-acting splicing regulatory elements, but there is a misconception that all exons have an equal dependency on splicing regulatory elements and thus a similar susceptibility to aberrant splicing. We investigated exonic mutations in ACADM exon 5 to experimentally examine their effect on splicing and found that 7 out of 11 tested mutations affected exon inclusion, demonstrating that this constitutive exon is particularly vulnerable to exonic splicing mutations. Employing ACADM exon 5 and 6 as models, we demonstrate that the balance between splicing enhancers and silencers, flanking intron length, and flanking splice site strength are important factors that determine exon definition and splicing efficiency of the exon in question. Our study shows that two constitutive exons in ACADM have different inherent vulnerabilities to exonic splicing mutations. This suggests that in silico prediction of potential pathogenic effects on splicing from exonic mutations may be improved by also considering the inherent vulnerability of the exon. Moreover, we show that single nucleotide polymorphism that affect either of two different exonic splicing silencers, located far apart in exon 5, all protect against both immediately flanking and more distant exonic splicing mutations.

Identification of SRSF10 as a regulator of SMN2 ISS-N1.

Understanding the splicing code can be challenging as several splicing factors bind to many splicing-regulatory elements. The SMN1 and SMN2 silencer element ISS-N1 is the target of the antisense oligonucleotide drug, Spinraza, which is the treatment against spinal muscular atrophy. However, limited knowledge about the nature of the splicing factors that bind to ISS-N1 and inhibit splicing exists. It is likely that the effect of Spinraza comes from blocking binding of these factors, but so far, an unbiased characterization has not been performed and only members of the hnRNP A1/A2 family have been identified by Western blot analysis and nuclear magnetic resonance to bind to this silencer. Employing an MS/MS-based approach and surface plasmon resonance imaging, we show for the first time that splicing factor SRSF10 binds to ISS-N1. Furthermore, using splice-switching oligonucleotides we modulated the splicing of the SRSF10 isoforms generating either the long or the short protein isoform of SRSF10 to regulate endogenous SMN2 exon 7 inclusion. We demonstrate that the isoforms of SRSF10 regulate SMN1 and SMN2 splicing with different strength correlating with the length of their RS domain. Our results suggest that the ratio between the SRSF10 isoforms is important for splicing regulation.

Down-regulation of CK2α correlates with decreased expression levels of DNA replication minichromosome maintenance protein complex (MCM) genes.

Protein kinase CK2 is a serine/threonine kinase composed of two catalytic subunits (CK2α and/or CK2α') and two regulatory subunits (CK2ß). It is implicated in every stage of the cell cycle and in the regulation of various intracellular pathways associated with health and disease states. The catalytic subunits have similar biochemical activity, however, their functions may differ significantly in cells and in vivo. In this regard, homozygous deletion of CK2α leads to embryonic lethality in mid-gestation potentially due to severely impaired cell proliferation. To determine the CK2α-dependent molecular mechanisms that control cell proliferation, we established a myoblast-derived cell line with inducible silencing of CK2α and carried out a comprehensive RNA-Seq analysis of gene expression. We report evidence that CK2α depletion causes delayed cell cycle progression through the S-phase and defective response to replication stress. Differential gene expression analysis revealed that the down-regulated genes were enriched in pathways implicated in cell cycle regulation, DNA replication and DNA damage repair. Interestingly, the genes coding for the minichromosome maintenance proteins (MCMs), which constitute the core of the replication origin recognition complex, were among the most significantly down-regulated genes. These findings were validated in cells and whole mouse embryos. Taken together, our study provides new evidence for a critical role of protein kinase CK2 in controlling DNA replication initiation and the expression levels of replicative DNA helicases, which ensure maintenance of proliferative potential and genome integrity in eukaryotic cells.

Next generation sequencing of RNA reveals novel targets of resveratrol with possible implications for Canavan disease.

Resveratrol (RSV) is a small compound first identified as an activator of sirtuin 1 (SIRT1), a key factor in mediating the effects of caloric restriction. Since then, RSV received great attention for its widespread beneficial effects on health and in connection to many diseases. RSV improves the metabolism and the mitochondrial function, and more recently it was shown to restore fatty acid ß-oxidation (FAO) capacities in patient fibroblasts harboring mutations with residual enzyme activity. Many of RSV's beneficial effects are mediated by the transcriptional coactivator PGC-1α, a direct target of SIRT1 and a master regulator of the mitochondrial fatty acid oxidation. Despite numerous studies RSV's mechanism of action is still not completely elucidated. Our aim was to investigate the effects of RSV on gene regulation on a wide scale, possibly to detect novel genes whose up-regulation by RSV may be of interest with respect to disease treatment. We performed Next Generation Sequencing of RNA on normal fibroblasts treated with RSV. To investigate whether the effects of RSV are mediated through SIRT1 we expanded the analysis to include SIRT1-knockdown fibroblasts. We identified the aspartoacylase (ASPA) gene, mutated in Canavan disease, to be strongly up-regulated by RSV in several cell lines, including Canavan disease fibroblasts. We further link RSV to the up-regulation of other genes involved in myelination including the glial specific transcription factors POU3F1, POU3F2, and myelin basic protein (MBP). We also observe a strong up-regulation by RSV of the riboflavin transporter gene SLC52a1. Mutations in SLC52a1 cause transient multiple acyl-CoA dehydrogenase deficiency (MADD). Our analysis of alternative splicing identified novel metabolically important genes affected by RSV, among which is particularly interesting the α subunit of the stimulatory G protein (Gsα), which regulates the cellular levels of cAMP through adenylyl cyclase. We conclude that in fibroblasts RSV stimulates the PGC-1α and p53 pathways, and up-regulates genes affecting the glucose metabolism, mitochondrial ß-oxidation, and mitochondrial biogenesis. We further confirm that RSV might be a relevant treatment in the correction of FAO deficiencies and we suggest that treatment in other metabolic disorders including Canavan disease and MADD might be also beneficial.

Blocking of an intronic splicing silencer completely rescues IKBKAP exon 20 splicing in familial dysautonomia patient cells.

Familial dysautonomia (FD) is a severe genetic disorder causing sensory and autonomic dysfunction. It is predominantly caused by a c.2204+6T>C mutation in the IKBKAP gene. This mutation decreases the 5' splice site strength of IKBKAP exon 20 leading to exon 20 skipping and decreased amounts of full-length IKAP protein. We identified a binding site for the splicing regulatory protein hnRNP A1 downstream of the IKBKAP exon 20 5'-splice site. We show that hnRNP A1 binds to this splicing regulatory element (SRE) and that two previously described inhibitory SREs inside IKBKAP exon 20 are also bound by hnRNP A1. Knockdown of hnRNP A1 in FD patient fibroblasts increases IKBKAP exon 20 inclusion demonstrating that hnRNP A1 is a negative regulator of IKBKAP exon 20 splicing. Furthermore, by mutating the SREs in an IKBKAP minigene we show that all three SREs cause hnRNP A1-mediated exon repression. We designed splice switching oligonucleotides (SSO) that blocks the intronic hnRNP A1 binding site, and demonstrate that this completely rescues splicing of IKBKAP exon 20 in FD patient fibroblasts and increases the amounts of IKAP protein. We propose that this may be developed into a potential new specific treatment of FD.

Corrigendum: Splicing factor 1 modulates dietary restriction and TORC1 pathway longevity in C. elegans.

This corrects the article DOI: 10.1038/nature20789.

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.

The ETFDH c.158A>G variation disrupts the balanced interplay of ESE- and ESS-binding proteins thereby causing missplicing and multiple Acyl-CoA dehydrogenation deficiency.

Multiple acyl-CoA dehydrogenation deficiency is a disorder of fatty acid and amino acid oxidation caused by defects of electron transfer flavoprotein (ETF) or its dehydrogenase (ETFDH). A clear relationship between genotype and phenotype makes genotyping of patients important not only diagnostically but also for prognosis and for assessment of treatment. In the present study, we show that a predicted benign ETFDH missense variation (c.158A>G/p.Lys53Arg) in exon 2 causes exon skipping and degradation of ETFDH protein in patient samples. Using splicing reporter minigenes and RNA pull-down of nuclear proteins, we show that the c.158A>G variation increases the strength of a preexisting exonic splicing silencer (ESS) motif UAGGGA. This ESS motif binds splice inhibitory hnRNP A1, hnRNP A2/B1, and hnRNP H proteins. Binding of these inhibitory proteins prevents binding of the positive splicing regulatory SRSF1 and SRSF5 proteins to nearby and overlapping exonic splicing enhancer elements and this causes exon skipping. We further suggest that binding of hnRNP proteins to UAGGGA is increased by triggering synergistic hnRNP H binding to GGG triplets located upstream and downsteam of the UAGGGA motif. A number of disease-causing exonic elements that induce exon skipping in other genes have a similar architecture as the one in ETFDH exon 2.

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

PKU is caused by mutations in PAH. A c.30C>G synonymous variation in exon 1, previously reported as neutral, was observed in two patients. The variation creates a GGG triplet, which is part of several exonic splicing silencer (ESS) motifs. Because the 5'-splice site of PAH exon 1 is intrinsically weak and therefore could be responsive to a new flanking ESS, we hypothesized that c.30C>G could cause aberrant mRNA splicing. We demonstrate that c.30C>G causes aberrant mRNA splicing in two different reporter minigenes, and that this is abolished if a preexisting flanking GGG triplet is disrupted. GGG triplets are part of the consensus motif bound by splicing-inhibitory hnRNPH proteins and we observed a dramatic increase in hnRNPH binding to c.30C>G PAH RNA. We conclude that c.30C>G creates a hnRNPH-binding ESS, which can disrupt mRNA splicing. A disease-causing mutation in HEXB, which has previously been associated with exon skipping in patients also creates a GGG triplet. We show that the mutant HEXB motif causes exon skipping of a reporter minigene and that this is also influenced by a flanking GGG triplet. We suggest that aberrant splicing caused by creation/abolishment of GGG triplets located together with a preexisting flanking GGG triplet, may be an underreported cause of human disease. It is important to recognize that exonic sequence changes may disrupt mRNA splicing. This is particularly important in PAH, since PKU patients harboring such mutations are unlikely to respond to therapy with 6R-tetrahydrobiopterin (BH(4)), despite the fact that the genetic code indicates otherwise.