A new laboratory evolution approach to select for constitutive acetic acid tolerance in Saccharomyces cerevisiae and identification of causal mutations
© The Author(s) 2016
Received: 1 April 2016
Accepted: 27 July 2016
Published: 12 August 2016
Acetic acid, released during hydrolysis of lignocellulosic feedstocks for second generation bioethanol production, inhibits yeast growth and alcoholic fermentation. Yeast biomass generated in a propagation step that precedes ethanol production should therefore express a high and constitutive level of acetic acid tolerance before introduction into lignocellulosic hydrolysates. However, earlier laboratory evolution strategies for increasing acetic acid tolerance of Saccharomyces cerevisiae, based on prolonged cultivation in the presence of acetic acid, selected for inducible rather than constitutive tolerance to this inhibitor.
Preadaptation in the presence of acetic acid was shown to strongly increase the fraction of yeast cells that could initiate growth in the presence of this inhibitor. Serial microaerobic batch cultivation, with alternating transfers to fresh medium with and without acetic acid, yielded evolved S. cerevisiae cultures with constitutive acetic acid tolerance. Single-cell lines isolated from five such evolution experiments after 50–55 transfers were selected for further study. An additional constitutively acetic acid tolerant mutant was selected after UV-mutagenesis. All six mutants showed an increased fraction of growing cells upon a transfer from a non-stressed condition to a medium containing acetic acid. Whole-genome sequencing identified six genes that contained (different) mutations in multiple acetic acid-tolerant mutants. Haploid segregation studies and expression of the mutant alleles in the unevolved ancestor strain identified causal mutations for the acquired acetic acid tolerance in four genes (ASG1, ADH3, SKS1 and GIS4). Effects of the mutations in ASG1, ADH3 and SKS1 on acetic acid tolerance were additive.
A novel laboratory evolution strategy based on alternating cultivation cycles in the presence and absence of acetic acid conferred a selective advantage to constitutively acetic acid-tolerant mutants and may be applicable for selection of constitutive tolerance to other stressors. Mutations in four genes (ASG1, ADH3, SKS1 and GIS4) were identified as causative for acetic acid tolerance. The laboratory evolution strategy as well as the identified mutations can contribute to improving acetic acid tolerance in industrial yeast strains.
KeywordsYeast Acetate Evolutionary engineering Bioethanol Inhibitors
Second generation bioethanol production uses lignocellulosic material from forestry residues, agricultural residues or energy crops as feedstocks. Use of these substrates is considered advantageous because of their abundance, availability, low cost and complementarity with food production [1–3]. Current processes for conversion of lignocellulosic biomass involve pretreatment of the substrate to disrupt its structure, followed by enzymatic hydrolysis to release monomeric sugars. The subsequent fermentation of these sugars to ethanol almost exclusively relies on the yeast Saccharomyces cerevisiae [4–7]. Hydrolysis not only releases fermentable sugars, but also furans, phenols and weak acids, whose presence can negatively affect yeast growth and ethanol production [8, 9].
Acetic acid, released during de-acetylation of hemicellulose, is the most abundant weak acid in lignocellulosic hydrolysates, at concentrations that can exceed 10 g/L . Inoculation into acetic acid containing media can cause a reduction in the specific growth rate and biomass yield of S. cerevisiae, as well as a substantially increased lag phase [11–13]. To minimize the risk of bacterial contamination and the requirement for base addition, alcoholic fermentation is preferably carried out at low pH, which enhances the toxicity of acetic acid [14–16]. At pH values below its dissociation constant (pKa = 4.76), acetic acid predominantly occurs in its protonated form, which can enter cells via passive diffusion [12, 16]. In the near-neutral cytosol, the acid dissociates into the acetate anion and a proton . To prevent intracellular acidification, yeast cells expel protons via the plasma membrane H+-ATPase Pma1 [11, 18–20]. The resulting diversion of ATP from cellular growth and maintenance is an important contributor to the inhibitory effects of acetic acid on S. cerevisiae . Additionally, intracellular accumulation of acetate anions can contribute to inhibition of specific cellular processes , osmotic stress, and in aerobic cultures, oxidative stress [22–24].
Upon transfer of S. cerevisiae cultures from a medium without acetic acid to a medium with an inhibitory concentration of acetic acid, only a small fraction of the population is able to resume growth . The fraction of cells able to grow decreases with increasing acetic acid concentration and is strain dependent [25, 26]. This culture heterogeneity contributes to a latency phase, defined as the time that elapses before growth is observed spectrophotometrically in an acetic acid containing medium. The fraction of a non-stressed population that is able to grow upon transfer to a medium with acetic acid is therefore a relevant measure for acetic acid tolerance of different strains . The actual length of the latency phase is further influenced by the time that this fraction of cells needs to start growing, and by their specific growth rate [20, 25].
To mitigate inhibitory effects of acetic acid in lignocellulosic hydrolysates, development of yeast strains with increased tolerance would be highly beneficial. Several previous studies aimed to understand and improve the response of S. cerevisiae to acetic acid. Transcriptome analysis identified a large number of genes that are differentially expressed upon exposure to acetic acid stress. These genes encode proteins involved in functions, such as transcription control, internal pH homeostasis, carbohydrate metabolism, cell-wall assembly, and biogenesis of mitochondria, ribosomes and the vacuole [27, 28]. However, despite clear progress, the underlying mechanisms of acetic acid tolerance remain incompletely understood , which hinders knowledge-based engineering strategies.
Natural and induced diversity of acetic acid tolerance among S. cerevisiae strains has also been explored to investigate its genetic basis. For example, generation of haploid segregants from crosses of highly tolerant and less tolerant strains, followed by quantitative-trait-locus (QTL) analysis, enabled the identification of multiple alleles that affect acetic acid tolerance in this yeast . Increased acetic acid tolerance has also been reported for strains carrying targeted genetic modifications, such as overexpression of TAL1, encoding a transaldolase , overexpression of PEP3, encoding a protein subunit involved in vacuolar biogenesis , overexpression of the transcription factor HAA1 , introduction of an artificial zinc-finger based transcription factor  and introduction of an ascorbic-acid production pathway . However, the levels of tolerance reached by these approaches are not sufficient for the high concentrations of acetic acid present in lignocellulosic hydrolysates at industrially relevant pH values.
Laboratory evolution allows the selection of specific phenotypes without the requirement for in-depth understanding of the underlying mechanisms [9, 36, 37] and has been successfully used to improve the tolerance of yeast to butanol  and other stresses, such as freezing–thawing cycles, increased temperature, ethanol, and oxidative stress . This approach has also been used to improve acetic acid tolerance in bacteria  and yeast. Wright et al.  used repeated batch cultivation at progressively increasing concentrations of acetic acid to select a S. cerevisiae strain with a higher acid–acid tolerance. However, after growing the evolved strain in the absence of acetic acid, the strain was no longer able to grow at high acetic acid concentrations. Although the strain was genetically stable, increased tolerance to acetic acid could only be shown after pre-cultivation in the presence of acetic acid. These observations indicated that the increased tolerance of the evolved strain was not constitutive but required induction by acetic acid . In another study, induction of acetic acid tolerance was shown to occur upon exposure of cells to low levels of acetic acid, resulting in a reversible acetic acid tolerant phenotype . In accordance with the proposed role of cytosolic pH in acetic acid tolerance, these adapted cells showed a less dramatic decrease in intracellular pH upon exposure to acetic acid compared to non-adapted cells upon exposure to acetic acid [11, 26]. Preadaptation to acetic acid was also shown to protect yeast cells from acetic acid-mediated programmed cell death . To benefit from such an inducible acetic acid tolerance in industrial processes, an adaptation step would have to be implemented prior to the fermentation, making the whole process more complex. Yeast strains with constitutively high-level tolerances, which do not require adaptation, would therefore clearly be preferable for industrial applications.
The aim of this study is to investigate how laboratory evolution and mutagenesis can be used to obtain S. cerevisiae strains with constitutive acetic acid tolerance. To maximize the selective pressure on constitutively tolerant strains and simultaneously minimize the competitiveness of cells with inducible acetic acid tolerance, laboratory evolution was performed by alternatingly growing S. cerevisiae cells in medium with and without acetic acid. In parallel, mutagenesis and screening was applied to obtain constitutively tolerant mutants without laboratory evolution. Subsequently, a combination of whole-genome sequencing and crossing and segregation studies was used to investigate the mutations underlying the constitutively acetic acid tolerant phenotype. The impact of these mutations on acetic acid tolerance was confirmed by reverse engineering the evolved phenotype into the unevolved reference strain.
Strains, media and growth conditions
Acetic acid tolerant mutants derived from CEN.PK113-7D through evolution/mutagenesis and crossing
Acetic acid tolerant mutant evolved from CEN.PK113-7D
Acetic acid tolerant mutant evolved from CEN.PK113-7D
Acetic acid tolerant mutant evolved from CEN.PK113-7D
Acetic acid tolerant mutant evolved from CEN.PK113-7D
Acetic acid tolerant mutant evolved from CEN.PK113-7D
Acetic acid tolerant UV-radiated mutant from CEN.PK113-7D
Diploid obtained from crossing MUT1A with IMK439
Diploid obtained from crossing MUT2B with IMK439
Diploid obtained from crossing MUT3E with IMK439
Diploid obtained from crossing HAT1E with IMK439
Diploid obtained from crossing HAT2A with IMK439
Engineered S. cerevisiae strains constructed in this study by replacement, in the reference strain CEN.PK113-7D, of the native alleles of ASG1, ADH3, GIS4 and/or SKS1 by mutated alleles identified in mutagenized and/or evolved acetic acid tolerant strains
Origin of mutation
MATa, asg1G1248A, adh3G416T
MATa, asg1G1248T, adh3T966G
MATa, asg1G1248T, adh3T966G, sks1G821T
MATa, asg1G2881C, adh3T201A
MATa, asg1G2881C, adh3T201A, sks1C617A
Repeated batch cultivation of S. cerevisiae CEN.PK113-7D was performed in a simple microaerobic cultivation system. Cultures were grown in 30 mL serum bottles, containing 25 mL medium and closed with a butyl-rubber seal. The seal was punctured with a needle connected to a syringe to enable release of overpressure generated by CO2 production. These cultures were incubated at 30 °C and shaken at 200 rpm. Serial transfer was performed by alternating cultivation in SMG without acetic acid with cultivation in SMG containing increasing concentrations of acetic acid, ranging from 10 to 18 g/L. Cultures with acetic acid were inoculated by direct transfer of 1 mL from the prior bottle without acetic acid. Cultures without acetic acid were inoculated by spinning down 0.1 mL of broth from the prior culture with acetic acid, which was washed with water and resuspended in SMG prior to inoculation. Five independent evolution lines were performed (MUT1, MUT2, MUT3, HAT1 and HAT2). The HAT1 and HAT2 evolution lines were started from CEN.PK113-7D stock cultures while the MUT1, MUT2 and MUT3 evolution lines were started with UV-mutagenized CEN.PK113-7D. Before UV mutagenesis, exponentially growing shake-flask cultures on YPD were washed with water and resuspended in SMG to an OD660 of 1.0. Subsequently, 25 mL of the resulting cell suspension was poured in a (9 cm diameter) petri dish, shaken at 100 rpm inside a laminar-flow cabinet and irradiated with a UV lamp (TUV 30 W T8, Philips, Eindhoven, The Netherlands) at a radiation peak of 253.7 nm. Suspensions used to initiate evolution lines MUT1, MUT2 and MUT3 were irradiated with UV doses yielding survival rates of 18, 22 and 16 %, respectively. Evolution experiments MUT1 and MUT2 involved 50 transfers, while the MUT3, HAT1 and HAT2 experiments involved 55, 54 and 51 transfers, respectively. At the end of the laboratory evolution experiments, single cell lines were isolated from each evolution line by streaking on SMG agar plates.
Mutagenesis and screening for acetic acid tolerant mutants
Saccharomyces cerevisiae CEN.PK113-7D was irradiated with UV light as described above, but at a dose resulting in a survival rate of 2 %. After irradiation, the cell suspension was transferred to a serum flask and grown in SMG in the dark for 2 days. 40 colonies were isolated on SMG agar plates with 9 g/L acetic acid and each transferred into one well of a 96-well microtiter plate containing 50 μL of SMG using a sterile pipet tip. Additionally, 8 wells were inoculated with the reference strain CEN.PK113-7D. The microtiter plate was then sealed with a gas impermeable seal (NUNC, Roskilde, Denmark) and incubated for 2 days in an orbital shaker (PHMP-4 Thermoshaker for microplates, Grant Instruments, Shepreth, UK) at 30 °C and 700 rpm. This initial plate was used to inoculate four 96-well plates containing 100 µL/well of SMG, two with 9 g/L acetic acid and two with 10 g/L acetic acid. These 96-well plates were sealed and incubated (30 °C, 700 rpm). After four days, the OD660 of all wells was measured and the five mutants with the highest OD660 values at each acetic acid concentrations were selected. Single-cell lines of mutants were isolated by picking a colony from YPD agar plates and growing them overnight in microaerobic shake-flasks containing 20 mL SMG, along with a similar culture of the reference strain CEN.PK113-7D. All cultures were then diluted with sterile water to an OD660 of 5.0 and eight 50 µL aliquots of each strain were distributed into a 96-well plate. This microtiter plate was used to inoculate a microtiter plate containing 100 µL/well of SMG with 10 g/L acetic acid using a sterile pin replicator. The plate was sealed, incubated for 4 days and the final OD660 was measured. The mutant with the highest average OD660 was named UV-E3 and further characterized for its acetic acid tolerance.
Analysis of acetic acid tolerance: OD660 in stationary phase
Shake-flask cultures were grown overnight until stationary phase in SMG and diluted with sterile water to an OD660 of 5. 96-well plates were filled with 100 μL/well of SMG containing acetic acid concentrations ranging from 0 to 15 g/L (12 different concentrations with 8 replicas each). The plates were inoculated with a sterile pin replicator, sealed with a gas impermeable seal (NUNC, Roskilde, Denmark) and incubated at 30 °C for 5 days without shaking. After incubation, cells were resuspended in a MS2 minishaker (IKA, Staufen, Germany), seals were removed and OD660 of each well was measured in a GENIos Pro micro plate spectrophotometer (Tecan, Männedorf, Switzerland). Average OD660 and standard deviation of replicate wells were then calculated. Where indicated in the results section, glucose concentrations in culture supernatants were obtained after centrifugation and 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.
Analysis of acetic acid tolerance: latency phase and growth rates
Overnight shake-flask cultures grown until stationary phase in SMG were used to inoculate 1.5-mL Eppendorf tubes containing 1 mL of SMG with acetic acid concentrations ranging from 0 to 15 g/L, at an initial OD660 of 0.2 for growth rate determination or of 0.1 for latency phase determination. 100 μL aliquots from each tube were distributed into different wells of a 96-well plate. The plates were sealed with a gas impermeable seal and incubated in a GENIos Pro micro plate spectrophotometer at 400 rpm for 5 days. OD660 was recorded at 15 min intervals. The duration of the latency phase was determined as the time after inoculation at which the OD660 increased above 0.12. Specific growth rates were estimated from exponential fits of the OD660 versus time.
Analysis of acetic acid tolerance: fraction of cells growing in the presence of acetic acid
Overnight shake-flask cultures grown until stationary phase in SMG were diluted with sterile water to an OD660 of 1. 100 μL of serial tenfold dilution of these cell suspensions were plated on SMG agar plates (pH 4.5) with and without 9 g/L acetic acid. The number of colony-forming units (CFU) per mL for each dilution was calculated after 5 days incubation in a Bactron X-2E anaerobic chamber (Sheldon Manufacturing Inc., Cornelius, OR) at 30 °C. For each strain, the fraction of cells able to grow on agar plates containing acetic acid was calculated by dividing the CFU/mL counts on agar plates with acetic acid by the CFU/mL count without acetic acid.
Whole-genome sequencing of the acetic acid tolerant mutants
Genomic DNA of the strains CEN.PK113-7D, HAT1E, HAT2A, MUT1A, MUT2B, MUT3E and UV-E3 was prepared as described previously . Libraries of 300-bp inserts were constructed and paired-end sequenced (101 base pair reads) using an Illumina HISeq 2500 sequencer (Baseclear BV, Leiden, The Netherlands). A minimum data quantity of 1 GB was generated for each strain, representing a minimum 83-fold average coverage. Sequence reads of each strain were mapped onto the CEN.PK113-7D genome  using the Burrows–Wheeler Alignment tool (BWA) and further processed using SAMtools [52–54]. Single-nucleotide variations were determined 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 variant quality, defined as the Phred-scaled probability that the mutation call is incorrect [55, 56], of at least 20 were kept and (iii) mutations with a depth of coverage <10X were discarded. Finally, the single-nucleotide variations were positioned by their genomic locations; coding effects were predicted and functionally annotated according to the CEN.PK113-7D sequence annotation . Raw sequencing data have been deposited as short-read archives (Bioproject: PRJNA313456/SRP070976, individual accession numbers in Additional file 1). The Magnolya algorithm was used to analyse copy number variation using Newbler (454 Life Sciences, Branford, CT) for co-assembly of sequence reads of CEN.PK113-7D and the sample of interest . Integer copy numbers of the assembled contigs were estimated using the Poisson mixture model (PMM) algorithm in Magnolya.
Crossing, sporulation and screening for acetic acid-tolerant haploid segregants
Mutants evolved for increased acetic acid tolerance were crossed with IMK439, a strain derived from CEN.PK113-1A (MATα, otherwise isogenic to CEN.PK113-7D) by replacing its URA3 gene with the kanMX marker . To this end, each mutant was grown together with IMK439 on SMG agar plates without uracil and with 200 µg/mL G418, on which only diploidized strains could grow. After sporulation, random spore isolation was performed as described previously  and spores were plated on SMG without uracil. Each segregant was then transferred into one well of a 96-well plate containing 50 µL of SMG using a sterile pipet tip. On every plate, the reference strain CEN.PK113-7D, the evolved acetic acid tolerant parent and the diploid strain from which the segregants were derived, were each inoculated in 8 wells as internal controls. After sealing the plate and overnight incubation at 30 °C without shaking, these plates were used to inoculate 96-well plates filled with 100 μL/well of SMG with acetic acid concentrations ranging from 9 to 15 g/L using a sterile pin replicator. After sealing with a gas impermeable seal, plates were incubated at 30 °C for 5 days without shaking, followed by measurement of OD660. For each mutant, ten segregants with the same tolerance as the evolved tolerant mutant and ten segregants with the same sensitivity as CEN.PK113-7D were selected for further analysis. These strains were grown overnight in microaerobic shake-flasks SMG, were diluted with water to an OD660 of 5 and eight 100 µL aliquots of each strain were distributed into a 96-well plate. From this plate, 96-well plates containing 100 µL SMG/well, containing acetic acid concentrations ranging from 8 to 15 g/L, were inoculated. Plates were sealed with a gas impermeable seal, incubated at 30 °C for 5 days without shaking, followed by measurement of OD660. Segregants for which the phenotype observed in the initial test was confirmed were selected for further study.
Sequencing of selected genes in haploid segregants
Different haploid segregants derived from the cross of the evolved strains with IMK439 were grown separately in SMG and the volume of aliquot equivalent to 1 mL with an OD660 of 1.0 was determined for each culture. For each evolution mutant, the aliquots of all (4–9) acetic acid-tolerant segregants were pooled (with a maximum of 6 segregants per pool), centrifuged, and after removal of supernatant, resuspended in water. This procedure was also applied to the (1 or 2) sensitive segregants of each evolved strain. Genomic DNA was extracted from the resulting cell suspensions with the YeaStar™ Genomic DNA kit (Zymo Research, Irvine, CA) and used as a template for amplification of the mutated genes identified in each evolved strain (primers listed in Additional file 2). Amplified alleles were Sanger sequenced (Baseclear) and presence of the wild-type, mutated allele or both were assessed with Clone Manager (Scientific & Educational Software, Morrisville, NC).
Reverse engineering of acetic acid tolerance
Inducibility of acetic acid tolerance in S. cerevisiae
To investigate induction of acetic acid tolerance by permissive concentrations of acetic acid, shake-flask cultures of S. cerevisiae CEN.PK113-7D cells were pre-grown on SMG medium pH 4.5 with or without 9 g/L acetic acid. These non-adapted and adapted cultures were used to inoculate microtiter plate cultures on SMG, supplemented with acetic acid at concentrations ranging from 0 to 15 g/L. When growth occurred, glucose was always fully consumed within 5 days. Final OD660, measured after 5 days, was therefore used as an indicator for biomass yield on glucose.
To investigate how prior adaptation to acetic acid affected the fraction of the population capable of initiating growth in acetic acid containing media, CEN.PK113-7D was grown until stationary phase in five sequential microaerobic shake-flask cultures: first in SMG medium without acetic acid, then in SMG medium with acetic acid at pH 4.5 to induce acetic acid tolerance, and then in three sequential cultures in SMG medium without acetic acid. At the end of each culture, tenfold dilutions of the cell broth were plated on agar with and without acetic acid and colony-forming units (CFU) were counted after 4–6 days of incubation to determine the fraction of growing cells. In the first culture, grown on SMG without acetic acid (pH 6.0), only 50 out of a million cells (a fraction of 0.00005) grew upon transfer to SMG agar with 9 g/L acetic acid (referred to as #1 AA−, Fig. 2d). In the stationary phase of the subsequent culture, grown on SMG and supplemented with 9 g/L acetic acid (pH 4.5), the fraction of growing cells had increased by more than 10,000-fold to 0.88 (referred to as #2 AA+ , Fig. 2d). Already in the stationary phase of the first of three subsequent shake-flasks on SMG without acetic acid (pH 6.0), this fraction had decreased to 0.00025 and did not change significantly during the two subsequent cultures in the same medium (#3 AA−, #4 AA−, #5 AA−, Fig. 2d). These results show that adaptation of S. cerevisiae to acetic acid induces an apparent tolerance of the overall culture, primarily by increasing the fraction of the population that is able to grow in acetic acid containing medium. However, after a single culture under non-stressed conditions, this fraction rapidly returns to within one order of magnitude of its original, non-adapted level.
Laboratory evolution for constitutive acetic acid tolerance
To investigate whether the tolerance observed in the laboratory evolution experiments was indeed constitutive, three additional transfers of all five evolved cultures were performed in SMG without acetic acid (pH 6.0). After the last culture, the fraction of cells able to grow on SMG agar plates with 9 g/L acetic acid was determined (Fig. 3b). As a reference, a culture of CEN.PK113-7D was pre-grown with 9 g/L acetic acid to induce acetic acid tolerance and subsequently transferred to new medium containing 18 g/L acetic acid prior to the determination of the fraction of growing cells after three additional cultures without acetic acid (Fig. 3b). After these three transfers, all five independently evolved cultures showed at least a 1000-fold higher fraction of the population that was able to grown on SMG agar with 9 g/L acetic acid (pH 4.5) compared to the unevolved reference strain CEN.PK113-7D (Fig. 3b).
Selection of a constitutively acetic acid tolerant mutant by mutagenesis and screening
In parallel to the laboratory evolution approach, a constitutively acetic acid tolerant mutant was isolated by a modified procedure for random mutagenesis and screening. Preliminary experiments based on UV mutagenesis followed by testing of random colonies for acetic acid tolerance failed to yield constitutively tolerant mutants. Indeed, we were unable to isolate mutants (based on single colonies on plates) with increased acetic acid tolerance directly after UV mutagenesis (data not shown). We assume this is due to the extremely low probability that a single cell has to resume proliferation in the presence of acetic acid, even if it contains a mutation that is beneficial for acetic acid tolerance. Therefore, after UV-irradiation, cells were first grown in SMG (pH 6.0) for 2 days without acetic acid before plating on SMG with 9 g/L acetic acid (pH 4.5). In total 40 randomly selected single-cell-colony isolates from this experiment were then pre-grown in SMG without acetic acid (pH 6.0), followed by screening of their tolerance in 96-well plates containing SMG with 9 and 10 g/L acetic acid (pH 4.5). The five UV-mutants that reached the highest OD660 in this screening were further characterized in 96-well plates (8 replicates per strain) containing 7, 8, 9 or 10 g/L acetic acid. In these experiments, strain UV-E3 consistently showed the highest final OD660. Its increased constitutive acetic acid tolerance relative to CEN.PK113-7D was confirmed by growing it in 96-well plates on SMG with acetic acid concentrations ranging between 0 and 15 g/L (Fig. 4). By allowing mutants to proliferate by prior cultivation in medium without acetic acid, the absolute number of cells harbouring specific mutations was increased. As even for beneficial mutants, only a fraction of cells is able to grow upon transfer to medium with acetic acid, the generally increased number of cells decreased the chance to lose a beneficial genotype due to single-cell stochasticity of growth. This modified procedure resulted in the selection of strain UV-E3 with improved constitutive acetic acid tolerance.
Genomic mutations in constitutively acetic acid tolerant mutants
Mutations present in multiple evolved S. cerevisiae strains and genetic analysis of their relevance for acetic acid tolerance
Nucleotide in THS
Nucleotide in SHS
Genetic analysis for identification of causal mutations
The observation that ASG1, ADH3, SKS1, GIS4, SAC6 and EUG1 were mutated in multiple acetic acid tolerant strains suggested that they might contribute to this acquired phenotype. To further investigate a possible role of these genes in tolerance, a genetic analysis was performed. The strains obtained from the laboratory evolution experiments (but not the UV-mutagenized strain UV-E3), were mated with IMK439, a MATα ura3Δ::KanMX strain congenic with CEN.PK113-7D. The resulting diploids were sporulated after which haploid segregants were tested for acetic acid tolerance. Subsequently, the presence or absence of mutated alleles was tested in haploid segregants with the same tolerance level as the mutants (tolerant segregants) and as CEN.PK113-7D (sensitive segregants). To this end, genomic DNA of all tolerant segregants derived from the same diploid was pooled (Additional file 9) and fragments containing the relevant SNPs in ASG1, ADH3, SKS1, GIS4, SAC6 and EUG1 were PCR amplified. The same was done for a pool of sensitive haploid segregants derived from each diploid. Sanger sequencing revealed that mutated and wild-type alleles of SAC6 and EUG1 were randomly distributed in the population of tolerant and sensitive haploid segregants from all five evolved strains (Table 3). Conversely, mutated alleles of ASG1 and GIS4 were exclusively found in tolerant segregants, while all sensitive segregants tested contained the wild-type allele (Table 3). It was therefore concluded that mutations in ASG1 and in GIS4, but not those in SAC6 and EUG1, contributed to the acquired acetic acid tolerance. For both ADH3 and SKS1, all the sensitive haploid segregants contained the wild-type alleles, but the tolerant haploid segregants did not exclusively contain mutated alleles. The pooled tolerant segregants from strain MUT2B contained both wild-type and mutated alleles of ADH3 and only the mutated allele of SKS1. Oppositely, for HAT2A both wild-type and mutated alleles of SKS1 and only the mutated allele of ADH3 were identified in the pooled tolerant segregants. This suggests that mutations in ADH3 and SKS1 contributed to acetic acid tolerance, but their impact may not be additive.
Reverse engineering of constitutive acetic acid tolerance
Acetic acid is a key inhibitor of yeast performance in lignocellulosic hydrolysates, but can also be produced by S. cerevisiae during cultivation under anaerobic conditions . Wild-type strains of S. cerevisiae have evolved an inducible tolerance against this organic acid, in which the HAA1 regulon, which includes the TPO2- and TPO3-encoded membrane transporters, plays an important role . In a previous study, in which S. cerevisiae was evolved in the laboratory through prolonged growth in the presence of progressively increasing acetic acid concentrations, the evolved cultures acquired an increased tolerance which, however, required prior induction by acetic acid . At first glance, such a ‘hyper-inducible’ acetic acid tolerance phenotype can be interpreted as reflecting an increased maximum expression level of the natural, inducible acetic acid tolerance mechanisms by all cells in the culture . Recent observations, however, suggest that the situation may be more complicated. Pronounced heterogeneity in the adaptation of genetically homogeneous microbial populations to changing environmental conditions has been described for many micro-organisms and culture parameters [68–70]. Recently, such a heterogeneous response was also observed upon the exposure of non-stressed S. cerevisiae cultures to acetic acid stress, which showed that only a small fraction of the population resumed growth [20, 25, 71]. For a given S. cerevisiae strain, the fraction of a non-stressed and non-induced population able to grow at a specific acetic acid concentration and pH was shown to be constant and an important determinant of acetic acid tolerance .
In this study, we show that preadaptation in the presence of acetic acid strongly increases the fraction of the population that is able to grow in acetic acid containing medium. In particular, preadaptation strongly reduced the duration of the latency phase upon transfer to acetic acid containing media, while growth rate and biomass yield remain unaltered. Preadaptation and culture heterogeneity are not explicitly included in other recent studies on selection and/or construction of acetic acid tolerant S. cerevisiae strains [30–35]. However, in those studies, pre-cultivation typically involved batch cultivation at low pH in glucose-containing media. It is conceivable that in these cultures, acetic acid produced by the yeast cultures was sufficient to induce acetic acid tolerance, an effect known as hormesis . In contrast, in the current industrial practice of second-generation bioethanol production, yeast biomass is typically pre-grown in a separate aerobic, sugar-limited fed-batch process. Conditions in this biomass propagation phase are especially chosen to optimize biomass yield, and therefore, to minimize formation of by-products such as ethanol and acetic acid. To ensure an efficient start of anaerobic fermentation of the subsequent acetic acid-containing lignocellulosic hydrolysates, inhibitor tolerance should already be fully expressed during the preceding ‘non-stressed’ biomass propagation phase.
The observation that acetic acid adaptation of S. cerevisiae cultures was rapidly lost upon growth in the absence of acetic acid at pH 6.0, provided the basis for a simple evolutionary engineering strategy. Prolonged cultivation in five independent evolution experiments with alternating stressed and non-stressed culture cycles reproducibly yielded evolved S. cerevisiae cultures that, after pre-cultivation in the absence of acetic acid, could immediately initiate growth at high acetic acid concentrations. As a side effect, this strategy adds a selective pressure against mutations that in addition to benefits for acetic acid tolerance results in decreased fitness during the off-phase, which would potentially decrease the industrial applicability. This simple ‘on–off’ strategy is likely to be applicable for the selection of constitutive tolerance to a variety of chemical or physical stresses, especially in situations where micro-organisms already harbour native, inducible tolerance mechanisms.
After initial identification of mutations by whole-genome resequencing of independent evolution lines and classical genetic analysis, causality of SNPs in ASG1, ADH3, SKS1 and GIS4 for acetic acid tolerance was confirmed by reverse engineering of an unevolved strain. Interestingly, three of these genes (ADH3, SKS1 and GIS4) were not previously associated with acetic acid tolerance in genome-wide expression studies [27, 28], screening of deletion libraries  or in a recent quantitative-trait-loci analysis . This difference can probably be explained from the focus, in this study, on constitutive tolerance rather than on the maximum permissive acetic acid concentration. Indeed, the acquired acetic acid tolerance of the evolved cultures was predominantly due to a strong increase of the fraction of the population that was able to initiate growth upon a transfer to acetic acid containing conditions.
Based on sequence homology and limited functional analysis, Asg1, Sks1 and Gis4 have been implicated in regulation, but their precise functions in acetic acid tolerance remain to be elucidated. Protein sequence information available from the Saccharomyces genome database  indicates that the mutations identified in this study did not affect amino-acid residues that are known or predicted to be post-translationally modified or to be in narrowly defined functional domains. However, deletion of ASG1 in S. cerevisiae was previously shown to result in increased susceptibility to acetic acid. Indeed, in Candida glabrata deletion of ASG1 resulted in reduced plasma membrane H+-ATPase activity, an inability to regulate intracellular pH and leading to growth inhibition in acidic environments. Identification of causal mutations in ADH3, which encodes the main mitochondrial isoenzyme of NAD+-dependent alcohol dehydrogenase, is puzzling, since it is difficult to imagine how its catalytic role could influence acetic acid tolerance . A recently reported role of Adh3 in cold tolerance of Saccharomyces yeasts  was tentatively attributed to its role in cellular redox balancing. Clearly, further molecular studies are required to elucidate how the four genes identified in this study affect the distribution of acetic acid-tolerant and -sensitive cells in cultures that have not previously been exposed to acetic acid. Recently, a link between the cytosolic pH and stochasticity of growth in the presence of acetic acid was identified . The fact that stochasticity of growth is determined at the single cell level precludes elucidation of the underlying mechanisms at a population level. In this light, recent developments on genomics of single cells in combination with microfluidic devices might prove to be valuable to investigate this cell-to-cell heterogeneity as well as the mechanisms of identified beneficial mutations [76–78]. Even without full characterization of their mechanistic roles in acetic acid tolerance, the genes and alleles identified in this study provide leads for the improvement of industrial strains of S. cerevisiae.
Alternating growth in the presence and absence of acetic acid provided a selective pressure for mutants with constitutive acetic acid tolerance, as reflected by an increased fraction of the population that was able to initiate growth in the presence of high concentrations of acetic acid without prior adaptation to this stressor. Awareness of the importance of culture heterogeneity and preadaptation proved valuable in the design of optimized laboratory evolution as well as mutagenesis strategies. Similar approaches are likely to be applicable to improve other characteristics of industrial micro-organisms that are inducible and/or involve culture heterogeneity. Mutations in three genes (ASG1, ADH3 and SKS1) were identified that individually, as well as synergistically, contributed to acetic acid tolerance. Additionally, alleles of GIS4 were identified that caused a similar level of constitutive acetic acid tolerance. Both strain improvement approaches, as well as the underlying mutations, can be used for industrial strain engineering and thereby stimulate further development of second generation processes for the production of fuels and chemicals.
DGR, SS, EN, JMGD, JTP and AJAvM conceived and designed the experiments. Laboratory evolution was performed by DGR. Mutagenesis and screening was done by MCvB and DGR. Quantification of constitutive and induced tolerance was done by ARGDV and DGR. Sequence assembly of the genome sequences was performed by MvdB. Crossing and sporulation of the mutants, as well as introduction of the evolved alleles into the parental strain, was done by SSG and DGR. The manuscript was written by DGR, ARGDV, JMGD, JTP and AJAvM. All authors read and approved the final version of the manuscript.
We thank Marijke A. H. Luttik and Erik de Hulster for their assistance throughout the project. We thank Prof. Rinji Akada for kindly providing plasmid pGG119.
The authors declare that they have no competing interests.
Availability of supporting data
Raw sequencing data have been deposited as short-read archives (Bioproject: PRJNA313456/SRP070976, individual accession numbers in additional file 1). Other datasets supporting the conclusions of this article are included within the article and its additional files.
This work was funded through the ERA-NET Scheme of the 6th EU Framework Program (INTACT).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Sarkar N, Ghosh SK, Bannerjee S, Aikat K. Bioethanol production from agricultural wastes: an overview. Renew Energ. 2012;37:19–27.View ArticleGoogle Scholar
- Demirbaş A. Bioethanol from cellulosic materials: a renewable motor fuel from biomass. Energ Sour. 2005;27:327–37.View ArticleGoogle Scholar
- Antizar-Ladislao B, Turrion-Gomez JL. Second-generation biofuels and local bioenergy systems. Biofuel Bioprod Biorefin. 2008;2:455–69.View ArticleGoogle Scholar
- Taherzadeh MJ, Karimi K. Pretreatment of lignocellulosic wastes to improve ethanol and biogas production: a review. Int J Mol Sci. 2008;9:1621–51.View ArticleGoogle Scholar
- Sun Y, Cheng J. Hydrolysis of lignocellulosic materials for ethanol production: a review. Bioresour Technol. 2002;83:1–11.View ArticleGoogle Scholar
- Laluce C, Schenberg AC, Gallardo JC, Coradello LF, Pombeiro-Sponchiado SR. Advances and developments in strategies to improve strains of Saccharomyces cerevisiae and processes to obtain the lignocellulosic ethanol–a review. Appl Biochem Biotechnol. 2012;166:1908–26.View ArticleGoogle Scholar
- Kumar P, Barrett DM, Delwiche MJ, Stroeve P. Methods for pretreatment of lignocellulosic biomass for efficient hydrolysis and biofuel production. Ind Eng Chem Res. 2009;48:3713–29.View ArticleGoogle Scholar
- Quintas C, Leyva JS, Sotoca R, Loureiro-Dias MC, Peinado JM. A model of the specific growth rate inhibition by weak acids in yeasts based on energy requirements. Int J Food Microbiol. 2005;100:125–30.View ArticleGoogle Scholar
- Almeida JRM, Modig T, Petersson A, Hähn-Hägerdal B, Lidén G, Gorwa-Grauslund MF. Increased tolerance and conversion of inhibitors in lignocellulosic hydrolysates by Saccharomyces cerevisiae. J Chem Technol Biot. 2007;82:340–9.View ArticleGoogle Scholar
- Taherzadeh MJ, Eklund R, Gustafsson L, Niklasson C, Lidén G. Characterization and fermentation of dilute-acid hydrolyzates from wood. Ind Eng Chem Res. 1997;36:4659–65.View ArticleGoogle Scholar
- Pampulha ME, Loureiro-Dias MC. Energetics of the effect of acetic acid on growth of Saccharomyces cerevisiae. FEMS Microbiol Lett. 2000;184:69–72.View ArticleGoogle Scholar
- Stratford M, Anslow PA. Comparison of the inhibitory action on Saccharomyces cerevisiae of weak-acid preservatives, uncouplers, and medium-chain fatty acids. FEMS Microbiol Lett. 1996;142:53–8.View ArticleGoogle Scholar
- Narendranath NV, Thomas KC, Ingledew WM. Effects of acetic acid and lactic acid on the growth of Saccharomyces cerevisiae in a minimal medium. J Ind Microbiol Biotechnol. 2001;26:171–7.View ArticleGoogle Scholar
- Thomas KC, Hynes SH, Ingledew WM. Influence of medium buffering capacity on inhibition of Saccharomyces cerevisiae growth by acetic and lactic acids. Appl Eniron Microbiol. 2002;68:1616–23.View ArticleGoogle Scholar
- Graves T, Narendranath NV, Dawson K, Power R. Effect of pH and lactic or acetic acid on ethanol productivity by Saccharomyces cerevisiae in corn mash. J Ind Microbiol Biotechnol. 2006;33:469–74.View ArticleGoogle Scholar
- Casal M, Cardoso H, Leão C. Mechanisms regulating the transport of acetic acid in Saccharomyces cerevisiae. Microbiology. 1996;142:1385–90.View ArticleGoogle Scholar
- Baronofsky JJ, Schreurs WJA, Kashket ER. Uncoupling by acetic acid limits growth of and acetogenesis by Clostridium thermoaceticum. Appl Environ Microbiol. 1984;48:1134–9.Google Scholar
- Carmelo V, Santos H, Sá-Correia I. Effect of extracellular acidification on the activity of plasma membrane ATPase and on the cytosolic and vacuolar pH of Saccharomyces cerevisiae. Biochim Biophys Acta. 1997;1325:63–70.View ArticleGoogle Scholar
- Holyoak CD, Stratford M, McMullin Z, Cole MB, Crimmins K, Brown AJ, Coote PJ. Activity of the plasma membrane H+-ATPase and optimal glycolytic flux are required for rapid adaptation and growth of Saccharomyces cerevisiae in the presence of the weak-acid preservative sorbic acid. Appl Environ Microbiol. 1996;62:3158–64.Google Scholar
- Viegas CA, Almeida PF, Cavaco M, Sá-Correia I. The H+-ATPase in the plasma membrane of Saccharomyces cerevisiae is activated during growth latency in octanoic acid-supplemented medium accompanying the decrease in intracellular pH and cell viability. Appl Environ Microbiol. 1998;64:779–83.Google Scholar
- Russell JB. Another explanation for the toxicity of fermentation acids at low pH—anion accumulation versus uncoupling. J Appl Bacteriol. 1992;73:363–70.View ArticleGoogle Scholar
- Piper P, Calderon CO, Hatzixanthis K, Mollapour M. Weak acid adaptation: the stress response that confers yeasts with resistance to organic acid food preservatives. Microbiology. 2001;147:2635–42.View ArticleGoogle Scholar
- Pampulha ME, Loureiro-Dias MC. Activity of glycolytic-enzymes of Saccharomyces cerevisiae in the presence of acetic acid. Appl Microbiol Biotechnol. 1990;34:375–80.View ArticleGoogle Scholar
- Giannattasio S, Guaragnella N, Zdralević M, Marra E. Molecular mechanisms of Saccharomyces cerevisiae stress adaptation and programmed cell death in response to acetic acid. Front Microbiol. 2013;4:33.View ArticleGoogle Scholar
- Swinnen S, Fernández-Niño M, González-Ramos D, van Maris AJA, Nevoigt E. The fraction of cells that resume growth after acetic acid addition is a strain-dependent parameter of acetic acid tolerance in Saccharomyces cerevisiae. FEMS Yeast Res. 2014;14:642–53.View ArticleGoogle Scholar
- Fernández-Niño M, Marquina M, Swinnen S, Rodríguez-Porrata B, Nevoigt E, Ariño J. The cytosolic pH of individual Saccharomyces cerevisiae cells is a key factor in acetic acid tolerance. Appl Environ Microbiol. 2015;81:7813–21.View ArticleGoogle Scholar
- Mira NP, Palma M, Guerreiro JF, Sá-Correia I. Genome-wide identification of Saccharomyces cerevisiae genes required for tolerance to acetic acid. Microb Cell Fact. 2010;9:79.View ArticleGoogle Scholar
- Li B-Z, Yuan Y-J. Transcriptome shifts in response to furfural and acetic acid in Saccharomyces cerevisiae. Appl Microbiol Biotechnol. 2010;86:1915–24.View ArticleGoogle Scholar
- Guo Z, Olsson L. Physiological response of Saccharomyces cerevisiae to weak acids present in lignocellulosic hydrolysate. FEMS Yeast Res. 2014;14:1234–48.View ArticleGoogle Scholar
- Meijnen JP, Randazzo P, Foulquié-Moreno MR, van den Brink J, Vandecruys P, Stojiljkovic M, Dumortier F, Zalar P, Boekhout T, Gunde-Cimerman N, et al. Polygenic analysis and targeted improvement of the complex trait of high acetic acid tolerance in the yeast Saccharomyces cerevisiae. Biotechnol Biofuels. 2016;9:5.View ArticleGoogle Scholar
- Hasunuma T, Sanda T, Yamada R, Yoshimura K, Ishii J, Kondo A. Metabolic pathway engineering based on metabolomics confers acetic and formic acid tolerance to a recombinant xylose-fermenting strain of Saccharomyces cerevisiae. Microb Cell Fact. 2011;10:1–13.View ArticleGoogle Scholar
- Ding J, Holzwarth G, Bradford CS, Cooley B, Yoshinaga AS, Patton-Vogt J, Abeliovich H, Penner MH, Bakalinsky AT. PEP3 overexpression shortens lag phase but does not alter growth rate in Saccharomyces cerevisiae exposed to acetic acid stress. Appl Microbiol Biotechnol. 2015;99:8667–80.View ArticleGoogle Scholar
- Tanaka K, Ishii Y, Ogawa J, Shima J. Enhancement of acetic acid tolerance in Saccharomyces cerevisiae by overexpression of the HAA1 Gene, encoding a transcriptional activator. Appl Environ Microbiol. 2012;78:8161–3.View ArticleGoogle Scholar
- Ma C, Wei X, Sun C, Zhang F, Xu J, Zhao X, Bai F. Improvement of acetic acid tolerance of Saccharomyces cerevisiae using a zinc-finger-based artificial transcription factor and identification of novel genes involved in acetic acid tolerance. Appl Microbiol Biotechnol. 2015;99:2441–9.View ArticleGoogle Scholar
- Martani F, Fossati T, Posteri R, Signori L, Porro D, Branduardi P. Different response to acetic acid stress in Saccharomyces cerevisiae wild-type and l-ascorbic acid-producing strains. Yeast. 2013;30:365–78.View ArticleGoogle Scholar
- Bailey JE, Shurlati A, Hatzimanikatis V, Lee K, Renner WA, Tsai PS. Inverse metabolic engineering: a strategy for directed genetic engineering of useful phenotypes. Biotechnol Bioeng. 1996;52:109–21.View ArticleGoogle Scholar
- Oud B, van Maris AJA, Daran JM, Pronk JT. Genome-wide analytical approaches for reverse metabolic engineering of industrially relevant phenotypes in yeast. FEMS Yeast Res. 2012;12:183–96.View ArticleGoogle Scholar
- González-Ramos D, van den Broek M, van Maris AJA, Pronk JT, Daran JM. Genome-scale analyses of butanol tolerance in Saccharomyces cerevisiae reveal an essential role of protein degradation. Biotechnol Biofuels. 2013;6:48.View ArticleGoogle Scholar
- Çakar ZP, Seker UOS, Tamerler C, Sonderegger M, Sauer U. Evolutionary engineering of multiple-stress resistant Saccharomyces cerevisiae. FEMS Yeast Res. 2005;5:569–78.View ArticleGoogle Scholar
- Steiner P, Sauer U. Long-term continuous evolution of acetate resistant Acetobacter aceti. Biotechnol Bioeng. 2003;84:40–4.View ArticleGoogle Scholar
- Wright J, Bellissimi E, de Hulster E, Wagner A, Pronk JT, van Maris AJA. Batch and continuous culture-based selection strategies for acetic acid tolerance in xylose-fermenting Saccharomyces cerevisiae. FEMS Yeast Res. 2011;11:299–306.View ArticleGoogle Scholar
- Sànchezinogué V, Narayanan V, Gorwa-Grauslund MF. Short-term adaptation improves the fermentation performance of Saccharomyces cerevisiae in the presence of acetic acid at low pH. Appl Microbiol Biotechnol. 2013;97:7517–25.View ArticleGoogle Scholar
- Sousa M, Duarte AM, Fernandes TR, Chaves SR, Pacheco A, Leão C, Côrte-Real M, Sousa MJ. Genome-wide identification of genes involved in the positive and negative regulation of acetic acid-induced programmed cell death in Saccharomyces cerevisiae. BMC Genom. 2013;14:838.View ArticleGoogle Scholar
- Nijkamp JF, vanden Broek K, Datema E, de Kok S, Bosman L, Luttik MAH, Daran-Lapujade P, Vongsangnak W, Nielsen J, Heijne WHM, et al. 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 ArticleGoogle Scholar
- Entian KD, Kötter P. Yeast genetic strain and plasmid collections. Method Microbiol. 2007;36:629–66.View ArticleGoogle Scholar
- Verduyn C, Postma E, Scheffers WA, Van Dijken JP. Physiology of Saccharomyces cerevisiae in anaerobic glucose-limited chemostat cultures. J Gen Microbiol. 1990;136:395–403.View ArticleGoogle Scholar
- Hensing MCM, Bangma KA, Raamsdonk LM, de Hulster E, Van Dijken JP, Pronk JT. Effects of cultivation conditions on the production of heterologous α-galactosidase by Kluyveromyces lactis. Appl Microbiol Biotechnol. 1995;43:58–64.View ArticleGoogle Scholar
- 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–17.View ArticleGoogle Scholar
- Akada R, Hirosawa I, Kawahata M, Hoshida H, Nishizawa Y. Sets of integrating plasmids and gene disruption cassettes containing improved counter-selection markers designed for repeated use in budding yeast. Yeast. 2002;19:393–402.View ArticleGoogle Scholar
- 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 ArticleGoogle Scholar
- Kok S, Nijkamp JF, Oud B, Roque FC, de Ridder D, Daran JM, Pronk JT, Maris AJA. Laboratory evolution of new lactate transporter genes in a jen1Δ mutant of Saccharomyces cerevisiae and their identification as ADY2 alleles by whole-genome resequencing and transcriptome analysis. FEMS Yeast Res. 2012;12:359–74.View ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25:1754–60.View ArticleGoogle Scholar
- 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–9.View ArticleGoogle Scholar
- Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010;26:589–95.View ArticleGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998;8:175–85.View ArticleGoogle Scholar
- Ewing B, Green P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998;8:186–94.View ArticleGoogle Scholar
- 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–903.View ArticleGoogle Scholar
- Nijkamp JF, van den Broek MA, Geertman JMA, Reinders MJT, Daran JM, de Ridder D. De novo detection of copy number variation by co-assembly. Bioinformatics. 2012;28:3195–202.View ArticleGoogle Scholar
- Marks VD, Ho Sui SJ, Erasmus D, vander Merwe GK, Brumm J, Wasserman WW, Bryan J, van Vuuren HJ. Dynamics of the yeast transcriptome during wine fermentation reveals a novel fermentation stress response. FEMS Yeast Res. 2008;8:35–52.View ArticleGoogle Scholar
- Coste AT, Ramsdale M, Ischer F, Sanglard D. Divergent functions of three Candida albicans zinc-cluster transcription factors (CTA4, ASG1 and CTF1) complementing pleiotropic drug resistance in Saccharomyces cerevisiae. Microbiology. 2008;154:1491–501.View ArticleGoogle Scholar
- Young ET, Pilgrim D. Isolation and DNA sequence of ADH3, a nuclear gene encoding the mitochondrial isozyme of alcohol dehydrogenase in Saccharomyces cerevisiae. Mol Cell Biol. 1985;5:3024–34.View ArticleGoogle Scholar
- Bakker BM, Bro C, Kötter P, Luttik MAH, Van Dijken JP, Pronk JT. The mitochondrial alcohol dehydrogenase Adh3p is involved in a redox shuttle in Saccharomyces cerevisiae. J Bacteriol. 2000;182:4730–7.View ArticleGoogle Scholar
- Vagnoli P, Bisson LF. The SKS1 gene of Saccharomyces cerevisiae is required for long-term adaptation of snf3 null strains to low glucose. Yeast. 1998;14:359–69.View ArticleGoogle Scholar
- Ye T, García-Salcedo R, Ramos J, Hohmann S. GIS4, a new component of the ion homeostasis system in the yeast Saccharomyces cerevisiae. Eukaryot Cell. 2006;5:1611–21.View ArticleGoogle Scholar
- La Rue J, Tokarz S, Lanker S. SCF Grr1-mediated ubiquitination of GIS4 modulates glucose response in yeast. J Mol Biol. 2005;349:685–98.View ArticleGoogle Scholar
- Fernandes AR, Mira NP, Vargas RC, Canelhas I, Sa-Correia I. Saccharomyces cerevisiae adaptation to weak acids involves the transcription factor Haa1p and Haa1p-regulated genes. Biochem Biophys Res Commun. 2005;337:95–103.View ArticleGoogle Scholar
- Chen Y, Sheng J, Jiang T, Stevens J, Feng X, Wei N. Transcriptional profiling reveals molecular basis and novel genetic targets for improved resistance to multiple fermentation inhibitors in Saccharomyces cerevisiae. Biotechnol Biofuels. 2016;9:1.View ArticleGoogle Scholar
- Grimbergen AJ, Siebring J, Solopova A, Kuipers OP. Microbial bet-hedging: the power of being different. Curr Opin Microbiol. 2015;25:67–72.View ArticleGoogle Scholar
- Abee T, Koomen J, Metselaar KI, Zwietering MH, Besten HMW. Impact of pathogen population heterogeneity and stress-resistant variants on food safety. Annu Rev Food Sci Technol. 2016;7:439.View ArticleGoogle Scholar
- van Heerden JH, Wortel MT, Bruggeman FJ, Heijnen JJ, Bollen YJM, Planqué R, Hulshof J, O’Toole TG, Wahl SA, Teusink B. Lost in transition: start-up of glycolysis yields subpopulations of nongrowing cells. Science. 2014;343:1245114.View ArticleGoogle Scholar
- Stratford M, Steels H, Nebe-von-Caron G, Novodvorska M, Hayer K, Archer DB. Extreme resistance to weak-acid preservatives in the spoilage yeast Zygosaccharomyces bailii. Int J Food Microbiol. 2013;166:126–34.View ArticleGoogle Scholar
- Stebbing ARD. Hormesis—the stimulation of growth by low levels of inhibitors. Sci Total Environ. 1982;22:213–34.View ArticleGoogle Scholar
- Ding J, Bierma J, Smith MR, Poliner E, Wolfe C, Hadduck AN, Zara S, Jirikovic M, van Zee K, Penner MH, et al. Acetic acid inhibits nutrient uptake in Saccharomyces cerevisiae: auxotrophy confounds the use of yeast deletion libraries for strain improvement. Appl Microbiol Biotechnol. 2013;97:7405–16.View ArticleGoogle Scholar
- Saccharomyces genome database http://www.yeastgenome.org/, Accessed March 2016.
- Paget CM, Schwartz JM, Delneri D. Environmental systems biology of cold-tolerant phenotype in Saccharomyces species adapted to grow at different temperatures. Mol Ecol. 2014;23:5241–57.View ArticleGoogle Scholar
- Wang D, Bodovitz S. Single cell analysis: the new frontier in ‘omics’. Trends Biotechnol. 2010;28:281–90.View ArticleGoogle Scholar
- Yin H, Marshall D. Microfluidics for single cell analysis. Curr Opin Biotechnol. 2012;23:110–9.View ArticleGoogle Scholar
- Woyke T, Jarett J. Function-driven single-cell genomics. Microb Biotechnol. 2015;8:38–9.View ArticleGoogle Scholar
- Kuijpers NGA, Solis-Escalante D, Bosman L, van den Broek M, Pronk JT, Daran JM, Daran-Lapujade P. A versatile, efficient strategy for assembly of multi-fragment expression vectors in Saccharomyces cerevisiae using 60 bp synthetic recombination sequences. Microb Cell Fact. 2013;12:47.View ArticleGoogle Scholar
- Güldener U, Heck S, Fiedler T, Beinhauer J, Hegemann JH. A new efficient gene disruption cassette for repeated use in budding yeast. Nucleic Acids Res. 1996;24:2519–24.View ArticleGoogle Scholar
- Solis-Escalante D, Kuijpers NGA, van der Linden FH, Pronk JT, Daran JM, Daran-Lapujade P. Efficient simultaneous excision of multiple selectable marker cassettes using I-SceI-induced double-strand DNA breaks in Saccharomyces cerevisiae. FEMS Yeast Res. 2014;14:741–54.View ArticleGoogle Scholar
- de Kok S, Yilmaz D, Suir E, Pronk JT, Daran J-M, van Maris AJ. Increasing free-energy (ATP) conservation in maltose-grown Saccharomyces cerevisiae by expression of a heterologous maltose phosphorylase. Metab Eng. 2011;13:518–26.View ArticleGoogle Scholar