Reducing phenolic off-flavors through CRISPR-based gene editing of the FDC1 gene in Saccharomyces cerevisiae x Saccharomyces eubayanus hybrid lager beer yeasts

Today’s beer market is challenged by a decreasing consumption of traditional beer styles and an increasing consumption of specialty beers. In particular, lager-type beers (pilsner), characterized by their refreshing and unique aroma and taste, yet very uniform, struggle with their sales. The development of novel variants of the common lager yeast, the interspecific hybrid Saccharomyces pastorianus, has been proposed as a possible solution to address the need of product diversification in lager beers. Previous efforts to generate new lager yeasts through hybridization of the ancestral parental species (S. cerevisiae and S. eubayanus) yielded strains with an aromatic profile distinct from the natural biodiversity. Unfortunately, next to the desired properties, these novel yeasts also inherited unwanted characteristics. Most notably is their phenolic off-flavor (POF) production, which hampers their direct application in the industrial production processes. Here, we describe a CRISPR-based gene editing strategy that allows the systematic and meticulous introduction of a natural occurring mutation in the FDC1 gene of genetically complex industrial S. cerevisiae strains, S. eubayanus yeasts and interspecific hybrids. The resulting cisgenic POF- variants show great potential for industrial application and diversifying the current lager beer portfolio.


Introduction
Although interspecific hybridization (i.e. hybridization between two different species) is believed to be rare in nature, next generation sequencing recently revealed the presence of several interspecific hybrid yeasts within the Saccharomyces yeast clade [1][2][3]. While some hybrids PLOS ONE | https://doi.org/10.1371/journal.pone.0209124 January 9, 2019 1 / 20 a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 have occasionally been isolated from natural habitats, most isolates were obtained from manmade industrial environments [4,5]. The best-known and most studied example of such interspecific hybrid is Saccharomyces pastorianus, the yeast species used for lager beer production. This hybrid yeast originated from a cross between the commonly used ale beer yeast Saccharomyces cerevisiae and the cold-tolerant species Saccharomyces eubayanus [6][7][8]. The resulting S. pastorianus hybrid combines the ability of S. cerevisiae to efficiently ferment sugars in beer wort with the cold-tolerance of S. eubayanus, making it the ideal yeast to perform lager beer fermentations which are typically performed at lower temperatures [6,9,10]. Besides S. pastorianus, other hybrid types are associated with industrial environments, such as S. cerevisiae and S. kudriavzevii hybrids (isolated from ale beer and wine fermentations), and S. cerevisiae X S. uvarum hybrids, which are sometimes isolated from wine and beer fermentations [1,[11][12][13].
The discovery of an increasing number of interspecific hybrids in industrial fermentation processes inspired researchers to mimic the hybridization in the lab, often with the aim of generating new variants that would expand the existing spectrum of industrial yeasts [8]. Hybridization between different strains of S. cerevisiae strains has proven an effective method to generate new variants with interesting industrial properties [14][15][16], and the ability to include non-cerevisiae strains in the breeding schemes further broadens the gene pool and thus the phenotypic diversity of the resulting hybrids. Over the past years, multiple reports describe hybridization between S. cerevisiae strains and S. uvarum [17][18][19], S. eubayanus [19][20][21][22][23][24], S. kudriavzevii [25,26] or S. arboricola [19]. The newly generated interspecific hybrids often show hybrid vigor in their fermentative capacity, broadened temperature tolerance and/or diversified metabolite and aroma profiles. Interspecific hybrids therefore open new routes to address changes in the global fermented beverage market, including the increasing demand for low alcohol, high-flavor and unique products [8,27].
While newly formed interspecific hybrids combine interesting characteristics of their respective parental species, they can also inherit undesired phenotypes that impede their direct implementation in industrial production processes. More specifically, the non-cerevisiae strains are less adapted to industrial fermentations and as a result they often perform poorly in the specific stress conditions imposed by industrial environments. Moreover, the wild yeasts often produce certain undesirable aroma compounds. Arguably the most important drawback of the wild species is the production of phenolic off-flavors (POF), most notably 4-vinyl guaiacol (4VG) [19,21,23]. In beer, 4VG is usually an undesirable yeast metabolite that imposes a very distinct spicy, clove-like flavor. It is produced by yeast through the bioconversion of ferulic acid, present in the endosperm of the malt and barley, to its decarboxylated derivative, 4VG [28,29]. The genetic underpinnings of this phenotype are well-described, and involve the action of Fdc1p and Pad1p [28,29]. The first decarboxylates ferulic acid, while the latter provides a prenylated flavin-mononucleotide (FMN) cofactor of Fdc1p, required for its function.
Several successful strategies have been described to obtain artificial interspecific yeast hybrids that are POF -. First, after the hybrid has been formed, an additional step consisting of a backcross to the POF -S. cerevisiae parent effectively removes the phenotype. However, newly formed interspecific hybrids are usually sterile due to the postzygotic barriers between members of the Saccharomyces clade, which limit interspecific hybrids to a vegetative lifestyle [14,30]. Interestingly though, recent publications have found a way to circumvent this hybrid sterility. A first approach is the generation of allotetraploid interspecific hybrids via rare mating of a diploid S. cerevisiae strain with a diploid non-cerevisiae yeast cell [18,20,31]. The resulting allotetraploid interspecific hybrids can form viable allodiploid spores which could be used for backcrossing with spores of its POFparental strain or with another yeast. This approach allowed Krogerus and coworkers to generate a POFinterspecific yeast hybrid, combining genetic material of three parental strains [20]. Nevertheless, this approach also has some major drawbacks. It relies on a rather complex breeding scheme, on the generation of auxotrophic mutants of the candidate parental yeasts and on two consecutive rounds of breeding that require extensive screening of the segregants.
A second approach involves direct modification of the non-cerevisiae parent to eliminate the POF phenotype using mutagenesis [32]. Segregants of the S. eubayanus parental strain are subjected to UV mutagenesis and subsequently screened to identify POFmutants that can be directly applied in a breeding scheme. However, this strategy also has substantial limitations, most notably the large screenings required to identify positive mutants and the risk of off-target mutations with potential undesired phenotypic effects. Despite the availability of a highthroughput screening method for POF production [33], identification of a POFmutant with no or a very limited number of mutations in other genes would require a screening setup with an even higher throughput.
The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and CRISPR associated protein (Cas9)-based genome editing technology offers a new tool to modify phenotypes of industrial Saccharomyces yeasts [34,35]. For example, a CRISPR-based genome editing strategy successfully reduced urea production in wine yeasts [36] and introduced the hops monoterpenes biosynthesis pathway in an ale brewing yeast, yielding yeast cells that produce certain hop aromas [37]. Moreover, a recent publication describes an optimized CRISPR strategy to alter the genome of the industrial S. pastorianus yeasts CBS1483 and W34/70, which allows to efficiently alter specific phenotypes such as ester production, by knocking out the responsible genes [38].
While CRISPR-based technologies demonstrate a broad spectrum of potential applications, the legislation surrounding the use of cisgenic gene-edited organisms differs widely between different countries across the planet. A recent EU ruling states that these organisms should follow the same guidelines as other genetically-modified organisms. By contrast, other countries, including Brazil, USA, Japan and Argentina, have installed specific guidelines for the use of CRISPR-based gene editing that in some cases allow such modified organisms to be used without following the GM legislation, which greatly increases their industrial potential [39][40][41].
In this study, we report a CRISPR-based gene editing strategy to develop cisgenic POFvariants of genetically complex industrial yeasts and interspecific hybrids by introducing a naturally occurring loss-of-function mutation in the FDC1 gene. We applied and optimized this strategy for S. cerevisiae strains with varying ploidy, a non-cerevisiae species (S. eubayanus), and newly developed interspecific lager yeasts. We show that our strategy allows the introduction of homozygous mutations, resulting in cisgenic mutants that lost the ability to produce POF without undesirable side effects.

Yeast strains used in this study
Yeast parental strains for the generation of interspecific hybrids were selected from a collection of 301 industrial and wild Saccharomyces yeasts, previously characterized by Steensels and coworkers [16]. S. cerevisiae strains SA003 and BE011 were selected based on their POF pheno-and genotype, as well as for their ability to form viable spores [16,42]. Additionally, S. eubayanus strains WL2022 (NPCC1286 [43]) and WL024 (NPCC1292 [43]) were selected based on their temperature tolerance and sporulation efficiency and viability. S. cerevisiae strains BE002, BE014, BE020 and BE074 were selected as candidate strains to test the proposed CRISPR-based genome editing strategy in industrially relevant S. cerevisiae strains. An overview of the used yeast strains is given in Table 1.

Interspecific hybrid generation through spore to spore mating and hybridstate confirmation
Interspecific yeast hybrids were generated by the previously described spore to spore mating technique [23]. First, parental strains were subjected to a random spore isolation protocol [15]. Next, a micromanipulator (MSM-singer instruments) was used to pair two spores on a YPD2% agar plate (2%[wt vol -1 ] Bacto peptone,1%[wt vol -1 ] yeast extract, 1.5%[wt vol -1 ] and 2%[wt vol -1 ] glucose), one from each parental strain. After six to eight hours of incubation at room temperature, the formation of a zygote ('Shmoo') was investigated. Possible hybrids were further purified by restreaking (3x) the strains on synthetic 12˚P wort medium (Light  ). PCR-confirmed interspecific hybrids were prolonged stored at -80˚C to ensure strain purity. Afterwards, generated hybrids were genetically stabilized according to the previously reported stabilization protocol [23].

CRISPR-Cas9 based gene editing
Description plasmid. A S. cerevisiae compatible version of the Candida albicans solo vector CRISPR system, previously described by Vyas and coworkers [45], was used as a platform for the CRISPR-based gene editing [46]. The S. cerevisiae and C. albicans codon optimized Cas9 endonuclease (CaCas9, where the use of the 'CUG' leucine codon is avoided, which is predominantly translated as serine by CTG clade species like C. albicans and S. cerevisiae [45]), as well as the single guide (sg)RNA were introduced into the yeast shuttle vector pR5416, which provides a CEN/ARS sequence for plasmid maintenance in yeast [47]. Both the CaCas9 and the sgRNA are preceded by a constitutive promoter (respectively TEF1p and SNR52p). Species specific sgRNA sequences (S. cerevisiae 5-GGCAAGTACTTACAAACGTA-3'; S. eubayanus 5'-GGCAAGTATTTGCAAACGTA-3') were cloned into the vector as described previously [45].
Repair templates. Double stranded (ds) DNA oligos were created as repair template for the homology directed repair of the induced DSB. Each of the repair templates are 100 nucleotides long, and are centered around the induced DSB in the yeast DNA and contained desired mutations in its sequence. The 100 nucleotides long repair templates were generated via a PCR fill-in of two 60bp long primers with an overlap of 20 nucleotides at their three prime end (primer sequence, see S10 Table).
Transformation. A standard lithium acetate-based yeast transformation protocol was used to transform both the CRISPR plasmid, as well as the repair template into the target strains [33]. Firstly, yeast was grown for one overnight in 5 mL YPD2% growth medium at 30˚C, 200 rpm, after which 1 mL of the pregrowth was transferred to 50 mL YPD2% growth medium and incubated for an extra 4hours (30˚C, 200 rpm). Next, the yeast cell culture was centrifuged (3 minutes at 3000 rpm) and cells were resuspended into 200 μL 0.1 M lithium acetate solution. After 10 minutes incubation at room temperature, 50 μL of the cell culture was mixed with 500 ng plasmid, in which the corresponding sgRNA was cloned and 5 to 25 μg (adjusted protocol) repair template DNA, 300 μL PLI (142 M Polyethylene glycol, 0.12 M lithium acetate, 0.01 M Tris (pH7.5) and 0.001M EDTA) and 5 μL salmon sperm DNA (1mg. mL -1 ) and incubated for 30 minutes at 42˚C. Cells were centrifuged (3 minutes at 3000 rpm) and resuspended in fresh YPD2%, after which cells were recuperated for one overnight on YPD2% agar plates at 30˚C. Selection of plasmid-containing cells was done via replica plating onto selective agar plates (YPD2% + clonat (0.2 μg.mL -1 ), followed by a 2 day incubation at 30˚C. Growing colonies were subjected to a second round of selection on YPD2% + clonat agar plates or immediately plated on YPD2% agar plates (adjusted protocol). After selection, strains were grown for three consecutive rounds onto YPD2% agar plates in order to induce plasmid loss, prior to long term storage at -80˚C.
Introduction of the correct mutation was determined via Sanger sequencing with speciesspecific primers (see S10 Table for an overview of used primers).

Absorbance-based POF measurement
The ability of the yeasts to produce POF was tested via the absorbance based detection method, described previously [33]. Yeasts were inoculated in 150 μL liquid YPD2% growth medium, supplemented with 100 mg.L -1 ferulic acid in a 96 well plate. In each plate, a POF -(W34/70) and a POF + (S. mikatae NCYC2888) control were included. 96-well plates were sealed with an aluminum sticker and incubated for 5 days at 30˚C, 200 rpm. After centrifugation (3 min, 3000 rpm), 100 μL of the supernatant was transferred to a new 96 well plate and remaining concentration of ferulic acid was measured at a wavelength of 325 nm (Tecan Infinite 200 PRO, Switzerland). Yeasts were regarded as POFif the absorbance at 325 nm was above the lower limit of the 90% confidence interval of the POFcontrol (W34/70).
Cell concentration of the pregrowths was measured (BioRad, TC20 automated cell counter, USA), and the calculated amount of cells was used to inoculate 150 mL of an 16˚P wort (17% [wt.vol -1 ] Light spray malt extract, Brewferm, Belgium, supplemented with 0.005 mg.L -1 zn 2+ , autoclaved for 10 minutes at 110˚C) to a starting concentration of 1.5 X 10 7 cells.mL -1 ).
The 250 mL bottles were equipped with a water lock and stirring bar after which they were incubated at 14˚C, agitated at 150 rpm. Weight loss was measured on a daily basis to track fermentation kinetics. Fermentations were stopped when the daily weight loss was equal or less than 0.05 g. Next, the fermentations were cooled on ice to prevent evaporation of the volatile compounds, and samples for chromatographic analysis (HS-GC-FID, Shimadzu corporation), ethanol (Alcolyzer beer ME, Anton Paar GmbH), sulfite and glycerol (Gallery Plus Beermaster, Thermo Scientific) measurements were taken. The leftover fermented medium was used for sensory analysis.

Data analysis and representation
All data analyses and visualization were performed in R [48]. Statistical analyses were conducted within the multicomp package (version 1.4-8 [49]). Figures were generated using the ggplot2 package (version 2.2.1 [50]).

Results
Our goal was to develop and test a CRISPR-based genome editing strategy aimed at modifying the POF phenotype of industrial yeasts, including polyploid and aneuploid yeasts as well as interspecific hybrids. More specifically, we aimed at introducing a loss-of-function SNP mutation in the FDC1 gene that occurs naturally in many domesticated industrial ale beer yeasts [42].

Identification of a CRISPR target region to modify the POF phenotype
The first step in CRISPR-based gene editing is finding an appropriate target region and designing the necessary sgRNA sequence [35]. This region needs to meet certain basic requirements. First, the region should harbor a loss-of-function mutation present in the natural biodiversity of Saccharomyces yeasts. Introducing such mutation in the same species or a closely related one is regarded as cisgenic modification, which will favor its industrial applicability. Second, the region should be highly conserved between different strains and even species, ensuring that the strategy works in several yeast strains, species and hybrids. Third, the region should contain a neighboring PAM sequence, essential for the correct identification by the Cas9 endonuclease of its target site. Fourth, the region should also be unique in the genome to avoid offtarget activity.
Analysis of the genome sequences of industrial POF -S. cerevisiae yeasts, as described by Gallone and coworkers [42], indicates that the disruption of either PAD1 or FDC1, inhibits POF production [28,33]. The majority (73.80%) of POFstrains from the 'Beer 1' lineage [42] share a C to T mutation at position 460 of the FDC1 coding sequence, that replaces a glutamine residue (CAA) by a stop codon (UAA) (Fig 1A). This premature stop codon likely disrupts the protein function, as both the dimerization domain and the catalytic pocket of the protein are not formed anymore (S3 Fig) [51]. Importantly, for the other species within the Saccharomyces species complex with a known FDC1 sequence, the same C to T mutation at position 460 of the FDC1 coding sequence introduces a similar stop codon (UAA) (Fig 1B).
Moreover, there is a high degree of conservation in the DNA sequence surrounding this mutation in all 156 previously sequenced S. cerevisiae strains [42]. Indeed, the sequence of the proposed sgRNA sequence is 100% identical for 155 out of the 156 sequenced S. cerevisiae strains. In addition, a PAM sequence (NGG) could be found in the close proximity of this mutation (nine nucleotides downstream of the mutation), which is crucial for the correct guidance of the CaCas9 endonuclease to its target (Fig 1B) [45]. Importantly, the same PAM sequence, situated nine nucleotides downstream of the target mutation site, can be found in the FDC1 sequence of all other sequenced members of the Saccharomyces species complex, allowing the design of a possible sgRNA sequence for these Saccharomyces species (Fig 1B).
Lastly, possible off-target reactions of the CRISPR system were assessed by blasting the newly designed species-specific 20 nucleotide guide sequences plus possible PAM sequence (NGG) against the genome of the 156 previously sequenced S. cerevisiae yeasts [42], as well as against the de-novo assembly of the S. eubayanus genome [53].This analysis shows that the proposed sgRNA sequences are species-and target unique as no other sequences with more than 85% similarity were detected. Moreover, mismatches or gaps in off-target sequences with the highest similarity occur in the first 13 bp immediately upstream of the PAM sequence, which has been shown to be sufficient to achieve a 100% off-target free gene editing in Saccharomyces yeasts [54].
Together, this indicates that the proposed region (Fig 1) is an appropriate candidate for effective cisgenic CRISPR-based engineering of the POF phenotype in pure and hybrid Saccharomyces species.

Evaluation of CRISPR efficiency in euploid S. cerevisiae and non-cerevisiae yeasts
Gene editing of industrial Saccharomyces yeasts is complicated by two main factors. First, many industrial yeasts are poly-and/or aneuploid, and therefore can contain multiple alleles of the target genes which all need to be modified by the CRISPR system. Analysis of the genomes of industrial ale yeasts shows that these strains have an average ploidy level of 3.52, with some yeasts showing a ploidy level above 4. As POF production is a dominant trait, all alleles of FDC1 or PAD1 need to be deactivated to affect the phenotype. Furthermore, CRISPR protocols are generally optimized for lab strains of S. cerevisiae, and their efficiency for editing non-cerevisiae or mixed genomes can be low [35,46,54].
To assess the efficiency of CRISPR-based gene editing in polyploid genomes, we introduced the desired nonsense mutation in FDC1 of different POF + S. cerevisiae strains with different ploidy levels. Besides a lab strain (haploid and diploid S288c), POF + industrial beer yeasts BE014 and BE020 (diploid) and BE002 and BE074 (triploid) were subjected to the CRISPR transformation. These strains contain two (BE014, BE020) or three (BE002, BE074) functional copies of the FDC1 gene, and at least one functional PAD1 gene [42].
Overall, a decrease in the efficiency of CRISPR gene editing was observed with increasing ploidy levels and genome complexity. The haploid and diploid strains showed a high success rate similar to that of the laboratory strains. Introduction of the mutations was observed in  [33,42]. Bars indicate the presence and the position of the mutation in Fdc1p based on S. cerevisiae s288C reference sequence. Height of the bars indicate the frequency of the mutation across the full collection of POF -S. cerevisiae strains considered. Distribution of the mutation is further dissected across S. cerevisiae subpopulations (circles-colors) and its proportion within each subpopulation is summarized (circle-size) (exact frequencies are reported in S1 Table). Type and position of the mutation in the coding sequence are annotated on top of each bar. Only POF -S. cerevisiae strains that harbor homozygous loss-of-function mutations were included in the analysis. (B) Natural FDC1 loss-of-function point mutation selected for the CRISPR-Cas9 gene editing procedure (yellow cone). The alignment represents a zoom-in of the targeted region in the coding sequence of FDC1 across six Saccharomyces species (-20nt, +20nt from the targeted point mutation). Colored boxes highlight specific areas of the targeted region: the targeted glutamine codon (CAA) that will be replaced by a stop codon (TAA) (yellow box), the PAM sequence (red box) and the guide RNA (green box). The bar-chart represents the alignment of consensus annotation for each position in the targeted region as calculated by Jalview [52].
To evaluate the efficiency of the strategy in non-cerevisiae species, we introduced the same C to T mutation in the FDC1 gene of two POF + diploid S. eubayanus yeasts (WL022 and WL024) (Fig 1). This gene editing was highly successful, with an efficiency of 100% (3/3 and 1/ 1 respectively) for both strains, yielding POFvariants that may be suitable for industrial application [8,32].
Evaluation of the newly formed mutants revealed that the introduction of the homozygous nonsense mutation effectively abolish the ability to produce the unwanted POF aroma (Fig 2).

Gene editing of interspecific hybrid yeasts can induce loss of chromosomal fragments
Next, the CRISPR-based gene editing was evaluated in interspecific hybrids. Therefore, we first generated POF + hybrids between S. cerevisiae and S. eubayanus, after which we tried to remove their POF production by introducing a nonsense mutation in their FDC1 genes using the developed CRISPR strategy.
Initially, the standard CRISPR-based gene editing strategy [46] was used to introduce the selected SNP mutation in the S. eubayanus derived FDC1 allele of the novel generated interspecific hybrids H1 and H2. The S. cerevisiae parent of H1 (BE011) is POFand its FDC1 carries a homozygous nonsense mutation in FDC1 p.W497 � [42]. The S. cerevisiae parent of H2 (SP003) is also POFand its FDC1 is heterozygous for the nonsense mutation p.Q154 � and homozygous for the insertion T>TA at position p.166 [42]. As these mutations differ from the one targeted in the developed CRISPR-based gene editing strategy, it allows to determine the species-specificity of the designed strategy towards the alleles derived from their S. eubayanus and S. cerevisiae parents.
For the first interspecific hybrid, H1, 30 POFvariants were obtained out of 32 tested (94%). Similarly, H2 yielded 21 POFvariants out of 24 (87.5%) (S1 Fig). However, control PCR reaction (primer pair SS_FWSE and SS_RVSE (S10 Table)), which amplifies 490 bp surrounding the DSB induced in the S. eubayanus derived FDC1 gene, did not yield an amplification product. Interestingly, further genetic characterization of this region revealed that the CRISPR editing induced the loss of the region downstream of the targeted region in the S. eubayanus Chr13 (S2 Fig). Genomes of artificial interspecific hybrids are notoriously unstable, and introduction of a double-stranded break likely caused partial loss of the respective chromosome [9,55,56]. Whereas this genomic rearrangement yields POFvariants, it is undesirable as it might have major unwanted pleiotropic consequences for other phenotypes. Therefore, the strategy was further optimized to eliminate this detrimental effect in interspecific hybrids. Increasing the concentration of repair template from 5 to 25 μg during the transformation stage and shortening the selection step for the presence of the CRISPR plasmid from two rounds to only one round of selection, allowed us to overcome this interspecific hybrid genome instability. This strategy yielded POFvariants (Fig 3), with a 12.9 (H1, 4/31) and 6.84% (H2; 8/117) efficiency. None of these variants showed partial chromosome loss. Moreover, none of these variants acquired the targeted mutation in the S. cerevisiae FDC1 allele (confirmed via Sanger sequencing), showing the specificity of the designed CRISPRbased genome editing strategy.

Gene edited variants do not show phenotypic side effects
To determine whether the gene editing procedure introduced any unwanted side effects to the fermentation performance of the yeasts, we evaluated their performance in lab-scale beer fermentations and compared the profiles to those of their respective wild types (WT).
Overall, profiles of the various gene edited variants were similar to those of their WT ancestral strains, with no significant differences in CO 2 production throughout the fermentation (Fig 4A and S2 Table), with the exception of strain WL022, where the mutant showed a faster fermentation at T2 (Anova-test; P-value <0.001) and T3 (P-value <0.01)). Furthermore, no  CRISPR mutants lose their ability to convert ferulic acid to 4VG. Yeast cultures were incubated with an excess of ferulic acid, the precursor for the POF aroma. A decrease in ferulic acid indicates POF aroma formation. Remaining ferulic acid was measured for two control strains (POFcontrol (green); POF + control (red)), S. cerevisiae and S. eubayanus parental strains (dark blue and orange, respectively), Interspecific hybrid H1 and H2 (dark purple) with four of their respective CRISPR variants (light purple). Error bars represent the standard deviation of two biological replicates. The dotted line represents the applied cut off value (lower border of the 90% confidence interval of the POFcontrol).
https://doi.org/10.1371/journal.pone.0209124.g003 At the end of the fermentation, 17 different flavor-active metabolites were evaluated. Apart from the targeted phenotype, namely the production of 4VG (see further), no differences between gene edited variants and the WT were observed (Anova and post-hoc Tukey test; Pvalue > 0.05; S3-S8 Tables), except for S. eubayanus WL024, where a significant difference was measured for one aroma compound, isoamyl alcohol (P-value = 0.004; average isoamyl alcohol production of WL024 is 227.34 mg.L -1 and 253.85 mg.L -1 for WL024_A). Principal component analysis (PCA), based on the production of the remaining 16 metabolites, shows a clear clustering of the different gene-edited variants with their respective WT ancestors. This further confirms that apart from the production of POF aroma, the overall flavor profile of the strains remained unchanged compared to that of the parental WT strains ( Fig 4B).
As expected, all mutants showed a dramatic decrease in 4VG production (Fig 4C, P-values ranging from <0.001 to 0.033; S3-S8 Tables), with the 4VG concentrations remaining below or around the reported flavor threshold of 0.3 mg.L -1 [29].
Subsequent sensorial analysis by a trained panel supported the above-mentioned measurements. Whereas BE014 and BE020 clearly produced clove like off-flavors, the fermentation products obtained with the respective gene edited variants were described as 'very fruity' and 'neutral'. Similarly, mutants of S.eubayanus strains WL022 and WL024 were described as 'slightly fruity' to 'fruity'. H1's aroma was defined as being 'slightly fruity' with 'phenolic' notes. All four gene edited variants of H1 were scored similar by the panel as being 'very fruity', highlighting the masking effect of 4VG on fruitiness. A similar trend was observed for H2. Interestingly, the aromatic contribution of all generated gene edited variants was still remarkably different compared to the aroma produced by the reference lager yeast W34/70, revealing the potential of these hybrids to broaden the aromatic diversity of lager beers, without introducing unwanted clove like phenolic off-flavors.

Discussion
Today's beer market is characterized by an increased demand for niche products and diversification [57]. These market forces have led to an increased interest into novel beer yeasts that can impart new aromas. However, some of the most interesting strains are characterized by the production of the undesirable aroma compound 4VG. Here we describe a new CRISPRbased gene editing strategy that allows to specifically modify 4VG production in various yeasts, including feral S. cerevisiae isolates, non-Saccharomyces strains and interspecific hybrids.
Compared to previous approaches aimed at modifying industrially-relevant phenotypes in yeast [36,38,54], this study introduces a naturally-occurring point mutation, rather than deleting the gene of interest, thereby generating cisgenic variants that are likely exempted from GM regulations within countries like Argentina, USA, Japan and Brazil [40,41,58].
Gene editing efficiency ranged from 100% for lab strains and diploid industrial S. cerevisiae strains to 10% for more complex, poly-and euploid industrial S. cerevisiae yeasts. We also report the first CIRPSR-based engineering of S. eubayanus, which was also highly efficient (100%). The reported gene editing efficiencies are in line with previous reported efficiencies in butyrate, 11 = isoamyl acetate, 12 = ethyl hexanoate, 13 = phenetyl alcohol, 14 = ethyl octanoate, 15 = phenetyl acetate and 16 = ethyl decanoate. (C) 4VG production (measured via HS-GC-FID; mg.L -1 ) of W34/70 (green), S. eubayanus WL022 and WL024 (dark orange) with their respective gene edited variants (light orange), S. cerevisiae BE014 and BE020 (dark blue) with their respective gene edited variants (light blue), as well as for artificial interspecific hybrids H1 and H2 (dark purple) with their gene edited variants (light purple). Error bars represent the standard deviation of two biological replicates. Dotted line represents the flavor threshold of 4VG in beer (0.3mg.L -1 ), as reported in [29].
Although CRISPR-based editing was previously reported to be an efficient route to modify S. pastorianus [38], the Cas9-induced double-stranded DNA break caused an unexpected and undesirable partial chromosome loss [35]. Off-target activity of CRISPR genome editing in haploid or homozygous Saccharomyces yeasts has been shown to be very rare [59]. More recently however, allele specific gene editing in artificial S. cerevisiae X S. eubayanus hybrids was reported to cause loss of heterozygosity (LOH). Specifically, the induced double-stranded break in the S. eubayanus derived locus was not repaired by the provided repair template via homology-directed repair, but rather repaired via LOH [60]. Although further research is needed, the fact that FDC1 is located in the subtelomeric region in S. cerevisiae (Chr IV) and S. euybayanus (Chr13) could favor partial chromosome loss over repair via LOH. One way to reduce such unwanted structural rearrangements is to design the repair template in such a way that the PAM site is inactivated [54]. This prevents continued cutting of the site by the Cas9 endonuclease after successful introduction of the desired mutation [61]. However, this strategy is not ideal for editing organisms that are targeted for food production, as in most cases the resulting mutant would not be cisgenic and thus subjected to GM laws. We therefore modified the gene editing protocol by increasing the amount of repair template and reducing the CRISPR-Cas9 endonuclease activity, which seemed to reduce the unwanted genomic rearrangements and increased repair via HDR. The proposed technique showed a seven to 12% efficiency in specifically introducing a single SNP in the subtelomeric S. eubayanus derived FDC1 gene in novel generated interspecific hybrids.
Overall, some of the generated POFvariants show great potential for industrial application, as their unique aroma profiles are no longer masked by 4VG. Additionally, the CRISPR gene editing strategy described in our study offers a general tool for tuning the characteristics of various aneuploid and non-Saccharomyces yeasts. Specifically, the combination of our gene editing protocol with the ever-increasing number of identified quantitative trait loci (QTL) represents a formidable opportunity to obtain superior industrial yeasts through gene editing [62].
Supporting information S1 Fig. POF phenotype CRISPR variants with standard protocol. Remaining ferulic acid was measured for two control strains (POFcontrol (green); POF + control (red)), Interspecific hybrid one and two (dark green) with their respective POF -(light green) and POF + (dark red) CRISPR variants, obtained via the standard CRISPR protocol. Error bars represent the standard deviation of two biological replicates. Dotted line represents the used cut off value (lower border of the 90% confidence interval of the POFcontrol). (EPS)

S2 Fig. Standard CRISPR protocol causes genetic instability in complex hybrid genomes.
Three primer pairs were designed, each targeting a different part of the S. eubayanus derived Chr13 in the generated interspecific hybrids H1 and H2. Primer pair one and two amplify a region 4kb and 8b upstream of the introduced DBS respectively. Primer pair three amplifies a region 1kb downstream of the introduced DBS on the S. eubayanus derived chromosome of our hybrids (visual representation not depicted to scale). Primer pair three does not yield any of the expected sized PCR product for gene edited variants generated with the standard CRISPR protocol (H1_1 and H2_1) but does for gene edited variants generated with the adjusted CRISPR protocol (H1_A and H2_A), proposed within this manuscript. (EPS) S3 Fig. 3D crystal structure of Fdc1p. Fdc1p acts as a dimer (blue and red form one monomer, cyan and orange, the second one).Blue and cyan parts of the molecules represent the N-terminal part of the protein that it still formed before the stop-gain mutation p.Q154, whilst the red and orange parts disappears when Q154 is replaced by a stop codon [51]. Two 4VG molecules, which are bound to the catalytic pockets of both monomers, are colored in green. . Gel electrophoreses shows both the cerevisiae and eubayanus specific PCR product for hybrids H1 and H2, as well as for S. pastorianus W34/70. The S. cerevisiae control strains WI009, BE001 and SP003 only show the S. cerevisiae specific PCR product, whereas S. eubayanus WL022 only shows the S. eubayanus specific PCR product. (TIF) S1 Table. Overview occurrence of natural loss-of-function mutations in the FDC1 gene across a collection of 76 POF -S. cerevisiae strains. Total frequency is the frequency of the mutations across the complete collection of POF -S. cerevisiae strain collection. Column four to nine show the frequency of each mutation across S. cerevisiae sub-populations, as described earlier [44]. Only POF -S. cerevisiae strains that harbor homozygous loss-of-function mutations were included in the analysis. (PDF) S2 Table. Statistical analysis of the weight loss measured during fermentation between gene edited variants and their respective WT (from time point 1 to end of the fermentation). P-values were obtained using ANOVA. All statistical analysis were conducted in R, with the multcomp package ( � P-value < 0.05; �� P-value <0.01; ��� P-values <0.001). "/" means fermentations were stopped before this time point. (PDF) S3 Table. Statistical analysis of the phenotypic behavior of H1 compared to the H1 geneedited variants. Column two represents the P-values obtained with ANOVA. Column three to twelve represent the obtained P-values of a post-hoc Tukey test. All statistical analyses were conducted in R, within the multcomp package ( � P-value < 0.05; �� P-value <0.01; ��� P-values <0.001). (PDF) S4 Table. Statistical analysis of the phenotypic behavior of H2 compared to the H2 geneedited variants. Column two represents the P-values obtained with ANOVA. Column three to twelve represent the obtained P-values of a post-hoc Tukey test. All statistical analyses were conducted in R, within the multcomp package ( � P-value < 0.05; �� P-value <0.01; ��� P-values <0.001). (PDF) S5 Table. Statistical analysis of the phenotypic behavior of BE014 compared to its geneedited variants. Column two represents the P-values obtained with ANOVA. Column three to twelve represent the obtained P-values of a post-hoc Tukey test. All statistical analyses were conducted in R, within the multcomp package ( � P-value < 0.05; �� P-value <0.01; ��� P-values <0.001). (PDF) S6 Table. Statistical analysis of the phenotypic behavior of BE020 compared to its geneedited variants. Column two represents the P-values obtained with ANOVA. Column three to twelve represent the obtained P-values of a post-hoc Tukey test. All statistical analyses were conducted in R, within the multcomp package ( � P-value < 0.05; �� P-value <0.01; ��� P-values <0.001). (PDF) S7 Table. Statistical analysis of the phenotypic behavior of WL022 compared to its geneedited variant. Column two represents the P-values obtained with ANOVA. Column three to twelve represent the obtained P-values of a post-hoc Tukey test. All statistical analyses were conducted in R, within the multcomp package ( � P-value < 0.05; �� P-value <0.01; ��� P-values <0.001). (PDF) S8 Table. Statistical analysis of the phenotypic behavior of WL024 compared to its geneedited variant. Column two represents the P-values obtained with ANOVA. Column three to twelve represent the obtained P-values of a post-hoc Tukey test. All statistical analyses were conducted in R, within the multcomp package ( � P-value < 0.05; �� P-value <0.01; ��� P-values <0.001). (PDF) S9 Table. Overview aroma and ethanol production from lab scale lager beer fermentation tests. Quantified yeast-related aroma compounds are represented as concentrations (mg.L -1 ), total weight loss as grams (g), ethanol production as volume percentage. Glycerol and SO 2 production are represented as concentrations (g.L -1 and mg.L -1 respectively). H 2 S production capacity is qualitatively indicated (+,+-, -). Lastly, the used score legend for flavors during sensory analysis was: VS = very slightly; S = slightly; V = very; N = neutral; FR = fruity; POF = cloves, phenolic; FRESH = fresh. (PDF) S10