Exon sequences at the splice junctions affect splicing fidelity and alternative splicing

Mechanisms and regulation of spliceosome-mediated pre-mRNA splicing in Saccharomyces cerevisiae

Pre-mRNA splicing, the removal of introns and ligation of flanking exons, is a crucial step in eukaryotic gene expression. The spliceosome, a macromolecular complex made up of five small nuclear RNAs (snRNAs) and dozens of proteins, assembles on introns via a complex pathway before catalyzing the two transesterification reactions necessary for splicing. All of these steps have the potential to be highly regulated to ensure correct mRNA isoform production for proper cellular function. While Saccharomyces cerevisiae (yeast) has a limited set of intron-containing genes, many of these genes are highly expressed, resulting in a large number of transcripts in a cell being spliced. As a result, splicing regulation is of critical importance for yeast. Just as in humans, yeast splicing can be influenced by protein components of the splicing machinery, structures and properties of the pre-mRNA itself, or by the action of trans-acting factors. It is likely that further analysis of the mechanisms and pathways of splicing regulation in yeast can reveal general principles applicable to other eukaryotes. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing.

Splicing factor Prp18p promotes genome-wide fidelity of consensus 3'-splice sites

The fidelity of splice site selection is critical for proper gene expression. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is challenging considering the low complexity of the 3'SS consensus sequence YAG. Here, we show that absence of the Prp18p splicing factor results in genome-wide activation of alternative 3'SS in S. cerevisiae, including highly unusual non-YAG sequences. Usage of these non-canonical 3'SS in the absence of Prp18p is enhanced by upstream poly(U) tracts and by their potential to interact with the first intronic nucleoside, allowing them to dock in the spliceosome active site instead of the normal 3'SS. The role of Prp18p in 3'SS fidelity is facilitated by interactions with Slu7p and Prp8p, but cannot be fulfilled by Slu7p, identifying a unique role for Prp18p in 3'SS fidelity. This fidelity function is synergized by the downstream proofreading activity of the Prp22p helicase, but is independent from another late splicing helicase, Prp43p. Our results show that spliceosomes exhibit remarkably relaxed 3'SS sequence usage in the absence of Prp18p and identify a network of spliceosomal interactions centered on Prp18p which are required to promote the fidelity of the recognition of consensus 3'SS sequences.

Predominant Liver Cystic Disease in a New Heterozygotic PKHD1 Variant: A Case Report

BACKGROUND The polycystic kidney and hepatic disease 1 (PKHD1) gene codes for fibrocystin-polyductin, a protein that takes part in cell-signaling for cell differentiation, especially in kidney tubules and bile ducts. A homozygous or compound heterozygous defect in this gene can cause autosomal recessive polycystic kidney disease (ARPKD). Polycystic liver disease (PCLD) can also be caused by single heterozygous variants in the PKHD1 gene. ARPKD presents with renal insufficiency and cystic dilatation of bile ducts, although disease is not expected with a single heterozygous mutation. PCLD presents with multiple cysts in the liver and dilated bile ducts as well, but with less of an impact on the kidneys than with ARPKD. Our purpose in publishing this report is to introduce an as-yet unknown variant to the body of genetic defects associated with ARPKD and PCLD, as well as to argue for the likely pathogenicity of the variant according to the prevailing criteria used for classifying gene variants. CASE REPORT We present a patient with a de novo PKHD1 variant currently classified as a variant of unknown significance manifesting with bilaterally enlarged cystic kidneys and echogenic cystic structures in the hepatic portal system, indicative of cystic disease. CONCLUSIONS Given this patient's liver and kidney presentation that does not fully align with either ARPKD or PCLD, the authors believe that the single heterozygous variant in this patient's PKHD1 gene is worthy of reporting. This new single heterozygous variant in PKHD1 gene causing cystic kidney and cystic hepatic disease in the patient should be considered 'likely pathogenic' according to the criteria set by the American College of Medical Genetics.

Spliceosomal mutations decouple 3' splice site fidelity from cellular fitness

The fidelity of splice site selection is thought to be critical for proper gene expression and cellular fitness. In particular, proper recognition of 3'-splice site (3'SS) sequences by the spliceosome is a daunting task considering the low complexity of the 3'SS consensus sequence YAG. Here we show that inactivating the near-essential splicing factor Prp18p results in a global activation of alternative 3'SS, many of which harbor sequences that highly diverge from the YAG consensus, including some highly unusual non-AG 3'SS. We show that the role of Prp18p in 3'SS fidelity is promoted by physical interactions with the essential splicing factors Slu7p and Prp8p and synergized by the proofreading activity of the Prp22p helicase. Strikingly, structure-guided point mutations that disrupt Prp18p-Slu7p and Prp18p-Prp8p interactions mimic the loss of 3'SS fidelity without any impact on cellular growth, suggesting that accumulation of incorrectly spliced transcripts does not have a major deleterious effect on cellular viability. These results show that spliceosomes exhibit remarkably relaxed fidelity in the absence of Prp18p, and that new 3'SS sampling can be achieved genome-wide without a major negative impact on cellular fitness, a feature that could be used during evolution to explore new productive alternative splice sites.

Adenosine Deaminase 2 Deficiency Caused by Biallele Variants Including Splicing Variant: The First Case in Korea

Deficiency of adenosine deaminase 2 (DADA2) is an autoinflammatory disease caused by pathogenic variants of the ADA2 gene and has similar clinical features to polyarteritis nodosa (PAN). We, herein, report a case of DADA2 in Korea that was diagnosed in a patient with childhood-onset PAN. The patient had a truncal ataxia and facial palsy caused by thalamic infarction at 34 months of age. Livedo reticularis with Raynaud phenomenon and abdominal pain with fever were followed. Radiologic examination showed multiple infarctions in brain and kidney. She was diagnosed with PAN using skin biopsy and angiography. She had severe hemorrhagic strokes despite medical treatments. Her disease activity was controlled after adding a tumor necrosis factor-α inhibitor. Molecular analysis revealed compound heterozygous pathogenic variants of ADA2 gene. This is the first case of DADA2 in Korea. Genetic analysis for ADA2 gene should be considered in patients with childhood-onset PAN.

Dynamic visualization of mRNA splicing variants with a transactivating reporter

Dynamic changes in intron sequences, with their loss and gain, are poorly detected due to the limited methods for the non-invasive monitoring of the pre-mRNA splicing process. Here, we describe the design of a two-step transcriptional activation (TSTA) reporter for the real-time imaging of the splicing process in living subjects. By taking advantage of the strong transactivating properties of the GAL4-VP16 fusion protein, which can target upstream activation sequence (UAS) elements to boost subsequent firefly luciferase reporter gene expression, we successfully and consistently detected the dynamic pre-mRNA splicing activity in response to exogenous splicing modulators in living cells and animals. Our findings provide a valuable tool for the high-throughput screening of splicing modulators, which could speed up the development of new drugs for the treatment of disordered splicing diseases.

Molecular Analysis of CYP27B1 Mutations in Vitamin D-Dependent Rickets Type 1A: c.590G > A (p.G197D) Missense Mutation Causes a RNA Splicing Error

Context

Vitamin D-dependent rickets type 1A (VDDR1A) is a rare autosomal recessively inherited disorder due to loss-of-function mutations in the CYP27B1 gene. CYP27B1 encodes an enzyme of 25-hydroxyvitamin D-1α-hydroxylase for converting inactive 25-OHD to biologically active 1,25-(OH)₂D.

Objective

To identify underlying genetic defects in patients with VDDR1A.

Methods

Twelve patients from 7 Turkish and 2 Saudi families were investigated. The coding exons and intron-exon boundaries of the CYP27B1 gene were amplified by Polymerase Chain Reaction (PCR) from peripheral lymphocyte DNA. PCR products were directly sequenced. The consequences of c.590G > A mutation were analyzed by in silico and functional analysis.

Results

CYP27B1 mutations were identified in all the patients. Two novel mutations were identified in two separate families: c.171delG (family 7) and c.398_400dupAAT (family 8). The intra-exon deletion of c.171delG resulted in a frameshift and premature stop codon 20 amino acids downstream from the mutation (p.L58Cfs^∗20). The intra-exon duplication of c.398_400dupAAT generated a premature stop codon at the mutation site (p.W134^∗). A missense c.590G > A (p.G197D) mutation was found in a patient from family 4 and caused a defect in pre-mRNA splicing. As a result, two populations of transcripts were detected: the majority of them with intron 3 retention (83%), and the minority (17%) being properly spliced transcripts with about 16% of wild-type enzymatic activity. The remaining nine patients from six families carried a previously reported c.1319_1325dupCCCACCC (F443Pfs^∗24) mutation. Clinically, all the patients need continued calcitriol treatment, which was consistent with inactivation of 25-hydroxy vitamin D1α-hydroxylase activity.

Conclusion

Two novel frameshift CYP27B1 mutations were identified and predicted to inactivate 25-hydroxyvitamin D-1α-hydroxylase. The loss of enzymatic activity by c.590G > A missense mutation was mainly caused by aberrant pre-mRNA splicing.

RNA-Primed Amplification for Noise-Suppressed Visualization of Single-Cell Splice Variants

Splice variants visualization is pivotal for a deeper understanding of cell growth and development. However, it remains technically challenging due to short lengths, similar sequences, and low abundance. The existing single-cell imaging strategies suffer from nonspecific amplification that causes considerable noise during visualization of the splice variants. Herein we develop a new RNA-primed amplification strategy for noise-suppressed visualization of single-cell splice variants. Block probes were designed to specifically identify the conjugated region of exons in mRNA, which was then digested by endonuclease and provided a hydroxyl group at the 3' terminal. The RNA target can act as primer to trigger rolling circle amplification, achieving visualization of splice variants with noise suppressed to nearly zero. We further explored the expression and distribution of BRCA1 splice variants in three breast cell lines, revealing cell-type specific mapping of this cancer suppressor gene.

RNA splicing process analysis for identifying antisense oligonucleotide inhibitors with padlock probe-based isothermal amplification

We report a highly sensitive method for quantifying the splicing products in different steps, enabling us to analyze the splicing process and identify ASO inhibitors.

RNA splicing, which mainly involves two transesterification steps, is a fundamental process of gene expression and its abnormal regulation contributes to serious genetic diseases. Antisense oligonucleotides (ASOs) are genetic control tools that can be used to specifically control genes through alteration of the RNA splicing pathway. Despite intensive research, how ASOs or various other factors influence the multiple processes of RNA splicing still remains obscure. This is largely due to an inability to analyze the splicing efficiency of each step in the RNA splicing process with high sensitivity. We addressed this limitation by introducing a padlock probe-based isothermal amplification assay to achieve quantification of the specific products in different splicing steps. With this amplified assay, the roles that ASOs play in RNA splicing inhibition in the first and second steps could be distinguished. We identified that 5′-ASO could block RNA splicing by inhibiting the first step, while 3′-ASO could block RNA splicing by inhibiting the second step. This method provides a versatile tool for assisting efficient ASO design and discovering new splicing modulators and therapeutic drugs.

Lessons from the canine Oxtr gene: populations, variants and functional aspects

Oxytocin receptor (OXTR) acts as a key behavioral modulator of the central nervous system, affecting social behavior, stress, affiliation and cognitive functions. Variants of the Oxtr gene are known to influence behavior both in animals and humans; however, canine Oxtr polymorphisms are less characterized in terms of possible relevance to function, selection criteria in breeding and domestication. In this report, we provide a detailed characterization of common variants of the canine Oxtr gene. In particular (1) novel polymorphisms were identified by direct sequencing of wolf and dog samples, (2) allelic distributions and pairwise linkage disequilibrium patterns of several canine populations were compared, (3) neighbor joining (NJ) tree based on common single nucleotide polymorphisms (SNPs) was constructed, (4) mRNA expression features were assessed, (5) a novel splice variant was detected and (6) in vitro functional assays were performed. Results indicate marked differences regarding Oxtr variations between purebred dogs of different breeds, free-ranging dog populations, wolf subspecies and golden jackals. This, together with existence of explicitly dog-specific alleles and data obtained from the NJ tree implies that Oxtr could indeed have been a target gene during domestication and selection for human preferred aspects of temperament and social behavior. This assumption is further supported by the present observations on gene expression patterns within the brain and luciferase reporter experiments, providing a molecular level link between certain canine Oxtr polymorphisms and differences in nervous system function and behavior.

The Fission Yeast Pre-mRNA-processing Factor 18 (prp18+) Has Intron-specific Splicing Functions with Links to G1-S Cell Cycle Progression

The fission yeast genome, which contains numerous short introns, is an apt model for studies on fungal splicing mechanisms and splicing by intron definition. Here we perform a domain analysis of the evolutionarily conserved Schizosaccharomyces pombe pre-mRNA-processing factor, SpPrp18. Our mutational and biophysical analyses of the C-terminal α-helical bundle reveal critical roles for the conserved region as well as helix five. We generate a novel conditional missense mutant, spprp18-5 To assess the role of SpPrp18, we performed global splicing analyses on cells depleted of prp18⁺ and the conditional spprp18-5 mutant, which show widespread but intron-specific defects. In the absence of functional SpPrp18, primer extension analyses on a tfIId⁺ intron 1-containing minitranscript show accumulated pre-mRNA, whereas the lariat intron-exon 2 splicing intermediate was undetectable. These phenotypes also occurred in cells lacking both SpPrp18 and SpDbr1 (lariat debranching enzyme), a genetic background suitable for detection of lariat RNAs. These data indicate a major precatalytic splicing arrest that is corroborated by the genetic interaction between spprp18-5 and spprp2-1, a mutant in the early acting U2AF59 protein. Interestingly, SpPrp18 depletion caused cell cycle arrest before S phase. The compromised splicing of transcripts coding for G₁-S regulators, such as Res2, a transcription factor, and Skp1, a regulated proteolysis factor, are shown. The cumulative effects of SpPrp18-dependent intron splicing partly explain the G₁ arrest upon the loss of SpPrp18. Our study using conditional depletion of spprp18⁺ and the spprp18-5 mutant uncovers an intron-specific splicing function and early spliceosomal interactions and suggests links with cell cycle progression.

Alternative Splicing in Next Generation Sequencing Data of Saccharomyces cerevisiae

mRNA splicing is required in about 4% of protein coding genes in Saccharomyces cerevisiae. The gene structure of those genes is simple, generally comprising two exons and one intron. In order to characterize the impact of alternative splicing on the S. cerevisiae transcriptome, we perform a systematic analysis of mRNA sequencing data. We find evidence of a pervasive use of alternative splice sites and detect several novel introns both within and outside protein coding regions. We also find a predominance of alternative splicing on the 3’ side of introns, a finding which is consistent with existing knowledge on conservation of exon-intron boundaries in S. cerevisiae. Some of the alternatively spliced transcripts allow for a translation into different protein products.

RNA-Binding Proteins: Splicing Factors and Disease

Pre-mRNA splicing is mediated by interactions of the Core Spliceosome and an array of accessory RNA binding proteins with cis-sequence elements. Splicing is a major regulatory component in higher eukaryotes. Disruptions in splicing are a major contributor to human disease. One in three hereditary disease alleles are believed to cause aberrant splicing. Hereditary disease alleles can alter splicing by disrupting a splicing element, creating a toxic RNA, or affecting splicing factors. One of the challenges of medical genetics is identifying causal variants from the thousands of possibilities discovered in a clinical sequencing experiment. Here we review the basic biochemistry of splicing, the mechanisms of splicing mutations, the methods for identifying splicing mutants, and the potential of therapeutic interventions.

Cwc21p promotes the second step conformation of the spliceosome and modulates 3' splice site selection

Pre-mRNA splicing involves two transesterification steps catalyzed by the spliceosome. How RNA substrates are positioned in each step and the molecular rearrangements involved, remain obscure. Here, we show that mutations in PRP16, PRP8, SNU114 and the U5 snRNA that affect this process interact genetically with CWC21, that encodes the yeast orthologue of the human SR protein, SRm300/SRRM2. Our microarray analysis shows changes in 3′ splice site selection at elevated temperature in a subset of introns in cwc21Δ cells. Considering all the available data, we propose a role for Cwc21p positioning the 3′ splice site at the transition to the second step conformation of the spliceosome, mediated through its interactions with the U5 snRNP. This suggests a mechanism whereby SRm300/SRRM2, might influence splice site selection in human cells.

TP53 mutations in human cancer: database reassessment and prospects for the next decade

More than 50% of human tumors carry TP53 gene mutations and in consequence more than 45,000 somatic and germline mutations have been gathered in the UMD TP53 database (http://p53.fr). Analyses of these mutations have been invaluable for bettering our knowledge on the structure-function relationships within the TP53 protein and the high degree of heterogeneity of the various TP53 mutants in human cancer. In this review, we discuss how with the release of the sequences of thousands of tumor genomes issued from high-throughput sequencing, the description of novel TP53 mutants is now reaching a plateau indicating that we are close to the full set of mutants that target the elusive tumor-suppressive activity of this protein. We performed an extensive and thorough analysis of the TP53 mutation database, focusing particularly on specific sets of mutations that were overlooked in the past because of their low frequencies, for example, synonymous mutations, splice mutations, or mutations-targeting residues subject to posttranslational modifications. We also discuss the evolution of the statistical methods used to differentiate TP53 passenger mutations and artifactual data from true mutations, a process vital to the release of an accurate TP53 mutation database that will in turn be an invaluable tool for both clinicians and researchers.

Widespread use of non-productive alternative splice sites in Saccharomyces cerevisiae

Saccharomyces cerevisiae has been used as a model system to investigate the mechanisms of pre-mRNA splicing but only a few examples of alternative splice site usage have been described in this organism. Using RNA-Seq analysis of nonsense-mediated mRNA decay (NMD) mutant strains, we show that many S. cerevisiae intron-containing genes exhibit usage of alternative splice sites, but many transcripts generated by splicing at these sites are non-functional because they introduce premature termination codons, leading to degradation by NMD. Analysis of splicing mutants combined with NMD inactivation revealed the role of specific splicing factors in governing the use of these alternative splice sites and identified novel functions for Prp17p in enhancing the use of branchpoint-proximal upstream 3′ splice sites and for Prp18p in suppressing the usage of a non-canonical AUG 3′-splice site in GCR1. The use of non-productive alternative splice sites can be increased in stress conditions in a promoter-dependent manner, contributing to the down-regulation of genes during stress. These results show that alternative splicing is frequent in S. cerevisiae but masked by RNA degradation and that the use of alternative splice sites in this organism is mostly aimed at controlling transcript levels rather than increasing proteome diversity.

Accurate gene expression requires the transfer of gene information from DNA to RNA. When DNA is transcribed into RNA, part of the RNA needs to be removed (spliced) to generate a proper copy of the genetic information. This process needs to be very accurate to preserve the genetic information that will be transferred into proteins. Our study shows that in baker's yeast, the splicing process does not always produce the correctly spliced products, as RNA splicing events frequently utilize incorrect splice sites. However, these deficient RNA molecules are eliminated from cells by a quality control mechanism to preserve the integrity of the genetic information. However, incorrect splicing is not useless, as it can be used to regulate the quantity of RNA that is generated.

Intronic features that determine the selection of the 3' splice site

Most eukaryotic primary transcripts include segments, or introns, that will be accurately removed during RNA biogenesis. This process, known as pre-messenger RNA splicing, is catalyzed by the spliceosome, accurately selecting a set of intronic marks from others apparently equivalent. This identification is critical, as incorrectly spliced RNAs can be toxic for the organism. One of these marks, the dinucleotide AG, signals the intronic 3' end, or 3' splice site (ss). In this review we will focus on those intronic features that have an impact on 3' ss selection. These include the location and type of neighboring sequences, and their distance to the 3' end. We will see that their interplay is needed to select the right intronic end, and that this can be modulated by additional intronic elements that contribute to alternative splicing, whereby diverse RNAs can be generated from identical precursors. This complexity, still poorly understood, is fundamental for the accuracy of gene expression. In addition, a clear knowledge of 3' ss selection is needed to fully decipher the coding potential of genomes.

CEF1/CDC5 alleles modulate transitions between catalytic conformations of the spliceosome

Conformational change within the spliceosome is required between the first and second catalytic steps of pre-mRNA splicing. A prior genetic screen for suppressors of an intron mutant that stalls between the two steps yielded both prp8 and non-prp8 alleles that suppressed second-step splicing defects. We have now identified the strongest non-prp8 suppressors as alleles of the NTC (Prp19 complex) component, CEF1. These cef1 alleles generally suppress second-step defects caused by a variety of intron mutations, mutations in U6 snRNA, or deletion of the second-step protein factor Prp17, and they can activate alternative 3' splice sites. Genetic and functional interactions between cef1 and prp8 alleles suggest that they modulate the same event(s) in the first-to-second-step transition, most likely by stabilization of the second-step spliceosome; in contrast, alleles of U6 snRNA that also alter this transition modulate a distinct event, most likely by stabilization of the first-step spliceosome. These results implicate a myb-like domain of Cef1/CDC5 in interactions that modulate conformational states of the spliceosome and suggest that alteration of these events affects splice site use, resulting in alternative splicing-like patterns in yeast.

Factors affecting splicing strength of yeast genes

Accurate and efficient splicing is of crucial importance for highly-transcribed intron-containing genes (ICGs) in rapidly replicating unicellular eukaryotes such as the budding yeast Saccharomyces cerevisiae. We characterize the 5' and 3' splice sites (ss) by position weight matrix scores (PWMSs), which is the highest for the consensus sequence and the lowest for splice sites differing most from the consensus sequence and used PWMS as a proxy for splicing strength. HAC1, which is known to be spliced by a nonspliceosomal mechanism, has the most negative PWMS for both its 5' ss and 3' ss. Several genes under strong splicing regulation and requiring additional splicing factors for their splicing also have small or negative PWMS values. Splicing strength is higher for highly transcribed ICGs than for lowly transcribed ICGs and higher for transcripts that bind strongly to spliceosomes than those that bind weakly. The 3' splice site features a prominent poly-U tract before the 3'AG. Our results suggest the potential of using PWMS as a screening tool for ICGs that are either spliced by a nonspliceosome mechanism or under strong splicing regulation in yeast and other fungal species.

Deciphering 3'ss selection in the yeast genome reveals an RNA thermosensor that mediates alternative splicing

Poor understanding of the spliceosomal mechanisms to select intronic 3' ends (3'ss) is a major obstacle to deciphering eukaryotic genomes. Here, we discern the rules for global 3'ss selection in yeast. We show that, in contrast to the uniformity of yeast splicing, the spliceosome uses all available 3'ss within a distance window from the intronic branch site (BS), and that in ∼70% of all possible 3'ss this is likely to be mediated by pre-mRNA structures. Our results reveal that one of these RNA folds acts as an RNA thermosensor, modulating alternative splicing in response to heat shock by controlling alternate 3'ss availability. Thus, our data point to a deeper role for the pre-mRNA in the control of its own fate, and to a simple mechanism for some alternative splicing.

The mechanism of the second step of pre-mRNA splicing

The molecular mechanisms of the second step of pre-mRNA splicing in yeast and higher eukaryotes are reviewed. The important elements in the pre-mRNA, the participating proteins, and the proposed secondary structures and roles of the snRNAs are described. The sequence of events in the second step is presented, focusing on the actions of the proteins in setting up and facilitating the second reaction. Mechanisms for avoiding errors in splicing are discussed.

Determinants of Nam8-dependent splicing of meiotic pre-mRNAs

Nam8, a component of yeast U1 snRNP, is optional for mitotic growth but required during meiosis, because Nam8 collaborates with Mer1 to promote splicing of essential meiotic mRNAs AMA1, MER2 and MER3. Here, we identify SPO22 and PCH2 as novel targets of Nam8-dependent meiotic splicing. Whereas SPO22 splicing is co-dependent on Mer1, PCH2 is not. The SPO22 intron has a non-consensus 5′ splice site (5′SS) that dictates its Nam8/Mer1-dependence. SPO22 splicing relies on Mer1 recognition, via its KH domain, of an intronic enhancer 5′-AYACCCUY. Mutagenesis of KH and the enhancer highlights Arg214 and Gln243 and the CCC triplet as essential for Mer1 activity. The Nam8-dependent PCH2 pre-mRNA has a consensus 5′SS and lacks a Mer1 enhancer. For PCH2, a long 5′ exon and a non-consensus intron branchpoint dictate Nam8-dependence. Our results implicate Nam8 in two distinct meiotic splicing regulons. Nam8 is composed of three RRM domains, flanked by N-terminal leader and C-terminal tail segments. The leader, tail and RRM1 are dispensable for splicing meiotic targets and unnecessary for vegetative Nam8 function in multiple synthetic lethal genetic backgrounds. Nam8 activity is enfeebled by alanine mutations in the putative RNA binding sites of the RRM2 and RRM3 domains.

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

Genetic predisposition to type 1 diabetes (T1D) has been associated with a chromosome 11 locus centered on the proinsulin gene (INS) and with differential steady-state levels of INS RNA from T1D-predisposing and -protective haplotypes. Here, we show that the haplotype-specific expression is determined by INS variants that control the splicing efficiency of intron 1. The adenine allele at IVS1-6 (rs689), which rapidly expanded in modern humans, renders the 3' splice site of this intron more dependent on the auxiliary factor of U2 small nuclear ribonucleoprotein (U2AF). This interaction required both zinc fingers of the 35-kD U2AF subunit (U2AF35) and was associated with repression of a competing 3' splice site in INS exon 2. Systematic mutagenesis of reporter constructs showed that intron 1 removal was facilitated by conserved guanosine-rich enhancers and identified additional splicing regulatory motifs in exon 2. Sequencing of intron 1 in primates revealed that relaxation of its 3' splice site in Hominidae coevolved with the introduction of a short upstream open reading frame, providing a more efficient coupled splicing and translation control. Depletion of SR proteins 9G8 and transformer-2 by RNA interference was associated with exon 2 skipping whereas depletion of SRp20 with increased representation of transcripts containing a cryptic 3' splice site in the last exon. Together, these findings reveal critical interactions underlying the allele-dependent INS expression and INS-mediated risk of T1D and suggest that the increased requirement for U2AF35 in higher primates may hinder thymic presentation of autoantigens encoded by transcripts with weak 3' splice sites.