Production of the Fragrance Geraniol in Peroxisomes of a Product-Tolerant Baker’s Yeast

Introduction

The metabolic engineering of microorganisms such as Saccharomyces cerevisiae is a promising tool for the eco-friendly production of aromas, fragrances, biofuels, and pharmaceutics (Zebec et al., 2016). The challenges lie in optimizing the metabolic flux and reducing the toxicity of the desired substances to the production host (Chen et al., 2020). Monoterpenoids, the smallest group of isoprenoids, are of great industrial interest because of their pleasant aroma and odor, and their potential applications as biofuels and pharmaceuticals (George et al., 2015). Geraniol, an acyclic monoterpenoid with a characteristic rose-like odor, is a constituent of plant essential oils like rose, citronella, or palmarosa oil and is widely used as a fragrance material in perfumes, cosmetics or household products. Further, it is commercially used as an insect repellent and natural pest control agent and it exhibits antimicrobial and antitumor properties (Chen and Viljoen, 2010). The industrial production of geraniol exceeds 1000 tons per year (Lapczynski et al., 2008). It is chemically synthesized from pinene or citral and only small amounts are gained from natural resources for the perfume industry by distillation of citronella oil (Surburg and Panten, 2006). The chemical manufacturing processes require multi-step protocols, the use of high temperatures and separation of isomers by fractional distillation. Further, due to seasonal variations and environmental conditions that affect essential oil production and composition, a consistent supply of biomass for the isolation of geraniol from natural resources is not guaranteed. Heterologous production in microorganisms in principle could provide an alternative eco-friendly and potentially cost-effective synthesis of geraniol (Maury et al., 2005). It benefits from energy and resource efficiency and small waste streams. However, so far the production in microorganisms is not industrially relevant, since only small yields are obtained due to low metabolic flux and the toxicity of geraniol to microorganisms such as Escherichia coli or S. cerevisiae (Liu et al., 2016; Zhao et al., 2016).

In this study, we evaluated two approaches for their suitability to be integrated into the engineering toolbox to increase geraniol product titers in S. cerevisiae cells for industrial biomanufacturers: increasing the tolerance of yeast cells to geraniol and compartmentalizing the geraniol producing enzymes into peroxisomes as production sites (Figure 1). Compartmentalization of biosynthetic enzymes offers the possibility of increasing the metabolic flux due to spatial proximity and increased local concentration of enzymes, improved substrate channeling and reduced metabolic crosstalk (DeLoache and Dueber, 2013). Peroxisomes, ubiquitous eukaryotic organelles with a single bilayer membrane, have already proven to be suitable compartments for the production of polyhydroxyalkanoates used for biodegradable plastics as well as for the production of fatty alcohols, alkanes, and olefins used as drop-in fuels (Poirier et al., 2001; DeLoache et al., 2016; Zhou et al., 2016). Recently, peroxisomes were successfully used as compartments for the production of the triterpene squalene (Liu et al., 2020). Here, we tested their suitability for monoterpenoid production. In a first step, S. cerevisiae cells with constantly high numbers of peroxisomes were generated by deleting the peroxin-encoding genes PEX30, PEX31, and PEX32, involved in peroxisome biogenesis (Vizeacoumar et al., 2004), as well as the pexophagy receptor-encoding gene ATG36, involved in peroxisome degradation (Motley et al., 2012a), in all possible combinations. Spot tests revealed that strains lacking the peroxins Pex30 and Pex32 separately and in combination showed an increased tolerance to geraniol. Comparative genomic analysis of these strains identified randomly inserted single-nucleotide polymorphisms (SNPs) resulting in a truncated α-arrestin-like adaptor Bul1, which seems to be responsible for the increase in geraniol tolerance. The double deletion strain pex30Δ/pex32Δ with high peroxisome numbers and high geraniol tolerance as well as the triple deletion strain pex30Δ/pex31Δ/atg36Δ with high peroxisome numbers but low geraniol tolerance were used for the insertion of the geraniol biosynthetic enzymes into peroxisomes and into cytoplasm as control. Whereas the compartmentalization into peroxisomes increased the geraniol titer by up to 13% in both strain backgrounds compared to cytoplasmic localization, the production was increased by up to 63% in the geraniol tolerant strain background pex30Δ/pex32Δ compared to the sensitive strain background pex30Δ/pex31Δ/atg36Δ. The combination of both approaches finally led to an 80% increase in geraniol titer in the tolerant strain with peroxisomal enzyme localization. Combined with other established tools to increase the metabolic flux (Zhao et al., 2017), the usage of geraniol tolerant strains and compartmentalization into peroxisomes are promising approaches to elevate geraniol titers to a level acceptable for industrial applications and to push ahead eco-friendly and cost-efficient production of this industrially valuable natural product.

FIGURE 1

Figure 1. Engineering of tolerant S. cerevisiae cells for compartmentalized geraniol production in peroxisomes. Peroxisome number is increased by deleting the genes for the peroxins Pex30 and Pex32. A truncated Bul1 α-arrestin-like adaptor protein, caused by a point-mutation in the corresponding gene, increases geraniol tolerance. Introduction of the geraniol biosynthetic enzymes into the peroxisomes leads to geraniol production.

Results

Single and Multiple Deletions of PEX30, PEX31, PEX32, and ATG36 Increase the Peroxisome Number of Saccharomyces cerevisiae Cells in Glucose Medium

Saccharomyces cerevisiae cells strictly regulate the peroxisome number in response to environmental signals by the interplay between biogenesis of new and degradation of present peroxisomes (Sibirny, 2016). At a molecular level, the regulation depends on the function of various PEX and ATG genes (Supplementary Figure S1). In order to obtain strains with a constant high number of peroxisomes, the genes for the peroxisome number controlling proteins Pex30, Pex31, and Pex32, as well as for the pexophagy receptor Atg36 were deleted in all possible combinations in LW2591Y, resulting in 15 deletion strains (Supplementary Table S1). A gene coding for the red fluorescence protein mCherry with the C-terminal, peroxisomal SKL amino acid target sequence (PTS1) was integrated (Till et al., 2012). The peroxisomes were quantified with fluorescence microscopy after growth in glucose-containing YPD medium (Figure 2A). All 15 deletion strains showed an increase in peroxisome number compared to the parental strain (Figure 2B and Supplementary Figure S2). The highest peroxisome number with a 5.6-fold increase was observed for pex30Δ/pex31Δ. The lowest number was found for pex31Δ with a 1.4-fold increase in comparison to the control.

FIGURE 2

Figure 2. Strains with all possible combinations of PEX30, PEX31, PEX32, and ATG36 deletions show increased peroxisome numbers. (A) Fluorescence microscopy of the reference S. cerevisiae strain LW2591Y, pex30Δ/pex32Δ, and pex30Δ/pex31Δ/atg36Δ with integrated mCherry-SKL under a constitutive promoter. Yeast strains were grown overnight in liquid YPD medium at 30°C, diluted and used for microscopy. Due to the reiterative recombination system site in the parental strain, all strains carry gfp integrated into the genome. Scale bar = 10 μm. (B) Quantification of the number of peroxisomes per cell detected via fluorescence microscopy in the parental strain LW2591Y and all PEX30, PEX31, PEX32, and ATG36 deletion combinations with introduced mCherry-SKL in YPD medium. For counting of the cells and peroxisomes in total 76 fluorescence images from two biological replicates were taken and the peroxisome/cell ratios were calculated. The 76 ratios are depicted in a Tukey’s box whisker plot. The p-value for each strain was calculated in comparison to LW2591Y (two-sample T-test). ***p ≤ 0.01; **p ≤ 0.05.

Yeast Strains Carrying a Truncated α-Arrestin-Like Adaptor Bul1 Show a Highly Increased Tolerance to Geraniol

Geraniol is toxic for yeast cells and the addition of 200 mg/l geraniol completely inhibits the growth of laboratory reference strains (Zhao et al., 2016). Therefore, the deletion strains were tested for growth deficiencies, which might cause problems when inserting the biosynthesis genes for geraniol. The strains were cultivated on synthetic complete (SC) medium and on yeast extract peptone medium with glucose (YPD) or oleate (YPO) as sole carbon source (Figure 3A). Oleate is metabolized by peroxisomal enzymes and is therefore an indicator for peroxisome functionality (Gurvitz and Rottensteiner, 2006). Whereas the single deletions showed no growth deficiencies on any medium, the double deletion pex32Δ/atg36Δ showed a diminished growth on SC and the triple deletion pex30Δ/pex32Δ/atg36Δ on SC and YPD. On oleate medium (YPO), the double deletions pex30Δ/atg36Δ, pex31Δ/pex32Δ, and pex32Δ/atg36Δ, the triple deletions pex30Δ/pex31Δ/pex32Δ and pex30Δ/pex32Δ/atg36Δ as well as the quadruple deletion pex30Δ/pex31Δ/pex32Δ/atg36Δ were diminished in growth, indicating a reduced functionality of the peroxisomes.

FIGURE 3

Figure 3. The S. cerevisiae strains pex30Δ, pex32Δ, and pex30Δ/pex32Δ with a truncated Bul1 protein are highly tolerant to geraniol. (A) Spot-tests of the reference strain LW2591Y and all single, double, triple, and quadruple deletion mutants of PEX30, PEX31, PEX32, and ATG36 on a synthetic complete medium (SC), a yeast extract peptone medium with glucose (YPD) or with oleate (YPO) as sole carbon-source and on YPD with 200 mg/l geraniol. A serial dilution was used (OD = 10^–1, 10^–2, 10^–3, 10^–4). The plates were incubated at 30°C for 3 days for SC, YPD, and YPD + geraniol and for 6 days for YPO medium. (B) Schematic representation of the truncated Bul1 proteins present in the tolerant strains pex30Δ, pex32Δ, and pex30Δ/pex32Δ. The wildtype Bul1 harbors an N-terminal conserved domain (Bul1_N) with the ubiquitin ligase binding motif PPSY and a C-terminal conserved domain (Bul1_C). pex30Δ carries a SNP from G to T at position 1579 bp of the BUL1 gene, leading to a stop-codon and a truncated Bul1 protein with 526 aa (Bul1^E527*). pex32Δ carries an additional adenine and pex30Δ/pex32Δ an additional thymine at position 2125 bp of the BUL1 gene, leading to a truncated Bul1 protein with 713 aa (Bul1^1–708+5aa). The nucleotide insertions induce a frameshift, leading to the insertion of an early stop-codon in addition to 5 mutated amino acids at the end of the truncated protein.

Toxic effects of geraniol were investigated by cultivation of all strains on YPD with 200 mg/l geraniol (Figure 3A). In agreement with the literature (Zhao et al., 2016), the reference strain did not grow. The deletion strains pex30Δ, pex32Δ, and pex30Δ/pex32Δ showed a highly increased tolerance to geraniol treatment. pex31Δ, atg36Δ, pex30Δ/pex31Δ, and pex30Δ/pex32Δ/atg36Δ showed weak growth, whereas all other mutants did not grow, similar to the reference strain. On plates with 250 mg/l geraniol, the tolerant strains pex30Δ, pex32Δ, and pex30Δ/pex32Δ did not grow anymore (Supplementary Figure S3). In order to determine the reason for the high tolerance in pex30Δ, pex32Δ, and pex30Δ/pex32Δ compared to the reference strain LW2591Y, a genome comparison was conducted, revealing the accumulation of frame shifts within the BUL1 gene resulting in premature stop-codons mutants. A SNP analysis of pex30Δ/pex32Δ showed a thymine insertion at position 2125 in the BUL1 gene, encoding an α-arrestin-like adaptor protein. The induced frameshift leads to the introduction of an early stop-codon and results in a truncated Bul1 protein with 713 aa missing the C-terminal domain (Figure 3B). Sanger sequencing of the BUL1 gene of all PEX/ATG deletion strains revealed that in pex32Δ an adenine nucleotide was inserted at the same position as in pex30Δ/pex32Δ leading to a truncated Bul1 with the same size. In pex30Δ, a nucleotide substitution from guanine to thymine at position 1579 bp of BUL1 leads to the introduction of a stop-codon, resulting in a truncated Bul1 protein with only 526 aa (Figure 3B). A spot-test of a BUL1 deletion strain in BY4741 background on geraniol-containing YPD medium verified that the loss of an intact Bul1 increases the geraniol tolerance (Supplementary Figure S4). pex30Δ and pex32Δ in BY4741 background carrying an intact BUL1 gene showed no increased geraniol tolerance, supporting the result that a loss-of-function Bul1 is responsible for the geraniol tolerance in the PEX strains in LW2591Y background.

In order to analyze and compare the influence of geraniol tolerance and compartmentalization into peroxisomes on the geraniol titer, two strains were selected for the introduction of the geraniol biosynthesis pathway that meet the following criteria: (1) one strain with and one strain without increased geraniol tolerance, (2) high peroxisome numbers, and (3) no growth deficiencies. Accordingly, pex30Δ/pex32Δ with high geraniol tolerance due to the truncated Bul1 protein and pex30Δ/pex31Δ/atg36Δ with low geraniol tolerance were chosen. Peroxisome numbers of both strains were further quantified with fluorescence microscopy as described before on glucose-containing SC medium, which corresponds to the cultivation conditions for geraniol extraction. After one and two days, both strains show a significant increase compared to the reference strain LW2591Y, although the overall peroxisome/cell ratio is reduced compared to YPD medium (Supplementary Figure S5).

Introduction of Geraniol Biosynthesis Pathway Genes Into Geraniol-Tolerant pex30Δ/pex32Δ and Geraniol-Sensitive pex30Δ/pex31Δ/atg36Δ Yeast Strains

Geraniol can be metabolized by S. cerevisiae into geranyl acetate through acetylation and into citronellol through reduction. Therefore, the genes ATF1, encoding an alcohol acetyltransferase, and OYE2, encoding an NADPH oxidoreductase, were deleted in pex30Δ/pex32Δ and pex30Δ/pex31Δ/atg36Δ (Supplementary Figure S6A; Steyer et al., 2013; Brown et al., 2015).

The biosynthetic precursor of geraniol is geranyl pyrophosphate (GPP). In yeast, GPP is naturally produced via the mevalonate pathway and needed for the production of sterols for the cell membrane (Oswald et al., 2007). To increase the flux through the mevalonate pathway, previously reported and well-established strategies were adapted (Figure 4). The following protein-encoding genes were integrated under the control of strong, constitutive promoters into the S. cerevisiae chromosomal genome of pex30Δ/pex32Δ and pex30Δ/pex31Δ/atg36Δ (Supplementary Figure S6A): an N-terminally truncated 3-hydroxy-3-methylglutaryl-coenzyme A (HMG-CoA) reductase 1 without feedback-regulation (trHMG1) (Donald et al., 1997; Rico et al., 2010); the isopentenyl diphosphate (IPP) isomerase Idi1 (Ignea et al., 2011); the negative regulator of RNA polymerase III Maf1, which represses tRNA biosynthesis from IPP (Liu et al., 2013); and the FPP-synthase from Gallus gallus carrying a point mutation for the N144W amino acid substitution (mFPS^N144W), which produces high amounts of GPP (Stanley Fernandez et al., 2000). The resulting strains were named Sen (“Sensitive,” in pex30Δ/pex31Δ/atg36Δ with low geraniol tolerance) and Tol (“Tolerant,” in pex30Δ/pex32Δ with high geraniol tolerance). We have integrated these genes on the chromosomes to enhance genetic stability. In comparison to an expression of the genes from high copy plasmids, this will most probably lead to reduced gene doses and less expression, which will most likely be accompanied by lower geraniol titers. However, it also leads to a higher comparability between the different strains, since exactly one additional copy of the gene is present in all strains. The expression level of the integrated genes was verified with semiquantitative reverse transcription-PCR (Supplementary Figure S7A). Whereas for trHMG1, MAF1, and mFPS^N144W the gene expression was increased in Tol and Sen compared to LW2591Y, the expression level of IDI1 was comparable to LW2591Y in both strains.

FIGURE 4

Figure 4. Schematic representation of the engineered biosynthetic pathway for geraniol in S. cerevisiae. Geranyl pyrophosphate (GPP) is produced in wildtype S. cerevisiae via the mevalonate pathway (dark gray box). One copy of the genes encoding the blue enzymes were introduced into the yeast genome. The genes for the green enzymes were introduced on a high-copy plasmid. For peroxisomal localization of the green enzymes, the peroxisomal target sequence 1 (PTS1), consisting of the amino acids SKL, was attached to the C-terminus. All genes were expressed under constitutive promoters. The genes encoding the red, crossed enzymes have been deleted from the genome.

Finally, the codon-optimized genes for the GPP-synthase (GPPS) from Abies grandis (AgGPPS2) and the geraniol-synthase (GES) from Ocimum basilicum were inserted under strong, constitutive promoters on a high copy plasmid in the strains Tol and Sen (Supplementary Figure S6B). The N-terminal chloroplast targeting sequences of both genes were removed. To compare cytoplasmic (Cyt) and peroxisomal (Per) production of geraniol, both genes were inserted with and without the PTS1 motif for peroxisomal localization. The successful expression of the inserted genes was verified with semiquantitative reverse transcription-PCR (Supplementary Figure S7B). As control strain, an empty vector without GPPS and GES (Neg) was used. The resulting six strains Sen_Cyt, Sen_Per, Sen_Neg, Tol_Cyt, Tol_Per, and Tol_Neg are listed in Table 1 and Supplementary Table S1.

TABLE 1

Table 1. Growth rate μ, maximum cell growth A and absolute geraniol titer of indicated S. cerevisiae strains.

GPP-synthase and geraniol-synthase with their peroxisomal SKL target sequence were both N-terminally tagged with mCherry and analyzed with fluorescence microscopy to examine whether they are accurately located within the cell (Supplementary Figure S8A). Both proteins were found to be localized to the peroxisomes. In case of mCherry-GES-SKL an additional signal in the cytoplasm is visible indicating either an incomplete transport into the peroxisomes or a partially instable mCherry-GES-SKL fusion protein. A Western analysis was performed using an antibody against mCherry (Supplementary Figure S8B). Whereas there is a single signal detectable for mCherry-GPPS-SKL, a band with the correct size is visible for mCherry-GES-SKL plus several smaller bands, indicating degraded protein versions. We assume that the degraded protein bands are responsible for the cytoplasmic signal in the fluorescence microscopy.

Peroxisomal Localization of GPP-Synthase and Geraniol-Synthase Reduces the Yeast Cell Growth Rate Compared to Cytosolic Localization

Geraniol is toxic for S. cerevisiae cells and reduces the growth rate. Growth curves for the geraniol-producing strains Sen_Cyt, Sen_Per, Tol_Cyt, and Tol_Per, as well as for the non-producing control strains Sen_Neg and Tol_Neg were compared (Figure 5A) and the growth rates μ as well as the maximum cell growth A were determined (Table 1).

FIGURE 5

Figure 5. Tolerance to geraniol combined with peroxisomal enzyme localization in yeast strains with high peroxisome numbers yields highest geraniol titer in this study. (A) Growth curves for indicated S. cerevisiae mutants. Data give the means with standard deviation from two biological with each three technical replicates. (B) Relative geraniol titer calculated from absolute values in Table 1. Data give the means with percentage difference from two biological with each two technical GCMS replicates.

The two control strains Sen_Neg and Tol_Neg grow in an equivalent manner with similar growth rates of 0.076 and 0.082, respectively, resulting in a similar maximum cell growth. All geraniol producing strains show a general growth deficiency. Growth rates vary greatly in the sensitive strains Sen_Cyt (μ = 0.059) and Sen_Per (μ = 0.035), resulting in the lowest determined maximum cell growths with A = 0.822 for Sen_Cyt and A = 0.778 for Sen_Per. The growth rates of the two tolerant strains Tol_Cyt (μ = 0.049) and Tol_Per (μ = 0.046) are similar, leading to similar maximum cell growth values with A = 0.851 and A = 0.854, respectively.

Geraniol Tolerance and High Peroxisome Numbers Combined With Peroxisomal Biosynthesis Improve Geraniol Titers up to 80%

For quantification of the geraniol titer, geraniol was extracted from Sen_Neg, Sen_Cyt, Sen_Per, Tol_Neg, Tol_Cyt, and Tol_Per and analyzed by gas chromatography coupled to mass spectrometry (GCMS). All strains except the two control strains with empty vector Sen_Neg and Tol_Neg produced geraniol in detectable amounts (Figure 5B, Table 1, and Supplementary Table S2).

The geraniol sensitive strains Sen_Cyt and Sen_Per produced the overall lowest geraniol levels. Sen_Cyt with cytoplasmic biosynthesis produced 1.53 mg/l geraniol (100%), whereas peroxisomal biosynthesis increased the geraniol titer to 110% in Sen_Per yielding 1.68 mg/l. The geraniol-tolerant strains Tol_Cyt and Tol_Per have an overall elevated geraniol production. In both strains the levels are increased to 159–163% compared to the corresponding Sen_Cyt and Sen_Per strains with 2.44 mg/l in Tol_Cyt and 2.75 mg/l in Tol_Per. Similar to Sen_Cyt and Sen_Per, the geraniol yield is increased to 113% when synthesized in the peroxisomes compared to the cytoplasm. The comparison of Sen_Cyt with 1.53 mg/l and Tol_Per with 2.75 mg/l reveals that the combination of peroxisomal biosynthesis with the usage of geraniol tolerant strains increased the geraniol titer by 80%.

In order to check whether the geraniol titer increase in strains with peroxisomal enzyme localization depends on the number of peroxisomes, we have inserted GPPS and GES with and without the SKL tag as well as the mevalonate pathway genes trHMG1, IDI1, MAF1, and mFPS^N144W in the parental strain LW2591Y, which has a low number of peroxisomes (Figure 2B and Supplementary Figures S2, S5). The resulting strains LW2591Y_Cyt (cytoplasmic GPPS and GES) and LW2591_Per (peroxisomal GPPS-SKL and GES-SKL) produced 1.47 and 0.95 mg/l, respectively (Supplementary Figure S9). This decrease between cytoplasmic and peroxisomal enzyme localization by 35% leads to the assumption that the increased geraniol titers in strains Tol_Per and Cyt_Per with peroxisomal enzyme localization depend on the number of peroxisomes.

Taken together, these results show that in the presence of high peroxisome numbers (i) strains with high geraniol tolerance show an up to 63% increased geraniol production compared to sensitive strains, (ii) peroxisomal biosynthesis increases the geraniol yield up to 13% compared to cytoplasmic production, and (iii) the combination of both approaches increased the geraniol titer by 80%.

Discussion

Geraniol is a highly valuable fragrance material for the cosmetics industry with further useful applications (Chen and Viljoen, 2010). Whereas the industrial chemical synthesis is energy consuming, its extraction from plants suffers from low yields and seasonal variations (Surburg and Panten, 2006). For an eco-friendly production with stable yields and probably low costs, the production of geraniol in microorganisms is an optimal solution (Maury et al., 2005). That said and despite much research effort, monoterpenoid production in microorganisms is limited by low metabolic flux through engineered monoterpenoid production pathways and product toxicity (Shah et al., 2013; Zhao et al., 2016). To circumvent both effects, this study applied compartmentalization of biosynthetic enzymes into peroxisomes in tolerant strains using geraniol production as an example, resulting in an 80% increase of product titers. For several years, peroxisomes have been considered as compartments for heterologous metabolic pathways (DeLoache et al., 2016). Peroxisomes are also naturally used as metabolite production sites, e.g., for the penicillin production in Aspergillus nidulans or for the AK-toxin biosynthesis in Alternaria alternata (Imazaki et al., 2010; Herr and Fischer, 2014), and increasing the number of peroxisomes can increase metabolite titer (Herr and Fischer, 2014). Since AK-toxin is a mevalonate pathway-derived terpenoid, the peroxisomal environment seems to be adequate for the biosynthesis of the monoterpenoid geraniol.

Yeast peroxisomes are generated and degraded depending on nutritional conditions (Sibirny, 2016). On glucose medium, low numbers of peroxisomes are present in a S. cerevisiae cell and excessive peroxisomes are degraded via an autophagic process called pexophagy coordinated by the specific pexophagy receptor Atg36 (Motley et al., 2012a,b; Till et al., 2012). In contrast, on oleate-containing medium, high numbers of peroxisomes are generated via fission or de novo proliferation from the endoplasmic reticulum depending on Pex30, Pex31, and Pex32 amongst others (Nuttall et al., 2011; Joshi et al., 2016). For the industrial production of geraniol or other compounds inside peroxisomes, constantly high peroxisome numbers as well as cultivation in cheap glucose medium would be beneficial. A previous study showed that PEX30, PEX31, and PEX32 deletions in different combinations increased the average number of peroxisomes per cell on oleate medium (Vizeacoumar et al., 2004). Here, we investigated peroxisome numbers during growth on glucose medium (YPD) and in combination with an ATG36 deletion (Figure 2B). In general, we confirmed the increase in peroxisome numbers in strains lacking Pex30, Pex31, or Pex32, but the amount differs between cells grown on oleate or glucose as carbon source. Whereas on oleate medium pex30Δ/pex31Δ/pex32Δ showed the highest peroxisome number (Vizeacoumar et al., 2004), on glucose medium the highest increase was achieved when simultaneously deleting PEX30 and PEX31 (Figure 2B). The additional deletion of ATG36 in most of the PEX deletions did not lead to an additional increase in peroxisome number but to a decrease compared to the PEX deletions alone (Figure 2B). A comparison of the peroxisome functionality of the PEX/ATG deletion library by spotting on oleate-containing medium revealed that 6 out of 15 deletion mutants showed defects, whereas all single deletions grow normally (Figure 3A), which is in accordance with the literature (Vizeacoumar et al., 2004). For pex31Δ/pex32Δ it was shown earlier that the localization of the GFP-tagged peroxisomal membrane protein Pex3 was altered, indicating that the protein import machinery of the peroxisome is not working properly and important enzymes for the function of the peroxisomes are not imported correctly (Zhou et al., 2016). This is in line with the observed growth defects on oleate-containing medium for pex31Δ/pex32Δ. However, mCherry-SKL was correctly imported in all PEX/ATG deletion strains in fluorescence microscopy for peroxisome quantification. Growth on glucose-containing SC medium led to a general decrease of peroxisomes numbers compared to geraniol-containing YPD medium (Supplementary Figure S5). Still, the deletion strains pex30Δ/pex32Δ and pex30Δ/pex31Δ/atg36Δ showed an increased peroxisome number compared to the parental strain LW2591Y, but the increase was not as high as for YPD medium (Figure 2B). Thus, the cultivation of strains in YPD medium would be preferable to obtain maximum peroxisome numbers. Since the GPPS and GES were expressed from plasmids, SC medium without uracil as selection marker was chosen to prevent plasmid loss.

A major drawback in engineering microorganisms as production platforms for monoterpenoids such as geraniol is the toxicity of them to the cells. At a concentration of 200 mg/l geraniol, growth of the laboratory strain CEN.PK102-5B is completely inhibited (Zhao et al., 2016). This is in line with our laboratory strains (Figure 3A and Supplementary Figure S4). The reason for its toxicity might be the membrane-disruptive features of geraniol (Bard et al., 1988; Carnesecchi et al., 2002). One way to circumvent toxicity is to increase the tolerance of the strains. High tolerance was achieved by SNPs in the BUL1 gene leading to truncated protein variants in pex30Δ, pex32Δ, and pex30Δ/pex32Δ (Figure 3A). Bul1 is an α-arrestin-like adaptor that links the E3 ubiquitin-ligase Rsp5 to its substrate for ubiquitination. This process is needed during endocytosis of plasma membrane transporters to adapt to changing external cues (Novoselova et al., 2012). Up to date, two ubiquitination targets of the Bul1/Rsp5 complex are identified: the monocarboxylate transmembrane transporter Jen1 and the general amino acid permease Gap1 (O’Donnell, 2012). Growth on either a preferred nitrogen source such as ammonium or the carbon source glucose induces Bul1/Rsp5-mediated ubiquitination of Jen1 or Gap1, respectively, and stimulates transporter inactivation via endocytosis. Consequently, a non-functional Bul1 protein prevents transporter degradation and might lead to increased tolerances to chemical compounds. Indeed, the BUL1 deletion strain is not only tolerant to geraniol, but also to a broad spectrum of other chemicals such as cisplatin, fluorine-containing anesthetics, camptothecin, tunicamycin, bleomycin, or rapamycin (Wolfe et al., 1999; Huang et al., 2005; Kapitzky et al., 2010). Since the truncated Bul1 strains as well as the deletion strain show no fitness defects on glucose medium but a high tolerance to geraniol (Figure 3A and Supplementary Figure S4), the non-functional Bul1 is a promising candidate for being incorporated as integral part into the geraniol engineering toolbox in yeast. Here, we gave the proof-of-principle and increased the geraniol titer by 63% in a strain with a Bul1 truncation (Figure 5B). Because of its manifold tolerances to chemicals, BUL1 deletions might also be considered for improving the heterologous production of other valuable chemical products in yeast. The reason for the weak geraniol tolerance in pex31Δ, atg36Δ, pex30Δ/pex31Δ, and pex30Δ/pex32Δ/atg36Δ could not be identified yet (Figure 3A). The genes for BUL1 were intact in these strains. Possibly, another not yet identified SNP is responsible for the increased tolerance. The future identification of this putative SNP might lead to the identification of a new candidate for further increasing geraniol tolerance in yeast.

For the compartmentalization of enzymes, a suitable targeting signal is necessary for the correct transport into the organelle. We have chosen the PTS1 for peroxisomal localization consisting of the three amino acids SKL attached to the C-terminus of the protein. There are three different common tags for peroxisomal localization. The SKL and the KANL tags are C-terminal PTS1 tags and bind the PEX5 receptor, whereas the N-terminal PTS2 tag binds the Pex7 receptor (Ehrenworth et al., 2017). Among those three tags, the SKL tag was shown to have the highest efficiency with 99%, which is almost independent of plasmid copy number (Ehrenworth et al., 2017). To verify the correct import of the SKL-tagged GPPS and GES in our study, we have tagged the proteins with an N-terminal mCherry protein and analyzed the localization with fluorescence microscopy (Supplementary Figure S8A). Both proteins were visible in the peroxisomes. For mCherry-GES-SKL, an additional signal in the cytoplasm was visible, which was assigned to degraded protein versions without the C-terminal SKL tag (Supplementary Figure S8B).

Compartmentalization of enzymes in general benefits from spatial proximity and increased local concentration of enzymes, improved substrate channeling and reduced metabolic crosstalk (DeLoache and Dueber, 2013). In our study, a 13% increase in geraniol production was observed when inserting the GPPS and GES into the peroxisomes in strains with high peroxisome numbers. In a control strain with low peroxisome numbers, the geraniol titer was decreased in strains with peroxisomal enzymes compared to cytosolic enzymes (Supplementary Figure S9). Although the increase of 13% in strains with high peroxisomes is small and other factors such as the different genetic background caused by the different PEX/ATG deletions must be considered, a trend can be seen that compartmentalization into peroxisomes is a promising approach to increase geraniol titers. For further optimizations of the titers, the whole MVA pathway can be introduced into the peroxisomes to further improve spatial proximity and local concentrations of precursors as described previously (Liu et al., 2020). The aim of this study was to test new methods for their suitability to be included in the engineering toolbox for geraniol production in yeast. Compared with the literature, our obtained titers are low. We have accepted the low titers in order to achieve a better comparability of the strains. Instead of overexpressing the genes for mevalonate pathway flux optimization on plasmids with several copies per cell, we have integrated only one copy of each gene in a defined spot on the chromosomes. This enables an even expression level and increases comparability between the strains, but the expression levels especially for IDI1 and HMG1 were low (Supplementary Figure S7A). The highest geraniol titer obtained from S. cerevisiae cultures is reported as 1.69 g/l (Zhao et al., 2017). This was achieved by combined deletion of the NADPH oxidoreductase gene OYE2 to reduce endogenous geraniol metabolism coupled with control of the GPP flux distribution by dynamic control of the FPP-synthase gene ERG20 in a leucine prototrophic strain. The strain, which carries the Valeriana officinalis GES, was cultivated with pure ethanol feeding in fed-batch fermentation (Zhao et al., 2017). Combinations of these previous approaches with the peroxisomal compartmentalization of geraniol biosynthetic enzymes contained in geraniol-tolerant S. cerevisiae strains might further enhance geraniol production. A comparison of the growth rates showed that strains with peroxisomal enzyme localization are reduced compared to the strains with cytosolic enzymes (Figure 5A and Table 1), which might be less optimal for industrial applications. Nevertheless, the tolerant strain with peroxisomal enzymes Tol_Per reaches a similar maximum cell growth after 30 h as the other strains with cytoplasmic enzymes and produced the highest geraniol titer in our study. Titers likely depend on overall available peroxisomal space (Herr and Fischer, 2014). Further increases in production might be reached using the methylotrophic yeasts Pichia pastoris and Hansenula polymorpha as both have giant peroxisomes that can account for more than 80% of the total cell volume (Veenhuis et al., 1978; Van Der Klei et al., 1991). These general chassis engineering approaches could be extended to other monoterpenoids many of which are important fragrances, flavors, and next-generation biofuels. Compartmentalization in product-tolerant strains will work alongside other strategies to increase product titers to industry production levels. Our work establishes that metabolic compartmentalization and, especially, the usage of product-tolerant strains are important strategies that should be considered for inclusion in the chassis engineering toolbox to increase product titers for industrial biomanufacture.

Materials and Methods

Growth Media and Growth Conditions

Saccharomyces cerevisiae strains were cultivated on either synthetic complete medium [SC; 0.2% yeast nitrogen base w/o amino acids and ammonium sulfate, 0.5% ammonium sulfate, 2% glucose, 5% (v/v) amino acid mix, and the appropriate supplements 0.002% histidine, 0.002% tryptophan, 0.002% uracil, 0.009% leucine], yeast extract peptone dextrose medium (YPD; 1% yeast extract, 2% bacto peptone, and 1% glucose), or an oleate-containing yeast extract peptone medium [YPO; 0.3% yeast extract, 0.5% bacto peptone, 0.12% (v/v) oleic acid, 0.2% (v/v) Tween-40, 0.5% KH₂PO₄ (pH = 6)]. The amino acid mix contains 0.04% adenine sulfate, 0.04% arginine-HCl, 0.2% aspartic acid, 0.2% glutamic acid, 0.06% isoleucine, 0.06% lysine-HCl, 0.04% methionine, 0.1% phenylalanine, 0.8% serine, 0.4% threonine, 0.06% tyrosine, and 0.3% valine. For recycling of the URA3-marker 0.05% 5-fluoroorotic acid was used. For solid plates 2% agar was used. For geraniol toxicity tests, 200 mg/l geraniol was added to the YPD medium.

Escherichia coli strains were cultivated on lysogeny broth medium (LB; 1% tryptone, 0.5% yeast extract, and 1% NaCl) supplemented with 100 μg/ml ampicillin.

Codon-Optimization and Synthesis of Genes

For gene synthesis, a geranyl diphosphate-synthase of Abies grandis (GenBank: AAN01134.1) without the chloroplast targeting sequence and with a C-terminal SKL-tag, a gene encoding a geraniol-synthase of Ocimum basilicum (UniProtKB/Swiss-Prot: Q6USK1.1) without the chloroplast targeting sequence and with a C-terminal SKL-tag, and a gene encoding the farnesyl pyrophosphate-synthase of Gallus gallus (UniProtKB/Swiss-Prot: P08836.2) with a N144W mutation (mFPS^N144W) were chosen. All genes were codon-optimized for S. cerevisiae with the GeneOptimizer^TM Software, synthesized and cloned with PmeI restriction sites into vector pMA-RQ (Thermo Fisher Scientific GENEART GmbH, Regensburg), resulting in plasmids pME4804–4806. In addition, the HIS3 promoter and terminator, the CYC1 promoter and terminator, and the TEF1 promoter and terminator were synthesized and cloned with SfiI restriction sites into the vectors pME4805, pME4806, and pME4804, respectively.

Plasmids Construction

All used and constructed plasmids are listed in Supplementary Table S3. All used primer sequences are given in Supplementary Table S4. E. coli strain DH5α [F^–, Φ80Δ(lacZ)M15^–1, Δ(lacZYA-argF) U169, recA1, endA1, hsdR17 (r_K^–, m_K⁺), supE44, λ^–, thi1, gyrA96, relA1] was used (Woodcock et al., 1989).

For targeting a protein to the peroxisomes, the peroxisomal target sequence 1 (PTS1) consisting of the amino acids SKL was fused to the C-terminus of the protein. For this, the corresponding coding sequence was attached to the reverse primer of the gene of interest.

For the construction of pME4803, P_GPD1 (amplified with JG965/1062 from pME4093), mCherry-SKL (amplified with JG1063/1064 from pME3772), and T_CYC1 (amplified with JG1065/969 from pME2787) were assembled together with pBluescript II KS (+) digested with EcoRV via the GeneArt^TM Seamless Cloning and Assembly Enzyme Mix Kit (Thermo Fisher Scientific).

For pME4804, the TEF1 promoter, mFPS^N144W, and the TEF1 terminator were synthesized and cloned into the pMA-RQ vector with additional SfiI sites (Thermo Fisher Scientific). In the same manner P_HIS3, GES-SKL, and T_HIS3 (pME4805) as well as P_CYC1, GPPS-SKL, and T_CYC1 (pME4806) were designed (Thermo Fisher Scientific).

For pME4807, the URA3 promoter and terminator were amplified from pME2787 with primers JG1228/1230 and JG1231/1229, respectively. The truncated HMG1 gene (trHMG1) was amplified with primers JG1232 and JG1233 from genomic DNA of LW2591Y. All three fragments and the EcoRV digested pBluescript II KS (+) vector were assembled with the GeneArt^TM Seamless Cloning and Assembly Kit (Thermo Fisher Scientific) resulting in plasmid pME4807.

For pME4808, the IDI1 gene was amplified with primers JG1224 and JG1225 from genomic DNA of LW2591Y and cloned with the GeneArt^TM Seamless Cloning and Assembly Kit (Thermo Fisher Scientific) into pME4806 digested with PmeI, resulting in pME4808, containing the IDI1 gene with the CYC1 promoter and terminator.

For pME4809, the MAF1 gene was amplified with primers JG1226 and JG1227 from genomic DNA of LW2591Y and cloned with the GeneArt^TM Seamless Cloning and Assembly Kit (Thermo Fisher Scientific) into pME4805 digested with PmeI, resulting in pME4809, containing the MAF1 gene with the HIS3 promoter and terminator.

For pME4810, P_GPD1 (amplified with JG1355/1172 from pME4803), GPPS-SKL (amplified with JG1220/1221 from pME4806), and T_CYC1 (amplified with 1137/1346 from pME4803) were assembled together with PvuII digested pESC_URA via the GeneArt^TM Seamless Cloning and Assembly Enzyme Mix Kit (Thermo Fisher Scientific). The plasmid carries an additional PmeI site downstream of T_CYC1. P_GPD1, GPPS (amplified with JG1220/1234 from pME4806), and T_CYC1 were assembled in the same way, resulting in pME4811.

For pME4812, P_TEF1 was amplified with JG1384/1385 from pME4804, GES-SKL was amplified with JG1269/1389 from pME4805, and T_TEF1 was amplified with JG1390/1388 from pME4804. All three fragments were assembled with pME4810 digested with PmeI via the GeneArt^TM Seamless Cloning and Assembly Enzyme Mix Kit (Thermo Fisher Scientific). In the same way P_TEF1, GES (amplified with JG1269/1386 from pME4805), and T_TEF1 were assembled with pME4811 digested with PmeI, resulting in pME4813 (Supplementary Figure S6B).

For pME4814, pESC_URA was digested with PvuII and the band with a size of 5.3 kb was religated with T4 DNA ligase (Thermo Fisher Scientific).

For pME5088, P_GPD1-mCherry (amplified with JG965/1811 from pME4803), and GPPS-SKL-T_CYC1 (amplified with JG1913/969 from pME4812) were assembled with pBluescript II KS (+) digested with EcoRV via the GeneArt^TM Seamless Cloning and Assembly Kit (Thermo Fisher Scientific).

For pME5089, mCherry was amplified with JG1910/1811 from pME4803 and GES-SKL was amplified with JG1911/1912 from pME4812. pME4804 was digested with PmeI and mCherry and GES-SKL were cloned between P_TEF1 and T_TEF1 via the GeneArt^TM Seamless Cloning and Assembly Kit (Thermo Fisher Scientific).

Saccharomyces cerevisiae Strain Construction

Saccharomyces cerevisiae strain LW2591Y [reiterative recombination parental acceptor strain; MAT-a inc; his3Δ200; leu2Δ, met15Δ, trp1Δ63, ura3Δ, P_PYK-GFP-HIS3-(HO cleavage site)] (Wingler and Cornish, 2011) was used as background strain. S. cerevisiae strains generated and used are listed in Supplementary Table S1.

Gene deletion was performed seamlessly according to Akada et al. (2006). Primers are listed in Supplementary Tables S4, S5. The 0.79 kb long 5′-region of the gene of interest was amplified from genomic DNA of LW2591Y. Additionally, the URA-marker was amplified from plasmid pME2787 with a forward primer, containing a 40 nt overhang matching the 3′-region of the gene of interest plus a minimum of 24 nt overhang matching the 5′-region of the gene of interest, and a reverse primer, containing a 40 nt overhang matching the 3′-region of the gene of interest. Both amplicons were fused via fusion-PCR (Szewczyk et al., 2007), resulting in a PCR-fragment, which can be used directly for transformation. For transformation, the protocol described in Wingler and Cornish (2011) was used.

For the introduction of P_GPD1-mCherry-SKL-T_CYC1, P_GPD1-mCherry-GPPS-SKL-T_CYC1, P_TEF1-mCherry-GES-SKL-T_TEF1, P_URA3-trHMG1-T_URA3, P_CYC1-IDI1-T_CYC1, P_HIS3-MAF1-T_HIS3, and P_TEF1-mFPS^N144W-T_TEF1 into the genome, the reiterative recombination system (Wingler and Cornish, 2011) was used. In a first step, the PCR-products for the transformation were amplified. For this, two PCRs with the primers and templates given in Supplementary Table S6 were performed. The template for PCR 2 is the amplicon of PCR 1. The amplicon of PCR 2 can be directly used for transformation together with the corresponding linearized plasmid also given in Supplementary Table S6. The transformation via electroporation, the induction, the selection, and the reiteration were performed according to the protocol described in Wingler and Cornish (2011).

SNP Analysis of Geraniol-Tolerant pex30Δ/pex32Δ Strain

Genomic DNA was extracted from overnight cultures of LW2591Y and pex30Δ/pex32Δ grown in YPD medium. The cells were lysed in 500 μl lysis buffer (30 mM Tris-HCl, 25.5 mM Na₂-EDTA, 0.5% SDS, pH = 8) with glass beads by shaking for 15 min. Proteins and cell debris were removed by the addition of 100 μl 8M potassium acetate and centrifugation for 15 min at 13,000 rpm. The DNA in the supernatant was precipitated with 300 μl isopropanol. After centrifugation at 13,000 rpm for 15 min and removal of the supernatant, the DNA was resuspended in 400 μl water with 3 μl RNAse A (10 mg/ml) and incubated at 37°C for 10 min. DNA was precipitated with 1 ml ethanol and 10 μl 4 M ammonium acetate, washed with ethanol, dried at 65°C and dissolved in 50 μl water at 65°C. For additional DNA purification, the DNeasy PowerClean Pro Cleanup Kit (Qiagen) was used according to the manufacturer’s instructions, if necessary. Concentration and purity of the isolated DNA was first checked with a Nanodrop ND-1000 (PeqLab) and exact concentration was determined using the Qubit^TM dsDNA HS Assay Kit as recommended by the manufacturer (Thermo Fisher Scientific). Illumina shotgun libraries were prepared using the NEBNext Ultra II DNA Library Prep with Beads (New England Biolabs) following the manufacturer’s protocol. To assess quality and size of the libraries, samples were run on an Agilent Bioanalyzer 2100 using an Agilent High Sensitivity DNA Kit as recommended by the manufacturer (Agilent Technologies). Concentrations of the libraries were determined using the Qubit^TM dsDNA HS Assay Kit as recommended by the manufacturer (Thermo Fisher Scientific). Sequencing was performed on a MiSeq system with the reagent kit v3 with 600 cycles (Illumina) as recommended by the manufacturer resulting in 6,729,204 paired-end reads for LW2591Y and 9,020,510 paired-end reads for RH3716 with average read lengths 260 and 269 bp, respectively. Reads were quality filtered with fastp v0.20.0 (Chen et al., 2018) with including overlap correction, removal of sequences with qualities <20, soft clipping with sliding window of 4, removal of sequences <50 bp, and Illumina adapter removal. Quality filtering resulted in 6,525,536 and 8,709,042 paired-end reads with average read lengths 257 and 267 bp. In addition, potential adapter remains were removed with cutadapt v2.5 (Martin, 2011) and potential phiX contamination was removed with bowtie2 v2.3.5.1 (Langmead and Salzberg, 2012) employing default parameters. SNP analysis of LW2591Y and RH3716 was performed with breseq v0.35.0 (Deatherage and Barrick, 2014) against the reference strain S288C (GCA_000146045.2) with default parameters. Sequence data were deposited in the NCBI Sequence Read Archive under the accession numbers SRR11285961 (RH3716) and SRR11285962 (LW2591Y) under Bioproject PRJNA611915.

Sanger Sequencing of PEX/ATG Deletion Strains

Genomic DNA of the PEX/ATG deletion strains was extracted according to the protocol described for SNP analysis and served as template for a PCR with primers JG1882/1883 and Phusion^TM DNA Polymerase (Thermo Fisher Scientific) to amplify the BUL1 gene. The amplified DNA was subjected to Sanger sequencing with primers JG1874 and JG1885 at Microsynth AG in Göttingen, Germany.

Semiquantitative Reverse Transcription-PCR

The strains were inoculated in 50 ml SC medium to an OD₆₀₀ = 0.1 from an overnight culture and grown at 30°C to an OD₆₀₀ ≈ 0.6. The cells were pelleted, washed, ground with liquid nitrogen and RNA was extracted with the RNeasy Plant Mini Kit (Qiagen). DNase digestion and subsequent cDNA synthesis with random hexamers was carried out using the QuantiTect Reverse Transcription Kit (Qiagen) on 0.8 mg of total RNA for each sample. Amplification was performed with the Phusion^TM High-Fidelity DNA Polymerase (Thermo Fisher Scientific) using 1 μl of cDNA and the following primers: JG1660/1233 for trHMG1, JG1659/1656 for IDI1, JG1301/1658 for MAF1, JG1304/1720 for mFPS^N144W, JG1515/1318 for GPPS, and JG1357/1320 for GES (Supplementary Table S4). Amplification conditions were as follows: 40 s at 98°C, 25 cycles of 15 s at 98°C, 30 s at 60°C, and 12 s at 72°C with an adjacent step at 72°C for 5 min. The amplicons were separated on a 1% agarose gel with the Thermo Scientific^TM GeneRuler^TM 1 kb DNA ladder SM0311 (Thermo Fisher Scientific) as size marker.

Peroxisome Quantification With Fluorescence Microscopy

The strains were grown for 1 or 2 days in liquid YPD or SC medium at 30°C, diluted approximately 1:10 with water and directly used for fluorescence microscopy. Per strain 76 (for YPD medium) and 60 (for SC medium) fluorescent images (38 and 20 per replicate, respectively) with 100× magnification were randomly taken using a Z1 microscope (Zeiss) equipped with a CSU-X1 A1 confocal scanner unit (Yokogawa), QuantEM:512SC digital camera (Photometrics) and SlideBook 6.0 software package (Intelligent Imaging Innovations). The number of cells (green fluorescence, GFP) and peroxisomes (red fluorescence, mCherry-SKL) were automatically counted from each image with the SlideBook 6.0 software package. The ratio of peroxisomes/cell was calculated from the number of cells and the number of peroxisomes per image, resulting in 76 and 60 values per strain, respectively. The values were depicted in a Tukey’s box-whisker-plot using R (R Development Core Team, 2011) and ggplot2 (Wickham, 2016).

The coding gene of GFP, which fluorescence was used for counting the cells, was integrated into the genome of the S. cerevisiae strains during the process of reiterative recombination (Supplementary Figure S6A; Wingler and Cornish, 2011). The fluorescence of mCherry, which was tagged with the peroxisomal target sequence 1, SKL, was used for counting the peroxisomes. In total, peroxisomes of a minimum of 2,131 cells per strain from two biological replicates for YPD medium and a minimum of 2,886 cells per strain from three biological replicates for SC medium were quantified. For statistical analysis and calculation of p-values, the two-tailed two-sample T-test was used.

Fluorescence Microscopy for Protein Localizations

The strains were grown overnight in liquid SC medium at 30°C, diluted approximately 1:10 with water and directly used for fluorescence microscopy. Images were taken with 100× magnification using a Z1 microscope (Zeiss) equipped with a CSU-X1 A1 confocal scanner unit (Yokogawa), QuantEM:512SC digital camera (Photometrics) and SlideBook 6.0 software package (Intelligent Imaging Innovations).

Western Blot Analysis

Forty microgram protein extract was separated by SDS-PAGE and transferred to a nitrocellulose membrane by electroblotting as described previously (Opitz et al., 2017). Total proteins were visualized with Ponceau Red (0.2% Ponceau S, 3% trichloroacetic acid). The membrane was blocked with Tris-buffered saline with Tween20 (TBS-T) with 5% skim milk powder for 1 h. α-RFP antibody (Abcam, ab62341), diluted 1:500 in blocking buffer, was used as primary and goat anti-rabbit IgG horseradish peroxidase conjugate (MoBiTec, G21234), diluted 1:1000 in blocking buffer, as secondary antibody.

Growth Tests

For growth tests on solid medium, the strains were grown to mid-log phase and diluted to an OD₆₀₀ of 0.1. After preparing serial 10-fold dilutions, the cells were spotted on SC medium, YPD medium, YPO medium or YPD containing geraniol, respectively. The plates were incubated for 3 days at 30°C. The experiments were repeated for at least two biological replicates (n = 2) with each three technical replicates. For growth tests in liquid medium, the strains were inoculated from an overnight culture to an OD₆₀₀ = 0.2 in 200 μl SC medium without histidine and uracil. Optical density measurements at 600 nm were performed in two biological (n = 2) and three technical replicates in 96-well plates for 40 h using a microplate reader with temperature control at 30°C (Infinite M200, Tecan) and the Magellan 7.2 software (Tecan). The growth rate μ and the maximum cell growth A were calculated using the grofit package for R (Kahm et al., 2010; R Development Core Team, 2011) and the spline fit. The two-tailed two-sample T-test was used to calculate the p-values.

Geraniol Extraction and Quantification With Gas Chromatography-Mass Spectrometry (GC-MS)

The strains were inoculated from an overnight liquid culture in 5 ml SC medium without histidine and uracil to an OD₆₀₀ = 0.1. The strains were cultivated for 2 days at 30°C with constant shaking in two biological replicates (n = 2). The OD₆₀₀ was measured and the culture was extracted with 1 ml dodecane. The dodecane phase was dried with MgSO₄. 200 μl of the dodecane phase were mixed with 4 μl 0.1% (w/v) menthol in dodecane as an internal standard and measured in a gas chromatography-mass spectrometry (GC-MS) system from Thermo Fisher Scientific, consisting of a TRACE GC Ultra gas chromatograph and a DSQ Quadrupole mass spectrometer with two technical replicates. The system was equipped with a TRI PLUS RSH autoinjector and a DB-WAX Ultra Inert column (30 m × 0.25 mm, 0.25 μm, Agilent Technologies). Samples were injected with a split ratio of 1:10 and helium as carrier gas at a flow rate of 1 ml/min. The PTV injector temperature was 250°C. The oven temperature gradient program started at 40°C for 1 min followed by heating to 230°C with a rate of 15°C/min and a hold for 3 min. The column effluent was introduced into the EI(+) ion source (200°C) of the MS. The mass scan range was m/z 50–300 (0.52 s). Geraniol was quantified using linear calibration curves from 5 × 10^–5% (w/v) to 0.01% (w/v) geraniol. Data are given in Supplementary Table S2.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/, SRR11285961 and https://www.ncbi.nlm.nih.gov/, SRR11285962.

Author Contributions

JG, ET, NS, and GB conceived the study. JG, HF, DS, TH, AP, ZZ, ET, and GB designed the experiments. JG, HF, MW, AP, and ZZ performed the experiments. JG, HF, DS, ZZ, and ET analyzed the data. JG, HF, DS, AP, NS, and GB wrote the manuscript. All authors read and approved the final manuscript.

Funding

The study was funded by the ERA-IB project “TERPENOSOME: Engineered compartments for monoterpenoid production using synthetic biology” (funding reference number 031A337A) and the Deutsche Forschungsgemeinschaft (DFG BR1502/19-1). ZZ, NS, and ET was funded by Biotechnology and Biological Sciences Research Council (BB/L027593/1) as part of the ERA-IB TERPENOSOME project. Funding for Open Access publication was supported by Göttingen University.

Conflict of Interest

TH was employed by the company Thermo Fisher Scientific.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Verena Große, Ute Neef, Verena Hofer-Preetz, and Mechthild Bömeke for their technical assistance, Maria Vinaixia for the data evaluation of the initial GC-MS experiments, and Aroma Chemical Services International GmbH for providing geraniol as reference. We also thank Virginia W. Cornish for providing LW2591Y and Helge Bode and Blagovesta Popova for providing plasmids.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2020.582052/full#supplementary-material

References

Akada, R., Kitagawa, T., Kaneko, S., Toyonaga, D., Ito, S., Kakihara, Y., et al. (2006). PCR-mediated seamless gene deletion and marker recycling in Saccharomyces cerevisiae. Yeast 23, 399–405. doi: 10.1002/yea.1365

PubMed Abstract | CrossRef Full Text | Google Scholar

Bard, M., Albrecht, M. R., Gupta, N., Guynn, C. J., and Stillwell, W. (1988). Geraniol interferes with membrane functions in strains of Candida and Saccharomyces. Lipids 23, 534–538. doi: 10.1007/bf02535593

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, S., Clastre, M., Courdavault, V., and O’Connor, S. E. (2015). De novo production of the plant-derived alkaloid strictosidine in yeast. Proc. Natl. Acad. Sci. U.S.A. 112, 3205–3210. doi: 10.1073/pnas.1423555112

PubMed Abstract | CrossRef Full Text | Google Scholar

Carnesecchi, S., Bradaia, A., Fischer, B., Coelho, D., Schöller-Guinard, M., Gosse, F., et al. (2002). Perturbation by geraniol of cell membrane permeability and signal transduction pathways in human colon cancer cells. J. Pharmacol. Exp. Ther. 303, 711–715. doi: 10.1124/jpet.102.039263

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, S., Zhou, Y., Chen, Y., and Gu, J. (2018). fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890. doi: 10.1093/bioinformatics/bty560

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., and Viljoen, A. M. (2010). Geraniol - a review of a commercially important fragrance material. South Afr. J. Bot. 76, 643–651. doi: 10.1016/j.sajb.2010.05.008

CrossRef Full Text | Google Scholar

Chen, X., Zhang, C., and Lindley, N. D. (2020). Metabolic engineering strategies for sustainable terpenoid flavor and fragrance synthesis. J. Agric. Food Chem. 23:acs.jafc.9b06203. doi: 10.1021/acs.jafc.9b06203

PubMed Abstract | CrossRef Full Text | Google Scholar

Deatherage, D. E., and Barrick, J. E. (2014). Identification of mutations in laboratory-evolved microbes from next-generation sequencing data using breseq. Methods Mol. Biol. 1151, 165–188. doi: 10.1007/978-1-4939-0554-6_12

CrossRef Full Text | Google Scholar

DeLoache, W. C., and Dueber, J. E. (2013). Compartmentalizing metabolic pathways in organelles. Nat. Biotechnol. 31, 320–321. doi: 10.1038/nbt.2549

PubMed Abstract | CrossRef Full Text | Google Scholar

DeLoache, W. C., Russ, Z. N., and Dueber, J. E. (2016). Towards repurposing the yeast peroxisome for compartmentalizing heterologous metabolic pathways. Nat. Commun. 7:11152. doi: 10.1038/ncomms11152

PubMed Abstract | CrossRef Full Text | Google Scholar

Donald, K. A., Hampton, R. Y., and Fritz, I. B. (1997). Effects of overproduction of the catalytic domain of 3-hydroxy-3-methylglutaryl coenzyme A reductase on squalene synthesis in Saccharomyces cerevisiae. Appl. Environ. Microbiol. 63, 3341–3344. doi: 10.1128/aem.63.9.3341-3344.1997

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehrenworth, A. M., Haines, M. A., Wong, A., and Peralta-Yahya, P. (2017). Quantifying the efficiency of Saccharomyces cerevisiae translocation tags. Biotechnol. Bioeng. 114, 2628–2636. doi: 10.1002/bit.26376

PubMed Abstract | CrossRef Full Text | Google Scholar

George, K. W., Alonso-Gutierrez, J., Keasling, J. D., and Lee, T. S. (2015). “Isoprenoid drugs, biofuels, and chemicals—artemisinin, farnesene, and beyond,” in Biotechnology of Isoprenoids. Advances in Biochemical Engineering/Biotechnology, Vol. 148, eds J. Schrader and J. Bohlmann (Cham: Springer), 355–389. doi: 10.1007/10_2014_288

CrossRef Full Text | Google Scholar

Gurvitz, A., and Rottensteiner, H. (2006). The biochemistry of oleate induction: transcriptional upregulation and peroxisome proliferation. Biochim. Biophys. Acta Mol. Cell Res. 1763, 1392–1402. doi: 10.1016/j.bbamcr.2006.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Herr, A., and Fischer, R. (2014). Improvement of Aspergillus nidulans penicillin production by targeting AcvA to peroxisomes. Metab. Eng. 25, 131–139. doi: 10.1016/J.YMBEN.2014.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, R. Y., Eddy, M., Vujcic, M., and Kowalski, D. (2005). Genome-wide screen identifies genes whose inactivation confer resistance to cisplatin in Saccharomyces cerevisiae. Cancer Res. 65, 5890–5897. doi: 10.1158/0008-5472.CAN-04-4093

PubMed Abstract | CrossRef Full Text | Google Scholar

Ignea, C., Cvetkovic, I., Loupassaki, S., Kefalas, P., Johnson, C. B., Kampranis, S. C., et al. (2011). Improving yeast strains using recyclable integration cassettes, for the production of plant terpenoids. Microb. Cell Fact 10:4. doi: 10.1186/1475-2859-10-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Imazaki, A., Tanaka, A., Harimoto, Y., Yamamoto, M., Akimitsu, K., Park, P., et al. (2010). Contribution of peroxisomes to secondary metabolism and pathogenicity in the fungal plant pathogen Alternaria alternata. Eukaryot. Cell 9, 682–694. doi: 10.1128/EC.00369-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Joshi, A. S., Huang, X., Choudhary, V., Levine, T. P., Hu, J., and Prinz, W. A. (2016). A family of membrane-shaping proteins at ER subdomains regulates pre-peroxisomal vesicle biogenesis. J. Cell Biol. 215, 515–529. doi: 10.1083/jcb.201602064

PubMed Abstract | CrossRef Full Text | Google Scholar

Kahm, M., Hasenbrink, G., Lichtenberg-Fraté, H., Ludwig, J., and Kschischo, M. (2010). grofit: fitting biological growth curves with R. J. Stat. Softw. 33, 1–21. doi: 10.18637/jss.v033.i07

CrossRef Full Text | Google Scholar

Kapitzky, L., Beltrao, P., Berens, T. J., Gassner, N., Zhou, C., Wüster, A., et al. (2010). Cross-species chemogenomic profiling reveals evolutionarily conserved drug mode of action. Mol. Syst. Biol. 6:451. doi: 10.1038/msb.2010.107

PubMed Abstract | CrossRef Full Text | Google Scholar

Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923

PubMed Abstract | CrossRef Full Text | Google Scholar

Lapczynski, A., Bhatia, S. P., Foxenberg, R. J., Letizia, C. S., and Api, A. M. (2008). Fragrance material review on geraniol. Food Chem. Toxicol. 46, S160–S170. doi: 10.1016/j.fct.2008.06.048

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, G. S., Li, T., Zhou, W., Jiang, M., Tao, X. Y., Liu, M., et al. (2020). The yeast peroxisome: a dynamic storage depot and subcellular factory for squalene overproduction. Metab. Eng. 57, 151–161. doi: 10.1016/j.ymben.2019.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Zhang, W., Du, G., Chen, J., and Zhou, J. (2013). Overproduction of geraniol by enhanced precursor supply in Saccharomyces cerevisiae. J. Biotechnol. 168, 446–451. doi: 10.1016/j.jbiotec.2013.10.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, W., Xu, X., Zhang, R., Cheng, T., Cao, Y., Li, X., et al. (2016). Engineering Escherichia coli for high-yield geraniol production with biotransformation of geranyl acetate to geraniol under fed-batch culture. Biotechnol. Biofuels 9:58. doi: 10.1186/s13068-016-0466-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10–12. doi: 10.14806/ej.17.1.200

CrossRef Full Text | Google Scholar

Maury, J., Asadollahi, M. A., Møller, K., Clark, A., and Nielsen, J. (2005). Microbial isoprenoid production: an example of green chemistry through metabolic engineering. Adv. Biochem. Eng. Biotechnol. 100, 19–51. doi: 10.1007/b136410

PubMed Abstract | CrossRef Full Text | Google Scholar

Motley, A. M., Nuttall, J. M., and Hettema, E. H. (2012a). Atg36: the Saccharomyces cerevisiae receptor for pexophagy. Autophagy 8, 1680–1681. doi: 10.4161/auto.21485

PubMed Abstract | CrossRef Full Text | Google Scholar

Motley, A. M., Nuttall, J. M., and Hettema, E. H. (2012b). Pex3-anchored Atg36 tags peroxisomes for degradation in Saccharomyces cerevisiae. EMBO J. 31, 2852–2868. doi: 10.1038/emboj.2012.151

PubMed Abstract | CrossRef Full Text | Google Scholar

Novoselova, T. V., Zahira, K., Rose, R. S., and Sullivan, J. A. (2012). Bul proteins, a nonredundant, antagonistic family of ubiquitin ligase regulatory proteins. Eukaryot. Cell 11, 463–470. doi: 10.1128/EC.00009-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Nuttall, J. M., Motley, A., and Hettema, E. H. (2011). Peroxisome biogenesis: recent advances. Curr. Opin. Cell Biol. 23, 421–426. doi: 10.1016/j.ceb.2011.05.005

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Donnell, A. F. (2012). The running of the Buls: control of permease trafficking by α-arrestins Bul1 and Bul2. Mol. Cell. Biol. 32, 4506–4509. doi: 10.1128/mcb.01176-12

PubMed Abstract | CrossRef Full Text | Google Scholar

Opitz, N., Schmitt, K., Hofer-Pretz, V., Neumann, B., Krebber, H., Braus, G. H., et al. (2017). Capturing the Asc1p/receptor for activated C Kinase 1 (RACK1) microenvironment at the head region of the 40s ribosome with quantitative BioID in yeast. Mol. Cell. Proteom. 16, 2199–2218. doi: 10.1074/mcp.M116.066654

PubMed Abstract | CrossRef Full Text | Google Scholar

Oswald, M., Fischer, M., Dirninger, N., and Karst, F. (2007). Monoterpenoid biosynthesis in Saccharomyces cerevisiae. FEMS Yeast Res. 7, 413–421. doi: 10.1111/j.1567-1364.2006.00172.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Poirier, Y., Erard, N., and Petétot, J. M. (2001). Synthesis of polyhydroxyalkanoate in the peroxisome of Saccharomyces cerevisiae by using intermediates of fatty acid beta-oxidation. Appl. Environ. Microbiol. 67, 5254–5260. doi: 10.1128/AEM.67.11.5254-5260.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

R Development Core Team (2011). R: A Language and Environment for Statistical Computing. Vienna: R Development Core Team.

Google Scholar

Rico, J., Pardo, E., and Orejas, M. (2010). Enhanced production of a plant monoterpene by overexpression of the 3-hydroxy-3-methylglutaryl coenzyme a reductase catalytic domain in saccharomyces cerevisiae. Appl. Environ. Microbiol. 76, 6449–6454. doi: 10.1128/AEM.02987-09

PubMed Abstract | CrossRef Full Text | Google Scholar

Shah, A. A., Wang, C., Chung, Y.-R., Kim, J.-Y., Choi, E.-S., and Kim, S.-W. (2013). Enhancement of geraniol resistance of Escherichia coli by MarA overexpression. J. Biosci. Bioeng. 115, 253–258. doi: 10.1016/J.JBIOSC.2012.10.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Sibirny, A. A. (2016). Yeast peroxisomes: structure, functions and biotechnological opportunities. FEMS Yeast Res. 16, 1–14. doi: 10.1093/femsyr/fow038

PubMed Abstract | CrossRef Full Text | Google Scholar

Stanley Fernandez, S. M., Kellogg, B. A., and Dale Poulter, C. (2000). Farnesyl diphosphate synthase. Altering the catalytic site to select for geranyl diphosphate activity. Biochemistry 39, 15316–15321. doi: 10.1021/bi0014305

PubMed Abstract | CrossRef Full Text | Google Scholar

Steyer, D., Erny, C., Claudel, P., Riveill, G., Karst, F., and Legras, J.-L. (2013). Genetic analysis of geraniol metabolism during fermentation. Food Microbiol. 33, 228–234. doi: 10.1016/j.fm.2012.09.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Surburg, H., and Panten, J. (2006). Common Fragrance and Flavor Materials, 5th Edn, Weinheim: Wiley.

Google Scholar

Szewczyk, E., Nayak, T., Oakley, C. E., Edgerton, H., Xiong, Y., Taheri-Talesh, N., et al. (2007). Fusion PCR and gene targeting in Aspergillus nidulans. Nat. Protoc. 1, 3111–3120. doi: 10.1038/nprot.2006.405

PubMed Abstract | CrossRef Full Text | Google Scholar

Till, A., Lakhani, R., Burnett, S. F., and Subramani, S. (2012). Pexophagy: the selective degradation of peroxisomes. Int. J. Cell Biol. 2012:512721. doi: 10.1155/2012/512721

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Der Klei, I. J., Harder, W., and Veenhuis, M. (1991). Biosynthesis and assembly of alcohol oxidase, a peroxisomal matrix protein in methylotrophic yeasts: a review. Yeast 7, 195–209. doi: 10.1002/yea.320070302

PubMed Abstract | CrossRef Full Text | Google Scholar

Veenhuis, M., Van Dijken, J. P., Pilon, S. A. F., and Harder, W. (1978). Development of crystalline peroxisomes in methanol-grown cells of the yeast Hansenula polymorpha and its relation to environmental conditions. Arch. Microbiol. 117, 153–163. doi: 10.1007/BF00402303

PubMed Abstract | CrossRef Full Text | Google Scholar

Vizeacoumar, F. J., Torres-Guzman, J. C., Bouard, D., Aitchison, J. D., and Rachubinski, R. A. (2004). Pex30p, Pex31p, and Pex32p form a family of peroxisomal integral membrane proteins regulating peroxisome size and number in Saccharomyces cerevisiae. Mol. Biol. Cell 15, 665–667. doi: 10.1091/mbc.E03

CrossRef Full Text | Google Scholar

Wickham, H. (2016). ggplot2: Elegant Graphics for Data Analysis. New York, NY: Springer-Verlag.

Google Scholar

Wingler, L. M., and Cornish, V. W. (2011). Reiterative recombination for the in vivo assembly of libraries of multigene pathways. Proc. Natl. Acad. Sci. U.S.A. 108, 15135–15140. doi: 10.1073/pnas.1100507108

PubMed Abstract | CrossRef Full Text | Google Scholar

Wolfe, D., Reiner, T., Keeley, J. L., Pizzini, M., and Keil, R. L. (1999). Ubiquitin metabolism affects cellular response to volatile anesthetics in yeast. Mol. Cell. Biol. 19, 8254–8262. doi: 10.1128/mcb.19.12.8254

PubMed Abstract | CrossRef Full Text | Google Scholar

Woodcock, D. M., Crowther, P. J., Doherty, J., Jefferson, S., Decruz, E., Noyer-Weidner, M., et al. (1989). Quantitative evaluation of Escherichia coli host strains for tolerance to cytosine methylation in plasmid and phage recombinants. Nucleic Acids Res. 17, 3469–3478. doi: 10.1093/nar/17.9.3469

PubMed Abstract | CrossRef Full Text | Google Scholar

Zebec, Z., Wilkes, J., Jervis, A. J., Scrutton, N. S., Takano, E., and Breitling, R. (2016). Towards synthesis of monoterpenes and derivatives using synthetic biology. Curr. Opin. Chem. Biol. 34, 37–43. doi: 10.1016/J.CBPA.2016.06.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Bao, X., Li, C., Shen, Y., and Hou, J. (2016). Improving monoterpene geraniol production through geranyl diphosphate synthesis regulation in Saccharomyces cerevisiae. Appl. Microbiol. Biotechnol. 100, 4561–4571. doi: 10.1007/s00253-016-7375-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, J., Li, C., Zhang, Y., Shen, Y., Hou, J., and Bao, X. (2017). Dynamic control of ERG20 expression combined with minimized endogenous downstream metabolism contributes to the improvement of geraniol production in Saccharomyces cerevisiae. Microb. Cell Fact 16:17. doi: 10.1186/s12934-017-0641-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y. J., Buijs, N. A., Zhu, Z., Gómez, D. O., Boonsombuti, A., Siewers, V., et al. (2016). Harnessing yeast peroxisomes for biosynthesis of fatty-acid-derived biofuels and chemicals with relieved side-pathway competition. J. Am. Chem. Soc. 138, 15368–15377. doi: 10.1021/jacs.6b07394

PubMed Abstract | CrossRef Full Text | Google Scholar