3.1 Yeast tolerance to Lithium Chloride is eroded by gene deletions of PEX11 and RIM20
Chemical-genetic approaches have the ability to allow the detailed annotation of the functional properties of various chemicals and bioactive compounds. They can shed light on a compound's primary mode of action as well as its secondary sites of interactions within a cell [24]. The sensitivity of mutant strains to a target compound is a powerful approach for determining the target pathways and the impact a compound has on a cell at the molecular level. In the current study, while investigating yeast gene deletion mutants that are sensitive to LiCl, we discovered two deletion mutants for PEX11 and RIM20, that exhibited increased sensitivity to LiCl (Figure 1A, B) compared to a control strain (Figure 1A). Using a spot test analysis, we demonstrate that when cells are grown in media containing galactose as the primary carbon source, deletion of PEX11 and RIM20 resulted in a substantial reduction in growth in the presence of 10 mM LiCl, indicating increased sensitivity. Deletion of TIF2 (eIF4A) was used as a positive control. Further we demonstrated that incorporating the deleted genes back into the corresponding gene deletion mutants reversed the observed growth reduction establishing a correlation between the observed phenotypes and the deleted gene. Colony count measurement analysis further validated these results by offering a quantitative viewpoint (Figure 1B). The formation of colonies observed in the presence of LiCl is compared and normalized to the number of colonies observed in the control strain. Lower number of colonies formed for gene deletion mutants indicates increased sensitivity to LiCl compared to the control strain. As shown in Figure 1B, deletion of TIF2, PEX11 and RIM20 under the influence of LiCl results in a significant reduction in colony formation. As before reintroduction of the deleted genes back into the corresponding mutant cells increased the number of colonies formed to levels comparable to that of the control strain. Additionally, we examined yeast strains' exposure to LiCl when glucose was used as the carbon source. As expected, no increased sensitivity for the mutant strains was observed when glucose was used (S1. Figure).
Previously, it was stated that overexpression of TIF2 reversed the toxicity of susceptible strains to LiCl. We inserted overexpression plasmids for our candidate genes into corresponding deletion strains spotted on media containing 10 mM LiCl to see whether they would similarly restore LiCl sensitivity. When the plasmids were integrated into the mutant strains, fitness was restored, suggesting that they may function similarly to eIF4A in the cell.
When galactose is used as the carbon source, LiCl inhibits PGM2 expression, resulting in the accumulation of galactose intermediate metabolites that are toxic to yeast cells [11]. GAL1 is the galactokinase in yeast that phosphorylates -D-galactose to -D-galactose-1-phosphate in the initial stages of galactose metabolism. By converting glucose-1-phosphate to glucose-6-phosphate, PGM2 promotes galactose entry into glycolysis. We investigated the impact of PEX11 and RIM20 on LiCl toxicity when galactose metabolism was disabled. For this purpose, we developed double gene deletions for PEX11 and RIM20 with the GAL1 gene. Not surprisingly, we observed GAL1 double mutant cells were no longer hypersensitive to LiCl treatment. PEX11 and RIM20 have not been previously linked to LiCl sensitivity or the underlying molecular mechanisms associated with it, rendering them intriguing gene candidates to investigate.
3.2 PGM2 Expression is Controlled at the Translational Level by PEX11 and RIM20.
PGM2 plays a key role in the sensitivity of yeast cells to LiCl[12], we investigated whether PEX11 and RIM20 will influence PGM2 expression. For this PGM2p was GFP-tagged, and western blot analysis was performed using anti-GFP antibodies to quantify protein content (Figure 2A). The deletion of PEX11 and RIM20 had no substantial effect on PGM2p protein levels in the absence of LiCl treatment. Interestingly, when cells were treated with LiCl, PEX11 and RIM20 gene deletion resulted in a significant reduction of PGM2p protein levels compared to the WT. To analyze the effect of PEX11 and RIM20 gene deletion on PGM2 transcription, the mRNA content of PGM2 was measured using qRT-PCR. Figure 2B indicates that deletion of PEX11 and RIM20 had similar levels of PGM2 mRNA levels as the WT in control media. However, PGM2 mRNA content increased in cells treated with LiCl, but there was no statistically relevant difference between mutant strains and WT. Thus, the deletion of PEX11 and RIM20 seems to have little influence on PGM2 mRNA content. As a result, it appears that PEX11 and RIM20 influence PGM2p at the level of protein synthesis. These findings are similar to the observations by Hajikarimlou et al., where in the presence of LiCl, deletion of YTA6 and YPR096C inhibited PGM2 expression at the translational level [16].
3.3 Deletion of PEX11 and RIM20 Influences Translation of β-Galactosidase Reporter mRNAs with a Hairpin Structure, but not Those Without
Several genes that affected PGM2 expression at the translation level, used the 5’-UTR of PGM2 mRNA to exert their activities. The 5’-UTR of PGM2 mRNA is predicted to have a structured region (S2. Figure) and in the absence of TIF2, a protein that unwinds mRNA structures during translation, PGM2 expression is reduced [14]. Next we intended to investigate whether PEX11 and RIM20 affect translation using the structured 5’-UTR of PGM2 mRNA. For this experiment, we utilized pPGM2 plasmid where the PGM2 5’-UTR is placed in front of a LacZ expression cassette in the p416 expression plasmid. The pPGM2 plasmid, as well as the parental p416 plasmid lacking a structural region in front of the LacZ gene were transformed into the deletion mutant strains for our candidate genes as well as the WT strain. β -galactosidase activity was quantified as a measure of translation (Figure 3A). No significant difference was observed for β -galactosidase activity derived from mRNAs lacking a 5’-UTR structure (p416) in different strains. β -galactosidase activity was significantly reduced however in pex11∆ and rim20∆ when the mRNA carried PGM2 5’-UTR suggesting a connection between the activity of PEX11 and RIM20, and the translation of structured mRNAs.
Next, we examined the effect of PEX11 and RIM20 on additional structured mRNAs. For this we utilized four additional constructs each carrying different structures at their 5’-UTR. pTAR carries a structure derived from HIV1 mRNA containing a ΔG= -57.9kcal/mol and pRTN carries the structured 5’-UTR of RTN4IP1 mRNA with ΔG= -29.8kcal/mol. Similarly, pBcell carries a structure derived from BCL-2 mRNA and a ΔG= -20kcal/mol. The last construct, p2hair contains a synthetic structure designed to have a high degree of complexity ΔG= -33kcal/mol. In our analysis. we observed that deletion of PEX11 and RIM20 significantly reduced the expression of all four highly structured mRNAs, comparing the WT suggesting a general role for these two genes in the translation of structured mRNAs (Figure 4).
3.4 Genetic Interaction Analysis Further Connects the Activity of PEX11 and RIM20 to Protein Biosynthesis
Genetic interaction (GI) analysis is based on the principle that parallel pathways allow for flexibility and tolerance to random harmful mutations, preserving cells viability and maintaining cell homeostasis [16]. A gene in one pathway might compensate for a gene in another, allowing the cell to survive. Consequently, when two genes in parallel pathways are deleted, cell fitness can be unexpectedly reduced (sickness) or even the cell dies (lethality). Therefore, when deletion of two genes result in a phenotype that is unexpected from the phenotypes of individual gene deletions, it is commonly said that the two genes are forming a genetic interaction. Due to the reduced fitness of double mutant this type of interaction is called negative genetic interaction (nGI). Many studies use nGIs to investigate gene function and pathway interactions [15, 16, 29, 37].
The high throughput analysis of GIs in yeast is performed by mating the two yeast mating types: α-mating type (Mat α) and a-mating type (Mat a). Mat "α" carries the target gene deletion and is crossed with an array of single gene deletions of Mat "a," mating type, to produce double gene deletions[28]. Colony size is used to measure the fitness of the strains[35, 38]. We utilised this method to investigate genetic connections between our query genes, PEX11 and RIM20, with approximately 1000 additional genes that includes approximately 700 genes associated with gene expression pathway and a random collection of 300 genes used as a control. (S1. Table).
In our analysis, we discovered numerous intriguing nGIs and several common gene hits between PEX11 and RIM20. Functional enrichment analysis of the hits revealed that a significant number of them are connected to relevant biological processes that includes translation regulation, ribosome biology,, mRNA catabolic processes, and protein synthesis and protein synthesis associated factors (Figure 5). For PEX11, we identified nGIs including DPH5, PCI8, EAP1, and TIF1 among others. The methyltransferase, DPH5 is engaged in the diphthamide biosynthesis pathway, responsible for assembly of translation elongation factor 2 (Eft1p or Eft2p[39]., PCI8 forms the subunit of translation initiation factor, eIF3b, contributing to the formation of 43S pre-initiation complex during translation initiation. It is a key factor in cellular signalling at all stages of protein synthesis, including elongation, termination, and ribosome recycling[40]. EAP1, is also involved in translation initiation process. It is an associated protein of translation initiation factor eIF4E, which forms a part of eIF4F complex. It could potentially enhance/inhibit the overall translation rate by facilitating mRNA degradation and promoting mRNA decapping[41].
Similarly, RIM20 interacted with several translation machinery associated genes including RPS23A, TOR1, and TMA19 (Figure 5). RPS23A codes for the ribosomal protein 28 (rp28) and forms a component of the 40S small ribosomal subunit[42]. Translation regulator protein TOR1 is a member of the TOR complex, which is known to regulate global translation rate through signal transduction. Under stress environments, the TOR-related proteins TOR1 and TOR2 control translation initiation and promote early G1 progression in yeast[43]. TMA19 codes for Translation Machinery Associated protein 19 that interacts with ribosomes during translation initiation. [44].
We also observed several common interactors between PEX11 and RIM20 including SLH1, and CAF20. SLH1 is a putative RNA helicase that is directly implicated in the translation inhibition of non-poly(A) mRNAs. [45]. CAF20 is a translational control phosphoprotein found in the mRNA cap-binding complex. It acts as a regulator of cap-dependent translational initiation [46].
Conditional nGIs are those interactions that are formed under a specific condition that includes the presence of a sub-inhibitory concentration of a bioactive compound, cold shock, heat shock, minimum media etc. They represent gene functional associations that are formed in response to a certain environment[30, 47]. For example, in the presence of DNA damage, the functions of certain genes may be altered, and it is the altered functions that are functionally related and hence form the basis of GIs in the presence of DNA damage[48]. We investigated nGIs for PEX11 and RIM20 using a mild sub-inhibitory concentration of LiCl (3 mM). The observed GIs in this case are those that are formed as a response to the presence of LiCl. As seen in Figure 6, we found new nGIs for our candidate genes PEX11 and RIM20. PEX11 interacted with several translation regulating genes including PBP1, and EBS1 among others. PBP1 is a component of glucose-deprived induced stress granulates that influences global translation rate through the TORC1 signalling pathway and autophagy. It is also known to have an active role in controlling mRNA polyadenylation by interacting with Pab1P, which influences initiation of translation during protein biosynthesis[49]. EBS1 codes for a protein involved in translation inhibition by interacting with cap binding proteins CDC33 and NAM7 upon glucose starvation. It also regulates translation initiation factor eIF4E, another cap binding protein that binds to 5’ end of the mRNA during the formation of 48S initiation complex in translation initiation step[50]. Interestingly, eIF4E is a component of the eIF4F complex, that comprises of eIF4A, eIF4G and eIF4E[51]. eIF4G is a scaffolding protein and eIF4A is a DEAD-box helicase protein that is actively involved in unwinding structured mRNAs during translation initiation. PEX11 interacting with EBS1 may indirectly influence helicase activity of eIF4A and potentially affect other helicases influencing global translation rate.
Similarly, RIM20 interacted with several news genes involved in translation control and regulation of translation. ETT1 is a nuclear protein that regulates overall translation rate and enhances/inhibits protein synthesis by influencing initiation and termination step during translation[52]. Under extreme stress conditions, PTP2 modulates phosphorylation in the MAPK signalling pathway, affecting mRNA decapping and mRNA stability influencing overall translation rate[53]. RPS19B is a protein component of the 40s ribosomal subunit that participates in the assembly of 43s pre-initiation complex by binding matured mRNA to small ribosomal subunit with translation initiation factors eIF1A and eIF3 during protein synthesis[54].
Notably, we also observed several common interactors between PEX11 and RIM20 including DHH1, TIF2, and SHE3 that are involved in translation control. DHH1 is a cytoplasmic DEAD-box helicase that regulates translation process influencing mRNA stability, mRNA degradation and polyadenylation at 5’ end of mRNA[55]. As explained above, translation initiation factor, TIF2 is another DEAD box RNA helicase that facilitates the binding and unwinding of mRNA during translation initiation[51]. SHE3 is a RNA binding protein that constitutes a part of the mRNA localization machinery and plays an active role in recruiting Myo4p-She3 complex that is involved in mRNA export / localization into the cytoplasm, influencing global translation rate[56].
Phenotypic suppression array (PSA) analysis focuses on another kind of interaction in which overexpression of one gene compensates for the absence of another [24, 33, 57, 58]. This is an important type of interaction as it can highlight certain functional overlaps between two genes. To this end we exposed the gene expression mutant arrays to 10 mM LiCl. Some of the investigated strains showed increased sensitivity. Then, by introducing PEX11 and RIM20 overexpression plasmids, we attempted to compensate for the observed sensitive phenotypes. Interestingly, the overexpression of either PEX11 or RIM20 compensated for the sensitivity of two gene deletions for RPS30A and SLF1 (Figure 6A). RPS30A encodes a cytosolic small ribosomal subunit protein that can also bind to mRNAs and contributes to the selection of the correct tRNA molecules[59]. SLF1 is an RNA binding protein that regulates mRNA translation, especially under stress conditions. [60].
The recovery of the sick phenotype in gene deletions for SLF1, RPS30A by introducing RIM20 and PEX11 overexpression plasmid was confirmed using spot test (Figure 6B). This was further verified by cell quantification analysis using colony count (Figure 6C). Using spot test and colony count analysis we verified that the overexpression of RIM20 and PEX11 compensated the sensitivity to LiCl that was seen for deletion strain for RPS30A and SLF1. Both SLF1 and RPS30A have reported genetic interactions with TIF2. Consequently, we also included the deletion strain for TIF2 in the spot test for analysing LiCl sensitivity compensation. Interestingly, when PEX11 and RIM20 are overexpressed in TIF2 deletion, LiCl sensitivity caused by TIF2 deletion was also recovered further connecting the activity of these genes.