Whole genome sequencing of Saccharomyces cerevisiae: from genotype to phenotype for improved metabolic engineering applications
- José Manuel Otero†1, 2, 5,
- Wanwipa Vongsangnak†1, 2, 8,
- Mohammad A Asadollahi1, 2, 6,
- Roberto Olivares-Hernandes1, 2,
- Jérôme Maury2, 7,
- Laurent Farinelli3,
- Loïc Barlocher3,
- Magne Østerås3,
- Michel Schalk4,
- Anthony Clark4 and
- Jens Nielsen1, 2Email author
© Otero et al; licensee BioMed Central Ltd. 2010
Received: 27 May 2010
Accepted: 22 December 2010
Published: 22 December 2010
The need for rapid and efficient microbial cell factory design and construction are possible through the enabling technology, metabolic engineering, which is now being facilitated by systems biology approaches. Metabolic engineering is often complimented by directed evolution, where selective pressure is applied to a partially genetically engineered strain to confer a desirable phenotype. The exact genetic modification or resulting genotype that leads to the improved phenotype is often not identified or understood to enable further metabolic engineering.
In this work we performed whole genome high-throughput sequencing and annotation can be used to identify single nucleotide polymorphisms (SNPs) between Saccharomyces cerevisiae strains S288c and CEN.PK113-7D. The yeast strain S288c was the first eukaryote sequenced, serving as the reference genome for the Saccharomyces Genome Database, while CEN.PK113-7D is a preferred laboratory strain for industrial biotechnology research. A total of 13,787 high-quality SNPs were detected between both strains (reference strain: S288c). Considering only metabolic genes (782 of 5,596 annotated genes), a total of 219 metabolism specific SNPs are distributed across 158 metabolic genes, with 85 of the SNPs being nonsynonymous (e.g., encoding amino acid modifications). Amongst metabolic SNPs detected, there was pathway enrichment in the galactose uptake pathway (GAL1, GAL10) and ergosterol biosynthetic pathway (ERG8, ERG9). Physiological characterization confirmed a strong deficiency in galactose uptake and metabolism in S288c compared to CEN.PK113-7D, and similarly, ergosterol content in CEN.PK113-7D was significantly higher in both glucose and galactose supplemented cultivations compared to S288c. Furthermore, DNA microarray profiling of S288c and CEN.PK113-7D in both glucose and galactose batch cultures did not provide a clear hypothesis for major phenotypes observed, suggesting that genotype to phenotype correlations are manifested post-transcriptionally or post-translationally either through protein concentration and/or function.
With an intensifying need for microbial cell factories that produce a wide array of target compounds, whole genome high-throughput sequencing and annotation for SNP detection can aid in better reducing and defining the metabolic landscape. This work demonstrates direct correlations between genotype and phenotype that provides clear and high-probability of success metabolic engineering targets. The genome sequence, annotation, and a SNP viewer of CEN.PK113-7D are deposited at http://www.sysbio.se/cenpk.
Metabolic engineering is the enabling technology for identification of targeted genetic modifications such as gene deletion, over-expression, or modulation. The genetic engineering implemented in a host microbial cell factory ideally will lead to re-direction of fluxes to enhance production or robustness of a given product or organism, respectively [1–5]. Metabolic engineering through systems biology has been complimented, and its application expanded in both scope and success. Systems biology is a multi-disciplinary approach to quantitative collection, analysis, and integration of whole genome scale data sets enabling construction of biologically relevant and often predictive mathematical models [6–8]. Genome sequencing of industrially relevant organisms, including S. cerevisiae strain S288c, the first eukaryote genome sequence reported, provided a framework for gene annotation through functional genomics. More relevant to metabolic engineering, an annotated genome sequence was a prerequisite for genome-scale metabolic network reconstructions [9, 10]. Such reconstructions offer a biochemical model describing the formation and depletion of each metabolite that by providing mass-balance boundary conditions makes possible constraint based simulations of how the metabolic network operates at different conditions. In simpler terms, using basic stoichiometry these models can be used to predict the relationships between gene functions in the cellular metabolic network. With nearly 14 years elapsing since the S. cerevisiae strain S288c genome sequence was made available, and more than 1,000 laboratories participating in functional genomics efforts, there are still 968 and 811 open reading frames (ORFs) classified as uncharacterized and dubious, respectively, according to the Saccharomyces Genome Database (SGD) [10, 11]. Furthermore, since 2003 there have been published five major S. cerevisiae genome-scale metabolic network reconstructions, with the most recent models encompassing between 13-14% genome coverage [12, 13]. The opportunity to further extend genotype to phenotype annotation is abundant.
Industrial biotechnology is dominated by efforts to confer a desirable phenotype onto strains using different methods of directed evolution and random mutagenesis, requiring screening and selection. This approach, while providing little to no mechanistic understanding of which specific genetic perturbations lead to improved strains so they could be further exploited, has proven to be commercially successful as illustrated by the more than 1,000 fold improvement in penicillin titer by Penicillium chrysogenum . As industrial biotechnology applications expand, and the desire to custom-engineer microbial cell factories with novel architecture for native and heterologous metabolic pathways increases, the necessity on a genome-wide level to understand direct genotype to phenotype relationships has rapidly increased. Within the same time period of approximately the last 10 years, the technologies and costs associated with whole genome sequencing have advanced and decreased, respectively. There are several excellent reviews of genome sequencing technologies, and their applications to functional genomics, strain engineering, and other investigatory biology efforts [5, 15–18]. Prior work, specifically focused on characterizing genome-wide analysis of nucleotide polymorphisms in S. cerevisiae utilized 25mers oligonucleotide microarrays (Affymetrix yeast tiling arrays) providing random and redundant coverage of the S. cerevisiae genome . This analysis included single nucleotide polymorphism (SNP) identification between S288c and the commonly used laboratory strain S. cerevisiae CEN.PK, where a total of 13,914 SNPs were identified. However, this approach is unable to identify the exact nucleotide substitution, and consequently whether the transcribed SNP results in an amino acid substitution, presumably required to confer a change in enzyme and/or protein function. More recently a collaborative project, the Saccharomyces Genome Resequencing Project (SGRP) between the Sanger Institute and Institute of Genetics, University of Nottingham, completed the ABI sequencing of haploids of 37 S. cerevisiae strains to a coverage of 1-3X. Furthermore, Illumina-Solexa genome sequencing of four of the 37 S. cerevisiae strains, one of which included S288c, was completed . This sequencing effort was focused on exploration of genomic variation in the context of evolution, thereby using multiple strains from different Saccharomyces species. It is a demonstration of a recent genome sequencing technology, referred to as Illumina-Solexa sequencing, compared to larger read methods such as Sanger or 454 sequencing. Illumina-Solexa sequencing is an ultra-high-throughput technology that performs sequencing-by-synthesis of random arrays of clonal DNA colonies attached to the surface of a flow-cell. The approach used in this study generates short, 35 base pair (bp) reads (currently, the technology limitations are 150 bp or 2 × 150 bp paired-reads), that must then be aligned to and assembled using a reference genome [21–23].
In this work we propose that high-throughput genome sequencing and annotation, integrated with a genome browser and SNP viewer of S. cerevisiae may serve as a commonplace tool, complementary to transcriptomics and physiological characterization, to extract direct genotype to phenotype information. Firstly, whole genome Illumina-Solexa sequencing of each strain was completed and then annotation was performed. To access annotated genome sequences and detected SNPs of CEN.PK113-7D, a genome browser and SNP viewer were developed. SNPs strictly related to metabolic genes were identified, characterized, and amino acid level analysis performed. In order to directly link genotype to phenotype cellular behavior was characterized in well-controlled batch fermentations on glucose and galactose, complimented with transcriptome analysis. More specifically, we demonstrate that S288c, the strain utilized for the publically available S. cerevisiae genome sequence, exhibits atypical S. cerevisiae behavior related to central carbon metabolism as compared to CEN.PK113-7D, a common laboratory strain for industrial biotechnology applications .
There were clear correlations between physiology and metabolic pathway enrichment of nonsynonymous SNPs observed, suggesting that sequencing and the annotated genome may assist in reducing the genetic target space for metabolic engineering applications. The analysis presented here serves as a foundation for comparative metabolic engineering SNP analysis, wherein the future reference strains may be compared to their metabolically engineered derivatives that use directed evolution in order to answer the age-old question: what changed in our strain that makes it a preferred microbial cell factory?
Genome Sequencing and Annotation
Illumina/Solexa genome sequencing and annotation results
No. of Reads
No. of Aligned Reads
Total Bases A (bp)
Calculated Average Coverage
Genome Percent Reference Coverage (%)
No. of chromosomes
Genome Browser and SNP Viewer
SNPs identification in metabolism
Genome sequencing and metabolic SNP detection
MAQ Software Determination
No. of Reads
Total No. of SNPs
Total No. of Non-Ambiguous SNPs
Total No. of Filtered SNPs A,B
Metabolism Focused Detection
Total No. of Metabolic Genes Considered C
Total No. of Metabolic Bases (bp)
No. of Aligned Reads
Total Gap Size (bp)
Metabolic Genome Percent Reference Coverage (%)
Total No. of Metabolic SNPs detected
Total No. Nonsynonymous Metabolic SNPs Detected
Percent of SNPs Detected Nonsynonymous (%)
Total No. of Metabolic Genes Containing SNP
Total No. of Metabolic Genes with Nonsynonymous SNP
Physiological characterization of S. cerevisiae strains S288c and CEN.PK113-7D
(n = 2)
(n = 2)
(n = 2)
(n = 2)
(n = 2)
μ-max (h -1 )
Carbon Recovery (%)
Specific Productivity or Consumption Rate A
-r Gluc or -r Gal
A similar characterization was performed using batch galactose supplemented fermentations. CEN.PK113-7D demonstrated a slight lag-phase compared to glucose fermentation; however, sustained a galactose specific growth rate of 0.27 h-1 and galactose uptake rate of 24.3 C-mmol g-DCW-1 h-1, representing a 34% and 77% reduction, respectively, compared to glucose (See Table 3). All extracellular metabolic specific productivity rates were significantly decreased (ethanol, acetate, and glycerol were 93%, 6.8%, and 88% reduced compared to glucose, respectively), with the exception of OUR, which was 47% higher on galactose compared to glucose, leading to an effectively lower RQ of 1.5 compared to 11.9 during glucose cultivation. Furthermore, given the significantly lower RQ during the exponential phase of galactose fermentation, relatively little ethanol was produced (2.7 g L-1), resulting in a short ERF phase (< 5h) (See Figure 5). Similarly, S288c was cultivated on galactose; however, a significant deficiency in the strain's ability to metabolize this carbon source was observed. A total of 25 h post-inoculation elapsed with no increase in biomass as compared to CEN.PK113-7D where after 6 h post-inoculation two cell doublings were observed. At 25 h post-inoculation a glucose bolus of 10 g L-1 was added to promote growth, and rapidly, glucose fermentation, a diauxic shift, and ethanol respiro-fermentation were observed (See Figure 5). Both co-consumption of galactose and ethanol, and a galactose only respiro-fermentative (GaRF) growth phase were observed. During co-consumption the specific growth rate was 0.14 h-1, while on galactose only the specific growth rate was 0.02 h-1. Similarly, the extracellular specific metabolite productivity rates were nearly zero when only galactose consumption was considered (See Table 3). Ethanol was consumed by 82 h post-inoculation, and in the period from 82 h to 128 h, only galactose consumption was observed, and biomass increased from 7.9 g-DCWL L-1 to 20.9 g-DCW L-1, representing a doubling time of 35 h compared to 2.6 h for CEN.PK113-7D.
Summary of Differential Gene Expression
Summary of Differential Expression (p adj < 0.01)
S288c v. CEN.PK113-7D
Glucose ( n = 2)
Galactose ( n = 2)
Total No. of Differentially Expressed Genes
No. of Genes LFC > 0
LFC ± SD
2.13 ± 1.41
1.81 ± 1.17
No. of Genes with SNPs Detected
No. of Genes with Nonsynonymous SNPs
No. of Genes LFC < 0
LFC ± SD
-2.12 ± 1.32
-1.53 ± 1.05
No. of Genes with SNPs Detected
No. of Genes with Nonsynonymous SNPs
For the condition S288c v. CEN.PK113-7D cultivated on glucose, the top 272 differentially expressed genes, ranked according to p adj value are characterized into GO process terms largely dominated by responses to stimuli and pheromone, with the dominant metabolic process categories being trehalose metabolism, steroid metabolism, and amino acid transport. Specific genes consistent with this categorization high in p adj value rank and lfc > 0 are GSY1 (glycogen synthase, lfc 2.0, p adj value rank 23) and for lfc < 0 is HMG1 (HMG-CoA reductase, lfc -1.7, p adj value rank 14). For the condition S288c v. CEN.PK113-7D cultivated on galactose, the top 501 differentially expressed genes, ranked according to p adj value are characterized into GO process terms response to simuli and stress, carbohydrate metabolism, and transport. Specific metabolic genes noteworthy in this category, high in p adj value rank amongst genes with lfc > 0 include MDH2 (malate dehydrogenase, lfc 2.8, p adj value rank 8), FBP1 (fructose-1,6-bisphosphatase, lfc 4.2, p adj value rank 15), GAD1 (glutamate decarboxylase, lfc 3.0, p adj value rank 30), GDH3(NADP+ dependent glutamate dehydrogenase, lfc 3.2, p adj value rank 32), GSY1 (lfc 1.4, p adj value rank 41), and ICL1 (isocitrate lyase, lfc 2.7, p adj value rank 54). Similarly, specific metabolic genes high in p adj value rank amongst genes with lfc < 0 include ARE2 (acyl-coA:sterol acetyltransferase, lfc -2.3, p adj value rank 10), and CYB5 (cytochrome b5, lfc -1.6, p adj value rank 47).
In an effort to further investigate if larger regulatory mechanisms could be identified the list of genes exhibiting significant differential expression were submitted to the Yeast Search for Transcriptional Regulators And Consensus Tracking (YEASTRACT) curated repository of associations between transcription factors and target genes in S. cerevisiae [25, 26]. The transcription factor, Tec1p, was identified as directly regulating 21.1% of the total submitted gene list (See Table 4, 272 genes, S288c glucose v. CEN.PK113-7D glucose), and was 1.7-fold higher expressed in CEN.PK113-7D compared to S288c (p adj value = 7.2 × 10-3). Tec1p was the only identified transcription factor to be significantly differentially expressed, and strongly regulates FLO11, a flocculin gene required for invasive growth, and pseudohyphal formation . The transcription factors regulating the highest percentage of the differentially expressed genes, yet not being differentially expressed themselves, were Sok2p and Ste12p, with 32.5% and 21.5%, respectively, of submitted genes being directly regulated. Sok2p and Ste12p are transcription factors negatively regulating pseudohyphal differentiation . A similar analysis was performed for galactose; however, similar results were obtained, with Sok2p and Ste12p directly regulating 23.1% and 17.4%, respectively, of the 501 differentially expressed genes (See Table 4). The transcription factors differentially expressed themselves were Msa1p (p adj value = 1.37 × 10-1) and Msa2p (p adj value = 1.65 × 10-1), putative G1-specific cell cycle transcription activators, and Usv1p, a putative zinc finger transcription factor regulating growth on non-fermentable carbon sources. USV1 expression was 2.2-fold higher in CEN.PK113-7D compared to S288c (p adj value = 3.6 × 10-3). Although relatively little is known about Usv1p, it has been shown to be induced post-diauxic shift, consistent with the deficiency in post-diauxic shift metabolism observed in S288c . With the exception of Usv1p, all transcription factors identified are more closely related to the significant difference in growth rates between strains rather than their respiro-fermentative metabolism. Metabolic SNPs identified and subsequent analysis did not identify clear correlations or pathway enrichment that could explain the lack of respiro-fermentative metabolism in S288c. Metabolic genes containing nonsynonymous SNPs in CEN.PK113-7D, significantly differentially expressed on galactose, and related to oxidative metabolism included ACS1 (Acetyl-CoA synthetase), GAD1 (Glutamate decarboxylase), YAT2 (Carnitine acetyltransferase), and CCP1 (Mitochondrial cytochrome-c peroxidase) which were 7.8-fold, 3.2-fold, 4.5-fold, and 2.1-fold higher in CEN.PK113-7D, respectively (See Additional file 1 Table S2).
The physiological characterization clearly suggests that S288c has a deficiency in metabolism of respiro-fermentative carbon sources, such as ethanol and galactose, when compared to CEN.PK113-7D. Inspection of the significantly differentially expressed genes between strains cultivated on glucose or galactose did not reveal an obvious gene cluster that would explain this significant physiological difference. This is supported both by the GO characterization and pathway expression mapping.
There were two central carbon metabolic pathways enriched with nonsynonymous SNPs that also correlated with significant differences in phenotype. S. cerevisiae CEN.PK113-7D exhibited significantly higher ergosterol content during growth on glucose, and to a lesser extent, galactose. This is consistent with previous work where CEN.PK2-1C had very high ergosterol/erg-ester (20.0 mg g-DCW-1) and triacylglycerols content (15.2 mg g-DCW-1) compared to 9 other S. cerevisiae strains, including FY169 (ergosterol/erg-ester content: 8.5 mg g-DCW-1; triacylglycerols content: 2.4 mg g-DCW-1) which is isogenic to S288c [29, 30]. The ergosterol biosynthetic pathway had significant nonsynonymous SNPs identified in ERG8 and ERG9, and silent SNPs identified in ERG20 and HMG1. Both ERG8 and ERG9 were not significantly differentially expressed, either in glucose or galactose, suggestive again that phenotypic observations, consistent with genome sequence variations, are not necessarily directly manifested at the transcriptome level. Both ERG8 (encodes phosophomevalonate kinase) and ERG9 (encodes squalene synthesase) are essential cytosolic enzymes in the biosynthetic pathway of isoprenoids and sterols (Δerg8 and Δerg9, were found to both be auxtrophic for ergosterol in the systematic deletion library), including ergosterol, from mevalonate [11, 31, 32]. The ergosterol biosynthetic pathway is highly regulated through feedback inhibition mechanisms and by several rate-controlling steps, including that catalyzed by HMG-CoA reductase, encoded by HMG1 [33, 34]. Under both glucose and galactose, HMG1 expression was significantly down-regulated in S288c compared to CEN.PK113-7D by 3.2-fold (p adj value = 3.3 × 10-4) and 1.8-fold (p adj value = 8.6 × 10-3), respectively, correlating with the significantly less ergosterol content in S288c cultivated on glucose and to a lesser extent, on galactose. Furthermore, ERG9 has been previously identified as also having a regulatory role , consistent with the hypothesis that a nonsynonymous SNP resulting in altered protein function could affect ergosterol synthesis. ERG8 on the other hand has not been explicitly shown to have a regulatory function, yet, when the specific activity of 0.06 μmol min-1 mg-1 is compared to other ergosterol synthetic enzymes such as ERG13 (2.1 in S. cerevisiae), ERG12 (0.77 in S. cerevisiae), ERG20 (5.22 in S. cerevisiae), and especially the known regulator HMG1/HMG2 (0.0035 in S. cerevisiae), it is suggestive that ERG8 is likely a rate limiting step [36–44]. There were a large number of nonsynonymous SNPs that encoded significant changes in amino acid classes, further suggestive that ERG8 is a strong metabolic engineering target for understanding the significantly higher ergosterol content in CEN.PK113-7D. Lastly, the observation that neither ERG8 nor ERG9 were differentially expressed under glucose or galactose, suggests their potential affect on phenotype is likely post-translational.
Similar to ergosterol biosynthesis, the galactose uptake pathway phenotype in S288c was vastly lower rate compared to CEN.PK113-7D, correlating with the nonsynonymous SNP enrichment in GAL1 and GAL10, and silent SNPs in GAL7. Neither GAL1 (encodes galactokinase) nor GAL10 (encodes UDP-glucose-4-epimerase) were significantly differentially expressed during growth on galactose; however, on glucose GAL1 was significantly up-regulated (p adj value = 9.7 × 10-4), 2.9-fold in CEN.PK113-7D. Both Δgal1 and Δgal10 mutants are unable to grow on galactose as sole carbon sources [45–47]. The significant number of nonsynonymous SNPs in both essential galactose genes suggests obvious targets for explanation of why S288c is incapable of galactose respiro-fermentative metabolism. Furthermore, it should be noted that S288c has been described as Δgal2 (See Additional file 1 Table S3), which may be ascribed to the presence of 4 SNPs between CEN.PK113-7D and S288c in the coding region, and 14 SNPs in the upstream region that result in a poor functional Gal2p. However, clearly, S288c is able to co-metabolize galactose with ethanol and this may be ascribed by galactose transport by e.g. the hexose transporters.
A further metabolic engineering benefit of whole genome sequencing was the detection of a nonsynonymous SNP resulting in a stop codon of PAD1 (encodes phenylacrylic acid decarboxylase). Pad1p is essential for decarboxylation of aromatic carboxylic acids conferring resistance to cinnamic acid, and a nonsynonymous SNP was detected at nucleotide position 294 (T to G), resulting in a stop codon (TAT → TAG) . Although Pad1 relevant phenotypes were not explored, the transcriptome response on glucose revealed significant differential expression of PAD1 (p adj value = 1.5 × 10-3), with 3.1-fold higher expression in S288c compared to CEN.PK113-7D. This is consistent with the stop codon detected in CEN.PK113-7D at position 294, noting that the total ORF genomic DNA sequence is 729 nucleotides, and therefore unlikely to be transcribed and detected.
Whole high-throughput genome sequencing of S. cerevisiae S288c and CEN.PK113-7D resulted in identification of 13,787 filtered SNPs in CEN.PK113-7D, with a total of 939 SNPs detected across 158 unique metabolic genes, 85 of which contained a total of 219 nonsynonymous SNPs. There were two central carbon metabolic pathways enriched with nonsynonymous SNPs that also correlated with significant differences in phenotype. S. cerevisiae CEN.PK113-7D exhibited significantly higher ergosterol content during growth on glucose, and to a lesser extent, galactose. The ergosterol biosynthetic pathway had significant nonsynonymous SNPs identified in ERG8 and ERG9, and silent SNPs identified in ERG20 and HMG1. The flux through the galactose uptake pathway was much lower in S288c compared with CEN.PK113-7D, correlating with the nonsynonymous SNP enrichment in GAL1 and GAL10, and silent SNPs in GAL7. More globally, the physiological characterization clearly suggests that S288c has a deficiency in metabolism of respiratory carbon sources, such as ethanol and galactose, when compared to CEN.PK113-7D. Inspection of the significantly differentially expressed genes between strains cultivated on glucose or galactose did not reveal an obvious gene cluster that would explain this significant physiological difference. In summary and perhaps not surprisingly, transcriptome analysis did not provide a clear hypothesis for major phenotypes observed, suggesting that genotype to phenotype correlations are manifested post-transcriptionally or post-translationally either through protein concentration and/or function. Clearly, future work must validate these correlations through genetic engineering of identified SNPs in either S288c or CEN.PK113-7D to see if desired phenotypes, such as increased galactose uptake or ergosterol synthesis in S288c, are observed. Future work must also expand on the metabolic SNP analysis presented to include all 13,787 SNPs, realizing phenotypic observations may not necessarily be linked directly to metabolic SNPs, but rather SNPs affecting larger regulatory mechanisms and networks, such as those governed by transcription factors. Certainly, as S. cerevisiae continues to be exploited, particularly for metabolic engineering applications, the integration of physiological characterization, transcriptome analysis, and metabolic SNP detection with high-throughput whole genome sequencing provides direct correlations between observed phenotypes and genotypes and offers high probability of success metabolic targets.
S. cerevisiae strain S288c (American Type Culture Collection, ATCC®) and strain CEN.PK113-7D (Scientific Research and Development, SRD GmbH were used in this study. Genotype of strain S288c was described by Mortimer et al, , Johnston et al , Goffeau et al , Cherry et al  and genotype of strain CEN.PK113-7D was described by Cherry et al , van Dijken et al . More information of S. cerevisiae strains was described in Additional file 1 Table S3.
Genome Sequencing, Annotation and Analysis
A standard 500 mL shake flask, supplemented with 10 g L-1 glucose and inoculated with a single colony of S. cerevisiae S288c or CEN.PK113-7D, was permitted to grow for 24-48 h at 30°C until visual inspection confirmed a high optical density. A total of 5 mL culture was aliquoted into 15 mL sterile tubes (one per extraction), centrifuged (4000 RCF) for 5 min, washed with 2 mL deionized water, and pelleted. Cell pellets were resuspended in 0.5 mL lysis buffer. Lysis buffer consisted of 0.1 M Tris pH 8.0, 50 mM EDTA, and 1% SDS final concentration. The lysis buffer suspension was transferred to a 1.5 mL FastPrep screw cap tube, to which 200 μL acid-washed glass beads (250-500 μm) and 25 μL 5 M NaCl was added. A FastPrep™ FP120 (QBiogene, Irwine, CA) was used for cell lysis, with two cycles of 20 s disruption and 1 min on ice. The resulting cell suspension was centrifuged (13,000 RCF) for 10 min., and the resulting clear liquid, approximately 350 μL, avoiding white cell debris and beads, was aspirated with a pipette and transferred to 1.5 mL microcentrifuge tubes. 400 μL chloroform (TE-saturated) was added to each tube, mixed, and a chloroform extraction performed. 1 mL 99% ethanol was added to the resulting suspension, mixed, centrifuged (13,000 RCF) for 6 min., ethanol decanted, and then resuspended in 70% ethanol. The resulting suspension was centrifuged (13,000 RCF) for 6 min., ethanol decanted, and pellet permitted to dry for 25-60 min. The pellet was then resuspended in 50 μL 2 mM Tris, incubated for 10 min. at 37°C, and stored at -20°C.
Illumina/Solexa Genome Sequencing and SNP Analysis
Isolated DNA from S. cerevisiae S288c and CEN.PK113-7D were processed to prepare short insert shotgun libraries at Fasteris SA (Geneva, Switzerland), utilizing the Solexa technology according to the manufacturer's recommendations (Illumina). The whole-genome sequencing was performed in 2007 on a Genome Analyzer "classic" instrument with sequencing kits version 1 and base calling on the Solexa Pipeline (version 0.2.2.5). SNP detection was performed using two independent approaches: mapping on the S288c reference sequence and de novo assembly. The MAQ software package http://maq.sourceforge.net was used to map the short base paired-reads (35 bp) to the reference genome S. cerevisiae S288c version 12.0, available at the Saccharomyces Genome Database (SGD). For consensus calling, the maq assemble -m command was used to call the consensus sequences from read mapping. The value m set to 1 specifies the maximum numbers of mismatches allowed for a filtered read. To detect high-quality SNP identification, the command maq.pl SNPfilter was set up. The threshold values applied to SNP detection for the CEN.PK113-7D sequence relative to the reference sequence were read depth (-d) > 5 for few reads, read depth (-D) < 255 for too many reads, uniqueness of the region (-w) < 1.5, quality of the read (-Q) > 50, and best read quality (-n) > 40. These threshold parameter values were tested such that the amount of coverage and proportion of genome with aligned sequences was maximized, and a graphical representation of SNPs was produced to confirm results. The de novo assembly has been performed using the EDENA software package http://www.genomic.ch/edena.php as described previously . The EDENA assembly results for both whole genome sequencing and for SNP detection in metabolic genes are presented in Additional file 1, Table S4. The de novo assembly used for SNP detection was essentially used for an independent verification of the SNPs detected by mapping on the reference sequence, as not all the genome was represented. Indeed a higher coverage, as well as longer insert libraries are required to achieve more complete de novo assembly. For purposes of the subsequent SNP analysis only the MAQ software results were used. The FASTA files of each genome sequence are available upon request.
All metabolic genes containing SNPs, both silent and nonsynonymous, were manipulated within the software BioEdit v7.08 http://www.mbio.ncsu.edu/BioEdit/bioedit.html. Specifically, the ORF genomics nucleotide sequence available on Saccharomyces Genome Database (SGD) http://www.yeastgenome.org were imported into BioEdit, and the sequences modified with the identified SNP, creating a new CEN.PK113-7D sequence for that ORF relative to the original S288c strain. Both the S288c and CEN.PK113-7D nucleotide sequences were then translated in fix full frames, and amino acid polymorphisms were identified, leading to the categorization of each SNP as either being silent or nonsynonymous. Subsequent physiological characterization of the gene and all relevant amino acid information from UnitProt were managed in a spreadsheet using Microsoft Excel. Multi-sequence Pfam alignments were performed using a custom BioPerl script and the UNIX operating environment. Calculations and characterization described in Additional file 1, Figure S3, related to amino acids, were then performed using Microsoft Excel.
Gene Prediction and Functional Annotation
Genes were predicted in the S. cerevisiae strain CEN.PK113-7D genome based on the homologies to known or putative genes in the public database of the S. cerevisiae S288c v12.0, sequence available at the Saccharomyces Genome Database (SGD). The statistical features of the genes were predicted by a combination of gene-finding software. To perform annotation, BLASTN and GlimmerHMM were applied. GlimmerHMM was used for prediction the gene structures of strain CEN.PK113-7D by training the gene models from the strain S288c. All of the predicted protein-coding genes were annotated by searching against the SGD database using BLASTP, followed by manual curation.
Development of Genome Browser and SNPs Viewer
The genome browser and SNP viewer were implemented by PHP source code and MySQL database management in our hosting server http://members.dot5hosting.com/controlpanel/. Currently, our MySQL database structure contains different features that correspond to chromosome location, gene position, gene function, GO process, metabolism, and SNPs.
A chemically defined minimal medium of composition 5.0 g L-1 (NH4)2SO4, 3.0 g L-1 KH2PO4, 0.5 g L-1 MgSO4.7H2O, 1.0 mL L-1 trace metal solution, 300 mg L-1 uracil, 0.05 g L-1 antifoam 204 (Sigma-Aldrich A-8311), and 1.0 mL L-1 vitamin solution was used for all shake flask and 2L well-controlled fermentations . The trace elment solution included 15 g L-1 EDTA, 0.45 g L-1 CaCl2.2H2O, 0.45 g L-1 ZnSO4 .7H2O, 0.3 g L-1 FeSO4.7H2O, 100 mg L-1 H3BO4, 1 g L-1 MnCl2.2H2O, 0.3 g L-1 CoCl2.6H2O, 0.3 g L-1 CuSO4.5H2O, 0.4 g L-1 NaMoO4.2H2O. The pH of the trace metal solution was adjusted to 4.0 with 2 M NaOH and heat sterilized. The vitamin solution included 50 mg L-1 d-biotin, 200 mg L-1 para- amino benzoic acid, 1 g L-1 nicotinic acid, 1 g L-1 Ca.pantothenate, 1 g L-1 pyridoxine HCl, 1 g L-1 thiamine HCl, and 25 mg L-1 m.inositol. The pH of the vitamin solution was adjusted to 6.5 with 2 M NaOH, sterile-filtered and the solution was stored at 4°C. The final formulated medium, excluding glucose and vitamin solution supplementation, is adjusted to pH 5.0 with 2 M NaOH and heat sterilized. For carbon-limited cultivations the sterilized medium is supplemented with 40 g L-1 glucose or 40 g L-1 galactose, heat sterilized separately, and 1.0 mL L-1 vitamin solution is added by sterile filtration (0.20 μm pore size Ministart®-Plus Sartorius AG, Goettingen, Germany).
Shake Flask Cultivations and Stirred Tank Fermentations
Shake flask cultivations were completed in 500 mL Erlenmeyer flasks with two diametrically opposed baffles and two side-necks with septums for sampling by syringe. Flasks were heat sterilized with 100 mL of medium, inoculated with a single colony, and incubated at 30°C with orbital shaking at 150 RPM. Stirred tank fermentations were completed in well-controlled, aerobic, 2.2L Braun Biotech Biostat B fermentation systems with a working volume of 2L (Sartorius AG, Goettingen, Germany). The temperature was controlled at 30°C. The fermenters were outfitted with two disk-turbine impellers rotating at 600 RPM. Dissolved oxygen was monitored with an autoclavable polarographic oxygen electrode (Mettler-Toledo, Columbus, OH). During aerobic cultivation the air sparging flow rate was 2 vvm. The pH was kept constant at 5.0 by automatic addition of 2 M KOH. Off-gas passed through a condenser to minimize the evaporation from the fermenter. The fermenters were inoculated from shake flask precultures to an initial OD600 0.01. Two biological replicates were performed for each fermentation condition.
The effluent fermentation gas was measured every 30 seconds for determination of O2(g) and CO2(g) concentrations by the off-gas analyzer Brüel and Kjær 1308 (Brüel & Kjær, Nærum, Denmark).
The optical density (OD) was determined at 600 nm using a Shimadzu UV mini 1240 spectrophotometer (Shidmazu Europe GmbH, Duisberg, Germany). Duplicate samples were diluted with deionized water to obtain OD600 measurements in the linear range of 0-0.4 OD600 Samples were always maintained at 4°C post-sampling until OD600 and dry cell weight (DCW) measurements were performed. DCW measurements were determined through the exponential phase, until stationary phase was confirmed according to OD600 and off-gas analysis. Nitrocellulose filters (0.45 μm Sartorius AG, Goettingen, Germany) were used. The filters were pre-dried in a microwave oven at 150W for 10 min., and cooled in a dessicator for 10 min. 5.0 mL of fermentation broth were filtered, followed by 10 mL DI water. Filters were then dried in a microwave oven for 20 min. at 150W, cooled for 15 min. in a desiccator, and the mass was determined.
Metabolite Concentration Determination
All fermentation samples were immediately filtered using a 0.45 μm syringe-filter (Sartorius AG, Goettingen, Germany) and stored at -20°C until further analysis. Glucose, ethanol, glycerol, acetate, succinate, pyruvate, fumarate, citrate, oxalate, and malate were determined by HPLC analysis using an Aminex HPX-87 H ion-exclusion column (Bio-Rad Laboratories, Hercules, CA). The column was maintained at 65°C and elution performed using 5 mM H2SO4 as the mobile phase at a flow rate of 0.6 mL min-1. Glucose, ethanol, glycerol, acetate, succinate, citrate, fumarate, malate, oxalate were detected on a Waters 410 differential refractometer detector (Shodex, Kawasaki, Japan), and acetate and pyruvate were detected on a Waters 468 absorbance detector set at 210 nm. Ergosterol measurements were made according to previous published methods .
RNA Sampling and Isolation
Samples for RNA isolation from two biological replicates from the late-exponential phase of glucose-limited and galactose-limited batch cultivations were taken by rapidly sampling 25 mL of culture into a 50 mL sterile Falcon tube with 40 mL of crushed ice in order to decrease the sample temperature to below 2°C in less than 10 s. Cells were immediately centrifuged (4000 RCF at 0°C for 2.5 min.), the supernatant discarded, and the pellet frozen in liquid nitrogen and it was stored at -80°C until total RNA extraction. Total RNA was extracted using the FastRNA Pro RED kit (QBiogene, Carlsbad, USA) according to manufacturer's instructions after partially thawing the samples on ice. RNA sample integrity and quality was determined prior to hybridization with an Agilent 2100 Bioanalyzer and RNA 6000 Nano LabChip kit according to the manufacturer's instruction (Agilent, Santa Clara, CA).
Probe Preparation and Hybridization to DNA Microarrays
mRNA extraction, cDNA synthesis, labeling, and array hybridization to Affymetrix Yeast Genome Y2.0 arrays were performed according to the manufacturer's recommendations (Affymetrix GeneChip® Expression Analysis Technical Manual, 2005-2006 Revision 2.0). Washing and staining of arrays were performed using the GeneChip Fluidics Station 450 and scanning with the Affymetrix GeneArray Scanner (Affymetrix, Santa Clara, CA).
Microarray Gene Transcription Analysis
Affymetrix Microarray Suite v5.0 was used to generate CEL files of the scanned DNA microarrays. These CEL files were then processed using the statistical language and environment R v5.3 (R Development Core Team, 2007, http://www.r-project.org), supplemented with Bioconductor v2.3 (Biconductor Development Core Team, 2008, http://www.bioconductor.org) packages Biobase, affy, gcrma, and limma . The probe intensities were normalized for background using the robust multiarray average (RMA) method only using perfect match (PM) probes after the raw image file of the DNA microarray was visually inspected for acceptable quality. Normalization was performed using the qspline method and gene expression values were calculated from PM probes with the median polish summary. Statistical analysis was applied to determine differentially expressed genes using the limma statistical package. Moderated t- tests between the sets of experiments were used for pair-wise comparisons. Empirical Bayesian statistics were used to moderate the standard errors within each gene and Benjamini-Hochberg's method was used to adjust for multi-testing. A cut-off value of adjusted p < 0.01 was used for statistical significance . Statistically significant differential gene expression lists were then submitted to the GO Term Finder (version 0.83) of the Saccharomyces Genome Database (SGD) for GO process, function, and component statistically significant identification (p < 0.01). Furthermore, the same differential gene expression lists were submitted to the Expression Viewer (Pathway Tools version 12.0 generated by SRI International on SGD) for metabolic pathway mapping and identification .
Genome sequence and annotated gene of S. cerevisiae CEN.PK113-7D were deposited at sysbio database http://www.sysbio.se/cenpk. Normalized gene expression data were deposited at the GEO database http://www.ncbi.nlm.nih.gov/geo/ with accession number GPL2529 (platform), GSM536874-GSM536881 (samples), and GSE21479 (series).
We would like to thank Merck Doctoral Fellowship for financial support to JMO and Chalmers University of Technology Foundation for financial support to WV. We would like to thank Kwanjeera Wanichthanarak for supporting on database development and Gionata Scalcinati for assisting on manuscript revision.
- Bailey JE: Toward a science of metabolic engineering. Science. 1991, 252 (5013): 1668-1675. 10.1126/science.2047876.PubMedView ArticleGoogle Scholar
- Nielsen J: Metabolic engineering. Appl Microbiol Biotechnol. 2001, 55 (3): 263-283. 10.1007/s002530000511.PubMedView ArticleGoogle Scholar
- Stephanopoulos G, Vallino JJ: Network rigidity and metabolic engineering in metabolite overproduction. Science. 1991, 252 (5013): 1675-1681. 10.1126/science.1904627.PubMedView ArticleGoogle Scholar
- Tyo KE, Alper HS, Stephanopoulos GN: Expanding the metabolic engineering toolbox: more options to engineer cells. Trends Biotechnol. 2007, 25 (3): 132-137. 10.1016/j.tibtech.2007.01.003.PubMedView ArticleGoogle Scholar
- Warner JR, Patnaik R, Gill RT: Genomics enabled approaches in strain engineering. Curr Opin Microbiol. 2009, 12 (3): 223-230. 10.1016/j.mib.2009.04.005.PubMedView ArticleGoogle Scholar
- Nielsen J, Jewett MC: Impact of systems biology on metabolic engineering of Saccharomyces cerevisiae. FEMS Yeast Res. 2008, 8 (1): 122-131. 10.1111/j.1567-1364.2007.00302.x.PubMedView ArticleGoogle Scholar
- Oliver SG: From genomes to systems: the path with yeast. Philos Trans R Soc Lond B Biol Sci. 2006, 361 (1467): 477-482. 10.1098/rstb.2005.1805.PubMed CentralPubMedView ArticleGoogle Scholar
- Westerhoff HV, Palssøn BO: The evolution of molecular biology into systems biology. Nat Biotechnol. 2004, 22 (10): 1249-1252. 10.1038/nbt1020.PubMedView ArticleGoogle Scholar
- Förster J, Famili I, Fu P, Palssøn , Nielsen J: Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. Genome Res. 2003, 13 (2): 244-53. 10.1101/gr.234503.PubMed CentralPubMedView ArticleGoogle Scholar
- Goffeau A, Barrell BG, Bussey H, Davis RW, Dujon B, Feldmann H, Galibert F, Hoheisel JD, Jacq C, Johnston M, Louis EJ, Mewes HW, Murakami Y, Philippsen P, Tettelin H, Oliver SG: Life with 6000 genes. Science. 1996, 274 (5287): 546-563. 10.1126/science.274.5287.546. 7PubMedView ArticleGoogle Scholar
- Cherry JM, Ball C, Weng S, Juvik G, Schmidt R, Adler C, Dunn B, Dwight S, Riles L, Mortimer RK, Botstein D: Genetic and physical maps of Saccharomyces cerevisiae. Nature. 1997, 387: 67-73. 10.1038/43025.PubMed CentralPubMedView ArticleGoogle Scholar
- Herrgård MJ, Swainston N, Dobson P, Dunn WB, Arga KY, Arvas M, Blüthgen N, Borger S, Costenoble R, Heinemann M, Hucka M, Le Novère N, Li P, Liebermeister W, Mo ML, Oliveira AP, Petranovic D, Pettifer S, Simeonidis E, Smallbone K, Spasić I, Weichart D, Brent R, Broomhead DS, Westerhoff HV, Kirdar B, Penttilä M, Klipp E, Palsson BØ, Sauer U, Oliver SG, Mendes P, Nielsen J, Kell DB: A consensus yeast metabolic network reconstruction obtained from a community approach to systems biology. Nat Biotechnol. 2008, 26 (10): 115-1160. 10.1038/nbt1492.Google Scholar
- Nookaew I, Jewett MC, Meechai A, Theammarongtham C, Laoteng K, Cheevadhanarak S, Nielsen J, Bhumriatana S: The genome-scale metabolic model iIN800 of Saccharomyces cerevisiae and its validation: a scaffold to query lipid metabolism. BMC Syst Biol. 2008, 2: 71-10.1186/1752-0509-2-71.PubMed CentralPubMedView ArticleGoogle Scholar
- Nielsen J: Physiological engineering aspects of Penicillium chrysogenum. Polyteknisk Forlag, Lyngby. 1995, 13-Google Scholar
- Kahvejian A, Quanckenbush J, Thompson JF: What would you do if you could sequence everything?. Nat Biotechnol. 2008, 26 (10): 1125-1133. 10.1038/nbt1494.PubMed CentralPubMedView ArticleGoogle Scholar
- Morozova O, Marra MA: Applications of next-generation sequencing technologies in functional genomics. Genomics. 2008, 92: 255-264. 10.1016/j.ygeno.2008.07.001.PubMedView ArticleGoogle Scholar
- Shendure J, Ji H: Next-generation DNA sequencing. Nat Biotechnol. 2008, 26 (10): 1135-1145. 10.1038/nbt1486.PubMedView ArticleGoogle Scholar
- Srivatsan A, Han Y, Peng J, Tehranchi AK, Gibbs R, Wang JD, Chen R: High-precision, whole-genome sequencing of laboratory strains facilitates genetic studies. PLOS Genetics. 2008, 4 (8): e1000139-10.1371/journal.pgen.1000139.PubMed CentralPubMedView ArticleGoogle Scholar
- Schacherer J, Ruderfer DM, Gresham D, Dolinski K, Botstein D, Kruglyak : Genome-wide analysis of nucleotide-level variation in commonly used Saccharomyces cerevisiae strains. PLOS One. 2007, 3: e322-10.1371/journal.pone.0000322.View ArticleGoogle Scholar
- Carter D: Saccharomyces Genome Resequencing Project: User Manual. Last accessed September 4, 2009 (Original publication: 2008), [http://www.sanger.ac.uk/Teams/Team118//sgrp/]
- Hernandez D, François P, Farinelli L, Osterås M, Schrenzel J: De novo bacterial genome sequencing: millions of very short reads assembled on a desktop computer. Genome Res. 2008, 18 (5): 802-809. 10.1101/gr.072033.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Nikolaev SI, Iseli C, Sharp AJ, Robyr D, Rougemont J, Gehrig C, Farinelli L, Antonarakis SE: Detection of genomic variation by selection of a 9 Mb DNA region and high throughput sequencing. PLOS One. 2009, 4 (8): e6659-10.1371/journal.pone.0006659.PubMed CentralPubMedView ArticleGoogle Scholar
- Rougemonet J, Amzallag A, Iseli C, Farinelli L, Xenarios I, Naef F: Probabilistic case calling of Solexa sequencing data. BMC Bioinformatics. 2008, 9: 431-10.1186/1471-2105-9-431.View ArticleGoogle Scholar
- van Dijken JP, Bauer J, Brambilla L, Duboc P, François JM, Gancedo C, Giuseppin ML, Heijnen JJ, Hoare M, Lange HC, Madden EA, Niederberger P, Nielsen J, Parrou JL, Petit T, Porro D, Reuss M, van Riel M, Rizzi M, Steensma HY, Verrips CT, Vindeløv J, Pronk JT: An interlaboratory comparison of physiological and genetic properties of four Saccharomyces cerevisiae strains. Enzyme Microb Technol. 2000, 26 (9-10): 706-714. 10.1016/S0141-0229(00)00162-9.PubMedView ArticleGoogle Scholar
- Monteiro PT, Mendes N, Teixiera MC, d-Orey S, Tenreiro S, Mira N, Pais H, Francisco AP, Carvalho AM, Lourenço A, Sá-Correia I, Oliveria AL, Freitas AT: YEASTRACT-DICOVERER: new tools to improve the analysis of transcriptional regulatory association in Saccharomyces cerevisiae. Nucl Acids Res. 2008, 36: D132-D136. 10.1093/nar/gkm976.PubMed CentralPubMedView ArticleGoogle Scholar
- Teixeria MC, Monteiro P, Jain P, Tenreiro S, Fernandes AR, Mira NP, Alenquer M, Freitas AT, Oliveira AL, Sá-Correia I: The YEASTRACT database: a tool for the analysis of transcription regulatory associations in Saccharomyces cerevisiae. Nucl Acids Res. 2006, 34: D446-D451. 10.1093/nar/gkj013.View ArticleGoogle Scholar
- Douglas LM, Li L, Yang Y, Dranginis AM: Expression and characterization of the flocculin Flo11/Muc1, a Saccharomyces cerevisiae mannoprotein with homotypic properties of adhesion. Eukaryot Cell. 2007, 6 (12): 2214-2221. 10.1128/EC.00284-06.PubMed CentralPubMedView ArticleGoogle Scholar
- McCammon MT, Epstein CB, Przybyla-Zawislak B, McAlister-Henn L, Butow RA: Global transcription analysis of Krebs tricarboxylic acid cycle mutants reveals an alternating pattern of gene expression and effects of hypoxic and oxidative genes. Mol Biol Cell. 2003, 14 (3): 958-972. 10.1091/mbc.E02-07-0422.PubMed CentralPubMedView ArticleGoogle Scholar
- Daum G, Tuller G, Nemec T, Hrastnik C, Balliano G, Cattel L, Milla P, Rocco F, Conzelmann A, Vionnet C, Kelly DE, Kelly S, Schweizer E, Schüller H-J, Hojad U, Greiner E, Finger K: Systematic analysis of yeast strains with possible defects in lipid metabolism. Yeast. 1999, 15 (7): 601-614. 10.1002/(SICI)1097-0061(199905)15:7<601::AID-YEA390>3.0.CO;2-N.PubMedView ArticleGoogle Scholar
- Winston F, Dollard C, Ricupero-Hovasse SL: Construction of a set of convenient Saccharomyces cerevisiae strains that are isogenic to S288c. Yeast. 1995, 11 (1): 53-55. 10.1002/yea.320110107.PubMedView ArticleGoogle Scholar
- Jennings SM, Tsay TH, Fisch TM, Robinson GW: Molecular cloning and characterization of the yeast gene for squalene synthesis. Proc Natl Acad Sci USA. 1991, 88: 6038-6042. 10.1073/pnas.88.14.6038.PubMed CentralPubMedView ArticleGoogle Scholar
- Tsay YH, Robinson GW: Cloning and characterization of ERG8, an essential gene of Saccharomyces cerevisiae that encodes phosphomevalonate kinase. Mol Cell Biol. 1991, 11 (2): 620-631.PubMed CentralPubMedView ArticleGoogle Scholar
- Basson ME, Thorsness M, Finer-Moore J, Stroud RM, Rine J: Structural and functional conservation between yeast and human 3-hydroxy-3-methglutaryl coenzyme A reductases, the rate-limiting enzyme of sterol biosynthesis. Mol Cell Biol. 1988, 8 (9): 3797-3808.PubMed CentralPubMedView ArticleGoogle Scholar
- Maury J, Asadollahi MA, Møller K, Clark A, Nielsen J: Microbial isoprenoid production: an example of green chemistry through metabolic engineering. Adv Biochem Engin/Biotechnol. 2005, 100: 19-51. full_text.View ArticleGoogle Scholar
- Grabowska D, Karst F, Szkopinska A: Effect of squalene synthase gene disruption on synthesis of polyprenols in Saccharomyces cerevisiae. FEBS Lett. 1998, 434 (3): 406-408. 10.1016/S0014-5793(98)01019-9.PubMedView ArticleGoogle Scholar
- Basson ME, Thorsness M, Rine J: Saccharomyces cerevisiae contains two functional genes encoding 3-hydroxy-3-methylglutaryl-coenzyme A reductase. Proc Natl Acad Sci USA. 1986, 83 (15): 5563-5567. 10.1073/pnas.83.15.5563.PubMed CentralPubMedView ArticleGoogle Scholar
- Bloch K, Chaykin S, Phillips AH, de Waard A: Mevalonic acid pyrophosphate and isopentenylpyrophosphate. J Biol Chem. 1959, 234: 2595-2604.PubMedGoogle Scholar
- Durr IF, Rudney H: The reduction of beta-hydroxy-beta-methyl-glutaryl coenzyme A to mevalonic acid. J Biol Chem. 1960, 235: 2572-2578.PubMedGoogle Scholar
- Eberhardt NL, Rilling HC: Prenyltransferase from Saccharomyces cerevisiae. Purification to homogeneity and molecular properties. J Biol Chem. 1975, 250 (3): 863-866.PubMedGoogle Scholar
- Gray JC, Kekwick RG: The inhibition of plant mevalonate kinase preparations by prenyl pyrophosphates. Biochim Biophys Acta. 1972, 279 (2): 290-296.PubMedView ArticleGoogle Scholar
- Middleton B, Tubbs PK: 3-Hydroxy-3-methylglutaryl-CoA synthase from baker's yeast. Methods Enzymol. 1975, 35: 173-177. 10.1016/0076-6879(75)35153-7.PubMedView ArticleGoogle Scholar
- Porter JW: Mevalonate kinase. Methods Enzymol. 1985, 110: 71-78. full_text.PubMedView ArticleGoogle Scholar
- Rilling HC: Eukaryotic prenyltransferases. Methods Enzymol. 1985, 110: 145-152. full_text.PubMedView ArticleGoogle Scholar
- Tchen TT: Mevalonic kinase: purification and properties. J Biol Chem. 1958, 233 (5): 1100-1103.PubMedGoogle Scholar
- Bhat PJ, Hopper JE: Overproduction of the GAL1 or GAL3 protein causes galactose-independent activation of the GAL4 protein: evidence for a new model of induction for the yeast GAL/MEL region. Mol Cell Biol. 1992, 12 (6): 2701-2707.PubMed CentralPubMedView ArticleGoogle Scholar
- Bhat PJ, Oh D, Hopper JE: Analysis of the GAL3 signal transduction pathway activating GAL4 protein-dependent transcription in Saccharomyces cerevisiae. Genetics. 1990, 125 (2): 281-291.PubMed CentralPubMedGoogle Scholar
- Douglas HC, Hawthorne DC: Enzymatic expression and genetic linkage of genes controlling galactose utilization in Saccharomyces. Genetics. 1964, 49: 837-844.PubMed CentralPubMedGoogle Scholar
- Clausen M, Lamb CJ, Megnet R, Doerner PW: PAD1 encodes phenylacrylic acid decarboxylase which confers resistance to cinnamic acid in Saccharomyces cerevisiae. Gene. 1994, 142 (1): 107-112. 10.1016/0378-1119(94)90363-8.PubMedView ArticleGoogle Scholar
- Mortimer RK, Johnston JR: Genealogy of principal strains of the yeast genetic stock center. Genetics. 1986, 113 (1): 35-43.PubMed CentralPubMedGoogle Scholar
- Johnston M, Andrew S, Brinkman R, Cooper J, Ding H, Dover J, Favello A, Fulton L, Gattung S, Geisel C, Kirsten J, Kucaba T, Hillier L, Jier M, Johnston L, Langston Y, Latreille P, Louis EJ, Macri C, Mardis E, Menezes S, Mouser L, Nhan M, Rifkin L, Riles L, Peter HS, Trevaskis E, Vaughan K, Vignati D, Wilcox L, Wohldman P, Waterston R, Wilson R, Vaudin M: Complete nucleotide sequence of Saccharomyces cerevisiae chromosome VIII. Science. 1994, 265 (5181): 2077-82. 10.1126/science.8091229.PubMedView ArticleGoogle Scholar
- Verudyn C, Postma E, Scheffers WA, Van Dijken JP: Effect of benzoic acid on metabolic fluxes in yeasts: a continuous culture study on the regulation of respiration and alcoholic fermentation. Yeast. 1992, 8: 501-517. 10.1002/yea.320080703.View ArticleGoogle Scholar
- Asadollahi MA, Maury J, Møller K, Nielsen KF, Schalk M, Clark A, Nielsen J: Production of plant sesquiterpenes in Saccharomyces cerevisiae: effect of ERG9 repression on sesquiterpene biosynthesis. Biotechnol Bioeng. 2008, 99 (3): 666-677. 10.1002/bit.21581.PubMedView ArticleGoogle Scholar
- Smyth GK: Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, Springer, New York, 397-420. full_text.View ArticleGoogle Scholar
- Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Statistical Applications in Genetics and Molecular Biology. 2000, 3: 1(3)Google 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.