- Research article
- Open access
- Published:
QTL mapping of the production of wine aroma compounds by yeast
BMC Genomics volume 13, Article number: 573 (2012)
Abstract
Background
Wine aroma results from the combination of numerous volatile compounds, some produced by yeast and others produced in the grapes and further metabolized by yeast. However, little is known about the consequences of the genetic variation of yeast on the production of these volatile metabolites, or on the metabolic pathways involved in the metabolism of grape compounds. As a tool to decipher how wine aroma develops, we analyzed, under two experimental conditions, the production of 44 compounds by a population of 30 segregants from a cross between a laboratory strain and an industrial strain genotyped at high density.
Results
We detected eight genomic regions explaining the diversity concerning 15 compounds, some produced de novo by yeast, such as nerolidol, ethyl esters and phenyl ethanol, and others derived from grape compounds such as citronellol, and cis-rose oxide. In three of these eight regions, we identified genes involved in the phenotype. Hemizygote comparison allowed the attribution of differences in the production of nerolidol and 2-phenyl ethanol to the PDR8 and ABZ1 genes, respectively. Deletion of a PLB2 gene confirmed its involvement in the production of ethyl esters. A comparison of allelic variants of PDR8 and ABZ1 in a set of available sequences revealed that both genes present a higher than expected number of non-synonymous mutations indicating possible balancing selection.
Conclusions
This study illustrates the value of QTL analysis for the analysis of metabolic traits, and in particular the production of wine aromas. It also identifies the particular role of the PDR8 gene in the production of farnesyldiphosphate derivatives, of ABZ1 in the production of numerous compounds and of PLB2 in ethyl ester synthesis. This work also provides a basis for elucidating the metabolism of various grape compounds, such as citronellol and cis-rose oxide.
Background
The wide diversity that can be observed among individuals of the same species is one of the most remarkable aspects of life. Deciphering the mechanisms explaining this phenotypic variety is among the major aims of evolutionists and geneticists. Quantitative genetics has been applied to untangle these issues, and over the last 30 years numerous studies have illustrated the power of these genetic approaches, and in particular quantitative trait locus (QTL) mapping, with the characterization of many genomic regions linked to or containing genes responsible for quantitative variations in a phenotype. These approaches have been extensively used in plant and cattle breeding programs; they have contributed to the understanding of resistance to several diseases [1] and also led to a significant improvement in crop yields and cattle breeding. Surprisingly quantitative genetic approaches have been applied only recently to budding yeast, initially to elucidate various complex mechanisms, including sporulation efficiency [2], thermotolerance [3, 4], and drug resistance [5]. Even more recently, this quantitative approach has been used to decipher complex traits [6, 7] at high resolution [8, 9]. It has also been applied successfully to the analysis of variations in gene expression [10, 11]. The QTL approach is now being used to study features important for the beverages industry, for example wine fermentation [12, 13], sake technological traits [14] and ethanol tolerance for ethanol production [15]. All these studies have implicated defective alleles, of for example AMN1[10], ASP1[16] or ABZ1[11], in the diversity of the phenotypes of segregants. Strain By4741 possesses a defective allele of AMN1 which leads to faster daughter cell separation; wine strain SB possesses a defective allele of the ASP1 gene involved in asparagine catabolism; and S288C possesses a defective allele of ABZ1 that codes for an enzyme which catalyzes the synthesis of 4-amino-4- deoxychorismate from chorismate, a step in the synthesis of paraminobenzoic acid. This defective allele of ABZ1 modulates the fermentation rate by controlling nitrogen utilization [11].
Wine aroma is complex and results from the blending of numerous compounds synthesized by vines, some of which are transformed by yeast, together with compounds directly produced by yeast as a result of its primary metabolism [17, 18]. The metabolic pathways leading to the synthesis of these yeast volatile compounds are numerous and incompletely described. The roles of some of the key genes, such as ATF1 for acetates and EEB1 for ethyl esters, has been demonstrated [19, 20]. Nevertheless, little is known about the factors explaining large strain-to-strain differences in the production of volatile compounds [21–23]. Holistic approaches [24] have given new insights into the roles of various key genes in the diversity of production of some volatile compounds. Further work from the same group highlighted how a few key players, such as transcription factors, may explain some of the differences between strains [25].
To analyze the differences in the production of wine aroma compounds linked to yeast strain diversity, we used QTL analysis with a population of 30 segregants arising from a cross between the laboratory strain S288C and 59A, a spore isolated from the industrial wine strain EC1118. This population of segregants has been genotyped with Affymetrix YGS98 microarrays to obtain a high density genetic map and was used for the first quantitative analysis of transcriptome variations during enological fermentation [11]. We tested this population of segregants for the production of aromatic compounds in two different experimental conditions: synthetic musts mimicking white and red wine fermentations. These analyses enabled us to detect the involvement of eight genomic regions in the production of various volatile compounds explaining 39% to 72% of the diversity. As examples, we characterized the role of two genes by hemizygote analysis and identified another candidate gene by analysis of the phenotype of a deleted mutant. Our findings provide new insights into the genetic architecture underlying the production of wine aroma by yeast.
Results
Each of the 30 segregants was tested in two experimental designs. In the first design (experiment A), white wine fermentation was simulated by fermentation at 20°C in medium with a low lipid content, whereas in the second design, mimicking red wine fermentation [11], the fermentations were run at 28°C in medium with a high lipid content (experiment B).
The fermentation kinetics of the 30 segregants presented significant diversity, from typical wine fermentations to clearly sluggish as observed for S288C (which presented the longest fermentation). Unlike the parental strains, several segregants presented a clear ability to flocculate. This resulted in large and significant diversity in the concentrations of volatile compounds at the end of the alcoholic fermentations. We measured a set of 27 compounds in experiment A, and 33 compounds in experiment B. We performed a principal component analysis to reduce the multidimensional data set of experiment A into three more informative dimensions (Figure 1). The first three axes explained 51% of the global variance (37.6% for axes 1 and 2 in Figure 1A and 33.7% for axes 1 and 3 in Figure 1B) and in this analysis the various compounds are grouped according to chemical family. Ethyl esters and medium chain fatty acids were correlated to the first axis, the various acetates correlated together with the second axis and 2-phenyl ethanol and isoamylalcohol were correlated to axis three. The representation of individual progeny strains in the factorial plan indicated a substantial diversity in the ability to produce volatile compounds. This was especially clear for acetates, as many strains were able to produce more acetates in the media than either of the two parents (S288C and 59A). In addition, some strains more than others metabolized geraniol into citronellol or into the high olfactive impact compound cis-rose oxide. This indicates that the characteristics of the yeast strain have a significant and variable impact on the grape aroma fraction. A similar picture was obtained from the analysis of experiment B.
Genetic analysis of the volatile compounds production
From experiment A, heritability was estimated to be greater than 70% for 21 of 27 compounds, which included the grape aroma compounds geraniol, linalool, citronellol and geraniol acetate.
To identify QTL for these technological features, we performed a linkage analysis with the previously reported genotypes for these progeny [11]. The concentrations of most compounds obtained for the population of segregants did not follow normal distributions (Additional file 1), so we performed linkage analysis with both parametric and non parametric models. We identified four and six regions involved in variations in the production of different compounds in experiments A and B, respectively (Tables 1 and 2, respectively). However, additional factors may contribute to wine aroma production: flocculation is one [26]; and the presence of the ABZ1-S288C allele, which is responsible for large variations in fermentation kinetics [11], may also have an indirect effect. To overcome the potential effects of these factors, we performed a second linkage analysis taking these two factors into account as covariables in the model. This enabled us (i) to improve the significance for some QTL detected after a simple scan (such as for ethyl octanoate), (ii) to detect a genetic effect for additional compounds (ethyl hexanoate) of one region already found and (iii) to detect three and one additional genomic regions in experiments A and B, respectively, for other compounds. The effects of flocculation and of ABZ1 allele on aroma production for each QTL are given as Additional file 2.
For several compounds we were unable to identify any QTL despite a high heritability. This were the case for instance for isoamyl-alcohol and its acetate ester. In the case of isoamyl-alcohol, this might be due to two isomeric compounds (3-methyl-1-butanol and 2-methyl-1-butanol) both being involved. Nevertheless, we detected several regions involved in the diversity of the production of various compounds in the acid, alcohols ethyl ester and isoprenoid chemical families. As a whole, these metabolic QTL (mQTL) explained between 43 and 73% of the metabolite variation.
One region on chromosome XII was identified in both experiments with high Lod score values and explained as much as 46 and 72% of the variations in the production of nerolidol. The same region was identified for other isoprenoids characterized only in design B. Another region, on chromosome XIV, was also detected in both experiments and was associated with various phenotypes: ethyl hexanoate in experiment A and several compounds (including 2-phenyl ethanol and ethyl octanoate) in experiment B. The ABZ1 gene which maps in this region has been reported to be responsible for variations in the rate of fermentation [11]. When ABZ1 polymorphism was used as a covariable, we detected other QTL for more compounds indicating that ABZ1 allelic variations may affect the production of several metabolites. These mQTL are scattered through the genome of strain S288C such that we did not detect a major region associated with all the compounds analyzed (Figures 2 and 3).
The other regions detected for several compounds were each found in only one experiment. These regions also explained a smaller part of the variations in the production of the various volatile metabolites despite a high heritability. Possibly, the production of most of these compounds is under multigenic control and the corresponding regions cannot be detected with such a small population of segregants.
In addition to compounds produced by yeast during alcoholic fermentation, we also studied the fate of geraniol, which is present in grape must at up to 3 mg/L [27]. We did not detect any region explaining variations in the metabolism of geraniol. However, we detected one QTL explaining 43% of the variations in the concentration of citronellol a compound produced from geraniol during alcoholic fermentation. The synthetic pathway for citronellol has not been clearly described and our results may indicate new targets to investigate. Two other QTL explained variations in the content of cis-rose oxide and in the ratio between cis and trans isomers of rose oxide. These QTLs may be technologically interesting as the two isomers of rose oxide present different olfactive thresholds.
Evaluation of the role of various candidate genes in the QTL
PDR8 is responsible for variations in nerolidol production during alcoholic fermentation
The major QTL responsible for variations in the concentration of nerolidol in experiment A and of nerolidol, farnesene and bisabolene in experiment B maps to a short region of 20 kb containing 26 ORFs. Nerolidol, farnesene and bisabolene are all derived from farnesyl diphosphate, an intermediate in isoprenoid and ergosterol biosynthesis (Figure 2): at acidic pH, the instability of the diphosphate group leads to the release of farnesol and its isomer nerolidol.
It seemed likely that the gene involved in the modulation of nerolidol, farnesene and bisabolene production is involved in ergosterol biosynthesis or in farnesol/nerolidol transport because of the size of these molecules. One of the genes mapping in this region is PDR8, a transcription factor that modulates the expression of 16 genes [28] including transporters (AZR1, PDR15, QDR2, YOR1), a gene of the ergosterol biosynthesis pathway (ERG8), and enzymes involved in oxido-reduction processes (CTT1, GTT2, YMR315w). This transcription factor was clearly a good candidate. The nucleotide sequences of the PDR8 genes in strains S288C and 59A show numerous single nucleotide differences. These SNPs generate five non-synonymous substitutions between the Pdr8p proteins in 59A and S288C.
To confirm the involvement of the PDR8 gene in the observed phenotype, we compared two reciprocal hemizygotes between S288C and 59A containing only one of the parental alleles. These hemizygotes presented the different phenotypes observed in the population of segregants, with the enological PDR8 allele of 59A leading to a lower production of nerolidol (Figure 3A). These results are in agreement with those obtained for the whole population (Figure 3B). However, the parental strain S288C produced less nerolidol than most of the segregants, indicating further interactions with the genetic background.
Characterization of PDR8 targets involved in the phenotype
To identify which of the targets of PDR8 explain the observed variations in nerolidol production, we measured the production of nerolidol by the corresponding 16 deletant strains in the By4741 background (Figure 4). The deletion of YOR1 led to an increase of nerolidol production whereas the deletion of four other PDR8 target genes (QDR2, PDR15, GPH1 and YMR135W) led to decreases of nerolidol production similar to that observed after the deletion of PDR8. The genes QDR2 and PDR15 encode transporters that may be involved in the export of nerolidol or derived compounds from the cell. The deletion of two other genes, GPH1 and YMR315W, resulted in a similar reductions in nerolidol production indicating other possible mechanisms. GPH1 is a glycogen phosphorylase required for the mobilization of glycogen, and YMR315W is an oxidoreductase enzyme thay may be involved in the reduction of farnesol (data not shown). ERG8, encoding mevalonate phosphate kinase, is an essential gene for isoprenoid and ergosterol biosynthesis, so it was not possible to conduct the appropriate tests with the deleted haploid strain. Deletion of only one copy of ERG8, in the diploid strain By4743, did not lead to any relevant change so we evaluated the effect of the overexpression of ERG8: no significant increase of the production of nerolidol was detected (data not shown). We did not detect any variation in the expression of QDR2, PDR15, GPH1 and YMR135W associated with the PDR8 allelic form reported in the experiment by Ambroset et al. [11], probably because of the high FDR rate. Therefore, we replaced the PDR8 allele in strain 59A and we compared the expression of these four genes between the strains containing each of the two alleles of PDR8. Quantative PCR (Figure 5) indicated that only QDR2 was more strongly expressed in the strain carrying the S288C-PDR8 allele.
ABZ1 allelic variations affect production of 2-phenylethanol and ethyl esters by yeast during fermentation
The variations in the concentrations of 2-phenyl ethanol and of ethyl hexanoate esters were linked to another mQTL corresponding to a 33 kb region of chromosome XIV. This region overlaps a region involved in differences in fermentation kinetics due to allelic variations of the ABZ1 gene [11]. There are five non synonymous mutations between the S288C and 59A alleles of ABZ1. We compared two reciprocal hemizygotes between strains S288C and 59A containing only one allele of each origin to confirm the role of this gene in 2-phenylethanol production. The hemizygote which carried the enological allele of 59A, produced more 2-phenylethanol than the hemizygote which carried the S288c allele (Figure 6) . The addition of 1 mg/L of p-aminobenzoic acid to the fermentation media suppressed the differences in the rates of fermentation of the two strains; it caused a reduction of only 15% of the difference in the production of 2-phenylethanol (Additional file 3: Table S3), but completely abolished the differences in 2-phenyl acetate production. Abz1p uses chorismate as a substrate, which is also one of the precursors of 2-phenylethanol synthesis. We tested for the effects of the two alleles on the concentration of the various compounds analyzed during mQTL analysis (Table 3). We observed significant effects on the concentrations of many volatile compounds, including ethyl esters, confirming the involvement of ABZ1 in their variations. These results also validate the use of ABZ1 as an additive covariable in the model used to search for mQTLs.
PLB2 allelic variations may affect ethyl ester production
A 60kb region of chromosome XIII was linked with variations in the production of ethyl esters, and we identified two candidate genes with two allelic forms in this region: PLB1 and PLB2. These genes code for phospholipase B which displays transacylase activity in vitro[29]. Plb1p in 59A presents some minor differences to that in S288C, whereas Plb2p of S288C carried a P378A substitution with respect to that in 59A. This proline residue is conserved in other Saccharomyces species and the mutation was not found in other available S. cerevisiae genome sequences. The ΔPLB2-By4741 strain produced much less octanoic ethyl ester than the control (Figure 7), and the difference was greater than that associated with deletion of EEB1, one of the key genes involved in the synthesis of decanoic ethyl ester [20]. These findings are consistent with the involvement of PLB2 in this phenotype. Deletion of PLB2 also led to a decrease in decanoic ethyl ester production and an increase of decanoic acid production, which was not observed after the deletion of EEB1[20].
Polymorphism of the various genes and adaptation
The two major QTL detected in this study, PDR8 and ABZ1, show substantial polymorphism with numerous differences between the allelic forms in S288C and 59A. We investigated whether the differences between the alleles originated from the introgression from a specific lineage by comparing the corresponding alleles from other yeast genome sequences. The phylogeny (Figure 8) reveals that the PDR8 allele of S288C is related to Malaysian or Asian alleles, and apparently one of the closest to its S. paradoxus ortholog; by contrast, the allele in 59A is a typical wine allele. As polymorphism may also result from specific adaptation, we performed a McDonald Kreitman test [34]. This test compares the ratio of nonsynonymous to synonymous polymorphism (intra species) to the ratio of nonsynonymous to synonymous divergence with the nearest species. This ratio is called the neutrality index (NI). An NI lower than one reflects a paucity of nonsynonymous polymorphism relative to nonsynonymous divergence, and is indicative of positive selection; an NI greater than one indicates negative selection of deletorious alleles driving divergence between species or balancing selection. This test was applied to a set of 15 PDR8 alleles from strains isolated from various substrates and NI was 2.30, indicating a significant excess of non-neutral mutations (p value=0.009). This suggests that PDR8 is subject to the accumulation of slightly deleterious mutations that are eliminated by negative selection during speciation, or alternatively that PDR8 presents substantial diversity that might be associated with balanced selection resulting from specific adaptation to different niches.
In contrast with PDR8, the overall phylogeny (Figure 9) revealed that the S288C ABZ1 sequence is related to copies from clinical isolate 322134S and bread strains YS2 and YS4. However, the S288C allelic form of ABZ1 is located at the end of a long branch such that it appears to be the result of the accumulation of numerous mutations. Similarly, the McDonald Kreitman test with a set of 15 ABZ1 sequences from strains isolated from various substrates indicated an excess of non-neutral mutations (NI = 3.00, p value<10-3).
Discussion
We report 13 regions linked to variations in the production of wine volatile compounds. This study is the first demonstration of the potential usefulness of QTL analysis for understanding the origin of the variations in the concentrations of wine aroma compounds and deciphering this “intricate lattice of chemical and biological interactions” [24]. It was not possible to detect QTL for all relevant compounds, despite high heritability. Presumably, the synthesis of many of these compounds is under multigenic control, such that the small size of our segregant population prevented exploration of their complexity. Until now, few key technological traits for alcoholic fermentation have been characterized [11, 13, 14, 16].
Several of the QTL found here are related to terpenoids, which constitute a large family of compounds. They include monoterpenes, which with their corresponding alcohols present useful properties, such as fragrances (in essential oils) or variety aroma (in wines), and even antimicrobial and cancer chemopreventive properties [35]. In yeast, these compounds are synthesized through the mevalonic acid pathway from acetyl-coA, which is converted to isopentenylpyrophosphate (IPP) and its isomer dimethylallyl pyrophosphate (DMAPP), the building blocks of isoprenoids. The main product of this pathway is ergosterol, and geranyldiphosphate and farnesyl diphosphate are intermediate metabolites.
We did not detect any QTL explaining variations in residual geraniol. However, one QTL explained some of the variation in the concentrations of citronellol; this QTL maps to a region of chromosome XIII containing several candidate genes. We did not find a candidate explaining the variations in the concentration of cis-rose oxide in the media or in the ratio between the cis and trans isomers. This compound is significant to wine-making because of its high odor activity [36] and it has been shown recently that yeast can produce cis-rose oxide in wine [37].
We demonstrate that the alleles of PDR8 found in S288C and 59A differently regulate the QDR2 gene responsible for the release of nerolidol into the media. Farnesol and its isomer nerolidol arise from farnesyl diphosphate instability at low pH, like that in the yeast vacuole or in the exocellular medium [38]. Therefore, it is possible that the transporter Qdr2p is responsible of the export of either farnesyl diphosphate or of nerolidol.
Enological strains have a PDR8 allele more divergent from S. paradoxus than that of S288C. The neutrality index [34] we calculated for this gene is incompatible with its neutral evolution: it presented a higher number of replacement polymorphisms than expected under neutral selection. This may be the result of adaptation to different niches or the results of isolation and multiple migrations as suggested by Aa et al. for SSU1[39]. The role of farnesol and nerolidol production by S. cerevisiae is not clear. Under the anaerobic conditions of wine fermentation ergosterol synthesis is blocked, however ERG20 expression correlates with fermentation speed [11] even in a fermentation medium containing ergosterol (experiment B). The synthesis of farnesol diphosphate is essential for (i) the synthesis of other compounds including dolichol which is necessary for cell wall assembly [40, 41], (ii) protein prenylation such as that of Skt5p [42] involved in chitin synthase activity, and (iii) ubiquinone synthesis which may be less important during fermentation. Furthermore, farnesol is a biologically active compound that at concentrations higher than 50 μM inhibits S. cerevisiae and C. albicans growth [43, 44] and at lower concentrations is involved in quorum sensing by C. albicans[45, 46].
In addition to the regions affecting terpenoid production, we also linked several other regions to variations in the concentrations of various volatiles. The ABZ1 gene in one of these regions seems to have the widest effect as its allelic variations affected 2-phenylethanol and ethyl ester synthesis in experiment B, and explained as much as 50% of the variations in the concentration of ethylhexanoate in experiment A. The lower production of 2-phenylethanol was only partially restored by the addition of p-aminobenzoic acid to the fermentation media, and this suggests that this phenotype is not solely the consequence of the substantial effect of the ABZ1-S288C allele on fermentation speed via its impact on nitrogen metabolism [11]. The phylogeny of ABZ1 was clearly different from that of PDR8: the ABZ1-S288C allele is located at the end of a long branch whose branch point is close to that of wine strains (relative to other origins). Again, the neutrality index [34] calculated for this gene was significantly higher than 1. Possibly, the allelic form of ABZ1 has accumulated several deleterious mutations leading to a loss of activity. This feature may explain the very particular phenotypic behavior observed for S288C [47].
The two experimental sets we report generated complementary results. We observed effects of PDR8 and ABZ1 allelic variations in both experiments. Nevertheless, the impact of ABZ1 polymorphism was much less pronounced in experiment A than B, as it explained only variation for ethylhexanoate production after correction for the flocculation effect. We also detected one additional region in experiment A only, and four in experiment B only; presumably the different sets of QTL identified in the two experiments reflect the effects of the different environmental conditions. By considering genes mapping in these regions, we identified PLB2 as possibly involved in the variations of ethyl ester content in experiment B. This gene may have a role complementary to EEB1 in the synthesis of ethyl esters during alcoholic fermentation.
Conclusion
This study shows that linkage analysis can give valuable information about the metabolic pathways involved in production of volatile compounds in yeast, even with only a small population of segregants. We identified the involvement of the PLB2 gene in the metabolism of ethyl esters, and evidenced the role of PDR8 in the release of nerolidol into the media via the regulation of QDR2 expression. We also showed that the weak activity of SC288C Abz1p allele leads to a lower production of many metabolites, including 2-phenyl ethanol, and that this effect was only partially relieved by supplementation with paminobenzoic acid. Other candidate genes are currently being evaluated (i.e. for citronellol synthesis). However, we could not find candidate genes in all regions detected, and despite a high heritability, we did not find any regions associated with the production diversity of many of the compounds considered. Possibly, a larger number of segregants is necessary for a more exhaustive analysis. Our results identify potential new targets for a marker-aided breeding strategy in yeast for the optimization of the production of volatile compounds during fermentation.
Interestingly, our genetic analysis revealed the particular evolution of the PDR8 gene. This may reflect a specific adaptation to wine fermentation conditions, but raises questions about the role of farnesol and nerolidol for S. cerevisiae during fermentation.
Methods
Strains, growth conditions, and fermentation conditions
The two parental Saccharomyces cerevisiae strains studied were the standard strain S288c (MATa; SUC2; gal2) and a haploid derivative of the industrial wine strain EC1118 (HO/ho), herein referred to as 59A (MATa; ho). This strain is prototrophic and has fermentation properties similar to the diploid strain EC1118. The population of 30 segregants obtained from these two parental strains used for QTL analysis have been genotyped after hybridization on high density olignonucleotide microarrays Affymetrix YGS98 oligoarrays.
The strains BY4742 (MATa; his3Δ1; leu2Δ 0; lys2Δ 0; ura3Δ 0) and BY4742Δ ABZ1 (Mata; his3Δ 1; leu2Δ 0; lys2Δ 0; ura3Δ 0; YNR033w::kanMX4), and BY4742 (MATa; his3Δ 1; leu2Δ 0; lys2Δ 0; ura3Δ 0) and BY4742Δ PDR8 (Mata; his3Δ 1; leu2Δ 0; lys2Δ 0; ura3Δ 0; YNR033w::kanMX4) were used for hemizygous constructions.
Allelic replacement at PDR8 in 59A was obtained in three steps: 1) deletion of PDR8 from 59A using the hphMX4 cassette for hygromycin resistance (pAG32). Primers for cassette amplification and verification were obtained form Euroscarf. 2) preparation of a replacement cassette containing PDR8-loxP-kanMX4-loxP by the insertion of loxP-KanMX4-loxP into the terminator of PDR8 in strain S288c (primers are given in Additional file 4: Table S4). 3) replacement of the hphMX4 cassette from 59A PDR8Δ::hph with the PDR8-loxP-kanMX4-loxP replacement cassette from S288c and selection on YPD containing G418 (200 μg.l-1). The loss of hphMX4 cassette was verified by PCR and the absence of growth on hygromycin.
YPD medium was used for precultures at 28°C for 24h in 125 mL flasks with shaking.
Synthetic MS300 medium, which mimics a natural must [48] and [11], was used for fermentation experiments. The first experimental design mimicked white wine fermentation (20°C, low lipid content and containing sitosterol; experiment A). Geraniol, one of the key aroma compounds found in Gewürztraminer wine, was added to study its metabolism during alcoholic fermentation. We also analyzed the production of volatile compounds during fermentation as described in Ambroset et al 2011, which differed by the higher lipid content of the synthetic must and fermentation temperature 28°C (Table 4; experiment B). In some experiments, p-aminobenzoic acid was added to the fermentation media at 1 mg/L to study the effect of ABZ1 alleles in the hemizygotes.
Fermentations were performed in 250 mL flasks equipped with airlocks to maintain anaerobiosis without stirring (design A) and in 1 L fermenters with constant stirring (design B). Small flask fermentations were weighed twice daily and stopped as soon as the daily loss was less than 1% of the expected total loss.
Volatile compounds analysis
Wine aroma compounds were analyzed by the Stir Bar Sorptive Extraction method [49] adapted to our laboratory conditions, with a 1 μL injection volume. The analyses were performed with an Agilent 6890N gas chomatograph equipped with an Agilent 7683 automatic liquid sampler coupled to an Agilent 5975B inert MSD (Agilent Technologies). The gas chomatograph was fitted with a DB-Wax capillary column (60 m × 0.32 mm i.d. × 0.50 μm film thickness, J&W Scientific) and helium was used as carrier gas (1 mL min-1 constant flow). The GC oven temperature was programmed without initial hold time at a rate of 2.7°C min-1 from 70°C to 235°C (hold 10 min). The injector was set to 250°C and used in pulsed splitless mode (25 psi for 0.50 min). The temperatures of the interface, MS ion source and quadrupole were 270°C, 230°C and 150°C, respectively. The mass spectrometer was operated in electron impact ionization mode (EI, 70 eV) and the masses were scanned over a m/z range of 29 – 300 amu. Agilent MSD chemStation software (G1701DA, Rev D.03.00) was used for instrument control and data processing. The mass spectra were compared with the Wiley’s library reference spectral bank
The following compounds were analyzed: isoamyl alcohol, isoamyl acetate, isobutanol, 2-3 butanediol, 2-phenylethyl acetate, ethyl hexanoate, ethyl octanoate, ethyl decanoate, ethyl dodecanoate, ethyl myristate, ethyl palmytate, ethyl laurate, 2-phenyl ethanol, hexanoic acid, octanoic acid, decanoic acid, dodecanoic acid, ethyl 9-decenoate, isoamyl octanoate, 2-phenyl ethyl hexanoate, 2-phenylethyl octanoate, 2-phenylethyl decanoate, ethyl 4 hydroxy butanoate, ethyl 3-hydroxydecanoate, ethyl 3-hydroxyoctanoate, nerolidol, farnesol, 2-3 dihydro farnesol, (E, Z)- or (Z, E)- farnesol (A), (E, E)- farnesol (A), farnesyl acetate, isoamyl octanoate, isoamyl decanoate, isoamyl dodecanoate, methyl oleate, trans-β-farnesene, (Z, E)-α-farnesene, α-bisabolene, ß-bisabolene, (E, E)-α-farnesene, α-terpineol, linalol, citronellol, geraniol, nerol, citronellyl acetate, geranyl acetate, neryl acetate, cis-rose oxide, trans-rose oxide.
Statistical and QTL analysis
Heritability was calculated according to the method of Brem et al. [50]. Statistical analyses were performed using R software version 2.13.1 [51]. QTL analysis was done for each phenotype of the two datasets (experiments A and B) using the genetic map of 1834 markers genotyped previously [11]. The distribution of each phenotype was verified using a Shapiro-Wilk normality test: the normality of the distribution was rejected for 17 of the 29 compounds analyzed in experiment A and 30 of 40 in experiment B, for a threshold of 0.05
As the distribution of most phenotypes was not normal and due to the small sample size, linkage analysis was performed using both parametric and non-parametric models to evaluate the robustness of the parametric model. The parametric model consists of a linkage analysis performed using a normal model with the Haley-Knott regression method implemented in the R/qtl package [52, 53]. As the results of the two analyses were concordant, only the normal analysis is presented.
To overcome the potential effects of flocculation [26] and of the presence of the ABZ1-S288C allele which provokes large variations in fermentation kinetics [11], we performed a second linkage analysis using a normal model with the Haley-Knott regression method, first with flocculation as an interactive covariate, and then with the ABZ1 marker (Chr 14, position 689.4 kb) as an additive covariate. For these regions, a significant effect was indeed observed for both flocculation and the ABZ1-specific markers.
For the three models and the two datasets analyzed, logarithm of odds (LOD) scores were computed for each marker every 2.5cM. An interval estimate of the location of each QTL was obtained as the 1-LOD support interval. The LOD significance threshold was estimated after permutation tests that were replicated 1000 times. The percentage of variance explained by each QTL was estimated from a drop-one-term analysis of results in the global model.
Q-PCR analysis of the expression of PDR8 targets after allelic replacement
Fermentations (900 mL of MS300 medium) were performed in triplicate with strains 59A and S288C-PDR8 59A, and cells were sampled when 70% of the glucose had been fermented. RNA was extracted with trizol as described previously [54]. cDNA was produced by reverse transcription and a 1 in 25 dilution of the resulting cDNA was used for the realtime PCR assays with gene-specific primers and Strategene’s Brilliant II SYBR Green QPCR Master Mix (Santa Clara, CA) and an ABI7300 QPCR machine. Expression levels were measured relative to those of UBC6 and SCR1, both giving similar results.
Sequence analysis and phylogeny
The comparison of the sequences of the 59A and S288C genomes and the differences between them can be found at http://genome.jouy.inra.fr/genyeastrait/[11].
To infer the evolutionary history of ABZ1 and PDR8, we collected their sequences from genomes available at SGD (http://www.yeastgenome.org/). All uncompleted or frameshift-containing sequences where discarded from this set. The phylogenies were inferred with MEGA [55] by the Maximum Likelihood method based on the Kimura 2-parameter model [56]. The trees with the highest log likelihood are shown. The trees are drawn to scale, with branch lengths proportional to the number of substitutions per site. The significance of the Neutrality Index [34] test was calculated using the http://bioinf3.uab.cat/mkt/MKT.asp website.
The list of the sequences used for the two analyses is given in supplementary data (Additional file 5).
Author’s contribution
DS: performed fermentation (experiment A), aroma and QTL analysis , candidate gene search and validation (PDR8, PLB2, …). DS wrote a first draft of the manuscript. CA produced the strain set and performed fermentations (experiment B). CB performed allelic replacement for S288c- PDR8 59A strain fermentations and QPCR expression analysis. PC, analyzed aroma compounds. PD built ABZ1 hemizygote strains and performed fermentations. IS, performed statistical analysis (QTL). JLL performed phylogenic analysis and tests. CE, BB, FK, JLL conceived the study ,designed and coordinated the research JLL wrote the manuscript. All authors analyzed the data. All authors read and approuved the final manuscript.
References
Young ND: QTL mapping and quantitative disease resistance in plants. Annu Rev Phytopathol. 1996, 34: 479-501. 10.1146/annurev.phyto.34.1.479.
Deutschbauer AM, Davis RW: Quantitative trait loci mapped to single-nucleotide resolution in yeast. Nat Genet. 2005, 37: 1333-1340. 10.1038/ng1674.
Sinha H, Nicholson BP, Steinmetz LM, McCusker JH: Complex genetic interactions in a quantitative trait locus. PLoS genetics. 2006, 2: e13-10.1371/journal.pgen.0020013.
Steinmetz LM, Sinha H, Richards DR, Spiegelman JI, Oefner PJ, McCusker JH, Davis RW: Dissecting the architecture of a quantitative trait locus in yeast. Nature. 2002, 416: 326-330. 10.1038/416326a.
Perlstein EO, Ruderfer DM, Roberts DC, Schreiber SL, Kruglyak L: Genetic basis of individual differences in the response to small-molecule drugs in yeast. Nat Genet. 2007, 39: 496-502. 10.1038/ng1991.
Demogines A, Smith E, Kruglyak L, Alani E: Identification and dissection of a complex DNA repair sensitivity phenotype in Baker’s yeast. PLoS genetics. 2008, 4: e1000123-10.1371/journal.pgen.1000123.
Ansel J, Bottin H, Rodriguez-Beltran C, Damon C, Nagarajan M, Fehrmann S, François J, Yvert G: Cell-to-cell stochastic variation in gene expression is a complex genetic trait. PLoS genetics. 2008, 4: e1000049-10.1371/journal.pgen.1000049.
Ehrenreich IM, Torabi N, Jia Y, Kent J, Martis S, Shapiro JA, Gresham D, Caudy AA, Kruglyak L: Dissection of genetically complex traits with extremely large pools of yeast segregants. Nature. 2010, 464: 1039-1042. 10.1038/nature08923.
Cubillos FA, Billi E, Zörgö E, Parts L, Fargier P, Omholt S, Blomberg A, Warringer J, Louis EJ, Liti G: Assessing the complex architecture of polygenic traits in diverged yeast populations. Mol Ecol. 2011, 20: 1401-1413. 10.1111/j.1365-294X.2011.05005.x.
Yvert G, Brem RB, Whittle J, Akey JM, Foss E, Smith EN, Mackelprang R, Kruglyak L: Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors. Nat Genet. 2003, 35: 57-64.
Ambroset C, Petit M, Brion C, Sanchez I, Delobel P, Guérin C, Chiapello H, Nicolas P, Bigey F, Dequin S, Blondin B: Deciphering the Molecular Basis of Wine Yeast Fermentation Traits Using a Combined Genetic and Genomic Approach. G3. 2011, 1: 263-281. 2011.
Marullo P, Yvert G, Bely M, Aigle M, Dubourdieu D: Efficient use of DNA molecular markers to construct industrial yeast strains. FEMS yeast research. 2007, 7: 1295-1306. 10.1111/j.1567-1364.2007.00281.x.
Marullo P, Mansour C, Dufour M, Albertin W, Sicard D, Bely M, Dubourdieu D: Genetic improvement of thermo-tolerance in wine Saccharomyces cerevisiae strains by a backcross approach. FEMS yeast research. 2009, 9: 1148-1160. 10.1111/j.1567-1364.2009.00550.x.
Katou T, Namise M, Kitagaki H, Akao T, Shimoi H: QTL mapping of sake brewing characteristics of yeast. J Biosci Bioeng. 2009, 107: 383-393. 10.1016/j.jbiosc.2008.12.014.
Hu XH, Wang MH, Tan T, Li JR, Yang H, Leach L, Zhang RM, Luo ZW: Genetic dissection of ethanol tolerance in the budding yeast Saccharomyces cerevisiae. Genetics. 2007, 175: 1479-1487.
Marullo P, Aigle M, Bely M, Masneuf-Pomarède I, Durrens P, Dubourdieu D, Yvert G: Single QTL mapping and nucleotide-level resolution of a physiologic trait in wine Saccharomyces cerevisiae strains. FEMS yeast research. 2007, 7: 941-952. 10.1111/j.1567-1364.2007.00252.x.
Loscos N, Hernandez-Orte P, Cacho J, Ferreira V: Release and formation of varietal aroma compounds during alcoholic fermentation from nonfloral grape odorless flavor precursors fractions. J Agric Food Chem. 2007, 55: 6674-6684. 10.1021/jf0702343.
Ferreira V, Lopez R, Cacho JF: Quantitative determination of the odorants of young red wines from different grape varieties. J Sci Food Agric. 2000, 80: 1659-1667. 10.1002/1097-0010(20000901)80:11<1659::AID-JSFA693>3.0.CO;2-6.
Lilly M, Lambrechts MG, Pretorius IS: Effect of increased yeast alcohol acetyltransferase activity on flavor profiles of wine and distillates. Appl Environ Microbiol. 2000, 66: 744-753. 10.1128/AEM.66.2.744-753.2000.
Saerens SMG, Verstrepen KJ, Van Laere SDM, Voet ARD, Van Dijck P, Delvaux FR, Thevelein JM: The Saccharomyces cerevisiae EHT1 and EEB1 genes encode novel enzymes with medium-chain fatty acid ethyl ester synthesis and hydrolysis capacity. J Biol Chem. 2006, 281: 4446-4456. 10.1074/jbc.M512028200.
Wondra M, Berovic M: Analyses of Aroma Components of Chardonnay Wine Fermented by Different Yeast Strains. Food Technology and Biotechnology. 2001, 39: 141-148.
Romano P: Function of yeast species and strains in wine flavour. International Journal of Food Microbiology. 2003, 86: 169-180. 10.1016/S0168-1605(03)00290-3.
Camarasa C, Sanchez I, Brial P, Bigey F, Dequin S: Phenotypic Landscape of Saccharomyces cerevisiae during Wine Fermentation: Evidence for Origin-Dependent Metabolic Traits. PLoS One. 2011, 6: e25147-10.1371/journal.pone.0025147.
Rossouw D, Naes T, Bauer FF: Linking gene regulation and the exo-metabolome: a comparative transcriptomics approach to identify genes that impact on the production of volatile aroma compounds in yeast. BMC Genomics. 2008, 9: 530-10.1186/1471-2164-9-530.
Rossouw D, Jacobson D, Bauer FF: Transcriptional regulation and the diversification of metabolism in wine yeast strains. Genetics. 2012, 190: 251-261. 10.1534/genetics.111.132720.
Pretorius IS: Gene Technology in Winemaking: New Approaches to an Ancient Art. Agric Conspec Sci. 2001, 66: 27-47.
Duchene E, Legras J-L, Karst F, Merdinoglu D, Claudel P, Jaegli N, Pelsy F: Variation of linalool and geraniol content within two pairs of aromatic and non-aromatic grapevine clones. Australian Journal of Grape and Wine Research. 2009, 15: 120-130. 10.1111/j.1755-0238.2008.00039.x.
Hikkel I, Lucau-Danila A, Delaveau T, Marc P, Devaux F, Jacq C: A general strategy to uncover transcription factor properties identifies a new regulator of drug resistance in yeast. J Biol Chem. 2003, 278: 11427-11432. 10.1074/jbc.M208549200.
Merkel O, Fido M, Mayr JA, Prüger H, Raab F, Zandonella G, Kohlwein SD, Paltauf F: Characterization and function in vivo of two novel phospholipases B/lysophospholipases from Saccharomyces cerevisiae. J Biol Chem. 1999, 274: 28121-28127. 10.1074/jbc.274.40.28121.
Akao T, Yashiro I, Hosoyama A, Kitagaki H, Horikawa H, Watanabe D, Akada R, Ando Y, Harashima S, Inoue T, Inoue Y, Kajiwara S, Kitamoto K, Kitamoto N, Kobayashi O, Kuhara S, Masubuchi T, Mizoguchi H, Nakao Y, Nakazato A, Namise M, Oba T, Ogata T, Ohta A, Sato M, Shibasaki S, Takatsume Y, Tanimoto S, Tsuboi H, Nishimura A, Yoda K, Ishikawa T, Iwashita K, Fujita N, Shimoi H: Whole-genome sequencing of sake yeast Saccharomyces cerevisiae Kyokai no. 7. DNA research: an international journal for rapid publication of reports on genes and genomes. 2011, 18: 423-434. 10.1093/dnares/dsr029.
Borneman AR, Desany BA, Riches D, Affourtit JP, Forgan AH, Pretorius IS, Egholm M, Chambers PJ: Whole-Genome Comparison Reveals Novel Genetic Elements That Characterize the Genome of Industrial Strains of Saccharomyces cerevisiae. PLoS Genetics. 2011, 7: e1001287-10.1371/journal.pgen.1001287.
Novo M, Bigey F, Beyne E, Galeote V, Gavory F, Mallet S, Cambon B, Legras J-L, Wincker P, Casaregola S, Dequin S: Eukaryote-to-eukaryote gene transfer events revealed by the genome sequence of the wine yeast Saccharomyces cerevisiae EC1118. Proc Natl Acad Sci. 2009, 106: 16333-16338. 10.1073/pnas.0904673106.
Argueso JL, Carazzolle MF, Mieczkowski PA, Duarte FM, Netto OVC, Missawa SK, Galzerani F, Costa GGL, Vidal RO, Noronha MF, Dominska M, Andrietta MGS, Andrietta SR, Cunha AF, Gomes LH, Tavares FCA, Alcarde AR, Dietrich FS, McCusker JH, Petes TD, Pereira GAG: Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production. Genome Res. 2009, 19: 2258-2270. 10.1101/gr.091777.109.
McDonald JH, Kreitman M: Adaptive protein evolution at the Adh locus in Drosophila. Nature. 1991, 351: 652-654. 10.1038/351652a0.
Oswald M, Fischer M, Dirninger N, Karst F: Monoterpenoid biosynthesis in Saccharomyces cerevisiae. FEMS yeast research. 2007, 7: 413-421. 10.1111/j.1567-1364.2006.00172.x.
Guth H: Quantitation and Sensory Studies of Character Impact Odorants of Different White Wine Varieties. J Agric Food Chem. 1997, 45: 3027-3032. 10.1021/jf970280a.
Koslitz S, Renaud L, Kohler M, Wüst M, Chem F, The C, However L: Stereoselective formation of the varietal aroma compound rose oxide during alcoholic fermentation. J Agric Food Chem. 2008, 56: 1371-1375. 10.1021/jf072518t.
Muramatsu M, Ohto C, Obata S, Sakuradani E, Shimizu S: Alkaline pH enhances farnesol production by Saccharomyces cerevisiae. J Biosci Bioeng. 2009, 108: 52-55. 10.1016/j.jbiosc.2009.02.012.
Aa E, Townsend JP, Adams RI, Nielsen KM, Taylor JW: Population structure and gene evolution in Saccharomyces cerevisiae. FEMS yeast research. 2006, 6: 702-715. 10.1111/j.1567-1364.2006.00059.x.
Juchimiuk M, Pasikowska M, Zatorska E, Laudy AE, Smoleńska-Sym G, Palamarczyk G: Defect in dolichol-dependent glycosylation increases sensitivity of Saccharomyces cerevisiae towards anti-fungal drugs. Yeast (Chichester, England). 2010, 27: 637-645. 10.1002/yea.1803.
Orłowski J, Machula K, Janik A, Zdebska E, Palamarczyk G: Dissecting the role of dolichol in cell wall assembly in the yeast mutants impaired in early glycosylation reactions. Yeast (Chichester, England). 2007, 24: 239-252. 10.1002/yea.1479.
Grabińska KA, Magnelli P, Robbins PW: Prenylation of Saccharomyces cerevisiae Chs4p Affects Chitin Synthase III activity and chitin chain length. Eukaryotic cell. 2007, 6: 328-336. 10.1128/EC.00203-06.
Fairn GD, Macdonald K, McMaster CR: A chemogenomic screen in Saccharomyces cerevisiae uncovers a primary role for the mitochondria in farnesol toxicity and its regulation by the Pkc1 pathway. J Biol Chem. 2007, 282: 4868-4874.
Langford ML, Hasim S, Nickerson KW, Atkin AL: Activity and toxicity of farnesol towards Candida albicans are dependent on growth conditions. Antimicrob Agents Chemother. 2010, 54: 940-942. 10.1128/AAC.01214-09.
Hornby JM, Jensen EC, Lisec AD, Tasto JJ, Jahnke B, Shoemaker R, Dussault P, Nickerson KW: Quorum sensing in the dimorphic fungus Candida albicans is mediated by farnesol. Appl Environ Microbiol. 2001, 67: 2982-2992. 10.1128/AEM.67.7.2982-2992.2001.
Palková Z, Váchová L: Life within a community: benefit to yeast long-term survival. FEMS Microbiol Rev. 2006, 30: 806-824. 10.1111/j.1574-6976.2006.00034.x.
Cubillos FA, Zia A, Gjuvsland A, TJ , Warringer J, Forsmark A, Durbin R, Omholt SW, Louis EJ, Liti G, Moses A, Blomberg A: Trait Variation in Yeast Is Defined by Population History. PLoS Genetics. 2011, 7: e1002111-10.1371/journal.pgen.1002111.
Bely M, Sablayrolles J-M, Barre P: Automatic detection of assimilable nitrogen deficiencies during alcoholic fermentation in oenological conditions. J Ferment Bioeng. 1990, 70: 246-252. 10.1016/0922-338X(90)90057-4.
Coelho E, Rocha S, Delgadillo I, Coimbra M: Headspace-SPME applied to varietal volatile components evolution during Vitis vinifera L. cv. “Baga” ripening. Anal Chim Acta. 2006, 563: 204-214. 10.1016/j.aca.2005.11.018.
Brem RB, Yvert G, Clinton R, Kruglyak L: Genetic dissection of transcriptional regulation in budding yeast. Science (New York, N.Y.). 2002, 296: 752-755. 10.1126/science.1069516.
R Core Team: R: A language and environment for statistical computing, Version 2.15.1. R Foundation for Statistical Computing. 2012
Broman KW, Wu H, Sen S, Churchill G: a.: R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003, 19: 889-890. 10.1093/bioinformatics/btg112.
Broman KW, Sen Ś: b.: A Guide to QTL Mapping with R/qtl. 2009, Germany/New York: Springer-Verlag, Berlin/Heidelberg
Legras J-L, Erny C, LeJeune C, Lollier M, Adolphe Y, Demuyter C, Delobel P, Blondin B, Karst F: Saccharomyces cerevisiae activates two different resistance mechanisms when exposed to octanoic and decanoic acids. Appl Environ Microbiol. 2010, 76: 7526-7535. 10.1128/AEM.01280-10.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Kimura M: A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16: 111-120. 10.1007/BF01731581.
Acknowledgments
This work has been supported by a grant from “Comité Interprofessionnel des Vins d’Alsace”, Colmar, France and from Region Alsace, France.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Electronic supplementary material
12864_2012_4741_MOESM1_ESM.pdf
Additional file 1: Distribution of the different phenotypes for the population of segregants and parental strains.(PDF 104 KB)
12864_2012_4741_MOESM2_ESM.pdf
Additional file 2: Impact of flocculation and ABZ1 on aroma production at the different QTL and variablity of fermentation length.(PDF 397 KB)
12864_2012_4741_MOESM3_ESM.pdf
Additional file 3: Table S3: Effect of the addition of p-aminobenzoic acid on the production of 2phenylethanol and 2-phenylethanolacetate. (PDF 40 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Steyer, D., Ambroset, C., Brion, C. et al. QTL mapping of the production of wine aroma compounds by yeast. BMC Genomics 13, 573 (2012). https://doi.org/10.1186/1471-2164-13-573
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1471-2164-13-573