<title>Background</title> <p>We engineered the first respiratory <italic>Saccharomyces cerevisiae </italic>strain, KOY.TM6*P, by integrating the gene encoding a chimeric hexose transporter, Tm6*, into the genome of an <italic>hxt </italic>null yeast [##REF##15071495##1##]. Subsequently we demonstrated the transferability of this respiratory phenotype to a yeast strain, V5 <italic>hxt1-7</italic>Δ, in which only <italic>HXT1-7 </italic>had been deleted [##REF##16204537##2##]. The resulting V5.TM6*P strain produced only minor amounts of ethanol when cultured on external glucose concentrations as high as 5% [##REF##16204537##2##]. Despite the fact that the glucose flux was reduced to 30% in the V5.TM6*P strain compared with that of its parental V5 strain, the V5.TM6*P strain produced biomass at a specific rate as high as 85% of that produced by the V5 wild-type strain, the yield itself increasing by 50% in the mutant compared to its parental strain.</p> <p>Having performed this initial physiological characterization of the V5.TM6*P strain, we wanted to examine more thoroughly the basis of its respiratory phenotype by comparing its transcriptome with that of its wild-type parent strain, V5. This approach has been widely used to obtain a global picture of the cellular responses of yeast to a wide range of physiological changes, including those experienced at the diauxic shift [##REF##15987902##3##, ####REF##15758028##4##, ##REF##9381177##5####9381177##5##]. Recently Ohlmeier and colleagues [##REF##14597615##6##] demonstrated that major changes at the transcriptional level are not reflected at the protein level. Indeed, focusing on the mitochondrial proteome following a diauxic shift, they showed that the levels of only 18 out of 253 identified proteins had changed (17 increased, and 1 decreased). Among them were proteins involved in the tricarboxylic acid (TCA) cycle (Sdh1, Sdh2, and Sdh4) and the respiratory chain (Cox4, Cyb2, and Qcr7). This seeming disagreement with their observations of heterogeneous changes in the transcriptome (where more than 4,000 changes were recorded) is also consistent with our own prior data where we showed that the protein patterns obtained from a 2D-PAGE analysis of the two strains (grown under the same conditions used in this study) were not substantially different [##REF##16204537##2##]. We did note, however, that the levels of the upper-part glycolytic intermediates and ATP did differ between the two strains, and that Cdc19 is present in lower amounts in V5.TM6*P. This discrepancy between proteome and transcriptome regulation has also been observed by others [##REF##15642185##7##] and is also consistent with the fact that there is an acknowledged lack of relationship between the levels of glycolytic proteins and the glycolytic flux [##REF##1533788##8##, ####REF##2528863##9##, ##REF##10862878##10####10862878##10##]. Nonetheless, researchers continue to depend on mRNA as an indicator of cellular state, which is a situation that will continue while methods for the global analysis of protein expression are improved [##REF##15642185##7##].</p> <p>In this respect, Brauer and colleagues [##REF##15987902##3##, ####REF##15758028##4##, ##REF##9381177##5####9381177##5##] have used transcriptome analysis to examine adjustments and remodelling of the metabolism of glucose-limited yeast cultures. In that study, the authors constructed metabolic models consistent with their cDNA array data. In this work, the V5.TM6*P strain was therefore characterized at the transcriptome level in comparison with the parental V5 strain in order to obtain a more global understanding of any genetic remodelling underlying the physiological differences between them. Both strains were cultured in 50 g/L glucose with six samples being taken at defined points between 38 and 6 g/L glucose for cDNA array analysis in two independent fermentations for each strain. Consequently, we generated a highly redundant dataset (n = 12) where 922 genes were either induced or repressed when comparing V5 and V5.TM6*P under all glucose concentrations. Using <italic>in vivo </italic>studies together with documented binding specificities of three transcription factors known to have a role in the regulation of respiratory genes (the Hap complex, Cat8 and Mig1), it was possible to correlate the up-regulation of all but 18 of the 146 genes induced by a factor ≥2.0 with the presence of one or more of their binding sites. This means that 75% of the transcriptional response of the genes in our dataset could be related to the activities of just three transcription factors, which in turn would account for 88% of the response of the induced set. Finally, we made the unexpected observation that both the induction and repression of many of the genes in our dataset had a clear glucose dependence in the parent V5 strain that was not present in V5.TM6*P. This important result indicates that the relief of glucose repression in the wild-type is already operable at much higher glucose concentrations than is currently widely accepted and that the phenotype of the respiratory strain is glucose insensitive by comparison.</p>
<title>Methods</title> <title>Yeast strains</title> <p>The auxotrophic wine yeast strain V5 (<italic>MAT</italic>a <italic>ura3 gal</italic>) and the V5 <italic>hxt1–hxt7</italic>Δ mutant strain (<italic>MAT</italic>a <italic>ura3 gal hxt514</italic>Δ::loxP <italic>hxt367</italic>Δ::lox P <italic>hxt2</italic>Δ::loxP) were kindly provided by the Professor Bruno Blondin, Montpellier, France. The construction of V5.TM6*P was previously described by Henricsson and colleagues [##REF##16204537##2##] where the expression cassette from KOY.TM6*P [##REF##15071495##1##] was transferred into the genome of the V5 <italic>hxt1–hxt7</italic>Δ strain by using the primers PROHXT3 (TCAAATGGCGGTGTAGTTTGAAAAG) and TERHXT7 (TTAAGTGACGGGCGATGAGTAAGAA). Transformants were selected by growth on glucose. The resulting <italic>ura3 </italic>strain is referred to as V5.TM6*. The auxotrophic V5 wild type and V5.TM6* used in this study were made prototrophic by integration of <italic>URA3 </italic>[##REF##10618490##37##] and are referred to as the V5 wild-type strain and V5.TM6*P, respectively.</p> <title>Growth conditions</title> <p>For both pre-cultures and main cultures, 5 × defined minimal medium [##REF##1523884##38##] with 5% glucose as the sole carbon and energy source was used. A two-step pre-cultivation was performed on a rotary shaker. First a 10 mL culture was grown for 72 h, and used to inoculate 100 mL medium, which after 30 h at 30°C was inoculated to give a final OD<sub>610 </sub>of 0.05 in a bioreactor (BRO2, Belach Bioteknik AB, Sweden) in a volume of 2.5 L. Polypropylene glycol P2000 was added as antifoam (100 <italic>μ</italic>L/L). The cultivation conditions were 30°C, 1,000 rpm, pH 5.0 with an airflow of 1.25 L/min. Gas evolution was monitored on-line (type CP460 O<sub>2</sub>/CO<sub>2</sub>, Belach Bioteknik AB). Cultivations were performed at least in duplicate.</p> <title>RNA preparation</title> <p>Yeast cells (2 – 20 mL) were directly transferred into twice their volume of ice-cold water. The cells were collected by centrifugation at 0°C for 5 min at 3000 × g, frozen in liquid nitrogen and stored at -80°C. Total RNA was then prepared using the RNeasy kit from Qiagen, following the manufacturer's instructions.</p> <title>cDNA array production</title> <p>PCR products for the yeast ORFs were generated from yeast DNA using the ORF-specific primer set from Research Genetics and then re-amplified with the Resgen universal forward and reverse primers (Catalogue number 40612). PCR products were prepared by Randy Strich (Fox Chase Cancer Centre, USA) as part of a collaboration with the Wistar Institute, Philadelphia, USA. PCR products were spotted with a BioRobotics Microgrid TASII (Genomic Solutions) by the Wistar Genomics Facility. Each array contained 6,319 individual yeast genes and 1,169 gene repeats. The complete list of arrayed genes on array YA04 is listed under GEO Series accession number GPL4423.</p> <title>RNA amplification, hybridization, and scanning</title> <p>Amplified RNA (aRNA) was prepared from 1 <italic>μ</italic>g of total RNA for each sample tested using the RiboAmp kit (Arcturus). Labelled targets were prepared from 1.6 <italic>μ</italic>g aRNA with Superscript II reverse transcriptase (Invitrogen), in the presence of 3,000–5,000 Ci/mM [α-<sup>33</sup>P] dCTP (Amersham Pharmacia Biotech), 1 mmol/L dATP, 1 mmol/L dTTP, 1 mmol/L dGTP, 67 ng/<italic>μ</italic>g oligo-dT (Promega Biosciences), and 0.65 × random decamer primers (Ambion). Labelled targets were hybridized to individual arrays at 42°C for 18 h in 3 mL Microhyb buffer (Invitrogen). Arrays were washed twice in 2 × standard saline citrate (SSC)/1% SDS solution for 30 min at 50°C, once in 0.5 × SSC/1% SDS and once in 0.1 × SSC/0.5% SDS for 30 min at 55°C. The arrays were exposed to phosphor screens (Amersham Biosciences) for 6 days and scanned in a Storm 820 PhosphorImager (Molecular Dynamics). Quantitation of each spot was assessed by Imagene 5.0 software (BioDiscovery Incorporated) using manual spot alignment, measuring median pixel intensity for each spot and subtracting the local background. The data have been deposited in the NCBI Gene Expression Omnibus and are accessible through GEO Series accession number GSE11799.</p> <title>cDNA array data analysis</title> <p>Observed values less than 0.1 (5% (8752 values) of the 176640 raw data values) were set to 0.1, in accordance with information from Wistar that variations in values below this threshold should be regarded as measurement noise. After converting to a logarithmic scale, average expression values were computed over spots with the same ORF, and arrays were scale normalized by subtracting the median for each of the 24 arrays. Averages were computed over repeated measurements with the same conditions producing 6 values at different glucose concentrations for each of the V5 (38.4, 35, 25.5, 24.4, 10.5, 7.4 g/L glucose) and V5.TM6*P (36.5, 34.8, 26.3, 25.1, 13.5, 5.6 g/L glucose) strains for each gene. T-tests were then performed, comparing the two strains, using Empirical Bayes robustification of variance estimates [##REF##16904203##39##,##UREF##0##40##].</p> <p>The T-tests were adjusted for multiple testing using the false discovery rate method (disregarding correlations between genes). Essentially, then, p-values were adjusted so that when selecting all genes with p-values less than a threshold q, a proportion of q false positives would be expected amongst these genes. The estimated effects are shown as factors in Additional file ##SUPPL##0##1##, as the linear modelling was performed on a log scale. The p-values in Additional file ##SUPPL##0##1## are the adjusted values, so that selecting genes with p-values &lt; 0.05 gives an expected false discovery rate of 5%. Calculations were performed using the BioConductor package limma [##UREF##0##40##]. To present the data, a tabular format was chosen over heatmaps in an attempt to make it as accessible as possible to the scientific community.</p> <p>When producing data for Figure ##FIG##3##4##, a similar linear model was used, but with a different design matrix: Essentially, for each gene, 2 straight lines were simultaneously fitted, one for each of the strains, to plots of the 6 log-expression values versus glucose concentration. Empirical Bayes robustification and adjustment for multiple testing was included even in this analysis, using the limma tool [##UREF##0##40##]. Lines with slopes significantly different from zero were found for 20 genes for the V5 strain, but for no genes from the V5.TM6*P strain.</p>
<title>Results</title> <p>We compared the transcriptional profiles of the V5 and V5.TM6*P strains cultured in 50 g/L glucose, with samples extracted between 38 and 6 g/L residual glucose. We identified 190 genes with adjusted p-values &lt; 0.05 that changed their expression levels either up by a factor of ≥2.0 (146 genes) or down by a factor of ≤0.5 (44 genes) when comparing the two strains at all three sampling points. These genes are presented by functional category of their gene products in Additional file ##SUPPL##0##1## and are discussed below. Additional file ##SUPPL##1##2## provides the information in Excel format and also includes other relevant literature data. A further 636 genes changed their expression levels either up by a factor of &lt;2.0 (317 genes) or down by a factor &gt; 0.5 (319 genes) and are presented as part of Additional file ##SUPPL##2##3##. The remaining 96 genes were either dubious, Ty-transposable elements, or were detected due to cross-hybridisation, as indicated in the legend to Additional file ##SUPPL##0##1##.</p> <title>The majority of genes that have an altered expression in V5.TM6*P compared with its parental strain also have an altered expression following diauxic shift</title> <p>We compared our dataset with three sets of expression data for yeast undergoing a diauxic shift, where the samples had been taken at least 2 hours after the glucose is depleted [##REF##15987902##3##, ####REF##15758028##4##, ##REF##9381177##5####9381177##5##]. Of the 190 genes in our dataset that changed their expression levels either up or down by a factor of 2 or more, 154 of these could be identified as changing in the same direction in at least one of these three datasets, with 124 having at least two matches. Only 8 genes (<italic>ERG3</italic>, <italic>RPS27A</italic>, <italic>YHB1</italic>, <italic>TRX1</italic>, <italic>YER188w</italic>, <italic>YKL066w</italic>, <italic>YAR068w </italic>and <italic>Y0L150c</italic>) had changes that were in the opposite direction in at least two of the datasets. This significant overlap between our data set and the data sets of the three previous studies is entirely consistent with the fact that there is a transition from a respiro-fermentative metabolism in V5 to a respiratory metabolism in V5.TM6*P.</p> <title>Glycolytic and gluconeogenic genes have substantially altered expression levels in V5.TM6*P</title> <p>Four of the most strongly induced genes in Additional file ##SUPPL##0##1## are the glycolytic genes, <italic>HXK1 </italic>(factor 9.0) and <italic>GLK1 </italic>(3.8), and the gluconeogenic genes, <italic>PCK1 </italic>(6.9) and <italic>FBP1 </italic>(5.3). Indeed, as shown in Figure ##FIG##0##1##, the induction in V5.TM6*P of several glycolytic and gluconeogenic genes could be consistent with these cells accommodating for a lack of glucose; the hypothesized mechanism behind the strain's reduced ethanol-producing, respiratory phenotype. The final step in glycolysis which produces pyruvate, is catalyzed by pyruvate kinase. In V5.TM6*P the pyruvate kinase genes, <italic>CDC19 </italic>and <italic>PYK2</italic>, respond in an opposing manner: <italic>CDC19 </italic>is repressed by a factor of 0.4 (consistent with our earlier proteome analysis where Cdc19 is present in lower amounts in V5.TM6*P than V5) while <italic>PYK2 </italic>is induced and is relieved of its glucose-repressed state as a result of the low glucose flux in the respiratory strain. Reduced ethanol production in V5.TM6*P compared to the parental strain is also consistent with the two most strongly repressed genes <italic>PDC1 </italic>(0.2) and <italic>PDC5 </italic>(0.2) encoding the main pyruvate decarboxylases.</p> <title>Alternative carbon source usage genes are induced in V5.TM6*P compared with V5, consistent with a reduced glucose flux</title> <p>An adaptation to enable growth on alternative carbon sources is seen in the induction of the plasma membrane transporter genes <italic>MAL31 </italic>and <italic>STL1 </italic>as well as the genes <italic>SUC2</italic>, <italic>GUT1</italic>, <italic>GUT2</italic>, <italic>ALD6 </italic>and <italic>BDH1 </italic>encoding proteins for the utilization of sucrose, glycerol, acetaldehyde and butanediol, respectively. <italic>GUT1 </italic>in particular is strongly induced by a factor of 6.2, and <italic>STL1 </italic>is the most strongly induced of all genes on the cDNA array (21.7), suggesting that V5.TM6* has undergone a substantial genetic adjustment compared to its parent. Furthermore, <italic>CYB2 </italic>(10.1) and <italic>DLD3 </italic>(0.5), whose gene products are involved in the conversion of lactate to pyruvate, and in the case of <italic>CYB2 </italic>feeding the released electrons into the electron transport chain, are oppositely regulated. Biochemical data and sequence patterns for <italic>YGL157w </italic>(0.3) [##REF##14574691##11##], <italic>YAL061w </italic>(2.9) and <italic>DSF1 </italic>(9.4, homologous to <sc>D</sc>-mannitol 2-dehydrogenase from <italic>Rhodobacter sphaeroides) </italic>suggest that these ORFs encode an oxidoreductase, a polyol dehydrogenase and a putative sugar dehydrogenase, respectively, which could be involved in channelling alternative carbon sources into the TCA cycle. The carbon-deficient state of V5.TM6*P [##REF##15071495##1##] could also be consistent with the induction of <italic>CRC1 </italic>and <italic>CAT2 </italic>whose gene products are involved in <italic>β</italic>-oxidation and carnitine-dependent transport of acetyl-CoA from peroxisomes to mitochondria as a way of using carbon from sources other than glucose [##REF##11329169##12##]. The induction of <italic>SFC1</italic>, <italic>JEN1</italic>, <italic>FBP1 </italic>and <italic>ATO2 </italic>seen in our study has previously been reported to be a specific identifier of carbon-limited growth [##REF##12414795##13##].</p> <title>All TCA- and glyoxylate cycle genes have increased expression in V5.TM6*P compared with V5</title> <p>The genes encoding the enzymes of the TCA cycle (<italic>CIT1</italic>, <italic>ACO1, IDH1</italic>, <italic>KGD1</italic>, <italic>KGD2</italic>, <italic>LSC2</italic>, <italic>SDH1-4</italic>, <italic>FUM1</italic>, <italic>MDH1</italic>) all have a significant (adjusted p-value &lt; 0.05) increase in expression when comparing V5 and V5.TM6*P. In addition, the gene encoding an isoform of aconitase (<italic>ACO2</italic>) has its expression level reduced by a factor of 0.5. Fermentation is almost completely abolished in the V5.TM6*P strain, and we have previously shown that the TCA cycle has a higher carbon flux [##REF##16204537##2##] to harness the electron flow into the respiratory chain and generate a proton gradient using F<sub>0</sub>F<sub>1</sub>-ATPase for ATP synthesis (several respiratory chain genes are induced, as discussed below). Interestingly, all genes of the TCA-cycle are induced by greater than a factor of 3.5 except for <italic>CIT1</italic>/<italic>CIT2 </italic>and <italic>KGD1</italic>/<italic>KGD2 </italic>where the factor change is 2.2 – 3.5 for each gene of the pair, presumably yielding the same overall transcriptional response. Consistent with this, Ohlmeier and colleagues have shown that Sdh1, 2 and 4 all have increased expression after a diauxic shift [##REF##14597615##6##].</p> <p>The glyoxylate cycle is anabolic, yielding gluconeogenic precursors, and all its genes (<italic>CIT2</italic>, <italic>ACO1, ICL1 </italic>(induced by a factor of 1.5; Additional file ##SUPPL##2##3##), <italic>MLS1 </italic>and <italic>MDH2</italic>) are induced in V5.TM6*P compared with V5 (Additional file ##SUPPL##0##1##, Fig. ##FIG##0##1##). This could be consistent with the correspondingly increased biomass yield of the V5.TM6*P strain (Table ##TAB##0##1##). We note that V5.TM6*P has a lower biomass yield than expected from a totally respiratory phenotype. In chemostat culture under aerobic, fully respiratory conditions, a yield of 0.5 g/g is typical at low dilution rates, whereas the observed low yield in V5.TM6*P has already been published using a range of sugars and concentrations [##REF##16204537##2##]. Furthermore, numerous experiments on KOY.TM6*P, our original respiratory strain based on CEN.PK2-1C, [##REF##15071495##1##] have yielded identical results. These observations are in agreement with our findings in Table ##TAB##0##1## of 0.32 g/g in V5.TM6*P and 0.21 g/g in V5 from batch cultures.</p> <title>Genes encoding respiratory chain enzymes are induced in V5.TM6*P compared with V5</title> <p>33 key genes encoding mitochondrial enzymes involved in respiration – as well as F<sub>0</sub>F<sub>1</sub>-ATPase and its chaperones – are all induced by a factor of 2 – 4 in V5.TM6*P (Additional file ##SUPPL##0##1##). The gene encoding Ndi1 (Fig. ##FIG##1##2##), which transfers electrons from NADH to ubiquinone (Q in Fig. ##FIG##1##2##) is induced the most (by a factor of 4.1), while all other components of the respiratory chain are induced by a factor of 2 – 2.5. Genes encoding heme synthesis and incorporation (<italic>HEM2</italic>, <italic>CYT2</italic>), complex assembly (<italic>MBA1</italic>) and heme degradation (<italic>HMX1</italic>) are also induced by a factor of 2 – 2.5. Cytochrome c oxidase also has different chaperones that are induced to form an active structure by proteolysis (<italic>COX20</italic>), as well as delivery of copper (<italic>COX17</italic>). In spite of the higher respiratory rate in V5.TM6*P, several genes involved in detoxification of oxygen radicals are repressed, <italic>e.g. </italic>genes encoding peroxiredoxin (<italic>AHP1</italic>) and thioredoxin (<italic>TRX2</italic>) are repressed by a factor of 0.4 – 0.5. Conversely, the capacity for nitrous oxide detoxification is increased by the induction of <italic>YHB1 </italic>by a factor of 2.4. It is possible that the respiratory phenotype of V5.TM6*P could induce a higher level of oxygen radical formation inside its mitochondria and thus oxidation of guanine could occur on the mitochondrial DNA. Such detrimental effects could be counteracted, <italic>via </italic>the observed induction of <italic>OGG1 </italic>by a factor of 2.2.</p> <title>Changes of expression in genes encoding mitochondrial and plasma membrane transporters are consistent with a respiratory phenotype in V5.TM6*P</title> <p>The genes encoding the hexose transporters, Hxt1 to Hxt7, have very high sequence identities [##REF##8929273##14##] and thus cross reactivity between them is expected on the array. This explains the observed signal for <italic>HXT1, 4 </italic>and <italic>6</italic>, which have been deleted in V5.TM6*P. <italic>HXT13</italic>, <italic>HXT15 </italic>and <italic>HXT16 </italic>are known to be repressed by high levels of glucose and are consequently induced as expected in V5.TM6*P. A carbon source deficiency triggers an up-regulation of genes encoding plasma membrane proteins responsible for the import of alternative carbon sources, as already noted: <italic>STL1 </italic>(21.7; glycerol); <italic>JEN1 </italic>(10.2; lactate), <italic>PUT4 </italic>(5.8; proline – the main nitrogen source in grape juice) and <italic>MAL31 </italic>(3.1; maltose). A gene associated with nutrient deficiency in yeast [##REF##8139573##15##], <italic>YGP1</italic>, is induced by a factor of 2.6, again supporting the notion that V5.TM6*P has undergone a genetic adjustment compared to its parent.</p> <p>It seems that V5.TM6*P also has an altered requirement for metal ions in that genes encoding three metal ion channels have altered expression: <italic>CTR3 </italic>(copper transport) is induced by a factor of 3.0, while <italic>ZRT1 </italic>(zinc transport) and <italic>ENB1 </italic>(iron transport) are repressed by a factor of 0.4 – 0.5. It is likely that copper import is increased so that copper can be transferred to cytochrome c oxidase, which is the main copper-requiring enzyme in Additional file ##SUPPL##0##1##, with 3 copper ions per complex and a copy number of 5,000–12,000 per cell. There are only 45 genes in Additional file ##SUPPL##0##1## with greater than 5,000 copies per cell and none of these has copper as a prosthetic group except cytochrome c oxidase.</p> <p>The amino acid uptake systems encoded by <italic>BAP3 </italic>and <italic>GNP1 </italic>are repressed while the proline permease encoded by <italic>PUT4 </italic>is induced. <italic>PUT4 </italic>has a carbon source responsive element (CSRE) [##REF##11024040##16##] that is expected to be induced at the low glucose flux experienced in the V5.TM6*P strain, while <italic>BAP3 </italic>and <italic>GNP1 </italic>are regulated by a sensor of extracellular amino acid concentration, where Ssy5 [##REF##15611869##17##] is a vital component. It is likely that in V5.TM6*P the repression of <italic>BAP3 </italic>and <italic>GNP1 </italic>is a result of the repression (by a factor 0.65; Additional file ##SUPPL##2##3##) of <italic>SSY5 </italic>seen on the cDNA array. Interestingly, induction of genes encoding the ammonia transport system in the V5.TM6*P strain (<italic>ATO2</italic>, <italic>MEP1</italic>) could suggest that the strain is experiencing general nutrient limitation. It has previously been reported for example that the induction of another <italic>MEP </italic>family member, <italic>MEP2</italic>, is associated with nitrogen limitation [##REF##12414795##13##].</p> <title>88% of genes induced in V5.TM6*P compared to its parent have transcription factor binding sites for one or more of the Hap complex, Cat8 or Mig1</title> <p>Genome-wide studies have yielded a wealth of information on which promoters are bound by different transcription factors [##REF##15343339##18##]. In order to rationalise these resources, we began by looking for biologically-verified data on transcription factors that might control the repression of the subset of 44 genes of Additional file ##SUPPL##0##1##. Since only limited information was publicly available, we did not pursue this part of the dataset further. In contrast, for the 146 induced genes in Additional file ##SUPPL##0##1##, our observation of the similarity between the transcriptional changes in the V5 to V5.TM6*P transition and those in the transition from glucose to ethanol growth lead us to analyse the roles of Hap4, Cat8 and Mig1 which have been previously been verified to be crucial to the expression of genes in the diauxic shift [##REF##11024040##16##,##REF##12537548##19##,##REF##8756637##20##]: Hap4 (in complex with Hap2, 3 and 5) [##REF##12537548##19##] and Cat8 (optionally in complex with Sip4 and Adr1) [##REF##14685767##21##] are transcriptional activators, whereas Mig1 (optionally with Mig2/3) is a transcriptional repressor that works in concert with Tup1-Cyc8 [##REF##7724528##22##].</p> <p>The Hap complex is comprised of proteins Hap2, 3, 4 and 5, and has separate subunits for DNA binding and transcriptional activation. DNA binding is mediated by Hap2, 3 and 5 [##REF##7828851##23##] while Hap4 functions as an activation domain. The expression of <italic>HAP4 </italic>is repressed by glucose and is responsible for regulation of Hap complex target genes. We began by scrutinising data for a <italic>HAP4 </italic>over-expression strain, in which 246 genes had their expression affected by at least a factor of 2 [##REF##12537548##19##], which is 5.3% of the genome. 56 of these genes were found in our dataset and are listed in Table ##TAB##1##2##. Only three of these (<italic>GLK1</italic>, <italic>CTR3 </italic>and <italic>ZRT1</italic>) are down-regulated in the <italic>HAP4 </italic>over-expression strain, whereas <italic>GLK1 </italic>and <italic>CTR3 </italic>are induced in the V5.TM6*P strain. Interestingly, a typical Hap complex binding site was not seen in the 500 nucleotide region upstream of the <italic>GLK1, CTR3 </italic>and <italic>ZRT1 </italic>start codon (where the first A in ATG is position 1) while all other genes in Table ##TAB##1##2## except <italic>AGX1</italic>, <italic>MNP1 </italic>and <italic>YNK1 </italic>were found to have at least one Hap complex binding site. These Hap complex-binding-site-containing genes comprise a significant subset (26% of our dataset of 190 genes in Additional file ##SUPPL##0##1##) with a verified biological dependence on Hap4 activation. Clearly, then, Hap4 has a statistically-significant role to play in the metabolism of V5.TM6*P as there is an association between Hap4 dependence and our dataset (χ<sup>2 </sup>= 182.75, degrees of freedom = 1, p = 0.00).</p> <p>Since six of the genes from the <italic>HAP4 </italic>over-expression study did not contain Hap complex binding sites, we did a further <italic>in silico </italic>analysis of the induced genes of Additional file ##SUPPL##0##1##. We noted that the genes of Table ##TAB##1##2##, contained motifs that could be sub-divided into 2 distinct families: <bold>CCAAT</bold>G (which we denoted as a Hap4_1 site) and (G/<bold>C)CAA</bold>(G/<bold>T</bold>)CAA (a Hap4_2 site). The bold sequences are the conserved sequences previously identified in the <italic>HAP4 </italic>over-expression strain study [##REF##12537548##19##]. The genes of Additional file ##SUPPL##0##1## were then examined 500 nucleotides upstream and 200 nucleotides downstream from the start codon using WebMOTIF [##REF##17584794##24##] for these two biologically-relevant sites since it has been previously demonstrated that Hap complex binding sites are statistically over-represented up to 400 nucleotides upstream from the start codon [##REF##12537548##19##]: we found that 112/146 (77%) induced genes and 127/190 (67%) of the complete dataset in Additional file ##SUPPL##0##1## had one or more Hap complex binding site (Additional files ##SUPPL##0##1## and ##SUPPL##3##4##). Our analysis included presumed binding sites such as CAAATC and CCAAAC which arose from the biologically-verified data on transcription factors (Table ##TAB##1##2##). For example, <italic>AQY2 </italic>is identified as Hap complex-dependent, and has two CCAAAC sites, but no CCAATNA sites. There are only two other genes that do not contain CCAATNA sites in Table ##TAB##1##2##, namely <italic>CRC1 </italic>and <italic>ALD6</italic>. In addition it was notable that Hap complex consensus sequences are predominantly found upstream of mitochondrial genes with 84/127 (66%) of such genes having at least a partial mitochondrial location. Furthermore, of the induced genes with a known cellular location and at least two Hap complex binding sites, 69/80 (86%) are localized to the mitochondria. It appears, then, that Hap4 activation is preferentially directed at mitochondrially-related functions especially within the TCA-cycle and respiratory chain, but cytosolic support via the glyoxylate cycle is also potentially regulated by Hap4 activation (Additional file ##SUPPL##0##1##). It was further noted that multiple Hap complex binding sites are associated with induction (Fig. ##FIG##2##3A##). A Fisher's exact test showed that there was a statistically-significant association between the change in gene expression and the number of Hap complex binding sites in the gene (p = 0.00). This supported the fact that genes with two or more Hap complex binding sites are likely to be induced, whereas those with one site may be either induced or repressed.</p> <p>In order to adopt a similar approach for Cat8, we examined available biological data for Cat8 and Adr1 (as they are co-regulators [##REF##15743812##25##]) from chromatin immunoprecipitations [##REF##15743812##25##], and mRNA expression ratios in wt/Δcat8 strains [##REF##11024040##16##]. We identified 20 genes out of 66 from these studies that were also in our dataset (Table ##TAB##2##3##). Only <italic>ATO2 </italic>was oppositely regulated when comparing our dataset. For all other genes, the change in expression followed the same trend of up-regulation (Table ##TAB##2##3##). As for the Hap4 analysis, we examined these 20 genes for Cat8 binding sites, and identified a typical Cat8 motif [##REF##14685767##21##] in 13 of them, nine being of the Cat8-Sip4 type, TCCATTSRTCCGR (Additional file ##SUPPL##3##4##). The remaining six induced genes (<italic>CYB2</italic>, <italic>GUT1</italic>, <italic>GUT2</italic>, <italic>LSC2</italic>, <italic>SUE1 </italic>and <italic>ODC1</italic>) were further scrutinised, and it was apparent that all except <italic>ODC1 </italic>had been identified from immunoprecipitation data for Adr1 alone, consistent with their lack of Cat8-Sip4 motif. The induced genes of Additional file ##SUPPL##0##1## were therefore searched from the start codon to 1,500 nucleotides upstream for the generic Cat8 motif, CC------CCG motif giving 54 genes (Additional file ##SUPPL##3##4##). The non-coding intergenic regions of these genes were then searched using WebMOTIF [##REF##17584794##24##] and <italic>HXK1</italic>, <italic>COX4</italic>, <italic>ODC1</italic>, <italic>OM14</italic>, <italic>OM45 </italic>and <italic>YOR019w </italic>were found to share the consensus sequence GCCSSTSS(W/Y))CMS in common with <italic>IDP2 </italic>and are thus also listed in Additional file ##SUPPL##3##4##. We noted that genes with these Cat8 binding sites are induced in V5.TM6*P by an average factor of 6.2 while the average for induced genes in Additional file ##SUPPL##0##1## is 3.3. The highly-induced, Cat8-site-containing genes include <italic>HXK1 </italic>(factor change 9.0), <italic>PCK1 </italic>(6.9), <italic>STL1 </italic>(21.7), <italic>JEN1 </italic>(10.3), <italic>ODC1 </italic>(6.3) and <italic>SFC1 </italic>(18.3).</p> <p>We observed that some strongly-induced genes in Additional file ##SUPPL##0##1## did not have Hap complex or Cat8 binding sites. Indeed, the expression of many genes such as <italic>SUC2 </italic>is already known to be repressed by both Mig1 and Mig2 in glucose media, although Mig1 is the main factor in this response [##REF##8756637##20##,##REF##9832517##26##]. We therefore examined biological data on gene expression in <italic>Δmig1 </italic>and <italic>Δmig2 </italic>strains [##REF##8756637##20##,##REF##9832517##26##]. In these studies, 11 genes were presented as being Mig dependent, of which 5 were found in Additional file ##SUPPL##0##1## (<italic>HXK1</italic>, <italic>EMI2</italic>, <italic>HXT13, HXT15 </italic>and <italic>DSF1</italic>) and 3 in Additional file ##SUPPL##2##3## (<italic>REG2</italic>, <italic>DOG2 </italic>and <italic>YLR042C</italic>). We then examined protein binding microarray data [##REF##15543148##27##] and tabulated the combined outputs as Table ##TAB##3##4##, which lists 15 Mig-dependent genes. When searching the induced set in Additional file ##SUPPL##0##1##<italic> in silico </italic>for genes containing Mig1 binding sites, 31 were identified, predominantly containing consensus sequences of the SUC2B type (CCCCGGAT), but also in some cases sites that were more homologous to SUC2A (sharing the consensus motif CCCC(G/A)(G/C)AT [##REF##8114729##28##]). We noted that the presence of a Mig1 site in a gene correlated with high induction (Fig. ##FIG##2##3B##) and that overall, our analysis provides a highly complete description of how three transcriptions factors might regulate the 146 induced genes of Additional file ##SUPPL##0##1## (Fig. ##FIG##2##3C##).</p> <p>We noted that <italic>HAP4 </italic>(factor change 1.7, Additional file ##SUPPL##2##3##) is induced and <italic>CAT8 </italic>(0.5, adjusted p value &lt; 0.06; Additional file ##SUPPL##4##5##), <italic>MIG2 </italic>(0.5, Additional file ##SUPPL##0##1##) and <italic>MIG3 </italic>(0.3, Additional file ##SUPPL##0##1##) are down-regulated in V5.TM6*P compared to V5. In our dataset, the specific change in array signal for <italic>MIG1 </italic>had an adjusted p-value &gt; 0.05 over the complete experimental range. The result for <italic>CAT8 </italic>was unexpected as it contains an upstream Mig1 binding site similar to <italic>HAP4 </italic>(which is induced) and Cat8 target genes such as <italic>FBP1 </italic>and <italic>PCK1 </italic>are themselves induced. It is possible that the transcript was not probed accurately on the array either due to cross-hybridisation between closely related DNA sequences or on account of its low base expression value (0.24). A future real time Q-PCR experiment would probe the individual transcript.</p> <title>Transcripts dependent on external glucose concentration are not found in the respiratory V5.TM6*P strain, even though several glucose- or carbon-source-associated transcripts are glucose dependent in the parental wild-type strain</title> <p>The transcriptome of the respiratory V5.TM6*P strain was not found to be responsive to changes in glucose concentration in the culture medium. However, in the wild-type parent strain (V5) the transcriptome was found to vary with glucose availability according to a least squares fit of the logarithmic array values for RNA extracted at the six different glucose concentrations. This dependence on external glucose concentrations is seen in the V5 strain for 20 genes (including those without functional annotations) with an adjusted p-value &lt; 0.05 (Fig. ##FIG##3##4##). The genes with the strongest dependence on external glucose concentrations were found to be <italic>HXK1</italic>, <italic>RGS2</italic>, <italic>ADH7</italic>, <italic>CHA4 </italic>and <italic>GLK1 </italic>(Fig. ##FIG##3##4##). The sugar kinases <italic>HXK1 </italic>and <italic>GLK1 </italic>are induced as expected as they have their maximal expression during growth on other carbon sources [##REF##11024040##16##]. The gene product of <italic>RGS2 </italic>inhibits Gpa2 in the PKA pathway, which is one of the pathways induced by glucose [##REF##12694616##29##]. <italic>ADH7 </italic>appears to be induced on carbon sources other than glucose [##REF##12423374##30##]. In a study by DeRisi and colleagues, batch yeast cultures were harvested for cDNA array analysis at gradually decreasing glucose concentration [##REF##9381177##5##]. Our examination of the supplementary data provided by the authors showed that as glucose was consumed between 18.7 and 7.5 g/L, three of the genes in Figure ##FIG##3##4## (<italic>GLK1</italic>, <italic>HXK1</italic>, and <italic>FAL1</italic>) were found in the top 22 genes with a changed expression from that study. Overall, the results in Figure ##FIG##3##4## show that glucose repression is gradually relieved on going from 37 g/L to 9 g/L glucose. Furthermore relief of repression is already apparent at a surprisingly high glucose concentration, in contrast to the assumption that it is typically triggered over a range of low glucose concentrations [##REF##12694616##29##], the trigger being very low for some genes [##REF##8793872##31##].</p> <p>In order to examine whether the glucose dependence of genes in V5 affected our data in Additional file ##SUPPL##0##1##, we excluded the data points at 13.5 and 5.6 g/L glucose for V5.TM6*P, and 10.5 and 7.4 g/L glucose for V5 (which gave the highest array signals) and compared the cDNA array data in this case with each other exactly as for the full data set. Following this procedure, only a small subset of genes that are highly induced had their factor changes increased substantially, but this did not affect the fact that they remain within the top 30 most induced genes in Additional file ##SUPPL##0##1##. They are <italic>HXK1 </italic>(9.0 in Additional file ##SUPPL##0##1##, 19.3 following re-calculation; the gene contains binding sites for the Hap complex, Cat8 and Mig1); <italic>STL1 </italic>(21.7, 27.5; Hap complex, Cat8); <italic>DSF1 </italic>(9.4, 15.1; Hap complex, Mig1); <italic>CYB2 </italic>(10.1, 13.4; Hap complex, Mig1) and <italic>FBP1 </italic>(5.3, 7.4; Cat8, Mig1). Fourteen of the seventeen genes that were most affected by this recalculation were found to contain a Mig1 binding site. It is therefore possible that Mig1 is also involved in relieving glucose repression in the 10 to 5 g/L range, and not only at the depletion of glucose at the diauxic shift. Overall, we found that the amplitudes of the factor changes varied by only 9.3% in this new calculation compared with the original one, which indicated that the glucose dependence of the V5 strain does not have a major influence on the way we generated the data presented in Additional file ##SUPPL##0##1##.</p>
<title>Discussion</title> <p>The main physiological change in both the diauxic shift and the transition from a respiro-fermentative to a respiratory metabolism is that glycolysis is channelled to support respiration: this occurs upon glucose depletion during the diauxic shift and throughout the growth curve in glucose medium in V5.TM6*P. Although previously-published genomic analyses of yeast's response to the diauxic shift transition have all been sampled under different conditions – rich (YPD) medium with 2% glucose [##REF##9381177##5##], minimal defined medium in chemostat cultures at different dilution rates [##REF##15758028##4##], and minimal defined media in batch cultures sampled 4 h after the diauxic shift [##REF##15987902##3##] – common themes emerge with changes on going from V5 to V5.TM6*P. Overall 81% of genes which are induced or repressed on going from V5 to V5.TM6*P in our study were found in at least one of these three datasets for the diauxic shift. Furthermore, we also find that for all proteins induced in a study of the diauxic shift by 2D-gel electrophoresis [##REF##14597615##6##], the corresponding genes are induced on going from V5 to V5.TM6*P except for the alcohol dehydrogenase, Adh2. There is a hence close agreement between induction of the following genes and the levels of the corresponding proteins: <italic>FBP1</italic>, <italic>ICL1</italic>, <italic>PCK1</italic>, <italic>IDP2</italic>, <italic>MLS1</italic>, <italic>DLD1</italic>, <italic>ALD6</italic>, <italic>SDH1 </italic>and <italic>CIT2</italic>. This strong correlation with previously-published data and our own observed correlation of data for Cdc19 and <italic>CDC19 </italic>supports the phenotypic similarity between the V5 to V5.TM6*P transition and the diauxic shift. It is worth noting that there is no clear correlation of our data in Additional file ##SUPPL##0##1## with any possible differences in growth rate between the two strains. In a study by Regenberg and colleagues [##REF##17105650##32##], 180 genes changed in expression by more than a factor of 4 when comparing growth rates. Of these 180, 3% of all genes are found in Additional file ##SUPPL##0##1## providing no significant overlap in data.</p> <p>It is also possible to see parallels between our data set and a study of nitrogen deprivation and stationary phase growth for wild-type yeast by Gasch and colleagues [##REF##11102521##33##]. This is supported by the fact that ten of the 15 genes that were found to be strongly induced upon nitrogen deprivation and growth in stationary phase in that study are also induced in the V5 to V5.TM6*P transition. Only two genes, <italic>YGR067c </italic>(likely to encode a transcription factor as it has high partial sequence homology with <italic>MIG3 </italic>and <italic>ADR1 </italic>and it fine-tunes the response to nutrient limitation) and <italic>ECM13</italic>, are not induced in the V5.TM6*P strain, but were highly induced in the Gasch study [##REF##11102521##33##] upon nitrogen deprivation and stationary phase growth. In addition, we find that the genes encoding Mep1 and Ato2, which are involved in ammonia transport are also strongly induced in V5.TM6*P in our study. Such an induction is seen in both the diauxic shift studies performed in defined media [##REF##15758028##4##,##REF##11406594##34##], but not in rich medium [##REF##9381177##5##], suggestive of nutrient limitation. In fact, V5.TM6*P also appears carbon limited, not least as a result of the fact that amongst the most strongly induced genes are <italic>FBP1 </italic>and <italic>PCK1</italic>, which encode gluconeogenic enzymes. <italic>STL1</italic>, <italic>MAL31 </italic>and <italic>SUC2 </italic>(encoding proteins involved in usage and/or transport of alternative carbon sources: glycerol, maltose and sucrose respectively) are induced in V5.TM6*P, also consistent with its low glycolytic flux and suggestive of an adaptation of metabolism ('discontinuous modelling) as seen by Brauer and colleagues [##REF##15758028##4##,##REF##11406594##34##] and supported by our previous work [##REF##15071495##1##]. Several of these genes including <italic>FBP1 </italic>and <italic>STL1</italic>, as well as <italic>JEN1</italic>, <italic>ODC1, SFC1 </italic>and <italic>ATO2 </italic>have already been identified as being either induced under carbon limitation or being specific indicators of carbon limitation in yeast cultures [##REF##12414795##13##], consistent with our own data.</p> <p>The analysis of transcription factor binding sites present in the genes of Additional file ##SUPPL##0##1## revealed that 127 genes had Hap complex binding sites, the majority of which encode mitochondrial proteins. Of genes that only had Hap complex binding sites present, thirteen had only one and of these only five were induced leading us to hypothesise that multiple Hap complex binding sites are correlated with induction of the Additional file ##SUPPL##0##1## genes (Fig. ##FIG##2##3A##). Interestingly, in a microarray screen of a <italic>hap1Δ </italic>mutant for Hap1-dependent genes [##REF##12112237##35##], 7/24 genes that were induced in the <italic>hap1Δ </italic>mutant are also induced in Additional file ##SUPPL##0##1## (<italic>CYB2</italic>, <italic>MLS1</italic>, <italic>HMX1</italic>, <italic>PUT4</italic>, <italic>STL1</italic>, <italic>HPA2</italic>, <italic>PUT1</italic>), and in each case also have a Hap complex binding site. In addition, 19 genes containing a Cat8 site and 31 genes containing a Mig1 site were predicted, with many of these genes being experimentally-verified in previous <italic>in vivo </italic>studies. Our analysis of the regulation of 88% of the induced genes of our unique respiratory strain (Fig. ##FIG##2##3##) is further supported by a novel computational method [##REF##17224918##36##] that predicts the Hap complex and Mig1 control the genes of Additional file ##SUPPL##0##1##; Cat8 is not included in that database.</p> <p>A quite unexpected result from comparing the transcriptional profiles of V5 with V5.TM6*P is the data set in Figure ##FIG##3##4##, showing genes from the parent strain that are clearly glucose dependent and that the glucose dependence observed begins at concentrations (6 – 36 g/L glucose) higher than the range generally regarded as the glucose-repressed region. Strain DBY7286 [##REF##9381177##5##], which is unrelated to V5, previously showed a gradual relief of glucose repression starting at 20 g/L glucose for <italic>GLK1</italic>, <italic>HXK1 </italic>and <italic>FAL1</italic>. These results for V5 and DBY7286 suggest that gradual relief is a general <italic>S. cerevisiae </italic>phenomenon. Future studies will show whether the gradient of the slope is different for different strains. That such glucose dependence requires an intact glucose uptake system is suggested by the fact that it is not seen in the V5.TM6*P strain: the array signal is not changed significantly for any of the genes in the V5.TM6*P strain at the different glucose concentrations sampled. This implies that in V5.TM6*P, glucose sensing is independent of external glucose concentration. The wild-type yeast, on the other hand, shows an extensive modulation of its expression profile at high glucose concentrations despite the fact that its physiological response is largely unaffected. The insensitivity of the V5.TM6*P strain therefore makes it an extremely valuable biotechnological tool as it can be cultured in a wide range of external glucose concentrations, whilst maintaining the same respiratory phenotype [##REF##16204537##2##].</p>
<title>Conclusion</title> <p>In this study, we have been able to characterize the transcriptome of a unique respiratory yeast strain. We have been able to identify highly complete collections of known genes in the TCA cycle, glyoxylate cycle and respiratory chain that are consistent with a respiratory metabolism. Our results suggest that there has been genetic remodelling predominantly through the activity of Hap4, Cat8 and Mig1, and that the gene expression profile of V5.TM6*P during growth on glucose resembles wild-type <italic>S. cerevisiae </italic>cells in the diauxic shift.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p>We previously described the first respiratory <italic>Saccharomyces cerevisiae </italic>strain, KOY.TM6*P, by integrating the gene encoding a chimeric hexose transporter, Tm6*, into the genome of an <italic>hxt </italic>null yeast. Subsequently we transferred this respiratory phenotype in the presence of up to 50 g/L glucose to a yeast strain, V5 <italic>hxt1-7</italic>Δ, in which only <italic>HXT1-7 </italic>had been deleted. In this study, we compared the transcriptome of the resultant strain, V5.TM6*P, with that of its wild-type parent, V5, at different glucose concentrations.</p> <title>Results</title> <p>cDNA array analyses revealed that alterations in gene expression that occur when transitioning from a respiro-fermentative (V5) to a respiratory (V5.TM6*P) strain, are very similar to those in cells undergoing a diauxic shift. We also undertook an analysis of transcription factor binding sites in our dataset by examining previously-published biological data for Hap4 (in complex with Hap2, 3, 5), Cat8 and Mig1, and used this in combination with verified binding consensus sequences to identify genes likely to be regulated by one or more of these. Of the induced genes in our dataset, 77% had binding sites for the Hap complex, with 72% having at least two. In addition, 13% were found to have a binding site for Cat8 and 21% had a binding site for Mig1. Unexpectedly, both the up- and down-regulation of many of the genes in our dataset had a clear glucose dependence in the parent V5 strain that was not present in V5.TM6*P. This indicates that the relief of glucose repression is already operable at much higher glucose concentrations than is widely accepted and suggests that glucose sensing might occur inside the cell.</p> <title>Conclusion</title> <p>Our dataset gives a remarkably complete view of the involvement of genes in the TCA cycle, glyoxylate cycle and respiratory chain in the expression of the phenotype of V5.TM6*P. Furthermore, 88% of the transcriptional response of the induced genes in our dataset can be related to the potential activities of just three proteins: Hap4, Cat8 and Mig1. Overall, our data support genetic remodelling in V5.TM6*P consistent with a respiratory metabolism which is insensitive to external glucose concentrations.</p>
<title>Authors' contributions</title> <p>NB initiated the study, prepared the RNA and participated in the data analysis and interpretation, particularly of the transcription factor binding sites. CF carried out the fermentations. PM performed the array analysis. MW performed the statistical analyses of the array data. CC and LS generated the arrays and performed the arraying. LG and CL participated in the experimental design and helped to draft the manuscript. RB participated in the study design, coordinated the data analysis and interpretation, and drafted the manuscript. All authors contributed to the final version of the manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was supported by the European Commission via contracts LSHG-CT-2004-504601 (E-MeP) and LSHG-CT-2006-037793 (OptiCryst) to RMB, and QLG2-CT-2002-00988 (SPINE) to LG and RMB. The Wistar Genomics facility is supported in part by NCI P30 CA10815-34S3. The BBSRC supports bioreactors in the RMB laboratory through an REI award.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>An overview of our dataset in the context of key yeast metabolic pathways</bold>. Genes with altered expression in V5.TM6*P compared with V5 are marked as follows: induced genes are denoted in large, bold font; down-regulated genes are crossed out. Genes that were previously found to be altered in studies of respiratory-state yeast [##REF##14597615##6##], but which are not found in our dataset are denoted in small, grey text. The chimeric Tm6*p transporter, which comprises the amino-terminal half of Hxt1 and the carboxy-terminal half of Hxt7, is responsible for the respiratory phenotype of V5.TM6*P at high glucose concentrations and is present in the plasma membrane. The grey boundaries represent the plasma membrane (top) and the mitochondrial membrane (left). <sup>†</sup>Note that <italic>ICL</italic>(1.5), <italic>NDE1 </italic>(1.8) and <italic>GPD1 </italic>(1.8) are induced in our study (see Additional file ##SUPPL##2##3##). <sup>‡ </sup><italic>LPD1 </italic>is a component of pyruvate dehydrogenase.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>An overview of our dataset in the context of the respiratory chain in yeast</bold>. Genes with altered expression in V5.TM6*P compared with V5 are marked as follows: genes induced by a factor ≥2.0 are denoted in large, bold font; genes induced by a factor &lt;2.0 are denoted in smaller, grey font.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>(A) The influence of multiple Hap complex sites on the change in gene expression of the genes in Additional file </bold>##SUPPL##0##1##. A Fisher's exact test was performed for genes containing predicted binding sites for the Hap complex and their association with induction or repression. Genes were grouped into either those with one Hap complex binding site or those with two or more according to Additional file ##SUPPL##0##1##. <bold>(B) The influence of Hap4, Cat8 and Mig1 on the magnitude of induction of the genes in Additional file </bold>##SUPPL##0##1##. The average fold-induction of genes in Additional file ##SUPPL##0##1##, grouped according to the transcription factor binding sites they contain, was calculated. It was observed that those containing binding sites for each of the Hap complex (where Hap4 is the activator), Cat8 and Mig1, had the highest average fold induction with Mig1 being a dominant factor in high induction. A bar for Cat8 is not included as only one gene (<italic>YOR019W</italic>, factor change 2.1) has a Cat8 site alone. <bold>(C) A Venn diagram of the distribution of binding sites for the Hap complex, Cat8 and Mig1 in the genes in Additional file </bold>##SUPPL##0##1##. 88% of the genes in Additional file ##SUPPL##0##1## contain a binding site for one or more of the Hap complex, Cat8 and Mig1. The range of induction for the 190 genes from Additional file ##SUPPL##0##1## is 2.0 to 21.7.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Transcript dependence on external glucose concentration for 20 genes in the V5 strain</bold>. A clear glucose-dependence (grey bars) was observed for all genes in the V5 strain with only a small variation in array signal (white bars) for the same genes in V5.TM6*P. Genes marked with a dot also have a significant change of expression on going from V5 to V5.TM6*P and as such are listed in Additional file ##SUPPL##0##1##. The factor change in gene expression was obtained as the quotient between the expression values at 37 and 9 g/L glucose given by the line fitted to the expression data for each gene. R<sup>2 </sup>values are for the linear fit for the V5 strain. The R<sup>2 </sup>values for V5.TM6*P were between 0.16 and 0.78 indicating no significant linear relationship.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Growth characteristics of V5 and V5.TM6*P at 30°C starting at 50 g/L glucose.</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left"><bold>V5</bold></td><td align="left"><bold>V5.TM6*P</bold></td></tr></thead><tbody><tr><td align="left"><bold>Generation time during glucose growth (h)</bold></td><td align="left">3.7 (0.2; n = 2)</td><td align="left">4.5 (0.3; n = 4)</td></tr><tr><td align="left"><bold>Biomass yield (g/g)</bold></td><td align="left">0.21 (0.02 ; n = 2)</td><td align="left">0.32 (0.01; n = 4)</td></tr><tr><td align="left"><bold>Glucose consumption (mmolg dry weight)<sup>-1</sup>h<sup>-1</sup></bold></td><td align="left">9.00 (0.70; n = 2)</td><td align="left">2.70 (0.20; n = 4)</td></tr><tr><td align="left"><bold>Ethanol yield (g/g)</bold></td><td align="left">0.33 (0.01; n = 2)</td><td align="left">0.05 (0.01; n = 4)</td></tr><tr><td align="left"><bold>Glycerol yield (g/g)</bold></td><td align="left">0.04 (0; n = 2)</td><td align="left">0</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Biologically-verified data on transcription factors that act on genes also found in Additional file ##SUPPL##0##1##: transcript dependence on <italic>HAP4 </italic>over-expression.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>ORF</bold></td><td align="left"><bold>GENE</bold></td><td align="left"><bold>FACTOR CHANGE IN GENE EXPRESSION FROM V5 TO V5.TM6*P</bold></td><td align="left"><bold>FACTOR CHANGE IN GENE EXPRESSION FROM WT TO A <italic>HAP4 </italic>OVER-EXPRESSION STRAIN </bold>[##REF##12537548##19##]</td></tr></thead><tbody><tr><td align="left" colspan="4"><bold>Glycolysis</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YCL040W</italic></td><td align="left"><italic>GLK1</italic></td><td align="left">3.8</td><td align="left">0.4</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Alternative carbon source utilization</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YML054C</italic></td><td align="left"><italic>CYB2</italic></td><td align="left">10.1</td><td align="left">5.0</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>TCA cycle</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YIL125W</italic></td><td align="left"><italic>KGD1</italic></td><td align="left">2.2</td><td align="left">2.2</td></tr><tr><td align="left"><italic>YDR148C</italic></td><td align="left"><italic>KGD2</italic></td><td align="left">3.5</td><td align="left">4.4</td></tr><tr><td align="left"><italic>YFL018C</italic></td><td align="left"><italic>LPD1</italic></td><td align="left">2.0</td><td align="left">2.2</td></tr><tr><td align="left"><italic>YKL085W</italic></td><td align="left"><italic>MDH1</italic></td><td align="left">3.6</td><td align="left">2.2</td></tr><tr><td align="left"><italic>YKL148C</italic></td><td align="left"><italic>SDH1</italic></td><td align="left">4.6</td><td align="left">4.7</td></tr><tr><td align="left"><italic>YLL041C</italic></td><td align="left"><italic>SDH2</italic></td><td align="left">4.1</td><td align="left">3.9</td></tr><tr><td align="left"><italic>YKL141W</italic></td><td align="left"><italic>SDH3</italic></td><td align="left">3.7</td><td align="left">2.3</td></tr><tr><td align="left"><italic>YDR178W</italic></td><td align="left"><italic>SDH4</italic></td><td align="left">5.3</td><td align="left">3.2</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Glyoxylate cycle</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YIR029W</italic></td><td align="left"><italic>DAL2</italic></td><td align="left">3.0</td><td align="left">3.0</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Respiratory chain</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YPL078C</italic></td><td align="left"><italic>ATP4</italic></td><td align="left">2.1</td><td align="left">2.0</td></tr><tr><td align="left"><italic>YKL016C</italic></td><td align="left"><italic>ATP7</italic></td><td align="left">2.3</td><td align="left">2.4</td></tr><tr><td align="left"><italic>YPR020W</italic></td><td align="left"><italic>ATP20</italic></td><td align="left">2.5</td><td align="left">3.4</td></tr><tr><td align="left"><italic>YGR174C</italic></td><td align="left"><italic>CBP4</italic></td><td align="left">2.8</td><td align="left">2.1</td></tr><tr><td align="left"><italic>YBL045C</italic></td><td align="left"><italic>COR1</italic></td><td align="left">2.0</td><td align="left">2.0</td></tr><tr><td align="left"><italic>YGL187C</italic></td><td align="left"><italic>COX4</italic></td><td align="left">2.0</td><td align="left">2.3</td></tr><tr><td align="left"><italic>YNL052W</italic></td><td align="left"><italic>COX5A</italic></td><td align="left">2.0</td><td align="left">2.1</td></tr><tr><td align="left"><italic>YHR051W</italic></td><td align="left"><italic>COX6</italic></td><td align="left">2.1</td><td align="left">2.2</td></tr><tr><td align="left"><italic>YLL009C</italic></td><td align="left"><italic>COX17</italic></td><td align="left">2.8</td><td align="left">4.7</td></tr><tr><td align="left"><italic>YOR065W</italic></td><td align="left"><italic>CYT1</italic></td><td align="left">2.5</td><td align="left">2.3</td></tr><tr><td align="left"><italic>YKL087C</italic></td><td align="left"><italic>CYT2</italic></td><td align="left">2.2</td><td align="left">2.8</td></tr><tr><td align="left"><italic>YIL098C</italic></td><td align="left"><italic>FMC1</italic></td><td align="left">2.5</td><td align="left">2.3</td></tr><tr><td align="left"><italic>YDL181W</italic></td><td align="left"><italic>INH1</italic></td><td align="left">2.2</td><td align="left">5.5</td></tr><tr><td align="left"><italic>YBR185C</italic></td><td align="left"><italic>MBA1</italic></td><td align="left">2.3</td><td align="left">2.3</td></tr><tr><td align="left"><italic>YML120C</italic></td><td align="left"><italic>NDI1</italic></td><td align="left">4.1</td><td align="left">4.4</td></tr><tr><td align="left"><italic>YEL024W</italic></td><td align="left"><italic>RIP1</italic></td><td align="left">2.3</td><td align="left">2.4</td></tr><tr><td align="left"><italic>YPR151C</italic></td><td align="left"><italic>SUE1</italic></td><td align="left">2.2</td><td align="left">3.0</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Plasma membrane transport</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YLL052C</italic></td><td align="left"><italic>AQY2</italic></td><td align="left">3.0</td><td align="left">3.0</td></tr><tr><td align="left"><italic>YLR411W</italic></td><td align="left"><italic>CTR3</italic></td><td align="left">3.0</td><td align="left">0.5</td></tr><tr><td align="left"><italic>YGL255W</italic></td><td align="left"><italic>ZRT1</italic></td><td align="left">0.5</td><td align="left">0.1</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Mitochondrial transport</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YNR002C</italic></td><td align="left"><italic>ATO2</italic></td><td align="left">3.8</td><td align="left">3.8</td></tr><tr><td align="left"><italic>YKL217W</italic></td><td align="left"><italic>JEN1</italic></td><td align="left">10.3</td><td align="left">2.9</td></tr><tr><td align="left"><italic>YPL134C</italic></td><td align="left"><italic>ODC1</italic></td><td align="left">6.3</td><td align="left">3.3</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Transcriptional regulation</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YOL071W</italic></td><td align="left"><italic>EMI5</italic></td><td align="left">2.1</td><td align="left">2.9</td></tr><tr><td align="left"><italic>YNL333W</italic></td><td align="left"><italic>SNZ2</italic></td><td align="left">3.0</td><td align="left">6.4</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Fatty acid metabolism</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YKL150W</italic></td><td align="left"><italic>MCR1</italic></td><td align="left">3.1</td><td align="left">2.8</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Biosynthesis</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YFL030W</italic></td><td align="left"><italic>AGX1</italic></td><td align="left">3.6</td><td align="left">7.8</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Ribosomal proteins in the mitochondria or cytosol</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YDR296W</italic></td><td align="left"><italic>MHR1</italic></td><td align="left">2.0</td><td align="left">2.5</td></tr><tr><td align="left"><italic>YGL068W</italic></td><td align="left"><italic>MNP1</italic></td><td align="left">2.1</td><td align="left">2.4</td></tr><tr><td align="left"><italic>YDR116C</italic></td><td align="left"><italic>MRPL1</italic></td><td align="left">2.6</td><td align="left">3.5</td></tr><tr><td align="left"><italic>YKL138C</italic></td><td align="left"><italic>MRPL31</italic></td><td align="left">2.2</td><td align="left">2.8</td></tr><tr><td align="left"><italic>YGR165W</italic></td><td align="left"><italic>MRPS35</italic></td><td align="left">2.0</td><td align="left">2.7</td></tr><tr><td align="left"><italic>YHR038W</italic></td><td align="left"><italic>RRF1 (FIL1)</italic></td><td align="left">2.2</td><td align="left">4.8</td></tr><tr><td align="left"><italic>YDR041W</italic></td><td align="left"><italic>RSM10</italic></td><td align="left">2.1</td><td align="left">2.4</td></tr><tr><td align="left"><italic>YDR175C</italic></td><td align="left"><italic>RSM24</italic></td><td align="left">2.4</td><td align="left">2.9</td></tr><tr><td align="left"><italic>YFR049W</italic></td><td align="left"><italic>YMR31</italic></td><td align="left">2.5</td><td align="left">2.6</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Other</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YBR262C</italic></td><td align="left"><italic>AIM5</italic></td><td align="left">2.3</td><td align="left">2.1</td></tr><tr><td align="left"><italic>YFR011C</italic></td><td align="left"><italic>AIM13</italic></td><td align="left">2.1</td><td align="left">3.4</td></tr><tr><td align="left"><italic>YLR168C</italic></td><td align="left"><italic>AIM30</italic></td><td align="left">2.6</td><td align="left">3.3</td></tr><tr><td align="left"><italic>YML087C</italic></td><td align="left"><italic>AIM33</italic></td><td align="left">4.2</td><td align="left">11.6</td></tr><tr><td align="left"><italic>YBR230C</italic></td><td align="left"><italic>OM14</italic></td><td align="left">3.3</td><td align="left">2.4</td></tr><tr><td align="left"><italic>YDL104C</italic></td><td align="left"><italic>QRI7</italic></td><td align="left">2.1</td><td align="left">2.7</td></tr><tr><td align="left"><italic>YOR187W</italic></td><td align="left"><italic>TUF1</italic></td><td align="left">2.3</td><td align="left">2.4</td></tr><tr><td align="left"><italic>YKL067W</italic></td><td align="left"><italic>YNK1</italic></td><td align="left">3.1</td><td align="left">2.1</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Not characterized</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YGR110W</italic></td><td align="left">N/A</td><td align="left">2.3</td><td align="left">2.0</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Biologically-verified data on transcription factors that act on genes also found in Additional file ##SUPPL##0##1##: transcript dependence on Cat8 and Adr1.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>ORF</bold></td><td align="left"><bold>GENE</bold></td><td align="right"><bold>FACTOR CHANGE IN GENE EXPRESSION FROM V5 TO V5.TM6*P<sup>a</sup></bold></td><td align="right"><bold>FACTOR CHANGE IN GENE EXPRESSION USING DATA FROM CHROMATIN IMMUNOPRECIPITATION FOR Cat8<sup>b</sup></bold>[##REF##15743812##25##]<bold>, mRNA EXPRESSION RATIOS FOR </bold><bold>WT</bold><bold>/Δcat8<sup>c</sup></bold>[##REF##11024040##16##]<bold> AND IMMUNOPRECIPITATION DATA FOR Adr1<sup>d </sup></bold>[##REF##15743812##25##].</td></tr></thead><tbody><tr><td align="left" colspan="4"><bold>Gluconeogenesis</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YLR377C</italic></td><td align="left"><italic>FBP1</italic></td><td align="right">5.3</td><td align="right">5.0<sup>b</sup>; MAX<sup>c</sup></td></tr><tr><td align="left"><italic>YKR097W</italic></td><td align="left"><italic>PCK1</italic></td><td align="right">6.9</td><td align="right">7.7<sup>b</sup>; 151<sup>c</sup></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Alternative carbon source utilization</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YPL061W</italic></td><td align="left"><italic>ALD6</italic></td><td align="right">2.3</td><td align="right">5.2<sup>b</sup>; 8<sup>c</sup></td></tr><tr><td align="left"><italic>YML054C</italic></td><td align="left"><italic>CYB2</italic></td><td align="right">10.1</td><td align="right">7.4<sup>d</sup></td></tr><tr><td align="left"><italic>YHL032C</italic></td><td align="left"><italic>GUT1</italic></td><td align="right">6.2</td><td align="right">3.7<sup>d</sup></td></tr><tr><td align="left"><italic>YIL155C</italic></td><td align="left"><italic>GUT2</italic></td><td align="right">3.3</td><td align="right">4.6<sup>d</sup></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>TCA cycle</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YGR244C</italic></td><td align="left"><italic>LSC2</italic></td><td align="right">2.9</td><td align="right">5.2<sup>d</sup></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Glyoxylate cycle</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YCR005C</italic></td><td align="left"><italic>CIT2</italic></td><td align="right">2.2</td><td align="right">1.7<sup>c</sup></td></tr><tr><td align="left"><italic>YLR174W</italic></td><td align="left"><italic>IDP2</italic></td><td align="right">4</td><td align="right">2.9<sup>b</sup>; 23<sup>c</sup></td></tr><tr><td align="left"><italic>YOL126C</italic></td><td align="left"><italic>MDH2</italic></td><td align="right">3</td><td align="right">12.7<sup>b</sup>; 3.6<sup>c</sup></td></tr><tr><td align="left"><italic>YNL117W</italic></td><td align="left"><italic>MLS1</italic></td><td align="right">3.5</td><td align="right">8.9<sup>b</sup>; 79<sup>c</sup></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Respiratory chain</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YPR151C</italic></td><td align="left"><italic>SUE1</italic></td><td align="right">3.7</td><td align="right">4.4<sup>d</sup></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Plasma membrane transport</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YOR348C</italic></td><td align="left"><italic>PUT4</italic></td><td align="right">5.8</td><td align="right">2.7<sup>b</sup>; 6.1<sup>c</sup>; 2.0<sup>d</sup></td></tr><tr><td align="left"><italic>YDR536W</italic></td><td align="left"><italic>STL1</italic></td><td align="right">21.7</td><td align="right">4.4<sup>c</sup></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Mitochondrial transport</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YNR002C</italic></td><td align="left"><italic>ATO2</italic></td><td align="right">3.8</td><td align="right">0.33<sup>c</sup></td></tr><tr><td align="left"><italic>YML042W</italic></td><td align="left"><italic>CAT2</italic></td><td align="right">2.6</td><td align="right">2.8<sup>b</sup>; 2.4<sup>c</sup></td></tr><tr><td align="left"><italic>YOR100C</italic></td><td align="left"><italic>CRC1</italic></td><td align="right">2.2</td><td align="right">8.1<sup>c</sup></td></tr><tr><td align="left"><italic>YKL217W</italic></td><td align="left"><italic>JEN1</italic></td><td align="right">10.3</td><td align="right">3.6<sup>b</sup>; 10<sup>c</sup>; 10.5<sup>d</sup></td></tr><tr><td align="left"><italic>YPL134C</italic></td><td align="left"><italic>ODC1</italic></td><td align="right">6.3</td><td align="right">3.2<sup>b</sup></td></tr><tr><td align="left"><italic>YJR095W</italic></td><td align="left"><italic>SFC1</italic></td><td align="right">18.3</td><td align="right">6.2<sup>b</sup>; MAX<sup>c</sup></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T4"><label>Table 4</label><caption><p>Biologically-verified data on transcription factors that act on genes also found in Additional file ##SUPPL##0##1##: transcript dependence on Mig1.</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left"><bold>ORF</bold></td><td align="left"><bold>GENE</bold></td><td align="right"><bold>FACTOR CHANGE IN GENE EXPRESSION FROM V5 TO V5.TM6*P</bold><sup>a</sup></td><td align="right"><bold>Mig1 BINDING VERIFIED USING DATA FROM PROTEIN-BINDING MICROARRAY<sup>b </sup></bold>[##REF##15543148##27##]<bold> AND <italic>IN VIVO</italic><sup>c </sup></bold>[##REF##9832517##26##]<bold> ANALYSES</bold></td></tr></thead><tbody><tr><td align="left" colspan="4"><bold>Glycolysis</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YFR053C</italic></td><td align="left"><italic>HXK1</italic></td><td align="right">9</td><td align="right">b,c</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Gluconeogenesis</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YLR377C</italic></td><td align="left"><italic>FBP1</italic></td><td align="right">5.3</td><td align="right">b</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Alternative carbon source utilization</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YHL032C</italic></td><td align="left"><italic>GUT1</italic></td><td align="right">6.2</td><td align="right">b</td></tr><tr><td align="left"><italic>YIL162W</italic></td><td align="left"><italic>SUC2</italic></td><td align="right">6.8</td><td align="right">c</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Glyoxylate cycle</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YCR005C</italic></td><td align="left"><italic>CIT2</italic></td><td align="right">2.2</td><td align="right">b</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Plasma membrane transport</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YEL069C</italic></td><td align="left"><italic>HXT13</italic></td><td align="right">3.9</td><td align="right">b,c</td></tr><tr><td align="left"><italic>YDL245C</italic></td><td align="left"><italic>HXT15</italic></td><td align="right">3.2</td><td align="right">c</td></tr><tr><td align="left"><italic>YBR298C</italic></td><td align="left"><italic>MAL31</italic></td><td align="right">3.1</td><td align="right">b</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Mitochondrial transport</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YNR002C</italic></td><td align="left"><italic>ATO2</italic></td><td align="right">3.8</td><td align="right">b</td></tr><tr><td align="left"><italic>YKL217W</italic></td><td align="left"><italic>JEN1</italic></td><td align="right">10.3</td><td align="right">b</td></tr><tr><td align="left"><italic>YPL134C</italic></td><td align="left"><italic>ODC1</italic></td><td align="right">6.3</td><td align="right">b</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Transcriptional regulation</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YDR516C</italic></td><td align="left"><italic>EMI2</italic></td><td align="right">3.6</td><td align="right">b,c</td></tr><tr><td align="left"><italic>YEL066W</italic></td><td align="left"><italic>HPA3</italic></td><td align="right">2.2</td><td align="right">b</td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left" colspan="4"><bold>Other</bold></td></tr><tr><td colspan="4"><hr/></td></tr><tr><td align="left"><italic>YEL070W</italic></td><td align="left"><italic>DSF1</italic></td><td align="right">9.4</td><td align="right">b,c</td></tr><tr><td align="left"><italic>YGR243W</italic></td><td align="left"><italic>FMP43</italic></td><td align="right">4.9</td><td align="right">b</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p><bold>Transcriptome analysis of V5.TM6*P compared with the parental V5 strain</bold>. Genes were tabulated if their expression was changed on going from the V5 parental strain to V5.TM6*P, as described in 'Materials and Methods'. All genes were then sorted according to their roles in yeast cellular physiology and alphabetically by gene name under each sub-heading. Note that <italic>HXT1-7 </italic>are deleted in the V5.TM6*P with a part of <italic>HXT1 </italic>and <italic>HXT7 </italic>reincorporated as the <italic>TM6</italic>* chimera. Due to overall high homology within the hexose transporter family we still observed changes on the cDNA array (<italic>HXT1 </italic>(<italic>YHR094C</italic>); factor 0.4; <italic><underline>HXT4 </underline></italic>(<italic>YHR092C</italic>); 0.2; <italic>HXT6 </italic>(<italic>YDR343C</italic>); 2.1; these are not listed in Additional file ##SUPPL##1##2##). As YLL053C is annotated as continuous with YLL052C (AQY2) in some strains, these ORFs are counted once here as YLL052C. Genes that were dubious were not included in Additional file ##SUPPL##0##1## and neither were the Ty-transposable elements, <italic>YGR161C-C </italic>(encoding TyA gag protein) and <italic>YNL054W-B </italic>(encoding TyB gag protein). The top 30 genes with the largest numerical factor change (both up- and down-regulated) are underlined. If no functional information was available, phenotypic data from a deletion mutant was entered under 'Function'. In the 'Factor Change' column, change is expressed as a factor, where that factor is x when a gene expressed with intensity '1' in V5 is expressed with intensity 'x' in V5.TM6*P. All T-tests were jointly adjusted for multiple testing using the false discovery rate method (disregarding correlations between genes): p-values were adjusted so that when selecting all genes with p-values less than a threshold q, a proportion of q false positives would be expected amongst these genes. The genes shown have p-values &lt; 0.05, and thus the expected false discovery rate is 5%. Transcription factors are listed if a binding site was predicted as described in the text. Data on 'Copies/cell' are from the yeast GFP fusion localization database and were collected from wild-type yeast grown on glucose. 'Protein Localisation' was extracted from the <italic>Saccharomyces </italic>Genome Database on 14<sup>th </sup>February 2008 and in each case has been manually curated by the site's curators unless the entry is underlined, indicating that the localization has been extracted from a genome wide study.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>This table provides the information from Additional file ##SUPPL##0##1##, but in Excel format, with base expression values as well as SAGE [##REF##9008165##41##] and other literature data also included. In the 'Factor Change' column, change is expressed as a factor, where that factor is x when a gene expressed with intensity '1' in V5 is expressed with intensity 'x' in V5.TM6*P. In Additional file ##SUPPL##0##1## we have removed dubious genes and Ty-elements, but all genes that were originally spotted on the array are included in Additional file ##SUPPL##1##2##. As for Additional file ##SUPPL##0##1##, <italic>YLL052C </italic>and <italic>YLL053C </italic>are treated as one data point.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional file 3</title><p>A further 636 genes changed their expression levels either up by a factor of &lt;2.0 (317 genes) or down by a factor &gt; 0.5 (319 genes) and are presented in Excel format.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional file 4</title><p>Binding sites are listed for the Hap complex, Cat8 and Mig1. The file lists exact position and sequences of the relevant sites.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S5"><caption><title>Additional file 5</title><p>A computation in which the datasets with the two lowest glucose values have been removed (V5: 10.5, 7.4 g/L glucose and V5.TM6*P: 13.5, 5.6 g/L glucose). Genes are only listed if they are not present in Additional files ##SUPPL##0##1##, ##SUPPL##1##2##, or ##SUPPL##2##3##. Overall, the amplitudes of the factor changes varied by only 9.3% in this new calculation compared with the original one, which indicated that the glucose dependence of the V5 strain does not have a major influence on the way we generated the data presented in Additional file ##SUPPL##0##1##.</p></caption></supplementary-material>
<table-wrap-foot><p>Yields were calculated as g product (dry weight, ethanol or glucose, respectively)/g of glucose consumed/catabolised for V5 and V5.TM6*P, cultured as described in 'Materials and Methods'. The numbers in parentheses are the standard error of the mean, together with the number of replicates (n) for each experiment. The seemingly low biomass yield for V5 was measured after consumption of all glucose. Yield increases as all ethanol and glycerol are consumed during subsequent respiratory growth.</p></table-wrap-foot> <table-wrap-foot><p>Genes were tabulated and sorted according to their roles in yeast cellular physiology, as described in the legend for Additional file ##SUPPL##0##1##. Note that the seven hexose transporter genes, <italic>HXT1</italic>-<italic>7</italic>, are deleted in the V5.TM6*P with the amino-terminal half of <italic>HXT1 </italic>and the carboxy-terminal half of <italic>HXT7 </italic>reincorporated as the <italic>TM6* </italic>chimera. In the 'Factor Change' column, change is expressed as a factor, where that factor is x when a gene expressed with intensity '1' in V5 is expressed with intensity 'x' in V5.TM6*P. All T-tests were jointly adjusted for multiple testing using the false discovery rate method (disregarding correlations between genes): p-values were adjusted so that when selecting all genes with p-values less than a threshold q, a proportion of q false positives would be expected amongst these genes. The genes shown have p-values &lt; 0.05, and thus the expected false discovery rate is 5%. The factor in the final column is taken from Lascaris <italic>et al </italic>[##REF##12537548##19##]. <sup>†</sup>Note that due to overall high homology within the hexose transporter family, changes were still observed on the cDNA array for <italic>HXT4 </italic>(<italic>YHR092C</italic>) even though it is deleted as already mentioned.</p></table-wrap-foot> <table-wrap-foot><p>The table is arranged as for Table 2. Factor changes denoted '<sup>a</sup>' are from Additional file ##SUPPL##0##1##, '<sup>b</sup>' refer to chromatin immunoprecipitation data for Cat8 [##REF##15743812##25##], '<sup>c</sup>' refer to mRNA expression ratios for wt/<italic>Δcat8 </italic>[##REF##11024040##16##] and '<sup>d</sup>' refer to chromatin immunoprecipitation data for Adr1 [##REF##15743812##25##].</p></table-wrap-foot> <table-wrap-foot><p>The table is arranged as for Table 2. Factor changes denoted '<sup>a</sup>' are from Additional file ##SUPPL##0##1##. Mig1 binding was verified using data from protein-binding microarray [##REF##15543148##27##] studies denoted '<sup>b' </sup>and <italic>in vivo </italic>[##REF##9832517##26##] analyses denoted '<sup>c' </sup>although no factor changes were indicated in those studies.</p></table-wrap-foot>
<graphic xlink:href="1471-2164-9-365-1"/> <graphic xlink:href="1471-2164-9-365-2"/> <graphic xlink:href="1471-2164-9-365-3"/> <graphic xlink:href="1471-2164-9-365-4"/>
<media xlink:href="1471-2164-9-365-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-365-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-365-S3.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-365-S4.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1471-2164-9-365-S5.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Smyth", "Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W"], "given-names": ["GK"], "article-title": ["Limma: linear models for microarray data"], "source": ["Bioinformatics and Computational Biology Solutions using R and Bioconductor"], "year": ["2005"], "publisher-name": ["New York , Springer"], "fpage": ["397"], "lpage": ["420"]}]
BMC Genomics. 2008 Jul 31; 9:365