- 1Programa de Biología Celular y Molecular, Instituto de Ciencias Biomédicas, Facultad de Medicina, Universidad de Chile, Santiago, Chile
- 2Departamento de Biología, Facultad de Química y Biología, Universidad de Santiago de Chile (USACH), Santiago, Chile
Background: Proteostasis refers to the processes that regulate the biogenesis, folding, trafficking, and degradation of proteins. Any alteration in these processes can lead to cell malfunction. Protein synthesis, a key proteostatic process, is highly-regulated at multiple levels to ensure adequate adaptation to environmental and physiological challenges such as different stressors, proteotoxic conditions and aging, among other factors. Because alterations in protein translation can lead to protein misfolding, examining how protein translation is regulated may also help to elucidate in part how proteostasis is controlled. Codon usage bias has been implicated in the fine-tuning of translation rate, as more-frequent codons might be read faster than their less-frequent counterparts. Thus, alterations in codon usage due to synonymous mutations may alter translation kinetics and thereby affect the folding of the nascent polypeptide, without altering its primary structure. To date, it has been difficult to predict the effect of synonymous mutations on protein folding and cellular fitness due to a scarcity of relevant data. Thus, the purpose of this work was to assess the effect of synonymous mutations in discrete regions of the gene that encodes the highly-expressed enzyme 3-phosphoglycerate kinase 1 (pgk1) in the fission yeast Schizosaccharomyces pombe.
Results: By means of systematic replacement of synonymous codons along pgk1, we found slightly-altered protein folding and activity in a region-specific manner. However, alterations in protein aggregation, heat stress as well as changes in proteasome activity occurred independently of the mutated region. Concomitantly, reduced mRNA levels of the chaperones Hsp9 and Hsp16 were observed.
Conclusion: Taken together, these data suggest that codon usage bias of the gene encoding this highly-expressed protein is an important regulator of protein function and proteostasis.
Introduction
Synonymous codons encoding the same amino acid, were long thought to be equivalent and interchangeable. Therefore, synonymous mutations were often called “silent.” However, a growing body of evidence has shown that synonymous codon choice in nature is not random and may affect multiple aspects of protein biogenesis in diverse organisms (Chen et al., 2004; Chaney and Clark, 2015; Chaney et al., 2017). Usage frequencies of synonymous codons vary greatly between species and within the same genome, a phenomenon called codon usage bias (CUB). In fast-growing unicellular organisms, the most-frequent codons (called optimal codons) are usually decoded by highly-expressed tRNAs, and the less-frequent codons (called rare or non-optimal codons) are decoded by less-expressed tRNAs (Ikemura, 1985; Kanaya et al., 1999; Quax et al., 2015).
Each organism has a defined CUB, which has been modeled by the co-evolution of the tRNA genes and the mutational rate of the coding genes. In fast-growing organisms, tRNA gene copy number is generally correlated with the abundance of the corresponding tRNA (Ikemura, 1985). In these organisms, CUB can lead to translation optimization or de-optimization of a group of related mRNAs, coordinating their expression (Zhou et al., 2016). Moreover, CUB has been invoked as a major factor in determining both mRNA and protein levels (Zhou et al., 2016). However, the extent to which translation efficiency of an mRNA relates codon bias to protein levels remains unclear.
In recent years, the issue of CUB has become even more relevant as it has been linked to many processes that affect gene expression, such as mRNA stability, translation accuracy as well as protein levels, folding and localization (reviewed in Drummond and Wilke, 2008; Zhang and Ignatova, 2009; Zhang et al., 2010; Subramaniam et al., 2013; Pechmann et al., 2014; Bali and Bebok, 2015; Buhr et al., 2016; Thommen et al., 2017; Bhattacharyya et al., 2018; Hanson and Coller, 2018). However, it is controversially discussed because while it has been established that codon usage is an important factor in estimating protein concentration in yeast (Hiraoka et al., 2009) and trypanosomatids (Jeacock et al., 2018), in humans it has been claimed that translation elongation speed is independent of CUB (Ingolia et al., 2011).
Of the processes mentioned above, the effect of codon choice on codon reading speed is one the least understood. Since protein folding is a co-translational process for many proteins, codon changes that affect codon reading speed can also alter the folding of the encoded protein, leading to abnormal protein functioning (reviewed in Hunt et al., 2014; Chaney and Clark, 2015; Brule and Grayhack, 2017; Rauscher and Ignatova, 2018; Walsh et al., 2020). For example, the synonymous substitution of five low-usage for high-usage codons in a linker region of the fatty acid-binding protein 1 gene in Echinococcus granulosus was shown to affect protein solubility (Cortazzo et al., 2002). Also, it was demonstrated that replacement of synonymous codons in the gene encoding Bos. taurus gamma-B crystallin affects the translation rate, generating different final folding forms of the protein when expressed in Escherichia coli (Buhr et al., 2016). Moreover, the synonymous polymorphisms found in a common haplotype [at 1236C > T (G412G), 2677G > T (A893S), and 3435C > T (I1145I)] in the human multidrug resistance 1 (MDR1) gene have been associated with multidrug resistance shown by cancer cells. This phenotype is attributable to an impairment of folding and activity of the protein, without altering mRNA and protein levels (Kimchi-Sarfaty et al., 2007; Fung and Gottesman, 2009). Recently, Walsh et al. (2020) demonstrated in E. coli that synonymous codon substitutions in the chloramphenicol acetyltransferase gene alter the translation elongation speed, affecting protein folding and leading to an active protein that is more susceptible to degradation. The effect of codon usage in different regions of the protein (e.g., linker region, protein domains, etc.) is still unknown, although it has been proposed that regions between protein domains are encoded by non-optimal codons, while structured domains are encoded by optimal codons (Zhou et al., 2015).
Given the above, and since alterations in the folding of one protein might affect the response of the cell to this abnormality and/or induce modifications in the folding or aggregation of other proteins, it is important to understand the contribution of codon selection on the regulation of protein homeostasis. To add to our knowledge of these processes, we studied the effects of synonymous mutations in the gene encoding 3-phosphoglycerate kinase (pgk1) from the fission yeast Schizosaccharomyces pombe on the properties of the encoded protein (Pgk1), proteostasis, and cell fitness. S. pombe is a model unicellular eukaryotic organism that has been used in the study of mechanisms involved in transcription, translation, and replication (Hoffman et al., 2015). Pgk1 (EC 2.7.2.3) is a highly-expressed monomeric protein that catalyzes the reversible conversion of 1,3-bisphosphoglycerate (1,3-BPGA) to 3-phosphoglycerate (3-PGA) in glycolysis (Banks et al., 1979; Bowler, 2013), and is one of the only two enzymes that generate ATP at the substrate level in this pathway. We selected this gene as a model because (i) almost all codons are optimal, (ii) CUB of pgk1 is conserved among several organisms, and (iii) the gene is highly-expressed. Previous work performed in pgk1 from S. cerevisiae showed that the replacement of a large number of codons with their synonymous counterparts affected mRNA stability and Pgk1 expression (Hoekema et al., 1987); however a detailed analysis of the contribution of each synonymous mutation to the observed results and their effect on yeast physiology remain unexplored.
In this work, we introduced synonymous mutations in short segments along the entire gene encoding Pgk1 and studied the effect of each mutated segment on the encoded protein and its impact on the organism. Depending on the mutated region, we observed subtle effects at the molecular level; however, independent of their location, all mutations altered protein aggregation without affecting Pgk1 levels. Also, some mutated pgk1 changed the response of the cell to heat and misfolding stress. The results obtained in the present study will help to elucidate the role of codon choice on protein aggregation and homeostasis, leading to a better understanding of the implications of synonymous mutations for cell physiology.
Results
The pgk1 gene shows a strong codon usage bias for optimal codons
The protein sequence of S. pombe Pgk1 is 68.36 and 67.07% identical to those of Saccharomyces cerevisiae (PDB: 1QPG_A) and Homo sapiens (PDB: 4AXX_A) respectively, with 99% coverage. Alignments of primary sequences and comparisons of secondary structures revealed a high degree of conservation between these three proteins (Supplementary Figure 1).
We calculated the codon adaptation index (CAI) and tRNA adaptation index (tAI) for the pgk1 gene from S. pombe (Sharp and Li, 1987; Dos Reis et al., 2004). The score measured by CAI ranges from 0 to 1, such that when the CAI value is higher (closer to 1), the gene is more likely to be highly-expressed. The CAI values revealed that most pgk1 codons are optimal, similar to those of genes that codify the most abundant proteins (e.g., ribosomal proteins) (Table 1; Mehdi, 2014; Quax et al., 2015). Only 37 of the 414 codons in the pgk1 sequence are non-optimal (9%), of which 15 are found proximal to the 150 codons encoding the carboxyl terminus of the protein. We found no obvious preference of non-optimal codons within specific regions of the protein (e.g., linkers between domains or inside domains). Additionally, the tAI values indicate that codons present in pgk1 are decoded by tRNAs with high gene copy number, a finding which is concordant with the CAI score. Thus, we conclude that pgk1 codon usage is optimized, as in most highly-expressed proteins. Altogether, these results indicate that pgk1 has a strong bias to use optimal codons, as found in other species, such as unicellular organisms with high proliferation rates (Supplementary Tables 3, 4). In this study, we did not include an analysis of higher eukaryotes genes, because neither global tAI nor CAI are reliable predictors of codon optimality in these multicellular organisms since the tRNA and mRNA expression levels are tissue-specific (Waldman et al., 2010).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1074741/fmicb-13-1074741-HTML/image_m/fmicb-13-1074741-t001.jpg)
Table 1. Codon adaptation index (CAI) and tRNA adaptation index (tAI) of Schizosaccharomyces pombe genes.
To test the potential role of the optimal codon bias in the expression and/or function of Pgk1, we randomly replaced wild-type optimal codons with synonymous non-optimal counterparts (based on the codon use of pgk1 and CUB of S. pombe) in 10 discrete segments covering the entire pgk1 gene (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1074741/fmicb-13-1074741-HTML/image_m/fmicb-13-1074741-g001.jpg)
Figure 1. Scheme of the mutated regions. The upper scheme presents the Pgk1 protein (in gray). The hatched rectangles indicate the protein domains predicted using the Pfam protein family database [3-phosphoglycerate (3-PGA) binding site and ATP-Mg2+ binding site]. Blue and green rectangles indicate α helices and β strands, respectively. Lower schemes indicate the different constructions, with the 7xHis tag at the C-terminus (His, red rectangle). The left column indicates the name and position of the mutated codons (WT correspond to the control without mutations, the numbers indicate the mutated region). Black framed rectangles indicate the mutated regions and yellow lines show the actual mutated codons.
Synonymous mutations located in the linker region alter growth rate
To assess the effect of synonymous codon mutations on the expression and function of Pgk1, 10 independent segments of pgk1 were mutated and replaced in the chromosome by homologous recombination (Figure 1). Thus, the mutated versions of pgk1 remain governed by the wild-type (WT) promoter and regulatory sequences, thereby maintaining the same genetic context as WT pgk1. We discarded the use of an ectopic vector for overexpression of pgk1 because higher levels of an artificially-produced protein may cause aggregation (Fink, 1998; Baneyx and Mujacic, 2004). The codons to be replaced were selected using the tabulations of codon usage provided in the study of Forsburg (Forsburg, 1994; Hiraoka et al., 2009). WT codons (optimal codons in general) were substituted by their non-optimal counterparts by random mutation at their third position in the selected regions, using the primers detailed in Supplementary Table 2. The positions of codon substitutions are indicated in Figure 1 (yellow bars). In the case of amino acids encoded by six synonymous codons (leucine, serine, and arginine), we also mutated the first position to randomly produce all codons with the lowest codon usage. As expected, tAI values for the described mutated sequences were lower than those of the control (Table 2). The DNA sequences of the mutated segments of pgk1 are shown in Table 2. To easily identify the WT and mutated Pgk1 proteins, we introduced a region coding for a His-tag at the 3′ end of each pgk1 gene (WT and mutants). The WT pgk1 plus a His-tag (Pgk1-His WT) was used as a control and named WT. This strain showed the same characteristics as the parental strain in terms of morphology and growth (data not shown). Mutant strains were named mut (xxx-yyy) where xxx and yyy are the position of the first and the last mutated codons in each region, respectively.
As a first approach to assess the effect of synonymous mutations, we determined the cellular fitness (measured as specific growth rate, μ index, under standard conditions) of the mutant strains. Five to ten clones of each mutated region were randomly selected and grown in minimal medium (EMM2). Growth curves of representative clones are shown in Figure 2 (see version of this figure in logaritmic scale in Supplementary Figure 4), whilst the μ index and generation times are shown in Table 3. The generation time of most mutant pgk1 strains (Table 3) was not affected, as seen in the cases of mut(60–85) and mut(112–144) (Figures 2A,B). However, mut(145–175) and mut(177–208) grew slower than the control (Figure 2C and Table 3). mut(321–352) possessed a μ index and generation time similar to the control but with an altered growth curve (Figure 2D and Table 3). No changes were observed in the mutants under microscopic observation or in terms of cell viability (data not shown). Altogether, the different mutations in discrete regions of pgk1 showed no effect on cell growth, except in the cases of mutants mut(145–175), mut(177–208) and to some extent mut(321–352).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1074741/fmicb-13-1074741-HTML/image_m/fmicb-13-1074741-g002.jpg)
Figure 2. Yeast growth in minimal medium is altered in certain strains harboring synonymous mutations. Yeasts were grown in minimal medium (EMM2) under standard conditions for 24 h at 30°C, and OD600 recorded each hour. Readings of (A) mut(60–85), (B) mut(112–142), (C) mut(145–175), and (D) mut(321–352) were plotted and compared to the control (WT). These results are representative of three independent experiments.
Synonymous mutations harbored by mut(145–175) and mut(177–208) are located in the linker region between the 3-PGA-and ADP-Mg2+-binding domains and the hinge region, respectively. In the case of mut(321–352), the mutations are located within the ADP-Mg2+-binding domain (Figure 1). For further analyses, we selected two mutant strains that showed altered growth [mut(145–175) and mut(321–352)] and two strains with no change in cell growth [mut(60–85) and mut(112–144)].
Synonymous mutations in pgk1 do not alter protein levels but affect protein aggregation and enzyme activity
To assess whether synonymous mutations alter the expression of Pgk1, we first measured the levels of pgk1 mRNA. We observed an increase (with low statistical significance) in the mRNA levels of mut(145–175) and mut(321–352), but not the other two mutant strains (Supplementary Figure 2A). However, the levels of Pgk1 protein in all mutant strains were not altered compared to WT (Supplementary Figure 2B). The increase in the mRNA levels for mut(145–175) and mut(321–352) was thus not reflected in an expected increase in the protein levels, indicating that the translation efficiency of these mRNAs may be lowered as a result of the synonymous mutations. One possible explanation for this finding could be that an altered secondary structure of mutant mRNAs might increase their stability but concomitantly decrease their translation efficiency (Hunt et al., 2014). To test whether these synonymous mutations could modify the folding of the encoded proteins, we analyzed their secondary structure using the circular dichroism (CD) spectra of the soluble fraction of Pgk1 (WT and mutants) (Figure 3A). CD analysis of the purified proteins showed a predominant alpha-helix pattern in the control (WT) sample (Greenfield, 2006). Synonymous mutations in codons 60–85, 112–144, or 321–352, but not in codons 145–175, resulted in slightly altered CD spectra. Of the mutants, the spectra of mut(112–144) was the most affected, indicating a probable loss of α-helix content (Greenfield, 2006). In addition, we performed a protease sensitivity assay to estimate the folding status (Angov et al., 2011), obtaining variable results. Pgk1-His mut(60–85) was insensitive to trypsin, but mut(112–142), mut(145–175), and mut(321–352) were slightly sensitive (Supplementary Figure 3). Taken together, these data suggest that synonymous mutations cause subtle changes in the folding of Pgk1 in a manner that is dependent on the mutated region.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1074741/fmicb-13-1074741-HTML/image_m/fmicb-13-1074741-g003.jpg)
Figure 3. Pgk1 folding and activity are slightly altered by synonymous mutations in pgk1. (A) Circular dichroism of the purified Pgk1 from WT and mutant strains. (B) Pgk1 specific activity assayed with 5 nM purified Pgk1 (WT and mutants), in an assay coupled to GAPDH. The results are representative of three independent assays. Data was analyzed by one-way ANOVA (*p < 0.05; **p < 0.02; ****p < 0.0001).
Given that the activity of a protein depends on its structure, we tested whether synonymous mutations produce alterations in the phosphorylation activity of 3-PGA by Pgk1. Activity of purified Pgk1-His from WT and mutant strains was analyzed in an enzymatic assay coupled to glyceraldehyde-3-phosphate dehydrogenase (GAPDH). Pgk1 from strains mut(60–85), mut(145–175), and mut(321–352) showed no change in enzymatic activity compared to WT. However, a marked 60–70% decrease in activity was observed for mut (112–144) (Figure 3B). This result correlates with the result that this strain has the most-affected CD spectrum (Figure 3A). Indeed, mut(112–144) harbors a large number of non-optimal substitutions (22 in total), 12 of which are found in tandem. The fact that these mutations are located in the linker between the two substrate-binding domains (Figure 1) led us to speculate that they might influence the translation rate of this segment, thus altering the folding of the protein.
Taken together, the data obtained thus far show that synonymous mutations cause slight structural variations in this soluble enzyme that may be reflected in changes in protein folding and activity. However, these alterations in the enzymatic activity of Pgk1 had no effect on the O2 consumption rate of the mutant strains (data not shown). This is probably because glycolysis is unaltered in the mutants since the reaction catalyzed by Pgk1 is not a limiting step in the pathway (Bolaños et al., 2008).
Mutations in the pgk1 gene increase protein aggregation
Because folding is a co-translational process for most proteins (Fedorov and Baldwin, 1995; Thommen et al., 2017), we examined whether synonymous mutations in the selected regions could induce Pgk1 aggregation. We determined that 1–2% of total WT Pgk1 is normally found as an aggregate (Figure 4A, quantified in Figure 4B). We observed that all synonymous mutations increased protein aggregation, reaching up to 10% in the case of mut (112–144) (Figure 4A, quantified in Figure 4B). Therefore, a reduced translation rate caused by the replacement of optimal codons by non-optimal counterparts, or a different secondary structure of the mRNA at the mutated sites might induce an altered co-translational protein folding regime in which perhaps the concomitant increase in exposure of buried hydrophobic amino acids might induce an increase in protein aggregation, irrespective of the mutated region (Nackley et al., 2006; Chen et al., 2014).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1074741/fmicb-13-1074741-HTML/image_m/fmicb-13-1074741-g004.jpg)
Figure 4. Synonymous mutations in pgk1 alter protein aggregation, response to heat stress and chaperone expression. (A) Western blot of WT and mutants (αHis) from aggregated (AF) and total (TF) protein fractions isolated from yeasts cultured in EMM2 during the late logarithmic phase in standard conditions. (B) Quantification of the aggregated fraction (aggregated Pgk1 normalized to total Pgk1) and expressed as fold-change compared to WT. Data was analyzed by one-way ANOVA (****p < 0.0001). (C) Response to heat stress. Yeasts were grown in EMM2 at 42°C for 2 h, then shifted to 37°C (chronic heat stress, indicated with orange arrow) for 24 h. Black circles and red squares represent WT and mutant strains, respectively. (D,E) Quantification by real-time PCR of the levels of chaperones Hsp16 (D) and Hsp19 (E) mRNAs from WT and mutants, grown under standard conditions (30°C for 24 h, no stress, black bars) or heat stress condition (2 h at 42°C and 24 h at 37°C, stress, gray bars). The results are representative of three independent experiments, and data was analyzed by two-way ANOVA (p < 0.005).
Synonymous mutations modify the response of cells to protein aggregation
Since WT pgk1 was replaced in the genome by the mutant variants, we envisioned that the cells of the latter should be chronically exposed to an increase in aggregated Pgk1. Thus, we analyzed whether this phenomenon affects protein homeostasis. We tested the cellular response to bortezomib, a proteasome inhibitor that increases protein aggregation (Canfield et al., 2006; Mathiassen et al., 2015). All selected strains harboring synonymous mutations (except mut145–175) grew at a similar rate as the control (vehicle of bortezomib, Figure 5A). Although bortezomib at 100 μM affected cell growth in the WT strain (compare Figure 5A with Figure 5B). See version of these figures in logaritmic scale in Supplementary Figure 6, all mutant strains grew better than this control (Figure 5B). Proteasome activity was measured in control and mutant cells (Table 4) revealing a slight increase (average 10%) in all mutants [except mut(112–144)].
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1074741/fmicb-13-1074741-HTML/image_m/fmicb-13-1074741-g005.jpg)
Figure 5. Response to the proteasome inhibitor bortezomib is improved by synonymous mutations in pgk1. Yeasts were incubated with vehicle (DMSO) (A) or 100 μM bortezomib (B) in EMM2 for 24 h in standard conditions. Blue squares, red triangles, green inverted triangles, purple diamonds, and orange circles represent WT, mut(60–85), mut(112–144), mut(145–175), and mut(321–352) respectively. Growth was monitored by recording OD600 at each hour.
Bortezomib is an inhibitor of the chymotrypsin-like subunit of the proteasome. We speculate that in the mutant cells, the chronic aggregation of Pgk1 might permanently increase the activity of the proteasome, giving rise to resistance of the cell to the drug, although other proteosomal functions might also be influenced by the mutant pgk1 (Canfield et al., 2006; Takeda et al., 2011). Taken together, we speculate that the chronic exposure of the cells to aggregated Pgk1 in the synonymous mutant strains might improve their adaptation to environmental stress.
Synonymous mutations alter the response of cells to stress
Since all the strains analyzed showed altered protein aggregation, and the response to proteasome inhibitor was opposite to that expected, we hypothesized that synonymous mutations might alter the cellular response to other stresses. To test whether synonymous mutations alter cell growth under heat stress, we exposed the cells to 42°C for 2 h followed by mild heat stress over a longer period to allow growth (37°C for 24 h) (Chung et al., 1998; Petersen and Russell, 2016). Under these conditions, most mutant strains showed a slight increase in growth rate during the logarithmic phase, yet all reached a higher OD600 at the stationary phase (Figure 4C). See a version of this figure in logaritmic scale in Supplementary Figure 5 than the control (WT) strain.
In the S. pombe genome, there are at least 17 members of HSP protein families that respond to different stress elicitors, where both chaperones Hsp9 and Hsp16 are highly expressed under heat treatment (Chen et al., 2003). For this reason, we measured mRNA levels of Hsp9 and Hsp16. As expected we observed that under heat stress, levels of Hsp16 mRNA increased up to 3-fold in WT control (Figure 4D). However Hsp9 mRNA showed no statistically significant increase (Figure 4E). Conversely, mRNA levels of both chaperones were lower in the mutant [except mut(60–85)] than in WT strains under both control and stress conditions (Figures 4D,E).
Taken together, as under heat shock the mutant strains grow better than WT, yet the expression of HSP mRNAs is the opposite as expected, we speculate that the exposure of the cells to aggregated Pgk1 might alter the heat shock stress response, reducing the levels of chaperone mRNA. Therefore, since mutant strains respond better than the WT to heat shock, grew better under proteasome inhibition and during the stationary phase (which is an inherently-stressful condition in yeast (Gray et al., 2004; Narayanaswamy et al., 2009; reviewed in Cotto and Morimoto, 1999), we speculate that the chronic exposure to aggregated Pgk1 alerts the mutant cells, enhancing their pre-adaptation to different stresses by an as yet unknown mechanism.
Discussion
Over the years, it has been postulated that the correct folding of a protein is dictated mainly by its primary structure. However, it is also known that several human diseases are caused by protein misfolding and aggregation, without any alterations in their primary amino acid sequence (Sauna and Kimchi-Sarfaty, 2011). Hence, other cellular processes that ensure proper protein folding must be considered. Since folding is a co-translational process for many proteins (Thommen et al., 2017), modification of the translation speed of mRNA by synonymous mutations or changes in the abundance or structure of tRNAs could lead to alternative folding, misfolding, and potential aggregation of the encoded protein. Thus codon usage has been proposed as an important determinant of not only protein folding but also efficiency and accuracy of expression and function. Since the influence of codon usage on protein expression, folding and cellular fitness is poorly understood, we studied the effect of synonymous mutations in the gene encoding the glycolytic enzyme Pgk1. To do so, we replaced the WT chromosomal copy of pgk1 with different mutant versions using a homologous recombination approach, meaning that the mutant proteins would be chronically expressed.
The pgk1 gene from S. pombe is rich in optimal codons (over 90%), as observed for many abundant proteins in other organisms such as S. cerevisiae, E. coli, and Candida albicans (Supplementary Tables 3, 4). It has been reported that codon usage of highly-expressed genes is conserved by natural selection since synonymous mutations could alter processes such as transcription and translation efficiency. For example, Glastad et al. (2015) described that synonymous substitutions may modify DNA methylation patterns of the pgk1 gene, altering transcription efficiency, and thus gene expression. Also, synonymous substitution could alter translation efficiency by modifying RNA structure, as found in the most common mutation in patients of cystic fibrosis, resulting in a translational pause which impairs the expression and proper folding of the encoded protein (Bartoszewski et al., 2010). Our studies of the effect of synonymous replacements in discrete segments of pgk1 did not reveal changes in protein levels (Supplementary Figure 1), although we did observe a subtle increase in the mRNA level in one of the mutated strains (Supplementary Figure 2A). This result is opposite to the findings of Hoekema et al. (1987) who found that synonymous mutations in the entire S. cerevisiae pgk1 gene caused a decrease in mRNA levels attributable to reduced translation speed, which also lowered protein levels. Therefore, we postulate that these effects may vary according to the context and the size of the mutated segments. Evidence supporting this hypothesis is clearly shown as we determined that cell growth is reduced in mutant strains [mut(145–175), mut(177–208) and, to a lesser extent, mut(209–233)] which possess mutations between protein domains (Table 3 and Figure 1). These results agree with previous observation indicating that CUB of regions between protein domains (intrinsically disordered regions) are conserved (low usage codons in general; Zhou et al., 2015). These observations could be explained because such discrete regions might be important in the timing of the folding of domains of the protein, given that the mis-timing of the exit of domains from the ribosome could lead to the abnormal interaction of protein segments. In other studies, differential effects depending on the extent and localization of synonymous mutations have been reported. For example, in experiments carried out in cell free extracts from Neurospora crassa (Yu et al., 2015), researchers found that synonymous mutations in a long segment of the luciferase gene (replacing rare codons for their optimal counterparts) increased luciferase activity (as a result of the increase in luciferase expression), although synonymous replacements on short segments within this region showed that not all mutations caused changes in protein expression. Hence, the position of synonymous codons could be determinant for the correct folding of the encoded protein. Moreover, it was found that mutations in different short segments of the luciferase gene produce a less marked effect compared with the mutation of all segments at the same time (Yu et al., 2015). Recently, Walsh et al. (2020) showed that synonymous substitutions in the gene encoding chloramphenicol acetyltransferase alter the co-translational folding of the protein. These changes produced a protein that was more susceptible to degradation, even though it was properly folded and active. Therefore, intermediates in the folding of the protein are probably more prone to degradation.
This evidence suggests that synonymous substitutions could produce subtle changes, depending on the mutated regions, impacting protein biogenesis in different ways (gene expression, mRNA stability, and/or translation rate), and producing different translational effects. It has been suggested that decreasing the translation elongation rate by the substitution of optimal codons with synonymous but rare codons, might give more time to nascent proteins to adopt a stable tertiary structure before the protein is completely synthesized by the ribosome (Hanson and Coller, 2018). This alteration of the translation rate can affect protein folding positively or negatively (Jacobson and Clark, 2016). Our results reveal that subtle changes in both secondary structure as well as activity of Pgk1 were observed in the synonymous mutants, in a region-specific manner (Figure 3). Nevertheless, protein aggregation increased in all mutant strains assayed (Figure 4A). As previously discussed, we consider that synonymous mutations in the regions studied could lead to the generation of different intermediates in the folding process, shaped by different translational kinetics that may alter the exit of the nascent peptide, exposing hydrophobic regions that increase protein aggregation (Horwich, 2002).
Although aggregation levels represent a small fraction (5–10%) of the total Pgk1 protein (Figure 4A), this enzyme is one of the most abundant proteins, with 310,000–560,000 molecules per cell in S. cerevisiae (Ghaemmaghami et al., 2003; Kulak et al., 2014). Based on such figures, the absolute number of misfolded molecules in the present study should be approximately 15,000–50,000 molecules per cell. Although protein chaperones are critical to the successful folding of many protein complex (Pechmann and Frydman, 2014; Bagchi et al., 2016; Kramer et al., 2019; Liu et al., 2019), we believe that in the mutant cells, the large amount of misfolded Pgk1 overcomes the function of chaperones. These numerous aggregated molecules may also interact with other proteins and interfere with various processes, for example increasing the malfunction or even triggering aggregation of other proteins interacting with misfolded stretches of Pgk1.
As the expression of mutant Pgk1 is chronic in these cells, this could lead to an adaptation to the stress caused by misfolded proteins, potentially causing the effects observed at the physiological level (Geiler-Samerotte et al., 2011) in response to different stresses, such as heat shock conditions (Figures 4C, 5). Although in the cells carrying mutant Pgk1, the levels of Hsp 16 and 9 chaperone mRNAs were lower than in the control, the response to bortezomib and the increased proteasome activity indicate possible compensatory mechanisms to improve the adaptation of the cells to misfolded proteins (Figures 4D,E). These contradictory observations could be the result of hormetic effects (Zemva et al., 2017), among other factors, involving compensatory events such as potential increases in transcription of those genes encoding proteasome proteases. Similarly, synonymous mutations of codons present in six highly-expressed genes in E. coli were shown to alter the overall translation efficiency of the rest of the transcriptome, due to the heightened demands on the translation machinery (Banks et al., 1979; Bowler, 2013; Hoffman et al., 2015). It has also been reported that aged organisms show depressed chaperone levels and activity (Norry and Loeschcke, 2003; Soti and Csermely, 2007) due to damaged proteins, impaired protein synthesis, and other alterations. Thus, S. pombe strains harboring the synonymous mutations could be more resistant to stresses due to a premature aging process produced by chronic Pgk1 aggregation.
Pgk1 is not exclusively a glycolytic enzyme; Pgk1 is also found in the nucleus and, during hypoxia, in mitochondria where it participates in DNA replication and regulation of the tricarboxylic acid cycle, respectively (Sun et al., 2015; Li et al., 2016). As several enzymes that participate in the phosphorylation of small metabolites, Pgk1 also function as a protein kinase in a number of fundamental cellular processes. This effect is known as moonlighting (reviewed in Lu and Hunter, 2018) Therefore, other Pgk1 functions are potentially affected by aggregation. Such multiple roles of Pgk1 must be considered in order to understand the effects of synonymous mutations in pgk1 on cellular response to stresses, as a measurement of cell fitness. It has been reported that in humans, Pgk1 deficiency is an X-linked metabolic disorder produced by a non-synonymous mutation in pgk1 that leads to altered function or even misfolding (Zoller, 1991). Therefore, the conserved CUB of pgk1 across the phylogenetic tree could be important in the regulation of pgk1 expression and function, where any synonymous substitution in the coding sequence could exert a pleiotropic effect, as observed in this work. Because the many functions of Pgk1, it would be very important to analyze the effect of synonymous mutations in pgk1 from different organisms as other yeast as well as more complex organisms such as mammals and even humans. The data obtained in this work led us to predict that synonymous mutations might exert changes in the physiology of these organisms at different levels.
Conclusion
The replacement of preferred codons by rare codons in discrete segments along the entire sequence of pgk1 from Schizosaccharomyces pombe generated subtle differences in the function of the expressed protein, particularly in the activity of the enzyme. However, irrespective of the location, every mutated segment induced aggregation of the protein. Therefore, the chronic expression of altered Pgk1 may make the cells permanently exposed to stress and thus more resistant to additional stresses.
Materials and methods
Calculation of codon and tRNA adaptation indexes
tRNA (tAI) and codon (CAI) adaptation indexes were calculated as described previously (Sharp and Li, 1987; Mehdi et al., 2014), using the data available in http://gtrnadb2009.ucsc.edu/ and http://www.kazusa.or.jp/codon/, respectively.
Yeast strains and media
The Schizosaccharomyces pombe 972 h- (ATCC 24843) strain was used for homologous recombination experiments. Primers and plasmids used in this study are listed in Supplementary Table 1. S. pombe was grown on YES medium (5 g/l yeast extract, 30 g/l glucose), YPD (10 g/l yeast extract, 20 g/l peptone, 20 g/l glucose), and EMM2 minimal medium (US Biological) supplemented with 20 g/l glucose.
Plasmid construction
The synonymous substitutions were performed based on S. pombe codon usage, replacing optimal codons with their non-optimal counterpart (according to Forsburg, 1994). The endogenous pgk1 gene was replaced by homologous recombination with the wild-type or mutated pgk1, both with a 7xHis-tag at the 3′ end. For homologous recombination, the flanking regions and pgk1 coding sequence were cloned in the pFA6a-KanMX6 vector (Addgene), as described below. Flanking regions were amplified from S. pombe genomic DNA (gDNA), using Herculase II Fusion DNA polymerase (Agilent Genomics) according to the manufacturer’s instructions. For the 5′ flanking region, amplification was performed in two rounds of PCR (Zoller, 1991) to eliminate an internal restriction site (NdeI), using 5′F and 5′R-intern and 5′F-intern and 5′R primers (Supplementary Table 1). For the 3′ flanking region, 3′F and 3′R primers were used (Supplementary Table 1). PCR products were cloned in pGEMT-easy (Promega, United States) and the recombinant plasmids were digested with PfoI and NdeI (for the 5′ flanking region) and SacI and EcoRI (for the 3′ flanking region). The digestion products were ligated into pFA6a-KanMX6 (first the 5′ flanking region, followed by the 3′ flanking region), that had been previously digested with the corresponding restriction enzymes. The final vector was named pFA6a-KanMX6-5′-3′ and confirmed by DNA sequencing. The WT pgk1 coding sequence was amplified from gDNA using the following primers: PGK1-F and PGK1-R (with a 7xHis tag sequence at the 3′ end) (Supplementary Table 1). The product was cloned into pGEMT-easy, and the construction was used as template for the amplification of the WT and mutant sequences. All pgk1 mutants were constructed as follows: Firstly, the amplification of the mutant pgk1 gene sequence was carried out in two separate PCR reactions using (i) PGK1-F primer and the corresponding reverse internal primer containing the mutated codons (for example, for the mutation of codons 112–142, the corresponding primer was 112–142-R) (Supplementary Table 2) to amplify one segment (in this case, 112–142-i) and (ii) PGK1-R and the corresponding forward primer containing the mutated codons (for the same example, 112–142-F) to amplify the second fragment (in this case, 112–142-ii). PCR products were purified and then used as template for two separate amplifications, where 112–142-i was amplified with the PGK1-F primer, and 112–142-ii was amplified with the PGK1-R primer. Secondly, both products were mixed and incubated in a thermocycler in the presence of 0.2 mM dNTPs, 1X Herculase II Fusion DNA polymerase buffer, and Herculase II Fusion DNA polymerase (30 cycles of: 30 s at 95°C, 30 s at 30°C; 1 min at 72°C). Thirdly, the final product was purified and digested with NdeI and BamHI, and ligated into pFA6a-KanMX6-5′-3′. Finally, E. coli JM109 strain was transformed with the corresponding construction by chemical transformation (Chung et al., 1989) and 3–5 clones were analyzed by DNA sequencing to corroborate the correct incorporation of mutations.
Homologous recombination
Purified plasmids containing the flanking regions plus the His-tagged WT and mutant coding sequences of pgk1 were used for amplification of the sequence employed for homologous recombination, using 5′F and 3′R primers (Supplementary Table 1). The corresponding PCR product (100 ng) was electroporated into 100 μl electrocompetent S. pombe prepared as described previously (Forsburg and Rhind, 2006). Yeasts were recovered in 1 M sorbitol at 30°C for 3 h, and then cells were pelleted, resuspended in 200 μl YE medium, and plated on YPD agar in the presence of 200 μg/ml G418 antibiotic (Sigma-Aldrich). Cells were grown for 4–5 days at 30°C and single colonies were picked and grown on YPD medium supplemented with 200 μg/ml G418 antibiotic. To corroborate the incorporation of synonymous mutations, the His-tagged pgk1 coding sequence (WT and mutants) were sequenced after amplification of the gDNA of 5–15 clones for each mutated region, using PGK1-F and PGK1-R primers. Additionally, to check that the recombination occurred in the endogenous pgk1 gene, we designed primers that amplify the locus segment outside the recombined segment (Supplementary Table 1), Forward_Out_PGKF and Reverse_Out_PGKR, combined with internal primers within the segment.
Growth assays
Yeasts were grown overnight in EMM2 medium at 30°C. Then, cells were pelleted (2,250 × g for 5 min) and resuspended in EMM2 at 0.05 OD600 (measured in an Infinite 200 PRO microplate plate reader, TECAN). Yeast cultures (200 μl) were inoculated in a 96-well plate (in triplicate), and cells were incubated for 24 h in constant agitation at 30°C in the plate reader, recording OD600 each hour. For thermal stress, cells were incubated at 42°C for 2 h, and then at 37°C for 22 h. For growth in the presence of the proteasome inhibitor Bortezomib (Santa Cruz Biotechnology), cells were prepared as previously described at a final OD600 of 0.05. Bortezomib was then added at a final concentration of 100 μM, or an equal quantity of DMSO (solvent) and cells were inoculated in 96-well plates in triplicate for each condition. Cells were grown for 24 h, and the OD600 was recorded. For the specific growth rate (μ index), the OD600 at the exponential phase of each strain was determined, and the values inserted into the following equation:
where k is the specific growth rate, N1 and N2 are the initial and terminal OD600 recordings, respectively, and t1 and t2 correspond to initial and terminal times, respectively (in hours). Thus, k is expressed in h−1. The generation (or doubling) time was calculated according to the following equation:
where generation time is expressed in hours.
Messenger RNA isolation and quantification
For pgk1 mRNA quantification, WT and mutant strains were grown to the late exponential phase (OD600, 1.0) in 20 ml EMM2 under standard conditions. For chaperone mRNA quantification, cells were grown in 20 ml EMM2 for 24 h under standard conditions or under heat stress, as described above. Yeasts were then pelleted (2,250 × g for 5 min) and washed once with sterilized water. The cell wall was disrupted with 0.12 μg zymolyase 20 T (US Biological) in 1 M sorbitol for 30 min at 37°C. Cells were pelleted and resuspended in 200 μl TRIzol (Thermo Scientific) and vortexed three times for 1 min each (intercalated with 1 min on ice). Then, 40 μl chloroform was added to the mix and immediately centrifuged at 12,000 × g for 20 min at 4°C. The supernatant was precipitated with 0.7 volumes of isopropanol at –80°C overnight. RNA was pelleted at 15,000 × g for 30 min and washed once with 80% ethanol. RNA was resuspended in 30 μl sterilized miliQ water and quantified in a nanodrop (BioTek), then visualized in an 1% agarose gel to check integrity. RNA (1 μg) was treated with DNAseI (Roche), according to the instructions provided in the manual, and then 500 ng RNA were used for cDNA synthesis followed by Real-Time PCR with the Brilliant II QRT-PCR, AffinityScript Two-Step Master Mix (Agilent), according to manufacturer instructions. The primers used for Real-Time PCR are listed in Supplementary Table 1 (PGK1RT-F and PGK1RT-R for pgk1 amplification). Quantifications were performed using actin as a control (ActRT-F and ActRT-R for actin amplification).
Isolation and quantification of total and aggregated protein
WT and mutant strains were grown to the exponential phase (OD600, 1.0) in 20 ml EMM2. Yeasts were harvested by centrifugation (2,250× g for 5 min), and the pellet was used to analyze protein aggregation as described previously (Rand and Grant, 2006). Briefly, spheroplasts were pelleted and resuspended in lysis buffer [50 mM potassium phosphate buffer, pH 7, 1 mM EDTA, 5% glycerol, 1 mM phenylmethylsulfonyl fluoride and Complete Mini protease inhibitor cocktail (Roche)]. Cell disruption was carried out by three vortex cycles (1 min of vortex and 1 min on ice) with 220 mg of acid-washed glass beads (Sigma-Aldrich G8772). Membrane proteins were removed by washing twice with 320 μl lysis buffer and 80 μl 10% NP-40 (Sigma-Aldrich), and the final aggregated protein extract was resuspended in 20 μl 1X loading buffer. Total and aggregated protein extracts were analyzed by Western blotting using an anti His-tag antibody (His Tag Antibody MAB050R-100, R&D Systems). Western blots against tubulin were performed as internal controls using T5168 monoclonal anti-a-tubulin clone B-5-1-2 (Sigma-Aldrich). Bands corresponding to WT, mutants, and tubulin (for total protein), or WT and mutants (for aggregated proteins) were quantified using ImageJ software.
Protein purification
Yeasts (WT and mutants) were grown at 30°C until the stationary phase, and then soluble proteins were isolated as described previously (Rand and Grant, 2006), scaled to a 160 ml culture (EMM2). Soluble protein was incubated for 2 h with 250 μl Ni-PentaTM Agarose 6 fast flow resin (Marvelgent Biosciences Inc), previously equilibrated with lysis buffer at 4°C. The resin was then washed with 20 volumes of washing buffer (20 mM HEPES pH 7.0, 300 mM NaCl, 1 mM β-mercaptoethanol) supplemented with 2 mM imidazole. His-tagged proteins were eluted with 10 volumes of elution buffer (washing buffer supplemented with 25 mM imidazole), and collected in 13 fractions (100 μl each). Fractions were analyzed by SDS-PAGE, and fractions 2–8 were selected to be dialyzed with the corresponding buffer for all yeasts assayed.
Circular dichroism
Purified proteins (WT and mutants) were dialyzed against 20 mM phosphate buffer, pH 7.0, at 4°C overnight. Purified proteins (300–500 mg/ml) were used for CD analyses in a Jasco J-1500 spectrometer. CD spectra were recorded between 200 and 250 nm, at a data pitch of 1.0 nm and a scanning speed of 50 nm/min. Three accumulations were obtained for each measurement, and the data were processed with the following equation:
Enzyme activity
Purified WT and mutant proteins were dialyzed against activity buffer (25 mM Tris–HCl, pH 7.4, 20% glycerol, 0.5 mM MgCl2, 1 mM phenylmethylsulfonyl fluoride) at 4°C overnight. The activity assay was carried out with 5 nM of each protein, in the presence of 7 mM 3-PGA, 3 mM ATP, 1 mM MgCl2, 50 mM Tris–HCl pH 7.4, and 0.5 mM NADH, in a final volume of 200 μl. The reaction mixture was placed in a 96-well plate and incubated at 25°C, measuring the absorbance at 340 nm in an Infinite 200 PRO microplate plate reader every 15 s.
Proteasome activity
Proteasome activity was measured with the Proteasome Activity Assay Kit (Fluorometric) (Abcam, ab0107921). Total protein extracts were obtained as described previously (Rand and Grant, 2006) from 5 ml culture of each strain (OD600 1.0). Then, 1 mg of each protein extract was mixed with Assay Buffer in a final volume of 100 μl, in independent wells of a 96-well dark plate (for fluorometric assays), in triplicate. Proteasome substrate (1 μl, 5 mM) was added to each well, and then proteasome activity was determined as the change in fluorescence, measuring excitation/emission at 350/440 nm in the Infinite 200 PRO microplate plate reader at 37°C for 25 min. A calibration curve was also generated, as indicated by the manufacturer, in order to determine the product concentration in the assays. One unit of proteasome activity corresponds to the amount of proteasome which generates 1 nmol of product per minute at 37°C.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
SM-R designed and performed most of the experiments and participated in writing the manuscript. LA carried out some experiments and discussed results. AK participated in the design and discussion of results. RF participated in the discussion and writing. GL performed some experiments and discussed results. OO provided the funding, designed the experimental strategy, discussed results, and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by FONDECYT, Chile Postdoctoral 2014 Grant 3150366 (SM-R) and FONDECYT Regular Grants 1110834 and 1190552 (OO), 1211386 (GL), 1191074 (AK) and Universidad de Chile and Universidad de Santiago de Chile. LA was the recipient of a fellowship for graduate studies from CONICYT, Chile.
Acknowledgments
We would like to thank Maira Rivera from the laboratory of Mauricio Baez at the Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile for assistance in the circular dichroism analysis (FONDEQUIP EQM140151).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.1074741/full#supplementary-material
References
Angov, E., Legler, P. M., and Mease, R. M. (2011). Adjustment of codon usage frequencies by codon harmonization improves protein expression and folding. Methods Mol. Biol. 705, 1–13. doi: 10.1007/978-1-61737-967-3_1
Bagchi, P., Inoue, T., and Tsai, B. (2016). EMC1-dependent stabilization drives membrane penetration of a partially destabilized non-enveloped virus. eLife 5:e21470. doi: 10.7554/eLife.21470
Bali, V., and Bebok, Z. (2015). Decoding mechanisms by which silent codon changes influence protein biogenesis and function. Int. J. Biochem. Cell Biol. 64, 58–74. doi: 10.1016/j.biocel.2015.03.011
Baneyx, F., and Mujacic, M. (2004). Recombinant protein folding and misfolding in Escherichia coli. Nat. Biotechnol. 22, 1399–1408. doi: 10.1038/nbt1029
Banks, R. D., Blake, C. C., Evans, P. R., Haser, R., Rice, D. W., Hardy, G. W., et al. (1979). Sequence, structure and activity of phosphoglycerate kinase: a possible hinge-bending enzyme. Nature 279, 773–777. doi: 10.1038/279773a0
Bartoszewski, R. A., Jablonsky, M., Bartoszewska, S., Stevenson, L., Dai, Q., Kappes, J., et al. (2010). A synonymous single nucleotide polymorphism in DeltaF508 CFTR alters the secondary structure of the mRNA and the expression of the mutant protein. J. Biol. Chem. 285, 28741–28748. doi: 10.1074/jbc.M110.154575
Bhattacharyya, S., Jacobs, W. M., Adkar, B. V., Yan, J., Zhang, W., and Shakhnovich, E. I. (2018). Accessibility of the Shine-Dalgarno sequence dictates N-terminal codon bias in E. coli. Mol. Cell 70, 894–905.e5. doi: 10.1016/j.molcel.2018.05.008
Bolaños, J. P., Delgado-Esteban, M., Herrero-Mendez, A., Fernandez-Fernandez, S., and Almeida, A. (2008). Regulation of glycolysis and pentose-phosphate pathway by nitric oxide: impact on neuronal survival. Biochim. Biophys. Acta 1777, 789–793. doi: 10.1016/j.bbabio.2008.04.011
Bowler, M. W. (2013). Conformational dynamics in phosphoglycerate kinase, an open and shut case? FEBS Lett. 587, 1878–1883. doi: 10.1016/j.febslet.2013.05.012
Brule, C. E., and Grayhack, E. J. (2017). Synonymous codons: choose wisely for expression. Trends Genet. 33, 283–297. doi: 10.1016/j.tig.2017.02.001
Buhr, F., Jha, S., Thommen, M., Mittelstaet, J., Kutz, F., Schwalbe, H., et al. (2016). Synonymous codons direct cotranslational folding toward different protein conformations. Mol. Cell 61, 341–351. doi: 10.1016/j.molcel.2016.01.008
Canfield, S. E., Zhu, K., Williams, S. A., and McConkey, D. J. (2006). Bortezomib inhibits docetaxel-induced apoptosis via a p21-dependent mechanism in human prostate cancer cells. Mol. Cancer Ther. 5, 2043–2050. doi: 10.1158/1535-7163.MCT-05-0437
Chaney, J. L., and Clark, P. L. (2015). Roles for synonymous codon usage in protein biogenesis. Annu. Rev. Biophys. 44, 143–166. doi: 10.1146/annurev-biophys-060414-034333
Chaney, J., Steele, A., Carmichael, R., Rodriguez, A., Specht, A. T., Ngo, K., et al. (2017). Widespread position-specific conservation of synonymous rare codons within coding sequences. PLoS Comput. Biol. 13:e1005531. doi: 10.1371/journal.pcbi.1005531
Chen, S. L., Lee, W., Hottes, A. K., Shapiro, L., and McAdams, H. H. (2004). Codon usage between genomes is constrained by genome-wide mutational processes. Proc. Natl. Acad. Sci. U. S. A. 101, 3480–3485. doi: 10.1073/pnas.0307827100
Chen, D., Toone, W. M., Mata, J., Lyne, R., Burns, G., Kivinen, K., et al. (2003). Global transcriptional responses of fission yeast to environmental stress. Mol. Biol. Cell 14, 214–229. doi: 10.1091/mbc.e02-08-0499
Chen, Y., Zhao, H., Schuck, P., and Wistow, G. (2014). Solution properties of γ-crystallins: compact structure and low frictional ratio are conserved properties of diverse γ-crystallins. Protein Sci. 23, 76–87. doi: 10.1002/pro.2395
Chung, K. S., Hoe, K. L., Kim, K. W., and Yoo, H. S. (1998). Isolation of a novel heat shock protein 70-like gene, pss1+ of Schizosaccharomyces pombe homologous to hsp110/SSE subfamily. Gene 210, 143–150. doi: 10.1016/s0378-1119(98)00061-4
Chung, C. T., Niemela, S. L., and Miller, R. H. (1989). One-step preparation of competent Escherichia coli: transformation and storage of bacterial cells in the same solution. Proc. Natl. Acad. Sci. U. S. A. 86, 2172–2175. doi: 10.1073/pnas.86.7.2172
Cortazzo, P., Cerveñansky, C., Marín, M., Reiss, C., Ehrlich, R., and Deana, A. (2002). Silent mutations affect in vivo protein folding in Escherichia coli. Biochem. Biophys. Res. Commun. 293, 537–541. doi: 10.1016/S0006-291X(02)00226-7
Cotto, J. J., and Morimoto, R. I. (1999). Stress-induced activation of the heat-shock response: cell and molecular biology of heat-shock factors. Biochem. Soc. Symp. 64, 105–118.
Dos Reis, M., Savva, R., and Wernisch, L. (2004). Solving the riddle of codon usage preferences: a test for translational selection. Nucleic Acids Res. 32, 5036–5044. doi: 10.1093/nar/gkh834
Drummond, D. A., and Wilke, C. O. (2008). Mistranslation-induced protein misfolding as a dominant constraint on coding-sequence evolution. Cells 134, 341–352. doi: 10.1016/j.cell.2008.05.042
Fedorov, A. N., and Baldwin, T. O. (1995). Contribution of cotranslational folding to the rate of formation of native protein structure. Proc. Natl. Acad. Sci. U. S. A. 92, 1227–1231. doi: 10.1073/pnas.92.4.1227
Fink, A. L. (1998). Protein aggregation: folding aggregates, inclusion bodies and amyloid. Fold. Des. 3, R9–R23. doi: 10.1016/S1359-0278(98)00002-9
Forsburg, S. L. (1994). Codon usage table for Schizosaccharomyces pombe. Yeast 10, 1045–1047. doi: 10.1002/yea.320100806
Forsburg, S. L., and Rhind, N. (2006). Basic methods for fission yeast. Yeast 23, 173–183. doi: 10.1002/yea.1347
Fung, K. L., and Gottesman, M. M. (2009). A synonymous polymorphism in a common MDR1 (ABCB1) haplotype shapes protein function. Biochim. Biophys. Acta 1794, 860–871. doi: 10.1016/j.bbapap.2009.02.014
Geiler-Samerotte, K. A., Dion, M. F., Budnik, B. A., Wang, S. M., Hartl, D. L., and Drummond, D. A. (2011). Misfolded proteins impose a dosage-dependent fitness cost and trigger a cytosolic unfolded protein response in yeast. Proc. Natl. Acad. Sci. U. S. A. 108, 680–685. doi: 10.1073/pnas.1017570108
Ghaemmaghami, S., Huh, W. K., Bower, K., Howson, R. W., Belle, A., Dephoure, N., et al. (2003). Global analysis of protein expression in yeast. Nature 425, 737–741. doi: 10.1038/nature02046
Glastad, K. M., Goodisman, M. A., Yi, S. V., and Hunt, B. G. (2015). Effects of DNA methylation and chromatin state on rates of molecular evolution in insects. G3 6, 357–363. doi: 10.1534/g3.115.023499
Gray, J. V., Petsko, G. A., Johnston, G. C., Ringe, D., Singer, R. A., and Werner-Washburne, M. (2004). Sleeping beauty: quiescence in Saccharomyces cerevisiae. Microbiol. Mol. Biol. Rev. 68, 187–206. doi: 10.1128/MMBR.68.2.187-206.2004
Greenfield, N. J. (2006). Using circular dichroism spectra to estimate protein secondary structure. Nat. Protoc. 1, 2876–2890. doi: 10.1038/nprot.2006.202
Hanson, G., and Coller, J. (2018). Codon optimality, bias and usage in translation and mRNA decay. Nat. Rev. Mol. Cell Biol. 19, 20–30. doi: 10.1038/nrm.2017.91
Hiraoka, Y., Kawamata, K., Haraguchi, T., and Chikashige, Y. (2009). Codon usage bias is correlated with gene expression levels in the fission yeast Schizosaccharomyces pombe. Genes Cells 14, 499–509. doi: 10.1111/j.1365-2443.2009.01284.x
Hoekema, A., Kastelein, R. A., Vasser, M., and de Boer, H. A. (1987). Codon replacement in the PGK1 gene of Saccharomyces cerevisiae: experimental approach to study the role of biased codon usage in gene expression. Mol. Cell. Biol. 7, 2914–2924. doi: 10.1128/mcb.7.8.2914-2924.1987
Hoffman, C. S., Wood, V., and Fantes, P. A. (2015). An ancient yeast for young geneticists: a primer on the Schizosaccharomyces pombe model system. Genetics 201, 403–423. doi: 10.1534/genetics.115.181503
Horwich, A. (2002). Protein aggregation in disease: a role for folding intermediates forming specific multimeric interactions. J. Clin. Invest. 110, 1221–1232. doi: 10.1172/JCI16781
Hunt, R. C., Simhadri, V. L., Iandoli, M., Sauna, Z. E., and Kimchi-Sarfaty, C. (2014). Exposing synonymous mutations. Trends Genet. 30, 308–321. doi: 10.1016/j.tig.2014.04.006
Ikemura, T. (1985). Codon usage and tRNA content in unicellular and multicellular organisms. Mol. Biol. Evol. 2, 13–34. doi: 10.1093/oxfordjournals.molbev.a040335
Ingolia, N. T., Lareau, L. F., and Weissman, J. S. (2011). Ribosome profiling of mouse embryonic stem cells reveals the complexity and dynamics of mammalian proteomes. Cells 147, 789–802. doi: 10.1016/j.cell.2011.10.002
Jacobson, G. N., and Clark, P. L. (2016). Quality over quantity: optimizing co-translational protein folding with non-‘optimal’ synonymous codons. Curr. Opin. Struct. Biol. 38, 102–110. doi: 10.1016/j.sbi.2016.06.002
Jeacock, L., Faria, J., and Horn, D. (2018). Codon usage bias controls mRNA and protein abundance in trypanosomatids. eLife 7:e32496. doi: 10.7554/eLife.32496
Kanaya, S., Yamada, Y., Kudo, Y., and Ikemura, T. (1999). Studies of codon usage and tRNA genes of 18 unicellular organisms and quantification of Bacillus subtilis tRNAs: gene expression level and species-specific diversity of codon usage based on multivariate analysis. Gene 238, 143–155. doi: 10.1016/s0378-1119(99)00225-5
Kimchi-Sarfaty, C., Oh, J. M., Kim, I. W., Sauna, Z. E., Calcagno, A. M., Ambudkar, S. V., et al. (2007). A “silent” polymorphism in the MDR1 gene changes substrate specificity. Science 315, 525–528. doi: 10.1126/science.1135308
Kramer, G., Shiber, A., and Bukau, B. (2019). Mechanisms of cotranslational maturation of newly synthesized proteins. Annu. Rev. Biochem. 88, 337–364. doi: 10.1146/annurev-biochem-013118-111717
Kulak, N. A., Pichler, G., Paron, I., Nagaraj, N., and Mann, M. (2014). Minimal, encapsulated proteomic-sample processing applied to copy-number estimation in eukaryotic cells. Nat. Methods 11, 319–324. doi: 10.1038/nmeth.2834
Li, X., Jiang, Y., Meisenhelder, J., Yang, W., Hawke, D. H., Zheng, Y., et al. (2016). Mitochondria-translocated PGK1 functions as a protein kinase to coordinate glycolysis and the TCA cycle in tumorigenesis. Mol. Cell 61, 705–719. doi: 10.1016/j.molcel.2016.02.009
Liu, K., Maciuba, K., and Kaiser, C. M. (2019). The ribosome cooperates with a chaperone to guide multi-domain protein folding. Mol. Cell 74, 310–319.e7. doi: 10.1016/j.molcel.2019.01.043
Lu, Z., and Hunter, T. (2018). Metabolic kinases moonlighting as protein kinases. Trends Biochem. Sci. 43, 301–310. doi: 10.1016/j.tibs.2018.01.006
Mathiassen, S. G., Larsen, I. B., Poulsen, E. G., Madsen, C. T., Papaleo, E., Lindorff-Larsen, K., et al. (2015). A two-step protein quality control pathway for a Misfolded DJ-1 variant in fission yeast. J. Biol. Chem. 290, 21141–21153. doi: 10.1074/jbc.M115.662312
Mehdi, A. M., Patrick, R., Bailey, T. L., and Bodén, M. (2014). Predicting the dynamics of protein abundance. Mol. Cell. Proteomics 13, 1330–1340. doi: 10.1074/mcp.M113.033076
Nackley, A. G., Shabalina, S. A., Tchivileva, I. E., Satterfield, K., Korchynskyi, O., Makarov, S. S., et al. (2006). Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure. Science 314, 1930–1933. doi: 10.1126/science.1131262
Narayanaswamy, R., Levy, M., Tsechansky, M., Stovall, G. M., O'Connell, J. D., Mirrielees, J., et al. (2009). Widespread reorganization of metabolic enzymes into reversible assemblies upon nutrient starvation. Proc. Natl. Acad. Sci. U. S. A. 106, 10147–10152. doi: 10.1073/pnas.0812771106
Norry, F. M., and Loeschcke, V. (2003). Heat-induced expression of a molecular chaperone decreases by selecting for long-lived individuals. Exp. Gerontol. 38, 673–681. doi: 10.1016/s0531-5565(03)00057-3
Pechmann, S., Chartron, J. W., and Frydman, J. (2014). Local slowdown of translation by nonoptimal codons promotes nascent-chain recognition by SRP in vivo. Nat. Struct. Mol. Biol. 21, 1100–1105. doi: 10.1038/nsmb.2919
Pechmann, S., and Frydman, J. (2014). Interplay between chaperones and protein disorder promotes the evolution of protein networks. PLoS Comput. Biol. 10:e1003674. doi: 10.1371/journal.pcbi.1003674
Petersen, J., and Russell, P. (2016). Growth and the environment of Schizosaccharomyces pombe. Cold Spring Harb. Protoc. 2016, 210–226. doi: 10.1101/pdb.top079764
Quax, T. E., Claassens, N. J., Söll, D., and van der Oost, J. (2015). Codon bias as a means to fine-tune gene expression. Mol. Cell 59, 149–161. doi: 10.1016/j.molcel.2015.05.035
Rand, J. D., and Grant, C. M. (2006). The thioredoxin system protects ribosomes against stress-induced aggregation. Mol. Biol. Cell 17, 387–401. doi: 10.1091/mbc.e05-06-0520
Rauscher, R., and Ignatova, Z. (2018). Timing during translation matters: synonymous mutations in human pathologies influence protein folding and function. Biochem. Soc. Trans. 46, 937–944. doi: 10.1042/BST20170422
Sauna, Z. E., and Kimchi-Sarfaty, C. (2011). Understanding the contribution of synonymous mutations to human disease. Nat. Rev. Genet. 12, 683–691. doi: 10.1038/nrg3051
Sharp, P. M., and Li, W. H. (1987). The codon adaptation index--a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Res. 15, 1281–1295. doi: 10.1093/nar/15.3.1281
Soti, C., and Csermely, P. (2007). Aging cellular networks: chaperones as major participants. Exp. Gerontol. 42, 113–119. doi: 10.1016/j.exger.2006.05.017
Subramaniam, A. R., Pan, T., and Cluzel, P. (2013). Environmental perturbations lift the degeneracy of the genetic code to regulate protein levels in bacteria. Proc. Natl. Acad. Sci. U. S. A. 110, 2419–2424. doi: 10.1073/pnas.1211077110
Sun, S., Liang, X., Zhang, X., Liu, T., Shi, Q., Song, Y., et al. (2015). Phosphoglycerate kinase-1 is a predictor of poor survival and a novel prognostic biomarker of chemoresistance to paclitaxel treatment in breast cancer. Br. J. Cancer 112, 1332–1339. doi: 10.1038/bjc.2015.114
Takeda, K., Mori, A., and Yanagida, M. (2011). Identification of genes affecting the toxicity of anti-cancer drug bortezomib by genome-wide screening in S. pombe. PLoS One 6:e22021. doi: 10.1371/journal.pone.0022021
Thommen, M., Holtkamp, W., and Rodnina, M. V. (2017). Co-translational protein folding: progress and methods. Curr. Opin. Struct. Biol. 42, 83–89. doi: 10.1016/j.sbi.2016.11.020
Waldman, Y. Y., Tuller, T., Shlomi, T., Sharan, R., and Ruppin, E. (2010). Translation efficiency in humans: tissue specificity, global optimization and differences between developmental stages. Nucleic Acids Res. 38, 2964–2974. doi: 10.1093/nar/gkq009
Walsh, I. M., Bowman, M. A., Soto Santarriaga, I. F., Rodriguez, A., and Clark, P. L. (2020). Synonymous codon substitutions perturb cotranslational protein folding in vivo and impair cell fitness. Proc. Natl. Acad. Sci. U. S. A. 117, 3528–3534. doi: 10.1073/pnas.1907126117
Yu, C. H., Dang, Y., Zhou, Z., Wu, C., Zhao, F., Sachs, M. S., et al. (2015). Codon usage influences the local rate of translation elongation to regulate co-translational protein folding. Mol. Cell 59, 744–754. doi: 10.1016/j.molcel.2015.07.018
Zemva, J., Fink, C. A., Fleming, T. H., Schmidt, L., Loft, A., Herzig, S., et al. (2017). Hormesis enables cells to handle accumulating toxic metabolites during increased energy flux. Redox Biol. 13, 674–686. doi: 10.1016/j.redox.2017.08.007
Zhang, G., and Ignatova, Z. (2009). Generic algorithm to predict the speed of translational elongation: implications for protein biogenesis. PLoS One 4:e5036. doi: 10.1371/journal.pone.0005036
Zhang, Z., Zhou, L., Hu, L., Zhu, Y., Xu, H., Liu, Y., et al. (2010). Nonsense-mediated decay targets have multiple sequence-related features that can inhibit translation. Mol. Syst. Biol. 6:442. doi: 10.1038/msb.2010.101
Zhou, Z., Dang, Y., Zhou, M., Li, L., Yu, C. H., Fu, J., et al. (2016). Codon usage is an important determinant of gene expression levels largely through its effects on transcription. Proc. Natl. Acad. Sci. U. S. A. 113, E6117–E6125. doi: 10.1073/pnas.1606724113
Zhou, M., Wang, T., Fu, J., Xiao, G., and Liu, Y. (2015). Nonoptimal codon usage influences protein structure in intrinsically disordered regions. Mol. Microbiol. 97, 974–987. doi: 10.1111/mmi.13079
Keywords: translation, codon usage bias, synonymous mutation, protein aggregation, proteostasis, stress tolerance
Citation: Moreira-Ramos S, Arias L, Flores R, Katz A, Levicán G and Orellana O (2023) Synonymous mutations in the phosphoglycerate kinase 1 gene induce an altered response to protein misfolding in Schizosaccharomyces pombe. Front. Microbiol. 13:1074741. doi: 10.3389/fmicb.2022.1074741
Edited by:
Haike Antelmann, Freie Universität Berlin, GermanyReviewed by:
Anna Glushakova, Lomonosov Moscow State University, RussiaFabian M. Commichau, University of Hohenheim, Germany
Copyright © 2023 Moreira-Ramos, Arias, Flores, Katz, Levicán and Orellana. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Omar Orellana, oorellan@uchile.cl