- Open Access
Metabolic engineering of Saccharomyces cerevisiae for high-level production of gastrodin from glucose
Microbial Cell Factories volume 19, Article number: 218 (2020)
The natural phenolic glycoside gastrodin is the major bioactive ingredient in the well-known Chinese herb Tianma and is widely used as a neuroprotective medicine in the clinic. Microbial production from sustainable resources is a promising method to replace plant extraction and chemical synthesis which were currently used in industrial gastrodin production. Saccharomyces cerevisiae is considered as an attractive host to produce natural plant products used in the food and pharmaceutical fields. In this work, we intended to explore the potential of S. cerevisiae as the host for high-level production of gastrodin from glucose.
Here, we first identified the plant-derived glucosyltransferase AsUGT to convert 4-hydroxybenzyl alcohol to gastrodin with high catalytic efficiency in yeast. Then, we engineered de novo production of gastrodin by overexpressing codon-optimized AsUGTsyn, the carboxylic acid reductase gene CARsyn from Nocardia species, the phosphopantetheinyl transferase gene PPTcg-1syn from Corynebacterium glutamicum, the chorismate pyruvate-lyase gene UbiCsyn from Escherichia coli, and the mutant ARO4K229L. Finally, we achieved an improved product titer by a chromosomal multiple-copy integration strategy and enhancement of metabolic flux toward the aglycon 4-hydroxybenzyl alcohol. The best optimized strain produced 2.1 g/L gastrodin in mineral medium with glucose as the sole carbon source by flask fermentation, which was 175 times higher than that of the original gastrodin-producing strain.
The de novo high-level production of gastrodin was first achieved. Instead of chemical synthesis or plants extraction, our work provides an alternative strategy for the industrial production of gastrodin by microbial fermentation from a sustainable resource.
The phenolic glycoside gastrodin is the major bioactive ingredient in the famous Chinese herb Tianma (the tuber of Gastrodia elata Blume), which has been widely used to treat nervous system diseases, such as ischemic stroke, epilepsy, dizziness, and dementia, for centuries in China [1, 2]. In addition to its neuroprotective effects , gastrodin has been proven to show memory improvement , anti-inflammation, and antioxidant effects [5, 6]. To date, gastrodin is produced mainly by chemical synthesis or direct extraction from Tianma, although different processing methods have also been explored and applied [7, 8]. However, total chemical synthesis is a complicated process and leads to severe environmental pollution due to the use of toxic phenols, phosphate and bromide. Direct extraction from Tianma is costly and time consuming because of the low gastrodin content in Tianma (below 0.7%) . Meanwhile, wild resources of Tianma are scarce, and its artificial cultivation is relatively complicated due to the nature of its symbiotic relationship with two compatible mycorrhizal fungi during seed germination and vegetative growth . The success of heterologous production of plant-derived natural products in microorganisms provides an alternative strategy for producing gastrodin in an economical and environmentally friendly manner .
The native biosynthetic pathway of gastrodin in plant has still not been documented. In our previous work , we first achieved microbial biosynthesis of gastrodin by constructing an artificial pathway in Escherichia coli, beginning with the conversion of chorismate to 4-hydroxybenzoic acid (4-HBA) by chorismate pyruvate-lyase (UbiC). Then, 4-HBA was transformed to 4-hydroxybenzyl aldehyde by carboxylic acid reductase (CAR), which is activated into its holo form by a phosphopantetheinyl transferase (sfp). The 4-hydroxybenzyl aldehyde was transformed to 4-hydroxybenzl alcohol by endogenous alcoho9l dehydrogenases (ADHs). Finally, 4-hydroxybenzl alcohol was converted to gastrodin by the glycosyltransferase UGT73B6 from Rhodiola. With further rational engineering of the glycosyltransferase, the resulting E. coli strain produces approximately 0.5 g/L gastrodin. Our previous results strongly suggest that the activity of UGTs is critical for the improvement of gastrodin biosynthesis. Unfortunately, the native glucosyltransferase for gastrodin biosynthesis is unknown. Tsai et al.  reported that a unigene encoding glycosyltransferase might participate in the gastrodin native biosynthetic pathway. However, its sequence has not yet been published. To date, UGT73B6 from Rhodiola is the only glycosyltransferase that has been characterized to be functional in gastrodin synthesis.
As common microbial platforms, E. coli and Saccharomyces cerevisiae have been engineered to produce aromatic amino acid-derived (AAA-derived), chorismate-derived or shikimate-derived products . As in E. coli, some successful strategies have also been used in yeast to enhance the production of these products to gram levels per liter [13, 14]. Liu et al.  constructed an S. cerevisiae platform strain with an AAA-derived p-coumaric acid yield of more than 10 g/L by rewiring carbon metabolism, achieving appropriate levels for large-scale production. Food-grade S. cerevisiae is considered an attractive host due to its superior capability of expressing enzymes requiring posttranslational modification, safety, robust fermentation and resistance to phage infection .
In this work, we explored the potential of S. cerevisiae as the host to produce gastrodin. We first identified a glucosyltransferase with high catalytic efficiency for gastrodin synthesis and constructed a gastrodin-producing yeast strain by pathway engineering. By a chromosomal integration strategy and metabolic flux enhancement toward the 4-hydroxybenzyl alcohol, we achieved high-level production of gastrodin. Engineered yeast strains could serve as the basis for creating an alternative method to produce gastrodin, instead of extraction from plant sources or chemical synthesis in the future.
Results and discussion
The glucosyltransferase AsUGT catalyzed gastrodin biosynthesis with high efficiency in yeast
To achieve heterologous biosynthesis of gastrodin in yeast, it is necessary to identify a compatible glucosyltransferase that recognizes 4-hydroxybenzyl alcohol as a substrate with high efficiency. It was reported that the glucosyltransferase AS from R. serpentina could accept phenols as substrates . Based on the structural similarity between phenols and 4-hydroxybenzyl alcohol, we hypothesized that AS might be able to catalyze the glucosylation of 4-hydroxybenzyl alcohol. To test this hypothesis, we overexpressed a codon-optimized AS (named AsUGTsyn in this work) in yeast and carried out feeding experiments with 4-hydroxybenzyl alcohol as the substrate. HPLC analyses of the metabolites in the fermentation broth supernatants produced by the engineered strain displayed a new peak at a retention time of 6.3 min (Fig. 1a). HPLC–MS analysis showed that the new compound had a molecular ion at m/z 304.1419 ([M + NH4]+), which was identical to that of gastrodin (Fig. 1b). Approximately 0.5 mM gastrodin was produced within 24 h, and the substrate conversion rate was 25%. Within 48 h, the substrate conversion rate reached up to 55%, with approximately 1.1 mM gastrodin produced.
Our previous study demonstrated that the glucosyltransferase UGT73B6FS, a mutant of UGT73B6, showed high catalytic efficiency in the conversion of 4-hydroxybenzyl alcohol to gastrodin in E. coli . To compare enzymatic efficiency between codon-optimized AsUGTsyn and UGT73B6FSy in yeast, ScTEF1 promoter-driven UGT73B6FSy was also introduced into S. cerevisiae BY4742. The resulting strain 4742-UGT73B6FSy was cultivated under the same conditions. As shown in Fig. 1c, strain 4742-UGT73B6FSy produced only trace amounts of gastrodin. The substrate conversion rate in strain 4742-AsUGTsyn was 33 times higher than that in strain 4742- UGT73B6FSy within 24 h. Taken together, these results clearly showed that AsUGT exhibited a higher catalytic efficiency than UGT73B6FS in yeast. Thus, AsUGT was used for heterologous gastrodin biosynthesis in subsequent experiments.
De novo biosynthesis of gastrodin from glucose was achieved in S. cerevisiae
As shown in our previous study, the artificial gastrodin biosynthesis pathway was extended from chorismate . In bacteria, chorismate lyase (UbiC) usually catalyzed chorismate to synthesize 4-HBA and the amounts of 4-HBA always low in cultures . In S. cerevisiae, a homolog of UbiC is missing, and the biosynthetic pathway of 4-HBA is unclear . There was no detectable 4-HBA accumulation in S. cerevisiae . Overexpression of the feedback inhibition resistance gene ARO4K229L is usually used to improve the chorismate pathway metabolic flux . In this work, we intended to overexpress ARO4K229L and UbiC from E. coli to preliminarily produce 4-HBA from chorismate in yeast. Coexpression of CAR from the Nocardia genus and PPTcg-1 from C. glutamicum could effectively reduce protocatechuic acid to protocatechuic aldehyde in S. cerevisiae . We tried to overexpress CAR from the Nocardia genus and PPTcg-1 from C. glutamicum in yeast to realize the reduction of 4-HBA to aldehyde. The subsequent reduction of aldehyde to 4-hydroxybenzylalcohol might be achieved by endogenous alcohol dehydrogenase (ADHs) in yeast. The de novo biosynthetic pathway of gastrodin in S. cerevisiae is shown in Fig. 2. We constructed 2 μ-based plasmids carrying strong constitutive promoters to realize target gene overexpression. The plasmid pCf302-CP and the plasmid pCf301-AUA were introduced into S. cerevisiae BY4742, and the resulting strain was named 4742-pGS. We observed two clear new compound peaks at retention times of 6.3 min (peak I) and 12.6 min (peak II) in the HPLC spectrum of the fermentation broth supernatant of 4742-pGS compared to that of control 4742-pCf (Fig. 3a). LC–MS analysis showed that peak I had a molecular ion at m/z 304.1419 ([M + NH4]+), which was identical to that of standard gastrodin (Fig. 3b). Peak II had a molecular ion at m/z 107.0515 ([M-H2O + H]+), which was identical to that of standard 4-hydroxybenzyl alcohol (Fig. 3c). The results clearly showed that gastrodin and 4-hydroxybenzylalcohol were produced in strain 4742-pGS, demonstrating that the pathway genes CARsyn, PPTcg-1syn and AsUGTsyn were functional in yeast. 4-Hydroxybenzylalcohol was synthesized from 4-HBA by CAR reduction, and subsequently, a portion of the 4-hydroxybenzylalcohol was glycosylated to gastrodin by AsUGT in S. cerevisiae BY4742. The overexpression of ARO4K229L and E. coli UbiC might improve the production of 4-HBA from chorismate.
Gastrodin production was improved by chromosomal integration of biosynthetic pathway genes in the rDNA locus
Stable integration and robust expression of foreign genes is critical for the success of heterologous biosynthesis of valuable compounds. It is challenging to stably express multiple genes in S. cerevisiae with the generally used CEN ori-based and 2 µm ori-based plasmids as tools [22, 23]. Moreover, plasmid-based pathway gene expression could result in growth and genetic instability . Chromosomal integration is an ideal method for introducing heterologous genes to achieve robust expression because of its structural stability . The rDNA encompasses approximately 140 copies of tandem repeats on chromosome XII , and δ sequences exist in over 100 copies on chromosome XV . The rDNA locus or δ sites are ideal for generating multicopy integration. Some genes associated with natural products, such as ginsenosides , glycyrrhetinic acid , and salidroside , have been successfully integrated into the rDNA locus or δ sites in yeast to achieve multiple-copy and stable expression by different biotechnological methods.
We aimed to simultaneously integrate CARsyn, PPTcg-1syn, AsUGTsyn, ubiCsyn and ARO4K229L into the yeast chromosomal rDNA locus in the expected order by one-step transformation. The arrangement of these genes in the rDNA locus is shown in Fig. 2. In S. cerevisiae, the knockout of chorismate mutase (ARO7) could block chorismate flux to prephenate and activate the shikimate pathway . Accordingly, we constructed the aro7Δ mutant in S. cerevisiae BY4742 as the parental strain. First, we PCR-amplified two DNA fragments containing the five target genes and URA3 auxotrophic marker from constructed plasmids in the work. Together with the two fragments rDNA-up and rDNA-down, these DNA fragments (approximately 14 kb) were simultaneously transformed into the strain S. cerevisiae aro7Δ by a one-step multiple-fragment yeast transformation protocol. The assembly order and integration of the DNA fragments into the rDNA locus were confirmed by diagnostic PCR and sequencing. Positive colonies were selected and cultured to determine gastrodin production, as shown in Additional file 1: Figure S1. The resulting strain with the highest concentration of gastrodin was named rGS3.
Another, we introduced the plasmids pCf302-CP and pCf301-AUA carrying CARsyn, PPTcg-1syn, AsUGTsyn, ubiCsyn and ARO4K229L into S. cerevisiae aro7Δ. The generated plasmid-based strain was named Δ7-pGS. Gastrodin production in rGS3 and Δ7-pGS was measured under the same fermentation conditions. Based on the production curves (Fig. 4a), gastrodin production in strain Δ7-pGS reached approximately 120 mg/L and was nearly stable after four days of fermentation. The gastrodin titer of strain rGS3 reached approximately 420 mg/L after 6 days of culture, which was 3.8 times higher than that of the plasmid-based strain Δ7-pGS. HPLC spectra are shown in Fig. 4b. The growth curves showed that there was a growth decline after two days of cultivation of Δ7-pGS, and the OD600 of strain rGS3 was 30% higher than that of Δ7-pGS at the end of fermentation (Fig. 4a). These results confirmed that the gene chromosomal integration strategy effectively improved gastrodin production and brought growth recovery compared to plasmid-based pathway gene expression.
Further improving the precursor supply increased gastrodin production
Enhancement of the precursor supply is a classic and effective strategy to improve final compound production. The initiation of chorismate biosynthesis requires two precursors: erythrose-4-phosphate (E4P) derived from the pentose phosphate pathway and phosphoenolpyruvate (PEP) derived from glycolysis. In S. cerevisiae, E4P and PEP were not abundant and were pathway bottlenecks based on the results of carbon tracing and 13C metabolic flux analysis . Many studies have revealed that the overexpression of transketolase (TktA) and phosphoenolpyruvate synthase (PpsA) effectively raised the intracellular pool of E4P and PEP [32, 33]. In yeast, the 3-deoxyD-arabino-heptulosonate-7-phosphate (DAHP) synthase (ARO4/ARO3) catalyzes the first step to condense PEP and E4P to form DAHP. The pentafunctional enzyme (ARO1) catalyzes five subsequent reactions from DAHP to 5-enolpyruvyl-shikimate-3-phosphate (EPSP). Chorismate synthase (ARO2) catalyzes the final conversion of EPSP to chorismate. Previous reports showed that co-overexpression of ARO2, ARO1 and ARO4 (or ARO3) could increase AAA-derived aromatic chemical production .
Given that a deficient supply of precursors might be the cause of the low amount of gastrodin, we intended to further enhance the carbon flux into 4-HBA precursor accumulation by combining the overexpression of ppsA, tktA, ARO1, and ARO2 using the multiple-copy integration strategy. First, we PCR-amplified three large DNA fragments containing ppsA, tktA, ARO1, ARO2, the mutated ARO4K229L and ubiCsyn from constructed plasmids in the work. Together with the two fragments δ DNA-up and δ DNA-down, these DNA fragments (approximately 21 kb) were simultaneously transformed into the strain S. cerevisiae aro7Δ by a one-step multiple-fragment yeast transformation protocol. The arrangement of these genes in δ sites is shown in Fig. 2. The assembly order and integration of the DNA fragments into the δ sites were confirmed by diagnostic PCR and sequencing. Positive colonies were selected and cultured. HPLC analysis of the culture of these colonies showed that there was a clear peak I at a retention time of 28.5 min, similar to the 4-HBA standard, while peak I was absent in S. cerevisiae aro7Δ (Fig. 5a). We collected the fraction of peak I for LC–MS in negative ion mode in a methanol–water (containing 20 mM NH4Ac) gradient system. LC–MS analysis showed that peak I had a molecular ion at m/z 137.0189 ([M-H]−), which was identical to that of standard 4-HBA (Fig. 5b), further confirming 4-HBA overproduction. Production of 4-HBA in some positive colonies was shown in Additional file 1: Figure S2. The colony with the highest 4-HBA production was named Δ7-HBA. Another, consistent with a previous report , there was no detectable 4-HBA accumulation in S. cerevisiae aro7Δ. These results demonstrated that the combined overexpression of ARO1, ARO2, mutated ARO4K229L, ppsA, tktA and ubiCsyn effectively enhanced metabolic flux toward 4-HBA biosynthesis in yeast.
Subsequently, these five DNA fragments were transformed into S. cerevisiae rGS3. Gastrodin production in the positive colonies was measured (Additional file 1: Figure S3), and the colony with the highest gastrodin production was named rGS-HBA. To evaluate the product titers and cell density of the engineered strains, we carried out fermentation in mineral medium with glucose as the sole carbon source by a 250-mL flask fermentation. The titers of gastrodin, 4-hydroxybenzyl alcohol and 4-HBA were determined after 6 days of fermentation. As shown in Fig. 6, the original gastrodin-producing strain 4742-pGS had 12 mg/L gastrodin, accumulated 8 mg/L 4-hydroxybenzyl alcohol and 120 mg/L 4-HBA. In strain 4742-pGS, the accumulation of 4-HBA showed that the overexpression of ARO4 K229L and E. coli UbiC indeed enhanced 4-HBA production. Compared to 4742-pGS, the production of all the three compounds was increased in Δ7-pGS, indicating that deletion of ARO7 activated the shikimate pathway and promoted the biosynthesis of gastrodin, 4-hydroxybenzyl alcohol and 4-HBA. The cell density of Δ7-pGS (OD600 = 5.8) was lower than that of 4742-pGS (OD600 = 9.0), probably because deletion of ARO7 affected strain growth. The plasmid-free strains rGS, Δ7-HBA and rGS-HBA showed similar cell densities, higher than that of plasmid-based Δ7-pGS. In addition to gastrodin, the total amount of 4-hydroxybenzylalcohol produced by plasmid-free rGS3 was 3 times that produced by the plasmid-based strain Δ7-pGS. These results further confirmed that the genes integrating expression in yeast were superior to plasmid-based expression in gastrodin production. With the cooverexpression of ppsA, tktA, ARO1, ARO2, ARO4K229L and ubiCsyn integrated in δ sites, 1.3 g/L 4-HBA accumulation was observed in Δ7-HBA. The enhanced pool of 4-HBA further increased the metabolic flux toward the aglycon 4-hydroxybenzyl alcohol, and consequently, high gastrodin production (2.1 g/L) was achieved in strain rGS-HBA, which was approximately 4.8 times that of S. cerevisiae rGS. With these rational engineering, the titer of gastrodin was improved 175 times higher than that of the original 4742-pGS.
The expression levels of target genes were detected between engineered strains by quantitative real-time PCR. Relative transcript levels were analyzed individually after normalization to the actin internal reference gene. Results were as shown in Additional file 1: File S1The analysis of expression levels of target genes. The changes of transcript expression levels of target genes showed good consistency with the production of gastrodin between engineered strains.
Here, we explored the potential of S. cerevisiae as the host for high-level production of gastrodin from glucose. The glucosyltransferase AsUGT was firstly identified to convert 4-hydroxybenzylalcohol to gastrodin with high catalytic efficiency in yeast. The de novo biosynthesis of gastrodin from glucose was realized by introducing the AsUGTsyn and other pathway genes into S. cerevisiae. With multiple-copy integration of the entire pathway genes in the rDNA locus and enhancement of metabolic flux toward 4-hydroxybenzyl alcohol, the best optimized strain produced 2.1 g/L gastrodin in mineral medium with glucose as the sole carbon source. To our knowledge, this is the highest reported titer of gastrodin produced by an engineered microorganism. Instead of chemical synthesis or plants extraction, this work provides a strong basis for the biotechnological production of gastrodin at an industrial scale for pharmaceutical applications.
Material and methods
Strains and medium
All plasmids and strains used in this work are shown in Table 1. E. coli DH5α was used for plasmid construction. E. coli strains were cultivated in Luria broth (LB). S. cerevisiae BY4742 (MATα his3Δ1 leu2Δ0 lys2Δ0 ura3Δ0) was used as the parental strain for all the yeast strains. Yeast strains were grown in SC dropout medium containing 6.7 g/L yeast nitrogen base, 2 g/L dropout mix, and 20 g/L glucose. SC dropout medium without uracil (SC-ura), SC dropout medium without leucine (SC-leu), and SC dropout medium without uracil and leucine (SC-ura, leu) were used in this work. Yeast strains were also cultured in YPD medium containing 10 g/L yeast extract, 20 g/L peptone and 20 g/L glucose. Standards, including gastrodin, 4-hydroxybenzl alcohol and 4-HBA, were purchased from Aladdin Chemistry Co., Ltd. (Shanghai, China).
Genes and enzymes
Glucosyltransferase AsUGT (Q9AR73)  from Rauvolfia serpentina (serpentwood), Rhodiola-derived glycosyltransferase UGT73B6 mutant gene UGT73B6FS , carboxylic acid reductase CAR (Q6RKB1.1) gene  from Nocardia iowensis, phosphopantetheinyl transferase PPTcg-1 (WP_003857486.1) gene  from Corynebacterium glutamicum, and chorismate pyruvate-lyase UbiC (AAY88959.1)  gene from E. coli BL21 (DE3) were synthesized with S. cerevisiae codon optimization by Generay Biotech Co., Ltd. (Shanghai, China). The feedback-resistant DAHP synthase mutant ARO4K229L from S. cerevisiae was generated from DAHP synthase ARO4 (NP_009808) by a site-directed mutagenesis method. The pentafunctional EPSP synthase gene ARO1 (Gene ID: 851705) and chorismate synthase gene ARO2 (GeneID: 852729) were amplified from S. cerevisiae BY4742. The phosphoenolpyruvate synthase gene ppsA (Gene ID: 946209) and transketolase gene tktA (Gene ID: 947420) were amplified from E. coli MG1655. The synthesized genes and sequences applied in this study are showed in Additional file 1: Table S2.
All the primers used in this study are listed in Additional file 1: Table S1.
The codon-optimized genes AsUGTsyn and UGT73B6FSy in S. cerevisiae were separately cloned into the pCf301  vector using NheI/AatII to create the plasmids pCf301- AsUGTsyn and pCf301- UGT73B6FSy, respectively. The codon-optimized genes CARsyn and PPTcg-1syn were respectively cloned into the pCf302  vector using NheI/AatII and HindIII/BamHI to create the pCf302-CP plasmids. The ARO4K229L, codon-optimized AsUGTsyn, and ubiCsyn genes were respectively cloned into the pCf301 vector using SpeI/BglII, NheI/AatII and HindIII/BamHI to create the pCf301-AUA plasmids. The genes ARO4K229L and ubiCsyn were also respectively cloned into the pCf302 vector using SpeI/BglII and HindIII/BamHI to create the pCf302-AU plasmids. The ppsA and tktA genes were respectively cloned into the pCf301 vector using NheI/AatII and HindIII/BamHI to create the pCf301-pt plasmids. The genes ARO1 and ARO2 were respectively cloned into the pCf301 vector using NheI/AatII and HindIII/BamHI to create the pCf301-ARO1/2 plasmids.
Chromosomal disruption of ARO7
Targeted chromosomal disruption of the chorismate mutase gene ARO7 in the strain S. cerevisiae BY4742 was performed by homologous recombination using antibiotic marker-containing disruption cassettes created by PCR, as described by Gueldener et al. . The loxP-flanked kanMX cassettes were amplified from pUG6. The antibiotic marker was removed using the Cre recombinase plasmid pSH65. Clones successfully carrying the ARO7 disruption cassette were further confirmed by colony PCR and named S. cerevisiae BY4742 aro7Δ.
The S. cerevisiae BY4742 strains respectively harboring pCf301, pCf301-UGT73B6 FSy or pCf301-AsUGTsyn were cultured in 3 mL of SC-ura medium at 30 °C and 220 r.p.m. for 24 h. Then, 2 mL of the culture was diluted into 50 mL of medium and shaken at 30 °C at 220 r.p.m. for 18 h. Then, the substrate 4-hydroxybenzyl alcohol was added into the cultures at a final concentration of 2 mM. After 24 h and 48 h, samples were taken and analyzed by HPLC or HPLC–MS. Three replicates were performed in experiments. The data were presented as means ± SD (standard deviation).
All measurements of optical density at 600 nm (OD600) were made using a SHIMADZU UV-1800 spectrophotometer. A single colony was first inoculated into 35-mL culture tubes containing 3 mL of SC dropout medium and grown at 30 °C and 220 r.p.m. for 18 h. Then, seed cultures were transferred into 50 mL of SC dropout medium to obtain an OD600 of 0.2 in 250-mL culture flasks. The strains were cultured at 30 °C and 220 r.p.m. for approximately 7 days. Samples were taken every 24 h and analyzed by HPLC. All of the data generated in this study were collected from three independent experiments. The data were presented as means ± SD (standard deviation).
One-step multiple-fragment yeast transformation
The procedure of multiple-fragment transformation was based on the method described by Jolanda van Leeuwen for rapid and efficient plasmid construction . A single colony was inoculated into 3 mL of liquid medium and incubated overnight at 200 r.p.m. and 30 °C. Then, 2 mL of the overnight culture was diluted into 50 mL of medium to obtain an OD600 of 0.2 and incubated at 30 °C and 200 r.p.m. The culture with an OD600 of 0.8 was centrifuged at 2000g for 3 min, and the cell pellet was suspended in 25 mL of sterile water. After centrifugation, the cells were suspended in 1.0 mL of 0.1 M lithium acetate (LiAc). The cell suspension was transferred to a 1.5-mL microcentrifuge tube and centrifuged for 15 s at 2000g. After discarding the supernatant, the cells were suspended in 320 μL of LiAc (0.1 M) as competent cells. Ten microliters of DNA fragment mixture (100–200 ng of each DNA fragment), 2 µL of ssDNA (10 mg/mL) and 12 µL of competent yeast cells were added into a 1.5-mL microcentrifuge tube and mixed gently. Then, 100 µL of the transformation buffer (800 µL of 50% PEG3350, 100 µL of 1 M LiAc, 100 µL of 10 × TE, and 50 µL of DMSO) was added into the 1.5-mL microcentrifuge tube and vortexed for 10 s. After incubation for 30 min at 25 °C, the transformation mixture was incubated at 42 °C for 15 min and then on ice for 5 min. The transformation mixture was centrifuged at 800g for 1 min, and the cells were suspended in 1 mL of YPD. The cell suspension was incubated at 30 °C for 2 h and was centrifuged at 2000g for 1 min. The cells were suspended in sterile water and plated on SC plates lacking leucine and/or uracil. Plasmid transformation was carried out using the standard LiAc method.
Chromosomal multiple-copy integration of the pathway genes for gastrodin biosynthesis
The CARsyn, PPTcg-1syn, AsUGTsyn, ubiCsyn and ARO4K229L, URA3 selectable marker, a total sequence of approximately 14 kb, were integrated into the ribosomal DNA (rDNA) locus. The two fragments rDNA-up (476 bp) and rDNA-down (326 bp) were PCR-amplified from the DNA of S. cerevisiae BY4742 using the primer set rU-F&R and rD-F&R. An approximately 5.8-kb fragment containing CARsyn and PPTcg-1syn was PCR-amplified from the pCf302-CP plasmids using the primer set r5.8-F & R. An approximately 7.4-kb fragment containing AsUGTsyn, ubiCsyn, ARO4K229L and URA3 was PCR-amplified from the pCf301-AUA plasmids using the primer set r7.4-F & R. The four DNA fragments, with 30–40 bp homologous regions designed between the adjacent fragments, were transferred into S. cerevisiae aro7Δ using the protocol for one-step multiple-fragment transformation mentioned above. The colonies were selected on the SC-ura plates and verified by PCR amplification and sequencing. Positive colonies were cultured, and the highest gastrodin-producing strain was designated rGS3.
The ARO1, ARO2, ppsA, tktA, ubiCsyn, ARO4K229L and LEU2 selectable marker, a total sequence length of approximately 21 kb, were integrated separately into the δ-sites of the strains S. cerevisiae aro7Δ and rGS3. The two fragments δ DNA-up (329 bp) and δ DNA-down (316 bp) were PCR-amplified from the DNA of S. cerevisiae BY4742 using the primer set δ U-F&R and δ D-F&R. A 7.7-kb fragment containing ARO1 and ARO2 was PCR-amplified from the pCf301-ARO1/2 plasmid using the primer set δ 7.7-F&R. An approximately 6.8-kb fragment containing LEU2, ubiCsyn, and ARO4K229L was PCR-amplified from the pCf302-AU plasmid using the primer set δ 6.8-F&R. A 6.2-kb fragment containing ppsA and tktA was amplified from the pCf301-pt plasmid using the primer set δ 6.2-F&R. The five DNA fragments, with 30–40 bp homologous regions designed between the adjacent fragments, were transformed into the strain S. cerevisiae aro7Δ using the protocol for one-step multiple-fragment transformation mentioned above. Positive colonies were selected on the SC-leu plates, and the strain with the highest 4-HBA production was designated strain Δ7-HBA. The five DNA fragments were also transformed into the strain rGS3. Colonies were selected on the SC-ura, leu plates, and the strain with the highest gastrodin production was designated strain rGS-HBA. The sequences of the rDNA and δ DNA fragments applied are shown in Additional file 1: Table S2.
The chemicals gastrodin, 4-hydroxybenzl alcohol and 4-HBA were analyzed by HPLC on a Shimadzu system with a UV detector. The column used was an Agela Innoval C18, 4.6 × 250 mm, with a particle size of 5 μm. A methanol–water (containing 0.1% trifluoroacetic acid) gradient system at a flow rate of 1 mL/min was used to separate the compounds. The HPLC program was as follows: 10% methanol for 40 min, 10% to 100% methanol for 1 min, 100% methanol for 10 min, 100% to 10% methanol for 1 min, and 10% methanol for 10 min. Products were detected at 224 nm or 201 nm. HPLC–MS was performed on an Agilent 1260 system with a 1260 Infinity UV detector and a Bruker micro Q-TOF II mass spectrometer equipped with an ESI ionization probe. The HPLC conditions were the same as described above. For HPLC–MS analysis of 4-HBA, the fraction of the peak was collected for MS identification in negative ion mode, and methanol–water (containing 20 mM NH4Ac) was used as the gradient system. The other conditions were the same as described above. To quantify the titers of gastrodin, 4-hydroxybenzl alcohol and 4-HBA, standard calibration curves were prepared with a series of known concentrations of the standard. R2 value for the standard curve was higher than 0.999. Samples were diluted as required for falling into the concentration range of the calibration curves. The titers were presented as means ± SD (standard deviation).
RNA isolation and quantitative real-time PCR
Total RNA from S. cerevisiae was isolated and purified using hot-acid-phenol method. To synthesize cDNA, about 2.0 μg of RNA was added to a FastKing gDNA Dispelling RT SuperMix (TIANGEN Biotech, Beijing, China). Each 20.0 μL reaction contained 10.0 μL of SuperReal PreMix Plus (SYBR green) (TIANGEN Biotech, Beijing, China). Real-time PCR was monitored and analyzed using Roche LightCycler® 96 System. Samples of Δ7-pGS and rGS3 were prepared after 2 days and 4 days fermentation culture. Samples of rGS3-PHB were prepared after 2 days. All samples were prepared in triplicate. The relative expression level of genes in engineered strain was calculated by the relative quantitative 2−ΔΔCT method  and ß-actin was used as the internal control gene Additional file 1: File S1.
Availability of data and materials
All data generated or analyzed during this study are included in this article and its Additional file 1.
Matias M, Silvestre S, Falcao A, Alves G. Gastrodia elata and epilepsy: Rationale and therapeutic potential. Phytomedicine. 2016;23:1511–26.
Chinese Pharmacopoeia Commission. Pharmacopoeia of the People’s Republic of China. Chinese Medical Science and Technology Press (Beijing). 2010; 54–55.
Liu Y, Gao J, Peng M, Meng H, Ma H, Cai P, Xu Y, Zhao Q, Si G. A review on central nervous system effects of gastrodin. Front Pharmacol. 2018;9:24.
Deng CK, Mu ZH, Miao YH, Liu YD, Zhou L, Huang YJ, Zhang F, Wang YY, Yang ZH, Qian ZY, et al. Gastrodin ameliorates motor learning deficits through preserving cerebellar long-term depression pathways in diabetic rats. Front Neurosci. 2019;13:1239.
Jiang T, Chu J, Chen H, Cheng H, Su J, Wang X, Cao Y, Tian S, Li Q. Gastrodin Inhibits H2O2-Induced Ferroptosis through Its Antioxidative Effect in Rat Glioma Cell Line C6. Biol Pharm Bull. 2020;43:480–7.
Peng Z, Wang S, Chen G, Cai M, Liu R, Deng J, Liu J, Zhang T, Tan Q, Hai C. Gastrodin alleviates cerebral ischemic damage in mice by improving anti-oxidant and anti-inflammation activities and inhibiting apoptosis pathway. Neurochem Res. 2015;40:661–73.
Hu M, Yan H, Fu Y, Jiang Y, Yao W, Yu S, Zhang L, Wu Q, Ding A, Shan M. Optimal Extraction Study of Gastrodin-Type Components from Gastrodia Elata Tubers by Response Surface Design with Integrated Phytochemical and Bioactivity Evaluation. Molecules. 2019;24:547.
Liu Y, Huang G. The content analysis of gastrodin and gastrodigenin obtained by different processing methods. J Chromatogr Sci. 2018;56:65–7.
Tsai CC, Wu KM, Chiang TY, Huang CY, Chou CH, Li SJ, Chiang YC. Comparative transcriptome analysis of Gastrodia elata (Orchidaceae) in response to fungus symbiosis to identify gastrodin biosynthesis-related genes. BMC Genomics. 2016;17:212.
Borodina I, Nielsen J. Advances in metabolic engineering of yeast Saccharomyces cerevisiae for production of chemicals. Biotechnol J. 2014;9:609–20.
Bai Y, Yin H, Bi H, Zhuang Y, Liu T, Ma Y. De novo biosynthesis of Gastrodin in Escherichia coli. Metab Eng. 2016;35:138–47.
Cao M, Gao M, Suastegui M, Mei Y, Shao Z. Building microbial factories for the production of aromatic amino acid pathway derivatives: from commodity chemicals to plant-sourced natural products. Metab Eng. 2020;58:94–132.
Gao S, Lyu Y, Zeng W, Du G, Zhou J, Chen J. Efficient Biosynthesis of (2S)-Naringenin from p-Coumaric Acid in Saccharomyces cerevisiae. J Agric Food Chem. 2020;68:1015–21.
Averesch NJH, Kromer JO. Metabolic engineering of the shikimate pathway for production of aromatics and derived compounds-present and future strain construction strategies. Front Bioeng Biotechnol. 2018;6:32.
Liu Q, Yu T, Li X, Chen Y, Campbell K, Nielsen J, Chen Y. Rewiring carbon metabolism in yeast for high level production of aromatic chemicals. Nat Commun. 2019;10:4976.
Nielsen J. Yeast systems biology. Model organism and cell factory. Biotechnol J. 2019;14:e1800421.
Hefner T, Arend J, Warzecha H, Siems K, Stockigt J. Arbutin synthase, a novel member of the NRD1beta glycosyltransferase family, is a unique multifunctional enzyme converting various natural products and xenobiotics. Bioorg Med Chem. 2002;10:1731–41.
Winter G, Averesch NJ, Nunez-Bernal D, Kromer JO. In vivo instability of chorismate causes substrate loss during fermentative production of aromatics. Yeast. 2014;31:333–41.
Marbois B, Xie LX, Choi S, Hirano K, Hyman K, Clarke CF. para-Aminobenzoic acid is a precursor in coenzyme Q6 biosynthesis in Saccharomyces cerevisiae. J Biol Chem. 2010;285:27827–38.
Kromer JO, Nunez-Bernal D, Averesch NJ, Hampe J, Varela J, Varela C. Production of aromatics in Saccharomyces cerevisiae–a feasibility study. J Biotechnol. 2013;163:184–93.
Hansen EH, Moller BL, Kock GR, Bunner CM, Kristensen C, Jensen OR, Okkels FT, Olsen CE, Motawia MS, Hansen J. De novo biosynthesis of vanillin in fission yeast (Schizosaccharomyces pombe) and baker’s yeast (Saccharomyces cerevisiae). Appl Environ Microbiol. 2009;75:2765–74.
Romanos MA, Scorer CA, Clare JJ. Foreign gene expression in yeast: a review. Yeast. 1992;8:423–88.
Erhart E, Hollenberg CP. The presence of a defective LEU2 gene on 2 mu DNA recombinant plasmids of Saccharomyces cerevisiae is responsible for curing and high copy number. J Bacteriol. 1983;156:625–35.
Englaender JA, Jones JA, Cress BF, Kuhlman TE, Linhardt RJ, Koffas MAG. Effect of genomic integration location on heterologous protein expression and metabolic engineering in E. coli. ACS Synth Biol. 2017;6:710–20.
Lv Y, Edwards H, Zhou J, Xu P. Combining 26s rDNA and the Cre-loxP System for Iterative Gene Integration and Efficient Marker Curation in Yarrowia lipolytica. ACS Synth Biol. 2019;8:568–76.
Wery J, Gutker D, Renniers AC, Verdoes JC, van Ooyen AJ. High copy number integration into the ribosomal DNA of the yeast Phaffia rhodozyma. Gene. 1997;184:89–97.
Parekh RN, Shaw MR, Wittrup KD. An integrating vector for tunable, high copy, stable integration into the dispersed Ty delta sites of Saccharomyces cerevisiae. Biotechnol Prog. 1996;12:16–21.
Dai Z, Liu Y, Zhang X, Shi M, Wang B, Wang D, Huang L, Zhang X. Metabolic engineering of Saccharomyces cerevisiae for production of ginsenosides. Metab Eng. 2013;20C:146–56.
Wang C, Su X, Sun M, Zhang M, Wu J, Xing J, Wang Y, Xue J, Liu X, Sun W, Chen S. Efficient production of glycyrrhetinic acid in metabolically engineered Saccharomyces cerevisiae via an integrated strategy. Microb Cell Fact. 2019;18:95.
Jiang JJ, Yin H, Wang SA, Zhuang YB, Liu SW, Liu T, Ma YH. Metabolic Engineering of Saccharomyces cerevisiae for High-Level Production of Salidroside from Glucose. J Agric Food Chem. 2018;66:4431–8.
Suastegui M, Guo W, Feng X, Shao Z. Investigating strain dependency in the production of aromatic compounds in Saccharomyces cerevisiae. Biotechnol Bioeng. 2016;113:2676–85.
Gao M, Cao M, Suastegui M, Walker J, Rodriguez Quiroz N, Wu Y, Tribby D, Okerlund A, Stanley L, Shanks JV, Shao Z. Innovating a nonconventional yeast platform for producing shikimate as the building block of high-value aromatics. ACS Synth Biol. 2017;6:29–38.
Koma D, Kishida T, Yoshida E, Ohashi H, Yamanaka H, Moriyoshi K, Nagamori E, Ohmoto T. Chromosome engineering to generate plasmid-free phenylalanine- and tyrosine-overproducing Escherichia coli which can be applied in generation of aromatic compound-producing strains. Appl Environ Microbiol. 2020. https://doi.org/10.1128/AEM.00525-20.
Siebert M, Severin K, Heide L. Formation of 4-Hydroxybenzoate in Escherichia Coli Characterization of the Ubic Gene and Its Encoded Enzyme Chorismate Pyruvate-Lyase. Microbiology-Sgm. 1994;140:897–904.
Gueldener U, Heinisch J, Koehler GJ, Voss D, Hegemann JH. A second set of loxP marker cassettes for Cre-mediated multiple gene knockouts in budding yeast. Nucleic Acids Res. 2002;30:e23.
Van Leeuwen J, Andrews B, Boone C, Tan G. Rapid and efficient plasmid construction by homologous recombination in yeast. Cold Spring Harb Protoc. 2015;2015:pdb prot085100.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3:1101–8.
This work was supported by grants from the National Natural Science Foundation of China (31400026, 31970065, U1902214), the Yunnan Key Research and Development Program (2019ZF011-2), and the Tianjin Synthetic Biotechnology Innovation Capacity Improvement Project (TSBICIP-KJGG-002). We thank Jie Zhang at Laboratory of Structural Biology (TIB, CAS) for technical help with RT-PCR experiments.
Ethics approval and consent to participate
Consent for publication
Part of the work has been included in a patent application by the Tianjin Institute of Industrial Biotechnology.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The analysis of expression levels of target genes. Figure S1. Titer of gastrodin produced by colonies with gene integration into the yeast chromosomal rDNA locus. Figure S2. Titer of 4-HBA accumulated in the culture of colonies with gene integration into yeast chromosomal δ sites. Figure S3. Titer of gastrodin produced by colonies with gene integration into yeast chromosomal δ sites. Figure S4. The relative transcript levels of target genes in rGS3 compared to Δ7-pGS at 2 day and 4 day. Figure S5. The relative transcript levels of target genes in rGS3-PHB compared to rGS3 at 2 day. Table S1. The primers used in this study. Table S2. The synthesized genes, sequence of rDNA and δ DNA fragment applied in this study. Table S3. The real-time PCR primers used in this study. Table S4. Information of genes used in this study.
About this article
Cite this article
Yin, H., Hu, T., Zhuang, Y. et al. Metabolic engineering of Saccharomyces cerevisiae for high-level production of gastrodin from glucose. Microb Cell Fact 19, 218 (2020). https://doi.org/10.1186/s12934-020-01476-0
- Saccharomyces cerevisiae
- Phenolic glycoside
- Chromosomal integration
- Metabolic engineering