Skip to main content
Download PDF

Boosting the biosynthesis of betulinic acid and related triterpenoids in Yarrowia lipolytica via multimodular metabolic engineering

Microbial Cell Factories volume 18, Article number: 77 (2019) Cite this article

Abstract

Background

Betulinic acid is a pentacyclic lupane-type triterpenoid and a potential antiviral and antitumor drug, but the amount of betulinic acid in plants is low and cannot meet the demand for this compound. Yarrowia lipolytica, as an oleaginous yeast, is a promising microbial cell factory for the production of highly hydrophobic compounds due to the ability of this organism to accumulate large amounts of lipids that can store hydrophobic products and supply sufficient precursors for terpene synthesis. However, engineering for the heterologous production of betulinic acid and related triterpenoids has not developed as systematically as that for the production of other terpenoids, thus the production of betulinic acid in microbes remains unsatisfactory.

Results

In this study, we applied a multimodular strategy to systematically improve the biosynthesis of betulinic acid and related triterpenoids in Y. lipolytica by engineering four functional modules, namely, the heterogenous CYP/CPR, MVA, acetyl-CoA generation, and redox cofactor supply modules. First, by screening 25 combinations of cytochrome P450 monooxygenases (CYPs) and NADPH-cytochrome P450 reductases (CPRs), each of which originated from 5 different sources, we selected two optimal betulinic acid-producing strains. Then, ERG1, ERG9, and HMG1 in the MVA module were overexpressed in the two strains, which dramatically increased betulinic acid production and resulted in a strain (YLJCC56) that exhibited the highest betulinic acid yield of 51.87 ± 2.77 mg/L. Then, we engineered the redox cofactor supply module by introducing NADPH- or NADH-generating enzymes and the acetyl-CoA generation module by directly overexpressing acetyl-CoA synthases or reinforcing the β-oxidation pathway, which further increased the total triterpenoid yield (the sum of the betulin, betulinic acid, betulinic aldehyde yields). Finally, we engineered these modules in combination, and the total triterpenoid yield reached 204.89 ± 11.56 mg/L (composed of 65.44% betulin, 23.71% betulinic acid and 10.85% betulinic aldehyde) in shake flask cultures.

Conclusions

Here, we systematically engineered Y. lipolytica and achieved, to the best of our knowledge, the highest betulinic acid and total triterpenoid yields reported in microbes. Our study provides a suitable reference for studies on heterologous exploitation of P450 enzymes and manipulation of triterpenoid production in Y. lipolytica.

Background

Betulinic acid is a pentacyclic lupane-type triterpenoid that widely exists in plants and possesses many pharmacological activities, such as antiviral (especially anti-HIV infection) [1, 2], antitumor (as a melanoma-specific cytotoxic agent) [3, 4] and other antimetabolic syndrome activities [5,6,7,8]. This compound can be acquired by purification of extracts from the outer bark of a variety of tree species [9]. Though this process has been optimized, problems such as low natural abundance, requirement of a large amount of land resources and long processing time are nonnegligible. The necessity for the development of an economical, ecofriendly and efficient production process has inspired researchers to use microorganisms for the production of betulinic acid [10].

The biosynthetic pathway of betulinic acid from glucose in Yarrowia lipolytica is shown in Fig. 1. In yeast, the biosynthetic precursors of betulinic acid are isopentenyl pyrophosphate (IPP) and dimethylallyl pyrophosphate (DMAPP), which are derived from acetyl-CoA via the mevalonic acid (MVA) pathway [11,12,13]. IPP then gives rise to higher order building blocks, namely, geranyl pyrophosphate (GPP) and farnesyl pyrophosphate (FPP), via the action of prenyltransferase enzymes. Two FPPs are condensed into squalene via a reaction catalyzed by squalene synthase (encoded by ERG9) and then converted to 2,3-oxidosqualene via a reaction catalyzed by squalene epoxidase (encoded by ERG1). The first step in the exogenous biosynthesis of betulinic acid from endogenous 2,3-oxidosqualene is catalyzed by lupeol synthase (LUS), which is followed by three sequential oxidations catalyzed by lupeol C-28 oxidases (cytochrome P450 monooxygenases, CYPs), leading to the generation of betulin, betulinic aldehyde, and betulinic acid. As increasing CYPs have been mined by bioinformatics analysis [14, 15], amino acid sequence alignment [16,17,18], and the combination approach of activation tagging and reverse genetic TILLING [19], heterologous betulinic acid production has been established in microbes [16, 20,21,22,23]. Li et al. demonstrated that attenuation of the competitive fatty acid biosynthesis pathway [23] and increasing the NADPH and oxygen supply [21] were beneficial for improvement of the production capability of betulinic acid in Saccharomyces cerevisiae. The titer of betulinic acid in yeast could be enhanced to 0.16 mg/L/OD600 by applying a multiple-approach strategy, which included using a novel CYP, employing a better yeast strain and upregulating the expression of key genes in the betulinic acid biosynthesis pathway [16]. Moreover, multicopy integration of AtLUP1 and CYP716Al1 with single-copy integration of ATR2 and overexpression of endogenous ERG1 and tHMG1 could improve betulinic acid production to ~ 28 mg/L in shake flask cultures in S. cerevisiae [22]. Furthermore, by coexpression of LUS and the fusion protein of CYP with CPR and multicopy integration of tHMG1 and ERG9, betulinic acid biosynthesis was achieved in Y. lipolytica and the yield was improved to 26.53 mg/L by substitution of glycerol for glucose as a carbon source [20]. Although the exogenous production pathway of betulinic acid has been introduced in yeast, no considerable yield was observed, and the engineering of heterologous production of betulinic acid was not performed as systematically as that of the production of other terpenoids.

Fig. 1
figure 1

Overview of the multimodular strategy for betulinic acid-related triterpenoid production in Y. lipolytica. All the genes that were engineered in this study are presented, and the betulinic acid biosynthesis pathway was divided into four modules: the red arrow represents the heterogenous CYP/CPR module (CYP is lupeol C-28 oxidase and CPR is NADPH-cytochrome P450 reductase. Both originate from the five sources). The yellow arrow represents the MVA module with 3 genes (ERG1/ERG9/HMG1) that are overexpressed separately or simultaneously. The green arrow represents the redox cofactor supply module with four introduced genes, namely, EMC, EMT, Rtme (encoding malic enzyme, which is responsible for NADPH generation) and Gapc (encoding glyceraldehyde-3-phosphate dehydrogenase, which is responsible for NADH generation). The blue arrow represents the acetyl-CoA generation module with seven endogenous genes overexpressed (ACL1 and ACL2, encoding ATP citrate lyase, which can directly increase acetyl-CoA levels; PXA1, MFE1, PEX10, POT1 and TGL3 in the β-oxidation pathway, which are responsible for the catabolism of fatty acids to generate acetyl-CoA)

Full size image

Multimodular metabolic engineering is an efficient strategy for addressing perturbances such as imbalance in pathway flux, module incompatibility or accumulation of intermediates caused by manipulation of endogenous genes and the introduction of heterologous pathways, and this method enables systematic manipulation of strains and metabolic pathways [24, 25]. This potential strategy has been widely implemented in Escherichia coli and S. cerevisiae for the production of value-added compounds. By partitioning the biosynthetic pathway of the diterpene taxadiene into two modules—an upstream methylerythritol-phosphate (MEP) module and a heterologous downstream terpenoid-forming module—and balancing the two modules by using various promoters and gene copy numbers in E. coli, production of the taxadiene was increased 15,000-fold compared to that in the parent strain and a yield of 1.02 ± 0.08 g/L achieved in fed-batch bioreactor fermentations [26]. Similarly, the production of miltiradiene (a kind of diterpenoid) in S. cerevisiae was improved from 0.5 to 12.5 mg/L by engineering MVA and heterologous synthase modules [27]. Recently, this strategy was also used to boost the production of a kind of triterpene—protopanaxadiol. By dividing the biosynthetic pathway of this compound into four modules—MVA, triterpene biosynthesis, sterol biosynthesis, and acetyl-CoA formation modules—and manipulating these modules simultaneously, protopanaxadiol production was enhanced to 66.55 mg/L/OD600 [28]. These studies demonstrated the potential of multimodular metabolic engineering for high-level production of terpenes and for addressing module incompatibility.

Yarrowia lipolytica, as a GRAS (generally recognized as safe) and oleaginous microorganism, has become a potential platform for the production of highly hydrophobic compounds [29, 30] because of the ability of this organism to accumulate large amounts of lipids that can store hydrophobic products [31, 32] and to supply sufficient intermediates and energy for terpene synthesis [33, 34]. In addition, the recent development of efficient and predictable synthetic biological tools in Y. lipolytica has contributed to the convenience of genetic manipulation in this organism [35,36,37,38]. Based on these factors, Y. lipolytica has been used to produce various hydrophobic terpenoids, such as linalool (monoterpene alcohol) [39], α-farnesene (sesquiterpene) [40], squalene (triterpene) [41], β-carotene, astaxanthin and lycopene (tetraterpene) [42,43,44,45]. In some cases, terpenoid production in Y. lipolytica was even higher than that in S. cerevisiae. For example, maximum production of β-carotene in Y. lipolytica was realized at 6.5 g/L [43], which demonstrated the potential of this strain for terpenoid production.

In this study, we used a multimodular metabolic engineering strategy to systematically engineer Y. lipolytica for production of the hydrophobic betulinic acid and related triterpenoids by dividing the biosynthetic pathway into four functional modules, including the heterogenous CYP/CPR, MVA, acetyl-CoA generation, and redox cofactor supply modules. By introducing exogenous genes or overexpressing endogenous key genes associated with these modules, we obtained a 51.87 ± 2.77 mg/L yield of betulinic acid and a 204.89 ± 11.56 mg/L yield of triterpenoids, both of which were the highest yields reported in shake flask cultures.

Results and discussion

Screening the CYP/CPR module for biosynthesis of betulinic acid in Y. lipolytica

Lupeol is the first heterologous intermediate metabolite from 2,3-oxidosqualene to betulinic acid, and this metabolite can be produced by the expression of LUSs. Subsequently, CYP716A subfamily members, which were identified as multifunctional triterpene C-28 oxidases [15], catalyze three sequential oxidations of lupeol to generate betulin, betulinic aldehyde and betulinic acid (Additional file 1: Figure S1). In this study, the RcLUS gene from Ricinus communis under the control of the TEFin promoter was integrated into the Ku70 locus (YALI0C08701 g) of Y. lipolytica to produce lupeol (Fig. 2a, Additional file 2: Figure S2). Because the native P450 enzymes of Y. lipolytica did not catalyze betulinic acid production (Additional file 3: Table S1), we tried to screen effective CYPs from plant species. Five CYP genes, namely, BPLO (CYP716A180) from Betula platyphylla [16], CrAO (CYP716AL1) from Catharanthus roseus [12], MtAO12 (CYP716A12) from Medicago truncatula [15], and VvAO15 (CYP716A15) and VvAO17 (CYP716A17) from Vitis vinifera [15], were expressed under the FBAin promoter and integrated into the Ku80 locus (YALI0E02068g), which generated strains YLJCC1, YLJCC7, YLJCC13, YLJCC19, YLJCC25. However, the desired product was not detected (Fig. 2a).

Fig. 2
figure 2

Screening of CYP/CPR for betulinic acid production. a Triterpenoid production by 30 strains that coexpressed LUS with CYP or LUS with CYP and CPR (both CYPs and CPRs originated from 5 different sources). b Triterpenoid production by strains that coexpressed LUS and fused CYP with CPR. N indicates that the encoded protein was located at the N-terminus of the fusion protein; C indicates that the encoded protein was located at the C-terminus of the fusion protein. L0 indicates that the two proteins were fused without any linker. L1, L2, L3 and L4 represent four different linkers, with amino acid sequences GGGS, GSG, GGGGS, and EAAAK, respectively. c Triterpenoid production of strains that coexpressed LUS and endoplasmic reticulum (ER)-targeted CYP and CPR. The ER-targeting sequence was fused to the C-terminus of both CYP and CPR. The red block in the first column of the gene forms represents the CYP/CPR module. Asterisk represents the strains used for subsequent optimization

Full size image

Previous studies have shown that in S. cerevisiae, recombinant plant CYPs work well with endogenous CPRs [21, 23]. However, our results and those of several other studies [20, 46, 47] indicated that in Y. lipolytica, plant CYPs might require the coexpression of heterogenous CPRs for exogenous terpenoid production. Accordingly, five CPR genes from different sources, namely, ATR2 from Arabidopsis thaliana, CrCPR from C. roseus, GuCPR from Glycyrrhiza uralensis, LjCPR from Lotus japonicus and MTR from M. truncatula, were expressed under the TEFin promoter and integrated into the IntF site (YALI0F3161413 to YALI0F3162449) [31]. Twenty-five strains were generated, and the fermentation results are shown in Fig. 2a (GC profiles and mass spectra of YLJCC6 are shown as an example in Additional file 4: Figure S3a and b). The results show that the combinations of CYPs from two sources (BPLO or VvAO15) and CPRs from three sources (ATR2, LjCPR, or MTR) could achieve the synthesis of betulinic acid (six strains were generated: YLJCC2, YLJCC5, YLJCC6, YLJCC20, YLJCC23, YLJCC24). Among the various CYPs, BPLO possesses relatively higher catalytic capacity for betulinic acid production in S. cerevisiae [16, 17]. In this study, BPLO also showed high catalytic capacity in Y. lipolytica. VvAO15 shared close evolutionary relationships and high homology with BPLO (Additional file 5: Figure S4a and b), which might be the reason why VvAO15 also performed well in Y. lipolytica. In addition, among the six strains, YLJCC5 and YLJCC6 showed the highest production of betulinic acid and total triterpenoids (Fig. 2a), which might be caused by the efficient interaction between BPLO and the related CPRs (that is LjCPR or MTR) [17, 48]. The titer of betulinic acid in YLJCC5 and YLJCC6 was 25.62 ± 2.20 mg/L and 32.33 ± 4.34 mg/L, respectively.

To further optimize the CYP/CPR module, we tried to fuse BPLO with LjCPR or MTR or to localize the proteins in the endoplasmic reticulum (ER) (Fig. 2b, c). We used two strategies for fusion of CYP and CPR. In the first strategy, the two proteins were fused without any linker and designated L0. In the second strategy, CYP and CPR were fused with four different linkers, designated L1, L2, L3 and L4. The amino acid sequences of these four linkers were GGGS, GSG, GGGGS, EAAAK, where the first three represent flexible linkers and the last one represents a rigid linker. The results showed that when CYP was fused to the N-terminus of the CPRs, the generated strains YLJCC31-35 and YLJCC41-45 showed a dramatic decline in the production of betulinic acid compared with the control strains, regardless of the strategy of fusion (with or without linker, or with any kind of linker). However, when CYP was fused to the C-terminus of the CPRs, the engineered strains barely produced any triterpenoids (Fig. 2b). In summary, the fusion of CYP with CPR could not further improve betulinic acid production. On the other hand, an ER-targeting peptide (KDEL) was fused to the C-terminus of the CYP (BPLO) and CPRs (LjCPR and MTR) [30] to generate strains YLJCC101 and YLJCC102, which coexpressed LUS, CYP-ER and CPR-ER (Fig. 2c). Although previous studies reported that the involvement of the ER could facilitate P450 function [49, 50], neither total triterpenoid production nor betulinic acid production was improved in the engineered strains, possibly because we used the full-length sequences of CYP and CPR, not truncated sequences, for heterologous expression, and these sequences already have signal peptides in the N-terminus of the proteins for membrane translocation [26, 51]. We fused a sfGFP [52] to the C-terminus of BPLO, LjCPR and MTR (without the ER-targeting peptide), and the confocal microscopy images confirmed our hypothesis. As shown in Additional file 6: Figure S5, the expressed proteins were already located in the membrane or organelle in Y. lipolytica due to the presence of the native signal peptide.

Finally, because no further improvement was observed upon fusion of CYP with CPR or upon addition of the ER-targeting peptide, the two best production strains from the above optimization, namely, YLJCC5 and YLJCC6 (Fig. 2a), were selected and subjected to further engineering.

Improving the production of triterpenoids by engineering the MVA module

Although the betulinic acid heterologous biosynthetic module was constructed in Y. lipolytica, the quantities of betulinic acid produced by strains YLJCC5 and YLJCC6 remained low. To increase the supply of the 2,3-oxidosqualene precursor, we engineered the MVA module with three structural genes, namely, ERG1, ERG9 and HMG1. These genes were overexpressed individually or in combination. Expression cassettes of these structural genes were integrated into the rDNA site of the Y. lipolytica genome (Fig. 3).

Fig. 3
figure 3

Engineering of the MVA module to enhance betulinic acid production. The genes ERG1, ERG9 and HMG1 in the MVA module were overexpressed individually or in combination. The light orange background represents strains generated from YLJCC5, and the reseda background represents strains generated from YLJCC6. Red and yellow blocks in the first column of the gene form represent the CYP/CPR and MVA modules, respectively. Asterisk represents the strains used in the final multimodular optimization

Full size image

First, squalene epoxidase (encoded by ERG1), which is responsible for the first oxygenation step in sterol biosynthesis and is recognized as the key enzyme in this pathway [53], was overexpressed under the control of the EXP promoter in strains YLJCC5 and YLJCC6, generating strains YLJCC51 and YLJCC54. As shown in Fig. 3, the total triterpenoid levels in YLJCC51 and YLJCC54 were slightly increased. Then, we further overexpressed squalene synthase (encoded by ERG9), which converts FPP to squalene and represents the branching point of synthetic pathways of triterpenoids, sterols and sesquiterpenoids. The resulting strains, namely, YLJCC52 and YLJCC55, exhibited a further increase in betulinic acid and total triterpenoid production (Fig. 3). Finally, with the addition of HMG1, which encodes a rate-limiting enzyme in the MVA pathway, the production of betulinic acid in the generated strains YLJCC53 and YLJCC56 improved 1.25- and 1.60-fold, respectively, compared with YLJCC5 and YLJCC6, and the total triterpenoid production increased by 1.34-fold and 1.91-fold, respectively (Fig. 3).

Previously, overexpression of ERG1 in S. cerevisiae increased the levels of the triterpenoid protopanaxadiol [54]. It has been demonstrated that overexpression of ERG9 enhanced triterpenoid accumulation in plants [55,56,57] and β-amyrin production in S. cerevisiae [58]. Overexpression of tHMG1 (encoding the N-terminal truncated version of HMG1) is a commonly used approach to improve the production of terpenoids because this method can overcome the tight regulation of HMG1 in host cells [40, 42, 59,60,61]. However, some data indicated that overexpression of homologous HMG1 was more effective than overexpression of tHMG1 for production of the monoterpene limonene in Y. lipolytica [62]. Recently, overexpression of endogenous HMG1 improved squalene synthesis 3.2-fold in Y. lipolytica [41]. Hence, we chose to overexpress HMG1 rather than tHMG1 in Y. lipolytica. Our results suggested that overexpression of ERG1 and ERG9 only slightly enhanced betulinic acid production in Y. lipolytica. Further overexpression of HMG1 led to a marked improvement in the yield of the product.

Improving triterpenoid production via engineering of the redox cofactor supply module

One of the key factors required for triterpenoid biosynthesis is sufficient redox supply [21]. Here, we tried to improve the NADPH or NADH supply by overexpression of malic enzyme or glyceraldehyde-3-phosphate (G3P) dehydrogenase, respectively. EMC, EMT and Rtme encoding malic enzyme and Gapc encoding G3P dehydrogenase (shown in Fig. 1) were integrated into the IntE site (YALI0E3673678–YALI0E3672531) [31] of YLJCC5 and YLJCC6 (Fig. 4), respectively, to generate strains YLJCC63-YLJCC70. The fermentation results are shown in Fig. 4. Compared with YLJCC5, the production of total triterpenoids and betulinic acid in YLJCC64 and YLJCC65 (harboring the introduced genes EMT and Rtme) increased slightly, although this increasing trend was not observed in YLJCC68 and YLJCC69 compared with YLJCC6. In contrast, when Gapc, which is responsible for NADH generation, was overexpressed, triterpenoid production decreased (YLJCC66). These results demonstrated that manipulation of the NADPH supply was more beneficial than manipulation of the NADH supply for betulinic acid production, which might be because NADPH is the direct redox supplier for endogenous MVA and the heterogenous CYP/CPR system. A previous study also demonstrated that improvement of the NADPH supply was an efficient method for improvement of betulinic acid production by introduction of mBDH1, which can increase cellular NADPH level with acetoin supplementation [21].

Fig. 4
figure 4

Engineering of the redox cofactor supply module to improve triterpenoid production. The EMC, EMT and Rtme genes that encode malic enzymes are responsible for NADPH generation. The gene Gapc, encoding glyceraldehyde-3-phosphate dehydrogenase, is responsible for NADH generation. The light orange background represents strains generated from YLJCC5, and the reseda background represents strains generated from YLJCC6. Red and green blocks in the first column of the gene form represent the CYP/CPR and redox cofactor supply modules, respectively. Asterisk represents the strains that exhibited increased triterpenoid production, and similar engineering was used in the final multimodular optimization

Full size image

Improving triterpenoid synthesis via engineering of the acetyl-CoA generation module

Acetyl-CoA is the biosynthetic precursor of lipids, terpenes and sterols. Y. lipolytica, which has been characterized as a nonconventional yeast, is able to accumulate large amounts of intracellular lipids, which seems to be a primary competition pathway for terpene production. To increase the carbon flux into the MVA pathway, we tried to recover acetyl-CoA from lipids by enhancing β-oxidation of fatty acids via overexpression of four structural genes, namely, the triacylglycerol lipase gene (TGL3), long-chain fatty acid transporter gene (PXA1), multifunctional β-oxidation enzyme gene (MFE1), 3-ketoacyl-CoA thiolase gene (POT1), and one regulatory gene, namely, the regulatory gene PEX10, which is necessary for correct peroxisomal biogenesis and morphology for the β-oxidation process (shown in Fig. 1). Moreover, we also overexpressed ATP citrate lyase (ACL1 and ACL2) to directly increase the synthesis of acetyl-CoA. In summary, the engineered strains YLJCC71-YLJCC84 were constructed, and the fermentation results are shown in Fig. 5.

Fig. 5
figure 5

Engineering of the acetyl-CoA generation module to increase triterpenoid production. ACL1 and ACL2 encode ATP citrate lyase, which can directly increase acetyl-CoA levels; PXA1, MFE1, PEX10, POT1, and TGL3 in the β-oxidation pathway are responsible for the catabolism of fatty acids to acetyl-CoA. The light orange background represents strains generated from YLJCC5, and the reseda background represents strains generated from YLJCC6. Red and blue blocks in the first column of gene form represent the CYP/CPR and acetyl-CoA generation modules, respectively. Asterisk represents the strains that exhibited increased triterpenoid production, and similar engineering was used in the final multimodular optimization

Full size image

Overexpression of the endogenous ATP citrate lyase gene (ylACL1) has been previously shown to increase the level of cytosolic acetyl-CoA 1.5-fold in Y. lipolytica [41, 63]. However, in this study, overexpression of ACL1 (YLJCC71) failed to improve the production of triterpenes, while overexpression of ACL2 (YLJCC72) increased the titers of total triterpenoids and betulinic acid 1.23- and 1.30-fold, respectively (Fig. 5). On the other hand, overexpression of all the genes in the β-oxidation pathway except PEX10 in YLJCC5 led to increased triterpene production. Among the generated strains, YLJCC76 (overexpressing POT1) exhibited the most significant enhancement of production with a triterpene yield of 66.77 ± 2.79 mg/L. Similarly, this elevated production did not occur in the strains generated from YLJCC6, which might be caused by module incompatibility.

In most cases, studies have focused on elimination of the competition pathway to reduce carbon loss. For example, α-santalene production increased upon knocking out of citrate synthase or malate synthase, which convert acetyl-CoA to citrate or C4 organic acids, respectively [64]. Fatty acid biosynthesis is also considered to be a competitive pathway for terpenoid synthesis because these pathways use the same precursor—acetyl-CoA. To reduce competition, the endogenous promoters of the key genes in the fatty acid pathway were replaced with relatively weak promoters, and the betulinic acid yield increased to 1.92 mg/L/OD in S. cerevisiae [23]. In this study, we applied a different strategy to address competition by enhancing fatty acid β-oxidation to recover acetyl-CoA. The results indicated that the strategy was effective for betulinic acid-related triterpenoid production in Y. lipolytica.

Further increasing the production of triterpenoids by multimodular optimization

By partitioning the betulinic acid biosynthetic pathway into four modules and engineering these modules separately, we identified the effective genes within each module. To boost triterpenoid production, we manipulated these modules in combination. Strains YLJCC53 and YLJCC56 (Fig. 3), which possessed the optimized heterogenous CYP/CPR module and MVA module, were used for the engineering process. EMT and Rtme in the redox cofactor supply module (Fig. 4) or ACL2, PXA1, MFE1, POT1, and TGL3 in the acetyl-CoA generation module (Fig. 5), which showed the capacity to increase triterpenoid production, were further overexpressed, generating strains YLJCC85-YLJCC98. The fermentation results for these strains are shown in Fig. 6. The strains generated from YLJCC53 showed observable improvement in betulin production, with a 1.8- to 2.9-fold increase in total triterpenoid titers compared with YLJCC53. Among these strains, the highest production was observed for strain YLJCC91, which exhibited a total triterpenoid yield of 204.89 ± 11.56 mg/L; the triterpenoids were composed of 65.44% betulin, 23.71% betulinic acid and 10.85% betulinic aldehyde. Although this dramatic increase was not observed in strains generated from YLJCC56, which might be due to the incompatibility of elevated redox or acetyl-CoA supply with the BPLO/MTR module, YLJCC56 exhibited the best betulinic acid production. Previously, the highest betulinic acid and triterpenoid yields in yeast in shake flask cultures were 28 mg/L and 147 mg/L, respectively [22]. Here, via systematic engineering, we obtained two high-yield strains. The first strain was YLJCC56 (Fig. 3), which produced 51.87 ± 2.77 mg/L betulinic acid. The second strain was YLJCC91 (Fig. 6b), which produced 204.89 ± 11.56 mg/L total triterpenoids. To the best of our knowledge, both titers are the highest in microbes in shake flask cultures to date. However, the yield of betulinic acid and total triterpenoids in this study still remained far from the theoretical yield, thus, more studies need to be carried out to improve the heterologous production of those valuable compounds.

Fig. 6
figure 6

Multimodular optimization to further improve triterpenoid production. a Triterpenoid production by strains in which the CYP/CRP, MVA and redox cofactor supply modules were optimized in combination. b Triterpenoid production by strains in which the CYP/CPR, MVA and acetyl-CoA generation modules were optimized in combination. The light orange background represents strains generated from YLJCC53, and the reseda background represents strains generated from YLJCC56. Red, yellow, green and blue blocks in the first column of the gene forms represent the CYP/CPR, MVA, redox cofactor supply, and acetyl-CoA generation modules, respectively. Asterisk represents the strain exhibiting the highest total triterpenoid production in this study

Full size image

To identify the impact of the engineering on cell growth and metabolism, we determined the dry cell weights (DCWs) and the levels of squalene and lipids in the parent strain 201249 and engineered strains marked with * (YLJCC5 and YLJCC6 in Fig. 2a; YLJCC53 and YLJCC56 in Fig. 3; YLJCC64 and YLJCC65 in Fig. 4; YLJCC72, YLJCC73, YLJCC74, YLJCC76, and YLJCC77 in Fig. 5; YLJCC91 in Fig. 6b). The results are presented in Fig. 7. Although the lipid levels and DCWs showed slight differences among the engineered strains, the change pattern was not directly correlated with triterpenoid production (the low biomass of the parent strain 201249 might result from its triple auxotroph characteristics [65]). However, the squalene levels showed a strong correlation with MVA module manipulation. After upregulation of the MVA module, the levels of squalene in YLJCC53 and YLJCC56 increased 1.54- and 1.68-fold, respectively, compared with those in the parent strain 201249. Based on YLJCC53, improvement of acetyl-CoA generation by overexpression of MFE1 (YLJCC91) further increased squalene accumulation 1.72-fold, which might represent an increased carbon flux to the MVA pathway and terpenoid synthesis.

Fig. 7
figure 7

Triterpenoid, lipid, DCW, and squalene production in the parent strain and engineered strains marked with Asterisk. Red, yellow, green and blue blocks in the first column of the gene form represent the CYP/CPR, MVA, redox cofactor supply, and acetyl-CoA generation modules, respectively

Full size image

Conclusions

Herein, we systematically engineered Y. lipolytica using a multimodular strategy to increase triterpenoid production. By matching CYPs and CPRs, both of which originated from five different sources, and coexpressing these proteins with LUS, we established betulinic acid synthesis in Y. lipolytica. Then, the MVA module, acetyl-CoA generation module, and redox cofactor supply module were optimized separately. The betulinic acid yield improved to 51.87 ± 2.77 mg/L when ERG1, ERG9, and HMG1 in the MVA module were overexpressed with RcLUS, BPLO and MTR. Finally, the optimal modules were integrated, and a 204.89 ± 11.56 mg/L total triterpenoid yield was realized when RcLUS, BPLO and LjCPR were coexpressed with the upregulated MVA (overexpressing ERG1, ERG9 and HMG1) and acetyl-CoA generation module (overexpressing MFE1). Both betulinic acid and total triterpenoid production achieved the highest reported yields in microbes in shake flask cultures. This study can serve as a foundation for the construction of efficient triterpene production strains of Y. lipolytica.

Methods

Strains, medium and culture conditions

TransT1 E. coli was used for plasmid construction and propagation, and cell growth was performed at 37 °C in Luria–Bertani (LB) broth with 100 mg/L ampicillin or 50 mg/L kanamycin and shaking at 220 rpm. The incipient Y. lipolytica strain ATCC 201249 was kindly provided by Professor Ying-Jin Yuan (Tianjin University, China). All the Y. lipolytica strains in this article (shown in Table 1) were cultivated at 28 °C with shaking at 250 rpm. Yeast extract-peptone-dextrose (YPD) medium, consisting of 20 g/L glucose, 20 g/L peptone, and 10 g/L yeast extract, was used for transformation, activation and preculture.

Table 1 Strains used in this work
Full size table

The proper auxotrophic medium used to grow transformants consisted of 20 g/L glucose H2O, 5 g/L (NH4)2SO4, 1.7 g/L yeast nitrogen base without amino acids, 2 g/L complete supplement mixture (CSM) lacking uracil (SC-Ura) or leucine (SC-Leu). YPD-hum medium containing 200 mg/L hygromycin was also used to screen the recombined strains, and 20 g/L agar was added for preparation of solid plates. Rich YPD medium, consisting of 50 g/L glucose, 20 g/L peptone, and 10 g/L yeast extract, was used for fermentation.

The seed cultures were first prepared in 25-mL polypropylene tubes containing 5 mL of YPD medium. After 24 h, the seed cultures were inoculated into a 250-mL shake flask containing 50 mL of rich YPD medium with an initial OD600 of 0.2 for 168 h. A 1/5 volumetric ratio of isopropyl myristate (IPM) was added to the shake flasks at 120 h to extract the products. All the titers presented in this research were based on the volume of the fermentation media. All assays were performed in at least triplicate (at least three independent assays for three different clones).

Plasmids and strain construction

All the plasmids used in this study are listed in Table 2. The integration plasmids pIntK, pIntF, pIntE, pK8FB, pLD01, pLD02in, pLD03in, prDNA-hph, pt1-R, pt2-R, and pt3-R were kindly provided by Professor Ying-Jin Yuan (Tianjin University, China). Five CYP-coding genes (BPLO from Betula platyphylla, CrAO from Catharanthus roseus, MtAO12 from Medicago truncatula, VvAO15 and VvAO17 from Vitis vinifera), five CPR-coding genes (ATR2 from Arabidopsis thaliana, CrCPR from Catharanthus roseus, GuCPR from Glycyrrhiza uralensis, LjCPR from Lotus japonicas and MTR from Medicago truncatula), four genes coding enzymes related to NADH or NADPH generation (EMC from Mucor circinelloides, EMT from Mortierella alpina, Gapc from Kluyveromyces lactis and Rtme from Rhodotorula toruloides NP11) and sfGFP were codon optimized and synthesized by GenScript (Nanjing, China) and listed in Additional file 7: Table S3. Three structural genes of MVA pathway (ERG1, ERG9 and HMG1), the long-chain fatty acid transporter coded gene PXA1, 3-ketoacyl-CoA thiolase coded gene POT1, multifunctional beta-oxidation enzyme coded gene MFE1, transcription factor coded gene PEX10, ATP citrate lyase subunit a coded gene ACL1, ATP citrate lyase subunit b coded gene ACL2, and triacyl -glycerol lipase coded gene TGL3 were amplified from genomic DNA of Y. lipolytica. All the genes were amplified with the proper sticky ends and incorporated into the corresponding integration plasmids by Golden Gate assembly or Gibson assembly.

Table 2 Plasmids used in this work
Full size table

Fusion protein-encoding genes were obtained by overlap extension polymerase chain reaction (PCR). All the primers used to construct fusion protein-encoding genes were synthesized by Genewiz Inc. (China) and are listed in Additional file 7: Table S2. We applied two types of fusion between CYP and CPR. In the first case, CYPs were fused with CPRs directly (marked as L0) and generated pK8FB-BPLO-LjCPR, pK8FB-BPLO-MTR, pK8FB-LjCPR-BPLO, and pK8FB-MTR-BPLO. In the other case, CYPs and CPRs were fused with different linkers, which contained linker1 (L1) “GGAGGAGGATCT”, linker2 (L2) “GGATCTGGA, linker3 (L3) “GGAGGAGGAGGATCT”, and linker4 (L4) “GAGGC CGCCGCCAAG” (the first three represents flexible linkers and the last one represents rigid linker), and generated pK8FB-BPLO-L1-LjCPR, pK8FB-BPLO-L2-LjCPR, pK8FB-BPLO-L3-LjCPR, pK8FB-BPLO-L4-LjCPR, pK8FB-LjCPR-L1-BPLO, pK8FB-LjCPR-L2-BPLO, pK8FB-LjCPR-L3-BPLO, pK8FB-BPLO-L1-MTR, pK8FB-BPLO-L2-MTR, pK8FB-BPLO-L3-MTR, pK8FB-BPLO-L4-MTR, pK8FB-MTR-L2-BPLO, pK8FB-MTR-L3-BPLO and pK8FB-MTR-L4-BPLO. All the fragments were amplified with the proper sticky ends and incorporated into the integration plasmid K8FB by Golden Gate assembly. The ER-targeting peptide (KDEL) and sfGFP were fused to the C-terminus of CYP and CPR to generate pK8FB-BPLO-ER, pIntE-LjCPR-ER, pIntE-MTR-ER, and pK8FB-BPLO-sfGFP, pK8FB-LjCPR-sfGFP, and pK8FB-MTR-sfGFP.

All the integrative fragments were transformed into Y. lipolytica using the Zymogen Frozen EZ Yeast Transformation Kit II (Zymo Research Corporation). Two micrograms of linearized DNA were used in the transformation reaction, and then, the cells were harvested and plated on selective solid plates. The Ura3 selection maker was recycled by incubating the transformants in 4 mL of YPD for 4 days and diluting 10,000 times on SC solid medium containing 1.2 mg/mL 5-fluoroorotic acid for 2 days. The genotypes of all the strains were confirmed by colony PCR using KOD FX DNA polymerase (Toyobo Co., Ltd.; Shanghai, China). The restriction endonucleases and ligase used for plasmid construction were purchased from NEB. Plasmid extraction and DNA fragment purification were performed using the TIANprep Mini Plasmid Kit and TIANgal Midi Purification Kit, respectively (TIANGEN BIOTECH, China). Gibson assembly was carried out using the ClonExpressII One Step Cloning Kit (Vazyme, China).

Lupeol, betulinic acid, betulin, and betulinic aldehyde extraction and analysis

The IPM phase was collected by centrifugation of 10 mL of fermentation broth mixture (10 min, 5000 rpm) followed by filtration with a 0.22-μm film. Then, 50 μL of the resulting extract was trimethylsilylated by mixing with 150 μL of N-methyl-N-(trimethylsilyl) trifluoracetamide (MSTFA) at 80 °C for 1 h and filtering with a 0.22-μm film. Appropriate concentrations of lupeol, betulin (Aladdin Industrial Corporation, US) and betulinic acid (MedChemExpress, US) were dissolved in IPM to prepare standard curves via GC–MS analysis. The methods used for betulinic acid detection were based on previous studies with some modification [15,16,17, 20]. A Rts-5 MS (30 m × 0.25 mm × 0.5 μm) GC column was used on a SHIMADZU GCMS-QP2020 system. A 1-μL sample was injected by AOC-20i+s with a split ratio of 10. The temperature of the injector was set at 250 °C. The initial oven temperature was maintained at 80 °C for 2 min, and then, the temperature was increased to 300 °C at a rate of 20 °C/min and held for 17 min. The GC system was operated under a constant pressure of 91 kPa. For the MS system, the temperature of the ion source was 230 °C, and full scan mode was used from 50 to 700 m/z.

Squalene extraction and analysis

Extraction and analysis of squalene were performed according to Guo XJ et al. [66] with some modification. Cells were harvested by centrifugation and washed 3 times with distilled water. Then, 1 mL of 3 M HCl was added for resuspension, followed by boiling in a 100 °C water bath for 5 min. The cell debris was washed with distilled water until the pH was neutral. NaOH solution was used to neutralize the residual HCl. Then, the cell pellet was resuspended with 1.5 M NaOH-methanol solution and incubated at 60 °C for 4 h. Then, 500 μL of hexyl hydride was added to the tubes, and the mixture was vortexed for 20 min. After centrifugation, the hexyl hydride layer was transferred to new tubes. This process was repeated once. The hexyl hydride layer was vacuum-dried, resuspended in 100 μL of MSTFA and incubated at 30 °C for 2 h. After filtering with a 0.22-μm film, the sample was used for GC–MS detection. A Rts-5 MS (30 m × 0.25 mm × 0.5 μm) GC column was used on the SHIMADZU GCMS-QP2020 system. A 1-μL sample was injected by AOC-20i + s with a split ratio of 5. The temperature of the injector was set at 280 °C. The initial oven temperature was maintained at 70 °C for 1 min, and then, the temperature was increased to 250 °C at a rate of 30 °C/min and held for 2 min. Then, the temperature was increased to 280 °C at a rate of 20 °C/min and held for 20 min. The GC system was operated under a constant pressure of 91 kPa. For the MS system, the temperature of the ion source was 220 °C, and full scan mode was used from 50 to 600 m/z.

Lipid extraction and analysis

Extraction and analysis of lipids were performed according to Wang et al. [67]. Cells were harvested by centrifugation and washed 3 times with distilled water. One milliliter of methanol containing 3 M HCl and 100 mg/L heptadecanoic acid (as an internal standard) was added to the tubes and vortexed for mixing. Then, these tubes were incubated at 70 °C for 3 h and inverted several times for 40 min. Then, these tubes were cooled naturally to the indoor temperature. Then, sodium chloride was added into the solutions until saturation, and the solutions were shaken for 1 min. Then, 500 μL of hexyl hydride was added into these tubes with shaking for 1 min. The hexyl hydride phase was collected by centrifugation, filtered with a 0.22-μm film and used for the detection of lipid production by GC–MS. A Rts-5 MS (30 m × 0.25 mm × 0.5 μm) GC column was used on a SHIMADZU GCMS-QP2020 system. A 1-μL sample was injected by AOC-20i+s with a split ratio of 50. The temperature of the injector was set at 260 °C. The initial oven temperature was maintained at 70 °C for 2 min, and then, the temperature was increased to 290 °C at a rate of 8 °C/min and held for 6 min. The GC system was operated under a constant pressure of 91 kPa. For the MS system, the temperature of the ion source was 220 °C, and full scan mode was used from 50 to 500 m/z.

Phylogenetic tree and homology analyses

Full-length amino acid sequences were obtained from GenBank and aligned by ClustalW. A phylogenetic tree was constructed using the neighbor-joining method with 1000 replicates using MEGA6 software. Amino acid sequence homology analysis was performed with BLASTP.

Confocal laser scanning microscopy

Cell imaging was performed using UltraView Vox (PerkingElmer, America) after washing with ddH2O. Excitation/emission wavelengths were 488 nm/530 nm for GFP.

References

  1. Fujioka T, Kashiwada Y, Kilkuskie RE, Cosentino LM, Ballas LM, Jiang JB, et al. Anti-AIDS agents, 11. betulinic acid and platanic acid as anti-HIV principles from Syzigium claviflorum, and the anti-HIV activity of structurally related triterpenoids. J Nat Prod. 1994;57(2):243–7.

    Article  CAS  PubMed  Google Scholar 

  2. Aiken C, Chen CH. Betulinic acid derivatives as HIV-1 antivirals. Trends Mol Med. 2005;11(1):31–6.

    Article  CAS  PubMed  Google Scholar 

  3. Pisha E, Chai H, Lee I-S, Chagwedera TE, Farnsworth NR, Cordell GA, et al. Discovery of betulinic acid as a selective inhibitor of human melanoma that functions by induction of apoptosis. Nat Med. 1995;1(10):1046–51.

    Article  CAS  PubMed  Google Scholar 

  4. Cichewicz RH, Kouzi SA. Chemistry, biological activity, and chemotherapeutic potential of betulinic acid for the prevention and treatment of cancer and HIV infection. Med Res Rev. 2004;24(1):90–114.

    Article  CAS  PubMed  Google Scholar 

  5. Silva FS, Oliveira PJ, Duarte MF. Oleanolic, ursolic, and betulinic acids as food supplements or pharmaceutical agents for type 2 diabetes: promise or illusion? J Agric Food Chem. 2016;64(15):2991–3008.

    Article  CAS  PubMed  Google Scholar 

  6. Gautam R, Jachak SM. Recent developments in anti-inflammatory natural products. Med Res Rev. 2009;29(5):767–820.

    Article  CAS  PubMed  Google Scholar 

  7. Kim J, Lee YS, Kim CS, Kim JS. Betulinic acid has an inhibitory effect on pancreatic lipase and induces adipocyte lipolysis. Phytother Res. 2012;26(7):1103–6.

    Article  CAS  PubMed  Google Scholar 

  8. Rios JL, Manez S. New pharmacological opportunities for betulinic acid. Planta Med. 2018;84(1):8–19.

    Article  CAS  PubMed  Google Scholar 

  9. Yogeeswari P, Sriram D. Betulinic acid and its derivatives: a review on their biological properties. Curr Med Chem. 2005;12(6):657–66.

    Article  CAS  PubMed  Google Scholar 

  10. Ajikumar PK, Tyo K, Carlsen S, Mucha O, Phon TH, Stephanopoulos G. Terpenoids: opportunities for biosynthesis of natural product drugs using engineered microorganisms. Mol Pharm. 2008;5(2):167–90.

    Article  CAS  PubMed  Google Scholar 

  11. Sacchettini JC, Poulter CD. Creating isoprenoid diversity. Science. 1997;277(5333):1788–9.

    Article  CAS  PubMed  Google Scholar 

  12. Huang L, Li J, Ye H, Li C, Wang H, Liu B, et al. Molecular characterization of the pentacyclic triterpenoid biosynthetic pathway in Catharanthus roseus. Planta. 2012;236(5):1571–81.

    Article  CAS  PubMed  Google Scholar 

  13. Miziorko HM. Enzymes of the mevalonate pathway of isoprenoid biosynthesis. Arch Biochem Biophys. 2011;505(2):131–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Yu L, Tang W, He W, Ma X, Vasseur L, Baxter SW, et al. Characterization and expression of the cytochrome P450 gene family in diamondback moth, Plutella xylostella (L.). Sci Rep. 2015;5:8952.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Fukushima EO, Seki H, Ohyama K, Ono E, Umemoto N, Mizutani M, et al. CYP716A subfamily members are multifunctional oxidases in triterpenoid biosynthesis. Plant Cell Physiol. 2011;52(12):2050–61.

    Article  CAS  PubMed  Google Scholar 

  16. Zhou C, Li J, Li C, Zhang Y. Improvement of betulinic acid biosynthesis in yeast employing multiple strategies. BMC Biotechnol. 2016;16(1):59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Suzuki H, Fukushima EO, Umemoto N, Ohyama K, Seki H, Muranaka T. Comparative analysis of CYP716A subfamily enzymes for the heterologous production of C-28 oxidized triterpenoids in transgenic yeast. Plant Biotechnol. 2018;35(2):131–9.

    Article  CAS  Google Scholar 

  18. Miettinen K, Pollier J, Buyst D, Arendt P, Csuk R, Sommerwerk S, et al. The ancient CYP716 family is a major contributor to the diversification of eudicot triterpenoid biosynthesis. Nat Commun. 2017;8:14153.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Carelli M, Biazzi E, Panara F, Tava A, Scaramelli L, Porceddu A, et al. Medicago truncatula CYP716A12 is a multifunctional oxidase involved in the biosynthesis of hemolytic saponins. Plant Cell. 2011;23(8):3070–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Sun J, Zhang C, Nan W, Li D, Ke D, Lu W. Glycerol improves heterologous biosynthesis of betulinic acid in engineered Yarrowia lipolytica. Chem Eng Sci. 2018. https://doi.org/10.1016/j.ces.2018.10.052.

    Article  Google Scholar 

  21. Li J, Zhang Y. Modulating betulinic acid production in Saccharomyces cerevisiae by managing the intracellular supplies of the co-factor NADPH and oxygen. J Biosci Bioeng. 2015;119(1):77–81.

    Article  CAS  PubMed  Google Scholar 

  22. Czarnotta E, Dianat M, Korf M, Granica F, Merz J, Maury J, et al. Fermentation and purification strategies for the production of betulinic acid and its lupane-type precursors in Saccharomyces cerevisiae. Biotechnol Bioeng. 2017;114(11):2528–38.

    Article  CAS  PubMed  Google Scholar 

  23. Li J, Zhang Y. Increase of betulinic acid production in Saccharomyces cerevisiae by balancing fatty acids and betulinic acid forming pathways. Appl Microbiol Biotechnol. 2014;98(7):3081–9.

    Article  CAS  PubMed  Google Scholar 

  24. Biggs BW, De Paepe B, Santos CNS, De Mey M, Ajikumar PK. Multivariate modular metabolic engineering for pathway and strain optimization. Curr Opin Biotechnol. 2014;29:156–62.

    Article  CAS  PubMed  Google Scholar 

  25. Li J, Meng H, Wang Y. Synbiological systems for complex natural products biosynthesis. Synth Syst Biotechnol. 2016;1(4):221–9.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Ajikumar PK, Xiao WH, Tyo KE, Wang Y, Simeon F, Leonard E, et al. Isoprenoid pathway optimization for taxol precursor overproduction in Escherichia coli. Science. 2010;330(6000):70–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zhou YJ, Gao W, Rong Q, Jin G, Chu H, Liu W, et al. Modular pathway engineering of diterpenoid synthases and the mevalonic acid pathway for miltiradiene production. J Am Chem Soc. 2012;134(6):3234–41.

    Article  CAS  PubMed  Google Scholar 

  28. Zhao F, Bai P, Nan W, Li D, Zhang C, Lu C, et al. A modular engineering strategy for high-level production of protopanaxadiol from ethanol by Saccharomyces cerevisiae. AIChE J. 2018. https://doi.org/10.1002/aic.16502.

    Article  Google Scholar 

  29. Sabirova JS, Haddouche R, Van Bogaert I, Mulaa F, Verstraete W, Timmis K, et al. The ‘LipoYeasts’ project: using the oleaginous yeast Yarrowia lipolytica in combination with specific bacterial genes for the bioconversion of lipids, fats and oils into high-value products. Microb Biotechnol. 2011;4(1):47–54.

    Article  CAS  PubMed  Google Scholar 

  30. Xu P, Qiao K, Ahn WS, Stephanopoulos G. Engineering Yarrowia lipolytica as a platform for synthesis of drop-in transportation fuels and oleochemicals. Proc Natl Acad Sci USA. 2016;113(39):10848–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Matthaus F, Ketelhot M, Gatter M, Barth G. Production of lycopene in the non-carotenoid-producing yeast Yarrowia lipolytica. Appl Environ Microbiol. 2014;80(5):1660–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Friedlander J, Tsakraklides V, Kamineni A, Greenhagen EH, Consiglio AL, MacEwen K, et al. Engineering of a high lipid producing Yarrowia lipolytica strain. Biotechnol Biofuels. 2016;9(1):77.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  33. Blank LM, Lehmbeck F, Sauer U. Metabolic-flux and network analysis in fourteen hemiascomycetous yeasts. FEMS Yeast Res. 2005;5(6–7):545–58.

    Article  CAS  PubMed  Google Scholar 

  34. Wasylenko TM, Ahn WS, Stephanopoulos G. The oxidative pentose phosphate pathway is the primary source of NADPH for lipid overproduction from glucose in Yarrowia lipolytica. Metab Eng. 2015;30:27–39.

    Article  CAS  PubMed  Google Scholar 

  35. Abdel-Mawgoud AM, Markham KA, Palmer CM, Liu N, Stephanopoulos G, Alper HS. Metabolic engineering in the host Yarrowia lipolytica. Metab Eng. 2018;50:192–208.

    Article  CAS  PubMed  Google Scholar 

  36. Zhang JL, Peng YZ, Liu D, Liu H, Cao YX, Li BZ, et al. Gene repression via multiplex gRNA strategy in Y. lipolytica. Microb Cell Fact. 2018;17(1):62.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  37. Bredeweg EL, Pomraning KR, Dai Z, Nielsen J, Kerkhoven EJ, Baker SE. A molecular genetic toolbox for Yarrowia lipolytica. Biotechnol Biofuels. 2017;10(1):2.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  38. Yu AQ, Pratomo N, Ng TK, Ling H, Cho HS, Leong SS, et al. Genetic engineering of an unconventional yeast for renewable biofuel and biochemical production. J Vis Exp. 2016;2016(115):54371.

    Google Scholar 

  39. Cao X, Wei LJ, Lin JY, Hua Q. Enhancing linalool production by engineering oleaginous yeast Yarrowia lipolytica. Bioresour Technol. 2017;245(Pt B):1641–4.

    Article  CAS  PubMed  Google Scholar 

  40. Yang X, Nambou K, Wei L, Hua Q. Heterologous production of α-farnesene in metabolically engineered strains of Yarrowia lipolytica. Bioresour Technol. 2016;216:1040–8.

    Article  CAS  PubMed  Google Scholar 

  41. Huang YY, Jian XX, Lv YB, Nian KQ, Gao Q, Chen J, et al. Enhanced squalene biosynthesis in Yarrowia lipolytica based on metabolically engineered acetyl-CoA metabolism. J Biotechnol. 2018;281:106–14.

    Article  CAS  PubMed  Google Scholar 

  42. Gao S, Tong Y, Zhu L, Ge M, Zhang Y, Chen D, et al. Iterative integration of multiple-copy pathway genes in Yarrowia lipolytica for heterologous β-carotene production. Metab Eng. 2017;41:192–201.

    Article  CAS  PubMed  Google Scholar 

  43. Larroude M, Celinska E, Back A, Thomas S, Nicaud JM, Ledesma-Amaro R. A synthetic biology approach to transform Yarrowia lipolytica into a competitive biotechnological producer of β-carotene. Biotechnol Bioeng. 2018;115(2):464–72.

    Article  CAS  PubMed  Google Scholar 

  44. Kildegaard KR, Adiego-Perez B, Domenech Belda D, Khangura JK, Holkenbrink C, Borodina I. Engineering of Yarrowia lipolytica for production of astaxanthin. Synth Syst Biotechnol. 2017;2(4):287–94.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Nambou K, Jian X, Zhang X, Wei L, Lou J, Madzak C, et al. Flux balance analysis inspired bioprocess upgrading for lycopene production by a metabolically engineered strain of Yarrowia lipolytica. Metabolites. 2015;5(4):794–813.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Guo X, Sun J, Li D, Lu W. Heterologous biosynthesis of (+)-nootkatone in unconventional yeast Yarrowia lipolytica. Biochem Eng J. 2018;137:125–31.

    Article  CAS  Google Scholar 

  47. Ma YR, Wang KF, Wang WJ, Ding Y, Shi TQ, Huang H, et al. Advances in the metabolic engineering of Yarrowia lipolytica for the production of terpenoids. Bioresour Technol. 2019;281:449–56.

    Article  CAS  PubMed  Google Scholar 

  48. Schuler MA. Plant cytochrome P450 monooxygenases. Crit Rev Plant Sci. 2010;15(3):235–84.

    Article  Google Scholar 

  49. Arendt P, Miettinen K, Pollier J, De Rycke R, Callewaert N, Goossens A. An endoplasmic reticulum-engineered yeast platform for overproduction of triterpenoids. Metab Eng. 2017;40:165–75.

    Article  CAS  PubMed  Google Scholar 

  50. Brignac-Huber LM, Park JW, Reed JR, Backes WL. Cytochrome P450 organization and function are modulated by endoplasmic reticulum phospholipid heterogeneity. Drug Metab Dispos. 2016;44(12):1859–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Bassard JE, Mutterer J, Duval F, Werck-Reichhart D. A novel method for monitoring the localization of cytochromes P450 and other endoplasmic reticulum membrane associated proteins: a tool for investigating the formation of metabolons. FEBS J. 2012;279(9):1576–83.

    Article  CAS  PubMed  Google Scholar 

  52. Pedelacq JD, Cabantous S, Tran T, Terwilliger TC, Waldo GS. Engineering and characterization of a superfolder green fluorescent protein. Nat Biotechnol. 2006;24(1):79–88.

    Article  CAS  PubMed  Google Scholar 

  53. Leber R, Zenz R, Schrottner K, Fuchsbichler S, Puhringer B, Turnowsky F. A novel sequence element is involved in the transcriptional regulation of expression of the ERG1 (squalene epoxidase) gene in Saccharomyces cerevisiae. Eur J Biochem. 2001;268(4):914–24.

    Article  CAS  PubMed  Google Scholar 

  54. Dai Z, Liu Y, Zhang X, Shi M, Wang B, Wang D, et al. Metabolic engineering of Saccharomyces cerevisiae for production of ginsenosides. Metab Eng. 2013;20:146–56.

    Article  CAS  PubMed  Google Scholar 

  55. Lee M-H, Jeong J-H, Seo J-W, Shin C-G, Kim Y-S, In J-G, et al. Enhanced triterpene and phytosterol biosynthesis in Panax ginseng overexpressing squalene synthase gene. Plant Cell Physiol. 2004;45(8):976–84.

    Article  CAS  PubMed  Google Scholar 

  56. Seo J-W, Jeong J-H, Shin C-G, Lo S-C, Han S-S, Yu K-W, et al. Overexpression of squalene synthase in Eleutherococcus senticosus increases phytosterol and triterpene accumulation. Phytochemistry. 2005;66(8):869–77.

    Article  CAS  PubMed  Google Scholar 

  57. Kim YS, Cho JH, Park S, Han JY, Back K, Choi YE. Gene regulation patterns in triterpene biosynthetic pathway driven by overexpression of squalene synthase and methyl jasmonate elicitation in Bupleurum falcatum. Planta. 2011;233(2):343–55.

    Article  CAS  PubMed  Google Scholar 

  58. Madsen KM, Udatha GD, Semba S, Otero JM, Koetter P, Nielsen J, et al. Linking genotype and phenotype of Saccharomyces cerevisiae strains reveals metabolic engineering targets and leads to triterpene hyper-producers. PLoS ONE. 2011;6(3):e14763.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. Jackson BE, Hart-Wells EA, Matsuda SP. Metabolic engineering to produce sesquiterpenes in yeast. Org Lett. 2003;5(10):1629–32.

    Article  CAS  PubMed  Google Scholar 

  60. Scalcinati G, Knuf C, Partow S, Chen Y, Maury J, Schalk M, et al. Dynamic control of gene expression in Saccharomyces cerevisiae engineered for the production of plant sesquitepene alpha-santalene in a fed-batch mode. Metab Eng. 2012;14(2):91–103.

    Article  CAS  PubMed  Google Scholar 

  61. Kirby J, Romanini DW, Paradise EM, Keasling JD. Engineering triterpene production in Saccharomyces cerevisiae-β-amyrin synthase from Artemisia annua. FEBS J. 2008;275(8):1852–9.

    Article  CAS  PubMed  Google Scholar 

  62. Cao X, Lv YB, Chen J, Imanaka T, Wei LJ, Hua Q. Metabolic engineering of oleaginous yeast Yarrowia lipolytica for limonene overproduction. Biotechnol Biofuels. 2016;9(1):214.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  63. Jingwen Z, Xiaoxia Y, Catherine M, Guocheng D, Jian C. Enhanced α-ketoglutarate production in Yarrowia lipolytica WSH-Z06 by alteration of the acetyl-CoA metabolism. J Biotechnol. 2012;161(3):257–64.

    Article  CAS  Google Scholar 

  64. Chen Y, Daviet L, Schalk M, Siewers V, Nielsen J. Establishing a platform cell factory through engineering of yeast acetyl-CoA metabolism. Metab Eng. 2013;15:48–54.

    Article  CAS  PubMed  Google Scholar 

  65. Corbacho I, Teixido F, Velazquez R, Hernandez LM, Olivero I. Standard YPD, even supplemented with extra nutrients, does not always compensate growth defects of Saccharomyces cerevisiae auxotrophic strains. Antonie Van Leeuwenhoek. 2011;99(3):591–600.

    Article  CAS  PubMed  Google Scholar 

  66. Guo XJ, Xiao WH, Wang Y, Yao MD, Zeng BX, Liu H, et al. Metabolic engineering of Saccharomyces cerevisiae for 7-dehydrocholesterol overproduction. Biotechnol Biofuels. 2018;11:192.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  67. Wang M, Liu GN, Liu H, Zhang L, Li BZ, Li X, et al. Engineering global transcription to tune lipophilic properties in Yarrowia lipolytica. Biotechnol Biofuels. 2018;11:115.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

Download references

Authors’ contributions

CCJ and JLZ performed the experiments. YXC, CCJ, JLZ and HS conceived of this study and drafted the manuscript. All authors read and approved the final manuscript.

Acknowledgements

The authors are grateful for the financial support from the Major Program of The National Natural Science Foundation of China (21621004) and the Young Elite Scientists Sponsorship Program by Tianjin (TJSQNTJ-2018-16).

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Data will be made available from the corresponding author on reasonable request.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Funding

Major Program of The National Natural Science Foundation of China (21621004) and Young Elite Scientists Sponsorship Program by Tianjin (TJSQNTJ-2018-16).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Author notes

  1. Cong-Cong Jin and Jin-Lai Zhang equally contributed to this work

Authors and Affiliations

  1. Frontier Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin, 300072, China

    Cong-Cong Jin, Jin-Lai Zhang, Hao Song & Ying-Xiu Cao

  2. Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), School of Chemical Engineering and Technology, Tianjin University, Tianjin, 300072, China

    Cong-Cong Jin, Jin-Lai Zhang, Hao Song & Ying-Xiu Cao

Authors
  1. Cong-Cong Jin
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Jin-Lai Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Hao Song
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Ying-Xiu Cao
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Ying-Xiu Cao.

Additional files

Additional file 1: Figure S1.

Structures of the endogenous compound 2,3-oxidosqualene and the heterogenous compounds lupeol, betulin, betulinic aldehyde and betulinic acid.

Additional file 2: Figure S2.

Production of lupeol in the parent strain 201249 and YLJCC0 (expressing RcLUS from R. communis).

Additional file 3: Table S1.

Triterpenoid production by strains with or without introduction of a plant CPR (LiCPR or MTR). No matter whether heterologous CPRs were expressed, no betulinic acid was observed, which indicated that the native P450 enzymes in Y. lipolytica did not catalyze betulinic acid production.

Additional file 4: Figure S3.

GC–MS profiles and mass spectra of chemical standards and fermentation extracts from YLJCC6. (a) GC–MS profiles of the (1) lupeol standard, (2) betulin standard, (3) betulinic acid standard, and (4) metabolite extracts of YLJCC6. (b) Mass spectra of lupeol and betulinic acid. Red represents the mass spectrum of the metabolite extracts of YLJCC6, and gray represents the mass spectrum of the standards.

Additional file 5: Figure S4.

(a) Phylogenetic tree of plant CYPs that were identified or deduced as lupeol c-28 oxidases constructed by MEGA6. (b) Homology analysis of CrAO, MtAO12, VvAO15 and VvAO17 with BPLO by BLASTP.

Additional file 6: Figure S5.

Confocal image of strains that expressed CYP or CPR fused with sfGFP. sfGFP was fused to the C-terminus of CYP or CPR.

Additional file 7: Table S2.

Oligonucleotides used in construct plasmids in this study. Table S3. The Codon-optimized gene sequences involved in this study.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jin, CC., Zhang, JL., Song, H. et al. Boosting the biosynthesis of betulinic acid and related triterpenoids in Yarrowia lipolytica via multimodular metabolic engineering. Microb Cell Fact 18, 77 (2019). https://doi.org/10.1186/s12934-019-1127-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12934-019-1127-8

Keywords

Microbial Cell Factories

ISSN: 1475-2859

Contact us