Synthetic control of a fitness tradeoff in yeast nitrogen metabolism
© Bayer et al; licensee BioMed Central Ltd. 2009
Received: 18 September 2008
Accepted: 02 January 2009
Published: 02 January 2009
Skip to main content
© Bayer et al; licensee BioMed Central Ltd. 2009
Received: 18 September 2008
Accepted: 02 January 2009
Published: 02 January 2009
Microbial communities are involved in many processes relevant to industrial and medical biotechnology, such as the formation of biofilms, lignocellulosic degradation, and hydrogen production. The manipulation of synthetic and natural microbial communities and their underlying ecological parameters, such as fitness, evolvability, and variation, is an increasingly important area of research for synthetic biology.
Here, we explored how synthetic control of an endogenous circuit can be used to regulate a tradeoff between fitness in resource abundant and resource limited environments in a population of Saccharomyces cerevisiae. We found that noise in the expression of a key enzyme in ammonia assimilation, Gdh1p, mediated a tradeoff between growth in low nitrogen environments and stress resistance in high ammonia environments. We implemented synthetic control of an endogenous Gdh1p regulatory network to construct an engineered strain in which the fitness of the population was tunable in response to an exogenously-added small molecule across a range of ammonia environments.
The ability to tune fitness and biological tradeoffs will be important components of future efforts to engineer microbial communities.
Many natural and man-made processes, such as lignocellulose digestion , wastewater treatment , environmental remediation , and biofilm formation  are mediated by consortia of microbes rather than a single organism. Often microbial consortia are composed of specialist strains that carry out individual metabolic reactions that benefit multiple community members, increase overall biochemical efficiency and buffer the community from environmental changes. In a recent example, a process involving two metabolic specialist strains of Escherichia coli was observed to efficiently convert xylose and glucose mixtures into fermentation products  more quickly than using a single generalist organism and adapted to changing concentrations of the two sugars by changing the relative abundance of each organism. The manipulation of existing microbial communities and the construction of synthetic communities will be increasingly important for engineering complex biological functions [6, 7].
Synthetic biologists are beginning to design microbial consortia using bacterial quorum sensing. A recent study demonstrated how communicating populations of E. coli can act as an AND gate, exhibiting a gene expression response only when both populations are present . Synthetic ecologies have been constructed where two types of bacteria act as 'predator' and 'prey', with each sensing the other by quorum sensing and dependent on the other for growth . The ability to rationally engineer fitness in a given environment will advance the widespread use of microbial communities for performing biotechnologically-relevant processes.
Underlying the design of microbial consortia is the understanding of and ability to control fitness, communication, and ecological strategies [10, 11] in a single organism. These tools and concepts can then be used to construct synthetic ecologies of interacting microbes useful in downstream engineering applications. Tradeoffs between fitness in different environments are well known in the ecology and engineering literature [12–15], and recent work has recapitulated such an ecological tradeoff by modulating the noise in expression of an antibiotic resistance gene . Yeast populations driving the expression of an antibiotic resistance gene from noisy promoters were better able to survive challenge by antibiotics, but exhibited a fitness disadvantage in media without antibiotic. Such tradeoffs between stress resistance and fitness in stable environments have been observed in many classic ecological studies ranging from prokaryotes to metazoa to plants [10, 13, 17]. The ability to tune the performance of a population of cells for a given environment, including growth, adaptability, and stress resistance, will be useful in future engineering efforts.
Phenotypic variation between members of a population has been shown to be an important parameter in determining how organisms respond to biotic or abiotic environmental challenges . Recent work has highlighted the prevalence of biological variability due to the fundamental limits of deterministic behavior at the cellular level [19–22]. Variability, or noise in gene expression, is a ubiquitous feature of the natural world and has been demonstrated to arise from the small number of molecules involved in cellular processes such as the levels of transcription factors, polymerases, and ribosomes . Noise has been shown to be critical in several biological processes, including determination of competence in Bacillus subtilis [23, 24], eye color vision development in Drosophila melanogaster , and viral latency in bacteriophages  and human pathogens .
Here, we examined the integration of synthetic regulation strategies with endogenous genetic networks to control 'ecological' parameters in a microbial population. Endogenous genetic networks dictate important ecological parameters, including fitness, phenotypic diversity , evolvability , and stress response, such that our ability to rationally manipulate these networks will be important to controlling more complex population- and consortia-level functions. We demonstrated that engineered strains differing only in the expression variability of an enzyme required for metabolizing ammonia, Gdh1p , displayed differences in fitness under 'normal' and stressful ammonia environments . A strain exhibiting Gdh1p expression variability greater than the wildtype strain demonstrated increased resistance to ammonia stress, but lower fitness than wildtype at normal ammonia concentrations. A strain exhibiting lower variability in Gdh1p expression than wildtype displayed the opposite fitness trends – lower than wildtype resistance to ammonia stress and similar fitness to wildtype under normal ammonia concentrations. Finally, we constructed an engineered strain in which the fitness tradeoff was controlled by exogenous addition of a small molecule by placing the endogenous Gdh1p regulator, Dal80p , under the transcriptional control of a galactose-titratable promoter system . Our results suggest that synthetic control of such fitness tradeoffs could be exploited to construct microbial populations and consortia with defined ecological behaviors.
The toxic effects of ammonia have been well documented in animals and plants ; however, recent work by Hess et al demonstrated that yeast are also susceptible to ammonia toxicity . Yeast possess mechanisms to excrete excess nitrogen in the form of amino acids, a rudimentary form of ammonia detoxification analogous to urea in mammals. The authors also demonstrated that ammonia toxicity is increased under low potassium ion levels, potentially due to the ability of ammonia ions to enter the cell via potassium channels. Therefore, while ammonia is an essential source of cellular nitrogen, environmental conditions, such as excess ammonia or low potassium, can be toxic to the cells. This paradoxical nature of ammonia (essential but toxic) provides a promising experimental system to examine the synthetic regulation of a fitness tradeoff between resource abundant and resource limited environments.
We used a previously reported direct competition assay  to examine the population fitness of the GDH1 promoter mutant clones under different ammonia concentrations. The relative contributions of noise and abundance in Gdh1p levels to the fitness tradeoff were determined by examining the fitness of sets of strains exhibiting constant abundance but different noise levels and constant noise but different abundance levels. To represent the relative contributions of fitness and stress resistance, we measured fitness via direct competition with a reference strain at both 17 and 556 mM ammonia. We defined a fitness term, Wenv, as the ratio of fitness at the higher ammonia concentration to fitness at the lower ammonia concentration, with Wenv for a wildtype strain containing the endogenous Gdh1p promoter equal to 1. Thus, clones with Wenv values greater than 1 exhibit enhanced growth at high ammonia versus low ammonia concentrations, while clones with Wenv values lower than 1 exhibit higher growth at low ammonia than high ammonia concentrations. Constant abundance-varying noise mutant sets with low, medium, and high abundance Gdh1p expression exhibited variations in Wenv that correlated with noise values (Fig. 2c). In contrast, Wenv did not show any correlation with abundance values under constant noise values (Fig. 2d). Our results demonstrate a correlation between Gdh1p noise and a fitness tradeoff across varying ammonia environments.
As previously indicated, ammonia toxicity is enhanced in environments with low amounts of potassium. To further examine the correlation between Gdh1p noise and stress resistance, we challenged the high and low noise strains in media with 600 mM ammonia and low potassium (17 mM), conditions at which ammonia is toxic to the cells. The high noise strain exhibited increases in CFUs at 30 minutes after ammonia challenge, but a decrease in viable cells at 60 minutes (Fig. 3c). In contrast, the low noise strain exhibited a slight loss of viable cells at 30 minutes after ammonia challenge that becomes more pronounced at 60 minutes. While it appeared that both strains were susceptible to the toxic effects of ammonia in low potassium media, the high noise strain displayed a significant lag period before the number of viable cells decreased. Such temporal differences in resistance may be critical in uncertain or fluctuating environments, where a high noise strain may be better able to buffer deleterious fluctuations in ammonia concentration.
We also examined the effects of noise in Gdh1p expression on fitness under 'physiological' ammonia environments. Yeast are commonly grown in 40 mM ammonia in the laboratory, while ammonia concentrations in the natural yeast habitat may be lower . Therefore, we measured population growth of the low and high noise strains across a range of low ammonia concentrations. The growth assays that examine increases in CFUs over time exhibited too much experimental error at these ammonia concentrations to reliably measure growth differences between the high and low noise strains and the wildtype strain. Therefore, we used the direct competition growth assay to measure fitness at multiple ammonia concentrations. The low noise strain exhibited similar fitness to wildtype at each assayed ammonia concentration (Fig. 3d). In contrast, the high noise strain exhibited significantly lower fitness than the wildtype and low noise strains at 10 mM and 20 mM ammonia, and similar fitness at 40 mM ammonia. Our results indicate that while the high noise strain exhibited increased resistance to ammonia stress, it also exhibited decreased fitness under lower ammonia concentrations, supporting the role of Gdh1p noise in affecting a tradeoff between stress resistance at high ammonia concentrations and fitness at low ammonia concentrations.
We measured the effects of varying Dal80p expression levels on the Gdh1p expression profile through flow cytometry analysis. Mean Gdh1p abundance did not change with varying Dal80p levels (Fig. 4c), indicating that combinatorial interactions at the GDH1 promoter or a complex regulatory network may be associated with this nitrogen assimilation pathway. However, the noise in Gdh1p expression decreased with increasing Dal80p levels (Fig. 4d). Specifically, in our engineered strain, low levels of Dal80p (under the absence of galactose) resulted in ~15% higher noise in Gdh1p expression than the wildtype strain, whereas high levels of Dal80p (in the presence of ~1% galactose) reduced noise to ~80% of the wildtype strain.
Variation between organisms in a population has long been recognized as an important parameter in predicting evolutionary dynamics. The past decade of research on the stochastic nature of gene expression has further highlighted the importance of variation on the functions of biological systems. Using a model of the interplay between variability and fitness tradeoffs , we discovered a similar tradeoff in yeast ammonia metabolism and achieved synthetic control of this tradeoff by manipulating noise through an endogenous regulatory circuit. We demonstrated that the level of noise in Gdh1p expression dictated the relative balance between resistance to toxic levels of ammonia and fitness in lower levels of ammonia. Furthermore, by examining the endogenous regulation of Gdh1p we discovered a convenient point in the circuit to regulate noise, and thus bring the fitness tradeoff under tunable control.
While our current work does not demonstrate a mechanism for stress resistance and fitness effects, similar studies in the literature may highlight routes for future investigation. Of particular interest is the mechanism by which increased noise in Gdh1p expression confers enhanced resistance to ammonia toxicity. One hypothesis is that in populations with larger Gdh1p distributions, the subset of the population with high Gdh1p expression is able to tolerate more ammonia by excreting more amino acids as described by Hess et al . In this situation, at least a subset of the population would be able to survive temporary excesses of ammonia. The tradeoff between stress resistance and fitness at lower ammonia concentrations may be due to the energetic cost and deleterious effects of large fluctuations in protein expression, as has been observed in a previous bacterial study . Such a mechanism may also explain why fitness tradeoffs were not observed with other nitrogen sources.
As potassium concentrations in the endogenous yeast habitat are likely lower than in laboratory conditions , the resistance to ammonia toxicity may be a significant contributor to survival of individuals and populations. Whether natural populations have taken advantage of modulating noise in enzyme expression in response to environments with excess ammonia, rather than manipulating amino acid excretion or other mechanisms, is an open and intriguing question, and is supported here by the observation that modulation of an endogenous transcriptional regulator modulates enzyme expression noise. Tuning adaptation to ammonia toxicity by synthetic manipulation of Gdh1p noise would result in populations that have low fitness in stable environments, but an enhanced ability to survive stressful periods. This phenotype is reminiscent of bacterial persistence, which has been shown to be driven by noise in gene expression . As we continue to link cellular processes with ecological parameters we will gain new insight into evolutionary and ecological processes such as adaptation, variation, and evolvability. For example, recent computational studies have predicted that the design of regulatory networks are determined by the fitness benefits of regulating noise in a population . In addition, the development of tools and concepts for manipulating ecological parameters will allow engineers to begin to more effectively build microbial consortia for potential applications in environmental remediation, energy, and therapeutics.
All manipulations were performed with the S288c background. Yeast were grown in synthetic complete media with the nitrogen source as specified. Primer sequences are provided in Additional file 1.
To construct mutant libraries of the GDH1 promoter, primers flanking 500 nucleotides upstream of the GDH1 coding region (1043500 – 1043050, chromosome XV) were used to amplify the fragment from yeast genomic DNA using KOD polymerase (Novagen) according to the manufacturer's instructions. The fragment was then diluted into mutagenic PCR buffer  (7 mM MgCl2, 0.5 mM MnCl2, 50 mM KCl, 10 mM Tris pH 8.3 with 1 mM dGTP, 0.2 mM dCTP, 0.2 mM dTTP, 0.2 mM dATP) and further amplified using Taq polymerase (Roche) according to the manufacturer's instructions. Separately, the leucine biosynthesis gene (LEU2) was amplified from pRS315  using KOD polymerase. The LEU2 gene fragment and the promoter library were then ethanol precipitated, resuspended in water, and PCR assembled together through overlapping primer sequences. The resulting large fragment was then transformed into yeast strains using a standard lithium acetate procedure . Transformants were selected in liquid synthetic complete leucine dropout media. The resulting library was grown to stationary phase and frozen in 15% glycerol at -80°C.
Integration of the GAL promoter was performed by amplifying the GAL1-10 promoter sequence from pRS314-Gal  using KOD polymerase. This fragment was PCR assembled with the leucine biosynthesis gene (LEU2) from pRS315  along with flanking homologous regions to the DAL80 upstream region (506000 – 504030 and 506500 – 506530 on chromosome XI) using KOD polymerase. The construct was transformed into yeast using a standard lithium acetate procedure and colonies were selected on synthetic complete leucine dropout agar plates. Integration was confirmed by colony PCR with primers flanking and internal to the integrated construct. Yeast DNA extraction was performed as previously described using the 'bust n' grab' method .
Cells were pelleted and frozen in liquid nitrogen. Pellets were resuspended in a 50 mM NaOAc (pH 5.2), 10 mM EDTA buffer. Cells were lysed by the addition of SDS to a final concentration of 1.6% and an equal volume of acid phenol. Solutions were kept at 65°C with intermittent vortexing for 10 min. After cooling on ice, the aqueous phase was extracted and further extraction was carried out with an equal volume of chloroform. RNA was further isolated and concentrated by use of RNeasy columns (Qiagen) according to manufacturer's instructions. Total RNA was quantified by OD260 readings. RNA samples were treated with DNase (Invitrogen) according to manufacturer's instructions. cDNA was synthesized using gene-specific primers and Superscript III reverse transcriptase (Invitrogen) according to the manufacturer's instructions. qRT-PCR was carried out on this cDNA using an iCycler iQ system (BioRAD). Samples were prepared using the iQ SYBR green supermix and primer pairs specific for different templates. Data were analyzed using the iCycler iQ software.
Yeast cells growing in exponential phase (OD600~0.5) were spun down, washed with sterile water, and resuspended in synthetic complete media with 600 mM ammonia. An aliquot of this culture was serially diluted (10-fold dilutions) and 50 μL of the dilutions were plated on YPD agar media. Cells were grown on the solid media for 2 days at 30°C and colonies were counted to measure the change in colony forming units (CFUs) over time.
Fitness was assayed by direct competition versus a common reference strain . The competitor and reference strain constitutively express different fluorescent proteins (GFP and CFP, respectively) from the ADH1 promoter integrated into the chromosome. The frequency of competitor and reference strains were quantified before and after the growth period by counting the numbers of GFP expressing cells to non-GFP expressing cells. Fitness (w) of the competitor strain is reported as the natural log of the change in frequency of the strain during the competitive growth period versus the change in frequency of the reference strain over the same growth period:
w = ln (Δfrequency of competitor strain/Δfrequency of reference strain)
All competitor strains were derivatives of the S288C background, while the reference strain was derived from the W303 background. These strains showed different electronic volume versus side-scatter distributions that can also be used to quantify population numbers, in good agreement with the values obtained from fluorescent measurements. Equal amounts of competitor and reference strain were mixed and grown in indicated liquid media for 3 generations (approximately 6 hours). The frequency of competitor and reference strain were quantified before and after the growth period by counting the numbers of GFP expressing cells to non-GFP expressing cells by flow cytometry using a Quanta SC flow cytometer (Beckman Coulter) equipped with the MPL system. Samples were excited with a 488 nm laser and GFP fluorescence was detected with a 525 nm bandpass filter. A gate was set above the non-GFP expressing cells in the Quanta analysis software to partition fluorescent from non-fluorescent cells. Samples of only reference or competitor strains and serial dilutions of ratios of competitor to reference strains were run in parallel as controls. 5,000 events were collected per sample.
Two gates were used to standardize each cell population for analysis using 'magnetic gating' in FlowJo flow cytometry analysis software (Tree Star, Inc.). The first gate isolated cells displaying regular morphology based on electronic volume and side-scatter, while the second gate removed non-fluorescent cells from the distribution. This gating method was compared against other methods previously described and the abundance and noise trends observed were consistent between methods [19, 34]. Noise was calculated as the square of the coefficient of variation (σ2/p2) of the distribution . Abundance was calculated as the mean of the distribution. 50,000 events were analyzed to calculate noise for each sample. Noise trends were similar when calculated as the coefficient of variation (σ/p) and the variance (σ2).
We thank D. Endy, M. Elowitz, J. Silberg, and J. Tabor for critical reading and comments on this manuscript. This work was supported by the Caltech Grubstake Program (grant to CDS), the National Institutes of Health (grant to CDS; training grant TSB; fellowship to KGH), and the Department of Defense (fellowship to CLB).
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.