Functional characterization of two novel splicing mutations of glucokinase gene associated with maturity-onset diabetes of the young type 2 (MODY2)

Two novel mutations in the glucokinase gene (GCK) have been identified in patients with maturity-onset diabetes of the young type-2 (MODY2), i.e., a C-for-G substitution at position ?1 of the acceptor splice site of intron 7 (c. 864-1G>C) and a synonymous c.666C>G substitution (GTC>GTG, p.V222V) at exon 6. An analysis of the splicing products obtained upon the transfection of human embryonic HEK293 cells with GCK minigene constructs carrying these mutations showed that both substitutions impaired normal splicing. As a result of c.864-1G>C, the usage of the normal acceptor site was blocked, which activated cryptic acceptor splice sites within intron 7 and generated several aberrant RNAs containing fragments of intron 7. The synonymous substitution c.666C>G created a novel donor splice site in exon 6, which results in the formation of an abnormal GCK mRNA with a 16-nucleotide deletion in exon 6. In vitro experiments on minigene splicing confirmed the inactivating effect of these mutations on glucokinase gene expression.     

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.     

Exonic Splicing Mutations Are More Prevalent than Currently Estimated and Can Be Predicted by Using In Silico Tools

The identification of a causal mutation is essential for molecular diagnosis and clinical management of many genetic disorders. However, even if next-generation exome sequencing has greatly improved the detection of nucleotide changes, the biological interpretation of most exonic variants remains challenging. Moreover, particular attention is typically given to protein-coding changes often neglecting the potential impact of exonic variants on RNA splicing. Here, we used the exon 10 of MLH1, a gene implicated in hereditary cancer, as a model system to assess the prevalence of RNA splicing mutations among all single-nucleotide variants identified in a given exon. We performed comprehensive minigene assays and analyzed patient’s RNA when available. Our study revealed a staggering number of splicing mutations in MLH1 exon 10 (77% of the 22 analyzed variants), including mutations directly affecting splice sites and, particularly, mutations altering potential splicing regulatory elements (ESRs). We then used this thoroughly characterized dataset, together with experimental data derived from previous studies on BRCA1, BRCA2, CFTR and NF1, to evaluate the predictive power of 3 in silico approaches recently described as promising tools for pinpointing ESR-mutations. Our results indicate that ΔtESRseq and ΔHZ_EI-based approaches not only discriminate which variants affect splicing, but also predict the direction and severity of the induced splicing defects. In contrast, the ΔΨ-based approach did not show a compelling predictive power. Our data indicates that exonic splicing mutations are more prevalent than currently appreciated and that they can now be predicted by using bioinformatics methods. These findings have implications for all genetically-caused diseases.     

Vitamin A Metabolite,All-trans-retinoic Acid,Mediates Alternative Splicing of Protein Kinase C δVIII (PKCδVIII) Isoform via Splicing Factor SC35

Vitamin A metabolite, all-trans-retinoic acid (RA), induces cell growth, differentiation, and apoptosis and has an emerging role in gene regulation and alternative splicing events. Protein kinase Cδ (PKCδ), a serine/threonine kinase, has a role in cell proliferation, differentiation, and apoptosis. We reported an alternatively spliced variant of human PKCδ, PKCδVIII that functions as a pro-survival protein (1). RA regulates the splicing and expression of PKCδVIII via utilization of a downstream 5′ splice site of exon 10 on PKCδ pre-mRNA. Here, we further elucidate the molecular mechanisms involved in RA regulation of alternative splicing of PKCδVIII mRNA. Overexpression and knockdown of the splicing factor SC35 (i.e. SRp30b) indicated that it is involved in PKCδVIII alternative splicing. To identify the cis-elements involved in 5′ splice site selection we cloned a minigene, which included PKCδ exon 10 and its flanking introns in the pSPL3 splicing vector. Alternative 5′ splice site utilization in the minigene was promoted by RA. Further, co-transfection of SC35 with PKCδ minigene promoted selection of 5′ splice site II. Mutation of the SC35 binding site in the PKCδ minigene abolished RA-mediated utilization of 5′ splice splice II. RNA binding assays demonstrated that the enhancer element downstream of PKCδ exon 10 is a SC35 cis-element. We conclude that SC35 is pivotal in RA-mediated PKCδ pre-mRNA alternative splicing. This study demonstrates how a nutrient, vitamin A, via its metabolite RA, regulates alternative splicing and thereby gene expression of the pro-survival protein PKCδVIII.     

Identification of a Novel Synonymous Mutation in the Human β -Ureidopropionase Gene UPB1 Affecting Pre-mRNA Splicing

β-Ureidopropionase is the third enzyme of the pyrimidine degradation pathway and it catalyzes the conversion of N-carbamyl-β-alanine and N-carbamyl-β-aminoisobutyric acid to β-alanine and β-aminoisobutyric acid, respectively, and ammonia and CO₂. To date, only 16 genetically confirmed patients with a complete ß-ureidopropionase deficiency have been reported. Here, we report the clinical, biochemical, and molecular analysis of a newly identified patient with β-ureidopropionase deficiency. Mutation analysis of the UPB1 gene showed that the patient was compound heterozygous for a novel synonymous mutation c.93C >T (p.Gly31Gly) in exon 1 and a previously described missense mutation c.977G >A (p.Arg326Gln) in exon 9. The in silico predicted effect of the synonymous mutation p.Gly31Gly on pre-mRNA splicing was investigated using a minigene approach. Wild-type and the mutated minigene constructs, containing the entire exon 1, intron 1, and exon 2 of UPB1, yielded different splicing products after expression in HEK293 cells. The c.93C >T (p.Gly31Gly) mutation resulted in altered pre-mRNA splicing of the UPB1 minigene construct and a deletion of the last 13 nucleotides of exon 1. This deletion (r.92_104delGCAAGGAACTCAG) results in a frame shift and the generation of a premature stop codon (p.Lys32SerfsX31). Using a minigene approach, we have thus identified the first synonymous mutation in the UPB1 gene, creating a cryptic splice-donor site affecting pre-mRNA splicing.     

A Novel Zinc Finger Gene ZNF468 with Two Co-Expressional Splice Variants, ZNF468.1 and ZNF468.2

A novel human zinc finger protein encoding gene ZNF468 was obtained from a fetal brain cDNA library. By BLAST-N analysis we found two different splice variants. We termed the two splice variants ZNF468.1 and ZNF468.2. By BLAST search against the human genome database, ZNF468 was mapped to 19q13.4. The ZNF468.1 cDNA has four exons, and the ZNF468.2 cDNA has one more, between the third and fourth exon. This extra exon creates a difference between the deduced protein N-termini of the two splice variants. The ZNF468.1 cDNA is 3906 bp in length, encoding a 522a a protein, and ZNF468.2 is 4024 bp, encoding a 469-aa-protein. Both proteins contain 11 C2H2-type zinc finger motifs at their C-termini. The N-terminus of the deduced protein of ZNF468.1 has a well-conserved Krüppel-associated box (KRAB) domain that consists of KRAB boxes A and B, whereas the protein of ZNF468.2 does not have the {KRAB} domain. Tissue distribution of the ZNF468 gene indicates that the two splice variants are widely expressed in normal human tissues, except in heart and brain, and they are also co-expressional.     

Cis requirements for alternative splicing of the cardiac troponin T pre-mRNA.

The cardiac troponin T (cTNT) pre-mRNA splices 17 exons contiguously but alternatively splices (includes or excludes) the fifth exon. Because both alternative splice products are processed from the same pre-mRNA species, the cTNT pre-mRNA must contain cis-acting sequences which specify exon 5 as an alternative exon. A cTNT minigene (SM-1) transfected into cultured cells produces mRNAs both including and excluding exon 5. The junctions of exons 4-5-6 and 4-6 in the cTNT minigene mRNAs are identical to those of endogenous cTNT mRNAs and no other exons are alternatively spliced. Thus, the SM-1 pre-mRNA is correctly alternatively spliced in transfected cells. To circumscribe the pre-mRNA regions which are required for the alternative nature of exon 5, we have constructed a systematic series of deletion mutants of SM-1. Transfection of this series demonstrates that a 1200 nt pre-mRNA region containing exons 4, 5, and 6 is sufficient to direct alternative splicing of exon 5. Within this region are two relatively large inverted repeats which potentially sequester the alternative exon via intramolecular base-pairing. Such sequestration of an alternative exon is consistent with models which propose pre-mRNA conformation as being determinative for alternative splicing of some pre-mRNAs. However, deletion mutants which remove the majority of each of the inverted repeats retain the ability to alternatively splice exon 5 demonstrating that neither is required for cTNT alternative splice site selection. Taken together, deletion analysis has limited cis elements required for alternative splicing to three small regions of the pre-mRNA containing exons 4, 5, and 6. In addition, the cTNT minigene pre-mRNA expresses both alternative splice products in a wide variety of cultured non-muscle cells as well as in cultured striated muscle cells, although expression of the cTNT pre-mRNA is normally restricted to striated muscle. This indicates that cis elements involved in defining the cTNT exon 5 as an alternative exon do not require muscle-specific factors in trans to function.     

Tissue-specific alternative splicing of rat brain fructose 6-phosphate 2-kinase/fructose 2,6-bisphosphatase

We have reported the occurrence of eight splice variants of rat brain fructose 6-phosphate 2-kinase/fructose 2,6-bisphosphatase (RB2K). In the present study, we quantified these splice variants in various tissues using a RNAse protection assay and found a tissue-specific pattern of alternative splicing of the RB2K gene. Splice variants containing exon F were specifically expressed in brain. Moreover, exons D and E were spliced in brain, skeletal muscle and heart. Consequently, eight, six, four and two splice variants were expressed in brain, skeletal muscle, heart and liver plus testis, respectively. These results suggest that distinct RB2K isoforms could be involved in regulation of glycolysis in a tissue-specific manner.     

Control of hnRNP A1 alternative splicing: an intron element represses use of the common 3' splice site

Alternative splicing of exon 7B in the hnRNP A1 pre-mRNA produces mRNAs encoding two proteins: hnRNP A1 and the less abundant A1B. We have reported the identification of several intron elements that contribute to exon 7B skipping. In this study, we report the activity of a novel element, conserved element 9 (CE9), located in the intron downstream of exon 7B. We show that multiple copies of CE9 inhibit exon 7B-exon 8 splicing in vitro. When CE9 is inserted between two competing 3' splice sites, a single copy of CE9 decreases splicing to the distal 3' splice site. Our in vivo results also support the conclusion that CE9 is a splicing modulator. First, inserting multiple copies of CE9 into an A1 minigene compromises the production of fully spliced products. Second, one copy of CE9 stimulates the inclusion of a short internal exon in a derivative of the human beta-globin gene. In this case, in vitro splicing assays suggest that CE9 decreases splicing of intron 1, an event that improves splicing of intron 2 and decreases skipping of the short internal exon. The ability of CE9 to act on heterologous substrates, combined with the results of a competition assay, suggest that the activity of CE9 is mediated by a trans-acting factor. Our results indicate that CE9 represses the use of the common 3' splice site in the hnRNP A1 alternative splicing unit.     

Functional analysis by minigene assay of putative splicing variants found in Bardet–Biedl syndrome patients

Bardet–Biedl syndrome (BBS) and Alström syndrome (ALMS) are rare diseases belonging to the group of ciliopathies. Although mutational screening studies of BBS/ALMS cohorts have been extensively reported, little is known about the functional effect of those changes. Thus, splicing variants are estimated to represent 15% of disease‐causing mutations, and there is growing evidence that many exonic changes are really splicing variants misclassified. In this study, we aimed to analyse for the first time several variants in BBS2, ARL6/BBS3, BBS4 and ALMS1 genes predicted to produce aberrant splicing by minigene assay. We found discordance between bioinformatics analysis and experimental data when comparing wild‐type and mutant constructs. Remarkably, we identified nonsense variants presumably resistant to nonsense‐mediated decay, even when a premature termination codon would be introduced in the second amino acid (p.(G2*) mutation in ARL6/BBS3 gene). As a whole, we report one of the first functional studies of BBS/ALMS1 variants using minigene assay, trying to elucidate their role in disease. Functional studies of variants identified in BBS and ALMS patients are essential for their proper classification and subsequent genetic counselling and could also be the start point for new therapeutic approaches, currently based only on symptomatic treatment.     

Novel alternative splicing of GABA receptor RDL exon 9 from Laodelphax striatellus modulates agonist potency

The resistance to dieldrin gene (RDL) encodes the primary subunit of the insect ionotropic γ-aminobutyric acid (GABA) receptor (GABAR), which is the target of phenylpyrazole and isoxazoline insecticides. The splice variants in exons 3 and 6 of RDL, which have been widely explored in many insects, modulate the agonist potency of the homomeric RDL GABAR and potentially play an important role in the development of insects. In the present study, four splice variants of exon 9 were identified in RDL of the small brown planthopper, Laodelphax striatellus (LsRDL), resulting in LsRDL-9a, LsRDL-9a’, LsRDL-9b, and LsRDL-9c. LsRDL-9a has one more amino acid (E, glutamic acid) compared with LsRDL-9a’, and LsRDL-9b lacked two amino acids and had seven different amino acids compared with LsRDL-9c. Two-electrode voltage-clamp recording on LsRDLs expressed in Xenopus oocytes showed that alternative splicing of exon 9 has significant impact on LsRDL sensitivity to the agonists GABA and β-alanine, whereas no significant difference was observed in the potencies of the non-competitive antagonists (NCAs) ethiprole and fluralaner on the splice variants. Our results suggest that alternative splicing of RDL exon 9 broadens functional capabilities of the GABAR in L. striatellus by influencing the action of GABA.     

Splicing of human chloride channel 1

Expression of chloride channel 1 (CLCN1/ClC-1) in skeletal muscle is driven by alternative splicing, a process regulated in part by RNA-binding protein families MBNL and CELF. Aberrant splicing of CLCN1 produces many mRNAs, which were translated into inactive proteins, resulting in myotonia in myotonic dystrophy (DM), a genetic disorder caused by the expansion of a CTG or CCTG repeat. This increase in abnormal splicing variants containing exons 6B, 7A or the insertion of a TAG stop codon just before exon 7 leads to a decrease in expression of the normal splice pattern. The majority of studies examining splicing in CLCN1 have been performed using mouse Clcn1, as have investigations into the activation and suppression of normal splicing variant expression by MBNL1-3 and CELF3–6, respectively. In contrast, examinations of human CLCN1 have been less common due to the greater complexity of splicing patterns. Here, we constructed a minigene containing CLCN1 exons 5–7 and established a novel assay system to quantify the expression of the normal splicing variant of CLCN1 using real-time RT-PCR. Antisense oligonucleotides could promote normal CLCN1 alternative splicing but the effective sequence was different from that of Clcn1. This result differs from previous reports using Clcn1, highlighting the effect of differences in splicing patterns between mice and humans.