The branchpoint residue is recognized during commitment complex formation before being bulged out of the U2 snRNA-pre-mRNA duplex

Aberrant splicing contributes to severe α-spectrin-linked congenital hemolytic anemia

The etiology of severe hemolytic anemia in most patients with recessive hereditary spherocytosis (rHS) and the related disorder hereditary pyropoikilocytosis (HPP) is unknown. Whole exome sequencing of DNA from probands of 24 rHS or HPP kindreds identified numerous mutations in erythrocyte membrane α-spectrin (SPTA1). Twenty-eight mutations were novel, with null alleles frequently found in trans to missense mutations. No mutations were identified in a third of SPTA1 alleles (17/48). Whole genome sequencing revealed linkage disequilibrium between the common rHS-linked α-spectrinBug Hill polymorphism and a rare intron 30 variant in all 17 mutation-negative alleles. In vitro minigene studies and in vivo splicing analyses revealed the intron 30 variant changes a weak alternate branch point (BP) to a strong BP. This change leads to increased utilization of an alternate 3' splice acceptor site, perturbing normal α-spectrin mRNA splicing and creating an elongated mRNA transcript. In vivo mRNA stability studies revealed the newly created termination codon in the elongated transcript activates nonsense mediated decay leading to spectrin deficiency. These results demonstrate a unique mechanism of human genetic disease contributes to the etiology of a third of cases of rHS, facilitating diagnosis and treatment of severe anemia, and identifying a new target for therapeutic manipulation.

A sequence-based, deep learning model accurately predicts RNA splicing branchpoints

Experimental detection of RNA splicing branchpoints is difficult. To date, high-confidence experimental annotations exist for 18% of 3' splice sites in the human genome. We develop a deep-learning-based branchpoint predictor, LaBranchoR, which predicts a correct branchpoint for at least 75% of 3' splice sites genome-wide. Detailed analysis of cases in which our predicted branchpoint deviates from experimental data suggests a correct branchpoint is predicted in over 90% of cases. We use our predicted branchpoints to identify a novel sequence element upstream of branchpoints consistent with extended U2 snRNA base-pairing, show an association between weak branchpoints and alternative splicing, and explore the effects of genetic variants on branchpoints. We provide genome-wide branchpoint annotations and in silico mutagenesis scores at http://bejerano.stanford.edu/labranchor.

Competition between the ATPase Prp5 and branch region-U2 snRNA pairing modulates the fidelity of spliceosome assembly

ATPase-facilitated steps during spliceosome function have been postulated to afford opportunities for kinetic proofreading. Spliceosome assembly requires the ATPase Prp5p, whose activity might thus impact fidelity during initial intron recognition. Using alanine mutations in S. cerevisiae Prp5p, we identified a suboptimal intron whose splicing could be improved by altered Prp5p activity and then, using this intron, screened for potent prp5 mutants. These prp5 alleles specifically alter branch region selectivity, with improved splicing in vivo of suboptimal substrates correlating with reduced ATPase activity in vitro for a series of mutants in ATPase motif III (SAT). Because these effects are abrogated by compensatory U2 snRNA mutations or other changes that increase branch region-U2 pairing, these results explicitly link a fidelity event with a defined physical structure, the branch region-U2 snRNA duplex, and provide strong evidence that progression of the splicing pathway requires branch region-U2 snRNA pairing prior to Prp5p-facilitated conformational change.

The intron 1 of HPV 16 has a suboptimal branch point at a guanosine

The branch point sequence (BPS) of intron 1 of the HPV-16 was determined via RT-PCR in a cell free system, using lariat intermediates obtained by in vitro splicing reactions. We used synthetic E6/E7 transcripts and HeLa nuclear protein extracts to obtain the splicing intermediates. Then, a divergent oligonucleotide primer set, pairing on the lariat RNA that encompassed the 2'-5' phosphodiester bond formed between the 5' end of the intron and the BPS, was used for cDNA synthesis and PCR amplification. Subsequent RT-PCR assays revealed four splicing intermediates, made up of a major intermediary corresponding to the BPS and four cryptic branched sequences. Only intermediates bound at the 5' end of the intron are probably the authentic branch point sequence, and all of them branch at guanosine 328 instead of the typical adenosine. Unusually, the BPS of intron 1 of HPV-16 is a suboptimal sequence (AGUGAGU) that differs from the eukaryotic consensus BPS, which correlates with the splicing profile observed for early transcripts of HPV-16 in tumors and tumor derived cell lines. The implications of this unusual branch point sequence for splicing of the HPV-16 pre-mRNA are discussed.

Cotranscriptional spliceosome assembly dynamics and the role of U1 snRNA:5'ss base pairing in yeast

To investigate the mechanism of spliceosome assembly in vivo, we performed chromatin immunoprecipitation (ChIP) analysis of U1, U2, and U5 small nuclear ribonucleoprotein particles (snRNPs) to intron-containing yeast (S. cerevisiae) genes. The snRNPs display patterns that indicate a cotranscriptional assembly model: U1 first, then U2, and the U4/U6*U5 tri-snRNP followed by U1 destabilization. cis-splicing mutations also support a role of U2 and/or the tri-snRNP in U1 destabilization. Moreover, they indicate that splicing efficiency has a major impact on cotranscriptional snRNP recruitment and suggest that cotranscriptional recruitment of U2 or the tri-snRNP is required to commit the pre-mRNA to splicing. Branchpoint (BP) mutations had a major effect on the U1 pattern, whereas 5' splice site (5'ss) mutations had a stronger effect on the U2 pattern. A 5'ss-U1 snRNA complementation experiment suggests that pairing between U1 and the 5'ss occurs after U1 recruitment and contributes to a specific U1:substrate conformation required for efficient U2 and tri-snRNP recruitment.

A novel intronic cis element, ISE/ISS-3, regulates rat fibroblast growth factor receptor 2 splicing through activation of an upstream exon and repression of a downstream exon containing a noncanonical branch point sequence

Mutually exclusive splicing of fibroblast growth factor receptor 2 (FGFR2) exons IIIb and IIIc yields two receptor isoforms, FGFR2-IIIb and -IIIc, with distinctly different ligand binding properties. Several RNA cis elements in the intron (intron 8) separating these exons have been described that are required for splicing regulation. Using a heterologous splicing reporter, we have identified a new regulatory element in this intron that confers cell-type-specific inclusion of an unrelated exon that mirrors its ability to promote cell-type-specific inclusion of exon IIIb. This element promoted inclusion of exon IIIb while at the same time silencing exon IIIc inclusion in cells expressing FGFR2-IIIb; hence, we have termed this element ISE/ISS-3 (for "intronic splicing enhancer-intronic splicing silencer 3"). Silencing of exon IIIc splicing by ISE/ISS-3 was shown to require a branch point sequence (BPS) using G as the primary branch nucleotide. Replacing a consensus BPS with A as the primary branch nucleotide resulted in constitutive splicing of exon IIIc. Our results suggest that the branch point sequence constitutes an important component that can contribute to the efficiency of exon definition of alternatively spliced cassette exons. Noncanonical branch points may thus facilitate cell-type-specific silencing of regulated exons by flanking cis elements.

Allosteric cascade of spliceosome activation

Introns are removed from precursor messenger RNAs in the cell nucleus by a large ribonucleoprotein complex called the spliceosome. The spliceosome contains five subcomplexes called snRNPs, each with one RNA and several protein components. Interactions of the snRNPs with each other and the intron are highly dynamic, changing in an ordered progression throughout the splicing process. This allosteric cascade of interactions is programmed into the RNA and protein components of the spliceosome, and is driven by a family of DExD/H-box RNA-dependent ATPases. The dependence of cascade progression on multiple intron-recognition events likely serves to enforce the accuracy of splicing. Here, the progression of the allosteric cascade from the first recognition event to the first catalytic step of splicing is reviewed.

Mutational analysis of a plant branchpoint and polypyrimidine tract required for constitutive splicing of a mini-exon

The branchpoint sequence and associated polypyrimidine tract are firmly established splicing signals in vertebrates. In plants, however, these signals have not been characterized in detail. The potato invertase mini-exon 2 (9 nt) requires a branchpoint sequence positioned around 50 nt upstream of the 5' splice site of the neighboring intron and a U11 element found adjacent to the branchpoint in the upstream intron (Simpson et al., RNA, 2000, 6:422-433). Utilizing the sensitivity of this plant splicing system, these elements have been characterized by systematic mutation and analysis of the effect on inclusion of the mini-exon. Mutation of the branchpoint sequence in all possible positions demonstrated that branchpoints matching the consensus, CURAY, were most efficient at supporting splicing. Branchpoint sequences that differed from this consensus were still able to permit mini-exon inclusion but at greatly reduced levels. Mutation of the downstream U11 element suggested that it functioned as a polypyrimidine tract rather than a UA-rich element, common to plant introns. The minimum sequence requirement of the polypyrimidine tract for efficient splicing was two closely positioned groups of uridines 3-4 nt long (<6 nt apart) that, within the context of the mini-exon system, required being close (<14 nt) to the branchpoint sequence. The functional characterization of the branchpoint sequence and polypyrimidine tract defines these sequences in plants for the first time, and firmly establishes polypyrimidine tracts as important signals in splicing of at least some plant introns.

Sequences upstream of the branch site are required to form helix II between U2 and U6 snRNA in a trans-splicing reaction

Three different base paired stems form between U2 and U6 snRNA over the course of the mRNA splicing reaction (helices I, II and III). One possible function of U2/U6 helix II is to facilitate subsequent U2/U6 helix I and III interactions, which participate directly in catalysis. Using an in vitro trans-splicing assay, we investigated the function of sequences located just upstream from the branch site (BS). We find that these upstream sequences are essential for stable binding of U2 to the branch region, and for U2/U6 helix II formation, but not for initial U2/BS pairing. We also show that non-functional upstream sequences cause U2 snRNA stem-loop IIa to be exposed to dimethylsulfate modification, perhaps reflecting a U2 snRNA conformational change and/or loss of SF3b proteins. Our data suggest that initial binding of U2 snRNP to the BS region must be stabilized by an interaction with upstream sequences before U2/U6 helix II can form or U2 stem-loop IIa can participate in spliceosome assembly.

Characterization of the effects of mutations in the putative branchpoint sequence of intron 4 on the splicing within the human lecithin:cholesterol acyltransferase gene

We have previously identified a point mutation (intervening sequence (IVS) 4: T --> C) in the branchpoint consensus sequence of intron 4 of the lecithin:cholesterol acyltransferase (LCAT) gene in patients with fish-eye disease. To investigate the possible mechanisms responsible for the defective splicing, we made a series of mutations in the branchpoint sequence and expressed these mutants in HEK-293 cells followed by the analysis of pre-mRNA splicing using reverse transcriptase-polymerase chain reaction as well as LCAT activity assay. The results reveal that 1) the mutation of the branchpoint adenosine to any other nucleotide completely abolishes splicing; 2) the insertion of a normal branch site into the intronic sequence of the natural (IVS4-22c) or the branchpoint (IVS4-20t) mutant completely restores splicing; 3) the natural mutation can be partially rescued by making a single nucleotide change (G --> A) within the branchpoint consensus sequence; and 4) other single base changes, particularly around the branchpoint adenosine residue, significantly decrease the efficiency of splicing and thus enzyme activity. Surprisingly, the nucleotide transversion at the last position of the branchpoint sequence (i.e. IVS4-25a or -25g) results in a 2.7-fold increase in splicing efficiency. Therefore, these observations clearly establish the functional significance of the branchpoint sequence of intron 4 for the splicing of the human LCAT mRNA precursors.

Characterization of human RNA splice signals by iterative functional selection of splice sites

An iterative in vitro splicing strategy was employed to select for optimal 3' splicing signals from a pool of pre-mRNAs containing randomized regions. Selection of functional branchpoint sequences in HeLa cell nuclear extract yielded a sequence motif that evolved from UAA after one round of splicing toward a UACUAAC consensus after seven rounds. A significant part of the selected sequences contained a conserved AAUAAAG motif that proved to be functional both as a polyadenylation signal and a branch site in a competitive manner. Characterization of the branchpoint in these clones to either the upstream or downstream adenosines of the AAUAAAG sequence revealed that the branching process proceeded efficiently but quite promiscuously. Surprisingly, the conserved guanosine, adjacent to the common AAUAAA polyadenylation motif, was found to be required only for polyadenylation. In an independent experiment, sequences surrounding an optimal branchpoint sequence were selected from two randomized 20-nt regions. The clones selected after six rounds of splicing revealed an extended polypyrimidine tract with a high frequency of UCCU motifs and a highly conserved YAG sequence in the extreme 3' end of the randomized insert. Mutating the 3' terminal guanosine of the intron strongly affects complex A formation, implying that the invariant AG is recognized early in spliceosome assembly.

A limited number of pseudouridine residues in the human atac spliceosomal UsnRNAs as compared to human major spliceosomal UsnRNAs

Two forms of spliceosomes were found in higher eukaryotes. The major form contains the U1, U2, U4, U5, and U6 snRNAs; the minor form contains the U11, U12, U4atac, U5, and U6atac snRNAs. Assembly and function of the major form are based on a complex dynamic of UsnRNA-UsnRNA and UsnRNA-pre-mRNA interactions, and the involved UsnRNA segments are highly posttranscriptionally modified in plants and vertebrates. To further characterize the minor form of spliceosomes, we looked for the psi residues in HeLa cells' U11, U12, U4atac, and U6atac snRNAs, using chemical approaches. Four psi residues were detected in total for these four atac UsnRNAs, compared to 20 in their counterparts of the major spliceosomes. The two psi residues detected in U12 are also found in U2 snRNA. One of them belongs to the branch-site-recognition sequence. It forms one of the base pairs that bulge out the A residue, responsible for the nucleophilic attack. Conservation of this strategic psi residue probably reflects a functional role. Another psi residue was detected in a U4atac snRNA segment involved in formation of helix II with U6atac. The fourth one was detected in the additional stem-loop structure present at the 3' end of U6atac snRNA. Differences in psi content of the atac and major UsnRNAs of human cells may participate in the differentiation of the two splicing systems. Based on secondary structure similarity, U2 and U12 snRNAs on the one hand and U4 and U4atac snRNAs on the other hand may share common psi synthases.

Luc7p, a novel yeast U1 snRNP protein with a role in 5' splice site recognition

The characterization of a novel yeast-splicing factor, Luc7p, is presented. The LUC7 gene was identified by a mutation that causes lethality in a yeast strain lacking the nuclear cap-binding complex (CBC). Luc7p is similar in sequence to metazoan proteins that have arginine-serine and arginine-glutamic acid repeat sequences characteristic of a family of splicing factors. We show that Luc7p is a component of yeast U1 snRNP and is essential for vegetative growth. The composition of yeast U1 snRNP is altered in luc7 mutant strains. Extracts of these strains are unable to support any of the defined steps of splicing unless recombinant Luc7p is added. Although the in vivo defect in splicing wild-type reporter introns in a luc7 mutant strain is comparatively mild, splicing of introns with nonconsensus 5' splice site or branchpoint sequences is more defective in the mutant strain than in wild-type strains. By use of reporters that have two competing 5' splice sites, a loss of efficient splicing to the cap proximal splice site is observed in luc7 cells, analogous to the defect seen in strains lacking CBC. CBC can be coprecipitated with U1 snRNP from wild-type, but not from luc7, yeast strains. These data suggest that the loss of Luc7p disrupts U1 snRNP-CBC interaction, and that this interaction contributes to normal 5' splice site recognition.

A conditional U5 snRNA mutation affecting pre-mRNA splicing and nuclear pre-mRNA retention identifies SSD1/SRK1 as a general splicing mutant suppressor

A combination of point mutations disrupting both stem 1 and stem 2 of U5 snRNA (U5AI) was found to confer a thermosensitive phenotype in vivo. In a strain expressing U5AI, pre-mRNA splicing was blocked before the first step through an inability of the mutant U5 snRNA to efficiently associate with the U4/U6 di-snRNP. Formation of early splicing complexes was not affected in extracts prepared from U5 snRNA mutant cells, while the capacity of these extracts to splice a pre-mRNA in vitro was greatly diminished. In addition, significant levels of a translation product derived from intron containing pre-mRNAs could be detected in vivo. The SSD1/SRK1 gene was identified as a multi-copy suppressor of the U5AI snRNA mutant. Single copy expression of SSD1/SRK1 was sufficient to suppress the thermosensitive phenotype, and high copy expression partially suppressed the splicing and U4/U6.U5 tri-snRNP assembly pheno-types. SSD1/SRK1 also suppressed thermosensitive mutations in the Prp18p and U1-70K proteins, while inhibiting growth of the cold sensitive U1-4U snRNA mutant at 30 degrees C. Thus we have identified SSD1/SRK1 as a general suppressor of splicing mutants.

Transient interaction of BBP/ScSF1 and Mud2 with the splicing machinery affects the kinetics of spliceosome assembly

Removal of introns from pre-mRNA is an essential step of gene expression. The splicing reaction is catalyzed in a large complex termed the spliceosome. Introns are recognized during the early steps of spliceosome assembly with the formation of commitment complexes. Intron recognition is mediated by the interaction of splicing factors with conserved sequences present in the pre-mRNA. BBP/SF1 participates in this recognition by interacting with the pre-mRNA branch point in both yeast and mammals. This protein, which is essential in yeast, also interacts with the U2AF65/Mud2 splicing factor. However, its precise role in splicing complex formation is still unclear. We have now analyzed the presence of BBP and Mud2 in yeast splicing complexes using supershift and coprecipitation assays. We found that BBP is present together with Mud2 in commitment complex 2 (CC2), but is not detectable in commitment complex 1 (CC1). Furthermore, genetic and biochemical depletion of BBP demonstrated that it is required for CC2 formation. In addition we observed that BBP and Mud2 are not detectable in pre-spliceosomes. These are the first commitment complex components that are shown to be released during or immediately after pre-spliceosome formation. Interestingly, depletion of BBP or disruption of MUD2 had no significant effect on pre-spliceosome formation and splicing in vitro but led to a transient accumulation of CC1. These observations support a model in which BBP and Mud2 are recycled during transition from CC2 to pre-spliceosome.

Pseudouridine mapping in the Saccharomyces cerevisiae spliceosomal U small nuclear RNAs (snRNAs) reveals that pseudouridine synthase pus1p exhibits a dual substrate specificity for U2 snRNA and tRNA

Pseudouridine (Psi) residues were localized in the Saccharomyces cerevisiae spliceosomal U small nuclear RNAs (UsnRNAs) by using the chemical mapping method. In contrast to vertebrate UsnRNAs, S. cerevisiae UsnRNAs contain only a few Psi residues, which are located in segments involved in intermolecular RNA-RNA or RNA-protein interactions. At these positions, UsnRNAs are universally modified. When yeast mutants disrupted for one of the several pseudouridine synthase genes (PUS1, PUS2, PUS3, and PUS4) or depleted in rRNA-pseudouridine synthase Cbf5p were tested for UsnRNA Psi content, only the loss of the Pus1p activity was found to affect Psi formation in spliceosomal UsnRNAs. Indeed, Psi44 formation in U2 snRNA was abolished. By using purified Pus1p enzyme and in vitro-produced U2 snRNA, Pus1p is shown here to catalyze Psi44 formation in the S. cerevisiae U2 snRNA. Thus, Pus1p is the first UsnRNA pseudouridine synthase characterized so far which exhibits a dual substrate specificity, acting on both tRNAs and U2 snRNA. As depletion of rRNA-pseudouridine synthase Cbf5p had no effect on UsnRNA Psi content, formation of Psi residues in S. cerevisiae UsnRNAs is not dependent on the Cbf5p-snoRNA guided mechanism.

Interaction of the U1 snRNP with nonconserved intronic sequences affects 5' splice site selection

Intron definition and splice site selection occur at an early stage during assembly of the spliceosome, the complex mediating pre-mRNA splicing. Association of U1 snRNP with the pre-mRNA is required for these early steps. We report here that the yeast U1 snRNP-specific protein Nam8p is a component of the commitment complexes, the first stable complexes assembled on pre-mRNA. In vitro and in vivo, Nam8p becomes indispensable for efficient 5' splice site recognition when this process is impaired as a result of the presence of noncanonical 5' splice sites or the absence of a cap structure. Nam8p stabilizes commitment complexes in the latter conditions. Consistent with this, Nam8p interacts with the pre-mRNA downstream of the 5' splice site, in a region of nonconserved sequence. Substitutions in this region affect splicing efficiency and alternative splice site choice in a Nam8p-dependent manner. Therefore, Nam8p is involved in a novel mechanism by which a snRNP component can affect splice site choice and regulate intron removal through its interaction with a nonconserved sequence. This supports a model where early 5' splice recognition results from a network of interactions established by the splicing machinery with various regions of the pre-mRNA.

A role for U2/U6 helix Ib in 5' splice site selection

Selection of pre-mRNA splice sites is a highly accurate process involving many trans-acting factors. Recently, we described a role for U6 snRNA position G52 in selection of the first intron nucleotide (+1G). Because some U2 alleles suppress U6-G52 mutations, we investigated whether the corresponding U2 snRNA region also influenced 5' splice site selection. Our results demonstrate that U2 snRNAs mutated at position U23, but not adjacent nucleotides, specifically affect 5' splice site cleavage. Furthermore, all U2 position U23 mutations are synthetic lethal with the thermosensitive U6-G52U allele. Interestingly, the U2-U23C substitution has an unprecedented hyperaccurate splicing phenotype in which cleavage of introns with a +1G substitution is reduced, whereas the strain grows with wild-type kinetics. U2 position U23 forms the first base pair with U6 position A59 in U2/U6 helix Ib. Restoration of the helical structure suppresses 5' splice site cleavage defects, showing an important role for the helix Ib structure in 5' splice site selection. U2/U6 helix Ib and helix II have recently been described as being functionally redundant. This report demonstrates a unique role for helix Ib in 5' splice site selection that is not shared with helix II.

Multiple developmental requirements of noisette, the Drosophila homolog of the U2 snRNP-associated polypeptide SP3a60

We report the cloning of the noisette gene (noi), which encodes the Drosophila melanogaster ortholog of a U2 snRNP-associated splicing factor, SF3a60 (SAP61) in humans and PRP9p in Saccharomyces cerevisiae. Antibodies raised against human SF3a60 recognized NOI in flies, showing a nuclear localization in all the stages examined, including the embryo, the dividing cells of imaginal discs, and the larval polyploid nuclei. NOI is expressed in somatic and germinal cells of both male and female gonads. By mobilization of P transposons, we have generated a large number of noi mutations. Complete loss of function resulted in lethality at the end of embryogenesis, without obvious morphological defects. Hypomorphic alleles revealed multiple roles of noi for the survival and differentiation of male germ cells, the differentiation of female germ cells, and the development of several adult structures.

Genetic interaction between U6 snRNA and the first intron nucleotide in Saccharomyces cerevisiae

Nuclear pre-mRNA splicing necessitates specific recognition of the pre-mRNA splice sites. It is known that 5' splice site selection requires base pairing of U6 snRNA with intron positions 4-6. However, no factor recognizing the highly conserved 5' splice site GU has yet been identified. We have tested if the known U6 snRNA-pre-mRNA interaction could be extended to include the first intron nucleotides and the conserved 50GAG52 sequence of U6 snRNA. We observe that some combinations of 5' splice site and U6 snRNA mutations produce a specific synthetic block to the first splicing step. In addition, the U6-G52U allele can switch between two competing 5' splice sites harboring different nucleotides following the cleavage site. These results indicate that U6 snRNA position 52 interacts with the first nucleotide of the intron before 5' splice site cleavage. Some combinations of U6 snRNA and pre-mRNA mutations also blocked the second splicing step, suggesting a role for the corresponding nucleotides in a proofreading step before exon ligation. From studies in diverse organisms, various functions have been ascribed to the conserved U6 snRNA 47ACAGAG52 sequence. Our results suggest that these discrepancies might reflect variations between different experimental systems and point to an important conserved role of this sequence in the splicing reaction.

Construction of an in vivo-regulated U6 snRNA transcription unit as a tool to study U6 function

U6 snRNA is the only spliceosomal snRNA transcribed by RNA polymerase III in yeast. We have constructed a regulated U6 snRNA transcription unit by introducing the binding site for the Escherichia coli lacI repressor protein in the U6 snRNA promoter. GAL-induced expression of lacI protein led to a decrease in U6 snRNA levels and blocked cell growth. lacI dissociation from the promoter, and consequent U6 snRNA transcription, could be induced by addition of IPTG and repression of lacI transcription. To test the usefulness of this system in studying spliceosomal U6 snRNA function, we conditionally expressed U6 snRNAs with a single base substitution in position A51. We demonstrate that expression of the U6-A51 mutations confers a strong dominant negative phenotype as shown by severe reductions in growth rate. In these strains, splicing of endogenous pre-mRNAs was blocked before the second step.