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

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 amongst 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 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

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 3

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

Missense, nonsense and translationally silent mutations can inactivate genes by altering the inclusion of mutant exons in mRNA, but their overall frequency amongst 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 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

Position-dependent repression and promotion of DQB1 intron 3 splicing by GGGG motifs

Alternative splicing of HLA-DQB1 exon 4 is allele-dependent and results in variable expression of soluble DQbeta. We have recently shown that differential inclusion of this exon in mature transcripts is largely due to intron 3 variants in the branch point sequence (BPS) and polypyrimidine tract. To identify additional regulatory cis-elements that contribute to haplotype-specific splicing of DQB1, we systematically examined the effect of guanosine (G) repeats on intron 3 removal. We found that the GGG or GGGG repeats generally improved splicing of DQB1 intron 3, except for those that were adjacent to the 5' splice site where they had the opposite effect. The most prominent splicing enhancement was conferred by GGGG motifs arranged in tandem upstream of the BPS. Replacement of a G-rich segment just 5' of the BPS with a series of random sequences markedly repressed splicing, whereas substitutions of a segment further upstream that lacked the G-rich elements and had the same size did not result in comparable splicing inhibition. Systematic mutagenesis of both suprabranch guanosine quadruplets (G4) revealed a key role of central G residues in splicing enhancement, whereas cytosines in these positions had the most prominent repressive effects. Together, these results show a significant role of tandem G4NG4 structures in splicing of both complete and truncated DQB1 intron 3, support position dependency of G repeats in splicing promotion and inhibition, and identify positively and negatively acting sequences that contribute to the haplotype-specific DQB1 expression

SERPING1 rs2511988 and age-related macular degeneration

Using a two-stage case-control protocol, Sarah Ennis and colleagues (Nov 22, p 1828)1 show an association between the SERPING1 gene and age-related macular degeneration. Although this exciting study provides substantive evidence for association with several closely linked intragenic and extragenic variants, it is unclear which biological mechanisms are affected by the population variability in SERPING1 and which variants have a major role in genetic predisposition to this disease.We inspected gene sequences flanking the newly identified variants and found that one of them (rs2511988) was highly likely to have functional consequences. This single-nucleotide substitution is located 20 base pairs upstream of the 3? splice site of exon 7 and is much closer to the SERPING1 coding sequence than the rs2511989 variant initially discovered after completion of their first screen. The rs2511988 variant (A/G) is adjacent to a branch point sequence (TGTTAAG; branch point is underlined), which has been predicted computationally.2Mutations or variants located in the vicinity of this key exon recognition signal have been shown to promote or inhibit splicing efficiency in several human precursor messenger RNAs.3 In yeast, adenine at the same position relative to the branch site is the preferred base4 and facilitates interactions with the branch point binding protein.5 Although human branch point sequences are much more degenerate than those in yeast,[2] and [3] the rs2511988 variant is likely to alter efficiency of intron 6 removal by interfering with assembly of protein-RNA complexes at the 3? splice site and, ultimately, SERPING1 expression level

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

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

Exonic splicing code and protein binding sites for calcium

Auxilliary splicing sequences in exons, known as enhancers (ESEs) and silencers (ESSs), have been subject to strong selection pressures at the RNA and protein level. The protein component of this splicing code is substantial, recently estimated at ∼50% of the total information within ESEs, but remains poorly understood. The ESE/ESS profiles were previously associated with the Irving-Williams (I-W) stability series for divalent metals, suggesting that the ESE/ESS evolution was shaped by metal binding sites. Here, we have examined splicing activities of exonic sequences that encode protein binding sites for Ca2+, a weak binder in the I-W affinity order. We found that predicted exon inclusion levels for the EFhand motifs and for Ca2+-binding residues in nonEFhand proteins were higher than for average exons. For canonical EF-hands, the increase was centred on the EF-hand chelation loop and, in particular, on Ca2+-coordinating residues, with a 1>12>3∼5>9 hierarchy in the 12-codon loop consensus and usage bias at codons 1 and 12. The same hierarchy but a lower increase was observed for noncanonical EFhands, except for S100 proteins. EF-hand loops preferentially accumulated exon splits in two clusters, one located in their N-terminal halves and the other around codon 12. Using splicing assays and published crosslinking and immunoprecipitation data, we identify candidate trans-acting factors that preferentially bind conserved GA-rich motifs encoding negatively charged amino acids in the loops. Together, these data provide evidence for the high capacity of codons for Ca2+-coordinating residues to be retained in mature transcripts, facilitating their exonlevel expansion during eukaryotic evolution