Genome-scale analyses of butanol tolerance in Saccharomyces cerevisiae reveal an essential role of protein degradation
© González-Ramos et al.; licensee BioMed Central Ltd. 2013
Received: 18 January 2013
Accepted: 22 March 2013
Published: 3 April 2013
n-Butanol and isobutanol produced from biomass-derived sugars are promising renewable transport fuels and solvents. Saccharomyces cerevisiae has been engineered for butanol production, but its high butanol sensitivity poses an upper limit to product titers that can be reached by further pathway engineering. A better understanding of the molecular basis of butanol stress and tolerance of S. cerevisiae is important for achieving improved tolerance.
By combining a screening of the haploid S. cerevisiae knock-out library, gene overexpression, and genome analysis of evolutionary engineered n-butanol-tolerant strains, we established that protein degradation plays an essential role in tolerance. Strains deleted in genes involved in the ubiquitin-proteasome system and in vacuolar degradation of damaged proteins showed hypersensitivity to n-butanol. Overexpression of YLR224W, encoding the subunit responsible for the recognition of damaged proteins of an ubiquitin ligase complex, resulted in a strain with a higher n-butanol tolerance. Two independently evolved n-butanol-tolerant strains carried different mutations in both RPN4 and RTG1, which encode transcription factors involved in the expression of proteasome and peroxisomal genes, respectively. Introduction of these mutated alleles in the reference strain increased butanol tolerance, confirming their relevance in the higher tolerance phenotype. The evolved strains, in addition to n-butanol, were also more tolerant to 2-butanol, isobutanol and 1-propanol, indicating a common molecular basis for sensitivity and tolerance to C3 and C4 alcohols.
This study shows that maintenance of protein integrity plays an essential role in butanol tolerance and demonstrates new promising targets to engineer S. cerevisiae for improved tolerance.
KeywordsSaccharomyces cerevisiae Butanol tolerance Evolutionary engineering Deletion collection screening Whole genome sequencing Proteasome Multivesicular bodies
Four-carbon alcohols, including n-butanol, produced from renewable biomass, are promising alternatives to ethanol as biofuels: they are less volatile, less hygroscopic and less corrosive than ethanol, and have higher energy content. They can also be added to gasoline as a fuel additive or even replace it completely without modification of the existing car engines [1, 2] and are compatible with the existing infrastructure for gasoline distribution .
n-Butanol can be naturally produced by some Clostridium species through the so called ABE fermentation, a process that yields acetone, butanol and ethanol . ‘Biobutanol’ has been historically produced through this process , however, it exhibits several drawbacks that limit its current application for large scale production. Clostridium species are strictly anaerobic and slow growing microorganisms . The low n-butanol yield and the production of by-products such as acetone and ethanol also affect butanol production negatively. Although attempts have been made to engineer Clostridium in order to overcome these limitations [3, 5, 6], its complex physiology and lack of efficient genetic tools make it difficult to engineer .
The limitations of clostridial n-butanol production were an incentive to develop and use alternative hosts. To this end, the genes encoding the n-butanol biosynthetic pathway from C. acetobutilicum have been expressed in different microorganisms. In engineered Escherichia coli strains, n-butanol titers of up to 600 mg/l were reached via this strategy [7–9]. Pseudomonas putida, Bacillus subtilis and Lactobacillus brevis have also been used as butanol-producing hosts, although the concentrations obtained were lower than those reached by E. coli[10, 11]. Similarly, S. cerevisiae has been engineered for n-butanol production by substituting the clostridial enzymes by isozymes from different microorganisms, however, the concentration obtained remained very low (2.5 mg/l) . Interestingly, S. cerevisiae can also produce n-butanol from an endogenous pathway that involves norvaline as an intermediate, which may offer new, interesting options for pathway engineering [13–15].
In addition to pathway engineering, product toxicity represents a major challenge in microbial butanol production [16, 17]. As an organic solvent, n-butanol tends to partition into biological membranes, thereby increasing their fluidity and changing their structures [17–19]. Membrane functions are severely hampered in the presence of n-butanol; cells lose the ability to maintain internal pH due to the increased proton permeability of the cytoplasmic membrane and inhibition of the membrane ATPase [20, 21]. The increase in membrane fluidity also results in a loss of intracellular molecules such as proteins, RNA and ATP  and glucose uptake is strongly impeded . In situ product recovery systems allow the efficient removal of n-butanol and thereby reduce toxicity to the cells [5, 22], however, they add substantial capital and operating costs to the process. These costs would be reduced with a higher operating concentration of n-butanol, and thus more tolerant strains would greatly improve process economics [3, 16, 17].
Improvement of n-butanol tolerance has been explored in both E. coli and S. cerevisiae. The overexpression in E. coli of genes related to iron transport and metabolism increase isobutanol tolerance as well as genes involved in membrane functions, amino acid transport, sugar transport and stress response, played a role in tolerance . In another study, genes related to the synthesis of glucosamine-6-phosphate (precursor of peptidoglycan and lipopolysaccharide), multidrug efflux system, degradation of L-cysteine and L-tryptophan, and galactitol metabolism were identified as potential targets for further engineering of butanol tolerance in E. coli. The overexpression in E. coli of GroEL/GroES, molecular chaperones that prevent protein aggregation under stress conditions and assist in protein folding, also resulted in strains with an improved butanol tolerance . However, the levels of tolerance obtained so far in E. coli through different strategies do not match its production capacity and therefore tolerance needs to be further investigated and improved.
S. cerevisiae is known for its high tolerance to alcohols [17, 26] and low pH , which is especially relevant when fermenting lignocellulosic material. Several groups attempted to increase n-butanol tolerance in S. cerevisiae by modifying genes involved in multidrug resistance , cell wall integrity , high osmolarity response , filamentous growth  and amino acid starvation . Some of the modified strains showed increased biomass yield or growth rate relative to their reference in the presence of specific n-butanol concentrations. However, the maximum n-butanol concentration that allows growth in S. cerevisiae has not been increased and the mechanisms of butanol toxicity remain largely unknown.
Experimental design and evaluation of n-butanol tolerance of S. cerevisiae strains BY4741 and CEN.PK113-7D
The reference strain BY4741 displayed a BSI value of 2. Strains that are more sensitive than this reference will have a BSI >2, while strains exhibiting a more tolerant phenotype with 1% n-butanol will have BSI values between 1 and 2.
Screening for genes involved in n-butanol tolerance
To identify the molecular mechanisms involved in n-butanol tolerance in S.cerevisiae, the haploid nonessential gene deletion collection (5154 strains), made in the BY4741 genetic background  was screened for growth in the absence and presence of 1% n-butanol. The BSI values of the entire knockout collection were compared to BY4741, and no strain revealed increased tolerance compared to the reference. After multiple replicate tests (n=8), a total of 105 deletion strains were confirmed to be more sensitive towards n-butanol (Additional file 1: Table S1). The 55 strains displaying BSI values between 2 and 10 were to a large majority still able to grow in the presence of 1% n-butanol (final OD with 1% butanol > 0.1) although they reached a lower OD660 than the reference (Additional file 1: Table S1). The 50 strains exhibiting a BSI value above 10 (Additional file 1: Table S1) were unable to grow in the presence of 1% n-butanol (final OD with 1% butanol < 0.08), indicating a hyper-sensitive phenotype. To exclude specific effects of the BY4741 genetic background, gene deletions resulting in absence of growth at 1% n-butanol were reintroduced in a S. cerevisiae strain of the CEN.PK lineage, with the exception of dubious open reading frames (YKL118W, YPL062W, YBL094C, YLR338W and YDR157W). Additionally, STP11 and DID4 were arbitrarily choosen and tested in CEN.PK113-7D, bringing the total number of tested deletions to 47 (Additional file 1: Table S1). Interestingly, three deletions did not yield transformants (GON7, UAF30 and OCH1). The genes UAF30 and GON7 have already been shown to be essential in the Sigma1278b strain . No such obvious explanation was found for our inability to obtain an och1Δ strain in the CEN.PK genetic background.
Set of 35 genes whose deletion results in n -butanol sensitivity in both BY4741 and CEN.PK113-7D genetic backgrounds
Alpha 3 subunit of the 20S proteasome. PRE9 encodes the only nonessential proteasome subunit.
Chaperone required for correct maturation of the 20S proteasome.
Subunit of the SCF ubiquitin ligase complexes, responsible of recognizing misfolded proteins.
Ubiquitin-specific protease and its cofactor, respectively. Responsible of the deubiquitination of proteins. Important to maintain the pull of free ubiquitin in the cells.
Components of the ESCRT-I, -II and –II complexes, respectively. Responsible of sorting ubiquitinated membrane proteins into Multivesicular Bodies (MVB) for their degradation in the vacuole.
Responsible of deubiquitination in the MBV. Important to maintain the pull of free ubiquitin.
Vps34 and Vps15 form a complex responsible for the synthesis of phosphatidylinositol 3-phosphate, involved in endosomal membrane trafficking and in the regulation of protein sorting.
Component of the conserved R2TP complex. Interacts with Hsp90 to mediate assembly large protein complexes such as box C/D snoRNPs and RNA polymerase II.
Transcription cofactor required for the unfolded protein response.
Subunits of the GET complex. Involved in the insertion of tail anchored proteins into the ER. Tail anchored proteins play a role in vesicular traffic and folding or degradation of membrane proteins.
Subunit F of the eight-subunit V1 peripheral membrane domain of vacuolar H+−ATPase, and peripheral membrane protein that is required for vacuolar H+−ATPase function, respectively.
Protein that regulates myosin function; involved in endocytosis.
6-phosphogluconate dehydrogenase, required for adaptation to oxidative stress.
Subunit of the alpha-1,6 mannosyltransferase complex. Involved in osmotic sensitivity.
Protein of unknown function, required for mitochondrial genome maintenance.
Subunits of the THO complex. Involved in transcription elongation by RNA polymerase II and in telomere maintenance.
Cytoskeletal protein binding protein required for assembly of the cortical actin cytoskeleton.
Epsilon-COP subunit of the coatomer; regulates retrograde Golgi-to-ER protein traffic.
Mitochondrial glutamyl-tRNA synthetase.
Non-essential kinetochore protein, subunit of the Ctf19 central kinetochore complex.
Cytosolic aldehyde dehydrogenase, required for conversion of acetaldehyde to acetate.
Member of the NineTeen Complex. Involved in splicing of nuclear RNAs via the spliceosome.
Regulatory subunit of Glc7p, involved in negative regulation of glucose-repressible genes.
Component of the RSC chromatin remodeling complex. involved in telomere maintenance.
Third subunit of DNA polymerase delta, involved in chromosomal DNA replication.
Cytoplasmic DExD/H-box helicase, stimulates mRNA decapping.
Actin-associated protein involved in cytoskeletal organization and cytokinesis.
Aspartic beta semi-aldehyde dehydrogenase, catalyzes the second step in the common pathway for methionine and threonine biosynthesis.
Evolutionary engineering for increased butanol tolerance
To further characterize the evolved mutants, we tested and compared the growth of the most tolerant strain, IMS0351, with that of the parental strain CEN.PK113-7D in anaerobic bioreactor cultures grown on synthetic medium with 1.5% n-butanol. Under these conditions, CEN.PK113-7D was unable to grow, while IMS0351 exhibited a growth rate of 0.085 h-1 (Figure 6C), thereby confirming the increased tolerance of the evolved strain. Although the n-butanol concentration decreased slightly due to the nitrogen sparging in these cutlures, the butanol concentration remained above 1.38% throughout the fermentation (Figure 6C).
Tolerance of the evolved strains to different alcohols
Identification and reverse engineering of the mutations present in evolved butanol-tolerant strains
Single nucleotide variation identified in the evolved strains IMS0344 and IMS0351 by whole-genome resequencing
Amino acid change
change from G to C 296 bp upstream of RPL10,
change from T to C 984 bp upstream of SST4.
The only gene identified in the screening of the deletion library whose overexpression led to increased butanol tolerance was YLR224W, which encodes one of 22 F-box proteins in S. cerevisiae (http://www.yeastgenome.org/). Together with Skp1, Cdc53 and Rbx1, individual F-box proteins form SCF ubiquitin-ligase complexes involved in ubiquitin-dependent protein degradation. F-box proteins contribute to the specificity of SCF complexes by, independently of the complex, aggregating to target proteins and subsequently binding to Skp1 [37, 62]. With the exception of Met30 and Cdc4, all S. cerevisiae F-box protein deletion mutants are viable. However, in our screen, YLR224W was the only gene encoding an F-box protein whose deletion led to increased butanol sensitivity. This suggests that Ylr224W specifically targets a (subset of) protein(s) whose degradation is important for butanol tolerance. A similar conclusion was drawn based on a previous observation that overexpression of YLR224W conferred increased tolerance to methylmercury . Overexpression of genes whose products are involved in proteasome-driven protein degradation, but act downstream of YLR224W, did not result in improved tolerance (Figure 10). This suggested that ubiquitination of YLR224W targets is the rate-limiting step in their degradation. Consistent with this notion, the effect of the overexpression of YLR224W on methylmercury tolerance was suppressed by the proteasome inhibitor MG132. . Identification of the protein targets of Ylr224W should provide additional insight in the molecular mechanisms of butanol sensitivity and tolerance in S. cerevisiae.
Screening of the haploid deletion collection of S. cerevisiae was limited to an exploration of the impact of loss-of-function mutations in non-essential genes on butanol tolerance. The laboratory evolution strategy allowed inspection of a much wider genetic landscape, which was not limited to non-essential genes but also widely expanding the type of mutations (i.e. nucleotide variation, small insertion-deletion and duplication ) resulting in proteins with modified catalytic or structural activities. This approach is especially relevant for the exploration of central cellular processes such as protein degradation, in which a substantial fraction of the genes involved cannot be studied via haploid deletion strains because their protein products are essential for growth.
In two independently evolved butanol-tolerant strains, mutations in RPN4 and RTG1 were shown to contribute to increased butanol tolerance. RPN4 encodes a transcription factor involved in the expression of proteasome genes  and is required for normal levels of intracellular proteolysis. Rpn4 is an unstable protein characterized by an extremely short half life . It is regulated by a negative feedback circuit; Rpn4 promotes expression of proteasome genes and is itself subject to proteasomal degradation . The distinct RPN4 alleles present in the two evolution lines might, for example, encode Rpn4 versions that are less susceptible to degradation by the proteasome or that are more active transcription factors. As a result of this change, the transcription of proteasome genes would be increased or maintained in time, thus increasing or extending the activity of the proteasome. Such a mechanism would be entirely consistent with the important role of protein degradation in butanol tolerance revealed by the screening of yeast deletion mutants.
RTG1 encodes a transcription factor involved in the communication between mitochondria and the nucleus  and it is also required for the expression of genes encoding peroxisomal proteins . Reverse engineering of the rtg1-1 allele in CEN.PK113-7D resulted in a strain with an increased butanol tolerance. Additionally, the strain defective in RTG1 in the knockout collection showed an increased butanol sensitivity relative to BY4741 (Additional file 1: Table S1), indicating its relevance for tolerance. However, the mechanism by which the mutated alleles present in the evolved strains confer tolerance is unclear and will necessitate further dedicated experimental work.
This study has identified several defined mutations that can be applied to increase tolerance to propanol and butanol isomers in S. cerevisiae. While the resulting significant increases in tolerance are modest from the viewpoint of industrial application, the methodology used in this study should be directly applicable for further improvement of tolerance to these and other compounds. In particular, this study demonstrates how combining whole-genome sequencing and classical yeast genetics can accelerate the analysis and reverse engineering of strains generated by laboratory evolution. Over the past few years, molecular characterization of evolutionary engineered yeast strains has been strongly facilitated by the fast developments in Next Generation Sequencing methodologies [35, 40, 41, 68]. However, establishing which of the observed mutations contribute to the selected phenotype still remained a laborious process, in which mutations have to be reversed engineered in a naive strain, either individually or in combination. Several solutions have been proposed to simplify the detection of biologically relevant mutations. First, availability of a high-quality reference genome sequence for the parental strain greatly contributes to reducing the number of false positives resulting from sequence comparison with evolved strains . In this study, direct comparison of the evolved strains to the assembled genome of the non-evolved parent led to fewer than 10 mutation calls in each evolutionary run. Secondly, the use of independent, parallel evolution lines may rapidly identify genes or processes that are affected by mutation in independent experiments and are therefore likely to contribute to the observed phenotype [35, 40, 69]. In this study, we used backcrossing, a classical approach in yeast genetics, to enrich and isolate mutations that contributed to butanol tolerance in evolved strains. The tremendous power of this approach is evident from the fact that, after only three cycles of back crossing, the selected F3 segregant IMS0344 retained only three mutations of which, after reverse engineering, two were shown to directly contribute to butanol tolerance of the original evolved strain.
Screening of a deletion mutant library of S. cerevisiae and analysis of butanol-tolerant strains isolated by whole genome sequencing revealed that protein degradation is a key process in the tolerance of S. cerevisiae to n-butanol. Mutations in three genes were shown to significantly increase butanol tolerance and can be applied in metabolic engineering strategies for butanol production with this yeast. Combination of whole-genome sequencing of evolved strains with repeated backcrossing to a native, isogenic strain was shown to be an extremely powerful approach to identify biologically relevant mutations. Further studies on the molecular mechanisms by which protein turnover affects butanol tolerance should provide further leads for strain improvement.
Strains, media and growth conditions
Saccharomyces cerevisiae strains used in this study
MATa his3Δ1 leu2Δ met15Δ ura3Δ
Yeast KO collection
MATa his3Δ1 leu2Δ met15Δ ura3Δ ORFΔ::KanMX
MATa ORF pr::KanMX-TPI1 pr::ORF
MATa rpn4-1 rtg1-1 ubr1-1 ssk2-1
MATa/MATα RPN4/rpn4-1 RTG1/rtg1-1 UBR1/ubr1-1 SSK2/ssk2-1 URA3/ura3Δ::KanMX
MATα rpn4-1 rtg1-1 ubr1-1 ssk2-1
MATa/MATα RPN4/rpn4-1 RTG1/rtg1-1 UBR1/ubr1-1 SSK2/ssk2-1 URA3/ura3Δ::KanMX
MATα rpn4-1 rtg1-1 ubr1-1 ssk2-1
MATa/MATα RPN4/rpn4-1 RTG1/rtg1-1 UBR1/ubr1-1 SSK2/ssk2-1 URA3/ura3Δ::KanMX
MATa rpn4-1 rtg1-1 ssk2-1
MATa rpn4-2 rtg1-2 ubr1-2 nma111-2 rpl10-2 sto1-2 sst4-2
MATa/MATα RPN4/rpn4-2 RTG1/rtg1-2 UBR1/ubr1-2 NMA111/nma111-2 RPL10/rpl10-2 STO1/sto1-2 SST4/sst4-2 URA3/ura3Δ::KanMX
MATa rpn4-2 rtg1-2 nma111-2 rpl10-2
MATa/MATα RPN4/rpn4-1 RTG1/rtgi-2 NMA111/nma111-2 RPL10/rpl10-2 URA3/ura3Δ::KanMX
MATa rpn4-2 rtg1-2 nma111-2 rpl10-2
MATa/MATα RPN4/rpn4-2 RTG1/rtg1-2 NMA111/nma111-2 RPL10/rpl10-2 URA3/ura3Δ::KanMX
MATa rpn4-2 rtg1-2 nma111-2
MATa RPN4 pr::KanMX RPN4 pr -rpn4-1
MATa RTG1 pr::KanMX RTG1 pr -rtg1-1
The crossing studies were done in solid synthetic medium  (containing 20 g/l agar) without uracil and supplemented with G418 (InvivoGen) at a concentration of 200 μg/ml. For sporulation of diploid strains, they were grown in YPA (yeast extract 10 g/l, peptone 20 g/l, potassium acetate 10 g/l) and transferred to the sporulation medium (20 g/l potassium acetate, pH 7.0) .
Screening of the deletion collection
Analysis of butanol tolerance in 96 well plates
The strains were grown in YPD for 24 hours at 30°C and 200 rpm. The OD660 of the cultures was measured and they were diluted to an OD660 of 5 with fresh YPD. 96 well plates were filled with 200 μl/well of synthetic medium containing butanol concentrations ranging from 0 to 1.9% (12 different concentrations, 8 replicas per condition). The plates containing synthetic medium with butanol were inoculated from the diluted YPD cultures with the use of a 96 well pin replicator. The plates were sealed with a gas impermeable tape (NUNC) and incubated at 30°C for 48 hours without shaking. After the incubation time, the cells were resuspended in an orbital shaker (MS3, IKA-Werke GmbH & Co), the sealing tape was removed and the OD660 of each well was measured in a GENIos Pro microplate spectrophotometer (Tecan). At least 2 plates were inoculated per strain and for each butanol concentration the average OD660 was calculated with the data from both plates (16 replicas).
A 96 well plate setup was also used to compare the growth of two different strains in the presence of increasing concentrations of butanol. The strains were grown in YPD for 24 hours at 30°C and 200 rpm. From these cultures, 12 aliquots of 1 ml of synthetic medium with butanol concentrations ranging from 0% to 1.9% (12 different concentrations) were inoculated at an OD660 of 0.1. Each inoculated aliquot was distributed in 4 wells of a 96 well plate (200 μl per well). The plates were sealed with a gas impermeable sealing tape (NUNC) and incubated at 30°C in a GENIos Pro microplate spectrophotometer (Tecan) for 48 hours. A measurement of the OD660 of each well was taken every 30 minutes.
Anaerobic serial transfer for evolutionary engineering
Two independent sequential batch cultures of the strain CEN.PK113-7D were cultivated in synthetic medium  containing 1% of butanol. The cultivation was carried out in shake flask, using closed bottles to prevent butanol evaporation. After each batch, the OD660 of the cultures was measured and a new shake flask was inoculated from the previous one to an OD660 of 0.1. When a higher final OD660 was observed, the butanol concentration of the following batch was increased. In the first evolution line the butanol concentration was increased up to 1.25% after 55 batches. After 83 batches a colony was isolated and its butanol tolerance tested. In the second evolution line the butanol concentration was increased up to 1.2% after 30 batches, to 1.3% after 42 batches and to 1.35% after 52 batches. After 63 batches a colony was isolated and its butanol tolerance tested. The strains evolved for butanol tolerance in the first and second cultures are IMS0351 and IMS0351, respectively.
Batch fermentation in bioreactor in the presence of butanol and isobutanol
Each strain was grown for 24 hours in synthetic medium containing 2.3 g/l urea as the sole nitrogen source and 1% n-butanol or isobutanol, respectively. Cultivation was performed in 2 l bioreactors (Applikon, Schiedam, The Netherlands) containing 1 l of synthetic medium with 20 g/l glucose, and supplemented with 0.01 g/l ergosterol and 0.42 g/l Tween 80 for anaerobic cultivation. Antifoam Emulsion C (Sigma-Aldrich, Zwijndrecht, The Netherlands) was prepared as a 20% (w/v) solution, sterilized and added into the bioreactor at a final concentration of 0.2 g/l. The bioreactors were inoculated at an OD600 of 0.1 and after inoculation; the right concentration of n-butanol or isobutanol was added. Cultures were stirred at 800 rpm and kept anaerobic by flowing N2 gas in the head space at a flow-rate of 0.1 l/min. The pH of the culture was kept at 5 by automated addition of 2.0 M KOH. Samples were taken over time to measure the OD of the cultures and determine butanol and glucose concentrations. Glucose, ethanol, glycerol and butanol concentrations were analysed via HPLC using an Aminex HPX-87H ion exchange column operated at 60°C with 5 mM H2SO4 as mobile phase at a flow rate of 0.6 ml/min.
Construction of the deletion and overexpression strains
In order to delete the target genes in CEN.PK113-7D, the deletion cassettes were prepared by direct amplification of the deleted gene in the corresponding strain from the BY knockout collection (Additional file 1: Table S2). The genomic DNA of the strains from the collection was extracted with the YeaStar™ Genomic DNA kit (Zymo Research, Irvine, CA) following the instructions of the manufacturer, and used as a template to amplify deletion cassettes. The deletion cassettes were amplified with the Phusion high-fidelity DNA polymerase (Thermo Scientific, Landsmeer, The Netherlands), with primers annealing at least 100 base pairs upstream and downstream of the KanMX marker, corresponding to the promoter and terminator regions of the genes, respectively (Additional file 1: Table S3). The deletion cassettes were transformed in CEN.PK113-7D by the lithium acetate method , inserted in the genome by homologous recombination and the selection of transformants was done in YPD plates containing 200 μg/ml G418. The confirmation of the deletion was done by PCR amplification with DreamTaq DNA Polymerase (Thermo Scientific), with a forward primer that anneals upstream of the insertion point and the reverse primer KanMX-DCR that anneals in the KanMX gene. In some cases, the confirmation was done with the forward primer KanMX-DCF (annealing in the KanMX gene) and a reverse primer annealing downstream of the insertion point. The amplification with these two primer pairs is only possible if the ORF of the gene has been successfully replaced by KanMX.
The overexpression of the genes was done by inserting upstream of their ORF the strong constitutive promoter from the gene TPI1. Overexpression cassettes containing the KanMX marker and the TPI1 promoter were amplified from the plasmid pUG6-TPI1 prom, with the Phusion high-fidelity DNA polymerase (Thermo Scientific) and the primer pair named with the suffixes –OF and –OR, stating for overexpression forward and reverse, respectively (Additional file 1: Table S4). The primers –OF contain 50 base pairs homologous to the promoter region of each gene and the –OR ones, homologous to the beginning of the ORF. The overexpression cassettes were transformed in CEN.PK113-7D by the lithium acetate method , inserted in the genome by homologous recombination and the selection of transformants was done in YPD plates containing G418. The confirmation of the insertion was done by PCR amplification DreamTaq DNA Polymerase (Thermo Scientific) and the forward primer TPI1prom-ICF, annealing in the promoter of TPI1 and the reverse primer named with the suffix –OCR, annealing in the ORF of each gene. The amplification with these two primers is only possible if the promoter from TPI1 has been successfully inserted upstream of the ORF. In both the confirmation of the deletion and overexpression of the genes, a negative control PCR was done with the primer pair –OCF and –DCR that amplifies only in the native allele.
Reverse engineering of the mutated rpn4-1 and rtg1-1 alleles
The gDNA of the strain IMS0344 was extracted with the YeaStar™ Genomic DNA kit (Zymo Research) following the instructions of the manufacturer, and used as a template to amplify by PCR the alleles rpn4-1 and rtg1-1. For the reverse engineering of the mutated alleles in CEN.PK113-7D, two different cassettes were used. One cassette consists of the marker gene KanMX, flanked in one side by 50 base pairs homologous to a region upstream of the target gene, and in the other side by a synthetic sequence to promote homologous recombination (SHR). The second cassette consists of the mutated gene flanked in one side by the same SHR sequence used in the first cassette (Additional file 1: Table S5). The primers used for the amplification of the cassettes and for the analysis of the transformants are listed in Additional file 1: Table S6. When cotransforming both cassettes into CEN.PK113-7D, they recombine through the SHR sequence creating one single cassette. The new cassette is inserted into the corresponding locus in the genome by homologous recombination, replacing the original allele by the mutated one. The insertion of the cassette in the right locus was checked by PCR (Additional files 1: Tables S5 and S6), where a PCR product is only obtained if the cassette has been successfully inserted. The replacement of the original allele by the mutated one was checked by Real Time-PCR coupled with a High Resolution Melting analysis of the amplified product. The primers amplify a fragment of a size between 80 to 200 base pairs that includes the region with the SNV. The PCR mix was made by using the Type-it® HRM™ kit following the instructions of the manufacturer and the amplification was done into a Rotor-Gene Q thermocycler (Qiagen, Venlo, The Netherlands). After the amplification step, the temperature was increased from 65°C to 95°C in steps of 0.1°C and the fluorescence was recorded. The fluorescence profile of the samples was compared with the ones from IMS0344 (positive control) and CEN.PK113-7D (negative control).
Crossing, sporulation and screening for butanol tolerance of the haploid segregants
Two ancestor strains with different mating types were used for the crossing; IMK439 (MATa) and IMK440 (MATα), derived from CEN.PK113-7D and CEN.PK113-1A respectively, in which the URA3 locus has been replaced by the KanMX marker. The cassettes for the deletion of URA3 was amplified from the plasmid pUG6-TPI1 prom with the primers URA3-KanMXF and URA3-KanMXR, and transformed into CEN.PK113-7D and CEN.PK113-1A. The transformants were selected in YPD plates containing G418 and the deletion of URA3 was confirmed by their inability of growing in the absence of uracil. This way, when crossing the evolved strains (MATa) together with IMK440 (MATα) in a plate containing synthetic medium without uracil and with G418, only diploid cells are able to grow. The diploid cells were sporulated, their spores selected and plated on synthetic medium without uracil, selecting only the haploid segregants containing the native URA3 locus . The haploid segregants were screened for butanol tolerance and the one with the same tolerance as the evolved strain was selected. The selected haploid segregant was crossed with IMK439 or IMK440 (depending on its mating type) and a new sporulation and spore isolation was performed.
Sequencing and analysis of the sequences
Genomic DNA of the strains CEN.PK113-7D, CEN.PK113-1A, IMS0344, IMS0346, IMS0348, IMS0350, IMS0351, IMS0353, IMS0355 and IMS0357 was prepared as described previously . Libraries of 350-bp insert were constructed and paired end sequenced (100 base pair reads) using an illlumina HISeq 2000 sequencer (Baseclear BV, Leiden, The Netherlands). A minimum data quantity of 950 Mb was generated for each strain, representing a minimum 80-fold coverage. The sequence reads were mapped onto CEN.PK113-7D genome  using Burrows–Wheeler Alignment tool (BWA) and further processed using SAMtools [76–78]. Single-nucleotide variations were extracted from the mapping using SAMtools’ varFilter. Default settings were used, except that the maximum read depth was set to 400X (−D400). To minimize false positive mutation calls, custom Perl scripts were used for further mutation filtering: i) mutation calls containing ambiguous bases in mapping consensus were filtered out, ii) only the single-nucleotide variations with a quality of at least 20 were kept (variant quality is defined as the Phred-scaled probability that the mutation call is incorrect [79, 80], iii) mutations with a. depth of coverage < 10X were discarded and iv) the mutations found in CEN.PK113-1A were subtracted from the list sequence Eventually, the single nucleotide variations were physically positioned and functionally annotated according to the CEN.PK113-7D sequence annotation . The raw sequencing data of strains CEN.PK113-1A, IMS0344, IMS0346, IMS0348, IMS0350, IMS0351, IMS0353, IMS0355 and IMS0357 have been deposited as sequence read archives (SRA, http://www.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?) under bioproject ID PRJNA191134.
Fisher’s exact test
This work was funded by the Kluyver Centre for Genomics of Industrial Fermentations. We are grateful to Filipa Dias Valada, who contributed to the construction and characterization of strains, to Erik de Hulster, for his help with the fermentations and to Bart Oud for his help in the analysis of the genome sequencing data. The research group of JTP is part of the Kluyver Centre for Genomics of Industrial Fermentation (http://www.kluyvercentre.nl/ ) is supported by The Netherlands Genomics Initiative (http://www.genomics.nl/). JTP and J-MD were also supported by the “Platform Green Synthetic Biology” programme (http://www.pgsb.nl/) funded by NGI. This research was also supported by BIOFLAVOUR, COST Action FA0907 (http://www.bioflavour.insa-toulouse.fr).
- Durre P: Biobutanol: an attractive biofuel. Biotechnol J 2007, 2:1525–1534.View Article
- Fischer CR, Klein-Marcuschamer D, Stephanopoulos G: Selection and optimization of microbial hosts for biofuels production. Metab Eng 2008, 10:295–304.View Article
- Lee SY, Park JH, Jang SH, Nielsen LK, Kim J, Jung KS: Fermentative butanol production by Clostridia . Biotechnol Bioeng 2008, 101:209–228.View Article
- Jones DT, Woods DR: Acetone-butanol fermentation revisited. Microbiol Rev 1986, 50:484–524.
- Zheng YN, Li LZ, Xian M, Ma YJ, Yang JM, Xu X, He DZ: Problems with the microbial production of butanol. J Ind Microbiol Biotechnol 2009, 36:1127–1138.View Article
- Durre P: Fermentative production of butanol–the academic perspective. Curr Opin Biotechnol 2011, 22:331–336.View Article
- Inui M, Suda M, Kimura S, Yasuda K, Suzuki H, Toda H, Yamamoto S, Okino S, Suzuki N, Yukawa H: Expression of Clostridium acetobutylicum butanol synthetic genes in Escherichia coli . Appl Microbiol Biotechnol 2008, 77:1305–1316.View Article
- Atsumi S, Cann AF, Connor MR, Shen CR, Smith KM, Brynildsen MP, Chou KJ, Hanai T, Liao JC: Metabolic engineering of Escherichia coli for 1-butanol production. Metab Eng 2008, 10:305–311.View Article
- Shen CR, Lan EI, Dekishima Y, Baez A, Cho KM, Liao JC: Driving forces enable high-titer anaerobic 1-butanol synthesis in Escherichia coli . Appl Environ Microbiol 2011, 77:2905–2915.View Article
- Nielsen DR, Leonard E, Yoon SH, Tseng HC, Yuan C, Prather KL: Engineering alternative butanol production platforms in heterologous bacteria. Metab Eng 2009, 11:262–273.View Article
- Berezina OV, Zakharova NV, Brandt A, Yarotsky SV, Schwarz WH, Zverlov VV: Reconstructing the clostridial n-butanol metabolic pathway in Lactobacillus brevis . Appl Microbiol Biotechnol 2010, 87:635–646.View Article
- Steen EJ, Chan R, Prasad N, Myers S, Petzold CJ, Redding A, Ouellet M, Keasling JD: Metabolic engineering of Saccharomyces cerevisiae for the production of n-butanol. Microb Cell Fact 2008, 7–36.
- Guymon JF, Ingraham JL, Crowell EA: The formation of n-propyl alcohol by Saccharomyces cerevisiae . Arch Biochem Biophys 1961, 95:163–168.View Article
- Ingraham JL, Guymon JF, Crowell EA: The pathway of formation of n-butyl and n-amyl alcohols by a mutant strain of Saccharomyces cerevisiae . Arch Biochem Biophys 1961, 95:169–175.View Article
- Romagnoli G, Luttik MA, Kotter P, Pronk JT, Daran JM: Substrate specificity of thiamine pyrophosphate-dependent 2-oxo-acid decarboxylases in Saccharomyces cerevisiae . Appl Environ Microbiol 2012, 78:7538–7548.View Article
- Knoshaug EP, Zhang M: Butanol tolerance in a selection of microorganisms. Appl Biochem Biotechnol 2009, 153:13–20.View Article
- Liu S, Qureshi N: How microbes tolerate ethanol and butanol. N Biotechnol 2009, 26:117–121.View Article
- Vollherbst-Schneck K, Sands JA, Montenecourt BS: Effect of butanol on lipid composition and fluidity of Clostridium acetobutylicum ATCC 824. Appl Environ Microbiol 1984, 47:193–194.
- Sardessai Y, Bhosle S: Tolerance of bacteria to organic solvents. Res Microbiol 2002, 153:263–268.View Article
- Bowles LK, Ellefson WL: Effects of butanol on Clostridium acetobutylicum . Appl Environ Microbiol 1985, 50:1165–1170.
- Ingram LO: Microbial tolerance to alcohols: role of the cell membrane. Trends Biotechnol 1986, 4:40–44.View Article
- Ennis BM, Marshall CT, Maddox IS, Paterson AHJ: Continuous product recovery by in-situ gas stripping/condensation during solvent production from whey permeate using Clostridium acetobutylicum . Biotechnol Lett 1986, 8:725–730.View Article
- Reyes LH, Almario MP, Kao KC: Genomic library screens for genes involved in n-butanol tolerance in Escherichia coli . PLoS One 2011, 6:e17678.View Article
- Atsumi S, Wu TY, Machado IM, Huang WC, Chen PY, Pellegrini M, Liao JC: Evolution, genomic analysis, and reconstruction of isobutanol tolerance in Escherichia coli . Mol Syst Biol 2010, 6:449.View Article
- Zingaro KA, Papoutsakis TE: GroESL overexpression imparts Escherichia coli tolerance to i-, n-, and 2-butanol, 1,2,4-butanetriol and ethanol with complex and unpredictable patterns. Metab Eng 2012, 15:196–205.View Article
- Casey GP, Ingledew WM: Ethanol tolerance in yeasts. Crit Rev Microbiol 1986, 13:219–280.View Article
- de Lucena RM, Elsztein C, Simoes DA, de Morais MAJ: Participation of CWI, HOG and Calcineurin pathways in the tolerance of Saccharomyces cerevisiae to low pH by inorganic acid. J Appl Microbiol 2012, 113:629–640.View Article
- van Dyk TK: New recombinant yeast cell useful for producing butanol, comprises at least one genetic modification that reduces pleiotropic drug resistant 5 phenotype activity. 2010. WO2010099524-A1
- Bramucci MG, Larossa RA, Smulski DR: New recombinant yeast cell useful for producing 1-butanol, 2-butanol and isobutanol, comprises butanol biosynthetic pathway and genetic modification increasing cell wall integrity pathway activity. 2010, US2010167364-A1.
- Bramucci MG, Larossa RA, Smulski DR: New recombinant yeast cell useful for producing 1-butanol, 2-butanol and isobutanol, comprises butanol biosynthetic pathway and genetic modification increasing activity of high osmolarity/glycerol response pathway. 2010, US2010167365-A1.
- Bramucci MG, Larossa RA, Singh M: New recombinant yeast cell useful for producing 1-butanol, 2-butanol and isobutanol, comprises butanol biosynthetic pathway and genetic modification which increases activity of nitrogen starvation-induced filamentous growth response. 2010, US2010167363-A1.
- Larossa RA: New recombinant yeast host cell comprises a genetic modification which reduces the response in the general control response to amino acid starvation, useful for producing butanol. 2009, US2009280546-A1.
- Winzeler EA, Shoemaker DD, Astromoff A, Liang H, Anderson K, Andre B, Bangham R, Benito R, Boeke JD, Bussey H, Chu AM, Connelly C, Davis K, Dietrich F, Dow SW, El BM, Foury F, Friend SH, Gentalen E, Giaever G, Hegemann JH, Jones T, Laub M, Liao H, Liebundguth N, Lockhart DJ, Lucau-Danila A, Lussier M, M’Rabet N, Menard P, Mittmann M, Pai C, Rebischung C, Revuelta JL, Riles L, Roberts CJ, Ross-MacDonald P, Scherens B, Snyder M, Sookhai-Mahadeo S, Storms RK, Veronneau S, Voet M, Volckaert G, Ward TR, Wysocki R, Yen GS, Yu K, Zimmermann K, Philippsen P, Johnston M, Davis RW: Functional characterization of the S. cerevisiae genome by gene deletion and parallel analysis. Science 1999, 285:901–906.View Article
- Hazelwood LA, Walsh MC, Pronk JT, Daran JM: Involvement of vacuolar sequestration and active transport in tolerance of Saccharomyces cerevisiae to hop iso-alpha-acids. Appl Environ Microbiol 2010, 76:318–328.View Article
- Oud B, van Maris AJ, Daran JM, Pronk JT: Genome-wide analytical approaches for reverse metabolic engineering of industrially relevant phenotypes in yeast. FEMS Yeast Res 2012, 12:183–196.View Article
- Dowell RD, Ryan O, Jansen A, Cheung D, Agarwala S, Danford T, Bernstein DA, Rolfe PA, Heisler LE, Chin B, Nislow C, Giaever G, Phillips PC, Fink GR, Gifford DK, Boone C: Genotype to phenotype: a complex problem. Science 2010, 328:469–469.View Article
- Kus BM, Caldon CE, Andorn-Broza R, Edwards AM: Functional interaction of 13 yeast SCF complexes with a set of yeast E2 enzymes in vitro. Proteins 2004, 54:455–467.View Article
- Hwang GW, Ishida Y, Naganuma A: Identification of F-box proteins that are involved in resistance to methylmercury in Saccharomyces cerevisiae . FEBS Lett 2006, 580:6813–6818.View Article
- Sauer U: Evolutionary engineering of industrially important microbial phenotypes. Adv Biochem Eng Biotechnol 2001, 73:129–169.
- de Kok S, Nijkamp JF, Oud B, Roque FC, de RD, Daran JM, Pronk JT, van Maris AJ: Laboratory evolution of new lactate transporter genes in a jen1Delta mutant of Saccharomyces cerevisiae and their identification as ADY2 alleles by whole-genome resequencing and transcriptome analysis. FEMS Yeast Res 2012, 12:359–374.View Article
- Hong KK, Nielsen J: Metabolic engineering of Saccharomyces cerevisiae : a key cell factory platform for future biorefineries. Cell Mol Life Sci 2012, 69:2671–2690.View Article
- Koppram R, Albers E, Olsson L: Evolutionary engineering strategies to enhance tolerance of xylose utilizing recombinant yeast to inhibitors derived from spruce biomass. Biotechnol Biofuels 2012, 5:32.View Article
- Nijkamp JF, van den Broek M, Datema E, de KS, Bosman L, Luttik MA, Daran-Lapujade P, Vongsangnak W, Nielsen J, Heijne WH, Klaassen P, Paddon CJ, Platt D, Kotter P, van Ham RC, Reinders MJ, Pronk JT, de RD, Daran JM: De novo sequencing, assembly and analysis of the genome of the laboratory strain Saccharomyces cerevisiae CEN.PK113–7D, a model for modern industrial biotechnology. Microb Cell Fact 2012, 11:36.View Article
- Emori Y, Tsukahara T, Kawasaki H, Ishiura S, Sugita H, Suzuki K: Molecular cloning and functional analysis of three subunits of yeast proteasome. Mol Cell Biol 1991, 11:344–353.
- Hochstrasser M: Ubiquitin-dependent protein degradation. Annu Rev Genet 1996, 30:405–439.View Article
- Maurizi MR: Proteasome assembly: biting the hand. Curr Biol 1998, 8:R453-R456.View Article
- Ramos PC, Hockendorff J, Johnson ES, Varshavsky A, Dohmen RJ: Ump1p is required for proper maturation of the 20S proteasome and becomes its substrate upon completion of the assembly. Cell 1998, 92:489–499.View Article
- Rothman JH, Howald I, Stevens TH: Characterization of genes required for protein sorting and vacuolar function in the yeast Saccharomyces cerevisiae. EMBO J 1989, 8:2057–2065.
- Babst M, Odorizzi G, Estepa EJ, Emr SD: Mammalian tumor susceptibility gene 101 (TSG101) and the yeast homologue, Vps23p, both function in late endosomal trafficking. Traffic 2000, 1:248–258.View Article
- Katzmann DJ, Babst M, Emr SD: Ubiquitin-dependent sorting into the multivesicular body pathway requires the function of a conserved endosomal protein sorting complex. ESCRT-I. Cell 2001, 106:145–155.View Article
- Babst M, Katzmann DJ, Snyder WB, Wendland B, Emr SD: Endosome-associated complex, ESCRT-II, recruits transport machinery for protein sorting at the multivesicular body. Dev Cell 2002, 3:283–289.View Article
- Babst M, Katzmann DJ, Estepa-Sabal EJ, Meerloo T, Emr SD: Escrt-III: an endosome-associated heterooligomeric protein complex required for mvb sorting. Dev Cell 2002, 3:271–282.View Article
- Schu PV, Takegawa K, Fry MJ, Stack JH, Waterfield MD, Emr SD: Phosphatidylinositol 3-kinase encoded by yeast VPS34 gene essential for protein sorting. Science 1993, 260:88–91.View Article
- Stack JH, Herman PK, Schu PV, Emr SD: A membrane-associated complex containing the Vps15 protein kinase and the Vps34 PI 3-kinase is essential for protein sorting to the yeast lysosome-like vacuole. EMBO J 1993, 12:2195–2204.
- Stack JH, DeWald DB, Takegawa K, Emr SD: Vesicle-mediated protein transport: regulatory interactions between the Vps15 protein kinase and the Vps34 PtdIns 3-kinase essential for protein sorting to the vacuole in yeast. J Cell Biol 1995, 129:321–334.View Article
- Wendland B, McCaffery JM, Xiao Q, Emr SD: A novel fluorescence-activated cell sorter-based screen for yeast endocytosis mutants identifies a yeast homologue of mammalian eps15. J Cell Biol 1996, 135:1485–1500.View Article
- Kihara A, Noda T, Ishihara N, Ohsumi Y: Two distinct Vps34 phosphatidylinositol 3-kinase complexes function in autophagy and carboxypeptidase Y sorting in Saccharomyces cerevisiae . J Cell Biol 2001, 152:519–530.View Article
- Strahl T, Thorner J: Synthesis and function of membrane phosphoinositides in budding yeast, Saccharomyces cerevisiae. Biochim Biophys Acta 2007, 1771:353–404.View Article
- Fujita K, Matsuyama A, Kobayashi Y, Iwahashi H: The genome-wide screening of yeast deletion mutants to identify the genes required for tolerance to ethanol and other alcohols. FEMS Yeast Res 2006, 6:744–750.View Article
- Auesukaree C, Damnernsawad A, Kruatrachue M, Pokethitiyook P, Boonchird C, Kaneko Y, Harashima S: Genome-wide identification of genes involved in tolerance to various environmental stresses in Saccharomyces cerevisiae . J Appl Genet 2009, 50:301–310.View Article
- Finley D, Ulrich HD, Sommer T, Kaiser P: The ubiquitin-proteasome system of Saccharomyces cerevisiae . Genetics 2012, 192:319–360.View Article
- Kato M, Kito K, Ota K, Ito T: Remodeling of the SCF complex-mediated ubiquitination system by compositional alteration of incorporated F-box proteins. Proteomics 2010, 10:115–123.View Article
- Mannhaupt G, Schnall R, Karpov V, Vetter I, Feldmann H: Rpn4p acts as a transcription factor by binding to PACE, a nonamer box found upstream of 26S proteasomal and other genes in yeast. FEBS Lett 1999, 450:27–34.View Article
- Xie Y, Varshavsky A: RPN4 is a ligand, substrate, and transcriptional regulator of the 26S proteasome: a negative feedback circuit. Proc Natl Acad Sci USA 2001, 98:3056–3061.View Article
- Wang X, Xu H, Ju D, Xie Y: Disruption of Rpn4-induced proteasome expression in Saccharomyces cerevisiae reduces cell viability under stressed conditions. Genetics 2008, 180:1945–1953.View Article
- Rothermel BA, Shyjan AW, Etheredge JL, Butow RA: Transactivation by Rtg1p, a basic helix-loop-helix protein that functions in communication between mitochondria and the nucleus in yeast. J Biol Chem 1995, 270:29476–29482.View Article
- Chelstowska A, Butow RA: RTG genes in yeast that function in communication between mitochondria and the nucleus are also required for expression of genes encoding peroxisomal proteins. J Biol Chem 1995, 270:18141–18146.View Article
- Otero JM, Vongsangnak W, Asadollahi MA, Olivares-Hernandes R, Maury J, Farinelli L, Barlocher L, Osteras M, Schalk M, Clark A, Nielsen J: Whole genome sequencing of Saccharomyces cerevisiae : from genotype to phenotype for improved metabolic engineering applications. BMC Genomics 2010, 11:723.View Article
- Oud B, Flores CL, Gancedo C, Zhang X, Trueheart J, Daran JM, Pronk JT, van Maris AJ: An internal deletion in MTH1 enables growth on glucose of pyruvate-decarboxylase negative, non-fermentative Saccharomyces cerevisiae . Microb Cell Fact 2012, 11:131.View Article
- Brachmann CB, Davies A, Cost GJ, Caputo E, Li J, Hieter P, Boeke JD: Designer deletion strains derived from Saccharomyces cerevisiae S288C: a useful set of strains and plasmids for PCR-mediated gene disruption and other applications. Yeast 1998, 14:115–132.View Article
- van Dijken JP, Bauer J, Brambilla L, Duboc P, Francois 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 N, Rizzi M, Steensma HY, Verrips CT, Vindelov J, Pronk JT: An interlaboratory comparison of physiological and genetic properties of four Saccharomyces cerevisiae strains. Enzyme Microb Technol 2000, 26:706–714.View Article
- Entian KD, Kotter P: Yeast genetic strain and plasmid collections. Method Microbiol 2007, 36:629–666.View Article
- Verduyn 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.View Article
- Bahalul M, Kaneti G, Kashi Y: Ether-zymolyase ascospore isolation procedure: an efficient protocol for ascospores isolation in Saccharomyces cerevisiae yeast. Yeast 2010, 27:999–1003.View Article
- Gietz RD, Schiestl RH: High-efficiency yeast transformation using the LiAc/SS carrier DNA/PEG method. Nat Protoc 2007, 2:31–34.View Article
- Li H, Durbin R: Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 2010, 26:589–595.View Article
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25:1754–1760.View Article
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25:2078–2079.View Article
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res 1998, 8:175–185.
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res 1998, 8:186–194.
- Odorizzi G, Katzmann DJ, Babst M, Audhya A, Emr SD: Bro1 is an endosome-associated protein that functions in the MVB pathway in Saccharomyces cerevisiae . J Cell Sci 2003, 116:1893–1903.View Article
- Eilbeck K, Lewis SE, Mungall CJ, Yandell M, Stein L, Durbin R, Ashburner M: The Sequence Ontology: a tool for the unification of genome annotations. Genome Biol 2005, 6:R44.View Article
- Knijnenburg TA, de Winde JH, Daran JM, Daran-Lapujade P, Pronk JT, Reinders MJ, Wessels LF: Exploiting combinatorial cultivation conditions to infer transcriptional regulation. BMC Genomics 2007, 8:25.View Article
- Knijnenburg TA, Daran JM, van den Broek MA, Daran-Lapujade PA, de Winde JH, Pronk JT, Reinders MJ, Wessels LF: Combinatorial effects of environmental parameters on transcriptional regulation in Saccharomyces cerevisiae : a quantitative analysis of a compendium of chemostat-based transcriptome data. BMC Genomics 2009, 10:53.View Article
- Kresnowati MT, van Winden WA, Almering MJ, ten Pierick A, Ras C, Knijnenburg TA, Daran-Lapujade P, Pronk JT, Heijnen JJ, Daran JM: When transcriptome meets metabolome: fast cellular responses of yeast to sudden relief of glucose limitation. Mol Syst Biol 2006, 2:49.View Article
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.