Rapid evolutionary adaptation to growth on an ‘unfamiliar’ carbon source
BMC Genomics volume 17, Article number: 674 (2016)
Cells constantly adapt to changes in their environment. When environment shifts between conditions that were previously encountered during the course of evolution, evolutionary-programmed responses are possible. Cells, however, may also encounter a new environment to which a novel response is required. To characterize the first steps in adaptation to a novel condition, we studied budding yeast growth on xylulose, a sugar that is very rarely found in the wild.
We previously reported that growth on xylulose induces the expression of amino acid biosynthesis genes in multiple natural yeast isolates. This induction occurs despite the presence of amino acids in the growth medium and is a unique response to xylulose, not triggered by naturally available carbon sources. Propagating these strains for ~300 generations on xylulose significantly improved their growth rate. Notably, the most significant change in gene expression was the loss of amino acid biosynthesis gene induction. Furthermore, the reduction in amino-acid biosynthesis gene expression on xylulose was tightly correlated with the improvement in growth rate, suggesting that internal depletion of amino-acids presented a major bottleneck limiting growth in xylulose.
We discuss the possible implications of our results for explaining how cells maintain the balance between supply and demand of amino acids during growth in evolutionary ‘familiar’ vs. ‘novel’ conditions.
Growing cells need to adjust the rate of biomass accumulation with the cell division time. Tight regulation of protein synthesis rate is particularly critical, as it consumes a large fraction of cellular resources. Indeed, studies over the past decades have shown that expression of translation-related genes, such as those required for the production of ribosomes, are modulated in proportion to cellular growth-rate in rapidly growing cells such as microorganisms or cancer cells [1, 2]. In the budding yeast, for example, expression of genes coding for ribosomal proteins (RP), proteins required for ribosome biogenesis (Ribi) or other translation related factors, strongly correlated with cellular growth rate in various environments, and can in fact be used to predict growth rate [3–7].
Amino acids provide the substrate for protein synthesis and their influx should similarly be coordinated with growth rate. Amino acids are taken up from the environment using specialized transporters. In addition, most amino-acids can be produced de-novo as part of cellular metabolism. Pathways involved in de-novo synthesis are induced at the transcriptional level when intra-cellular amino acids are depleted [8–11]. In particular, depletion of amino acids leads to the accumulation of uncharged tRNAs, leading to increased translation of GCN4, the major transcription factor inducing the expression of dozens of genes required for amino acid biosynthesis [12–16]. This pathway is known to be readily induced under conditions of amino acid limitation. However, expression of this pathway is not modulated in response to environmental changes as ribosome-associated genes are, and does not appear to correlate with growth rate. Specifically, shifting cells between different carbon sources rapidly modulated the expression of RP and Ribi genes, but did not significantly alter the expression of genes required for amino acid biosynthesis. This indicates that the capacity of translation machinery (ribosomes), rather than amino acid de-novo production, is adjusted according to the metabolic capacity upon varying conditions, thereby ensuring optimal protein production.
We recently reported a notable exception to this behavior while studying the transcription profile of different budding yeast strains upon shift from glucose to xylulose, the only pentose sugar budding yeast are able to metabolize . Xylulose is phosphorylated to xylulose-5-phosphate, an intermediate metabolite of the pentose phosphate pathway and thus a substrate for cellular metabolism. However, xylulose is not naturally available in significant amounts outside of living cells [18, 19] and can therefore be viewed as an ‘unfamiliar’ or ‘novel’ carbon source for which yeast did not undergo evolutionary adaptation. Inter-strain diversity in growth rate and ethanol production was significantly higher during growth on xylulose than that observed during growth on other ‘familiar’ carbon sources such as glucose and galactose. While this large diversity could be a sign of local adaptation , it could also reflect the rarity by which this carbon source is encountered in nature.
Measuring expression profiles of twenty four budding yeast strains growing on xylulose we noted that no correlation was observed between growth rate and expression of the gene groups that were shown to be tightly correlated with growth rate during growth on different carbon sources (e.g. RP and Ribi). Instead, xylulose strongly induced expression of amino acid biosynthesis genes. Furthermore, the level of induction was significantly (anti)correlated with growth rate on xylulose. Searching for additional conditions that activate this gene group in a growth rate-correlated manner, we identified a large number of gene deletion mutants which may similarly represent conditions not optimized by evolution . Based on this observation, we hypothesized that familiarizing cells to such ‘novel’ conditions by evolutionary adaptation might restore the typical expression profile.
To test this prediction, we evolved twelve wild-type budding yeast strains on xylulose for ~300 generations and compared their gene expression profiles before and after this evolution process. As predicted, we found that induction of amino acid biosynthesis genes was practically eliminated in the evolutionarily-optimized cells. In addition, the degree of growth improvement was correlated with the change in expression of this gene group. We discuss the implications of our results for explaining how cells maintain the balance between amino acid supply and demand for optimizing protein synthesis and growth rate in ‘familiar’ vs. ‘novel’ conditions.
Microevolution improved growth on xylulose of multiple strains
We previously reported the growth rate and gene expression of a collection of twelve wild type S. cerevisiae strains, isolated from diverse geographical origins and natural habitats, on multiple carbon sources including xylulose . All of these strains could grow on xylulose, although their growth rates and fermentation yields varied greatly. With the aim of optimizing growth on xylulose, we subjected these twelve strains to microevolution by serial dilutions on media containing xylulose as the sole carbon source (Fig. 1a). This process was carried out for a total of ~300 generations.
Measuring the growth curves of the evolved strains revealed that the strains underwent different types of growth improvements (Fig. 1b). For example, while some strains improved growth rate during logarithmic growth, others improved the saturation density, improving the overall yield. In attempt to describe the effect of evolution using a single metric, we first considered separately the logarithmic growth rate and the final OD level reached in saturation (Table 1). However, the evolutionary effect with respect to these two parameters was not consistent across all strains and we therefore examined for a different metric.
To capture the different aspects of growth improvement using a single parameter, we defined an effective growth rate as the maximal OD level reached during growth divided by the time elapsed until saturation (following inoculation in xylulose media). By this measure, eleven of the twelve strains exhibited improved growth, while one maintained the same effective growth rate as its parent strain (Fig. 1c and Table 1).
Evolved strains eliminated the AA-biosynthesis gene induction upon shift from glucose to xylulose
As described above, we previously noted that growth on xylulose induces expression of AA-biosynthesis genes. We asked whether this induction is maintained also in our evolutionarily-optimized strains. Using RNA-Seq, we characterized the transcription profiles of all strains in our collection, before and after evolution, during logarithmic growth on xylulose and on glucose. This was performed in two biological repeats, totaling ninety six genome-wide expression measurements. For each strain (before and after evolution), we defined the (log) ratio of gene expression in xylulose vs. glucose. We then systematically defined the groups of genes whose level of expression significantly changed upon shift from glucose to xylulose. This was done by first identifying all the genes whose level of expression changed significantly and then examining for functional enrichments by projecting this set of varying genes on pre-defined functional gene groups . Thus, groups of genes were identified which were either induced or reduced in expression on xylulose compared to glucose, in pre-evolved wild type strains and in the evolved populations.
Consistent with our prediction, the most significant difference in the pattern of induced expression before and after evolution was that of the AA-biosynthesis gene group (Fig. 2a, b). These genes were significantly induced in the ancestral strains (prior to evolution), but this induction was practically absent following evolution. This result was independent of how the group was pre-defined, either as a co-regulated module or by the targets of GCN4, the AA-biosynthesis master regulator (Additional file 1). The observed change in expression pattern was significantly larger than the inter-strain variation and was observed exclusively during growth on xylulose, but not on glucose (Additional file 2).
We further characterized the expression changes in gene groups associated with ribosome production (RP and Ribi), protein synthesis and the environmental stress response (ESR). Expression of these gene groups typically correlates with cellular growth rate. Expression of RP, Ribi and protein synthesis genes was up-regulated in xylulose compared to glucose following evolution (Fig. 2b). In addition, genes induced as part of the known environmental stress response [6, 23] were down-regulated in the evolved populations, while genes which are normally reduced in expression as a response to environmental stress were up-regulated. These results are consistent with the fact that the expression of these gene groups is correlated with growth rate and demonstrate adaptation at the level of gene expression. Here too, the observed effect resulted from changes in expression during growth on xylulose, rather than on glucose (Additional file 2).
Interestingly, in evolved strains, genes whose expression is typically induced in the presence of environmental stress were down-regulated on xylulose, compared to glucose (Fig. 2b). This down-regulation resulted from both reduced expression of these genes on xylulose and also increased expression during growth on glucose (Additional file 2). This may indicate a tradeoff in the efficiency by which glucose and xylulose are utilized, whereby adaptation to growth on xylulose was acquired, to some extent, at the expense of growth on glucose.
To determine whether the above changes in gene expression profiles are reflected in expression of regulatory elements, we examined for expression changes among a large set of gene groups characterized by the transcription factors associated with each group (n = 69). Notably, the most significant changes in expression on xylulose following evolution were the down-regulation of GCN4 and up-regulation of FHL1 and SFP1 (Fig. 2c). GCN4 is the master regulator of amino acid biosynthesis genes and its down-regulation following evolution is coherent with reduction in expression of this group of genes. Both FHL1 and SFP1 are known to be involved in the regulation of RP gene expression. SFP1 is also associated with Ribi gene expression regulation.
The reduction in AA-biosynthesis gene expression correlates with effective growth rate improvement
The expression changes in the aforementioned gene groups as a result of evolutionary adaptation to growth on xylulose suggest that the extent of the change might be proportional to the degree of adaptation. To test whether this is the case, we examined the correlation, or lack thereof, of the expression changes of these genes groups with the degree of improvement in growth on xylulose, as defined by the ratio between effective growth rates after and before evolution. We found that the change in expression level (evolved/non-evolved) of the amino acid biosynthesis regulator – GCN4 – was anti-correlated with growth improvement (pv = 0.006, Fig. 3). In addition, the expression changes of ribosomal protein genes were (positively) correlated with growth improvement on xylulose (pv = 0.017). This was independent on how the group was defined; thus expression level change in amino acid biosynthesis genes was significantly anti-correlated with the level of growth improvement on xylulose (pv =0.008) while that of protein synthesis genes was positively correlated (pv = 0.009, Additional file 3). No correlation was observed between growth improvement and expression changes of other gene groups (Additional file 3).
Discussion and conclusions
We previously reported that cells growing on xylulose show a distinct expression pattern characterized by an induction of AA-biosynthesis group in a manner that was inversely correlated with the growth rate of the cell. This induction occurred despite the fact that the media was supplemented with all required amino-acids. Searching through a large dataset , we detected a similar general induction of this group of genes in a group of mutant strains deleted of a variety of genes , most of which are not related to AA metabolism (Fig. 4a). We reasoned that the two situations are related by the fact that both represent conditions that are not evolutionarily optimized: xylulose is not naturally available outside cells, so budding yeast were unlikely to frequently encounter it during the course of evolution. Similarly, deleting any particular gene modifies the internal state of the cell in a way that it did not undergo evolutionary optimization for. This indicated that induction of the AA-biosynthesis gene group may be a consequence of such non-optimized internal state.
In the preset study we examine this hypothesis by tested whether this induction is lost once the cells have had a chance to (evolutionarily) adapt to this new condition. Indeed, not only was the induction lost upon such adaptation, resulting in reduced expression of this group specifically in xylulose (Fig. 4a), but the improvement in effective growth rate was correlated with the extent to which AA-biosynthesis expression was reduced. While these results do not demonstrate a direct causal relationship between evolutionary improvement and the observed changes in gene expression, especially considering that at least two strains display similar improvement but different expression pattern, they strongly suggest that amino-acids biosynthesis forms a major bottleneck for growth in conditions not subject to evolutionary adaptation, as we postulated.
In order to achieve a balanced cellular level of amino acids, the incoming flux of amino acids must match the rate by witch proteins are synthesized. AA-biosynthesis is induced by accumulation of uncharged tRNAs, reflecting a demand for amino-acids that is larger than their supply. By contrast, expression of ribosomal-associated genes is largely dependent on signaling pathways such as TOR and PKA. This implies that in order to maintain balance between supply and demand, a tuning mechanism between these pathways is required, which we postulate is achieved primarily through evolutionary adaptation.
We propose a model comprised of two possible scenarios relating to this process (Fig. 4b). In an evolutionary ‘familiar’ environment, evolution had optimized signaling pathways to adjust the level of translation machinery with the expected rate of amino acid supply. However, in the absence of such evolutionary tuning, as may be exemplified by growth on xylulose or in the presence of genetic mutations, a transient imbalance occurs which triggers an internal correction mechanism activated through the induction of the GCN4-module. Subsequent evolutionary adaptation to the ‘unfamiliar’ environment will improve the tuning mechanism, thereby reducing GCN4-dependent gene expression.
The need for internal adjustment compromises cellular growth. Thus, according to this model and confirmed by our results, the first step in evolutionary adaptation to ‘unfamiliar’ conditions would be to optimize the balance between supply and demand of amino acids by adjusting the level of translation machinery, thereby abrogating the need for the internal ‘backup’ control of amino-acid biosynthesis adjustments.
Biotechnological applications often rely on engineering cells to either utilize new substrates or produce novel products. For example, a major limitation of using budding yeast as an economically viable solution for generating biofuel is its inability to ferment xylose, motivating engineering efforts for establishing a xylose-utilizing budding yeast strain. These engineered strains may experience a similar situation to the one we describe above, in that they are put in a condition not optimized by evolution. Our results suggest that co-optimizing amino-acid metabolism may be useful for improving growth performance of these strains.
Strains and media
A complete description of wild type yeast strains used in this study can be found in . YP-xylulose medium was prepared as described in . Briefly, D-xylose was enzymatically converted to D-xylulose using xylose isomerase (Sigma-Aldrich). This yielded a 70:30 mixture of xylose:xylulose, as measured by HPLC (Agilent 1200 series). YP media was supplemented with the pentose mixture to a final concentration of 2 g/L xylulose. Although growth on xylose was previously reported for certain wild type S. cerevisiae strains , in our experiments xylose was not utilized by any of the strains during the growth period, as determined by HPLC analysis of the growth media after reaching saturation.
Evolution and growth measurements
All laboratory evolution experiments were carried out by serial dilution in 24-well plates with breathable seals. Plates were grown on an orbital shaker at 600 rpm in 30 °C. Each evolving strain was grown on 1.2 ml of YP-xylulose medium until reaching stationary phase and then diluted daily by a factor of 1:64 into fresh media (6 generations per dilution). Growth of wild type and evolved populations was characterized using a Hamilton robotic system. Cultures were grown in 96-well plates with 600 rpm shaking in 30 °C. Mineral oil was added to each well to prevent evaporation. OD600 measurements were taken during growth using a Tecan F500 plate reader.
RNA sequencing and normalization
Total mRNA of each of the twelve strains from our wild type yeast strain collection, both before and after evolution, was sequenced, during growth on glucose and during growth on xylulose. Two repeats were performed for each strain and condition, thus totaling 96 expression measurements. Prior to RNA extraction, populations of each strain before and after evolution were grown on YP media supplemented with either glucose or xylulose. Cells were harvested at mid log phase and frozen immediately in liquid nitrogen. For total RNA extraction, cell lysis was performed in a 96 deep-well plate by adding 450 μl of lysis buffer containing 1 M sorbitol (Sigma-Aldrich), 100 mM EDTA and 0.45 μl lyticase (10 IU/μl). The plate was then incubated in 30 °C for 30′ and centrifuged for 10′ at 2500 rpm after which the supernatant was removed. RNA extraction proceeded by following the NucleoSpin® 96 RNA kit (Macherey-Nagel) protocol, substituting β-mercaptoethanol with DTT. RNA extracts were subsequently fragmented and selected for polyadenylation using PolydT beads. Barcoded cDNA was prepared using Smartscribe and unique molecular identifier (UMI)-containing primers , and subjected to RNase H treatment. A second adaptor was ligated, and cDNA was amplified by PCR using KAPA HiFi DNA polymerase for 15 cycles. cDNA was sequenced using an Illumina HiSeq2000 platform with 50 bps single reads. Reads were mapped to the S. cerevisiae S288c reference genome (SGD R64) using Bowtie. Reads mapping to rRNA were filtered out. Normalization for PCR bias was carried out using UMI scoring . The total number of mapped reads per sample was normalized independently to 106.
The expression value of each gene in each strain, either before or after evolution, was averaged over two repeats. Except for the individual amino-acid genes expression comparison (Fig. 2a), all expression values of gene groups before and after evolution were obtained by averaging the number of mapped reads per gene in the gene group, per strain, and then averaging across all strains either before or after evolution. For comparison of amino-acid biosynthesis genes expression between populations before and after evolution (Fig. 2a), the mean expression level of each gene before and after evolution was first calculated by averaging over all strains and the two mean values were normalized by the total number of reads per gene across all samples.
Environmental stress response
Unique molecular identifier
Airoldi EM, Huttenhower C, Gresham D, Lu C, Caudy AA, Dunham MJ, Broach JR, Botstein D, Troyanskaya OG. Predicting cellular growth from gene expression signatures. PLoS Comput Biol. 2009;5:e1000257.
Castrillo JI, Zeef LA, Hoyle DC, Zhang N, Hayes A, Gardner DCJ, Cornell MJ, Petty J, Hakes L, Wardleworth L, Rash B, Brown M, Dunn WB, Broadhurst D, O’Donoghue K, Hester SS, Dunkley TPJ, Hart SR, Swainston N, Li P, Gaskell SJ, Paton NW, Lilley KS, Kell DB, Oliver SG. Growth control of the eukaryote cell: a systems biology study in yeast. J Biol. 2007;6:4.
Brauer MJ, Huttenhower C, Airoldi EM, Rosenstein R, Matese JC, Gresham D, Boer VM, Troyanskaya OG, Botstein D. Coordination of growth rate, cell cycle, stress response, and metabolic activity in yeast. Mol Biol Cell. 2008;19:352–67.
Regenberg B, Grotkjaer T, Winther O, Fausbøll A, Akesson M, Bro C, Hansen LK, Brunak S, Nielsen J. Growth-rate regulated genes have profound impact on interpretation of transcriptome profiling in Saccharomyces cerevisiae. Genome Biol. 2006;7:R107.
Berry DB, Gasch AP. Stress-activated genomic expression changes serve a preparative role for impending stress in yeast. Mol Biol Cell. 2008;19:4580–7.
Gasch AP, Spellman PT, Kao CM, Carmel-Harel O, Eisen MB, Storz G, Botstein D, Brown PO. Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2000;11:4241–57.
Lai L-C, Kissinger MT, Burke PV, Kwast KE. Comparison of the transcriptomic “stress response” evoked by antimycin A and oxygen deprivation in Saccharomyces cerevisiae. BMC Genomics. 2008;9:627.
Schürch A, Miozzari J, Hütter R. Regulation of tryptophan biosynthesis in Saccharomyces cerevisiae: mode of action of 5-methyl-tryptophan and 5-methyl-tryptophan-sensitive mutants. J Bacteriol. 1974;117:1131–40.
Wolfner M, Yep D, Messenguy F, Fink GR. Integration of amino acid biosynthesis into the cell cycle of Saccharomyces cerevisiae. J Mol Biol. 1975;96:273–90.
Jia MH, Larossa RA, Lee JM, Rafalski A, Derose E, Gonye G, Xue Z. Global expression profiling of yeast treated with an inhibitor of amino acid biosynthesis, sulfometuron methyl. Physiol Genomics. 2000;3:83–92.
Natarajan K, Meyer MR, Jackson BM, Slade D, Roberts C, Hinnebusch AG, Marton MJ. Transcriptional profiling shows that Gcn4p is a master regulator of gene expression during amino acid starvation in yeast. Mol Cell Biol. 2001;21:4347–68.
Hinnebusch AG. Translational regulation of GCN4 and the general amino acid control of yeast. Annu Rev Microbiol. 2005;59:407–50.
Wek SA, Zhu S, Wek RC. The histidyl-tRNA synthetase-related sequence in the eIF-2 alpha protein kinase GCN2 interacts with tRNA and is required for activation in response to starvation for different amino acids. Mol Cell Biol. 1995;15:4497–506.
Zhu S, Sobolev AY, Wek RC. Histidyl-tRNA synthetase-related sequences in GCN2 protein kinase regulate in vitro phosphorylation of eIF-2. J Biol Chem. 1996;271:24989–94.
Wek RC, Jackson BM, Hinnebusch AG. Juxtaposition of domains homologous to protein kinases and histidyl-tRNA synthetases in GCN2 protein suggests a mechanism for coupling GCN4 expression to amino acid availability. Proc Natl Acad Sci. 1989;86:4579–83.
Dever TE, Feng L, Wek RC, Cigan AM, Donahue TF, Hinnebusch AG. Phosphorylation of initiation factor 2α by protein kinase GCN2 mediates gene-specific translational control of GCN4 in yeast. Cell. 1992;68:585–96.
Tamari Z, Rosin D, Voichek Y, Barkai N. Coordination of gene expression and growth-rate in natural populations of budding yeast. PLoS One. 2014;9:e88801.
Gunther WR, Wang Y, Ji Y, Michaelis VK, Hunt ST, Griffin RG, Román-Leshkov Y. Sn-Beta zeolites with borate salts catalyse the epimerization of carbohydrates via an intramolecular carbon shift. Nat Commun. 2012;3:1109.
Granström TB, Takata G, Tokuda M, Izumori K. Izumoring. J Biosci Bioeng. 2004;97:89–94.
Warringer J, Zörgö E, Cubillos FA, Zia A, Gjuvsland A, Simpson JT, Forsmark A, Durbin R, Omholt SW, Louis EJ, Liti G, Moses A, Blomberg A. Trait variation in yeast is defined by population history. PLoS Genet. 2011;7:e1002111.
Hughes TR, Marton MJ, Jones AR, Roberts CJ, Stoughton R, Armour CD, Bennett HA, Coffey E, Dai H, He YD, Kidd MJ, King AM, Meyer MR, Slade D, Lum PY, Stepaniants SB, Shoemaker DD, Gachotte D, Chakraburtty K, Simon J, Bard M, Friend SH. Functional discovery via a compendium of expression profiles. Cell. 2000;102:109–26.
Ihmels J, Friedlander G, Bergmann S, Sarig O, Ziv Y, Barkai N. Revealing modular organization in the yeast transcriptional network. Nat Genet. 2002;31:370–7.
Causton HC, Ren B, Koh SS, Harbison CT, Kanin E, Jennings EG, Lee TI, True HL, Lander ES, Young RA. Remodeling of yeast genome expression in response to environmental changes. Mol Biol Cell. 2001;12:323–37.
Olsson L, Lindén T, Hahn-Hägerdal B. A rapid chromatographic method for the production of preparative amounts of xylulose. Enzyme Microb Technol. 1994;16:388–94.
Wenger JW, Schwartz K, Sherlock G. Bulk segregant analysis by high-throughput sequencing reveals a novel xylose utilization gene from Saccharomyces cerevisiae. PLoS Genet. 2010;6:e1000942.
Kivioja T, Vähärautio A, Karlsson K, Bonke M, Enge M, Linnarsson S, Taipale J. Counting absolute numbers of molecules using unique molecular identifiers. Nat Methods. 2012;9:72–4.
We would like to thank Dr. Justin Fay for providing wild type S. cerevisiae strains. This work was supported by the NIH (P50GM068763), the European Research Council, the Israel Science Foundation, Azrieli award for innovative student proposals and the Helen and Martin Kimmel Award for Innovative Investigations. NB is the incumbent of the Lorna Greenberg Scherzer Professorial Chair, and the Minerva Foundation. YP acknowledges support of the Minerva Center for Live Emulation of Genome Evolution in the Lab. ZT is supported by the Israeli Ministry of Science and Technology. AHY is supported by the Human Frontier Science Program.
NIH, European Research Council, Israel Science Foundation, Helen and Martin Kimmel Award for Innovative Investigations, Minerva Center for Live Emulation of Genome Evolution in the Lab, Israeli Ministry of Science and Technology and Azrieli award for innovative student proposals.
Availability of data and materials
Expression data generated during the current study are available in the GEO repository, under accession number GSE76077.
ZT performed RNA-Seq experiments and data analysis. AY performed evolutionary experiments. All authors conceived the study. ZT and NB wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
About this article
Cite this article
Tamari, Z., Yona, A.H., Pilpel, Y. et al. Rapid evolutionary adaptation to growth on an ‘unfamiliar’ carbon source. BMC Genomics 17, 674 (2016). https://doi.org/10.1186/s12864-016-3010-x
- Amino-acid biosynthesis