Design and Application of Temperature Sensitive Mutants in Essential Factors of RNA Splicing and RNA Interference Pathway in Schizosaccharomyces Pombe

Abstract

Gene deletions are a powerful method to uncover the cellular functions of a given gene in living systems. A limitation to this methodology is that it is not applicable to essential genes. Even for non-essential genes, gene knockouts cause complete absence of gene product thereby limiting genetic analysis of the biological pathway. Alternatives to gene deletions are mutants that are conditional, for e.g, temperature sensitive (ts) mutants are robust tools to understand temporal and spatial functions of genes. By definition, products of such mutants have near normal activity at a lower temperature or near-optimal growth temperature which is called as the permissive temperature and reduced activity at a higher, non-optimal temperature called as the non-permissive temperature. Generation of ts alleles in genes of interest is often time consuming as it requires screening a large population of mutants to identify those that are conditional. Often many essential proteins do not yield ts such alleles even after saturation mutagenesis and extensive screening (Harris et al., 1992; Varadarajan et al., 1996). The limited availability of such mutants in many essential genes prompted us to adopt a biophysical approach to design temperature-sensitive missense mutants in an essential gene of fission yeast. Several studies report that mutations in buried or solvent-inaccessible amino acids cause extensive changes in the thermal stability of proteins and specific substitutions create temperature-sensitive mutants (Rennell et al., 1991; Sandberg et al., 1995). We used the above approach to generate conditional mutants in the fission yeast gene spprp18+encoding an essential predicted second splicing factor based on its homology with human and S. cerevisiae proteins. We have used a missense mutant coupled with a conditional expression system to elucidate the cellular functions of spprp18+. Further, we have employed the same biophysical principle to generate a missense mutant in spago1+ RNA silencing factor that is non-essential for viability but has critical functions in the RNAi pathway of fission yeast. Fission yeast pre-mRNA splicing: cellular functions for the protein factor SpPrp18 Pre-mRNA splicing is an evolutionarily conserved process that excises introns from nascent transcripts. Splicing reactions are catalyzed by the large ribonuclear protein machinery called the spliceosome and occur by two invariant trans-esterification reactions (reviewed in Ruby and Abelson, 1991; Moore et al., 1993). The RNA-RNA, RNA–protein and protein-protein interactions in an assembly of such a large protein complex are numerous and highly dynamic in nature. These interactions in in vitro splicing reactions show ordered recruitment of essential small nuclear ribonucleic particles snRNPs and non–snRNP components on pre-mRNA cis-elements. Further these trans acting factors recognize and poise the catalytic sites in proximity to identify and excise introns. The precision of the process is remarkable given the diversity in architecture for exons and introns in eukaryotic genes (reviewed in Burge et al., 1999; Will and Luhrmann, 2006). Many spliceosomal protein components are conserved across various organisms, yet introns have diverse features with large variations in primary sequence. We hypothesize that co-evolution of splicing factor functions occurs with changes in gene and intron architectures and argue for alternative spliceosomal interactions for spliceosomal proteins that thus enabling splicing of the divergent introns. In vitro biochemical and genetic studies in S. cerevisiae and biochemical studies with human cell lines have indicated that ScPRP18 and its human homolog hPRP18 function during the second catalytic reaction. In S. cerevisiae, ScPrp18 is non-essential for viability at growth temperatures <30°C (Vijayraghavan et al., 1989; Vijayraghavan and Abelson, 1990; Horowitz and Abelson, 1993b). The concerted action of ScSlu7 - ScPrp18 heteromeric complex is essential for proper 3’ss definition during the second catalytic reaction (Zhang and Schwer, 1997; James et al., 2002). These in vitro studies also hinted at a possible intron -specific requirement for ScPrp18 and ScSlu7 factors as they were dispensable for splicing of intron variants made in modified ACT1 intron containing transcripts (Brys and Schwer, 1996; Zhang and Schwer, 1997). A short spacing distance between branch point adenosine to 3’splice site rendered the substrate independent of Prp18 and Slu7 for the second step (Brys and Schwer, 1996; Zhang and Schwer, 1997). Extensive mutational analyses of budding yeast ScPrp18 identified two functional domains and suggested separate roles during splicing (Bacikova and Horowitz, 2002; James et al., 2002). Fission yeast with its genome harboring multiple introns and degenerate splice signals has recently emerged as a unique model to study relationships between splicing factors and their role in genomes with short introns. Previously, studies in our lab had initiated genetic and mutational analysis of S. pombe Prp18, the predicted homolog of budding yeast Prp18. Genetic analysis showed its essentiality, but a set of missense mutants based on studies of budding yeast ScPrp18 (Bacikova and Horowitz, 2002) gave either inactive null or entirely wild type phenotype for the fission yeast protein. In this study, we have extended our previous mutational analysis of fission yeast Prp18 by adopting biophysical and computational approaches to generate temperature-sensitive mutants. A missense mutant was used to understand the splicing functions and interactions of SpPrp18 and the findings are summarized below. Fission yeast SpPrp18 is an essential splicing factor with transcript-specific functions and links efficient splicing with cell cycle progression We initiated our analysis of SpPrp18 by adopting a biophysical approach to generate ts mutants. We used the PREDBUR algorithm to predict a set of buried residues, which when mutated could result in a temperature-sensitive phenotype that complements the null allele at permissive temperature. These predictions are based upon two biophysical properties of amino acids: 1) Hydrophobicity, which is calculated in a window of seven amino acids 2) Hydrophobic moment, which is calculated in a sliding window of nine amino acids in a given protein sequence. Several studies correlate these properties to protein stability and function (Varadarajan et al., 1996). One of the buried residue mutants V194R, in helix 1 of SpPrp18 conferred weak temperature- sensitivity and strong cold-sensitivity even when the protein was over expressed from a plasmid. Through semi-quantitative RT-PCR we showed splicing-defects for tfIId+ intron1 in these cells even when grown at permissive temperature. The primary phenotype was the accumulation of pre-mRNA. Further, we showed this splicing arrest is co-related with reduced levels of SpPrp18 protein, linking protein stability and splicing function. Next we examined the effects of this mutation on function by further reduction of protein levels. This was done by integrating the expression cassette nmt81:spprp18+/spprp18V194R at the leu1 chromosomal locus and by metabolic depletion of the integrated allele. Through RT-PCRs we demonstrated that depletion of wild type or missense protein has intron specific splicing defects. These findings showed its non-global and possibly substrate-specific splicing function. In the affected introns, precursor accumulation is the major phenotype, confirming prior data from our lab that hinted at its likely early splicing role. This contrasts with the second step splicing role of the human or budding yeast Prp18 proteins. Previous data from our lab showed loss of physical interaction between SpPrp18 and SpSlu7 by co-immunoprecipitation studies. This again differs from the strong and functionally important ScPrp18 and ScSlu7 interaction seen in budding yeast. We show the absence of charged residues in SpSlu7 interaction region formed by SpPrp18 helix1 and helix2 which can explain the altered associations for SpPrp18 in fission yeast. Importantly, as the V194R mutation in helix 1 shows splicing defects even at permissive temperature, the data indicate a critical role for helix 1 for splicing interactions, possibly one that bridges or stabilizes the proposed weak association of SpPrp18-SpSlu7 with a yet unknown splicing factor. We also investigated the effects of mutations in other helices; surprisingly we recovered only mutations with very subtle growth phenotypes and very mild splicing defects. Not surprisingly, stop codon at L239 residue predicted to form a truncated protein lacking helices 3, 4 and 5 conferred recessive but null phenotype implicating essential functions for other helices. Other amino acid substitutions at L239 position had near wild type phenotype at 30°C and 37°C. Helix 3 buried residue mutant I259A conferred strong cold-sensitivity when over expressed from plasmid, but semi quantitative analysis indicated no splicing defects for intron1 in the constitutively expressed transcript tfIId+. These findings indicate cold sensitivity either arises due to compromised splicing of yet unknown transcripts or that over-expressed protein has near wild type activity. We find mutations in the helix 5 buried residues L324 also conferred near WT phenotype. Earlier studies in the lab found that substitution of surface residues KR that are in helix 5 with alanine lead to null phenotypes (Piyush Khandelia and Usha Vijayraghavan unpublished data). We report stable expression of all of these mutant proteins; L239A, L239P, L239G, I259A, I259V, L324F, L324A as determined by our immunoblot analysis at 30°C and 37°C. The mild phenotypes of many buried residues can be attributed to orientation of their functional groups into a protein cavity between the helices. Lastly, our microscopic cellular and biochemical analysis of cellular phenotypes of spprp18 mutant provided a novel and direct role of this factor in G1-S transition of cell cycle. Our RT-PCR data suggest spprp18+ is required for efficient splicing of several intron containing transcripts involved in G1-S transition and subsequent activation of MBF complex (MluI cell cycle box-binding factor complex) during S-phase and shows a mechanistic link between cell cycle progression and splicing. A tool to study links between RNA interference, centromeric non-coding RNA transcription and heterochromatin formation S.pombe possesses fully functional RNA interference machinery with a single copy for essential RNAi genes ago1+, dcr1+ and rdp1+. Deletion of any of these genes causes loss of heterochromatinzation with abnormal cytokinesis, cell-cycle deregulation and mating defects (Volpe et al., 2002). In S.pombe, exogenous or endogenously generated dsRNA’s from transcription of centromeric repeats are processed by the RNaseIII enzyme dicer to form siRNA. These siRNA’s are loaded in Ago1 to form minimal RNA induced silencing complex (RISC) complex or specialized transcription machinery complex RNA induced transcriptional silencing (RITS) complex and target chromatin or complementary mRNAs for silencing. Thus as in other eukaryotes, fission yeast cells deploy RNAi mediated silencing machinery to regulate gene-expression and influence chromatin status. Several recent studies point to emerging new roles of RNAi and its association with other RNA processes (Woolcock et al., 2011; Bayane et al., 2008; Kallgren et al., 2014). Many recent reports suggest physical interactions of RISC or RITS and RNA dependent RNA polymerase complex (RDRC) with either some factors of the spliceosomal machinery, heterochromatin machinery (CLRC complex) and the exosome mediated RNA degradation machinery (Bayne et al., 2008 and Chinen et al., 2010 ; Hiriart et al., 2012; Buhler et al., 2008; Bayne et al., 2010 ). Thus we presume conditional alleles in spago1+ will facilitate future studies to probe the genetic network between these complexes as most analyses thus far rely on ago1∆ allele or have been based on proteomic pull down analyses of RISC or RITS complexes. In this study, we employed biophysical and modeling approaches described earlier to generate temperature sensitive mutants in spago1+ and spdcr1+. We tested several mutants for their ability to repress two reporter genes in a conditional manner. Our modeling studies on SpAgo1 PAZ domain indicated structural similarities with human Ago1 PAZ domain. We created site-directed missense mutants at predicted buried residues or in catalytic residues. We also analyzed the effects of random amino acid replacements in specific predicted buried or catalytic residues of SpAgoI. These ago1 mutants were screened as pools for their effects on silencing of GFP or of ura4+ reporter genes. These assays assessed post transcriptional gene silencing (PTGS) or transcriptional gene silencing (TGS) activity of these mutants. We obtained three temperature sensitive SpAgo1 mutants V324G, V324S and L215V while the V324E replacement was a null allele. Based upon our modeling, a likely explanation for the phenotype of these mutants is structural distortion or mis-orientation of the functional groups caused due to these mutations, which affect activity in a temperature dependent manner. This distortion in the PAZ domain may affect binding of siRNA and thereby lead to heterochromatin formation defects that we observed. Our data on the SpAgo1 V324 mutant shows conditional centromeric heterochromatin formation confirmed by semi quantitative RT-PCR for dh transcripts levels that shows temperature dependent increase in these transcripts. We find reduced H3K9Me2 levels at dh locus by chromatin immunoprecipitation (ChIP) assay, linking the association of siRNAs for establishment of heterochromatin at this loci. The data on PTGS of GFP transcripts show SpAgo1 V324G mutation has decreased slicing activity as semi-quantitative RT-PCR for GFP transcripts show increased levels at non permissive temperature. These studies point out the importance of siRNA binding to the PAZ domain and its effect on slicing activity of SpAgo1. The mutations in Y292 showed residue loss of centromeric heterochromatin formation phenotype. Thus, we ascribe critical siRNA binding and 3’ end recognition functions to this residue of SpAgo1. These studies point out functional and structural conservation across hAgo1 and SpAgo1. Adopting the aforementioned biophysical mutational approach, we generated mutants in spdcr1+ and screened for those with conditional activity. Our modeling studies on SpDcr1 helicase domain shows it adopts the conserved helicase domain structure seen for other DEAD Box helicases. Our data on mutational analysis of a conserved buried residue I143 in the walker motif B created inactive protein. The data confirm critical functions for dicer in generation of siRNAs and also in recognition of dsRNA ends. Mutants in buried residues L1130 and I1228 of RNase IIIb domain were inactive and the proximity of these residues to the catalytic core suggest that the critical structural alignment of catalytic residues is indispensable for carrying out dsRNA cleavage to generate siRNAs. We also attribute critical catalytic functions to SpDcr1 D1185 residue for generation of siRNA and heterochromatin formation as measured by our transcriptional gene silencing assay. Our studies employing biophysical and computational approaches to design temperature-sensitive mutants have been successfully applied to an essential splicing factor SpPrp18, which was refractory for ts mutants by other methods. Using a missense mutant, we showed its intron-specific splicing function for subsets of transcripts and deduced that its ubiquitous splicing role is arguable. We have uncovered a link between the splicing substrates of SpPrp18 and direct evidence of splicing based cell cycle regulation, thus providing a mechanistic link to the cell cycle arrest seen in some splicing factor mutants. The same methodology was applied to another important biological pathway, the RNAi machinery, where central factors SpAgoI and SpDcrI were examined We report the first instance of conditional gene silencing tool by designing Ago1 ts mutants which will be useful for future studies of the global interaction network between RNAi and other RNA processing events.