Identification of functional single nucleotide polymorphisms in the branchpoint site

Background

The human genome contains millions of single nucleotide polymorphisms (SNPs); many of these SNPs are intronic and have unknown functional significance. SNPs occurring within intron branchpoint sites, especially at the adenine (A), would presumably affect splicing; however, this has not been systematically studied. We employed a splicing prediction tool to identify human intron branchpoint sites and screened dbSNP for identifying SNPs located in the predicted sites to generate a genome-wide branchpoint site SNP database.

Results

We identified 600 SNPs located within branchpoint sites; among which, 216 showed a change in A. After scoring the SNPs by counting the As in the ±?10 nucleotide region, only four SNPs were identified without additional As (rs13296170, rs12769205, rs75434223, and rs67785924). Using minigene constructs, we examined the effects of these SNPs on splicing. The three SNPs (rs13296170, rs12769205, and rs75434223) with nucleotide substitution at the A position resulted in abnormal splicing (exon skipping and/or intron inclusion). However, rs67785924, a 5-bp deletion that abolished the branchpoint A nucleotide, exhibited normal RNA splicing pattern, presumably using two of the downstream As as alternative branchpoints. The influence of additional As on splicing was further confirmed by studying rs2733532, which contains three additional As in the ±?10 nucleotide region.

Conclusions

We generated a high-confidence genome-wide branchpoint site SNP database, experimentally verified the importance of A in the branchpoint, and suggested that other nearby As can protect branchpoint A substitution from abnormal splicing.

     

Role of six single nucleotide polymorphisms,risk factors in coronary disease,in OLR1 alternative splicing

The OLR1 gene encodes the oxidized low-density lipoprotein receptor (LOX-1), which is responsible for the cellular uptake of oxidized LDL (Ox-LDL), foam cell formation in atheroma plaques and atherosclerotic plaque rupture. Alternative splicing (AS) of OLR1 exon 5 generates two protein isoforms with antagonistic functions in Ox-LDL uptake. Previous work identified six single nucleotide polymorphisms (SNPs) in linkage disequilibrium that influence the inclusion levels of OLR1 exon 5 and correlate with the risk of cardiovascular disease. Here we use minigenes to recapitulate the effects of two allelic series (Low- and High-Risk) on OLR1 AS and identify one SNP in intron 4 (rs3736234) as the main contributor to the differences in exon 5 inclusion, while the other SNPs in the allelic series attenuate the drastic effects of this key SNP. Bioinformatic, proteomic, mutational and functional high-throughput analyses allowed us to define regulatory sequence motifs and identify SR protein family members (SRSF1, SRSF2) and HMGA1 as factors involved in the regulation of OLR1 AS. Our results suggest that antagonism between SRSF1 and SRSF2/HMGA1, and differential recognition of their regulatory motifs depending on the identity of the rs3736234 polymorphism, influence OLR1 exon 5 inclusion and the efficiency of Ox-LDL uptake, with potential implications for atherosclerosis and coronary disease.     

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.     

Resequencing Three Candidate Genes for Major Depressive Disorder in a Dutch Cohort

Major depressive disorder (MDD) is a psychiatric disorder, characterized by periods of low mood of more than two weeks, loss of interest in normally enjoyable activities and behavioral changes. MDD is a complex disorder and does not have a single genetic cause. In 2009 a genome wide association study (GWAS) was performed on the Dutch GAIN-MDD cohort. Many of the top signals of this GWAS mapped to a region spanning the gene PCLO, and the non-synonymous coding single nucleotide polymorphism (SNP) rs2522833 in the PCLO gene became genome wide significant after post-hoc analysis. We performed resequencing of PCLO, GRM7, and SLC6A4 in 50 control samples from the GAIN-MDD cohort, to detect new genomic variants. Subsequently, we genotyped these variants in the entire GAIN-MDD cohort and performed association analysis to investigate if rs2522833 is the causal variant or simply in linkage disequilibrium with a more associated variant. GRM7 and SLC6A4 are both candidate genes for MDD from literature. We aimed to gather more evidence that rs2522833 is indeed the causal variant in the GAIN-MDD cohort or to find a previously undetected common variant in either PCLO, GRM7, or SLC6A4 with a higher association in this cohort. After next generation sequencing and association analysis we excluded the possibility of an undetected common variant to be more associated. For neither PCLO nor GRM7 we found a more associated variant. For SLC6A4, we found a new SNP that showed a lower P-value (P = 0.07) than in the GAIN-MDD GWAS (P = 0.09). However, no evidence for genome-wide significance was found. Although we did not take into account rare variants, we conclude that our results provide further support for the hypothesis that the non-synonymous coding SNP rs2522833 in the PCLO gene is indeed likely to be the causal variant in the GAIN-MDD cohort.     

Functional Analysis of Mutations in Exon 9 of NF1 Reveals the Presence of Several Elements Regulating Splicing

Neurofibromatosis type 1 (NF1) is one of the most common human hereditary disorders, predisposing individuals to the development of benign and malignant tumors in the nervous system, as well as other clinical manifestations. NF1 is caused by heterozygous mutations in the NF1 gene and around 25% of the pathogenic changes affect pre-mRNA splicing. Since the molecular mechanisms affected by these mutations are poorly understood, we have analyzed the splicing mutations identified in exon 9 of NF1, which is particularly prone to such changes, to better define the possible splicing regulatory elements. Using a minigene approach, we studied the effect of five splicing mutations in this exon described in patients. These highlighted three regulatory motifs within the exon. An in vivo splicing analysis of an extensive collection of changes generated in the minigene demonstrated that the CG motif at c.910-911 is critical for the recognition of exon 9. We also found that the GC motif at c.945-946 is involved in exon recognition through SRSF2 and that this motif is part of a Composite Exon Splicing Regulatory Element made up of physically overlapping enhancer and silencer elements. Finally, through an in vivo splicing analysis and in vitro binding assays, we demonstrated that the c.1007G>A mutation creates an Exonic Splicing Silencer element that binds the hnRNPA1 protein. The complexity of the splicing regulatory elements present in exon 9 is most likely responsible for the fact that mutations in this region represent 25% of all exonic changes that affect splicing in the NF1 gene.     

MBNL and CELF proteins regulate alternative splicing of the skeletal muscle chloride channel CLCN1

The expression and function of the skeletal muscle chloride channel CLCN1/ClC-1 is regulated by alternative splicing. Inclusion of the CLCN1 exon 7A is aberrantly elevated in myotonic dystrophy (DM), a genetic disorder caused by the expansion of a CTG or CCTG repeat. Increased exon 7A inclusion leads to a reduction in CLCN1 function, which can be causative of myotonia. Two RNA-binding protein families—muscleblind-like (MBNL) and CUG-BP and ETR-3-like factor (CELF) proteins—are thought to mediate the splicing misregulation in DM. Here, we have identified multiple factors that regulate the alternative splicing of a mouse Clcn1 minigene. The inclusion of exon 7A was repressed by MBNL proteins while promoted by an expanded CUG repeat or CELF4, but not by CUG-BP. Mutation analyses suggested that exon 7A and its flanking region mediate the effect of MBNL1, whereas another distinct region in intron 6 mediates that of CELF4. An exonic splicing enhancer essential for the inclusion of exon 7A was identified at the 5′ end of this exon, which might be inhibited by MBNL1. Collectively, these results provide a mechanistic model for the regulation of Clcn1 splicing, and reveal novel regulatory properties of MBNL and CELF proteins.     

Position-dependent effects of hnRNP A1/A2 in SMN1/2 exon7 splicing

Heterogeneous nuclear ribonucleoprotein A1 and A2 (hnRNP A1/2) is a ubiquitously expressed RNA binding protein known to bind intronic or exonic splicing silencer. Binding of hnRNP A1/2 to survival of motor neuron gene (SMN1/2) exon 7 and flanking sequences strongly inhibits the inclusion of exon 7, which causes spinal muscular atrophy, a common genetic disorder. However, the role of hnRNP A1/2 on the side away from exon 7 is unclear. Here using antisense oligonucleotides, we fished an intronic splicing enhancer (ISE) near the 3′-splice site (SS) of intron 7 of SMN1/2. Mutagenesis identified the efficient motif of the ISE as “UAGUAGG”, coupled with RNA pull down and protein overexpression, we proved that hnRNP A1/2 binding to the ISE promotes the inclusion of SMN1/2 exon 7. Using MS2-tethering array and “UAGGGU” motif walking, we further uncovered that effects of hnRNP A1/2 on SMN1/2 exon 7 splicing are position-dependent: exon 7 inclusion is inhibited when hnRNP A1/2 binds proximal to the 5′SS of intron 7, promoted when its binds proximal to the 3′SS. These data provide new insights into the splicing regulatory mechanism of SMN1/2.     

An Alu-derived intronic splicing enhancer facilitates intronic processing and modulates aberrant splicing in ATM

We have previously reported a natural GTAA deletion within an intronic splicing processing element (ISPE) of the ataxia telangiectasia mutated (ATM) gene that disrupts a non-canonical U1 snRNP interaction and activates the excision of the upstream portion of the intron. The resulting pre-mRNA splicing intermediate is then processed to a cryptic exon, whose aberrant inclusion in the final mRNA is responsible for ataxia telangiectasia. We show here that the last 40 bases of a downstream intronic antisense Alu repeat are required for the activation of the cryptic exon by the ISPE deletion. Evaluation of the pre-mRNA splicing intermediate by a hybrid minigene assay indicates that the identified intronic splicing enhancer represents a novel class of enhancers that facilitates processing of splicing intermediates possibly by recruiting U1 snRNP to defective donor sites. In the absence of this element, the splicing intermediate accumulates and is not further processed to generate the cryptic exon. Our results indicate that Alu-derived sequences can provide intronic splicing regulatory elements that facilitate pre-mRNA processing and potentially affect the severity of disease-causing splicing mutations.     

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.     

Antagonistic factors control the unproductive splicing of SC35 terminal intron

Alternative splicing is regulated in part by variations in the relative concentrations of a variety of factors, including serine/arginine-rich (SR) proteins. The SR protein SC35 self-regulates its expression by stimulating unproductive splicing events in the 3′ untranslated region of its own pre-mRNA. Using various minigene constructs containing the terminal retained intron and flanking exons, we identified in the highly conserved last exon a number of exonic splicing enhancer elements responding specifically to SC35, and showed an inverse correlation between affinity of SC35 and enhancer strength. The enhancer region, which is included in a long stem loop, also contains repressor elements, and is recognized by other RNA-binding proteins, notably hnRNP H protein and TAR DNA binding protein (TDP-43). Finally, in vitro and in cellulo experiments indicated that hnRNP H and TDP-43 antagonize the binding of SC35 to the terminal exon and specifically repress the use of SC35 terminal 3′ splice site. Our study provides new information about the molecular mechanisms of SC35-mediated splicing activation. It also highlights the existence of a complex network of self- and cross-regulatory mechanisms between splicing regulators, which controls their homeostasis and offers many ways of modulating their concentration in response to the cellular environment.     

Defective pre-mRNA splicing in PKD1 due to presumed missense and synonymous mutations causing autosomal dominant polycystic disease

Autosomal dominant polycystic kidney disease is the most common human monogenic disorder and is caused by mutations in the PKD1 or PKD2 genes. Most patients with the disease present mutations in PKD1, and a considerable number of these alterations are single base substitutions within the coding sequence that are usually predicted to lead to missense or synonymous mutations. There is growing evidence that some of these mutations can be detrimental by affecting the pre-mRNA splicing process. The aim of our study was to test PKD1 mutations, described as missense or synonymous in the literature or databases, for their effects on exon inclusion. Bioinformatics tools were used to select mutations with a potential effect on pre-mRNA splicing. Mutations were experimentally tested using minigene assays. Exons and adjacent intronic sequences were PCR-amplified and cloned in the splicing reporter minigene, and selected mutations were introduced by site-directed mutagenesis. Minigenes were transfected into kidney derived cell lines. RNA from cultured cells was analyzed by RT-PCR and DNA sequencing. Analysis of thirty-three PKD1 exonic mutations revealed three mutations that induce splicing defects. The substitution c.11156G > A, previously predicted as missense mutation p.R3719Q, abolished the donor splice site of intron 38 and resulted in the incorporation of exon 38 with 117 bp of intron 38 and skipping of exon 39. Two synonymous variants, c.327A > T (p.G109G) and c.11257C > A (p.R3753R), generated strong donor splice sites within exons 3 and 39 respectively, resulting in incorporation of incomplete exons. These three nucleotide substitutions represent the first PKD1 exonic mutations that induce aberrant mRNAs. Our results strengthen the importance to evaluate the consequences of presumed missense and synonymous mutations at the mRNA level.