RNA-sequencing of a mouse-model of spinal muscular atrophy reveals tissue-wide changes in splicing of U12-dependent introns.

Spinal Muscular Atrophy (SMA) is a neuromuscular disorder caused by insufficient levels of the Survival of Motor Neuron (SMN) protein. SMN is expressed ubiquitously and functions in RNA processing pathways that include trafficking of mRNA and assembly of snRNP complexes. Importantly, SMA severity is correlated with decreased snRNP assembly activity. In particular, the minor spliceosomal snRNPs are affected, and some U12-dependent introns have been reported to be aberrantly spliced in patient cells and animal models. SMA is characterized by loss of motor neurons, but the underlying mechanism is largely unknown. It is likely that aberrant splicing of genes expressed in motor neurons is involved in SMA pathogenesis, but increasing evidence indicates that pathologies also exist in other tissues. We present here a comprehensive RNA-seq study that covers multiple tissues in an SMA mouse model. We show elevated U12-intron retention in all examined tissues from SMA mice, and that U12-dependent intron retention is induced upon siRNA knock-down of SMN in HeLa cells. Furthermore, we show that retention of U12-dependent introns is mitigated by ASO treatment of SMA mice and that many transcriptional changes are reversed. Finally, we report on missplicing of several Ca2+ channel genes that may explain disrupted Ca2+ homeostasis in SMA and activation of Cdk5.

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.

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.

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.

Absence of an intron splicing silencer in porcine Smn1 intron 7 confers immunity to the exon skipping mutation in human SMN2.

Spinal Muscular Atrophy is caused by homozygous loss of SMN1. All patients retain at least one copy of SMN2 which produces an identical protein but at lower levels due to a silent mutation in exon 7 which results in predominant exclusion of the exon. Therapies targeting the splicing of SMN2 exon 7 have been in development for several years, and their efficacy has been measured using either in vitro cellular assays or in vivo small animal models such as mice. In this study we evaluated the potential for constructing a mini-pig animal model by introducing minimal changes in the endogenous porcine Smn1 gene to maintain the native genomic structure and regulation. We found that while a Smn2-like mutation can be introduced in the porcine Smn1 gene and can diminish the function of the ESE, it would not recapitulate the splicing pattern seen in human SMN2 due to absence of a functional ISS immediately downstream of exon 7. We investigated the ISS region and show here that the porcine ISS is inactive due to disruption of a proximal hnRNP A1 binding site, while a distal hnRNP A1 binding site remains functional but is unable to maintain the functionality of the ISS as a whole.

Identification of six novel PTH1R mutations in families with a history of primary failure of tooth eruption.

Primary Failure of tooth Eruption (PFE) is a non-syndromic disorder which can be caused by mutations in the parathyroid hormone receptor 1 gene (PTH1R). Traditionally, the disorder has been identified clinically based on post-emergent failure of eruption of permanent molars. However, patients with PTH1R mutations will not benefit from surgical and/or orthodontic treatment and it is therefore clinically important to establish whether a given failure of tooth eruption is caused by a PTH1R defect or not. We analyzed the PTH1R gene in six patients clinically diagnosed with PFE, all of which had undergone surgical and/or orthodontic interventions, and identified novel PTH1R mutations in all. Four of the six mutations were predicted to abolish correct mRNA maturation either through introduction of premature stop codons (c.947C>A and c.1082G>A), or by altering correct mRNA splicing (c.544-26_544-23del and c.989G>T). The latter was validated by transfection of minigenes. The six novel mutations expand the mutation spectrum for PFE from eight to 14 pathogenic mutations. Loss-of-function mutations in PTH1R are also associated with recessively inherited Blomstrand chondrodysplasia. We compiled all published PTH1R mutations and identified a mutational overlap between Blomstrand chondrodysplasia and PFE. The results suggest that a genetic approach to preclinical diagnosis will have important implication for surgical and orthodontic treatment of patients with failure of tooth eruption.

CUGBP1 and MBNL1 preferentially bind to 3' UTRs and facilitate mRNA decay.

CUGBP1 and MBNL1 are developmentally regulated RNA-binding proteins that are causally associated with myotonic dystrophy type 1. We globally determined the in vivo RNA-binding sites of CUGBP1 and MBNL1. Interestingly, CUGBP1 and MBNL1 are both preferentially bound to 39 UTRs. Analysis of CUGBP1- and MBNL1-bound 39 UTRs demonstrated that both factors mediate accelerated mRNA decay and temporal profiles of expression arrays supported this. Role of CUGBP1 on accelerated mRNA decay has been previously reported, but the similar function of MBNL1 has not been reported to date. It is well established that CUGBP1 and MBNL1 regulate alternative splicing. Screening by exon array and validation by RT-PCR revealed position dependence of CUGBP1- and MBNL1-binding sites on the resulting alternative splicing pattern. This study suggests that regulation of CUGBP1 and MBNL1 is essential for accurate control of destabilization of a broad spectrum of mRNAs as well as of alternative splicing events.