Allele-specific recognition of the 3' splice site of INS intron 1

Understanding Insulin in the Age of Precision Medicine and Big Data: Under-Explored Nature of Genomics

Insulin is amongst the human genome's most well-studied genes/proteins due to its connection to metabolic health. Within this article, we review literature and data to build a knowledge base of Insulin (INS) genetics that influence transcription, transcript processing, translation, hormone maturation, secretion, receptor binding, and metabolism while highlighting the future needs of insulin research. The INS gene region has 2076 unique variants from population genetics. Several variants are found near the transcriptional start site, enhancers, and following the INS transcripts that might influence the readthrough fusion transcript INS-IGF2. This INS-IGF2 transcript splice site was confirmed within hundreds of pancreatic RNAseq samples, lacks drift based on human genome sequencing, and has possible elevated expression due to viral regulation within the liver. Moreover, a rare, poorly characterized African population-enriched variant of INS-IGF2 results in a loss of the stop codon. INS transcript UTR variants rs689 and rs3842753, associated with type 1 diabetes, are found in many pancreatic RNAseq datasets with an elevation of the 3'UTR alternatively spliced INS transcript. Finally, by combining literature, evolutionary profiling, and structural biology, we map rare missense variants that influence preproinsulin translation, proinsulin processing, dimer/hexamer secretory storage, receptor activation, and C-peptide detection for quasi-insulin blood measurements.

Restriction of an intron size en route to endothermy

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

Cancer-Associated Substitutions in RNA Recognition Motifs of PUF60 and U2AF65 Reveal Residues Required for Correct Folding and 3' Splice-Site Selection

U2AF65 (U2AF2) and PUF60 (PUF60) are splicing factors important for recruitment of the U2 small nuclear ribonucleoprotein to lariat branch points and selection of 3′ splice sites (3′ss). Both proteins preferentially bind uridine-rich sequences upstream of 3′ss via their RNA recognition motifs (RRMs). Here, we examined 36 RRM substitutions reported in cancer patients to identify variants that alter 3′ss selection, RNA binding and protein properties. Employing PUF60- and U2AF65-dependent 3′ss previously identified by RNA-seq of depleted cells, we found that 43% (10/23) and 15% (2/13) of independent RRM mutations in U2AF65 and PUF60, respectively, conferred splicing defects. At least three RRM mutations increased skipping of internal U2AF2 (~9%, 2/23) or PUF60 (~8%, 1/13) exons, indicating that cancer-associated RRM mutations can have both cis- and trans-acting effects on splicing. We also report residues required for correct folding/stability of each protein and map functional RRM substitutions on to existing high-resolution structures of U2AF65 and PUF60. These results identify new RRM residues critical for 3′ss selection and provide relatively simple tools to detect clonal RRM mutations that enhance the mRNA isoform diversity.

PUF60-activated exons uncover altered 3' splice-site selection by germline missense mutations in a single RRM

PUF60 is a splicing factor that binds uridine (U)-rich tracts and facilitates association of the U2 small nuclear ribonucleoprotein with primary transcripts. PUF60 deficiency (PD) causes a developmental delay coupled with intellectual disability and spinal, cardiac, ocular and renal defects, but PD pathogenesis is not understood. Using RNA-Seq, we identify human PUF60-regulated exons and show that PUF60 preferentially acts as their activator. PUF60-activated internal exons are enriched for Us upstream of their 3′ splice sites (3′ss), are preceded by longer AG dinucleotide exclusion zones and more distant branch sites, with a higher probability of unpaired interactions across a typical branch site location as compared to control exons. In contrast, PUF60-repressed exons show U-depletion with lower estimates of RNA single-strandedness. We also describe PUF60-regulated, alternatively spliced isoforms encoding other U-bound splicing factors, including PUF60 partners, suggesting that they are co-regulated in the cell, and identify PUF60-regulated exons derived from transposed elements. PD-associated amino-acid substitutions, even within a single RNA recognition motif (RRM), altered selection of competing 3′ss and branch points of a PUF60-dependent exon and the 3′ss choice was also influenced by alternative splicing of PUF60. Finally, we propose that differential distribution of RNA processing steps detected in cells lacking PUF60 and the PUF60-paralog RBM39 is due to the RBM39 RS domain interactions. Together, these results provide new insights into regulation of exon usage by the 3′ss organization and reveal that germline mutation heterogeneity in RRMs can enhance phenotypic variability at the level of splice-site and branch-site selection.

AtU2AF65b functions in abscisic acid mediated flowering via regulating the precursor messenger RNA splicing of ABI5 and FLC in Arabidopsis

In mammalians and yeast, the splicing factor U2AF65/Mud2p functions in precursor messenger RNA (pre-mRNA) processing. Arabidopsis AtU2AF65b encodes a putative U2AF65 but its specific functions in plants are unknown. This paper examines the function of AtU2AF65b as a negative regulator of flowering time in Arabidopsis. We investigated the expression and function of AtU2AF65b in abscisic acid (ABA)-regulated flowering as well as the transcript abundance and pre-mRNA splicing of flowering-related genes in the knock-out mutants of AtU2AF65b. The atu2af65b mutants show early-flowering phenotype under both long-day and short-day conditions. The transcript accumulation of the flowering repressor gene FLOWERING LOCUS C (FLC) is reduced in the shoot apex of atu2af65b, due to both increased intron retention and reduced transcription activation. Reduced transcription of FLC results, at least partially, from the abnormal splicing and reduced transcript abundance of ABSCISIC ACID-INSENSITIVE 5 (ABI5), which encodes an activator of FLC in ABA-regulated flowering signaling. Additionally, the expression of AtU2AF65b is promoted by ABA. Transition to flowering and splicing of FLC and ABI5 in the atu2af65b mutants are compromised during ABA-induced flowering. ABA-responsive AtU2AF65b functions in the pre-mRNA splicing of FLC and ABI5 in shoot apex, whereby AtU2AF65b is involved in ABA-mediated flowering transition in Arabidopsis.

Thioflavin T Monitoring of Guanine Quadruplex Formation in the rs689-Dependent INS Intron 1

The human proinsulin gene (INS) contains a thymine-to-adenine variant (rs689) located in the 3′ splice site (3′ ss) recognition motif of the first intron. The adenine at rs689 is strongly associated with type 1 diabetes. By weakening the polypyrimidine tract, the adenine allele reduces the efficiency of intron 1 splicing, which can be ameliorated by antisense oligonucleotides blocking a splicing silencer located upstream of the 3′ ss. The silencer is surrounded by guanine-rich tracts that may form guanine quadruplexes (G4s) and modulate the accessibility of the silencer. Here, we employed thioflavin T (ThT) to monitor G4 formation in synthetic DNAs and RNAs derived from INS intron 1. We show that the antisense target is surrounded by ThT-positive segments in each direction, with oligoribonucleotides exhibiting consistently higher fluorescence than their DNA counterparts. The signal was reduced for ThT-positive oligonucleotides that were extended into the silencer, indicating that flanking G4s have a potential to mask target accessibility. Real-time monitoring of ThT fluorescence during INS transcription in vitro revealed a negative correlation with ex vivo splicing activities of corresponding INS constructs. Together, these results provide a better characterization of antisense targets in INS primary transcripts for restorative strategies designed to improve the INS splicing defect associated with type 1 diabetes.

U2AF65 enhances milk synthesis and growth of bovine mammary epithelial cells by positively regulating the mTOR-SREBP-1c signalling pathway

U2 snRNP auxiliary factor 65 kDa (U2AF65) is a splicing factor that promotes prespliceosome assembly. The function of U2AF65 in alternative splicing has been identified; however, the essential physiological role of U2AF65 remains poorly understood. In this study, we investigated the regulatory role of U2AF65 in milk synthesis and growth of bovine mammary epithelial cells (BMECs). Our results showed that U2AF65 localizes in the nucleus. Treatment with amino acids (Met and Leu) and hormones (prolactin and β-estradiol) upregulated the expression of U2AF65 in these cells. U2AF65 overexpression increased the synthesis of β-casein, triglycerides, and lactose; increased cell viability; and promoted proliferation of BMECs. Furthermore, our results showed that U2AF65 positively regulated mTOR phosphorylation and expression of mature mRNA of mTOR and SREBP-1c. Collectively, our findings demonstrate that U2AF65 regulates the mRNA expression of signalling molecules (mTOR and SREBP-1c) involved in milk synthesis and growth of BMECs, possibly via controlling the splicing and maturation of these mRNAs. U2 snRNP auxiliary factor 65 kDa (U2AF65) is a splicing factor that promotes prespliceosome assembly. The essential physiological role of U2AF65 remains poorly understood. In the present study, we confirmed that U2AF65 functions as a positive regulator of milk synthesis in and proliferation of bovine mammary epithelial cells via the mTOR-SREBP-1c signalling pathway. Therefore, our study uncovers the regulatory role of U2AF65 in milk synthesis and cell proliferation.

Disease-specific biases in alternative splicing and tissue-specific dysregulation revealed by multitissue profiling of lymphocyte gene expression in type 1 diabetes

Genome-wide association studies (GWAS) have identified multiple, shared allelic associations with many autoimmune diseases. However, the pathogenic contributions of variants residing in risk loci remain unresolved. The location of the majority of shared disease-associated variants in noncoding regions suggests they contribute to risk of autoimmunity through effects on gene expression in the immune system. In the current study, we test this hypothesis by applying RNA sequencing to CD4+, CD8+, and CD19+ lymphocyte populations isolated from 81 subjects with type 1 diabetes (T1D). We characterize and compare the expression patterns across these cell types for three gene sets: all genes, the set of genes implicated in autoimmune disease risk by GWAS, and the subset of these genes specifically implicated in T1D. We performed RNA sequencing and aligned the reads to both the human reference genome and a catalog of all possible splicing events developed from the genome, thereby providing a comprehensive evaluation of the roles of gene expression and alternative splicing (AS) in autoimmunity. Autoimmune candidate genes displayed greater expression specificity in the three lymphocyte populations relative to other genes, with significantly increased levels of splicing events, particularly those predicted to have substantial effects on protein isoform structure and function (e.g., intron retention, exon skipping). The majority of single-nucleotide polymorphisms within T1D-associated loci were also associated with one or more cis-expression quantitative trait loci (cis-eQTLs) and/or splicing eQTLs. Our findings highlight a substantial, and previously underrecognized, role for AS in the pathogenesis of autoimmune disorders and particularly for T1D.

Stoichiometries of U2AF35, U2AF65 and U2 snRNP reveal new early spliceosome assembly pathways

The selection of 3΄ splice sites (3΄ss) is an essential early step in mammalian RNA splicing reactions, but the processes involved are unknown. We have used single molecule methods to test whether the major components implicated in selection, the proteins U2AF35 and U2AF65 and the U2 snRNP, are able to recognize alternative candidate sites or are restricted to one pre-specified site. In the presence of adenosine triphosphate (ATP), all three components bind in a 1:1 stoichiometry with a 3΄ss. Pre-mRNA molecules with two alternative 3΄ss can be bound concurrently by two molecules of U2AF or two U2 snRNPs, so none of the components are restricted. However, concurrent occupancy inhibits splicing. Stoichiometric binding requires conditions consistent with coalescence of the 5΄ and 3΄ sites in a complex (I, initial), but if this cannot form the components show unrestricted and stochastic association. In the absence of ATP, when complex E forms, U2 snRNP association is unrestricted. However, if protein dephosphorylation is prevented, an I-like complex forms with stoichiometric association of U2 snRNPs and the U2 snRNA is base-paired to the pre-mRNA. Complex I differs from complex A in that the formation of complex A is associated with the loss of U2AF65 and 35.

Wild-Type U2AF1 Antagonizes the Splicing Program Characteristic of U2AF1-Mutant Tumors and Is Required for Cell Survival

We have asked how the common S34F mutation in the splicing factor U2AF1 regulates alternative splicing in lung cancer, and why wild-type U2AF1 is retained in cancers with this mutation. A human lung epithelial cell line was genetically modified so that U2AF1S34F is expressed from one of the two endogenous U2AF1 loci. By altering levels of mutant or wild-type U2AF1 in this cell line and by analyzing published data on human lung adenocarcinomas, we show that S34F-associated changes in alternative splicing are proportional to the ratio of S34F:wild-type gene products and not to absolute levels of either the mutant or wild-type factor. Preferential recognition of specific 3' splice sites in S34F-expressing cells is largely explained by differential in vitro RNA-binding affinities of mutant versus wild-type U2AF1 for those same 3' splice sites. Finally, we show that lung adenocarcinoma cell lines bearing U2AF1 mutations do not require the mutant protein for growth in vitro or in vivo. In contrast, wild-type U2AF1 is required for survival, regardless of whether cells carry the U2AF1S34F allele. Our results provide mechanistic explanations of the magnitude of splicing changes observed in U2AF1-mutant cells and why tumors harboring U2AF1 mutations always retain an expressed copy of the wild-type allele.

Upregulation of INS-IGF2 read-through expression and identification of a novel INS-IGF2 splice variant in insulinomas

Fusion transcripts arising from the combination of exons residing on neighboring genes on the same chromosome may give rise to chimeric or novel proteins. Such read-through transcripts have been detected in different cancers where they may be of pathogenetic interest. In this study, we describe for the first time the expression of a read-through transcript in insulinomas, a functioning neuroendocrine pancreatic neoplasm. The read-through transcript INS-IGF2, composed of exons from the two genes proinsulin precursor (INS) and insulin‑like growth factor 2 (IGF2), both mapping to chromosomal subband 11p15.5, was highly expressed in the two insulinomas analyzed. More precisely, version 2 of the INS-IGF2 transcript was expressed, indicating possible expression of the chimeric INS-IGF2 protein. We further identified a novel splice variant of the INS-IGF2 read-through transcript in one of the insulinomas, composed of exon 1 of INS3 and exons of IGF2. In the same tumor, we found high expression of INS3 and the presence of the A allele at SNP rs689. SNP rs689 has been previously described to regulate splicing of the INS transcript, indicating that this regulatory mechanism also affects splicing of INS-IGF2. The identification of the INS-IGF2 read-through transcript specifically in tumor tissue but not in normal pancreatic tissue suggests that high expression of INS-IGF2 could be neoplasia‑specific. These results may have potential clinical applications given that the read-through transcript could be used as a biomarker in insulinoma patients.

Alternative splicing of U2AF1 reveals a shared repression mechanism for duplicated exons

The auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF) facilitates branch point (BP) recognition and formation of lariat introns. The gene for the 35-kD subunit of U2AF gives rise to two protein isoforms (termed U2AF35a and U2AF35b) that are encoded by alternatively spliced exons 3 and Ab, respectively. The splicing recognition sequences of exon 3 are less favorable than exon Ab, yet U2AF35a expression is higher than U2AF35b across tissues. We show that U2AF35b repression is facilitated by weak, closely spaced BPs next to a long polypyrimidine tract of exon Ab. Each BP lacked canonical uridines at position -2 relative to the BP adenines, with efficient U2 base-pairing interactions predicted only for shifted registers reminiscent of programmed ribosomal frameshifting. The BP cluster was compensated by interactions involving unpaired cytosines in an upstream, EvoFold-predicted stem loop (termed ESL) that binds FUBP1/2. Exon Ab inclusion correlated with predicted free energies of mutant ESLs, suggesting that the ESL operates as a conserved rheostat between long inverted repeats upstream of each exon. The isoform-specific U2AF35 expression was U2AF65-dependent, required interactions between the U2AF-homology motif (UHM) and the α6 helix of U2AF35, and was fine-tuned by exon Ab/3 variants. Finally, we identify tandem homologous exons regulated by U2AF and show that their preferential responses to U2AF65-related proteins and SRSF3 are associated with unpaired pre-mRNA segments upstream of U2AF-repressed 3′ss. These results provide new insights into tissue-specific subfunctionalization of duplicated exons in vertebrate evolution and expand the repertoire of exon repression mechanisms that control alternative splicing.

Exon-centric regulation of ATM expression is population-dependent and amenable to antisense modification by pseudoexon targeting

ATM is an important cancer susceptibility gene that encodes a critical apical kinase of the DNA damage response (DDR) pathway. We show that a key nonsense-mediated RNA decay switch exon (NSE) in ATM is repressed by U2AF, PUF60 and hnRNPA1. The NSE activation was haplotype-specific and was most promoted by cytosine at rs609621 in the NSE 3' splice-site (3'ss), which is predominant in high cancer risk populations. NSE levels were deregulated in leukemias and were influenced by the identity of U2AF35 residue 34. We also identify splice-switching oligonucleotides (SSOs) that exploit competition of adjacent pseudoexons to modulate NSE levels. The U2AF-regulated exon usage in the ATM signalling pathway was centred on the MRN/ATM-CHEK2-CDC25-cdc2/cyclin-B axis and preferentially involved transcripts implicated in cancer-associated gene fusions and chromosomal translocations. These results reveal important links between 3'ss control and ATM-dependent responses to double-strand DNA breaks, demonstrate functional plasticity of intronic variants and illustrate versatility of intronic SSOs that target pseudo-3'ss to modify gene expression.

The role of short RNA loops in recognition of a single-hairpin exon derived from a mammalian-wide interspersed repeat

Splice-site selection is controlled by secondary structure through sequestration or approximation of splicing signals in primary transcripts but the exact role of even the simplest and most prevalent structural motifs in exon recognition remains poorly understood. Here we took advantage of a single-hairpin exon that was activated in a mammalian-wide interspersed repeat (MIR) by a mutation stabilizing a terminal triloop, with splice sites positioned close to each other in a lower stem of the hairpin. We first show that the MIR exon inclusion in mRNA correlated inversely with hairpin stabilities. Employing a systematic manipulation of unpaired regions without altering splice-site configuration, we demonstrate a high correlation between exon inclusion of terminal tri- and tetraloop mutants and matching tri-/tetramers in splicing silencers/enhancers. Loop-specific exon inclusion levels and enhancer/silencer associations were preserved across primate cell lines, in 4 hybrid transcripts and also in the context of a distinct stem, but only if its loop-closing base pairs were shared with the MIR hairpin. Unlike terminal loops, splicing activities of internal loop mutants were predicted by their intramolecular Watson-Crick interactions with the antiparallel strand of the MIR hairpin rather than by frequencies of corresponding trinucleotides in splicing silencers/enhancers. We also show that splicing outcome of oligonucleotides targeting the MIR exon depend on the identity of the triloop adjacent to their antisense target. Finally, we identify proteins regulating MIR exon recognition and reveal a distinct requirement of adjacent exons for C-terminal extensions of Tra2α and Tra2β RNA recognition motifs.

Identification of U2AF(35)-dependent exons by RNA-Seq reveals a link between 3' splice-site organization and activity of U2AF-related proteins

The auxiliary factor of U2 small nuclear RNA (U2AF) is a heterodimer consisting of 65- and 35-kD proteins that bind the polypyrimidine tract (PPT) and AG dinucleotides at the 3′ splice site (3′ss). The gene encoding U2AF35 (U2AF1) is alternatively spliced, giving rise to two isoforms U2AF35a and U2AF35b. Here, we knocked down U2AF35 and each isoform and characterized transcriptomes of HEK293 cells with varying U2AF35/U2AF65 and U2AF35a/b ratios. Depletion of both isoforms preferentially modified alternative RNA processing events without widespread failure to recognize 3′ss or constitutive exons. Over a third of differentially used exons were terminal, resulting largely from the use of known alternative polyadenylation (APA) sites. Intronic APA sites activated in depleted cultures were mostly proximal whereas tandem 3′UTR APA was biased toward distal sites. Exons upregulated in depleted cells were preceded by longer AG exclusion zones and PPTs than downregulated or control exons and were largely activated by PUF60 and repressed by CAPERα. The U2AF(35) repression and activation was associated with a significant interchange in the average probabilities to form single-stranded RNA in the optimal PPT and branch site locations and sequences further upstream. Although most differentially used exons were responsive to both U2AF subunits and their inclusion correlated with U2AF levels, a small number of transcripts exhibited distinct responses to U2AF35a and U2AF35b, supporting the existence of isoform-specific interactions. These results provide new insights into function of U2AF and U2AF35 in alternative RNA processing.

Optimal antisense target reducing INS intron 1 retention is adjacent to a parallel G quadruplex

Splice-switching oligonucleotides (SSOs) have been widely used to inhibit exon usage but antisense strategies that promote removal of entire introns to increase splicing-mediated gene expression have not been developed. Here we show reduction of INS intron 1 retention by SSOs that bind transcripts derived from a human haplotype expressing low levels of proinsulin. This haplotype is tagged by a polypyrimidine tract variant rs689 that decreases the efficiency of intron 1 splicing and increases the relative abundance of mRNAs with extended 5' untranslated region (5' UTR), which curtails translation. Co-expression of haplotype-specific reporter constructs with SSOs bound to splicing regulatory motifs and decoy splice sites in primary transcripts revealed a motif that significantly reduced intron 1-containing mRNAs. Using an antisense microwalk at a single nucleotide resolution, the optimal target was mapped to a splicing silencer containing two pseudoacceptor sites sandwiched between predicted RNA guanine (G) quadruplex structures. Circular dichroism spectroscopy and nuclear magnetic resonance of synthetic G-rich oligoribonucleotide tracts derived from this region showed formation of a stable parallel 2-quartet G-quadruplex on the 3' side of the antisense retention target and an equilibrium between quadruplexes and stable hairpin-loop structures bound by optimal SSOs. This region interacts with heterogeneous nuclear ribonucleoproteins F and H that may interfere with conformational transitions involving the antisense target. The SSO-assisted promotion of weak intron removal from the 5' UTR through competing noncanonical and canonical RNA structures may facilitate development of novel strategies to enhance gene expression.

BRCA1 exon 11 a model of long exon splicing regulation

BRCA1 exon 11 is one of the biggest human exons, spanning 3426 bases. This gene is potentially involved in DNA repair as well as cell growth and cell cycle control. Exon 11 is regulated at the splicing level producing three main different combinations of BRCA1 mature transcripts; one including the whole of exon 11 (full isoform), one skipping the entire exon (D11 isoform), and one including only 117 base pairs of exon 11 (D11q isoform). Using minigene and deletion analyses, we have previously described important splicing regulatory sequences located at the beginning of this exon (5' end). We have now found additional important sequences located at its 3' end. In particular, we describe the presence of a strong splicing enhancer adjacent to the downstream 5' splice site, which minimizes competition from an upstream 5' splice site and so ensures long exon inclusion. Analyses of the proteins binding these RNA sequences have revealed that Tra2beta and hnRNP L are involved in the regulation of BRCA1 exon 11 by influencing the recognition of donor sites. Interestingly, BRCA1 exon 11 carrying deletion of the regulatory sequences bound by these factors also showed unexpected responses to up- or downregulation of these regulatory proteins, suggesting that they can also bind elsewhere in this large exon and elicit different effects on its recognition.   The identification of sequences and proteins relevant for the regulation of BRCA1 exon 11 now provides better knowledge on how this exon is recognized and may represent an important step toward understanding how large exons are regulated.

Recurrent mutations in the U2AF1 splicing factor in myelodysplastic syndromes

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders that often progress to chemotherapy-resistant secondary acute myeloid leukemia (sAML). We used whole-genome sequencing to perform an unbiased comprehensive screen to discover the somatic mutations in a sample from an individual with sAML and genotyped the loci containing these mutations in the matched MDS sample. Here we show that a missense mutation affecting the serine at codon 34 (Ser34) in U2AF1 was recurrently present in 13 out of 150 (8.7%) subjects with de novo MDS, and we found suggestive evidence of an increased risk of progression to sAML associated with this mutation. U2AF1 is a U2 auxiliary factor protein that recognizes the AG splice acceptor dinucleotide at the 3' end of introns, and the alterations in U2AF1 are located in highly conserved zinc fingers of this protein. Mutant U2AF1 promotes enhanced splicing and exon skipping in reporter assays in vitro. This previously unidentified, recurrent mutation in U2AF1 implicates altered pre-mRNA splicing as a potential mechanism for MDS pathogenesis.

Compensatory signals associated with the activation of human GC 5' splice sites

GC 5′ splice sites (5′ss) are present in ∼1% of human introns, but factors promoting their efficient selection are poorly understood. Here, we describe a case of X-linked agammaglobulinemia resulting from a GC 5′ss activated by a mutation in BTK intron 3. This GC 5′ss was intrinsically weak, yet it was selected in >90% primary transcripts in the presence of a strong and intact natural GT counterpart. We show that efficient selection of this GC 5′ss required a high density of GAA/CAA-containing splicing enhancers in the exonized segment and was promoted by SR proteins 9G8, Tra2β and SC35. The GC 5′ss was efficiently inhibited by splice-switching oligonucleotides targeting either the GC 5′ss itself or the enhancer. Comprehensive analysis of natural GC-AG introns and previously reported pathogenic GC 5′ss showed that their efficient activation was facilitated by higher densities of splicing enhancers and lower densities of silencers than their GT 5′ss equivalents. Removal of the GC-AG introns was promoted to a minor extent by the splice-site strength of adjacent exons and inhibited by flanking Alu repeats, with the first downstream Alus located on average at a longer distance from the GC 5′ss than other transposable elements. These results provide new insights into the splicing code that governs selection of noncanonical splice sites.

Prediction of single-nucleotide substitutions that result in exon skipping: identification of a splicing silencer in BRCA1 exon 6

Missense, nonsense, and translationally silent mutations can inactivate genes by altering the inclusion of mutant exons in mRNA, but their overall frequency among disease-causing exonic substitutions is unknown. Here, we have tested missense and silent mutations deposited in the BRCA1 mutation databases of unclassified variants for their effects on exon inclusion. Analysis of 21 BRCA1 variants using minigene assays revealed a single exon-skipping mutation c.231G>T. Comprehensive mutagenesis of an adjacent 12-nt segment showed that this silent mutation resulted in a higher level of exon skipping than the 35 other single-nucleotide substitutions. Exon inclusion levels of mutant constructs correlated significantly with predicted splicing enhancers/silencers, prompting the development of two online utilities freely available at http://www.dbass.org.uk. EX-SKIP quickly estimates which allele is more susceptible to exon skipping, whereas HOT-SKIP examines all possible mutations at each exon position and identifies candidate exon-skipping positions/substitutions. We demonstrate that the distribution of exon-skipping and disease-associated substitutions previously identified in coding regions was biased toward top-ranking HOT-SKIP mutations. Finally, we show that proteins 9G8, SC35, SF2/ASF, Tra2, and hnRNP A1 were associated with significant alterations of BRCA1 exon 6 inclusion in the mRNA. Together, these results facilitate prediction of exonic substitutions that reduce exon inclusion in mature transcripts.

AG-dependent 3'-splice sites are predisposed to aberrant splicing due to a mutation at the first nucleotide of an exon

In pre-mRNA splicing, a conserved AG/G at the 3′-splice site is recognized by U2AF³⁵. A disease-causing mutation abrogating the G nucleotide at the first position of an exon (E⁺¹) causes exon skipping in GH1, FECH and EYA1, but not in LPL or HEXA. Knockdown of U2AF³⁵ enhanced exon skipping in GH1 and FECH. RNA-EMSA revealed that wild-type FECH requires U2AF³⁵ but wild-type LPL does not. A series of artificial mutations in the polypyrimidine tracts of GH1, FECH, EYA1, LPL and HEXA disclosed that a stretch of at least 10–15 pyrimidines is required to ensure normal splicing in the presence of a mutation at E⁺¹. Analysis of nine other disease-causing mutations at E⁺¹ detected five splicing mutations. Our studies suggest that a mutation at the AG-dependent 3′-splice site that requires U2AF³⁵ for spliceosome assembly causes exon skipping, whereas one at the AG-independent 3′-splice site that does not require U2AF³⁵ gives rise to normal splicing. The AG-dependence of the 3′-splice site that we analyzed in disease-causing mutations at E⁺¹ potentially helps identify yet unrecognized splicing mutations at E⁺¹.