Computational analysis of GAL pathway pinpoints mechanisms underlying natural variation

Quantitative traits are measurable phenotypes that show continuous variation over a wide phenotypic range. Enormous effort has recently been put into determining the genetic influences on a variety of quantitative traits with mixed success. We identified a quantitative trait in a tractable model system, the GAL pathway in yeast, which controls the uptake and metabolism of the sugar galactose. GAL pathway activation depends both on galactose concentration and on the concentrations of competing, preferred sugars such as glucose. Natural yeast isolates show substantial variation in the behavior of the pathway. All studied yeast strains exhibit bimodal responses relative to external galactose concentration, i.e. a set of galactose concentrations existed at which both GAL-induced and GAL-repressed subpopulations were observed. However, these concentrations differed in different strains. We built a mechanistic model of the GAL pathway and identified parameters that are plausible candidates for capturing the phenotypic features of a set of strains including standard lab strains, natural variants, and mutants. In silico perturbation of these parameters identified variation in the intracellular galactose sensor, Gal3p, the negative feedback node within the GAL regulatory network, Gal80p, and the hexose transporters, HXT, as the main sources of the bimodal range variation. We were able to switch the phenotype of individual yeast strains in silico by tuning parameters related to these three elements. Determining the basis for these behavioral differences may give insight into how the GAL pathway processes information, and into the evolution of nutrient metabolism preferences in different strains. More generally, our method of identifying the key parameters that explain phenotypic variation in this system should be generally applicable to other quantitative traits.     

The regulatory roles of the galactose permease and kinase in the induction response of the GAL network in Saccharomyces cerevisiae

The GAL genetic switch of Saccharomyces cerevisiae exhibits an ultrasensitive response to the inducer galactose as well as the "all-or-none" behavior characteristic of many eukaryotic regulatory networks. We have constructed a strain that allows intermediate levels of gene expression from a tunable GAL1 promoter at both the population and the single cell level by altering the regulation of the galactose permease Gal2p. Similar modifications to other feedback loops regulating the Gal80p repressor and the Gal3p signaling protein did not result in similarly tuned responses, indicating that the level of inducer transport is unique in its ability to control the switch response of the network. In addition, removal of the Gal1p galactokinase from the network resulted in a regimed response due to the dual role of this enzyme in galactose catabolism and transport. These two activities have competing effects on the response of the network to galactose such that the transport effects of Gal1p are dominant at low galactose concentrations, whereas its catabolic effects are dominant at high galactose concentrations. In addition, flow cytometry analysis revealed the unexpected phenomenon of multiple populations in the gal1delta strains, which were not present in the isogenic GAL1 background. This result indicates that Gal1p may play a previously undescribed role in the stability of the GAL network response.     

Characterizing the memory of the GAL regulatory network in Saccharomyces cerevisiae

Genetic regulatory networks respond dynamically to perturbations in the intracellular and extracellular environments of an organism. The GAL system in the yeast Saccharomyces cerevisiae has evolved to utilize galactose as an alternative carbon and energy source, in the absence of glucose in the environment. We present a dynamic model for GAL system in Saccharomyces cerevisiae, which includes a novel mechanism for Gal3p activation upon induction with galactose. The modification enables the model to simulate the experimental observation that in absence of galactose, oversynthesis of Gal3p can also induce the GAL system. We then characterize the memory of the GAL system as the domain of attraction of the steady states.     

Estimating the Stochastic Bifurcation Structure of Cellular Networks

High throughput measurement of gene expression at single-cell resolution, combined with systematic perturbation of environmental or cellular variables, provides information that can be used to generate novel insight into the properties of gene regulatory networks by linking cellular responses to external parameters. In dynamical systems theory, this information is the subject of bifurcation analysis, which establishes how system-level behaviour changes as a function of parameter values within a given deterministic mathematical model. Since cellular networks are inherently noisy, we generalize the traditional bifurcation diagram of deterministic systems theory to stochastic dynamical systems. We demonstrate how statistical methods for density estimation, in particular, mixture density and conditional mixture density estimators, can be employed to establish empirical bifurcation diagrams describing the bistable genetic switch network controlling galactose utilization in yeast Saccharomyces cerevisiae. These approaches allow us to make novel qualitative and quantitative observations about the switching behavior of the galactose network, and provide a framework that might be useful to extract information needed for the development of quantitative network models.     

Evidence for Gal3p's cytoplasmic location and Gal80p's dual cytoplasmic-nuclear location implicates new mechanisms for controlling Gal4p activity in Saccharomyces cerevisiae

Genetics and in vitro studies have shown that the direct interaction between Gal3p and Gal80p plays a central role in galactose-dependent Gal4p-mediated GAL gene expression in the yeast Saccharomyces cerevisiae. Precisely how Gal3p-Gal80p interaction effects induction is not clear. It has been assumed that Gal3p interacts with Gal80p in the nucleus upon galactose addition to release Gal80p inhibition of Gal4p. Although Gal80p has been shown to possess nuclear localization signal (NLS) peptides, the subcellular distribution of neither Gal80p nor Gal3p was previously determined. Here we report that Gal3p is located in the cytoplasm and apparently excluded from the nucleus. We show that Gal80p is located in both the cytoplasm and the nucleus. Converting Gal80p into a nucleus-localized protein (NLS-Gal80p) by exogenous NLS addition impairs GAL gene induction. The impaired induction can be partially suppressed by targeting Gal3p to the nucleus (NLS-Gal3p). We document a very rapid association between NLS-Gal3p and Gal80p in vivo in response to galactose, illustrating that the nuclear import of Gal80p is very rapid and efficient. We also demonstrate that nucleus-localized NLS-Gal80p can move out of the nucleus and shuttle between nuclei in yeast heterokaryons. These results are the first indication that the subcellular distribution dynamics of the Gal3 and Gal80 proteins play a role in regulating Gal4p-mediated GAL gene expression in vivo.     

Galactose induction in yeast involves association of Gal80p with Gal1p or Gal3p

Gal1p carries out two functions in the galactose pathway of yeast. It activates Gal4p by interacting with Gal80p – a function that can also served by Gal3p – and it catalyzes the formation of galactose-1-phosphate. Recently, we and others have presented biochemical evidence for complex formation between Gal1p and Gal80p. Here, we extend these data and present genetic evidence for an interaction between Gal1p and Gal80p in vivo, using a two-hybrid assay. Interaction between Gal1p and Gal80p depends on the presence of galactose, but not on the catalytic activity of Gal1p. A new class of Kluyveromyces lactis mutants was isolated, designated Klgal1-m, which have lost the derepressing activity but retain galactokinase activity, indicating that the two Gal1p activities are functionally independent. The KlGal1-m proteins are defective in their ability to interact with Gal80p in a two-hybrid assay. The locations of gal1-m mutations identify putative interaction sites in Gal1p and Gal80p. A dominant mutation, KlGAL1-d, leads to a high level of constitutive expression of genes of the galactose pathway. The behavior of chimeric proteins consisting of Gal3p and KlGal1p sequences indicates that both the N-terminal and C-terminal halves of KlGal1p are involved in specific interaction with KlGal80p. Received: 12 November 1998 / Accepted: 18 December 1998     

Analysis of the galactose signal transduction pathway in Saccharomyces cerevisiae: interaction between Gal3p and Gal80p.

The GAL3 gene plays a critical role in galactose induction of the GAL genes that encode galactose- metabolizing enzymes in Saccharomyces cerevisiae. Defects in GAL3 result in a long delay in GAL gene induction, and overproduction of Gal3p causes constitutive expression of GAL. Here we demonstrate that concomitant overproduction of the negative regulator, Gal80p, and Gal3p suppresses this constitutive GAL expression. This interplay between Gal80p and Gal3p is direct, as tagged Gal3p coimmunoprecipitated with Gal80p. The amount of coprecipitated Gal80p increased when GAL80 yeast cells were grown in the presence of galactose. When both GAL80 and GAL3 were overexpressed, the amount of coprecipitated Gal80p was not affected by galactose. Tagged gal3 mutant proteins bound to purified Gal80p, but only poorly in comparison with the wild type, suggesting that formation of the Gal80p-Gal3p complex depends on the normal function of Gal3p. Gal3p appeared larger in Western blots (immunoblots) than predicted by the published nucleic acid sequence. Reexamination of the DNA sequence of GAL3 revealed several mistakes, including an extension at the 3' end of another predicted 97 amino acids.     

Glucose-induced monoubiquitination of the Saccharomyces cerevisiae galactose transporter is sufficient to signal its internalization

In Saccharomyces cerevisiae, the addition of glucose to cells growing on galactose induces internalization of the galactose transporter Gal2p and its subsequent proteolysis in the vacuole. Here we report that the essential step in Gal2p down-regulation is its ubiquitination through the Ubc1p-Ubc4p-Ubc5p triad of ubiquitin-conjugating enzymes and Npi1/Rsp5p ubiquitin-protein ligase. Moreover, Gal2p appears to be stabilized in mutant cells defective in the ubiquitin-hydrolase Npi2p/Doa4p, and the mutant phenotype can be reversed by overexpression of ubiquitin. An analysis of the fate of Gal2p in cells overexpressing wild-type ubiquitin as well as its variants incompetent to form polyubiquitin chains showed that monoubiquitination of Gal2p is sufficient to signal internalization of the protein into the endocytic pathway.     

Multi-component gene network design as a survival strategy in diverse environments

Background

Gene-environment interactions are often mediated though gene networks in which gene expression products interact with other network components to dictate network activity levels, which in turn determines the fitness of the host cell in specific environments. Even though a gene network is the right context for studying gene-environment interactions, we have little understanding on how systematic genetic perturbations affects fitness in the context of a gene network.

Results

Here we examine the effect of combinatorial gene dosage alterations on gene network activity and cellular fitness. Using the galactose utilization pathway as a model network in diploid yeast, we reduce the copy number of four regulatory genes (GAL2, GAL3, GAL4, GAL80) from two to one, and measure the activity of the perturbed networks. We integrate these results with competitive fitness measurements made in six different rationally-designed environments containing different galactose concentrations representing the natural induction spectrum of the galactose network. In the lowest galactose environment, we find a nonlinear relationship between gene expression and fitness while high galactose environments lead to a linear relationship between the two with a saturation regime reached at a sufficiently high galactose concentration. We further uncover environment-specific relevance of the different network components for dictating the relationship between the network activity and organismal fitness, indicating that none of the network components are redundant.

Conclusions

These results provide experimental support to the hypothesis that dynamic changes in the environment throughout natural evolution is key to structuring natural gene networks in a multi-component fashion, which robustly provides protection against population extinction in different environments.