- Research article
- Open Access
Multiway real-time PCR gene expression profiling in yeast Saccharomyces cerevisiae reveals altered transcriptional response of ADH-genes to glucose stimuli
© Ståhlberg et al; licensee BioMed Central Ltd. 2008
- Received: 17 September 2007
- Accepted: 16 April 2008
- Published: 16 April 2008
The large sensitivity, high reproducibility and essentially unlimited dynamic range of real-time PCR to measure gene expression in complex samples provides the opportunity for powerful multivariate and multiway studies of biological phenomena. In multiway studies samples are characterized by their expression profiles to monitor changes over time, effect of treatment, drug dosage etc. Here we perform a multiway study of the temporal response of four yeast Saccharomyces cerevisiae strains with different glucose uptake rates upon altered metabolic conditions.
We measured the expression of 18 genes as function of time after addition of glucose to four strains of yeast grown in ethanol. The data are analyzed by matrix-augmented PCA, which is a generalization of PCA for 3-way data, and the results are confirmed by hierarchical clustering and clustering by Kohonen self-organizing map. Our approach identifies gene groups that respond similarly to the change of nutrient, and genes that behave differently in mutant strains. Of particular interest is our finding that ADH4 and ADH6 show a behavior typical of glucose-induced genes, while ADH3 and ADH5 are repressed after glucose addition.
Multiway real-time PCR gene expression profiling is a powerful technique which can be utilized to characterize functions of new genes by, for example, comparing their temporal response after perturbation in different genetic variants of the studied subject. The technique also identifies genes that show perturbed expression in specific strains.
- Principal Component Analysis
- Hierarchical Cluster
- Ethanol Production
- Glucose Addition
- Score Vector
The extraordinary sensitivity and virtually unlimited dynamic range of real-time PCR makes it the preferred technology for quantitative gene expression profiling. Using microarray technology expression of entire genomes can be measured, identifying candidates for expression profiling. After validation of these genes on representative samples by real-time PCR, eliminating any false leads and possibly complementing with other genes, powerful panels of expression markers can be identified. The recent development of high throughput real-time PCR platforms [1, 2] will spur the development further. To extract maximum information from profiling experiments using such panels, methods to pre-process and process the gene expression data are needed.
The addition of glucose to Saccharomyces cerevisiae cells grown in ethanol causes an extensive reprogramming of gene expression and metabolism, making it a suitable model system to understand gene regulation. In this system glucose consumption rate correlates with glucose repression [3, 4]. We have previously reported on a series of strains, in which glucose uptake is mediated by different native and chimeric hexose transporters, which display a wide range of glucose uptake rates [3, 5, 6]. These strains are useful for investigating the effects of different glycolytic rates on glucose-induced signaling pathways. Many glucose induced and glucose repressed genes have been extensively studied, but several genes believed to be associated with metabolism remain poorly characterized. The alcohol dehydrogenases (ADH) are such group. Its first two members, ADH1 and ADH2, have well known characteristics [7–10], while the functions of ADH3-6 are poorly understood [11–16].
We have previously shown that glucose uptake is the rate limiting step for glycolytic flux in strains expressing a series of individual glucose transporters with reduced transport capacity . Ethanol production, also under aerobic conditions, in Saccharomyces cerevisiae is believed to be a result of overflow metabolism where rate limitation of the TCA cycle results in a flux towards ethanol production . In this study we were interested in exploring transcriptional responses of some of the less characterized ADH-genes to better understand their regulations under conditions of different glycolytic rates.
Here, four yeast strains representing the full range of glycolytic rates; namely, wild-type (high glycolytic rate), HXT-HXT7 (medium glycolytic rate), HXT-TM6* (low glycolytic rate), and HXT-null (no glucose uptake) were selected to study the responses of metabolic genes. In previous works we have shown that ethanol production rate correlates to glycolytic rate , which pointed us in the direction of the ADH-genes and in particular the less studied ADH3, ADH4, ADH5 and ADH6. Here we address whether a decrease in the rate of ethanol production rate can be attributed in part to the roles of the ADH-genes. The study is a 3-way design, with the three ways being (i) genotype, (ii) gene, and (iii) time. We use a variant of geNorm  and Normfinder  to find suitable reference genes for normalization along all three ways, and we present suitable pre-processing of the data for analysis. Finally, the data are analyzed by augmented principal component analysis (PCA), which is a variant of PCA for 3-way studies . The classification obtained by the augmented PCA is verified by hierarchical clustering and clustering by the self-organizing map (SOM) . The analyses classify the genes into five groups with characteristic temporal profiles, based on which functional similarities between the ADH3-6 genes and previously more characterized genes can be found.
Experimental setup and gene selection
In this study four yeast strains KOY.PK2-1C83 (wild-type), KOY.HXT7P (HXT-HXT7), KOY.TM6*P (HXT-TM6*) and KOY.VW100P (HXT-null) were used. In the HXT-HXT7, HXT-TM6* and HXT-n ull strains, the genes HXT1-HXT7, GAL2, STL1 and three maltose permeases (see Methods) have been deleted. This results in a strain unable to take up glucose as measured by uptake of radiolabeled C13-glucose . Into this HXT-null strain we introduced either HXT7 (HXT-HXT7) or TM6* (HXT-TM6*). The four strains used have high, intermediate, low and zero glucose consumption rates, i.e. glycolytic rates [3, 5, 6]. They were grown on ethanol and exposed to glucose to follow the transcriptional responses of a series of selected genes. Samples were collected before and during one hour after the pulse to study the short term response of the genes.
Validation of reference genes
where CT0 is the CT-value measured immediately before glucose addition (0 minutes). RQ was then converted to fold changes (FC) with log2 base:
FC = log2(RQ)
The average expression of each gene in each strain is now zero and its standard deviation is one. The autoscaled expression profiles for the yeast wild-type strain are shown in Figure 4B (RQ and FCAS for the other strains are shown in Additional data 1).
Classification of the genes' expression profiles with principal component analysis
The measured data consisted of eight time points (samplings) measured for fifteen genes of interest in each of the four strains. Hence, the total number of data points was: 8 (samplings) × 15 (genes) × 4 (strains) = 480. Clearly, to visualize such large amounts of data and to unravel patterns efficient methods are called for. The classical scatter plot is an intuitive way to visualize how genes are expressed in different samples. These are typically 2-dimensional plots where the two axes indicate the genes' expressions in the two samples. It is possible to indicate genes' expressions for three genes in a 3-dimensional scatter plot. But this is the limit, since we have no convenient way to plot data in more than three dimensions. To deal with higher order data multivariate biostatistical tools are required to reduce the number of dimensions without loss of essential information. A classical, widely used tool is Principal Components Analysis (PCA). PCA allows scientists to study many variables simultaneously. It reflects how the original variables are correlated, and also how the samples are grouped. Principal Components (PCs) are mathematical constructs that can be interpreted as linear combinations of the studied variables with the following important properties:
(i) The PCs are orthogonal. Once a PC is linked to the behavior of one or several genes, one can be reasonably sure that this correlation is unique and these genes do not correlate substantially with other PCs. The numerical coefficients, ranging from -1 to +1, given to each gene in each PC are called loadings and reflect how important the gene is to define this PC.
(ii) The PCs are sequential. This means that the first, and by definition the most significant, PC can be interpreted as the line in the original multidimensional space of all the samples that best fits the expression data and, hence, explains most of the observed variability and accounts for most of the information. The second most significant PC can be visualized as a vector perpendicular to the first PC that fits the expression data best, and accounts for most of the variability that is not accounted for by the first PC. Additional PCs are defined analogously. One can extract PCs until a certain percentage of all the information, let's say 80%, is accounted for, and then discard the remaining PCs, which will mainly reflect uncorrelated information and, hence, the experimental noise. In most cases it is not practical to calculate more than three PC's, since, as already mentioned, we are limited to make 3-dimensional plots. Once the PCs have been calculated, the samples can be located in this new space using the scores, which specify the location of each sample on each PC. The original data can now be presented in a scores scatter plot to reveal groups among the samples or a loadings scatter plot to reveal groups among the genes. Many times PC1-PC2 scatter plots are sufficient, but one can also construct PC1-PC2-PC3 scatter plots.
Figure 4C shows the two most significant scores vectors for the autoscaled wild-type (high glycolytic rate) yeast data. As the temporal expression profiles, the scores are functions of time, and they reflect the main features of all the measured profiles. The most significant PC score vector (PC1) reflects a rapid decrease in expression that levels off after some 20 minutes. The PC2 score vector is signifying extreme expression at an intermediate time of about 10–15 minutes. The expression profile of every gene can be approximated as a linear combination of the two scores vectors, and can be visualized in a scatter plot based on the weights (loadings) of the linear combination (Figure 4D). For wild-type yeast we see that the induced (blue) and repressed (red) genes are clearly distinguished by PC1: induced genes have negative PC1 loadings while repressed genes have positive. Among the repressed genes, we find that SUC2 is located off the cluster's center, indicating that the SUC2 profile may be somewhat different from that of the bulk of the repressed genes. In Figures 4A and 4B, we indeed see that one of the red temporal expression profiles shows a different behavior from the rest. This is SUC2. While the expression of the bulk of the repressed genes reaches minimum at about 20 minutes and then saturates, SUC2 expression goes through a minimum and thereafter it slowly increases. Indeed, the PC2 loadings sort genes based on their tendency to show extreme expression at an intermediate time of about 10–15 minutes (the PC2 scores vector in Figure 4C). Genes showing local maximum expression at an intermediate time are characterized by a positive PC2 loading, while genes exhibiting local minimum have negative PC2 loading. Genes that do not show extreme expression at intermediate time points have PC2 loadings around zero. The most negative PC2 loading is found for CYC1 (figure 4D), which shows maximal repression after about 15 minutes after glucose addition (Figure 4B). Positive PC2 loadings are found for HSP12, ADH4 and ADH6, which all have a local maximum in their temporal expression profile. The specific roles for Adh4 and Adh6 are not fully understood. ADH4 has previously been reported not to be expressed in laboratory strains or to affect ethanol production . ADH6 has a high specificity towards long chain aliphatic and bulky substrates and has been suggested to participate in the production of fusel alcohols . Fusel alcohols are produced mostly during fermentation, which could explain its induction after glucose addition. The induced genes, with negative PC1 loadings, and most of the repressed genes, with positive PC1 loadings (SUC2 being an exception), have PC2 loadings close to zero, indicating that their temporal expression profiles are unimodal. ADH3 and ADH5 are located close to each other in the loadings scatter plot, within the cluster of the repressed genes. They have somewhat lower PC2 loadings than the bulk of the repressed genes, but higher than SUC2 (Figure 4D). This suggests that ADH3 and ADH5 respond similarly to glucose addition in wild-type yeast, and their expression profiles are characteristic of repressed genes. Repression of ADH3 after glucose addition is consistent with previous studies, which report lower Adh3 activity during respiratory growth when compared to fermentative growth . Adh5 has been reported to be involved in ethanol formation but its function is only apparent in an adh1adh3 strain . More detailed studies are needed to understand differences in responses between ADH4/6 and ADH3/5.
From the eigenvalues of the PCA it was calculated that 95.6 % of the variance in all the measured data is accounted for by the first two PC's, reinforcing the usefulness of the PC1 vs. PC2 loadings scatter plot. The third PC of the wild-type strain accounts for an additional 3% of the variability in the data. The third score vector has the shape of the derivative of the second vector (Figure 4C), and sorts the genes based on when they reach intermediate extreme expression. Genes that reach extreme expression after 10–15 minutes obtain positive PC3 loading, while genes reaching extreme expression within 10 minutes obtain negative PC3 loading. The genes in wild-type yeast are clustered based on all three PC's in the PC1 vs. PC2 vs. PC3 loadings scatter plot in Figure 4E. In the 3D loadings scatter plot, MIG1 separates from the other induced genes, because its expression reaches maximum at an earlier time. The 3-dimensional plot accounts for 98.6 % of all the variability, which is essentially all biologically relevant information; remaining variability is mainly experimental noise. This can be verified by comparing the wild-type temporal expression profiles with those reproduced from the three main PC's. The agreement is excellent evidencing that the three PC's have indeed picked up all the important features of the genes' expression profiles in wild-type yeast (data not shown). These features are:
1) Expression either increases or decreases
2) Expression reaches an extreme negative or positive level at an intermediate time point from which it recovers
3) Extreme expression is reached before or after 10 minutes.
PCs calculated for HXT-HXT7 (medium glycolytic rate) had similar features as those for the wild-type strain. In the PC1 vs. PC2 loadings scatter plot (92.6% of the initial variance, Figure 4F), the PC1 still differentiates between induced (left) and down-regulated (right) genes. The repressed genes cluster more tightly in the HXT-HXT7 strain compared to wild-type, and the cluster also contains ADH3 and HSP12, which in wild-type had a distinct location. This is because their expression in wild-type initially increases, goes through a maximum (Figure 4A–B) and then decreases. In HXT-HXT7 they are instead repressed from start. ADH5, which was found close to the cluster of repressed genes in wild-type, is now distant from it having large negative PC2 loading. This suggests that its expression is affected in HXT-HXT7 such that it goes through a minimum. CYC1, which too has negative PC2 score, also goes through a minimum in its temporal expression profile in the HXT-HXT7 strain. In addition, the induced genes separate into two subclusters: ADH1, TPI1 and PDC1 in one and PGK1, MIG1, ADH4 and ADH6 in the other (Figure 4F). The separation is mainly along PC2, indicating that the subgroups differ in whether expression goes through (in this case) a maximum or if it levels off.
For HXT-TM6* (low glycolytic rate)the contribution from noise to the data was larger due to the overall less efficient glucose response and the two main PC's account for only 85 % of the total variance. PC1 loadings still reflect whether genes are up or down-regulated, while the PC2 loadings are less well defined (Figure 4G). The repressed genes still form a cluster, with SUC2 having a somewhat higher PC2 score. HSP12 is located within the cluster of repressed genes, while ADH3 and ADH5, which were within this cluster in wild-type strain, are separated due to negative PC2 loadings and they are now closer to CYC1. The induced genes, but PGK1, form a tight cluster containing also ADH4. ADH6 is found at similar PC1 score as the induced genes, suggesting it too has similar expression, although its positive PC2 score indicates its expression profile is somewhat shifted in time. PGK1 is located around zero loadings for both PC1 and PC2. This means that the first two PCs are not very useful to describe its temporal profile and we can only conclude it is different from that of the other induced genes. Inspecting the raw data (see Additional data 1) we find that PGK1 expression in HXT-TM6* is essentially unaffected by glucose addition and remains at a constant level.
For the HXT-null (no glucose uptake) strain the changes in expression upon glucose addition were small, and the PC loadings are quite different from those of the other strains. Therefore, the genes' locations in the HXT-null scatter plot cannot be compared to the previous ones. Anyway not much remains of the groups seen for the other strains, reflecting the relatively weak response in HXT-null when glucose is added. The most spectacular observation is that MIG1, which is a typical glucose-induced gene, is here clearly differentiated from the other glucose induced genes (Figure 4H and Additional data 1C). Inspection of the temporal profiles indicates that MIG1 is the only gene in this group that is still induced upon glucose addition; the expression of the other genes decreases when glucose is added. HSP12, which is repressed after the glucose pulse has previously been reported to be down-regulated at very low glucose concentrations , but our results rather suggest that the signal is extracellular. Our observation that MIG1 is still derepressed in HXT-null confirms that it is regulated through extracellular sensing, which previously has been suggested by Kaniak et al. .
I. PC1 << 0, PC2 ≈ 0; genes induced upon glucose addition.
II. PC1 >> 0, PC2 ≈ 0; genes repressed upon glucose addition.
III. PC1 ≈ 0; PC2 >> 2; genes with expression profiles that go through a maximum.
IV. PC1 ≈ 0; PC2 << 2; genes with expression profiles that go through a minimum.
V. PC1 ≈ 0; PC2 ≈ 0; genes with no regulation upon glucose addition.
Most of the ADH1, PGK1, TPI1, PDC1 and MIG1 genes are found in area I, indicating they are induced upon glucose addition. Exceptions are PGK1 and MIG1 in the HXT-HXT7 strain, which go through a maximum, and PGK1 in HXT-TM6*, which shows no response. Most FBP1, ADH2, MDH2 and SUC2 genes are in area II. Exceptions are SUC2 in wild-type, where it shows a local minimum in the temporal profile, and in HXT-TM6*, where it shows a transient induction similar to that observed in growing cells close to glucose depletion . The less well understood ADH3-6 genes show the following behavior: ADH4 and ADH6 are found in areas I and III co-localized with genes induced by glucose, while ADH3 and ADH5 are found in areas II and IV co-localized with genes repressed by glucose. Hence, we conclude this is their general response to glucose. Several of the ADH3-6 genes in some strains have PC2 scores significantly different from zero indicating that their temporal profiles may go through local minima/maxima. To verify these conclusions the entire experiment was repeated and analyzed separately by matrix-augmented PCA (Additional data 3). Only small differences in scores were seen and the genes were grouped the same way as in Figure 5
Confirmation with hierarchical clustering and Kohonen self-organizing maps
While PCA is a very robust approach to classify samples based on multivariate and multiway measurements and an excellent tool to unravel variable patterns, there are also other techniques for unsupervised clustering. The most common is hierarchical clustering although, recently, the Kohonen self-organizing map (SOM) is gaining attention. While PCA always yields the same unbiased result for a certain set of data, hierarchical clustering requires the user to select a distance or similarity measure and also to define how distances between groups shall be measured. Once these decisions are made, hierarchical clustering also gives the same result every time for certain data. SOM, on the other hand, is based on a particular type of artificial neural networks that can be used to create an organized map of expression profiles by treating the raw information in the experimental data using a chain of successive comparisons. The goal is to create a map where adjacent areas correspond to similar samples, alternatively to genes with similar expressions. A SOM consists of an array of unconnected artificial neurons. There are several options to organize them but in most cases they are arranged in a convenient square so that each neuron is adjacent to other neurons (termed neighbors). The underlying idea is to assign sets of neurons (regions in the map) to a distinct class of samples. This is achieved through an iterative process. Briefly, each neuron in the map is defined as a set of weights (a vector of values) that equals the number of variables, in this case temporal profiles, measured for each gene . Hence, the weight vectors can be interpreted as artificial expression profiles. Creating the SOM consists of adjusting the weights of the neuron during a training phase. First the map is initiated by assigning small random numbers to the weights. Then a gene is selected randomly and its measured temporal expression profile is compared to the weights of the SOM's neurons. The neuron with the most similar weights is identified and information about the gene's expression profile is added to its weights. Furthermore, the information is also added to the neighboring neurons. The process is then repeated with a new gene over and over until a stable SOM is obtained. This SOM will have neurons with weights that reflect all genes' expression profiles and weights of neighboring neurons will be similar. In a final step each gene is placed in the neuron of the SOM that has the most similar weight to its temporal profile. As a consequence genes with similar profiles will be found in the same neuron or in close-lying neurons.
The possible advantages of hierarchical clustering and SOM when compared to PCA are that the former methods are intuitive and use all the information in the data, while PCA is based on coordinate transformation and reduction of dimensionality, which are less obvious operations and also always loses some of the variability in the measured data. The loss of limited amount of information is usually not serious, at least if the first 2–3 PCs have large eigenvalues and, hence, account for a large percentage of the total variation. In fact, reduction of dimensionality may lead to improved signal to noise ratio in the measured data, since biological variability is systematic and is mainly contained in the initial PCs, while experimental noise is random and dominates the PCs that are discarded . A disadvantage of hierarchical and Kohonen-SOM clustering is that they only yield groups, and any further biochemical interpretation has to be extracted studying the original data matrices, while in PCA there is a relation between the measured data and the scores and loadings.
Classifications by PCA, hierarchical clustering, and SOM are based on quite different assumptions and mathematical models, which makes the methods complementary. Hierarchical clustering and SOM create groups based on some measure of similarity between the samples, which is calculated directly from the experimental data. Additional criteria are required to construct groups (hierarchical clustering) or define regions for groups (SOMs). In contrast, PCA is based on the variables and its main objective is to reveal patterns by calculating a set of abstract factors (the PCs). The number of PCs is much lower than the number of variables, which makes interpretation easier. It eliminates most of the random variation but also some systematic, biologically relevant, information may be missed. SOM takes account of all variation in the data, but at the expense of not having a linear scale. This makes SOM more suitable for validation than first hand classification. Hierarchical clustering creates clusters sequentially by inspecting subpopulations of the data. Once a sample is entered into a group it cannot be extracted again (i.e. hierarchical). Hence, even though all information in the measurement is considered, it is not considered all at once. Therefore, the final clusters may depend on subtle differences between samples' expression profiles. For this reason results of hierarchical clustering should be confirmed by an independent method. For the data presented here and also for the independent replicated experiment in Additional data 3, classifications by PCA, hierarchical clustering and SOM are highly consistent, suggesting that the conclusions reached based on the results are valid.
Strains and growth conditions
Strains used in this study were KOY.2-1C83 (MATa MAL2-8c SUC2), KOY.VW100P (MATa MAL2-8c SUC2 hxt17Δ ura3-52 gal2Δ::loxP stl1Δ::loxP agt1Δ::loxP ydl247wΔ::loxP yjr160cΔ::loxP hxt12Δ::loxP hxt15Δ::loxP hxt16Δ::loxP hxt14Δ::loxP hxt12Δ::loxP hxt9Δ::loxP hxt11Δ::loxP hxt10Δ::loxP hxt8Δ::loxP hxt514:loxP hxt2Δ::loxP hxt367Δ::loxP carrying an integration cassette at the former HXT367 site containing the truncated, constitutive promoter of HXT7  (the KlURA3 open reading frame for counter selection, and the HXT7 terminator), KOY.HXT7P (same as KOY.VW100P but klURA3, in the integration cassette, was replaced by HXT7 and ura3-52::URA3) and KOY.TM6*P (same as KOY.VW100P but klURA3 in the integration cassette was replaced by TM6* and ura3-52::URA3) (3-4). Detailed strain information has been reported previously [3, 5, 6]. Cells were grown in 1% ethanol and 5× concentration of minimal media  until an optical density (OD) of 1–1.5 at 610 nm. Glucose was pulsed to a final concentration of 5% (w/v) and samples for RNA extraction were taken at 0 min just before the glucose pulse and then at 1, 5, 10, 20, 30 and 60 min after the pulse.
RNA extraction, cDNA synthesis and quantitative real-time PCR
RNA was extracted using phenol/chloroform extraction. RNA was dissolved in water and 100 μg was DNase treated on RNeasy columns (QIAGEN) as described by the manufacturer. 1 μg of the DNA free RNA was used in the reverse transcriptase reaction (Superscript II, Invitrogen) using pd(T)12–18 (Amersham Bioscience) as primers. Quantitative real-time assays were performed in an iCycler (BIORAD). PCR products were validated by agarose gel electrophoresis and melting curve analysis. Other settings are specified separately in the results. All real-time PCR experiments were run as duplicates. All data analyses, including normalization and data pre-treatment were performed with GenEx (version 4.1.7, MultiD Analyses). Data Supplement 1 and 3 require GenEx for visual inspection of the data.
This research was in part supported by start-up grant to MK from the institute of Biotechnology of the Czech Academy of science (Institutional research Concept AV0Z50520701). AS was supported by a postdoctoral fellowship award from the Swedish Research Council. KE was supported by Lena Gustafsson through European Commission contract BIO4-CT98-0562, Swedish Research Council (621-2001-1988) and by Stefan Hohmann through European Commission project AMPKIN (LSHG-CT2005-51818).
- Elbing K, Ståhlberg A, Hohmann S, Gustafsson L: Transcriptional responses to glucose at different glycolytic rates in Saccharomyces cerevisiae. Eur J Biochem. 2004, 271: 4855-4864. 10.1111/j.1432-1033.2004.04451.x.PubMedView ArticleGoogle Scholar
- Ye L, Kruckeberg AL, Berden JA, van Dam K: Growth and glucose repression are controlled by glucose transport in Saccharomyces cerevisiae cells containing only one glucose transporter. J Bacteriol. 1999, 181: 4673-4675.PubMedPubMed CentralGoogle Scholar
- Elbing K, Larsson C, Bill RM, Albers E, Snoep JL, Boles E, Hohmann S, Gustafsson L: Role of hexose transport in control of glycolytic flux in Saccharomyces cerevisiae. Appl Environ Microbiol. 2004, 70: 5323-5330. 10.1128/AEM.70.9.5323-5330.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Otterstedt K, Larsson C, Bill RM, Ståhlberg A, Boles E, Hohmann S, Gustafsson L: Switching the mode of metabolism in the yeast Saccharomyces cerevisiae. EMBO Rep. 2004, 5: 532-537. 10.1038/sj.embor.7400132.PubMedPubMed CentralView ArticleGoogle Scholar
- Leskovac V, Trivic S, Pericin D: The three zinc-containing alcohol dehydrogenases from baker's yeast, Saccharomyces cerevisiae. FEMS Yeast Res. 2002, 2: 481-494.PubMedGoogle Scholar
- Denis CL, Ferguson J, Young ET: mRNA levels for the fermentative alcohol dehydrogenase of Saccharomyces cerevisiae decrease upon growth on a nonfermentable carbon source. J Biol Chem. 1983, 258: 1165-1171.PubMedGoogle Scholar
- Ciriacy M: Genetics of alcohol dehydrogenase in Saccharomyces cerevisiae. II. Two loci controlling synthesis of the glucose-repressible ADH II. Mol Gen Genet. 1975, 138: 157-164. 10.1007/BF02428119.PubMedView ArticleGoogle Scholar
- Voronkova V, Kacherovsky N, Tachibana C, Yu D, Young ET: Snf1-dependent and Snf1-independent pathways of constitutive ADH2 expression in Saccharomyces cerevisiae. Genetics. 2006, 172: 2123-2138. 10.1534/genetics.105.048231.PubMedPubMed CentralView ArticleGoogle Scholar
- Drewke C, Ciriacy M: Overexpression, purification and properties of alcohol dehydrogenase IV from Saccharomyces cerevisiae. Biochim Biophys Acta. 1988, 950: 54-60.PubMedView ArticleGoogle Scholar
- Larroy C, Rosario Fernandez M, Gonzalez E, Pares X, Biosca JA: Properties and functional significance of Saccharomyces cerevisiae ADHVI. Chem Biol Interact. 2003, 1: 229-238. 10.1016/S0009-2797(02)00166-7.View ArticleGoogle Scholar
- Smith MG, Des Etages SG, Snyder M: Microbial synergy via an ethanol-triggered pathway. Mol Cell Biol. 2004, 24: 3874-3884. 10.1128/MCB.24.9.3874-3884.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Young ET, Pilgrim D: Isolation and DNA sequence of ADH3, a nuclear gene encoding the mitochondrial isozyme of alcohol dehydrogenase in Saccharomyces cerevisiae. Mol Cell Biol. 1985, 5: 3024-3034.PubMedPubMed CentralView ArticleGoogle Scholar
- Petersson A, Almeida JR, Modig T, Karhumaa K, Hahn-Hägerdal B, Gorwa-Grauslund MF, Lidén G: A 5-hydroxymethyl furfural reducing enzyme encoded by the Saccharomyces cerevisiae ADH6 gene conveys HMF tolerance. Yeast. 2006, 23: 455-464. 10.1002/yea.1370.PubMedView ArticleGoogle Scholar
- Larroy C, Fernández MR, González E, Parés X, Biosca JA: Characterization of the Saccharomyces cerevisiae YMR318C (ADH6) gene product as a broad specificity NADPH-dependent alcohol dehydrogenase: relevance in aldehyde reduction. Biochem J. 2002, 361: 163-172. 10.1042/0264-6021:3610163.PubMedPubMed CentralView ArticleGoogle Scholar
- Kappeli O: Regulation of carbon metabolism in Saccharomyces cerevisiae and related yeasts. Adv Microbial Physiol. 1986, 28: 181-208.View ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: reasearch0034-10.1186/gb-2002-3-7-research0034.View ArticleGoogle Scholar
- Andersen CL, Ledet-Jensen J, Ørntoft T: Normalization of real-time quantitative RT-PCR data: a model based variance estimation approach to identify genes suited for normalization – applied to bladder- and colon-cancer data-sets. Cancer Res. 2004, 64: 5245-5250. 10.1158/0008-5472.CAN-04-0496.PubMedView ArticleGoogle Scholar
- Jolliffe IT: Principal Component Analysis. 2002, Springer Series in Statistics, Springer, 2Google Scholar
- Kubista M, Sjögreen B, Forootan A, Sindelka R, Jonak J, Andrade JM: Real-time PCR gene expression profiling. European Pharmaceutical Reviews. 2007, 1: 56-60.Google Scholar
- DeRisi JL, Iyer VR, Brown PO: Exploring the metabolic and genetic control of gene expression on a genomic scale. Science. 1997, 278: 680-686. 10.1126/science.278.5338.680.PubMedView ArticleGoogle Scholar
- Wieczorke R, Krampe S, Weierstall T, Freidel K, Hollenberg CP, Boles E: Concurrent knock-out of at least 20 transporter genes is required to block uptake of hexoses in Saccharomyces cerevisiae. FEBS Lett. 1999, 464: 123-128. 10.1016/S0014-5793(99)01698-1.PubMedView ArticleGoogle Scholar
- Bu Y, Schmidt MC: Identification of cis-acting elements in the SUC2 promoter of Saccharomyces cerevisiae required for activation of transcription. Nucleic Acids Res. 1998, 26: 1002-1009. 10.1093/nar/26.4.1002.PubMedPubMed CentralView ArticleGoogle Scholar
- Rubenstein EM, McCartney RR, Schmidt MC: Regulatory domains of Snf1-activating kinases determine pathway specificity. Eukaryot Cell. 2006, 5: 620-627. 10.1128/EC.5.4.620-627.2006.PubMedPubMed CentralView ArticleGoogle Scholar
- Treitel MA, Kuchin S, Carlson M: Snf1 protein kinase regulates phosphorylation of the Mig1 repressor in Saccharomyces cerevisiae. Mol Cell Biol. 1998, 18: 6273-6280.PubMedPubMed CentralView ArticleGoogle Scholar
- Ozcan S, Vallier LG, Flick JS, Carlson M, Johnston M: Expression of the SUC2 gene of Saccharomyces cerevisiae is induced by low levels of glucose. Yeast. 1997, 13: 127-137. 10.1002/(SICI)1097-0061(199702)13:2<127::AID-YEA68>3.0.CO;2-#.PubMedView ArticleGoogle Scholar
- Schuller HJ: Transcriptional control of nonfermentative metabolism in the yeast Saccharomyces cerevisiae. Curr Genet. 2003, 43: 139-160.PubMedGoogle Scholar
- Treitel MA, Carlson M: Repression by SSN6-TUP1 is directed by MIG1, a repressor/activator protein. Proc Natl Acad Sci USA. 1995, 92: 3132-3136. 10.1073/pnas.92.8.3132.PubMedPubMed CentralView ArticleGoogle Scholar
- Guarente L, Lalonde B, Gifford P, Alani E: Distinctly regulated tandem upstream activation sites mediate catabolite repression of the CYC1 gene of Saccharomyces cerevisiae. Cell. 1984, 36: 503-511. 10.1016/0092-8674(84)90243-5.PubMedView ArticleGoogle Scholar
- Praekelt UM, Meacock PA: HSP12, a new small heat shock gene of Saccharomyces cerevisiae : analysis of structure, regulation and function. Mol Gen Genet. 1990, 223: 97-106. 10.1007/BF00315801.PubMedView ArticleGoogle Scholar
- Varela JC, Praekelt UM, Meacock PA, Planta RJ, Mager WH: The Saccharomyces cerevisiae HSP12 gene is activated by the high-osmolarity glycerol pathway and negatively regulated by protein kinase A. Mol Cell Biol. 1995, 15: 6232-6245.PubMedPubMed CentralView ArticleGoogle Scholar
- de Groot E, Bebelman JP, Mager WH, Planta RJ: Very low amounts of glucose cause repression of the stress-responsive gene HSP12 in Saccharomyces cerevisiae. Microbiology. 2000, 146: 367-375.PubMedView ArticleGoogle Scholar
- Kubista M, Andrade JM, Bengtsson M, Forootan A, Jonák J, Lind K, Sindelka R, Sjöback R, Sjögreen B, Strömbom L, Ståhlberg A, Zoric N: The Real-time Polymerase Chain Reaction. Mol Aspects Med. 2006, 27: 95-125. 10.1016/j.mam.2005.12.007.PubMedView ArticleGoogle Scholar
- Drewke C, Thielen J, Ciriacy M: Ethanol formation in adh0 mutants reveals the existence of a novel acetaldehyde-reducing activity in Saccharomyces cerevisiae. J Bacteriol. 1990, 172: 3909-3917.PubMedPubMed CentralGoogle Scholar
- Kaniak A, Xue Z, Macool D, Kim JH, Johnston M: Regulatory network connecting two glucose signal transduction pathways in Saccharomyces cerevisiae. Eukaryot Cell. 2004, 3: 221-231. 10.1128/EC.3.1.221-231.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Felipe-Sotelo M, Andrade JM, Carlosena A, Tauler R: Temporal characterisation of river waters in urban and semi-urban areas using physico-chemical parameters and chemometric methods. Anal Chim Acta. 2007, 583: 128-137. 10.1016/j.aca.2006.10.011.PubMedView ArticleGoogle Scholar
- Melssen W, Wehrens R, Buydens L: Supervised Kohonen networks for classification problems. Chemom Intell Lab Syst. 2006, 83: 99-113. 10.1016/j.chemolab.2006.02.003.View ArticleGoogle Scholar
- Wise BM, Kowalski BR: Process Chemometrics. Process Analytical Chemistry. Edited by: McLennan F, Kowalski BR. 1995, London: Blackie Academic Press, 257-312.Google Scholar
- D'haeseleer P, Liang S, Somogyi R: Genetic network inference: from co-expression clustering to reverse engineering. Bioinformatics. 2000, 16: 707-726. 10.1093/bioinformatics/16.8.707.PubMedView ArticleGoogle Scholar
- Hauf J, Zimmermann FK, Muller S: Simultaneous genomic overexpression of seven glycolytic enzymes in the yeast Saccharomyces cerevisiae. Enzyme Microb Technol. 2000, 26: 688-698. 10.1016/S0141-0229(00)00160-5.PubMedView ArticleGoogle Scholar
- Verduyn C, Stouthamer AH, Scheffers WA, van Dijken JP: A theoretical evaluation of growth yields of yeasts. Antonie Van Leeuwenhoek. 1991, 59: 49-63. 10.1007/BF00582119.PubMedView ArticleGoogle Scholar
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.