Abstract
Pre-mRNA splicing is catalyzed in two steps: 5’ splice site (SS) cleavage and exon ligation. A number of proteins transiently associate with spliceosomes to specifically impact these steps (1st and 2nd step factors). We recently identified Fyv6 (FAM192A in humans) as a 2nd step factor in S. cerevisiae; however, we did not determine how widespread Fyv6’s impact is on the transcriptome. To answer this question, we have used RNA-seq to analyze changes in splicing. These results show that loss of Fyv6 results in activation of non-consensus, branch point (BP) proximal 3’ SS transcriptome-wide. To identify the molecular basis of these observations, we determined a high-resolution cryo-EM structure of a yeast product complex spliceosome containing Fyv6 at 2.3 Å. The structure reveals that Fyv6 is the only 2nd step factor that contacts the Prp22 ATPase and that Fyv6 binding is mutually exclusive with that of the 1st step factor Yju2. We then use this structure to dissect Fyv6 functional domains and interpret results of a genetic screen for fyv61 suppressor mutations. The combined transcriptomic, structural, and genetic studies allow us to propose a model in which Yju2/Fyv6 exchange facilitates exon ligation and Fyv6 promotes usage of consensus, BP distal 3’ SS.
eLife assessment
This important study addresses how 3' splice site choice is modulated by the conserved spliceosome-associated protein Fyv6. The authors provide compelling evidence Fyv6 functions to enable selection of 3' splice sites distal to a branch point and in doing so antagonizes more proximal, suboptimal 3' splice sites. The study would be improved through a more nuanced discussion of alternative possibilities and models, for instance in discussing the phenotypic impact of Fyv6 deletion.
Introduction
Precursor messenger RNA (pre-mRNA) splicing is catalyzed by a large macromolecular complex called the spliceosome. Spliceosomes are composed of 5 small nuclear ribonucleoproteins (snRNPs), each composed of a small nuclear RNA (snRNA; U1, U2, U4, U5, and U6) and several different protein splicing factors. The snRNPs and dozens of other proteins assemble de novo at each intron to form spliceosomes. Spliceosomes are highly dynamic and form different complexes as proteins and snRNAs join and leave. Many of these complexes have been characterized biochemically, genetically, and by cryo-electron microscopy (cryo-EM) (Plaschka et al., 2019). The splicing reaction itself is carried out in two sequential transesterification reactions (Fig. 1A). First, the 5’ splice site (5’ SS) is cleaved during formation of an intron lariat (1st step). Second, the intron lariat is removed simultaneously with exon ligation at the 3’ SS (2nd step).
![](https://prod--epp.elifesciences.org/iiif/2/100449%2Fv1%2Fcontent%2F592262v2_fig1.tif/full/max/0/default.jpg)
Widespread splicing changes due to loss of Fyv6.
A) Pre-mRNA splicing reaction steps and products. B) Growth of upf11 and fyv61upf11 yeast at 16, 23, 30, and 37°C. Images were taken on the indicated days. C) Number of alternative splicing events observed in a fyv61upf11 strain relative to a strain containing FYV6 after a temperature shift to 16, 30, and 37°C. D) Venn diagram of shared alternative splicing events at 16 (blue), 30 (green), and 37°C (red). E) Box plots of the log10 of the ratios of fyv61 upf11 to WT (upf11) FAnS for RPGs (red) and non-RPGs (black).
The integrity of the genetic information contained within a pre-mRNA depends on correct identification of the 5’ and 3’ SS by the splicing machinery since a single nucleotide shift in either site could destroy a protein reading frame. In yeast, several DExD/H-box ATPases function to limit usage of suboptimal SS by enhancing splicing fidelity (Chung et al., 2023; Semlow and Staley, 2012). Two of these ATPases, Prp16 and Prp22, impact the fidelity of the 1st and 2nd catalytic steps, respectively. In addition to these roles, Prp16 is also required for remodeling of the spliceosome to permit the 1st-to-2nd step transition (Schwer and Guthrie, 1992), and Prp22 is essential for releasing the mRNA product (Company et al., 1991). While Prp16 does not need to be present during the 1st step (Chung et al., 2023), Prp22 is required to be present for exon ligation if the branch point (BP) to 3’ SS distance is > 21 nucleotides (nt) (Schwer and Gross, 1998). Despite many biochemical, single molecule, and structural studies, how spliceosomes promote usage of BP distal, Prp22-dependent 3’ SS has remained elusive.
Identification of the correct 3’ SS is an especially challenging problem given that the consensus sequence in Saccharomyces cerevisiae (yeast) and humans is just three nucleotides (YAG, Y=U or C). Cryo-EM structures of spliceosome product (P) complexes have revealed how these nucleotides can be recognized within the spliceosome active site (Bai et al., 2017; Fica et al., 2019; Liu et al., 2017; Wilkinson et al., 2017; Zhang et al., 2019). The AG 3’ SS dinucleotide is recognized by formation of non-Watson Crick base-pairing interactions with the 5’ SS +1G and BP adenosine, while the pyrimidine at the -3 position of the 3’ SS is recognized by the Prp8 protein. Given this short consensus sequence, multiple different splicing factors also help to ensure that the correct 3’ SS is utilized, including Cwc21 and the 2nd step factors Slu7, Prp18, and Prp22 (Crotti et al., 2007; Frank and Guthrie, 1992; Gautam et al., 2015; Kawashima et al., 2014; Roy et al., 2023; Semlow et al., 2016). It is believed that Slu7 and Prp18 help to stabilize 3’ SS docking to the active site and that this contributes to 2nd step efficiency. It has been proposed that Prp22 antagonizes 3’ SS docking to permit sampling of different 3’ SS to ensure that optimal sequences are used as a proofreading mechanism (Mayas et al., 2006; Semlow et al., 2016). Importantly, identification of the correct 3’ SS depends not just on identification of the YAG sequence but also involves choosing which YAG sequence to use. Prp18 appears to aid in selection of BP distal 3’ SS and avoidance of BP proximal (and often non-consensus) 3’ SS. This observation could be due to Prp18 imposing a BP to 3’ SS distance constraint, enforcing use of canonical YAG 3’ SS sequences, or a combination of the two activities.
Recently, Fyv6 was identified as a novel 2nd step factor in yeast (Lipinski et al., 2023). Fyv6 is a homolog of the human protein FAM192A, which was identified in pre-C* structures of the human spliceosome (Zhan et al., 2022). While the impact of FAM192A on human spliceosome activity has not yet been well-characterized, loss of Fyv6 in yeast decreases 2nd step splicing efficiency in vitro and results in use of an alternative 3’ SS in the SUS1 pre-mRNA in vivo. It is likely that unassigned cryo-EM densities in yeast C* (just prior to exon ligation) and P complexes correspond to Fyv6 since these densities are in analogous locations to that for FAM192A in the human pre-C* structure (Zhan et al., 2022). However, the resolutions of the yeast C* and P complex structures make unambiguous assignment of Fyv6 amino acid side chains and interactions difficult.
Here, we use a combination of transcriptomic, structural, biochemical, and genetic assays to elucidate how Fyv6 controls 3’ SS usage in yeast. RNA-seq analysis reveals widespread activation of alternative, BP proximal nonconsensus 3’ SS in the absence of Fyv6 consistent with biochemical assays that show Fyv6 facilitates splicing of 3’ SS located > 21 nt from the BP. To further elucidate Fyv6 function, we determined a 2.3 Å cryo-EM structure of the yeast P complex spliceosome in which Fyv6 can be modeled. The structure reveals interactions between Fyv6 and Prp22 as well as a domain in Syf1 that interacts with either Fyv6 or the 1st step factor Yju2. 3D classification allowed us to identify two additional conformational states of P complex and provides insights into the coupling of 3’ SS active site docking and the presence of 2nd step factors. We then used this structural data to dissect the functional domains of Fyv6, interpret the results of a genetic screen for suppressors of fyv6Δ, and to probe Fyv6 and Prp22 interactions. Combined, we propose a model in which Fyv6 is recruited to the spliceosome via interactions with the NTC component Syf1 in order to promote usage of BP distal, and Prp22-dependent, 3’ SS.
Results
Deletion of FYV6 Results in Widespread Use of Alternative 3′ Splice Sites
Our previous work showed that loss of FYV6 causes changes in 3’ SS usage in the SUS1 pre-mRNA. To probe for changes transcriptome-wide, we used nonsense mediated decay (NMD)-deficient (upf11, WT) yeast cells expressing or lacking Fyv6. By suppressing NMD, we hoped to limit degradation of alternatively spliced mRNA isoforms generated in the absence of Fyv6 (Sayani et al., 2008). Consistent with results from fyv61 strains, the fyv61 upf11 double mutant strain also showed reduced growth at 30°C and cold-sensitive (cs) and temperature-sensitive (ts) phenotypes relative to a upf11 control strain (Fig. 1B). We then used RNA-Seq to analyze the isolated RNAs from these strains. As expected, many more changes in splicing were observed in the fyv61 upf11 strain relative to the fyv61 strain (Supplemental Fig. S1A, B; Supplemental Tables S1, S2). This is consistent with many mRNA isoforms generated due to loss of Fyv6 being substrates for NMD. Loss of Fyv6 resulted in a number of changes in gene expression with many non-intron containing genes also being up- or down-regulated (Supplemental Fig. S1C).
We compared and classified changes in splicing between strains containing or lacking FYV6 under permissive growth conditions (30°C) or after the yeast had been shifted to a non-permissive temperature (16 or 37°C) for 1 h (Fig. 1C). In this case, we only considered events resulting in at least a 10% change in the percent spliced in (PSI) value when FYV6 is deleted. At all temperatures, we see changes in alternative 3’ SS usage in the fyv61 strain relative to the FYV6 background. For example, we detected use of the alternative 3’ SS in the first intron of the SUS1 transcript as we previously reported (Supplemental Fig. S1D; Lipinski et al., 2023). We also saw increased use of a cryptic, nonconsensus GAG 3’ SS in the first, but not second, intron of the YOS1 transcript in the fyv61 data sets and confirmed this result by RT-PCR (Supplemental Fig. S1E, F). We observed several cases of alternative 5’ SS. However, the majority of these alternative 5’ SS were used with equal efficiencies in the WT and fyv61 strains (Supplemental Fig. S1G).
The strains temperature shifted to 16°C had the largest number of changes due to FYV6 deletion and the strains shifted to 37°C the least. The temperature shift to 16° caused an increase in the numbers of both alternative 3’ SS and retained intron events detected (Fig. 1C). The changes in alternative 3’ SS used at each temperature were mostly unique with few events detected under all three conditions (Fig. 1D). Collectively, we observed that 61 different RNAs changed splicing patterns due to loss of Fyv6. For many of these mRNAs, we detected multiple different 3’ SS being used.
To investigate these splicing events more closely, we analyzed changes in usage of alterative splice junctions by calculating the Fraction of Annotated Splicing (FAnS) in the presence and absence of Fyv6 (Roy et al., 2023). FAnS ratios report the relative abundance of an alternative splicing event in relation to the main spliced isoform. We calculated FAnS values for all exon junction reads corresponding to usage of a canonical 5’ SS and an alternative 3’ SS for WT and fyv6Δ strains at each temperature. We then calculated the ratios of the fyv6Δ/WT FAnS values for each transcript. These results confirm changes in alternative 3’ SS usage for both ribosomal protein gene (RPG) and non-RPG transcripts at all conditions (Fig. 1E). We could not detect any significant differences between RPG and non-RPG splicing with this analysis, suggesting that differences in gene expression and splicing efficiencies between these two classes are not correlated with alternative 3’ SS usage. Combined these results show that Fyv6 plays a critical role in defining mRNA isoform production in yeast, particularly under non-optimal growth conditions such as temperature stress.
Fyv6 Facilitates Usage of Branch Point Distal 3′ Splice Sites
We next analyzed the features of the 3’ SS impacted by Fyv6 loss. We used RNAs collected from yeast grown at 30°C and sequenced the libraries at higher depth than those used for the temperature-shift experiment described in Fig. 1 (∼400 million vs. 200 million reads; Supplemental Table S2). We first identified alternative 3’ SS used in each strain, annotated each as being upstream (5’) or downstream (3’) of the canonical 3’ SS (Grate and Ares, 2002), and then calculated the FAnS ratio for each site. This analysis shows that the majority of the alternative 3’ SS arising from Fyv6 deletion are upstream (Fig. 2A, pink points above the diagonal). When these alternative 3’ SS are mapped according to their nucleotide distance from the canonical site, most of these are found within ∼40 nt upstream (Fig. 2B, purple violin). Plotting these alternative 3’ SS sites relative to their distance from the annotated BP reveals that the highest density of sites is < 20 nt from the BP for transcripts with the highest FAnS ratios (Fig. 2C, purple violin). It is important to note that we did not map BP in Fyv6 deletion strains, and it is possible that in some cases a shift in 3’ SS usage could also coincide with a shift in the BP. Indeed, a small number of alternative 3’ SS detected in the Fyv6 data sets are located upstream of the annotated BP, suggesting activation of an alternative BP and 3’ SS. However, we believe that the changes are predominantly due to usage of alternative 3’ SS located between the annotated BP and canonical 3’ SS.
![](https://prod--epp.elifesciences.org/iiif/2/100449%2Fv1%2Fcontent%2F592262v2_fig2.tif/full/max/0/default.jpg)
Loss of Fyv6 activates BP proximal, non-consensus 3’ SS.
A) Plot of the log10 of FAnS values for changes in alternative 3’ SS usage in fyv6Δ upf1Δ and upf1Δ strains. Points are colored based on whether they are upstream (pink) or downstream (blue) to the canonical 3’ SS. B) Violin plots of the ratios of fyv6Δupf1Δ to the upf1Δ FAnS values based on the distances between the alternative and canonical 3’ SS. FAnS ratios >0 are indicative of the site being upregulated in fyv6Δ. C) Violin plots of the ratios of fyv6Δupf1Δ to the upf1Δ FAnS values based on the distances between the alternative 3’ SS and the annotated BP. FAnS ratios >0 are indicative of the site being upregulated in fyv6Δ. D) Diagram of ACT1-CUP1 reporter showing BP-3’ SS distances. E) Representative primer extension analysis of RNA products generated from splicing of the ACT1-CUP1 reporter in the presence (WT) or absence of Fyv6 (fyv6Δ). Bands for fully spliced mRNA and lariat intermediate are indicated. U6 snRNA was detected as a loading control. The * indicates an unknown product present in every lane. F) Quantification of the primer extension results from N=3 replicates represented by the ratio of band intensities for mRNA/(mRNA+lariat intermediate). Bars represent the average ratio of the replicates ±SD. Means between WT and fyv6Δ groups for each reporter were compared with an unpaired Welch’s two-tailed t-test. Significance is indicated: n.s. no significance; ** p<0.01; *** p<0.001. G) Sequence logos of alternative 3’ SS with FAnS > 0 sorted by either upstream or downstream of the canonical 3’ SS compared against the canonical 3’ SS for genes with FanS > 0. H) The log10 of the FAnS ratio for alternative 3’ SS sorted by the sequences of the 3’ SS.
The above results suggest that Fyv6 is important for splicing of 3’ SS located distant from the BP. To test this directly and systematically, we created a series of ACT1-CUP1 reporters with various BP to 3’ SS distances (9-50 nt, Fig. 2D, Supplemental Table S3) based on sequences used in previous studies of 3’ SS selection (Brys and Schwer, 1996; Frank and Guthrie, 1992; Schwer and Gross, 1998). We transformed strains containing or lacking Fyv6 and without debranchase (Dbr1) to limit degradation of lariat intermediates. We then detected RNA products by primer extension and calculated the 2nd step (exon ligation) efficiency for each reporter (Fig. 2E, F). While loss of Fyv6 has minimal impact on exon ligation when BP to 3’ SS distances are short (9-15 nt), exon ligation efficiency decreases significantly at distances of 21 nt or greater. In addition, we do not see evidence of cryptic BP usage in these assays. This result is consistent with the RNA-seq analysis and supports a function for Fyv6 in facilitating splicing at BP distal 3’ SS.
Finally, we analyzed the sequence features of the alternative 3’ SS activated by loss of Fyv6 (Fig. 2G, H). For sites located upstream of the canonical 3’ SS, we observed much greater sequence variability relative to the canonical site. We detected increased use of highly variable 3’ SS including those with atypical HAU (H=A,C,U) and BG (B=C,U,G) motifs and confirmed generation of new RNA isoforms for several of these by RT-PCR (Fig. 2H, Supplemental Fig. S2). Increased use of variable 3’ SS was also observed when the 2nd step factor PRP18 was deleted (Roy et al., 2023). However, the subset of RNAs most impacted by splicing factor deletion appear to be different. For example, we were not able to detect changes in 3’ SS usage for the UBC12, MAF1, MUD1, PHO85, SPT14, or YCL002C transcripts in the fyv61 data sets as was reported for prp181 (not shown). Prp18 and Fyv6 likely have transcript-specific effects on splicing outcomes. Combined our results strongly support Fyv6 as a 2nd step splicing factor that facilitates usage of consensus sites located distal to the BP in vivo.
A High-Resolution Spliceosome Structure Reveals Fyv6 Interactions
Previous structures of the yeast P-complex spliceosome were solved with resolutions at the core ranging from 3.3 Å to 3.7 Å (Bai et al., 2017; Liu et al., 2017; Wilkinson et al., 2017), and with peripheral components such as the Prp22 helicase, U2 snRNP, and Prp19 Complex (NTC) at much lower resolutions (5-10 Å) that precluded detailed investigation. To improve the resolution, we purified P complex as previously described by stalling spliceosome disassembly with a dominant negative mutant of Prp22 (S635A) defective in mRNA release (Schwer and Meszaros, 2000; Wilkinson et al., 2017). However, by collecting a much larger cryo-EM dataset on a more modern electron detector, we were able to resolve the structure to 2.3 Å within the catalytic core, as well as visualize peripheral components at resolutions from 3.0 Å – 3.7 Å (Fig. 3A, B; Supplemental Fig. S3-S6; Supplemental Table S4). This is the highest resolution spliceosome structure to date. The high-quality density at the active site shows the positions of bases unambiguously, confirming the manner of the non-Watson-Crick base pairing that mediates 3’ SS recognition (Liu et al., 2017; Wilkinson et al., 2017). Monovalent ions bound at the active site were previously shown to be important during the 1st step of splicing (Wilkinson et al., 2021), and our data indicate that the positions of these ions are preserved after exon ligation. Additionally, the density also shows a metal ion, likely potassium, that bridges the base pair between +1G of the 5’ SS and -1G of the 3’ SS (Fig. 3C).
![](https://prod--epp.elifesciences.org/iiif/2/100449%2Fv1%2Fcontent%2F592262v2_fig3.tif/full/max/0/default.jpg)
Cryo-EM structure of the yeast P complex spliceosome at 2.3 Å resolution.
A) A composite density map for P complex showing focused refinements of the Prp22, NTC, U2 snRNP, U5 Sm ring, and Cwc22 N-terminal domain regions. B) Overall model for the P complex spliceosome. C) Cryo-EM density for state I of P complex (above, low-pass filtered) and for the active site and 3’ SS (below, sharpened). D,E) As for C) but for states II (D) and III (E).
The 2nd-step splicing factor Prp22 is resolved at 3.0 Å resolution, allowing atomic interpretation for the first time (Fig. 3B and later figures). Bases 13 – 21 of the 3’ exon are visible between the two RecA domains of Prp22 consistent with biochemical footprinting data (Schwer, 2008). The S635A mutation used to stall P complex is also visible; however, the origin of defective mRNA release due to this mutation is not apparent (Schwer and Meszaros, 2000). The loop containing S635A is in a similar conformation in this structure and in the structure of WT human Prp22 (Felisberto-Rodrigues et al., 2019) (Supplemental Fig. S6). Part of the extensive N-terminal domain of Prp22 is resolved and constitutes a long helix that bridges between the Prp22 RecA2 domain and the C-terminal domain of Cwc22; this helix was unassigned in previous P-complex structures (Fig. 4B).
![](https://prod--epp.elifesciences.org/iiif/2/100449%2Fv1%2Fcontent%2F592262v2_fig4.tif/full/max/0/default.jpg)
Fyv6 and its interactors in the P complex cryo-EM structure.
A) Cryo-EM density segmented around Fyv6. B) Structure of Prp22. C) Structure of the NTC within P complex. D) IP6 site and cryo-EM density in NTC. E) Interaction of the hook domain of Fyv6 with the Prp22 RecA2 domain. F) Interaction of Fyv6 with Slu7 and the region of the intron between the BP and 3’ SS. G) Interaction of Fyv6 with Syf1 and an analogous view of the interaction of Syf1 with Yju2 in C complex (Wilkinson et al., 2021).
To improve the resolution of the peripheral, flexible regions, we used a new data-driven regularization strategy (Kimanius et al., 2023) that allows focused refinements of much smaller domains than previously possible (Supplemental Fig. S4, S6). Combined with AlphaFold2-assisted modelling, we were able to obtain more accurate models for the U2 snRNP, NTC, U5 Sm ring, and Cwc22 N-terminal domains. Our model for the NTC within the P complex is similar to our previous model for the NTC within C and Ci complexes (Wilkinson et al., 2021) (Fig. 4C). Notably however, we observed a star-shaped density coordinated by several conserved lysines from Clf1 and Ntc20, which we were able to model as inositol hexakisphosphate (IP6) (Fig. 4D). The same density, at lower resolution, is visible in the C complex (Wilkinson et al., 2021). A distinct, separate IP6 molecule had previously been observed within all catalytic spliceosomes between the Bact and P-complex stages coordinated by Prp8 and (for C* and P-complexes) Slu7 (Supplemental Fig. S6). We do not believe that the NTC-bound IP6 is essential for splicing since deletion of the Clf1 tetratricopeptide repeat that coordinates IP6 does not result in growth or splicing phenotypes in yeast and Ntc20 is a non-essential splicing factor (Chen et al, 2001; Chung et al., 1999). Nonetheless, it is possible that IP6 may regulate splicing in some manner via interactions with both Prp8 and the NTC.
We were able to unambiguously assign three, connected long α-helices to Fyv6 (Fig. 4A). These helices were visible in previous cryo-EM reconstructions of yeast C* and P complexes but were either unassigned or misassigned. In the cases of misassigned densities, Fyv6 was previously attributed to the C-terminal domain of the 1st step factor Yju2 (Wilkinson et al., 2017). Yju2 also has an elongated helical architecture and adopts a similar position on C complex (Wilkinson et al., 2021) and the post-P complex Intron Lariat Spliceosome (ILS) (Wan et al., 2017). We closely inspected the density from previously determined spliceosome structures and concluded that the Yju2 C-terminus is still a better fit than Fyv6 for densities in B* and C complexes, whereas Fyv6 is a better fit than Yju2 for densities in C* and P complexes (Supplemental Fig. S7). Therefore, this position of Fyv6 seems to be characteristic of a 2nd step conformation of the spliceosome and consistent with the position of its distant human homolog FAM192A in the human pre-C* complex (Zhan et al., 2022).
Fyv6 makes a multitude of interactions with essential splicing factors (Fig. 3B; 4E, F, G). One conserved region of Fyv6 is the N-terminus with a FVSE motif (Supplemental Fig. S8). This motif forms a ‘hook’ that interacts with a hydrophobic patch on the surface of the Prp22 RecA2 domain, where it also forms hydrogen bonds by β-sheet augmentation (Fig. 4E). The hook is followed by three α-helices. The start of helix 2 sits on top of the intron between the BP and the active site docked 3’ SS, where it may act as a steric block to prevent 3’ SS undocking and promote exon ligation (Fig. 4F) (Liu et al., 2017). Helix 2 also contains a patch of conserved arginines that form salt bridges with Asp317 and Glu318 of Slu7 (Fig. 4F). The linker between helices 2 and 3 interacts with the NTC component Cef1, and helix 3 interacts with NTC component Syf1 (Fig. 4G). The binding interface with Syf1 overlaps with the Syf1/Yju2 interaction in the C complex (Wilkinson et al., 2021), suggesting these interactions are mutually exclusive (Fig. 4G). Fyv6 thus provides a direct link between the Prp22 ATPase, NTC components, and other essential splicing factors over a distance of ∼100 Å.
By 3D classification, we were able to resolve a total of three states of the P-complex spliceosome (Fig. 3C-E). State I, represented by 51% of particles, is the most complete and has all the above-described factors visible, including a stably-docked 3’ SS. State II, represented by 20% of particles, is similar to State I but lacks density for the N-terminal half of the 2nd step factor Slu7. The endonuclease-like domain of Prp8, to which this domain of Slu7 binds, is correspondingly shifted. Interestingly, this state retains strong density for the 3’ SS, suggesting the N-terminal regions of Slu7 are not involved in maintaining 3’ SS docking after exon ligation. Finally, State III, containing 29% of particles, entirely lacks density for the 2nd step factors Slu7, Prp18, and Fyv6. Prp17 and the RNase H domain of Prp8 are still present but are weaker than in States I and II, suggesting they are more flexible in State III. Most of Prp22 can still be observed with weaker density except for the C-terminal domain that interacts with Prp8, which lacks density entirely. Importantly, State III also lacks density for the 3’ SS but has density for the 3’ exon, suggesting it is indeed post-catalytic and represents a state after loss of 2nd step factors and undocking of the 3’ SS from the catalytic core. The coincidence of loss of Fyv6 with loss of Prp18 and Slu7 reinforces the notion that these three factors act together and are important for 3’ SS docking and stabilization of Prp17, Prp22, and the RNase H domain of Prp8.
The Fyv6/Syf1 Interaction is Critical for Suppressing fyv61 Phenotypes
To determine the functionally critical regions and interactions of Fyv6, we used the cryo-EM structure to design five truncation mutants (Fig. 5A). The first three truncations remove amino acids from the N-terminus of Fyv6: deletion of the first 16 amino acids, which are not resolved in the cryo-EM structure (Δ1-16); deletion of the first 23 amino acids, which removes the conserved hook region of Fyv6 that interacts with the Prp22 RecA2 domain (Δ1-23); and deletion of amino acids 1 to 51, which additionally removes the first α-helix and interactions with Prp22, Prp8, and Slu7 (Δ1-51). The remaining two mutants are C-terminal truncations. In the first, amino acids 134 to 173, which do not have clear density in the cryo-EM structure and are not modeled but which do contain a predicted nuclear localization signal (NLS), are deleted (Δ134-173) (Nguyen Ba et al., 2009). In the second, amino acids 103 to 173 are deleted (Δ103-173), which additionally removes the region of the protein that interacts with Syf1.
![](https://prod--epp.elifesciences.org/iiif/2/100449%2Fv1%2Fcontent%2F592262v2_fig5.tif/full/max/0/default.jpg)
Structure-based analysis of Fyv6 domain function.
A) Diagram of Fyv6 protein structure, protein interactors, and protein truncations. NLS= nuclear localization signal B) Spot dilution assay of strains with Fyv6 truncations on –Trp DO plates at different temperatures. Plates were imaged after the number of days indicated. C) Representative gel image of SUS1 RT-PCR in strains with an empty vector or expressing the indicated Fyv6 variant. (+RT reactions contain reverse transcriptase; -RT control reactions do not contain reverse transcriptase.) D) Diagram of Syf1 protein structure and truncations. E) Spot dilution assay of strains with Syf1 truncations on –Trp DO plates at different temperatures. Plates were imaged after the number of days indicated. F) Representative gel image of SUS1 RT-PCR in strains with an empty vector or expressing the indicated Syf1 variant. G) Temperature growth assay of strains with Syf1 truncations in a fyv6Δ background on YPD plates. Plates were imaged after the number of days indicated. H) Representative gel image of RT-PCR of SUS1 in strains with plasmids containing WT Syf1 or a Syf1 truncation in a fyv6Δ background.
We tested each truncation mutant for its ability to suppress the cs and ts phenotypes of the fyv6Δ strain when expressed from a plasmid under the native promoter (Fig. 5B). All of the N-terminal truncation mutations could suppress the cs or ts phenotypes to some degree. Surprisingly, a strain with the complete deletion of the conserved hook domain that interacts with Prp22 (Δ1-23) grew similarly to the strain expressing WT Fyv6. We wondered if splicing changes could still be present even if the growth phenotype was suppressed; so, we assayed SUS1 isoforms by RT-PCR. We previously showed that loss of Fyv6 causes use of an alternative 3’ SS in the first intron of this pre-mRNA (Lipinski et al., 2023). In addition to suppression of the fyv6Δ growth phenotype, the N-terminal truncation mutants also showed no evidence for changes in SUS1 splicing (Fig. 5C). The N-terminal region of Fyv6 composed of the hook domain and the first α-helix is not required for suppression of either the temperature sensitivity or splicing phenotypes observed when Fyv6 is lost despite its conservation and interactions with critical splicing factors (Prp8, Prp22, and Slu7).
For the C-terminal truncations, the Δ134-173 truncation was also able to suppress the cs and ts phenotypes, indicating that the predicted NLS is also not essential (Fig. 5B). However, this mutant grew more poorly at 16°C than did the Δ1-16 and Δ1-23 truncations and showed some evidence of alternative 3’ SS usage in SUS1 (Fig. 5C). The largest effects were seen when the C-terminus was truncated further. The Δ103-173 truncation strain phenocopied the fyv6Δ strain and grew poorly at 37°C and was dead at 23° or 16°C (Fig. 5B). Additionally, RT-PCR showed increased use of the alternative 3’ SS in SUS1 with this truncation (Fig. 5C). To determine if this phenotype could be due to loss of protein expression, we assayed expression of the epitope-tagged proteins by western blot (Supplemental Fig. S9A). In all cases, we could detect protein; however, the abundance of the N-terminally epitope-tagged Δ103-173 variant was much lower than the others. Therefore, we also constructed a C-terminally epitope-tagged version of Δ103-173 (Supplemental Fig. S9B). The C-terminally-tagged protein expressed at much higher levels but also failed to suppress the cs and ts phenotypes. These data suggest that the C-terminal region of Fyv6 that includes a portion of α-helix 3 and the Syf1 interaction domain is critical for function.
To obtain additional evidence for the importance of the Fyv6/Syf1 interaction, we created a syf1Δ shuffle strain containing WT Fyv6 and expressed C-terminal Syf1 truncations (Fig. 5D). We hypothesized that loss of the Fyv6-interacting domain on Syf1 would phenocopy fyv6Δ. We observed just this effect. Loss of the Fyv6-interacting region in Syf1 Δ778-859 and Δ634-859 mutants resulted in either modest (Δ778-859) or severe (Δ634-859) cs and ts phenotypes. In addition, these Syf1 truncation mutants also triggered use of the alternative 3’ SS in SUS1 (Fig. 5F). Interestingly, if the Δ778-859 and Δ817-859 Syf1 mutants are expressed in a fyv6Δ background, it results in suppression of the fyv6Δ ts phenotype (Fig. 5G; c.f. Fig. 5B). However, the SUS1 alternative 3’ SS is still used (Fig. 5H). This indicates that ts phenotype suppression and triggering usage of the alternative 3’ SS in SUS1 are separable processes.
As mentioned above, the Syf1/Fyv6 interaction is mutually exclusive with the Syf1/Yju2 interaction that occurs during 5’ SS cleavage and is reformed in the ILS spliceosome. One possible explanation for suppression of the ts phenotype in fyv6Δ/syf1(Δ778-859 or Δ817-859) strains is that Fyv6 and Yju2 are competing with one another and binding of either protein to the wrong spliceosome complex results in growth and splicing defects. When the Yju2/Fyv6 binding site on Syf1 is removed, the observed phenotypes could be due to failure of Yju2 to be properly recruited to the P complex spliceosome for displacement of Fyv6. However, recruitment occurs correctly when Fyv6 is deleted since it no longer needs to be displaced. Similarly, phenotypes due to loss of Fyv6 could arise from retention of Yju2 during exon ligation. A prediction of this model is that failure to properly remove or recruit Yju2 and Fyv6 can result in similar phenotypes. Further work will be needed to understand the binding dynamics of these factors.
Mutations in Many Different Splicing Factors Can Suppress fyv6Δ Phenotypes
To gain additional insights into Fyv6 function, we carried out a genetic screen to detect spontaneously arising suppressors that can correct temperature-dependent growth defects (Fig. 6A). This method has previously been used to illuminate many aspects of spliceosome biochemistry and non-silent mutations within splicing factors are likely causative for suppression (Montemayor et al. 2014; Kuhn and Brow, 2000; Brow, 2019).
![](https://prod--epp.elifesciences.org/iiif/2/100449%2Fv1%2Fcontent%2F592262v2_fig6.tif/full/max/0/default.jpg)
Identification of novel suppressors of fyv61 temperature sensitivity.
A) Workflow for the suppressor screen. B) P complex spliceosome structure with fyv61 suppressor mutations labelled. (See Supplemental Table S5) C) Close-up view of Cef1 and Prp8 suppressor mutations in P complex. D) Spot dilution temperature growth assay with Cef1 mutants on YPD plates. Plates were imaged after the number of days indicated. E) Close-up view of Prp8 mutations in P complex. F) Spot dilution temperature growth assay with Prp8 mutants on YPD plates. Plates were imaged after the number of days indicated.
We plated saturated fyv6Δ yeast cultures onto 20 different YPD plates and incubated them at 37°C until colonies arose. Within several days, colonies of various sizes appeared and were isolated. A subset of picked colonies were then tested for suppression of the ts phenotype. All of these colonies were able to suppress growth defects at 37°C, and many also suppressed the cs phenotype at 23°C (Supplemental Fig. S10). We then performed whole genome sequencing on the largest of the picked colonies as these were expected to be the strongest suppressors. In 19 of 20 sequenced strains, we were able to identify at least one non-silent mutation in a splicing factor (Fig. 6B, Supplemental Table S5). When we carried out a similar experiment aimed at identifying cs suppressors, we were unable to identify any likely causative mutations in splicing factors. Therefore, we focused our analysis on the ts suppressor mutations.
Collectively, the identified mutations tended to fall into two categories: gain of positive charge or substitution with a bulkier side chain. Gain of positive charge mutations were found in Prp8 (G1868R) as well as the second step factors Slu7 (A24R, and N28K), and Prp22 (I1133R). Structurally, amino acids within Prp8 and Prp22 are proximal to RNA within the spliceosome and addition of charged amino acid side chains may stabilize interactions with the phosphodiester RNA backbone. The N-terminus of Slu7 containing A24 and N28 is not resolved in any spliceosome structure to our knowledge and normally contains many negatively charged residues. The majority of identified suppressors involve a mutation to a bulkier side chain including several in Prp8 (S1594F/Y, V1862L), Prp17 (A251V, G386W), Cef1 (A37V), and Clf1 (T33I). The increase in hydrophobic surface area provided by these mutations may help to stabilize interactions with nearby proteins, compensating for structural destabilization due to the absence of Fyv6.
We isolated one suppressor mutation within the U2 snRNA (A25G) near where the C-terminal domain of Prp22 inserts towards the active site and located beneath the Fyv6 binding site. It has previously been observed that mutation of this specific U2 snRNA nucleotide is not deleterious to growth (McPheeters and Abelson, 1992), and is therefore well-tolerated by yeast. The mechanism of suppression by this mutation is not clear; however, it is interesting to note that we also observed suppressor mutations in the U2 snRNP component Rse1 (Supplemental Table S5) which is released from the spliceosome prior to catalysis. The Rse1 mutations always occurred coincident with mutations in other splicing factors. It is possible that suppressors in factors present prior to or after exon ligation are able to suppress (or enhance suppression of) fyv6Δ by changing how pre-mRNAs compete for a limited set of splicing factors or by reducing blocks to pre-spliceosome assembly caused by accumulation of stalled complexes (Mendoza-Ochoa et al., 2019; Munding et al., 2013). Together these data indicate that mutations in many different splicing factors including second-step factors and the U2 snRNA can be isolated from fyv6Δ suppressor strains.
Multiple Second-step Splicing Defects can be Rescued by the Same Suppressor
We decided to study several of these mutations in greater detail and to confirm that they were sufficient for fyv6Δ ts phenotype suppression. In the case of Cef1, we obtained several mutants located near the interface between Cef1 and the 2nd step factor Prp17: A37V, M175I, and Q193P (Fig. 6C). This interface has only been observed in structures of C* and P-complex spliceosomes, suggesting its importance for exon ligation. In addition, it has previously been reported that another mutant of Cef1 A37, A37P, is able to suppress splicing defects due to a BS adenine to guanosine substitution (BS-G) which causes stalling after 5’ SS cleavage (Query and Konarska, 2012). Therefore, we wondered if other Cef1 mutants capable of suppressing defects due to BS-G could also suppress fyv6Δ (Cef1 A37P, V36R, and S48R). We created fyv6Δ shuffle strains for Cef1 and tested each mutant’s ability to suppress cs and ts phenotypes (Fig. 6D). All of these Cef1 mutants, both those isolated from our suppressor screen and those identified as BS-G suppressors, could at least weakly suppress the cs phenotype at 23°C with our isolated M175I mutant being a particularly strong suppressor. In addition, all of the mutants except Cef1 A37V could suppress the ts phenotype at 37°C to varying degrees. These results confirm that single point mutants in Cef1 can suppress fyv6Δ growth phenotypes and that at least some BS-G suppressors in Cef1 also suppress fyv6Δ. A common mechanism could be that these mutations stabilize the exon ligation conformation of the spliceosome to favor the 2nd step either when a BS-G is present or Fyv6 is absent.
In the case of Prp8, we tested five different suppressor mutations in prp8 shuffle strains (Fig 6E,F). Identification of mutants in Prp8 that suppress fyv6Δ is not surprising since many different alleles of Prp8 are known to modulate the conformational equilibrium between the 1st and 2nd steps (Fica and Nagai, 2017; Liu et al., 2007; Query and Konarska, 2004). We previously discovered that one Prp8 2nd step allele, P986T, could suppress a fyv6Δ growth defect (Lipinski et al., 2023).
We identified two suppressors within the Prp8 α-finger domain at position S1584 and another three within the Rnase H-like domain (Fig. 6E, Supplemental Table S5) (Galej et al. 2013). S1584 is located adjacent to a residue that crosslinks with the BP (aa 1585) (Galej et al., 2013). Proximity to the BP suggests that mutations at S1584 might destabilize RNA bound within the active site after branching, alleviating the stall caused by fyv6Δ after the 1st step or facilitating formation of the 2nd step conformation.
Suppressors of other splicing defects have been previously identified within the Prp8 RNase H-like domain, two of which are located at the same position as fyv6Δ suppressors. V1862A/Y/D and T1982A were identified as suppressors of, respectively, the U4-cs1 allele (Kuhn and Brow, 2000; Kuhn et al., 1999) or the 5’ SS U2A mutant (Siatecka et al., 1999). V1862 is located within the RNase H-like domain β-finger (Fig. 6E), while T1982 is located at the base of the β-finger (Fig. 6C). We identified fyv6Δ suppressors at both V1862 (V1862L) and T1982 (T1982S) as well as a third suppressor (G1868R), also located in the β-finger. Several 2nd step Prp8 alleles are also located proximal to G1868 (N1869D and V1870N) (Query and Konarska, 2004).
Four of the identified Prp8 suppressor mutations could partially suppress either or both of the cs and ts fyv6Δ phenotypes (Fig. 6F). Prp8 S1574Y, V1862L, G1868R, and T1982S suppressors were able to support growth at 37°C. Prp8 V1862L and T1982S also partially suppressed the growth defects at 16 and 23°C. Even though we isolated the Prp8 S1584F mutant in our screen and this substitution is very similar to S1584Y, we were not able observe suppression due to this mutant in these strains. Mutations within the Prp8 β-finger are able to compensate for loss of Fyv6 and these mutations are located near known Prp8 2nd step alleles. As with Cef1 mutants, it is possible that fyv6Δ suppression is originating from stabilization of the 2nd step conformation.
Finally, given our observations with the Cef1 and Prp8 mutants, we wondered about the generality of 2nd step splicing factor suppression. Several dominant mutants in the 2nd step factor Prp18 were previously shown to suppress ts phenotypes due to various mutations in Slu7 (Aronova et al., 2007). We wondered if these Prp18 mutants would also suppress the ts (or cs) phenotype due to Fyv6 loss. Indeed, the Prp18 V191A suppressor of slu7 alleles was able to improve growth at 37°C of fyv6Δ yeast (Supplemental Fig. S11). Together our results from the suppressor analysis support a model in which a given 2nd step suppressor mutation may be able to restore splicing activity lost by a variety of mechanisms.
A Prp22-Dependent Splicing Stall is Relieved by Fyv6 Deletion
The final fyv6Δ suppressor mutation that we studied was the Prp22 I1133R mutant. This mutation is located within the C-terminal domain of Prp22 that enters the spliceosome active site and is located proximal to U2 snRNA nucleotide A31 (Fig. 7A). Fyv6 bridges over this Prp22-U2 snRNA interaction. The I1133R mutation might stabilize Prp22/U2 snRNA contacts that can compensate for Fyv6 loss and promote exon ligation. The Prp22 I1133R mutant is able to weakly suppress both the cs and ts phenotypes in fyv6Δ yeast (Fig. 7B, C). Given its location in the C-terminal domain, we wondered if it or fyv6Δ would show synthetic genetic interactions or epistasis with mutations in the helicase domain. We tested two helicase domain mutations, R805A and G810A, both of which are known to impact Prp22-dependent proofreading (Mayas et al., 2006). The Prp22 R805A mutant is by itself cs (Schwer and Meszaros, 2000) and showed a synthetic lethal interaction with fyv6Δ (Fig. 7B, 7C). This was only slightly rescued by simultaneously introducing I1133R into Prp22. In contrast, Prp22 G810A mutant is viable in the presence and absence of Fyv6 as well as in combination with the I1133R suppressor mutation (Fig. 7B, 7C). Prp22 G810A by itself does not suppress the ts fyv6Δ phenotype; however, the Prp22 G810A/I1133R double mutant is a stronger suppressor of the ts phenotype than either mutant alone. These additive genetic interactions suggest that the Prp22 fyv6Δ suppressor mutation (I1133R) and the helicase domain mutations (R805A and G810A) impact distinct Prp22 functions. These functions could correspond to Prp22’s role in promoting exon ligation during the 2nd step (the C-terminal domain mutation) and its role in promoting ATP-dependent mRNA release and proofreading (the helicase domain mutations).
![](https://prod--epp.elifesciences.org/iiif/2/100449%2Fv1%2Fcontent%2F592262v2_fig7.tif/full/max/0/default.jpg)
Genetic interactions between Prp22 and Fyv6.
A) Location of the Prp22 I1133R fyv6Δ suppressor in the P complex structure. B) Spot dilution growth assay with Prp22 mutants in WT or fyv6Δ backgrounds grown on –Trp DO or –Trp+5-FOA plates. Plates were imaged after 3 days at 30°C. C) Spot dilution growth assay with strains containing Prp22 mutants on YPD plates at different temperatures. Plates were imaged after the number of days indicated. D) Representative gel of SUS1 RT-PCR products in strains with Prp22 mutants in the presence or absence of Fyv6. E) Quantitation of band intensities for each SUS1 splice isoform as a fraction of total product in the lane. Bars indicate average ± SD of N=3 replicates. F) ACT1-CUP1 reporter construct with the location of the 3’ SS substitution noted. G) Copper tolerance for each ACT1-CUP1 reporter with each Prp22 allele in strains with (▪) or without (○) Fyv6 present. Shown is a representative, single replicate of N=3. H) Model for the role of Fyv6 during the 2nd step. Fyv6 interacts with Syf1 and Prp22 and promotes usage of BP distal 3’ SS. When Fyv6 is absent, BP proximal 3’ SS are favored. It is possible that lack of Fyv6 also results in the presence of Yju2 during exon ligation and relaxing of the requirement for Prp22 at this step.
Recently, it was reported that the Prp22 G810A mutant could activate usage of some alternative 3’ SS that are also activated by PRP18 deletion (Roy et al., 2023). We therefore wondered if any of the Prp22 mutants could activate alternative 3’ SS also activated by FYV6 deletion. We isolated RNA from these strains and analyzed SUS1 splicing by RT-PCR (Fig. 7D, E). None of the Prp22 mutants activated usage of the SUS1 alternative 3’ SS or enhanced usage of this site beyond what is observed upon Fyv6 deletion. However, we noticed that the Prp22 helicase mutants, R805A and G810A, both caused an accumulation of unspliced pre-mRNA. This is consistent with these mutations causing a block during pre-mRNA splicing and having a different effect on splicing than Prp22 I1133R. Surprisingly, the block due to the G810A mutation was relieved in the fyv6Δ strains and little pre-mRNA accumulation was observed in the absence of Fyv6. This result suggests that Fyv6 could act as a negative regulator of splicing of the SUS1 transcript when a Prp22 helicase domain mutant is present. Removing Fyv6 restores splicing by activation of the SUS1 BP proximal 3’ SS.
To further explore the connections between Fyv6 and Prp22 fidelity, we used the ACT1-CUP1 assay in which yeast growth in the presence of increasing concentrations of Cu2+ is proportional to the extent of splicing of a reporter RNA (Fig. 7F, G; Supplemental Fig. S12, S13). Since the Prp22 R805A mutant grew poorly, we only assayed the Prp22 G810A, I1133R, and G810A/I1133R (double) mutants in the presence and absence of Fyv6. Consistent with previous results, we found that fyv6Δ is deleterious for splicing of the WT ACT1-CUP1 reporter containing consensus splice sites. Copper tolerance was improved when fyv6Δ was combined with Prp22 I1133R. Unexpectedly, Prp22 I1133R was actually deleterious when Fyv6 was present for this reporter. If this Prp22 mutant and Fyv6 both promote the 2nd step, then their combination may have an additive, negative effect on splicing due to over-stabilization of this spliceosome conformation.
Finally, we assayed strains containing a 3’ SS UAG to gAG reporter since it was previously shown that Prp22 helicase domain mutants could increase usage of nonconsensus 3’ SS due to faulty proofreading activity (Fig. 7F, G). In agreement with this, we observed a dramatic increase in copper tolerance when Prp22 contained both the C-terminal domain and helicase mutants (double). This suggests that splicing of nonconsensus 3’ SS splice sites is enhanced both by stabilization of the 2nd step conformation (I1133R) and perturbing mRNA release (G810A). This effect on Prp22 is Fyv6-dependent since loss of Fyv6 results in very poor copper tolerance of strains containing the gAG reporter, likely since this 3’ SS is located far from the BP (38 nt). These results show that fyv6Δ can be suppressed with a Prp22 C-terminal domain mutant and that mutations in the Prp22 C-terminal domain can have distinct functional consequences from those located in the helicase domain. In the case of SUS1, stalling of splicing due to the Prp22 G810A mutant can be relieved by removing Fyv6 and activation of a BP proximal 3’ SS.
Discussion
Here, we establish Fyv6 as a 2nd step splicing factor in yeast that promotes usage of BP-distal 3’ SS. Structural, genetic, and biochemical studies reveal critical contacts and interactions between Fyv6 and the splicing machinery. From these data, we propose that Yju2/Syf1 interactions need to be disrupted for proper Fyv6 recruitment and that Fyv6 can enforce the Prp22-dependence of at least some splicing events. Combined, our data allow us to propose a model for Fyv6’s role during the 2nd step (Fig. 7H). In this model, Fyv6 is recruited to the C* complex after Prp16-dependent release of Yju2, which exposes the binding site for Fyv6 on Syf1. Once bound, Fyv6 (possibly in concert with Slu7 and Prp18) facilitates docking of BP distal 3’ SS into the active site. In the absence of Fyv6, the use of BP distal 3’ SS becomes more difficult and BP proximal 3’ SS are favored. Under these conditions, the requirement for Prp22 during exon ligation is reduced and Yju2 may also potentially rebind Syf1. One or both of these factors could contribute to formation of a lower fidelity spliceosome with relaxed sequence constraints for exon ligation.
Structural Insights into the Recruitment and Release of 2nd step Factors
A key principle of spliceosome structural biochemistry is the presence of mutually exclusive interactions that coordinate the stepwise and directional progression of splicing (Townsend et al., 2020). Our cryo-EM structure of the P complex spliceosome suggests another such interaction is involved in the 1st to 2nd-step transition during which Fyv6 displaces Yju2. Yju2 and Fyv6 occupy overlapping binding sites on the splicing factor Syf1 and, in fact, unassigned density in spliceosome P complexes due to Fyv6 was previously mis-interpreted as Yju2 due to their structural similarity (Wilkinson et al., 2017; Wilkinson et al., 2021). We propose that Prp16-promoted spliceosome remodeling results in release of Yju2 from the spliceosome C complex and permits mutually exclusive binding of Fyv6. Fyv6 remains bound in the P complex but then dissociates to permit re-binding of Yju2 in the ILS complex.
Previously, Chiang and Cheng noted that while the N-terminal half of Yju2 was able to restore 1st step activity to Yju2-depleted extracts, it is much more stably associated if a second polypeptide containing the C-terminal, Syf1-interacting half was included (Chiang and Cheng, 2013). It is possible that competition between Yju2 and Fyv6 for Syf1 binding can explain these results. If only the N-terminal half of Yju2 is present, Fyv6 could associate prematurely with Syf1 and result in destabilization of the 1st step conformation in favor of 2nd step interactions. This may then result in weaker binding for the N-terminal half of Yju2. Failure to properly recruit and exchange the Yju2/Fyv6 interaction on Syf1 may also be the origin of some of the yeast growth and splicing phenotypes we observe (Fig. 5).
The need for Yju2/Fyv6 exchange hints at two distinct functions for Fyv6: a role in Yju2 displacement to prevent its re-binding during exon ligation and a role in stabilization of the 2nd step spliceosome conformation. These two functions can potentially result in the distinct cs and ts phenotypes in yeast. While more work is needed to dissect these phenotypes, it is possible that the cs phenotype results from increased stability of the Yju2/Syf1 interaction at lower temperatures and the need for Fyv6 to act as a competitor of this interaction. The ts phenotype on the other hand may result from destabilization of the 2nd step conformation of the spliceosome active site due to loss of Fyv6. It is possible that loss of Fyv6 results in destabilization due to removal of a buttress that holds the RNase H-like domain of Prp8 in place (along with Prp18), changes in how the C-terminal domain of Prp22 interacts with the spliceosome active site, removal of an impediment to reversal of the U2 snRNA/BS duplex conformational change that occurs during the C to C* complex transition, or a combination of these factors.
The high quality of our cryo-EM data combined with 3D classification also provide new insights into the requirements for 3’ SS docking into the active site and potentially for how the P complex spliceosome is disassembled into the ILS complex. It was previously speculated that the N-terminal domain of Slu7 may be critical for 3’ SS docking by reducing the entropic cost of this process when the BP to 3’ SS distance is long (Wilkinson et al., 2017). In our state II P complex structure, the 3’ SS remains docked but this domain of Slu7 is not observed. This suggests that either Slu7 N-terminal domain interactions with Cwc22 and Prp8 are not essential for 3’ SS docking into the active site for exon ligation and/or that these interactions are not essential for maintaining 3’ SS docking within P complex once they form in C*.
In the state III P complex structure, densities for the 2nd step factors Fyv6, Prp18, and Slu7 are entirely missing and the 3’ SS is no longer docked into the active site. This result suggests that 3’ SS undocking is coupled to 2nd step factor release. Prp22 is still present in this structure, and it is possible that state III represents a conformation occurring after exon ligation but prior to Prp22-dependent mRNA release. In this case, the state III structure implies that Prp22 ATPase activity is not needed for dissociation of Fyv6, Prp18, and Slu7. It is also interesting to note that the C-terminal domain of Prp22 can no longer be resolved in the state III structure. What exactly activates the Prp22 ATPase to promote proofreading or mRNA release is not known. It is conceivable that undocking of the 3’ SS from the active site can be sensed by the Prp22 C-terminal domain and this contributes to ATPase activation. This model differs from that previously proposed for Prp22-enabled 3’ SS sampling in which that ATPase activity of Prp22 is needed to undock the 3’ SS (Semlow et al., 2016). However, sampling could involve multiple steps: ATP-independent 3’ SS undocking and ATP-dependent entry into the sampling state. Whether or not the 2nd step factors must release or alter conformation to enable sampling is unknown. Regardless of the exact mechanism, it is important to recognize that the P complex structures determined here were all obtained using a Prp22 mutant defective in mRNA release. The chemical and kinetic competencies of these states will require future investigations.
Multiple, Non-Redundant 2nd Step Factors Promote Use of BP Distal 3’ SS
The challenge of 3’ SS selection involves not just selection of the correct nucleotide sequence (YAG) but also identifying which YAG sequence to use. Here, we define a function for Fyv6 in facilitating use of BP distal 3’ SS, especially those located > 21 nt away from the branch. This function is similar to those of other 2nd step factors important for efficient splicing when the BP-3’ SS distance is >12 nt (Slu7, Prp18) or >21 nt (Prp22) (Brys and Schwer, 1996; Schwer and Gross, 1998; Zhang and Schwer, 1997). Human Slu7 has additionally been reported to be important for preventing suppression of the correct 3’ SS in favor of up- or downstream alternative 3’ SS (Chua and Reed, 1999). Our previous study using competing 3’ SS in the ACT1-CUP1 intron did not provide evidence for suppression of the correct 3’ SS in the absence of Fyv6, suggesting that this function is not shared between human Slu7 and yeast Fyv6. Nonetheless, it is important to recognize that the splicing machinery employs multiple factors that enforce usage of BP distal 3’ SS. These functions are not redundant as changes in 3’ SS are still observed in the absence of Fyv6 but presence of Prp18 and vice versa (Roy et al., 2023).
Intron architecture can also contribute to 3’ SS selection since canonical “AG” 3’ SS are depleted in humans in the region between the BP and annotated 3’ SS (the AG exclusion zone, or AGEZ) (Corvelo et al., 2010; Gooding et al., 2006). AG sites are also typically depleted in the corresponding region in yeast introns (data not shown). This would suggest that annotated 3’ SS are also the most optimal sites located downstream of the BP. However, an important assumption of this model is that non-AG 3’ SS can also be effectively proofread and not used. Deletion of either Prp18 or Fyv6 leads to activation of not just BP proximal 3’ SS but also those with highly divergent sequences from the YAG consensus (Roy et al., 2023). This suggests that an AGEZ alone is insufficient for correct 3’ SS usage and that mechanisms to ensure preferential use of YAG sites are critical.
It is likely that Fyv6 and Prp18 enforce usage of YAG sites by common as well as factor-specific mechanisms. Both proteins stabilize the conformation of the Prp8 RNase H-like domain in P complex. Numerous splice site suppressor mutations map to this domain including those that enhance splicing of non-YAG 3’ SS (Galej et al., 2013). Fyv6 and Prp18 may be essential for establishing a “high fidelity” conformation of the RNase H-like domain that enforces use of canonical YAG 3’ SS and preventing activation of cryptic sites within the AGEZ. In support of this, we note that in our state III P complex structure that lacks both Fyv6 and Prp18 and a stably docked 3’ SS, the RNase H-like domain is also more flexible.
In addition to this possibility, Fyv6 and Prp18 may also contribute uniquely to 3’ SS fidelity and selection. The conserved loop region of Prp18 inserts into the spliceosome active site and contacts the 3’ SS. Partial deletion of this region of Prp18 phenocopies at least some of the changes in 3’ SS selection when Prp18 is deleted entirely (Crotti et al., 2007; Roy et al., 2023). This indicates that even when the Prp8 RNaseH-like domain is still stabilized by Prp18 (and Fyv6) disruptions to 3’ SS selection can still occur. Unlike Prp18, Fyv6 does not contact the spliceosome active site. However, it is the only 2nd step factor that contacts exon ligation fidelity factor Prp22. While these contacts are not essential for suppressing use of an alternative 3’ in SUS1, a mutation in Prp22 can nonetheless suppress a fyv61 phenotype and loss of Fyv6 relieves a Prp22 helicase mutant-dependent stall in SUS1 splicing. These observations suggest that Fyv6 and Prp22 are functionally connected to one another.
Some 2nd Step Factors May Enforce Usage of Prp22 during Exon Ligation
We were particularly struck by the similarities in the BP-3’ SS distances for the requirement of Prp22 during exon ligation and for promotion of BP distal 3’ SS usage by Fyv6, both >21 nt (Schwer and Gross, 1998). In addition, the alternative, BP proximal 3’ SS that are activated when Fyv6 is deleted show much less sequence conservation relative to the annotated 3’ SS and include sequences (GAG, BG) known to be proofread and rejected by Prp22. Together, this suggests that one function of Fyv6 is to enforce Prp22-dependence of exon ligation by stabilizing a spliceosome conformation that more easily enables usage of BP distal 3’ SS (Fig. 7H). As these sites require Prp22 for exon ligation, it is also likely that they are subject to Prp22-dependent proofreading, and this may enforce accurate exon ligation at the YAG 3’ SS consensus. When Fyv6 is absent, many alternative, BP proximal 3’ SS are activated. These sites could have a relaxed requirement for Prp22 and proofreading. While the Fyv6-dependence of Prp22 occupancy is beyond the scope of the current manuscript, loss of Prp18 (which also activates usage of non-canonical, BP-proximal 3’ SS) does result in lower occupancy of Prp22 on the ACT1 intron (Strittmatter et al., 2021). This suggests that Prp22 occupancy, and possible involvement in the exon ligation step, can be influenced by 2nd step factors.
There is precedent for nonessential splicing factors enforcing use of a proofreading ATPase. Cus2 is not essential for pre-spliceosome formation (Yan et al., 1998). However, when Cus2 is present then ATPase activity of Prp5 is essential for disruption of the Hsh155/Cus2 interaction (Perriman et al., 2003; Talkish et al., 2019). In this case, Cus2 is not essential for Prp5 proofreading of BP selection, and the connection between the ATP hydrolysis activity of Prp5, removal of Cus2, and proofreading is unclear (Liang and Cheng, 2015; Talkish et al., 2019; Xu and Query, 2007). It is possible that enforcement of Prp5 ATPase activity by Cus2 is important for other aspects of pre-spliceosome formation in vivo and the impact of Prp5 and Cus2 on BP selection has not yet been studied transcriptome-wide using high-resolution, sequencing-based approaches. The C-terminal domain of the 1st step splicing factor Yju2 (the same domain that competes with Fyv6 for Syf1 binding) is also not essential for 5’ SS cleavage. In the absence of this domain, the N-terminal domain of Yju2 is sufficient for promoting the 1st step reaction, and the spliceosome can transition to the 2nd step and carry out exon ligation without Prp16, albeit inefficiently (Chiang and Cheng, 2013). Importantly, addition of the C-terminal domain on a separate polypeptide reduces Prp16-independent activity. This suggests that the Yju2 C-terminal domain enforces the use of the Prp16 ATPase. As we propose with Fyv6 and Prp22, the C-terminal domain of Yju2 may not just help to promote the 1st to 2nd step transition by enforcing Prp16-dependence but could also facilitate Prp16-dependent proofreading during 5’ SS cleavage (Koodathingal et al., 2010; Semlow et al., 2016). General operating principles of the yeast splicing machinery may not just include proofreading but also use of splicing factors that enforce usage of the proofreading ATPases.
Materials and methods
The yeast strains and plasmids used in this study are described in Supplemental Tables S6 and S7. Yeast transformation, plasmid shuffling/5-FOA selection, and growth were carried out using standard techniques and media (Treco and Lundblad 1993; Sikorski and Boeke 1991).
Cryo-electron microscopy
Purification of dominant-negative Prp22
Yeast Prp22 with a dominant-negative S635A mutation was cloned into pRS424 and pRS426 vectors after an N-terminal CBP-His-TEV tag (Galej et al., 2013). The protein was expressed in 24 L of BCY123 yeast cells essentially as described (Galej et al., 2013). Recombinant Prp22 was then purified first by Calmodulin-Sepharose in a buffer containing 50 mM Tris-HCl pH 8.5, 350 mM NaCl, 2 mM CaCl2, 1 mM Mg acetate, 1 mM imidazole, 0.1% NP-40 substitute, 5 mM β-mercaptoethanol and eluted with 3 mM EGTA. The protein was dialyzed into 20 mM Tris-HCl pH 7.4, 350 mM NaCl, 10 mM imidazole, 0.1% NP-40 substitute and further purified over Ni-NTA resin at 500 mM NaCl and eluted with 200 mM imidazole. The eluted protein was treated with His-tagged TEV protease during dialysis against 20 mM HEPES-KOH pH 7.9, 300 mM KCl, 0.2 mM EDTA, 20% glycerol, 0.5 mM DTT, before removal of the protease with Ni-NTA resin. The final protein was frozen in liquid nitrogen and stored at -80 °C.
Preparation and purification of P complex for cryo-EM
P complex was purified as described in (Wilkinson et al., 2017). Brr2-TAPS yeast were grown in a 120 L fermenter and splicing extract was prepared by the liquid nitrogen method essentially as previously described (Lin et al., 1985). A pre-mRNA substrate consisting of 20 nt 5’-exon, 95 nt intron, and 32 nt 3’-exon from the UBC4 pre-mRNA followed by 3xMS2 stem loops was generated by run-off in vitro transcription from an EcoRI-linearized plasmid template. The RNA product was labelled at the 3’ end with fluorescein-5-thiosemicarbazide (Wu et al., 1996). A 180 mL in vitro splicing reaction was performed containing 60 mM potassium phosphate pH 7, 2 mM ATP, 2.5 mM MgCl2, 3% w/v PEG8000, 3 nM pre-mRNA, 37.5 nM MS2-MBP protein, and 40% (v/v) splicing extract treated with dominant negative Prp22 S635A. Reactions were incubated for 30 min at 23°C, then incubated for 20 min with 5 mM of a DNA oligonucleotide complementary to the 3’-exon (oligo sequence 5’-ATGAAGTAGGTGGAT-3’) to induce cleavage of the 3’-MS2 tag by the endogenous RNase H activity of the splicing extract. The reaction mixture was centrifuged through a 40% glycerol cushion in buffer A (20 mM HEPES, pH 7.9, 75 mM KCl, 0.25 mM EDTA). The cushion was collected and applied to amylose resin in the presence of 0.025% NP-40 substitute. After overnight incubation at 4°C the resin was washed and eluted with buffer A containing 5% glycerol, 0.01% NP-40 substitute and 12 mM maltose. Fractions containing spliceosomes were concentrated to 0.1 mg/mL then dialyzed against buffer A for 3 h.
For cryo-EM grid preparation, a freshly glow-discharged Cu300 R2/2 holey carbon grid with a 2-nm layer of amorphous carbon (Quantifoil) was mounted in the chamber of a Vitrobot Mark IV (Thermo Fisher Scientific) maintained at 4°C and 100% humidity. Spliceosomes (3 µL) were applied and after 60 s was manually blotted using Ø55 grade 595 filter paper (Ted Pella), followed by application of another 3 μL of spliceosomes, waiting 60 s, before blotting again and plunging into liquid ethane.
Cryo-EM data collection
Cryo-EM data were collected as two datasets from two grids made from the same sample. Dataset 1 was acquired using a Thermo Scientific Titan Krios cryo-TEM (LMB Krios 1) and dataset 2 was acquired using a Thermo Scientific Titan Krios G3i cryo-TEM (LMB Krios 3). Both used a K3 direct detector (Gatan) operated in super-resolution mode with 2-fold binning, and an energy filter with slit width of 20 eV. Micrographs were collected automatically using EPU in AFIS mode, yielding a total of 51,113 movies at 81,000x or 130,000x magnification with a real nominal pixel size of 0.93 Å or 0.669 Å, with defocus ranging from -1.3 μm to -3.1 μm and a total fluence per micrograph of 40 e-/Å2, fractionated into 40 frames.
Cryo-EM data processing
All cryo-EM data were processed using RELION-5.0 (Kimanius et al., 2023). Both datasets were initially processed separately. Movies were corrected for motion using the RELION implementation of MotionCor2, with 4×4 patches and dose-weighting. CTF parameters were estimated using CTFFIND-4.1. Particle picking was done using Topaz with the general model (Bepler et al., 2019), yielding 1,820,457 particles. 3D classification was performed on 4-fold (dataset 1) or 5-fold (dataset 2) binned particles using a P-complex reference map (EMD-10140), and 299,741 (dataset 1) or 599,123 (dataset 2) were selected and refined to 2.99 Å or 2.78 Å resolution respectively. After Bayesian Polishing and CTF refinement (per-particle defocus, anisotropic magnification, beam tilt, trefoil, and 4th order aberrations) the resolution improved to 2.67 Å or 2.32 Å respectively. To improve the general density quality, 3D classification without alignment was performed with a soft mask around the whole complex with T=15. Particles with good density, without selecting for Fyv6 occupancy, were selected (151,443 and 257,501 for datasets 1 and 2) and refined to 2.66 Å or 2.27 Å respectively. The particles were then merged, refined, and then refined for anisotropic magnification. The overall scaling factors from the “anisotropic” magnification matrix was then used to calculate the corrected pixel size for dataset 1 as 0.94057 Å, keeping dataset 2 constant at 0.669 Å per pixel. Defocus values for dataset 1 were then scaled by the square of 0.94057 / 0.93, with an empirical constant correction of -10.06 Å (Wilkinson et al., 2019). Finally, another round of per-particle defocus refinement and anisotropic magnification refinement was performed to reduce errors in defocus correction and dataset scaling. The combined particles refined to 2.24 Å resolution. Resolution is reported using the gold-standard Fourier Shell Correlation with 0.143 cutoff.
The refined particles after merging were subjected to 3D classification without alignment, with a soft mask around the entire spliceosome, with T=15. The resulting four classes were used to define States I, II and III. To improve the densities for the periphery of the spliceosome that displayed local flexibility, focused refinements were used with Blush regularization in RELION-5 to reduce overfitting (Kimanius et al., 2023). The Prp22 density was improved using a soft-mask just around the Prp22 helicase domain for refinement of particles from States I and II, yielding a 2.96 Å map. Densities for the U2 snRNP, NTC, U5 Sm ring, and Cwc22 N-terminal domain were each improved first by 3D classification without alignment of particles from all three states, selecting particles with strong density (72,869 for the U2 snRNP, 260,409 for the NTC, 113,561 for the U5 Sm ring, 201,984 for Cwc22-NTD), performing signal subtraction to remove density for the rest of the spliceosome (recentering on the mask center of mass and re-boxing to 300 pixels, or 192 pixels for Cwc22-NTD), and finally 3D refinement starting with 1.8 degree local angular searches, giving 3.42 Å for the U2 snRNP, 3.72 Å for the NTC, 3.06 Å for the U5 Sm, and 3.57 Å for Cwc22-NTD. In all cases, control refinements without Blush regularization did not produce interpretable maps (Supplemental Fig. S6).
Model building
The model was built using Coot (Casañal et al., 2020) and ISOLDE (Croll, 2018) starting with PDB 6EXN (Wilkinson et al., 2017) as an initial model for the core and PDB 7B9V (Wilkinson et al., 2021) as an initial model for the periphery. The core was improved by manual fixing in Coot, while the peripheral U2 snRNP and NTC were improved by applying torsion and distance restraints to reference models produced by AlphaFold2. The model was refined using PHENIX real_space_refine (Liebschner et al., 2019), just performing one macro-cycle. Figures were generated using UCSF ChimeraX.
Plasmid cloning
Fyv6 plasmids were made by PCR amplification from gDNA of the FYV6 coding sequence plus ∼250 nt upstream and downstream, restriction digest with NotI and SalI, and ligation into pRS414 or pRS416. An N-terminal 1x FLAG tag was added onto Fyv6 in the plasmids by PCR mutagenesis. Fyv6 truncations were cloned by Genewiz (Azenta Life Sciences) from pAAH1572. Restriction digest and ligation were used to place Fyv6 constructs into the pRS413 backbone.
Syf1 plasmids were made by PCR amplification from gDNA of the SYF1 coding sequence plus ∼275 nt upstream and downstream, restriction digest, and ligation into pRS414 or pRS416.
Cef1 mutant plasmids, Prp22 I1133R mutant plasmid, Prp8 mutant plasmids, ACT1-CUP1 BP-3’ SS distance reporters, and Syf1 truncation plasmids were prepared by PCR-based site directed mutagenesis of pAAH1611, pAAH1042, pAAH1440, pAAH0470, and pAAH1625 respectively. All plasmids were fully sequenced to confirm their nucleotide sequence.
Yeast strain creation
The DBR1 and FYV6 genes were deleted by replacement with a nourseothricin resistance cassette (NatMX4) or hygromycin resistance cassette (hphMX4), respectively, through homologous recombination (see Supplemental Table S7; Goldstein and McCusker, 1999). The SYF1 and CEF1 genes were deleted by replacement with a kanamycin resistance cassette (KanMX) through homologous recombination in strains transformed with pRS416-Syf1 or pRS316-Cef1 respectively. Gene deletion was confirmed by colony PCR and/or genomic DNA extraction from the strains and PCR amplification of the corresponding genomic locus.
RNA isolation
Yeast were grown in YPD media until OD600 0.5-0.8. For RNA-sequencing of temperature-shifted samples, cultures were shifted to either 37°C or 16°C or remained at 30°C and grown for 1 hr. A 15 mL volume of culture was harvested via vacuum filtration with mixed cellulose ester filter paper (Whatman, 47 mm diameter) in a Buchner funnel. Filter papers were separated into 5 mL tubes and flash frozen. Filters were stored at –80°C until RNA isolation. Filters were washed with 1 mL TRIzol to collect cells. For all other experiments involving RNA isolation, 10 OD600 units of OD600 0.5-0.8 yeast were harvested via centrifugation.
Cell walls were disrupted by incubation with Zymolyase (50 U/mL, Zymo Research) in Y1 buffer (1 M sorbitol, 100 mM EDTA pH 8.0, 13 mM β-mercaptoethanol) (primer extension, RNA-seq of non-temperature-shifted samples) or vortexing with silica disruption beads in TRIzol (RT-PCR, RNA-seq of temperature-shifted samples). Total RNA was then isolated using TRIzol reagent (Thermo Fisher Scientific) according to manufacturer’s instructions. Samples were cleaned using a 50 μg Monarch RNA Cleanup Kit (New England Biolabs) and treated with TURBO DNase (Thermo Fisher Scientific) to remove residual contaminating DNA. For RNA-sequencing, quality of RNA samples was assessed prior to library preparation by NanoDrop (concentration, A260/A280, and A260/A230), Qubit (RNA High Sensitivity and RNA IQ), and TapeStation (RNA ScreenTape Analysis). High quality samples showing minimal rRNA degradation with RIN >7 and A260/A280 ratios of ∼2.0 were selected for library preparation.
Library preparation and RNA sequencing
Library preparation was conducted by the University of Wisconsin-Madison Biotechnology Center Gene Expression Center. mRNA was selected from total RNA samples through polyA enrichment with the TruSeq Stranded mRNA kit (Illumina). Sequencing was performed with an NovaSeq6000 sequencer (Illumina).
Bioinformatic analysis of RNA-seq datasets
FASTQC was used for quality control of reads pre- and post-trimming (Andrews, 2010). Trimming of reads was accomplished using fastp (Chen et al., 2018) prior to mapping to the SacCer3 genome (Ensembl, R64-1-1) with STAR (Dobin et al., 2013). All samples were aligned with STAR in a first pass. Novel junctions from all samples (WT and fyv6Δ) were combined and filtered for likely false-positive junctions (see Supplemental Methods). STAR was run in a second pass with the additional input of the filtered set of novel junctions. Indexing of bam files was accomplished with SAMtools (Li et al., 2009). Quantitation of read counts aligned to the SacCer3 transcriptome was accomplished with Salmon (Patro et al., 2017). Differential gene expression analysis was conducted with DeSeq2 (Love et al., 2014). Junction reads within genes were counted with featureCounts (Liao et al., 2014) and used to calculate fraction of Annotation Splicing (FAnS) based on methods in Roy et al., 2023 (see Supplemental Methods). Canonical branch points were used as defined in the Ares Intron Database (Grate and Ares, 2002). Differential splicing analysis was conducted with SpliceWiz (Wong et al., 2024).
Primer extension
RNA was isolated using the protocol above. IR700 dye conjugated probes (Integrated DNA Technologies) were used for primer extension with SuperScript III (Thermo Fisher) of the ACT1-CUP1 reporters (2 pmol yAC6: /5IRD700/GGCACTCATGACCTTC) and U6 snRNA (0.4 pmol yU6: /5IRD700/GAACTGCTGATCATGTCTG) (Carrocci et al., 2017; van der Feltz et al., 2021). Primer extension products were visualized on a 7% (w/v) denaturing polyacrylamide gel run at 35 W for 80 min at RT. Gels were imaged with an Amersham Typhoon NIR laser scanner (Cytiva), and band intensities were quantified with Image J (version 1.53v, 2022).
RT-PCR
Yeast cultures were inoculated from cultures grown to stationary phase overnight and grown until OD600 = 0.7-0.9. RNA was isolated and depleted of contaminating DNA using the MasterPure Yeast RNA Purification Kit (LGC Biosearch Technologies) protocol as previously described (Carrocci et al., 2017) (RT-PCR shown in Fig. 4C) or using the RNA isolation protocol above (all others).
RT-PCR reactions were set up using the Access RT-PCR system (Promega Corporation) following the kit procedure with 100 ng of input RNA per 25 μl reaction. Primers to amplify SUS1 were SUS1-exon1 5’-TGGATACTGCGCAATTAAAGAGTC-3’ and SUS1-exon3 5’-TCATTGTGTATCTACAATCTCTTCAAG-3’. Reaction products were separated on 2% (w/v) agarose-TBE gels and imaged using ethidium bromide fluorescence.
Selection of fyv6Δ suppressors
Strain yAAH3353 was previously described (Lipinski et al. 2023). Twenty single colonies of yAAH3353 were selected from a YPD plate and used to inoculate 20 x 5mL YPD liquid culture. Cultures were grown for 48 h at 30°C. A 1 mL volume of culture was spun down, resuspended in 200 µL YPD liquid medium, and spread on a YPD plate containing 40 mg/L adenine hemisulfate (to prevent ade2 reversion). Duplicate plates for each overnight were placed in incubators at either 18°C or 37°C. Many colonies of various sizes appeared on each plate by 9 days for 37°C and 12 days for 18°C. Colonies were picked after that time and annotated first by growth temperature, then by overnight number, and finally roughly by colony size with 01 generally being the largest colony on a plate (ex: 370101 from a plate grown at 37°C inoculated from overnight 01 and the largest colony selected from the plate). Picked colonies were grown overnight in 5mL YPD at 30°C and frozen stocks were made. Anywhere from 0 to ∼15 colonies were picked per plate, covering the variety of colony sizes that arose.
Sequencing of yeast genomic DNA
Yeast were grown to saturated stationary phase in 5mL YPD and ∼10 OD units were sent as a frozen yeast pellet to Azenta or SeqCoast for genomic DNA extraction, library, preparation, and sequencing with 30x/400 mbp coverage.
Bioinformatic analysis of WGS data
Azenta and SeqCoast utilized the SacCer3 genome (Ensembl, R64-1-1) for variant calling and provided a list of >6,000 SNPs. Since the strain is not a perfect match for this reference, non-unique SNPs were filtered from each dataset by taking the union of two sets with each set obtained from colonies isolated on different plates (to reduce the chance of computationally eliminating the same variant SNP). Filtering reduced the dataset to <100 unique variants per colony. Few unique SNPs were within protein- or snRNA-coding genes. Of SNPs within protein-coding regions, single amino acid changes were frequently found within splicing factors (see Supplemental Table S4). Datasets were reexamined upon identification of suppressor mutations for suppressor mutation duplication that were potentially filtered out during comparison of SNP datasets.
Yeast temperature growth assays
Yeast were grown to stationary phase overnight in the appropriate liquid medium indicated in the figure legends. The cultures were first diluted to OD600 0.5, serially diluted 1:10, 1:100, and 1:1000, and stamped onto agar plates containing the appropriate growth medium. The plates were incubated at 16°, 23°, 30°, and 37°C, and plates were imaged after the number of days indicated in each figure.
ACT1-CUP1 copper tolerance assays
Yeast strains expressing ACT1-CUP1 reporters were grown to stationary phase in -Leu DO media, diluted to OD600 = 0.5 in 10% (v/v) glycerol, and spotted onto -Leu DO plates containing 0 to 2.5 mM CuSO4 (Lesser and Guthrie 1993; Carrocci et al., 2018). Plates were scored and imaged after 48 h of growth at 30°C for WT strains and after 72 h of growth at 30°C for fyv6Δ strains due to differential growth between strains.
Acknowledgements
We thank David Brow for his advice on the suppressor screen. We thank Justin Mabin for his advice on RNA-sequencing data analysis. MEW is grateful to Feng Zhang for funding and support, and to the labs of Kiyoshi Nagai and Kelly Nguyen for reagents and support. The authors utilized the University of Wisconsin – Madison Biotechnology Center Gene Expression Center (Research Resource Identifier – RRID:SCR_017757) for preparation of NGS libraries for RNA-sequencing and the Biochemistry Computational Research Facility (BCRF) for analysis of RNA-sequencing datasets. This work was supported by grants from the National Institutes of Health (R35 GM136261 to AAH), NIH Biotechnology Training Program (T32 GM135066) and National Science Foundation Graduate Research Fellowship Program (Grant No. DGE-1747503) fellowships to KAS, a Helen Hay Whitney Postdoctoral Fellowship (to MEW). This study was supported by the MRC Laboratory of Molecular Biology electron microscopy facility. We thank Guillaume Chanfreau for advice and Florian Heyd and Tucker Carrocci for critical reading of the manuscript.
Contributions
AAH, KAS, KAL, and MEW conceptualized the project; MEW collected and analyzed the cryo-EM data; KAL and NJZ prepared RNA for sequencing; KAL, NJZ, and AFG analyzed the RNA-sequencing data; AAH, KAS, NJZ, and KAL created yeast strains; KAS, KAL, NJZ, and AFG performed genetic and biochemical assays and analyzed data; KAS, KAL, MEW, and AAH wrote the manuscript.
Data and materials availability
RNA-seq data has been uploaded to NCBI SRA with BioProject ID PRJNA1113593. The cryo-EM maps will be deposited in the Electron Microscopy Data Bank and the atomic model to the Protein Data Bank upon manuscript submission. Raw cryo-EM data will be available from EMPIAR.
Competing interest statement
AAH is a member of the scientific advisory board and carrying out sponsored research for Remix Therapeutics.
References
- FastQC: a quality control tool for high throughput sequence data
- Functional interactions between Prp8, Prp18, Slu7, and U5 snRNA during the second step of pre-mRNA splicingRNA 13:1437–1444https://doi.org/10.1261/rna.572807
- Structure of the Post-catalytic Spliceosome from Saccharomyces cerevisiaeCell 171:1589–1598https://doi.org/10.1016/j.cell.2017.10.038
- Positive-unlabeled convolutional neural networks for particle picking in cryo-electron micrographsNature Methods 16https://doi.org/10.1038/s41592-019-0575-8
- An Allosteric Network for Spliceosome Activation Revealed by High-Throughput Suppressor Analysis in Saccharomyces cerevisiaeGenetics 212:111–124https://doi.org/10.1534/genetics.119.301922
- Requirement for SLU7 in yeast pre-mRNA splicing is dictated by the distance between the branchpoint and the 3’ splice siteRNA 2:707–717
- Functional analysis of Hsh155/SF3b1 interactions with the U2 snRNA/branch site duplexRNA 24:1028–1040https://doi.org/10.1261/rna.065664.118
- SF3b1 mutations associated with myelodysplastic syndromes alter the fidelity of branchsite selection in yeastNucleic Acids Research 45:4837–4852https://doi.org/10.1093/nar/gkw1349
- Current developments in Coot for macromolecular model building of Electron Cryo-microscopy and Crystallographic DataProtein Science 29:1055–1064https://doi.org/10.1002/pro.3791
- Identification and Characterization of Two Novel Components of The Prp19p-associated Complex, Ntc30p and Ntc20p *Journal of Biological Chemistry 276:488–494https://doi.org/10.1074/jbc.M006958200
- fastp: An ultra-fast all-in-one FASTQ preprocessorBioinformatics 34:i884–i890https://doi.org/10.1093/bioinformatics/bty560
- A Weak Spliceosome-Binding Domain of Yju2 Functions in the First Step and Bypasses Prp16 in the Second Step of SplicingMolecular and Cellular Biology 33:1746–1755https://doi.org/10.1128/MCB.00035-13
- The RNA splicing factor hSlu7 is required for correct 3′ splice-site choiceNature 402:207–210https://doi.org/10.1038/46086
- An ATP-independent role for Prp16 in promoting aberrant splicingNucleic Acids Research 51:10815–10828https://doi.org/10.1093/nar/gkad861
- Yeast ortholog of the Drosophila crooked neck protein promotes spliceosome assembly through stable U4/U6.U5 snRNP additionRNA 5:1042–1054
- Requirement of the RNA helicase-like protein PRP22 for release of messenger RNA from spliceosomesNature 349:487–493https://doi.org/10.1038/349487a0
- Genome-Wide Association between Branch Point Properties and Alternative SplicingPLOS Computational Biology 6https://doi.org/10.1371/journal.pcbi.1001016
- ISOLDE: A physically realistic environment for model building into low-resolution electron-density mapsActa Crystallographica Section D: Structural Biology 74:519–530https://doi.org/10.1107/S2059798318002425
- The Prp18 protein stabilizes the interaction of both exons with the U5 snRNA during the second step of pre-mRNA splicingGenes & Development 21:1204–1216https://doi.org/10.1101/gad.1538207
- STAR: Ultrafast universal RNA-seq alignerBioinformatics 29:15–21https://doi.org/10.1093/bioinformatics/bts635
- Structural and functional characterisation of human RNA helicase DHX8 provides insights into the mechanism of RNA-stimulated ADP releaseBiochemical Journal 476:2521–2543https://doi.org/10.1042/BCJ20190383
- Saccharomyces cerevisiae Ecm2 modulates the catalytic steps of pre-mRNA splicingRNA 27:591–603https://doi.org/10.1261/rna.077727.120
- Cryo-electron microscopy snapshots of the spliceosome: Structural insights into a dynamic ribonucleoprotein machineNature Structural & Molecular Biology 24https://doi.org/10.1038/nsmb.3463
- A human postcatalytic spliceosome structure reveals essential roles of metazoan factors for exon ligationScience 363:710–714https://doi.org/10.1126/science.aaw5569
- An essential splicing factor, SLU7, mediates 3’ splice site choice in yeastGenes & Development 6:2112–2124https://doi.org/10.1101/gad.6.11.2112
- Crystal structure of Prp8 reveals active site cavity of the spliceosomeNature 493:638–643https://doi.org/10.1038/nature11843
- Cwc21p promotes the second step conformation of the spliceosome and modulates 3′ splice site selectionNucleic Acids Research 43:3309–3317https://doi.org/10.1093/nar/gkv159
- Three new dominant drug resistance cassettes for gene disruption in Saccharomyces cerevisiaeYeast 15:1541–1553https://doi.org/10.1002/(SICI)1097-0061(199910)15:14<1541::AID-YEA476>3.0.CO;2-K
- A class of human exons with predicted distant branch points revealed by analysis of AG dinucleotide exclusion zonesGenome Biology 7https://doi.org/10.1186/gb-2006-7-1-r1
- Methods in Enzymology:380–392https://doi.org/10.1016/S0076-6879(02)50975-7
- Widespread Use of Non-productive Alternative Splice Sites in Saccharomyces cerevisiaePLOS Genetics 10https://doi.org/10.1371/journal.pgen.1004249
- Data-driven regularisation lowers the size barrier of cryo-EM structure determinationbioRxiv https://doi.org/10.1101/2023.10.23.563586
- The DEAH Box ATPases Prp16 and Prp43 Cooperate to Proofread 5′ Splice Site Cleavage during Pre-mRNA SplicingMolecular Cell 39:385–395https://doi.org/10.1016/j.molcel.2010.07.014
- Suppressors of a Cold-Sensitive Mutation in Yeast U4 RNA Define Five Domains in the Splicing Factor Prp8 That Influence Spliceosome ActivationGenetics 155:1667–1682https://doi.org/10.1093/genetics/155.4.1667
- Splicing Factor Prp8 Governs U4/U6 RNA Unwinding during Activation of the SpliceosomeMolecular Cell 3:65–75https://doi.org/10.1016/S1097-2765(00)80175-6
- Mutational analysis of pre-mRNA splicing in Saccharomyces cerevisiae using a sensitive new reporter gene, CUP1Genetics 133:851–863https://doi.org/10.1093/genetics/133.4.851
- The Sequence Alignment/Map format and SAMtoolsBioinformatics 25:2078–2079https://doi.org/10.1093/bioinformatics/btp352
- A novel mechanism for Prp5 function in prespliceosome formation and proofreading the branch site sequenceGenes & Development 29:81–93https://doi.org/10.1101/gad.253708.114
- featureCounts: An efficient general purpose program for assigning sequence reads to genomic featuresBioinformatics 30:923–930https://doi.org/10.1093/bioinformatics/btt656
- Macromolecular structure determination using X-rays, neutrons and electrons: Recent developments in PhenixActa Crystallographica. Section D, Structural Biology 75:861–877https://doi.org/10.1107/S2059798319011471
- Yeast mRNA splicing in vitroJournal of Biological Chemistry 260:14780–14792https://doi.org/10.1016/S0021-9258(17)38641-6
- Biochemical and Genetic Evidence Supports Fyv6 as a Second-Step Splicing Factor in Saccharomyces cerevisiaeRNA, rna 79607https://doi.org/10.1261/rna.079607.123
- Opposing classes of prp8 alleles modulate the transition between the catalytic steps of pre-mRNA splicingNature Structural & Molecular Biology 14https://doi.org/10.1038/nsmb1240
- Structure of the yeast spliceosomal postcatalytic P complexScience 358https://doi.org/10.1126/science.aar3462
- Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2Genome Biology 15https://doi.org/10.1186/s13059-014-0550-8
- Exon ligation is proofread by the DExD/H-box ATPase Prp22pNature Structural & Molecular Biology 13https://doi.org/10.1038/nsmb1093
- Mutational analysis of the yeast U2 snRNA suggests a structural similarity to the catalytic core of group I intronsCell 71:819–831https://doi.org/10.1016/0092-8674(92)90557-s
- Blocking late stages of splicing quickly limits pre-spliceosome assembly in vivoRNA Biology 16:1775–1784https://doi.org/10.1080/15476286.2019.1657788
- Core structure of the U6 small nuclear ribonucleoprotein at 1.7-Å resolutionNature Structural & Molecular Biology 21https://doi.org/10.1038/nsmb.2832
- Yeast vectors for the controlled expression of heterologous proteins in different genetic backgroundsGene 156:119–122https://doi.org/10.1016/0378-1119(95)00037-7
- Competition between Pre-mRNAs for the Splicing Machinery Drives Global Regulation of SplicingMolecular Cell 51:338–348https://doi.org/10.1016/j.molcel.2013.06.012
- NLStradamus: A simple Hidden Markov Model for nuclear localization signal predictionBMC Bioinformatics 10https://doi.org/10.1186/1471-2105-10-202
- Salmon provides fast and bias-aware quantification of transcript expressionNature Methods 14:417–419https://doi.org/10.1038/nmeth.4197
- ATP requirement for Prp5p function is determined by Cus2p and the structure of U2 small nuclear RNAProceedings of the National Academy of Sciences 100:13857–13862https://doi.org/10.1073/pnas.2036312100
- Structural Basis of Nuclear pre-mRNA Splicing: Lessons from YeastCold Spring Harbor Perspectives in Biology 11https://doi.org/10.1101/cshperspect.a032391
- Suppression of Multiple Substrate Mutations by Spliceosomal prp8 Alleles Suggests Functional Correlations with Ribosomal Ambiguity MutantsMolecular Cell 14:343–354https://doi.org/10.1016/S1097-2765(04)00217-5
- CEF1/CDC5 alleles modulate transitions between catalytic conformations of the spliceosomeRNA 18:1001–1013https://doi.org/10.1261/rna.029421.111
- Splicing factor Prp18p promotes genome-wide fidelity of consensus 3′-splice sitesNucleic Acids Research https://doi.org/10.1093/nar/gkad968
- Widespread Impact of Nonsense-Mediated mRNA Decay on the Yeast IntronomeMolecular Cell 31:360–370https://doi.org/10.1016/j.molcel.2008.07.005
- A Conformational Rearrangement in the Spliceosome Sets the Stage for Prp22-Dependent mRNA ReleaseMolecular Cell 30https://doi.org/10.1016/j.molcel.2008.05.003
- Prp22, a DExH-box RNA helicase, plays two distinct roles in yeast pre-mRNA splicingThe EMBO Journal 17:2086–2094https://doi.org/10.1093/emboj/17.7.2086
- A conformational rearrangement in the spliceosome is dependent on PRP16 and ATP hydrolysisThe EMBO Journal 11:5033–5039https://doi.org/10.1002/j.1460-2075.1992.tb05610.x
- RNA helicase dynamics in pre-mRNA splicingThe EMBO Journal 19:6582–6591https://doi.org/10.1093/emboj/19.23.6582
- Spliceosomal DEAH-Box ATPases Remodel Pre-mRNA to Activate Alternative Splice SitesCell 164https://doi.org/10.1016/j.cell.2016.01.025
- Staying on message: Ensuring fidelity in pre-mRNA splicingTrends in Biochemical Sciences 37:263–273https://doi.org/10.1016/j.tibs.2012.04.001
- Functional interactions of Prp8 with both splice sites at the spliceosomal catalytic centerGenes & Development 13:1983–1993
- In Vitro Mutagenesis and Plasmid Shuffling: From Cloned Gene to Mutant YeastMethods in Enzymology 194:302–318https://doi.org/10.1016/0076-6879(91)94023-6
- psiCLIP reveals dynamic RNA binding by DEAH-box helicases before and after exon ligationNature Communications 12https://doi.org/10.1038/s41467-021-21745-9
- Cus2 enforces the first ATP-dependent step of splicing by binding to yeast SF3b1 through a UHM–ULM interactionRNA 25:1020–1037https://doi.org/10.1261/rna.070649.119
- Mechanism of protein-guided folding of the active site U2/U6 RNA during spliceosome activationScience 370https://doi.org/10.1126/science.abc3753
- Preparation of Yeast MediaCurrent Protocols in Molecular Biology 23:13–1https://doi.org/10.1002/0471142727.mb1301s23
- Structure of an Intron Lariat Spliceosome from Saccharomyces cerevisiaeCell 171:120–132https://doi.org/10.1016/j.cell.2017.08.029
- Structural basis for conformational equilibrium of the catalytic spliceosomeMolecular Cell https://doi.org/10.1016/j.molcel.2021.02.021
- Postcatalytic spliceosome structure reveals mechanism of 3′–splice site selectionScience 358:1283–1288https://doi.org/10.1126/science.aar3729
- Methods for merging data sets in electron cryo-microscopyActa Crystallographica Section D: Structural Biology 75:782–791https://doi.org/10.1107/S2059798319010519
- SpliceWiz: Interactive analysis and visualization of alternative splicing in RBriefings in Bioinformatics 25https://doi.org/10.1093/bib/bbad468
- A Fluorescence-Labeling Method for Sequencing Small RNA on Polyacrylamide GelNucleic Acids Research 24:3472–3473https://doi.org/10.1093/nar/24.17.3472
- Competition between the ATPase Prp5 and Branch Region-U2 snRNA Pairing Modulates the Fidelity of Spliceosome AssemblyMolecular Cell 28:838–849https://doi.org/10.1016/j.molcel.2007.09.022
- CUS2, a Yeast Homolog of Human Tat-SF1, Rescues Function of Misfolded U2 through an Unusual RNA Recognition MotifMolecular and Cellular Biology 18:5000–5009https://doi.org/10.1128/MCB.18.9.5000
- Mechanism of exon ligation by human spliceosomeMolecular Cell 82https://doi.org/10.1016/j.molcel.2022.05.021
- Functional and physical interaction between the yeast splicing factors Slu7 and Prp18Nucleic Acids Research 25:2146–2152https://doi.org/10.1093/nar/25.11.2146
- Structures of the human spliceosomes before and after release of the ligated exonCell Research 29https://doi.org/10.1038/s41422-019-0143-x
Article and author information
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
Copyright
© 2024, Senn et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.