High production of valencene in Saccharomyces cerevisiae through metabolic engineering
Microbial Cell Factories volume 18, Article number: 195 (2019)
The biological synthesis of high value compounds in industry through metabolically engineered microorganism factories has received increasing attention in recent years. Valencene is a high value ingredient in the flavor and fragrance industry, but the low concentration in nature and high cost of extraction limits its application. Saccharomyces cerevisiae, generally recognized as safe, is one of the most commonly used gene expression hosts. Construction of S. cerevisiae cell factory to achieve high production of valencene will be attractive.
Valencene was successfully biosynthesized after introducing valencene synthase into S. cerevisiae BJ5464. A significant increase in valencene yield was observed after down-regulation or knock-out of squalene synthesis and other inhibiting factors (such as erg9, rox1) in mevalonate (MVA) pathway using a recyclable CRISPR/Cas9 system constructed in this study through the introduction of Cre/loxP. To increase the supplement of the precursor farnesyl pyrophosphate (FPP), all the genes of FPP upstream in MVA pathway were overexpressed in yeast genome. Furthermore, valencene expression cassettes containing different promoters and terminators were compared, and PHXT7-VS-TTPI1 was found to have excellent performance in valencene production. Finally, after fed-batch fermentation in 3 L bioreactor, valencene production titer reached 539.3 mg/L with about 160-fold improvement compared to the initial titer, which is the highest reported valencene yield.
This study achieved high production of valencene in S. cerevisiae through metabolic engineering and optimization of expression cassette, providing good example of microbial overproduction of valuable chemical products. The construction of recyclable plasmid was useful for multiple gene editing as well.
Valencene, a natural sesquiterpene, possesses good biological activity and is found in various citrus species, such as in the essential oil of Valencia orange . When used in flavoring and fragrances, valencene imparts a woody citrus characteristic and can be economically used as an additive in food and drinks . Besides, valencene can be further derivatized into many economically useful sesquiterpenes, such as its oxidation product nootkatone which can inhibit proliferation of cancer cells. Valencene has broad application prospects in food, cosmetics and pharmaceutical industries . However, it is not economically feasible to extract valencene in large-scale from natural sources due to its low concentration in citrus fruits and high cost of purification . Constructing microbial cell factories to carry out the production of valencene is more attractive on account of the short growth cycle and low cost.
Previous studies have achieved high valencene yield through expression of valencene synthase gene (CnVS) in genetically engineered Rhodobacter sphaeroides strains, including expression of mevalonate operon from Paracoccus zeaxanthinifaciens . Frohwitter et al. constructed recombinant Corynebacterium glutamicum containing valencene synthase from Nootka cypress and genes erg20 and ispA were overexpressed, resulting in valencene yield of 2.41 mg/L . In previous studies, Saccharomyces cerevisiae has been considered an ideal host for metabolic engineering of valencene production due to its favorable physiological properties, and optimized fermentations produced maximal levels of 20 mg/L valencene [4, 7, 8]. However, the yield of valencene in S. cerevisiae is still unable to meet the industrial demand, which indicates the need for achieving higher yields through metabolic engineering.
All terpenes are biosynthesized through mevalonate pathway (MVA) in yeast, originating from acetyl-CoA, and the intermediate product farnesyl diphosphate (FPP) is the direct precursor of valencene. FPP is a key metabolic point in MVA pathway toward the metabolism branches, such as the synthesis of squalene with squalene synthase encoded by erg9  and synthesis of GGPP with geranylgeranyl diphosphate synthase encoded by bts1 , leading to synthesis of isoprenes . It is reported that rox1 is a transcriptional factor inhibiting the expression of hypoxia-induced genes in MVA pathway and ergosterol biosynthesis [12, 13]. Moreover, some distantly located genetic loci might have potential interactions with some terpenoid pathway, such as ypl062w and yjl064w [7, 14]. Recent research pointed out that the ypl062w functioned as an important promoter for ald6 whose expression level was negatively correlated with terpenoid productivity . Apart from pathway engineering, the selection of appropriate expression cassettes including promoter and terminator is also important for the optimization of valencene production . This is because the transcription level is determined by the promoter strength, and also because promoters behave differently under different growth conditions. Stephanopoulos and coworkers systematically tested promoters of different strengths and plasmid copy-numbers, and were able to maximize taxadiene production to approximately 1 g/L with minimal accumulation of any toxic intermediate [17,18,19]. Thus, it is attractive to test the yield of valencene using different expression cassettes.
In this study, S. cerevisiae was selected as the host to construct cell factory and achieve the overproduction of valencene. A series of metabolic engineering strategies were performed in MVA pathway (Fig. 1 ). Firstly, a recyclable CRISPR/Cas9 system was constructed to achieve multiple genome editing to accumulate the FPP pool toward the synthesis of valencene. Secondly, to enhance the metabolic flux to precursor FPP, the key genes of FPP synthetic pathway in MVA pathway were overexpressed in yeast genome. Thirdly, valencene synthase expression cassette was optimized through constructing a promoter-terminator library. Using the above combinatorial engineering strategies, this work offers a good reference to increase the heterologous expression performance of valuable compounds in S. cerevisiae through metabolic engineering.
Results and discussion
Construction of recyclable plasmid and application in CRISPR/Cas9 system
CRISPR/Cas9 system, the essential components of which are gRNA expression plasmid and Cas9 protein, is the most widely used genome editing technology for prokaryotes and eukaryotes due to its efficient genome editing ability [20,21,22]. However, when using CRISPR/Cas9 system to perform multiple genetic modifications, such as gene knock-out, knock-in or tagging, there are bottlenecks because of the limited selection markers on each gRNA expression plasmid. To solve this problem, some marker recycling systems are used currently, such as FLP/FRT system [23, 24], Cre/loxP system , counter-selection system [26, 27] and I-SceI .
In this study, a recyclable plasmid mediated by Cre/loxP system was constructed. First, under the control of galactose-inducible promoter GAL1, the Cre recombinase gene was introduced into plasmid p426, and subsequently, two loxP direct repeats were placed flanking 2μ replication origin of plasmid p426, resulting in the recyclable plasmid P426-CL (Additional file 1: Fig S1). When the strain containing P426-CL was cultivated in galactose medium, expression of Cre recombinase was induced. This mediated the excision of the 2μ replication origin on P426-CL, leading to the lack of replication capacity. The plasmid P426-CL was lost within cells in the passage and the selection marker of URA could be repeatedly used, which meant genome editing could proceed again in strains that had been edited by recyclable gRNA expression plasmid. As the existence of selection marker URA is important for the efficiency of subsequent gene editing, it is necessary to verify the plasmid loss efficiency. The efficiency results are shown in Table 1. It can be seen that, after 36 h culturing in galactose medium, the ratio of plasmid-free cells was observed to be nearly 100% (Table 1), which was higher than that obtained using Cre/loxP in yeast Pichia (65%) . Thus, the engineered yeast strain could be applied to a new round of genome editing using the gRNA vector (based on P426-CL) targeting another gene of interest (knock-in or knock-out). Ryan et al.  reported a scarless and marker-free genome editing method based on CRISPR-Cas9 in S. cerevisiae. However, the plasmid that coexpresses the Cas9 endonuclease and guide RNA (sgRNA) expression cassette will still be remained in the host cell, and the selection marker on the plasmid could not be recycled. In the study of Stovicek et al. , researchers focused on the development of the genetic manipulation on unrelated prototrophic polyploidy S. cerevisiae strains with the efficiency ranging between 65% and 78%. Although it is claimed that marker-free recombinants were isolated through natural loss of plasmid, the loss efficiency was not stable and satisfactory. The high loss rate of P426-CL under galactose induction gave us a guarantee on the recyclable CRISPR/Cas9 system.
Based on recyclable plasmid P426-CL, five gRNA expression plasmids (Additional file 1: Table S1) were constructed to target different gene loci (bts1, rox1, erg9, ypl062w and yjl064w) (Additional file 1: Table S2) that may influence valencene production in FPP upstream. Then, 31 mutant combinations were obtained through repeating use of these five gRNA expression plasmids and recycling of selection marker (Fig. 2). To evaluate the efficiency of mutant combinations, 10 colonies of each genome editing strain were selected for genomic PCR and sequencing. The results (Additional file 1: Fig S2) showed that the mutant efficiency of almost all mutant combinations was 100%, which showed a higher genome editing success rate than that reported by Jakočiūnas et al. , especially in multiple sites editing. Only five single gRNA expression plasmids were required to achieve 31 mutant combinations instead of constructing many gRNA expression plasmids for each mutant, which may save a lot of work . More practically, with the help of recyclable CRISPR/Cas9 system, new genetic modifications can be made directly based on the optimized strains instead of repeating the previous strategies first.
Regulation of FPP downstream for valencene production
After successful multiple genome editing of five targeted gene loci that would influence valencene synthesis, the valencene yields of wild-type strain and 31 mutant strains (Additional file 1: Table S3) were examined through fermentation in shake flasks. Valencene synthase from Callitropsis nootkatensis (CnVS) was introduced into the recombination strain. The detection of valencene was performed by GC-FID and GC–MS (Fig. 3). From the fermentation results (Fig. 4), it was found that all mutant strains obtained higher yields of valencene compared to wild type and the highest strain gave 5.7 mg/g DCW, which was a 3.2-fold increase. As shown in previous study , the single knock-out of rox1 gene improved valencene yield more than other single gene mutant. Moreover, all the 16 mutant strains with knock-out of rox1 showed significant enhancement on valencene yield. Besides, erg9 repression through the deletion of 45-bp (coordinates are − 220 to − 175) in the range of upstream activating cis-element (UAS) in erg9 promoter resulted in obvious increase in yield of valencene in both single and multiple mutants. For the ergosterol synthesis charged by erg9 is necessary to S. cerevisiae survival, researchers often tried to downregulate the erg9 expression instead of blocking the pathway thoroughly. It has been demonstrated that the alteration of UAS could obviously reduce the ERG9 expression . When the deletion of rox1 and alteration of UAS in erg9 promoter were combined, the increase of valencene production was more significant, indicating the effect of two mutations were synergistic.
Among these 31 strains, higher multiple genome editing generally showed higher increase in yield of valencene than lower multiple genome editing. For instance, the quadruple mutant rox1;erg9;ypl062w;yjl064w (2.8-fold) showed higher increase in yield compared to triple mutant erg9;ypl062w;yjl064w (twofold). However, some combinatorial mutants did not show this pattern. For example, the quadruple mutant bts1;rox1;erg9;yjl064w (1.9-fold) showed much lower increase in yield than triple mutant rox1;erg9;yjl064w (2.7-fold), and bts1;rox1;ypl062w;yjl064w (1.8-fold) showed lower increase than triple mutant rox1;erg9;ypl062w (2.9-fold). This showed that the interaction between the five genes was not always positive under different genetic backgrounds.
The other three single mutants, bts1, ypl062w, and yjl064w, did not show significant increase in yield of valencene, even though knock-out of bts1 was supposed to restrict metabolic flow to isoprenoid and promote valencene synthesis. The squalene synthase encoded by ERG9 has a Km value of 2.5 μM for FPP and a kcat number of 0.53/s , while geranylgeranyl diphosphate (GGPP) synthase encoded by BTS1 has a Km value of 3.2 μM for FPP and a kcat number of 0.025/s . The higher Km value for FPP and low turnover number of GGPP synthase indicate its low capacity. This may be the reason why the knock-out of bts1 was less effective than down-regulation of erg9.
These findings suggested that the gene alteration of erg9, rox1 and bts1 was beneficial to valencene production. While, knock-out of ypl062w and yjl064w showed weak effect on valencene synthesis in S. cerevisiae.
Overexpression of FPP synthetic pathway in MVA pathway
FPP is derived from the MVA pathway as precursor of valencene synthesis. However, the native FPP synthetic pathway has a low synthetic capacity that limits the production of valencene . The overexpression of endogenous truncated 3-hydroxy-3-methylglutaryl coenzyme A reductase (tHMG1), a key enzyme in MVA pathway, has been shown to increase the biosynthesis of FPP . However, few studies focus on the overexpression of other key genes in FPP synthetic pathway.
To enhance FPP synthesis, every gene in FPP synthetic pathway (erg10, erg13, tHMG1, erg12, erg8, erg19, erg20 and idi1) was integrated into the genome of BJM-03, in which the ERG9 expression was down-regulated. Aforementioned results indicated that alteration in ypl062w and yjl064w showed no obvious benefit on valencene production. While the dpp1 and lpp1 were reported to be responsible for the hydrolysis of FPP and their knock-out could enhance FPP supplementation . Based on this information, the integration sites for were finally selected as bts1, rox1, dpp1 and lpp1. Through the selection marker recyclable plasmid P426-CL, BJM-33 were obtained after four rounds of gene knock-out/integration (Fig. 5a).
The recombinant strains were cultivated in shake flask and valencene yield was determined by GC-FID. The result showed that strain BJM-33 increased valencene production by 5.2-fold to 9.4 mg/g DCW compared with wild-type strain BJM-00 (Fig. 5b). However, compared to strain BJM-15 which possessed the best valencene production in 31 mutant strains (Fig. 4), the yield only increased by 1.7-fold, which was not in accordance with expectation. In a previous study, geraniol production was increased by 1.4-fold through overexpression of only tHMG1 . In other works, augmentation of the native yeast genes in MVA pathway increased nerolidol production by 3.5-fold . The capacity of valencene synthase was suspected to be the limiting factor of increasing valencene production rather than the precursor FPP pool. The amount of FPP may have exceeded the conversion capacity of valencene synthase. So increasing valencene synthase expression may be necessary for high valencene yield.
Screening of promoter-terminator pairs for valencene production
The above work focused on the regulation of MVA pathway, which showed some increase in yield of valencene. However, choosing a suitable CnVS expression cassette containing different promoters and terminators would influence the production of target compounds as well. In this study, a promoter-terminator library consisting of 7 combinations was constructed, and CnVS gene was introduced using restriction enzymes (labeled as PADH1-VS-TADH1, PCYC1-VS-TCYC1, PFBA1-VS-TCYC1, PHXT7-VS-TTPI1, PPGK1-VS-TENO2, PTDH3-VS-TPGI1, PTEF1-VS-TFBA1). The recombination strains containing different CnVS expression cassettes were cultivated in shake flask to evaluate the change in valencene yield compared to strain BJM-00 (PTDH3-VS-TADH1). According to the results shown in Fig. 6a, the strengths of the 8 promoters could be ranked in the following order: PHXT7 > PFBA1 > PTEF1 > PCYC1 > PPGK1 > PTDH3 > PADH1.
The above result is contrary to the result reported by Partow et al. that PTEF1 showed strongest activity . Moreover, Jie Sun et al.  compared 14 constitutive promoters from S. cerevisiae and found that PTEF1 was the strongest and PPGI1 was the weakest. The inconsistent result implied that the activity of a promoter may not be the same in different expression cassettes or in different product’s heterologous expression. Interestingly, PTDH3-VS-TADH1 and PTDH3-VS-TPGI1, consisting of the same promoter but different terminators, also showed different valencene yields (1.8 mg/g DCW and 1.7 mg/g DCW), which indicated that both promoter and terminator could influence valencene synthase transcription level. As for promoter HXT7, Partow et al.  found that the activity of PHXT7 was low when glucose concentration was high and it increased along with consumption of glucose. Previous work also showed that PHXT7 was a high-affinity hexose transporter which was highly expressed at low glucose concentration (< 4.4 mM) [41, 42].
In order to combine the advantages of FPP supply improvement and expression cassette optimization, the best performing plasmid YM10 containing PHXT7-VS-TTPI1, was transformed into the strain BJM-42 and BJM-43, generating BJM-44 and BJM-45, respectively. After 48 h cultivation, valencene production was dramatically increased as 22.7 mg/g DCW in BJM-45, which is increased by 2.7 fold compared with BJM-44 (Fig. 6b). The result showed that overexpression of FPP upstream in MVA pathway had a large influence on valencene yield with an optimized expression cassette. The high performance of BJM-45 also indicated the effectiveness of combinatorial engineering strategies in valencene synthesis, which might be meaningful for biotechnological production of other valuable chemicals. The increase in yield caused by optimization of CnVS expression cassette also implied that the activity or expression of CnVS may be the main limiting factor. Modification of the enzyme and searching for new valencene synthase would be useful strategies for further improvement.
The fed-batch fermentation of valencene overproduction
To evaluate the production performance of the engineered strain BJM-45, fed-batch fermentation was performed with SD-glucose/ΔLeu medium in shake flask. After fermentation for 192 h with 6 times supplement of glucose, the highest valencene yield of 89.7 mg/L was obtained (Fig. 7a). The result showed that the strain basically stopped growing after 48 h and the residual sugar concentration in the medium kept increasing. The low dissolved oxygen in shake flask may cause this problem and restrict the high production of valencene.
Fed-batch fermentation was performed at a 1 L scale medium in 3 L bioreactor to determine the best performance of optimized strain BJM-45. After 136 h with 4 times supplement of glucose, the highest valencene titer of 539.3 mg/L was obtained, which was sixfold higher than that of fed-batch fermentation in shake flask (Fig. 7b). It was also much higher than the yields reported previously in engineered S. cerevisiae strains [4, 43] or in other valencene producers, such as Corynebacterium glutamicum  and Schizophyllum commune  (Table 2).
The course profiles in bioreactor (Fig. 7b) showed that glucose was almost entirely consumed and the biomass increased by twofold compared to shake flask. However, the cell density was still relatively low in fed-batch experiment and far below high-density fermentation (30 g/L DCW), as the highest CDW only reached 6 g/L throughout the fermentation. Additionally, the yeast cells at the later period of fermentation were collected for microscopic observation and it was found that most of the cells were dead (Additional file 1: Fig S3). In future, increasing the cell density in fed-batch fermentation and improving the cytoactivity of yeast cells at the later period of fermentation would be useful measures. Intracellular and extracellular metabolomics analysis would be an efficient way to find suitable biomarkers for batch media or for optimization of feed solution. It is expected that valencene production by the engineered strain in this study could be further improved by continuous efforts in both metabolic engineering and fermentation optimization.
In this study, the high production of high value compound valencene was achieved through metabolic engineering using modified genome editing technology. Valencene synthase from Callitropsis nootkatensis (CnVS) was introduced into yeast, achieving valencene biosynthesis. A recyclable plasmid mediated by Cre/loxP system was constructed and applied in CRISPR/Cas9 system to perform multiple genetic modifications. The FPP pool toward valencene synthesis was enhanced through knock-out or down-regulate of gene loci of pathway branch or inhibiting factor. The valencene yield was increased to some extent through these strategies, but the metabolic flux in MVA pathway still needs further study and precise regulation. The effect of overexpression of FPP synthetic pathway into yeast genome was estimated through shake flask fermentation., and some increase was obtained. Moreover, it was found that the expression cassette PHXT7-VS-TTPI1 showed significant influence on valencene yield. To achieve further optimization of expression cassette, more regulation work is required in the future. The best valencene titer (539.3 mg/L) were achieved in fed-batch of 3 L bioreactors, which was highest production reported so far. Specific valencene yield was improved by approximately 160-fold compared to the initial strain. This study provides a good reference for microbial overproduction of valuable chemical products through combinatorial pathway engineering and optimization of expression cassette.
Materials and methods
Chemicals and media
All the heterologous genes and primers used in this study were obtained from Sangon (Shanghai, China). PrimerSTAR Max Premix used for fragment cloning was purchased from TaKaRa (Dalian, China). DNA purification and plasmid isolation kits were purchased from Tiangen (Beijing, China). T4 DNA ligase and restriction endonuclease were purchased from Thermo fermentas (MD, USA). The S.c EasyComp Transformation kit used for yeast transformation was provided by Invitrogen (CA, USA). Standards of valencene were obtained from Sigma (Sigma-Aldrich, MO, USA). KOD FX, High Success-rate DNA polymerase, was used in colony PCR to screen positive clones and was purchased from TOYOBO (Japan).
Escherichia coli DH5α were grown in Luria–Bertani (LB) medium containing 10 g/L tryptone, 10 g/L NaCl and 5 g/L yeast extract. Agar was added into the medium with concentration of 20 g/L if solid plate was needed. The yeast strains were grown in Synthetic Dropout (SD) medium containing 6.2 g/L Do Supplement (Clontech, USA), 67 g/L Yeast Nitrogen Base without amino acids (YNB) with 20 g/L glucose or 20 g/L galactose as carbon source (denoted as SD-glucose and SD-galactose, respectively), minus the auxotrophic amino acids complemented by recombinant plasmids. SD-glucose medium supplemented with 1 g/L 5-fluoroorotic acid (5-FOA) was used for selecting yeast colonies which showed loss of Orotidine 5′-phosphate decarboxylase (ODCase, encoded by URA3) activity.
Construction of plasmids and strains
Construction of recyclable plasmid
The cre gene was PCR amplified with primer pair Cre-F and Cre-R (Additional file 1: Table S4) from plasmid pUC57 synthesized by Sangon (Shanghai). The Cre recombinase expression cassette was constructed under the control of the inducible Gal1 promoter in the plasmid YEplac181 between SmaI site and BamHI site. The expression fragment PGAL1-Cre-TCYC1 was amplified by PCR using primer pair of GAL1-F and CYC1-R and cloned into plasmid p426-PSNR52-gRNA.CAN1.Y-TSUP4 (shortened as P426, Addgene (USA)) using ClonExpress II One Step Cloning Kit (Vazyme, China), resulting in the plasmid p426-Cre. To introduce the loxP direct repeats into both sides of 2μ replication origin of plasmid, the p426-Cre was linearized using primer pair PB-F and Ori-R. Here, both primers contained part of loxP direct repeat sequence as homologous arms. Then, the linearized fragment was self-ligated using ClonExpress II One Step Cloning Kit and one loxP direct repeat was introduced into 3′ side of 2μ replication origin of the plasmid. Another loxP direct repeat sequence was inserted using primer pair of Ori-F and PB-R in the same way, and plasmid P426-CL was obtained.
Valencene expression plasmid construction
The valencene synthase gene from Callitropsis nootkatensis (CnVS) was synthesized by Sangon (Shanghai) with codon optimization and inserted into plasmid YEplac181 (LEU2 selection marker) under the control of TDH3 promoter between SmaI and BamHI restriction sites, yielding plasmid pYM06. The plasmid was transformed into yeast strain BJ5464 using standard protocols.
Knock-out and knock-down strains construction
To improve the valencene production, 5 gene loci that could influence metabolic flux toward valencene synthesis were chosen for metabolic engineering. Among these 5 gene loci, rox1 [12, 13], bts1 [43, 46], ypl062w and yjl064w  were knocked out, and erg9 [48, 49] was down-regulated. Then, 31 combinations of these 5 gene loci from single mutants to quintuple mutants were designed.
To obtain these 31 mutant constructions, 5 gRNA expression plasmids, P426-bts1, P426-rox1, P426-erg9, P426-ypl062w, and P426-yjl064w were constructed based on recyclable plasmid P426-CL. Taking P426-bts1 as an example, the bts1 gene sequence was acquired from NCBI, and CRISPy tool [11, 50] (http://staff.biosustain.dtu.dk/laeb/crispy_cenpk/) was used to choose the specific gRNA sequence (Additional file 1: Table S5) targeting bts1 gene in yeast BJ5464 genome. Subsequently, two primers bts1-F and bts1-R were designed, both containing gRNA sequence as homologous arms. Two fragments named bts1-1 and bts1-2 were PCR amplified with primer pairs tong-F/bts1-R and bts1-F/tong-R, respectively using P426-CL as template. Two fragments were assembled together using ClonExpress II One Step Cloning Kit and plasmid P426-bts1 containing bts1 gRNA sequence was generated. The other four gRNA expression plasmids were obtained in the same way. These gRNA expression plasmids were transformed into yeast strain containing Cas9 expression plasmid (Addgene, USA) with 90-bp donor DNA (Additional file 1: Table S6) to construct single mutant strains.
Regarding the multiple mutants construction, triple mutation bts1;rox1;erg9 was shown as an example. First, single mutant strain BJM-01 was cultivated in SD-galactose medium for 36 h to induce the expression of Cre recombinase. Then, yeast solution was plated in SD-glucose/5-FOA solid plate to select the positive colony that lost gRNA expression plasmid P426-bts1 containing URA3 selection marker. After recycling of selection marker, another gRNA expression plasmid P426-rox1 and appropriate donor DNA were transformed into the above positive strain, resulting in double mutant strain BJM-14. Similarly, the URA3 selection marker in BJM-14 was recycled through galactose induction and gRNA expression plasmid P426-erg9 was introduced into the strain, resulting in triple mutation strain BJM-24. Other multiple mutants were achieved through the reuse of 5 gRNA expression plasmids as described above (Fig. 2).
Overexpression of FPP biosynthetic pathway genes
To further enhance the supply of FPP, the FPP biosynthesis pathway genes (erg10, erg13, thmg1, erg12, erg8, erg19, idi1, erg20) were overexpressed in yeast genome through the plasmid P426-CL constructed in this study (Fig. 5-a). All the genes were amplified using S. cerevisiae genomic DNA as template. Four integration sites of bts1, rox1, dpp1 and lpp1 were chosen for the aforementioned 8 genes overexpression. The ypl062w and yjl064w were replaced by dpp1 and lpp1, because in this study, the alterations on them did not show positive effects on valencene production as other researchers’ report [7, 14]. While knock-out of dpp1 and lpp1 could reduce the FPP hydrolysis and increase the substrate supplementation for valencene synthesis .
A bidirectional expression vector containing PTDH3-TADH1 and PTEF1-TCYC1 was first constructed based on the plasmid YEplac181, denoted as pYM01. Genes of erg20 and idi1 were cloned from yeast genome. Two fragments named backone-1 and backone-2 were amplified with primer pairs cas1-F/181-bone-R and 181-bone-2-F/181-bone-2-R using pYM01 as template. These four fragments were assembled into a recombinant plasmid named pYM02 using ClonExpress MultiS One Step Cloning Kit (Vazyme, China) (Fig. 5a). The bidirectional expression fragment, TADH1-ERG20-PTDH3-PTEF1-IDI1-TCYC1, amplified from plasmid pYM02 was transformed into strain BJM-03 as a donor DNA with gRNA expression plasmid P426-rox1. After knock-out of rox1 gene loci and knock-in of fragment TADH1-ERG20-PTDH3-PTEF1-IDI1-TCYC1, the overexpression of erg20 gene and idi1 gene was achieved. Similar method was used to obtain overexpression of other remaining genes, resulting in strain BJM-33.
Construction of different CnVS expression cassettes
To optimize the expression cassette of valencene synthase (CnVS), a promoter-terminator library was constructed. A series of promoters (ADH1, CYC1, HXT7, FBA1, PGK1, TDH3, and TEF1) and terminators (ADH1, CYC1, TPI1, PGI1, and FBA1) from S. cerevisiae genome were chosen. The fragments of promoters and terminators were amplified by PCR with primer pairs containing restriction sites. The resulting fragments were assembled into plasmid YEplac181 after digestion by different restriction enzymes, generating the promoter-terminator library. Then, valencene synthase gene was inserted between promoter and terminator using XbaI and PstI to obtain different expression cassettes. All recombinant plasmids were transformed into wild-type strain BJ5464 to examine the valencene production.
Plasmid loss efficiency determination
The recyclable plasmid P426-CL containing the selection marker URA3 was transformed into S. cerevisiae BJ5464. Single colony of recombinant yeast was selected and inoculated into SD-galactose medium to induce Cre recombinase expression. After incubation at 30 °C for 36 h with shaking at 200 rpm, colony-forming unit per milliliter (CFU/mL) was determined using haemocytometer through microscopic counts. Then, 0.2 mL of yeast culture with different dilution factors was plated onto SD-glucose agar plate without uracil. Cells harboring the plasmid P426-CL grew out and colony number were recorded.
The plasmid loss efficiency was calculated using the following formula:
A1, cell density of yeast culture, CFU/mL, A2, number of colonies per plate, D, dilution factor.
Triplicate measurements were performed for each set.
Cultivation in shake flasks
The yields of valencene-producing strains constructed in this study were evaluated in shake flask through two-phase flask cultivation. SD-glucose medium was used in the entire cultivation process. Strains were recovered from glycerol stocks by streaking on SD-glucose agar plates and pre-cultured in SD-glucose medium to exponential phase (cell density OD600 between 1 and 3). Two-phase flask cultivation was initiated by inoculating pre-cultured seed broth to OD600 = 0.05 with 10 mL SD-glucose medium in 50 mL flasks. After adding 2 mL n-dodecane to extract valencene, flasks were sealed with breathable microporous film. Flask cultivation was performed at 30 °C and 200 rpm.
Fed-batch fermentation of valencene in bioreactors
Strain BJM-45 was selected for fed-batch fermentation. Seed cultures were prepared by inoculating 250 µL of glycerol-stock into a test tube containing 5 mL SD-glucose and culturing at 30 °C for 24 h to an OD600 of 3–5. Then, 3–4 mL of precultures were inoculated into a 500 mL shake-flask containing 100 mL SD-glucose/and subcultured for an additional 24 h at 30 °C to an OD600 of 5–6. Seed cultures were transferred into a 3 L bioreactor containing 1 L SD-glucose batch medium at a 5–10% (v/v) inoculum. Valencene was extracted using 200 mL n-dodecane. Fermentation was carried out at 30 °C with an air flow rate of 0.8 vvm. The dissolved oxygen level was kept at 60% by adjusting the agitation speed from 250 to 500 rpm and pH was controlled at 6.0 by automatic addition of 1 M NaOH. After 36 h of fermentation when the initial glucose was fully consumed, concentrated glucose solution (400 g/L) was fed into bioreactors every 24 h to maintain the growth of strain. The fermentation was ended at 135 h when the cells were basically inactive (Additional file 1: Fig S3).
Cell density in the medium was determined by measuring absorbance at 600 nm using an ultraviolet–visible spectrophotometer (Thermo Fisher Scientific GENESYS 10, Germany). The dry cell weight (DCW) was measured by analytical balance (Sartorius, Beijing) after the yeast solution was dried in a drying oven.
Extraction and analysis of valencene
The culture of valencene fermentation was covered with 20% n-dodecane to extract valencene secreted in the medium. To determine content of valencene, 500 μL sample collected from n-dodecane organic phase in the two-phase fermentation medium was mixed with 500 μL ethyl acetate and 2 μL isolongifolene. The resulting mixture was passed through a 0.22 μm sterile filter. Valencene concentration was detected by gas chromatography (GC, Agilent Technologies, 7890) equipped with a flame ionization detector (FID). An HP-5 column (30 m × 0.32 mm; 0.25 μm film thickness; Agilent Technologies) was used, with nitrogen as the carrier gas. Samples were run using the following analysis method: GC oven temperature program was maintained at 100 °C for 10 min, followed by a ramp to 200 °C at 10 °C/min, and then maintained at 200 °C for 8 min.
A GC/MS-HP7890 system (Agilent, USA) equipped with a 5975C series mass selective detector (MSD) was used to perform the GC/MS analysis. The same GC oven temperature program were used as described above for GC-FID detection. MS data were recorded at 70 eV (EI), m/z (rel. intensity in %) as TIC, total ion current. The data were collected in full scan mode (m/z 50–650). Compounds in the samples were identified by comparison of retention time and mass in the GC/MS spectra to the commercially available chemical standards and mass spectral data of the NIST Standard Reference Database.
Analysis of metabolites by HPLC
Glucose consumption and ethanol accumulation were determined by analyzing the glucose and ethanol concentration in the medium using high-performance liquid chromatography (HPLC) as previously described .
Availability of data and materials
The datasets used or analysed during the current study are available from the corresponding author on reasonable request.
Hunter GLK, Brogden WB Jr. Analysis of the terpene and sesquiterpene hydrocarbons in some citrus oils. J Food Sci. 1965;30:383–7.
Ignea C, Trikka FA, Kourtzelis I, Argiriou A, Kanellis AK, Kampranis SC, Makris AM. Positive genetic interactors of HMG2 identify a new set of genetic perturbations for improving sesquiterpene production in Saccharomyces cerevisiae. Microb Cell Fact. 2012;11:162.
Kolwek J, Behrens C, Linke D, Krings U, Berger RG. Cell-free one-pot conversion of (+)-valencene to (+)-nootkatone by a unique dye-decolorizing peroxidase combined with a laccase from Funalia trogii. J Ind Microbiol Biot. 2018;45:89–101.
Takahashi S, Yeo Y, Greenhagen BT, McMullin T, Song L, Maurina-Brunker J, Rosson R, Noel JP, Chappell J. Metabolic engineering of sesquiterpene metabolism in yeast. Biotechnol Bioeng. 2007;97:170–81.
Beekwilder J, van Houwelingen A, Cankar K, van Dijk AD, de Jong RM, Stoopen G, Bouwmeester H, Achkar J, Sonke T, Bosch D. Valencene synthase from the heartwood of Nootka cypress (Callitropsis nootkatensis) for biotechnological production of valencene. Plant Biotechnol J. 2014;12:174–82.
Frohwitter J, Heider SA, Peters-Wendisch P, Beekwilder J, Wendisch VF. Production of the sesquiterpene (+)-valencene by metabolically engineered Corynebacterium glutamicum. J Biotechnol. 2014;191:205–13.
Ozaydin B, Burd H, Lee TS, Keasling JD. Carotenoid-based phenotypic screen of the yeast deletion collection reveals new genes with roles in isoprenoid production. Metab Eng. 2013;15:174–83.
Asadollahi MA, Maury J, Moller K, Nielsen KF, Schalk M, Clark A, Nielsen J. Production of plant sesquiterpenes in Saccharomyces cerevisiae: effect of ERG9 repression on sesquiterpene biosynthesis. Biotechnol Bioeng. 2008;99:666–77.
Jennings SM, Tsay YH, Fisch TM, Robinson GW. Molecular cloning and characterization of the yeast gene for squalene synthetase. P Natl Acad Sci USA. 1991;88:6038–42.
Jiang Y, Proteau P, Poulter D, Ferro-Novick S. BTS1 encodes a geranylgeranyl diphosphate synthase in Saccharomyces cervisiae. J Biol Chem. 1995;270:21793–9.
Jakočiūnas T, Bonde I, Herrgård M, Harrison SJ, Kristensen M, Pedersen LE, Jensen MK, Keasling JD. Multiplex metabolic pathway engineering using CRISPR/Cas9 in Saccharomyces cerevisiae. Metab Eng. 2015;28:213–22.
Henry KW, Nickels JT, Edlind TD. ROX1 and ERGR regulation in Saccharomyces cerevisiae: implications for antifungal susceptibility. Eukaryot Cell. 2002;1:1041–4.
Montañés FM, Pascual-Ahuir A, Proft M. Repression of ergosterol biosynthesis is essential for stress resistance and is mediated by the Hog1 MAP kinase and the Mot3 and Rox1 transcription factors. Mol Microbiol. 2011;79:1008–23.
Chen Y, Xiao W, Wang Y, Liu H, Li X, Yuan Y. Lycopene overproduction in Saccharomyces cerevisiae through combining pathway engineering with host engineering. Microb Cell Fact. 2016;15:113.
Chen Y, Wang Y, Liu M, Qu J, Yao M, Li B, Ding M, Liu H, Xiao W, Yuan Y. Primary and secondary metabolic effects of a key gene deletion (ΔYPL062W) in metabolically engineered terpenoid-producing Saccharomyces cerevisiae. Appl Environ Microbiol. 2019;85:e01990–2018.
Sun J, Shao Z, Zhao H, Nair N, Wen F, Xu J-H, Zhao H. Cloning and characterization of a panel of constitutive promoters for applications in pathway engineering in Saccharomyces cerevisiae. Biotechnol Bioeng. 2012;109:2082–92.
Alper H, Fischer C, Nevoigt E, Stephanopoulos G. Tuning genetic control through promoter engineering. P Natl Acad Sci USA. 2006;102:12678–83.
Nevoigt E, Kohnke J, Fischer CR, Alper H, Stahl U, Stephanopoulos G. Engineering of promoter replacement cassettes for fine-tuning of gene expression in Saccharomyces cerevisiae. Appl Environ Microb. 2006;72:5266–73.
Ajikumar PK, Xiao W-H, Tyo KEJ, Wang Y, Simeon F, Leonard E, Mucha O, Phon TH, Pfeifer B, Stephanopoulos G. Isoprenoid pathway optimization for taxol precursor overproduction in Escherichia coli. Science. 2010;330:70–4.
Sakuma T, Nishikawa A, Kume S, Chayama K, Yamamoto T. Multiplex genome engineering in human cells using all-in-one CRISPR/Cas9 vector system. Sci Rep. 2014;4:5400.
Iwata S, Yoshina S, Suehiro Y, Hori S, Mitani S. Engineering new balancer chromosomes in C. elegans via CRISPR/Cas9. Sci Rep. 2016;6:33840.
Cong L, Ran FA, Cox D, Lin S, Barretto R, Habib N, Hsu PD, Wu X, Jiang W, Marraffini LA, Zhang F. Multiplex genome engineering using CRISPR-Cas systems. Science. 2013;339:819–23.
Kilby NJ, Snaith MR, Murray JAH. Site-specific recombinases: tools for genome engineering. Trends Genet. 1994;9:413–21.
Kopke K, Hoff B, Kuck U. Application of the Saccharomyces cerevisiae FLP/FRT recombination system in filamentous fungi for marker recycling and construction of knockout strains devoid of heterologous genes. Appl Environ Microb. 2010;76:4664–74.
Zhang DX, Lu HL, Liao X, St Leger RJ, Nuss DL. Simple and efficient recycling of fungal selectable marker genes with the Cre-loxP recombination system via anastomosis. Fungal Genet Biol. 2013;61:1–8.
Yang J, Jiang W, Yang S. mazF as a counter-selectable marker for unmarked genetic modification of Pichia pastoris. FEMS Yeast Res. 2009;9:600–9.
Akada R, Hirosawa I, Kawahata M, Hoshida H, Nishizawa Y. Sets of integrating plasmids and gene disruption cassettes containing improved counter-selection markers designed for repeated use in budding yeast. Yeast. 2002;19:393–402.
Solis-Escalante D, Kuijpers NG, van der Linden FH, Pronk JT, Daran JM, Daran-Lapujade P. Efficient simultaneous excision of multiple selectable marker cassettes using I-SceI-induced double-strand DNA breaks in Saccharomyces cerevisiae. FEMS Yeast Res. 2014;14:741–54.
Li C, Lin Y, Zheng X, Yuan Q, Pang N, Liao X, Huang Y, Zhang X, Liang S. Recycling of a selectable marker with a self-excisable plasmid in Pichia pastoris. Sci Rep. 2017;7:11113.
Ryan OW, Poddar S, Cate JH. CRISPR-Cas9 Genome Engineering in Saccharomyces cerevisiae Cells. Cold Spring Harb Protoc. 2016. https://doi.org/10.1101/pdb.prot086827.
Stovicek V, Borodina I, Forster J. CRISPR–Cas system enables fast and simple genome editing of industrial Saccharomyces cerevisiae strains. Metab Eng Commun. 2015;2:13–22.
Kennedy MA, Bard M. Positive and negative regulation of squalene synthase (ERG9), an ergosterol biosynthetic gene, in Saccharomyces cerevisiae. Biochim Biophys Acta. 2001;1517:177–89.
Lograsso PV, Soltis DA, Boettcher BR. Overexpression, purification, and kinetic characterization of a carboxyl-terminal-truncated yeast squalene synthetase. Arch Biochem Biophys. 1993;307:193–9.
Chang TH, Guo RT, Ko TP, Wang AH, Liang PH. Crystal structure of type-III geranylgeranyl pyrophosphate synthase from Saccharomyces cerevisiae and the mechanism of product chain length determination. J Biol Chem. 2006;281:14991–5000.
Vickers CE, Behrendorff JBYH, Bongers M, Brennan TCR, Bruschi M, Nielsen LK. Production of industrially relevant isoprenoid compounds in engineered microbes. In: Kamm B, editor. Microbiology monographs. Berlin: Springer; 2015. p. 303–34.
Ma T, Shi B, Ye Z, Li X, Liu M, Chen Y, Xia J, Nielsen J, Deng Z, Liu T. Lipid engineering combined with systematic metabolic engineering of Saccharomyces cerevisiae for high-yield production of lycopene. Metab Eng. 2018;52:134–42.
Faulkner A, Chen X, Rush J, Horazdovsky B, Waechter CJ, Carman GM, Sternweis PC. The LPP1 and DPP1 gene products account for most of the isoprenoid phosphate phosphatase activities in Saccharomyces cerevisiae. J Biol Chem. 1999;274:14831–7.
Liu J, Zhang W, Du G, Chen J, Zhou J. Overproduction of geraniol by enhanced precursor supply in Saccharomyces cerevisiae. J Biotechnol. 2013;168:446–51.
Peng B, Plan MR, Chrysanthopoulos P, Hodson MP, Nielsen LK, Vickers CE. A squalene synthase protein degradation method for improved sesquiterpene production in Saccharomyces cerevisiae. Metab Eng. 2017;39:209–19.
Partow S, Siewers V, Bjorn S, Nielsen J, Maury J. Characterization of different promoters for designing a new expression vector in Saccharomyces cerevisiae. Yeast. 2010;27:955–64.
Reifenberger E, Freidel K, Ciriacy M. Identification of novel HXT genes in Saccharomyces cerevisiae reveals the impact of individual hexose transporters on qlycolytic flux. Mol Microbiol. 1995;16:157–67.
Sedlak M, Ho NWY. Characterization of the effectiveness of hexose transporters for transporting xylose during glucose and xylose co-fermentation by a recombinant Saccharomyces yeast. Yeast. 2004;21:671–84.
Asadollahi MA, Maury J, Schalk M, Clark A, Nielsen J. Enhancement of farnesyl diphosphate pool as direct precursor of sesquiterpenes through metabolic engineering of the mevalonate pathway in Saccharomyces cerevisiae. Biotechnol Bioeng. 2010;106:86–96.
Binder D, Frohwitter J, Mahr R, Bier C, Grunberger A, Loeschcke A, Peters-Wendisch P, Kohlheyer D, Pietruszka J, Frunzke J, et al. Light-controlled cell factories: employing photocaged isopropyl-β-d-thiogalactopyranoside for light-mediated optimization of lac promoter-based gene expression and (+)-valencene biosynthesis in Corynebacterium glutamicum. Appl Environ Microb. 2016;82:6141–9.
Scholtmeijer K, Cankar K, Beekwilder J, Wösten BHA, Lugones LG, Bosch D. Production of (+)-valencene in the mushroom-forming fungus S. commune. Appl Microbiol Biotechnol. 2014;98:5059–68.
Verwaal R, Wang J, Meijnen JP, Visser H, Sandmann G, van den Berg JA, van Ooyen AJ. High-level production of beta-carotene in Saccharomyces cerevisiae by successive transformation with carotenogenic genes from Xanthophyllomyces dendrorhous. Appl Environ Microb. 2007;73:4342–50.
Giaever G, Chu AM, Ni L, Connelly C, Riles L, Véronneau S, Dow S, Lucau-Danila A, Anderson K, André B, et al. Functional profiling of the Saccharomyces cerevisiae genome. Nature. 2002;418:387.
Paradise EM, Kirby J, Chan R, Keasling JD. Redirection of flux through the FPP branch-point in Saccharomyces cerevisiae by down-regulating squalene synthase. Biotechnol Bioeng. 2008;100:371–8.
Ro D-K, Paradise EM, Ouellet M, Fisher KJ, Newman KL, Ndungu JM, Ho KA, Eachus RA, Ham TS, Kirby J, et al. Production of the antimalarial drug precursor artemisinic acid in engineered yeast. Nature. 2006;440:940.
Blin K, Pedersen LE, Weber T, Lee SY. CRISPy-web: an online resource to design sgRNAs for CRISPR applications. Synth Syst Biotechnol. 2016;1:118–21.
Zhu M, Fan W, Cha Y, Yang X, Lai Z, Li S, Wang X. Dynamic cell responses in Thermoanaerobacterium sp. under hyperosmotic stress. Sci Rep. 2017;7:10088.
Wriessnegger T, Augustin P, Engleder M, Leitner E, Muller M, Kaluzna I, Schurmann M, Mink D, Zellnig G, Schwab H, Pichler H. Production of the sesquiterpenoid (+)-nootkatone by metabolic engineering of Pichia pastoris. Metab Eng. 2014;24:18–29.
Farhi M, Marhevka E, Masci T, Marcos E, Eyal Y, Ovadis M, Abeliovich H, Vainstein A. Harnessing yeast subcellular compartments for the production of plant terpenoids. Metab Eng. 2011;13:474–81.
This work was financial supported by the Natural Science Foundation of China (Grants 21878104, 31800063 and U1701243), National Key R&D Program of China (2018YFA0901504), the Project on the Integration of Industry, Education, and Research of Guangzhou, China (Grants 201704020183 and 201903010086) and the Science and Technology Planning Project of Guangdong Province, China (Grant 2017A010105019).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Chen, H., Zhu, C., Zhu, M. et al. High production of valencene in Saccharomyces cerevisiae through metabolic engineering. Microb Cell Fact 18, 195 (2019). https://doi.org/10.1186/s12934-019-1246-2