Production of Triterpene Ginsenoside Compound K in the Non

Feb 13, 2019 - Compound K (CK) is a rare, tetracyclic, triterpenoid compound with important medical properties, such as antitumor and anti-inflammator...
2 downloads 0 Views 758KB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Biotechnology and Biological Transformations

Production of triterpene ginsenoside compound K in the nonconventional yeast Yarrowia lipolytica Dashuai Li, Yufen Wu, Chuanbo Zhang, Jie Sun, Zhijiang Zhou, and Wenyu Lu J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.9b00009 • Publication Date (Web): 13 Feb 2019 Downloaded from http://pubs.acs.org on February 14, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 37

Journal of Agricultural and Food Chemistry

Production of triterpene ginsenoside compound K in the nonconventional yeast Yarrowia lipolytica Dashuai Lia#, Yufen Wua#, Chuanbo Zhanga#, Jie Suna, Zhijiang Zhoua, Wenyu Lua, b, c*

(a) School of Chemical Engineering and Technology, Tianjin University, Tianjin, People’s Republic of China (b) Key Laboratory of System Bioengineering (Tianjin University), Ministry of Education, Tianjin, People’s Republic of China (c) SynBio Research Platform, Collaborative Innovation Center of Chemical Science and Engineering (Tianjin), Tianjin, People’s Republic of China #these

authors contributed equally to this work.

*Corresponding

author: Wenyu Lu

Tel: +86-022-27892132; Fax: +86-022-27892132. E-mails: [email protected]

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

1

Abstract

2

Compound K (CK) is a rare, tetracyclic, triterpenoid compound with important

3

medical properties, such as antitumor and anti-inflammatory activity. Herein, an

4

efficient biosynthetic pathway of CK was constructed in metabolically engineered

5

Yarrowia lipolytica for the first time, and the engineered strain, YL-CK0, produced 5.1

6

mg/L of CK. The production of CK was further increased by 5.96-fold to 30.4 mg/L

7

with overexpression of key genes in the MVA pathway and fusion of cytochrome P450

8

monooxygenase (PPDS) and NADPH-P450 reductase. Finally, 161.8 mg/L CK

9

production was achieved by fed-batch fermentation in a 5-L fermenter using the strain

10

YL-MVA-CK. To the best of our knowledge, this is the first report on heterologous CK

11

synthesis with the highest titer in Y. lipolytica. This study demonstrates the feasibility

12

of producing high-value triterpenoid compounds using Y. lipolytica as a platform.

13 14

Keywords:

15

Yarrowia lipolytica; compound K; cytochrome P450 monooxygenase; metabolic

16

engineering

17

2

ACS Paragon Plus Environment

Page 2 of 37

Page 3 of 37

Journal of Agricultural and Food Chemistry

18

1. Introduction

19

Triterpenoids, a member of large isoprene families, are C30 natural compounds

20

with important medical values, such as anticancer, antitumor, cytotoxic, antiviral, and

21

antibiotic properties1,

22

belongs to dammarane-type ginsenosides. CK, as one of the active ingredients in the

23

valuable Chinese medicine, ginseng, has been reported to exhibit strong antitumor

24

activity and to have an effective role in treating heart disease and depression, improving

25

immunity, and enhancing learning ability4-7. Phytoextraction, enzymatic conversion,

26

and chemical synthesis have been applied to produce CK8, 9, 10. However, these methods

27

were limited owing to the shortage of raw materials and the complexity of traditional

28

preparation processes. On one hand, CK is rarely accumulated in ginseng plants whose

29

growth is affected by climate and seasons11, 12. On the other hand, CK production via

30

enzymatic conversion and chemical synthesis led to high production costs and

31

complicated by-products13, so it is necessary to explore other methods for CK

32

production.

2, 3.

Compound K (CK), a tetracyclic triterpenoid compound,

33

Synthetic biology has made significant progress in producing natural products via

34

different microorganism hosts, such as ginsenoside Rh2 production in Saccharomyces

35

cerevisiae14, β-carotene production in Escherichia coli15, and dammarenediol-II

36

production in Pichia pastoris or E. coli16, 17. In a previous study, the biosynthesis of CK

37

has been successfully realized in S. cerevisiae by introducing an exogenous CK

38

synthesis pathway, including the dammarenediol-II synthase (DS), cytochrome P450

39

enzyme (PPDS), and UDP-glycosyltransferase (UGT71A28) from P. ginseng, together 3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

40

with an NADPH-cytochrome P450 reductase ATR2-1 from Arabidopsis thaliana. The

41

engineered S. cerevisiae produced only 1.4 mg/L of CK18. The CK pathway was also

42

introduced into the tobacco genome by Gwak19, and the transgenic plant yielded 1.55–

43

2.64 μg/g DCW of CK. However, the production volume cannot meet the consumer

44

demand.

45

An unconventional yeast, Y. lipolytica, generally regarded as safe (GRAS), has

46

been modified to produce valuable compounds such as fatty acids, oils, biofuels, and

47

few terpenoids. Due to its convenient genetic manipulation, the robust acetyl-CoA

48

synthesis, and energy supply system, Y. lipolytica is considered as a potential host for

49

heterologous terpenoids synthesis20-23. Like S. cerevisiae, Y. lipolytica contains a native

50

mevalonate (MVA) pathway that provides the prerequisite compounds, IPP and

51

DMAPP24, which further demonstrates that Y. lipolytica could serve as a better natural

52

host for terpenoids synthesis, as S. cerevisiae lacks acetyl-CoA25, 26. In recent years,

53

engineering Y. lipolytica has been applied to synthesis of terpenoids, such as

54

monoterpene (limonene), sesquiterpenoid (α-farnesene), and tetraterpene (β-carotene,

55

lycopene). Various metabolic strategies have been applied to enhance the yield of these

56

terpenoids, such as increasing the cytosolic acetyl-CoA supply, optimizing the

57

mevalonate (MVA) pathway, and balancing redox cofactors. The overexpression of

58

HMG1 and ERG12 genes, involved in the MVA pathway, can effectively increase the

59

production of limonene27. The sesquiterpene α-farnesene yield was increased by 20.8-

60

fold by tHMG1, IDI, and ERG20 overexpression28. Besides, HMG1 overexpression has

61

also been reported to increase the production of tetraterpene, lycopene, in Y. lipolytica29. 4

ACS Paragon Plus Environment

Page 4 of 37

Page 5 of 37

Journal of Agricultural and Food Chemistry

62

However, metabolic engineering Y. lipolytica for triterpenoid synthesis, especially for

63

triterpenoid saponins, has been rarely reported. Moreover, Y. lipolytica can use a variety

64

of cheap hydrophobic or hydrophilic carbon sources30, such as fatty acids and waste oil,

65

which makes it more suitable for industrial applications31.

66

In this study, Y. lipolytica was selected as the chassis to de novo synthesize

67

triterpenoid saponin CK. Multi-step metabolic engineering strategies were applied to

68

increase the production of CK to 30.4 mg/L. The CK production was further increased

69

to 161.8 ± 9.3 mg/L via fed-batch fermentation in the 5-L fermenter. To our knowledge,

70

this is the first report on the production of triterpenoid saponin in Y. lipolytica and the

71

highest titer of CK in microbial cell factories, which implys the feasibility for

72

production of other triterpenoids in Y. lipolytica.

73

2. Materials and Methods

74

2.1 Strains, Plasmids, and Medium

75

Y. lipolytica ATCC 201249 was used as the parental strain for strain engineering

76

and plasmid pINA1269 was used as the vector for gene expression32, 33. Y. lipolytica

77

ATCC 201249 and plasmid pINA1269 were kindly provided by Professor Yingjin Yuan

78

(School of Chemical Engineering and Technology, Tianjin University). E. coli DH5α

79

was grown in LB medium (1% tryptone, 0.5% yeast extract, and 1% NaCl) with

80

ampicillin (100 mg/L) at 37 °C for plasmid construction. All Y. lipolytica strains were

81

cultivated in YPD medium (2% glucose, 2% peptone, and 1% yeast extract) at 30 °C

82

and 200 rpm for 4 days. SC medium (0.67% yeast nitrogen base, 2% glucose, and 2%

5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

83

agar) lacking leucine or uracil was used for screening Y. lipolytica transformants. Feed

84

solution (containing: glucose, 400 g/L; KH2PO4, 9 g/L; MgSO4· 7H2O, 5.12 g/L; K2SO4,

85

3.5 g/L; Na2SO4, 0.28 g/L; adenine, 0.5 g/L; uracil, 0.6 g/L; lysine, 1.2 g/L) at pH 5.5

86

was used for fed-batch fermentation.

87

2.2 Plasmid construction

88

Primers used for plasmids and expression cassettes constructing are listed in Table

89

S1. Promoters and terminators were amplified from the genomic DNA of Y. lipolytica

90

ATCC 201249. The DS (GenBank: AB265170.1), PPDS (GenBank: JN604537.1),

91

UGT1 (GenBank: AIE12479.1), and ATR1 (GenBank: BT008426. 1) genes were

92

synthesized and cloned into pUC57 plasmids by GENEWIZ (Suzhou, China) with

93

codon optimization for Y. lipolytica (Codon optimized sequences are showed in Table

94

S2.). Other genes used for this study were all amplified from the genomic DNA of Y.

95

lipolytica ATCC 201249. ClonExpress® II One Step Cloning Kit (Vazyme, Nanjing)

96

was used for inserting PCR fragments into proper restriction sites of the parental

97

plasmid. Figure S1 illustrates the construction processes of pINA1269-HE9-20 from

98

pINA1269. For single gene insertion, primer pairs PmlI-gene-F and PmlI-gene-R were

99

used for gene amplification by PCR, and the PCR fragment was then cloned into the

100

PmlI site of pINA1269 resulting in pINA1269-gene. Subsequently, ERG1, ERG9,

101

ERG12, ERG19, and ERG20 expression cassettes were amplified from pINA1269-

102

ERG1, pINA1269-ERG9, pINA1269-ERG12, pINA1269-ERG19, and pINA1269-

103

ERG20 with primers P-X-T-SalI-F and P-X-T-SalI-R, and inserted into the SalI site of

104

pINA1269-tHMG1 resulting in pINA1269-HE1, pINA1269-HE9, pINA1269-HE12, 6

ACS Paragon Plus Environment

Page 6 of 37

Page 7 of 37

Journal of Agricultural and Food Chemistry

105

pINA1269-HE19, and pINA1269-HE20, respectively. The ERG20 gene expression

106

cassette was amplified from pINA1269-ERG20 using the primers P-X-T-ClaI-F and P-

107

X-T-ClaI-R; the fragments were inserted into the ClaI site of pINA1269-HE9 to get

108

pINA1269-HE9-20. The pINA1269-HE1-9 was constructed by inserting ERG9

109

expression cassette (hp4d-ERG9-XPR2t) into the ClaI site of pINA1269-HE1, and

110

pINA1269-HE9-UGT1 was constructed by inserting UGT1 expression cassette (hp4d-

111

UGT1-XPR2t) into the ClaI site of pINA1269-HE9. Constructed plasmids are shown in

112

Table 1.

113

2.3 Strain construction

114

The pINA1269 derived plasmids were digested with Not I, and the gene expression

115

cassettes used for genome integration were constructed by fusion PCR as shown in

116

Figure S2. The linearized plasmids and gene expression cassettes were transformed to

117

Y. lipolytica by the LiAc/ssDNA/PEG method according to our previously reported

118

methods34, 35, 36. Constructed strains are shown in Table 2.

119

2.4 Copy number assay

120

Gene copy number assays were performed by qPCR with absolute quantification

121

method37. ERG5 gene was selected as a reference gene. Open reading frames of ERG5,

122

OpDS, PPDS, and ATR1 were amplified and quantified to generate standard curves.

123

Genome DNA of different strains were extracted using a TIANamp Yeast DNA kit

124

(Tiangen, Beijing, China). TaqMan probe-based determinations were performed as

125

described by Zhang38. All experiments were performed in triplicates. Primers used for

7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

126

qPCR are listed in Table S3.

127

2.5 Flask and 5-L bioreactor fermentation

128

Y. lipolytica strains were inoculated into 5 mL of YPD medium as the pre-culture

129

for 24 h at 30 °C and 200 rpm and the culture was transferred to 250-mL shake flasks

130

loaded with 50 mL YPD medium for 5 days under the same cultivation conditions. All

131

shake flask fermentation experiments were carried out in 3 parallel experimental groups.

132

For batch fermentation in the 5-L bioreactor (Bailun, Shanghai, China), 5 mL of pre-

133

culture was transferred to 150 mL of YPD medium in 500-mL shake flasks, and

134

cultivated for 24 h at 30 °C and 220 rpm. The seed medium was inoculated into 5-L

135

bioreactor with an initial OD600 of 0.5 in 2 L of YPD medium. The temperature was

136

maintained at 30 °C. The pH was maintained at 6.8 by automatically adding 5 M

137

ammonia water or 3 M H2SO4. The rotating speed was set at 500 rpm with an air flow

138

rate of 2 vvm. When the glucose concentration was lower than 2 g/L, the feed solution

139

was added to maintain the glucose concentration at below 5 g/L.

140

2.6 Metabolite extraction and analysis

141

Cell growth was determined by measuring the optical density at 600 nm (OD600)

142

using a UV-VIS spectrophotometer. Glucose concentration was measured via a

143

bioanalyzer (SBA-40C, Shandong Academy of Sciences, China). Squalene, 2, 3-

144

oxidosqualene, dammarenediol-II, and PPD were extracted with ethyl acetate,

145

determined by LC/APCI/MS, and quantified by HPLC following the protocols

146

previously reported39.

8

ACS Paragon Plus Environment

Page 8 of 37

Page 9 of 37

Journal of Agricultural and Food Chemistry

147

CK was extracted with n-butanol. Five hundred microliters of n-butanol and 0.5 g

148

SiO2 (0.5 mm diameter) were added into 1 mL of culture broth, and the mixture was

149

agitated by vortexing for 20 min. Then, the mixture was centrifuged at 11564 (× g) rpm

150

for 10 min, and the n-butanol phase was collected for HPLC analysis. HPLC analysis

151

was carried out on an Elite (Dalian, China) HPLC system equipped with an Elite P230II

152

pump and an UV230II detector at 203 nm. The mobile phase was 80% acetonitrile at a

153

flow rate of 1 mL/min for CK. Standards which were purchased from Meilun

154

Biotechnology Co., Ltd (Dalian, China). The LC-MS analysis of DMD, PPD, and CK

155

is shown in Figure S3.

156

3. Results and Discussion

157

3.1 Production of DMD and PPD in Y. lipolytica ATCC 201249

158

PPD, precursor of CK, cannot be directly generated in Y. lipolytica. Although Y.

159

lipolytica has a native MVA pathway, but it cannot directly synthesize CK due to the

160

lack of dammarenediol-II synthase (DS), cytochrome P450 enzyme (PPDS), and UDP-

161

glycosyltransferase (UGT1) (Figure 1).

162

The plasmid pINA1269 (pINA1269-DS) containing a DS gene derived from

163

Panax ginseng was transformed into Y. lipolytica ATCC 201249 resulting in strain YL-

164

DS. However, DMD was not detected in strain YL-DS. It has been reported that Y.

165

lipolytica has a special codon preference for some heterologous genes27, 28. Here, plant-

166

derived DS gene might not be correctly expressed in Y. lipolytica. To solve this problem,

167

a codon-optimized DS gene was integrated into single (pINA1269-OpDS) and muti-

9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

168

copy sites (rDNA)33 of the Y. lipolytica genome resulting in YL-OpDS0 and YL-OpDS1,

169

respectively. As shown in Figure 2, DMD was detected in both strains and the yield of

170

DMD in YL-OpDS1 (17.9 ± 2.3 mg/L) was 2.2-fold higher than that in YL-OpDS0 (7.9

171

± 1.3 mg/L). The levels of DMD precursors, squalene and 2, 3-oxidosqualene, also

172

decreased in YL-OpDS1. It has been reported that the expression level of the first

173

exogenous gene that competed with precursors directly affected the yield of the target

174

product during heterologous steroid synthesis40. Then, we inserted TEF1p-OpDS-

175

TXPR2t, EXP1p-PPDS-MIG1t, and GPD1p-ATR1-LIP2t expression modules into the

176

Y. lipolytica ATCC 201249 multi-copy rDNA sites obtaining strain YL-PPDS0. The

177

YL-PPDS0 strain produced 12.9 mg/L of PPD which was mostly detected in

178

extracellular (10.4 mg/L) medium. However, DMD was only detected in intracellular

179

(4.8 mg/L) compartment. These results demonstrated Y. lipolytica has the ability to

180

express heterologous P450 systems, which indicated the feasibility of synthesizing

181

tetracyclic triterpenoid products in Y. lipolytica.

182

However, the production of PPD was still low, presumably due to the following

183

reasons. On one hand, the precursor flux is low in the endogenous MVA pathway, such

184

as acetyl-CoA, IPP, and FPP. On the other hand, low PPD production may be due to

185

the compatibility issues between the heterologous genes and the endogenous

186

metabolisms, such as the endogenous P450 enzymes competed with the introduced

187

PPDS for the electrons transmitted by ATR1. In fact, DMD was not fully converted

188

with a conversion rate of 72%. Therefore, we speculated that increasing the flux of

189

precursors in the triterpene synthesis pathway and optimizing the heterologous P450 10

ACS Paragon Plus Environment

Page 10 of 37

Page 11 of 37

Journal of Agricultural and Food Chemistry

190

oxidation system could improve PPD production.

191

3.2 Overexpression of key genes of MVA pathway to improve the production of

192

PPD in Y. lipolytica

193

It has been reported that β-carotene production increased by 134% when truncated

194

hydroxymethylglutaryl-CoA reductase (tHMGR) from Y. lipolytica was overexpressed.

195

Moreover, additional two tHMGR1 copy numbers effectively increased β-carotene

196

production in Y. lipolytica41. Overexpression of hydroxymethylglutaryl-CoA reductase

197

(HMGR1) was beneficial to lycopene production in Y. lipolytica29, and to beta-carotene

198

production in S. cerevisiae42. Overexpression of enzymes, tHMG1 and ERG9, in the

199

MVA pathway increased betulinic acid production in Y. lipolytica36. Therefore, to

200

evaluate the influences on the production of PPD in Y. lipolytica, the MVA pathway

201

(Figure 1) genes were analyzed in this study. Firstly, MVA pathway genes including

202

ERG1, ERG8, ERG9, ERG10, ERG12, ERG13, ERG19, ERG20, IDI, HMG1, and

203

tHMG1 (The sequence encoding the first 511 amino acids at the N-terminal of HMG1

204

was deleted) were expressed in YL-PPDS0 using the single copy vector pINA1269

205

(Figure 3A). Except for YL-PPDS9 that overexpressed IDI, the synthesis of rest of the

206

PPDs were improved in strains that overexpressed other genes. Particularly, PPD

207

production in YL-PPDS11 (tHMG1 overexpression) increased by 230% compared to

208

that in YL-PPDS0. Overexpression of ERG1, ERG9, ERG12, ERG19, ERG20, and

209

HMG1 increased PPD production by 84%, 110%, 96%, 38%, 87%, and 180%,

210

respectively. The results indicated that HMG1 is also an important rate-limiting enzyme

211

for triterpenoid synthesis in Y. lipolytica, followed by ERG9 and ERG20. In addition, 11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

212

these results also corroborated with previous studies that overexpressed key enzymes

213

in the MVA pathway could significantly enhance the yield of terpenoids in the Y.

214

lipolytica or S. cerevisiae43, 44, 45. In this study, both overexpressed HMG1 and tHMG1

215

increased the titer of PPD, but the effect of tHMG1 was more prominent. The

216

hydroxymethylglutaryl-CoA reductase HMG1 contains an N-terminal transmembrane

217

region and a C-terminal catalytically active region, which respectively play a role in

218

transmembrane localization and catalysis. HMG1 has been recognized as the most

219

important rate-limiting enzyme in the MVA pathway, and is an important regulatory

220

point in the metabolism of terpenoids in the cytoplasm. The increase in HMG1

221

expression level was beneficial to increase the supply of mevalonate in the MVA

222

pathway. Overexpression of HMG1 resulted in feedback inhibition of the mevalonate

223

metabolic pathway, during which downstream products activate HMG1 on the

224

endoplasmic reticulum for degradation. The tHMG1 (deleted N-terminal) avoids the

225

self-degradation mediated by its N-terminal domain, which enhances protein stability

226

and increases the metabolic flow of mevalonate46, 47. Therefore, the tHMG1 was used

227

for later strain construction instead of HMG1.

228

In order to further increase the titer of PPD, ERG1, ERG9, ERG12, ERG19, and

229

ERG20 were selected to co-express with tHMG1 in YL-PPDS0 respectively resulting

230

in strains YL-PPDS12-16 (Figure 3B). The highest PPD production of 65.5 mg/L was

231

achieved by overexpressing tHMG1 and ERG9 (YL-PPDS13). However, the PPD

232

production in YL-PPDS14 and YL-PPDS15 strains did not significantly increase

233

compared to its parental strain, YL-PPDS11. Therefore, ERG20 or ERG1 was co-

12

ACS Paragon Plus Environment

Page 12 of 37

Page 13 of 37

Journal of Agricultural and Food Chemistry

234

expressed with tHMG1 and ERG9 to obtain strain YL-PPDS17 and YL-PPDS18,

235

respectively. The production of PPD in YL-PPDS17 increased by 5-fold compared to

236

the original strain (YL-PPDS0), reaching 82.5 mg/L. The above results suggested that

237

the synergy between tHMG1 and ERG9 was the most effective for synthesizing high

238

levels of PPD, followed by ERG20. Therefore, the strategy of metabolic engineering to

239

overexpress key genes in the MVA pathway exhibited a significant effect on increasing

240

the metabolic flux of the precursors and promoting the accumulation of products.

241

As shown in Table 3, when the upstream MVA pathway was overexpressed,

242

although the PPD production was enhanced, DMD and endogenous sterol accumulation

243

was also increased. It is possible to explain why the ergosterol content increased as the

244

PPD synthesis yield increased. Therefore, it was speculated that DMD could not be

245

completely transformed, which was the key to strengthen the potentiality of

246

heterologous triterpenes. Finding a more suitable CPR (ATR1) for P450 (PPDS) was

247

an effective way to improve the P450 catalytic efficiency. In previous studies, we found

248

that PPDS catalytic efficiency plays an important role in balancing the heterologous

249

pathway in S. cerevisiae, and PPDS-ATR1 fusion could increase the catalytic efficiency

250

in S. cerevisiae39. Therefore, the fusion protein strategy was implemented in Y.

251

lipolytica. Strains YL-fuPPDS1 and YL-fuPPDS2 were obtained by PPDS and

252

t46AtCPR1 fusion with linker1 (GSTSSGSSG) and linker2 (GSTSSG)34,

253

fuPPDS2 accumulated 109.6 mg/L of PPD and 2.3 mg/L of DMD with a conversion

254

rate of 98%, which was in accordance with the effects in S. cerevisiae (Figure 3C). The

255

artificially constructed PPDS-ATR1 fusion protein improved the efficiency of electron

13

ACS Paragon Plus Environment

39.

YL-

Journal of Agricultural and Food Chemistry

256

transfer and the activity of P450 oxidation system in Y. lipolytica. This indicated that

257

fusion of key enzymes in the pathway may reduce loss of intermediates by diffusion,

258

degradation, or conversion of the intermediate through competitive pathways48.

259

3.3 Construction of CK synthesis pathway in Y. lipolytica

260

Protopanaxadiol (PPD) is a direct precursor compound of CK. In order to produce

261

CK in Y. lipolytica, the plasmid pINA1269-UGT1 carrying the codon-optimized

262

glycosyltransferase gene was integrated into the genome of YL-PPDS0 to obtain the

263

strain YL-CK0. Then, UGT1 was co-expressed with tHMG1, ERG20, and ERG9

264

resulting in YL-CK1. As shown in Figure 4, the production of CK in YL-CK1 was 21.7

265

mg/L which was 4.2-fold higher than that in YL-CK0 (5.1 mg/L). The results showed

266

that enhancing the metabolic flux of the precursor was important for increase in CK

267

production via the overexpression of the MVA pathway. Expressing the key genes of

268

MVA pathway, including tHMG1, ERG20, and ERG9, and genes involved in CK

269

synthesis pathway, including OpDS, PPDS, ATR1, and UGT1, into the Y. lipolytica

270

genome muti-copy sites obtained the strain YL-MVA-CK (Figure S2). YL-MVA-CK

271

accumulated 30.4 mg/L of CK which was 5.96-fold higher than that of YL-CK0. It

272

exhibited a significant increase compared with previous reported 1.4 mg/L yield in S.

273

cerevisiae18. To our knowledge, this is the first report of CK biosynthesis in oleaginous

274

yeast Y. lipolytica. Multi-copy integration strategy in key genes of MVA pathway and

275

fusion strategy between PPDS and ATR1 effectively increased the yield of CK.

276

However, PPD was highly accumulated in the final CK producing strain YL-MVA-CK.

277

Thus, enhancement in the UDP-glucose supply and glycosyltransferase activity could 14

ACS Paragon Plus Environment

Page 14 of 37

Page 15 of 37

Journal of Agricultural and Food Chemistry

278

be the feasible methods to solve this problem.

279

3.4 Copy number assay

280

In order to verify the effect of copy number on the results, we used qPCR to detect

281

the copy numbers of DS, PPDS, and ATR1 in screened strains (YL-OpDS1, YL-PPDS0,

282

YL-fuPPDS1, YL-fuPPDS2, and YL-MVA-CK). Based on the results of qPCR (Figure

283

S4), we found that the copy numbers of the coding genes DS, PPDS, and ATR1 at muti-

284

copy sites were more than one in different engineered strains, which indicated that the

285

multi-copy integration strategy has a positive impact on improving substrate conversion

286

rate. In Saccharomyces cerevisiae, the conversion of DMD and the reactive oxygen

287

species (ROS) level increased with increasing copy numbers of PPDS and ATR1, but

288

three copies of PPDS-ATR1 fusion proteins could convert 96.8% DMD to PPD39. In

289

this study, two copies of PPDS-ATR1 fusion proteins were sufficient to convert 98%

290

DMD to PPD in Y. lipolytica. Compared with the starting strain Y. lipolytica 201249,

291

there was no significant difference in the growth of strains YL-fuPPDS1 and YL-

292

fuPPDS2. Therefore, in this study, we believed that the ROS level caused by two copy

293

numbers has no significant effect on cell growth and the synthesis of product. Further

294

research is needed to explore the effect of more copy number or ROS level on the

295

synthesis of PPD and CK.

296

3.5 Fed-batch fermentation for CK production in the 5-L bioreactor

297

According to the growth curve and sugar consumption characteristics of the

298

engineered strain, YL-MVA-CK, in flaks, fed-batch fermentation was performed in a

299

5-L bioreactor with a working volume of 2.5 L. The results of biomass and CK 15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

300

production during fermentation course are shown in Figure 5. In this study, after 72

301

hours of fermentation, we found a CK titer of 161.8 mg/L from 100 g (total glucose)/L

302

(yield = 1.61 mg/g sugars, productivity = 2.25 mg/L/h) feed. This was the highest yield

303

reported. In the next 48 hours, the production of CK had mostly increased. It can be

304

considered that the short fermentation cycle for the production of CK in Y. lipolytica

305

has important industrial value.

306

Conflicts of interest

307 308

The authors declare no conflicts of interest. Acknowledgments

309

This work was financially supported by the National Natural Science Foundation

310

of China (No. 21878220) and the Major Research Plan of Tianjin [No.

311

16YFXTSF00460].

312

Supporting information

313

Table S1. All the primers used in this study.

314

Table S2. Codon optimized sequences of DS, PPDS, ATR1 and UGT1 for Y.lipolytica.

315

Table S3. Primers for qPCR.

316

Figure S1. The diagram of plasmid p1269-HE9-20 construction process.

317

Figure S2. The gene expression cassettes used in engineered strain YL-MVA-CK.

318

Figure S3. LC/MS analysis of DMD, PPD and CK.

319

Figure S4. Copy number of OpDS, PPDS and ATR1 in different engineered strains.

16

ACS Paragon Plus Environment

Page 16 of 37

Page 17 of 37

Journal of Agricultural and Food Chemistry

320

References

321

1.

Kim, Y. S.; Kim, D. S.; Kim, S. I., Ginsenoside Rh2 and Rh3 induce differentiation

322

of HL-60 cells into granulocytes: modulation of protein kinase C isoforms during

323

differentiation by ginsenoside Rh2. The international journal of biochemistry &

324

cell biology 1998, 30, 327-338.

325

2.

Cheng, Y.; SHEN, L. h.; ZHANG, J. t., Anti‐amnestic and anti‐aging effects of

326

ginsenoside Rg1 and Rb1 and its mechanism of action. Acta Pharmacologica

327

Sinica 2005, 26, 143-149.

328

3.

Kim, Y. H.; Chung, C. B.; Kim, J. G.; Ko, K. I.; Park, S. H.; Kim, J.-H.; Eom, S.

329

Y.; Kim, Y. S.; Hwang, Y.-I.; Kim, K. H., Anti-wrinkle activity of ziyuglycoside

330

I isolated from a Sanguisorba officinalis root extract and its application as a

331

cosmeceutical ingredient. Bioscience, biotechnology, and biochemistry 2008, 72,

332

303-311.

333

4.

biosynthesis in plants. Annual Review of Plant Biology 2014, 65, 225-257.

334 335

Thimmappa, R.; Geisler, K.; Louveau, T.; O'Maille, P.; Osbourn, A., Triterpene

5.

KARIKURA, M.; TANIzAwA, H.; HIRATA, T.; MIYASE, T.; TAKINO, Y.,

336

Studies on absorption, distribution, excretion and metabolism of ginseng saponins.

337

VIII. Isotope labeling of ginsenoside Rb2. Chemical and pharmaceutical bulletin

338

1992, 40, 2458-2460.

339

6.

in plant triterpene cyclization. Current opinion in plant biology 2006, 9, 305-314.

340 341

Phillips, D. R.; Rasbery, J. M.; Bartel, B.; Matsuda, S. P., Biosynthetic diversity

7.

Ko, S.-R.; Choi, K.-J.; Uchida, K.; Suzuki, Y., Enzymatic preparation of 17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

342

ginsenosides Rg2, Rh1, and F1 from protopanaxatriol-type ginseng saponin

343

mixture. Planta medica 2003, 69, 285-286.

344

8.

Bae, E.-A.; Choo, M.-K.; Park, E.-K.; Park, S.-Y.; Shin, H.-Y.; Kim, D.-H.,

345

Metabolism of ginsenoside Rc by human intestinal bacteria and its related

346

antiallergic activity. Biological and Pharmaceutical Bulletin 2002, 25, 743-747.

347

9.

Zhou, W.; Yan, Q.; Li, J. Y.; Zhang, X. C.; Zhou, P., Biotransformation of Panax

348

notoginseng saponins into ginsenoside compound K production by Paecilomyces

349

bainier sp. 229. Journal of applied microbiology 2008, 104, 699-706.

350

10. Yoo, M.-H.; Yeom, S.-J.; Park, C.-S.; Lee, K.-W.; Oh, D.-K., Production of

351

aglycon protopanaxadiol via compound K by a thermostable β-glycosidase from

352

Pyrococcus furiosus. Applied microbiology and biotechnology 2011, 89, 1019-

353

1028.

354

11. Hasegawa, H., Proof of the mysterious efficacy of ginseng: basic and clinical trials:

355

metabolic activation of ginsenoside: deglycosylation by intestinal bacteria and

356

esterification with fatty acid. Journal of pharmacological sciences 2004, 95, 153-

357

157.

358 359

12. Sivakumar, G.; Yu, K.; Paek, K., Biosafe Ginseng: A Novel Source for Human Well‐Being. Engineering in Life Sciences 2005, 5, 527-533.

360

13. Dürr, C.; Hoffmeister, D.; Wohlert, S. E.; Ichinose, K.; Weber, M.; von Mulert,

361

U.; Thorson, J. S.; Bechthold, A., The glycosyltransferase UrdGT2 catalyzes both

362

C‐and O‐glycosidic sugar transfers. Angewandte Chemie International Edition

363

2004, 43, 2962-2965.

18

ACS Paragon Plus Environment

Page 18 of 37

Page 19 of 37

Journal of Agricultural and Food Chemistry

364

14. Zhuang, Y.; Yang, G.-Y.; Chen, X.; Liu, Q.; Zhang, X.; Deng, Z.; Feng, Y.,

365

Biosynthesis of plant-derived ginsenoside Rh2 in yeast via repurposing a key

366

promiscuous microbial enzyme. Metabolic engineering 2017, 42, 25-32.

367

15. Yoon, S.-H.; Lee, S.-H.; Das, A.; Ryu, H.-K.; Jang, H.-J.; Kim, J.-Y.; Oh, D.-K.;

368

Keasling, J. D.; Kim, S.-W., Combinatorial expression of bacterial whole

369

mevalonate pathway for the production of β-carotene in E. coli. Journal of

370

biotechnology 2009, 140, 218-226.

371

16. Liu, X.-B.; Liu, M.; Tao, X.-Y.; Zhang, Z.-X.; Wang, F.-Q.; Wei, D.-Z., Metabolic

372

engineering of Pichia pastoris for the production of dammarenediol-II. Journal of

373

biotechnology 2015, 216, 47-55.

374

17. Li, D.; Zhang, Q.; Zhou, Z.; Zhao, F.; Lu, W., Heterologous biosynthesis of

375

triterpenoid dammarenediol-II in engineered Escherichia coli. Biotechnology

376

letters 2016, 38, 603-609.

377

18. Yan, X.; Fan, Y.; Wei, W.; Wang, P.; Liu, Q.; Wei, Y.; Zhang, L.; Zhao, G.; Yue,

378

J.; Zhou, Z., Production of bioactive ginsenoside compound K in metabolically

379

engineered yeast. Cell research 2014, 24, 770.

380

19. Gwak, Y. S.; Han, J. Y.; Adhikari, P. B.; Ahn, C. H.; Choi, Y. E., Heterologous

381

production of a ginsenoside saponin (compound K) and its precursors in transgenic

382

tobacco impairs the vegetative and reproductive growth. Planta 2017, 245, 1105-

383

1119.

384

20. Groenewald, M.; Boekhout, T.; Neuvéglise, C.; Gaillardin, C.; Van Dijck, P. W.;

385

Wyss, M., Yarrowia lipolytica: safety assessment of an oleaginous yeast with a

19

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

386

great industrial potential. Critical reviews in microbiology 2014, 40, 187-206.

387

21. Senam, S.; Barth, G., Detection and activation of the Ylt1 retrotransposon in the

388

yeast Yarrowia lipolytica. In Non-Conventional Yeasts in Genetics, Biochemistry

389

and Biotechnology, Springer: 2003; pp 357-364.

390

22. Blank, L. M.; Lehmbeck, F.; Sauer, U., Metabolic-flux and network analysis in

391

fourteen hemiascomycetous yeasts. FEMS yeast research 2005, 5, 545-558.

392

23. Wasylenko, T. M.; Ahn, W. S.; Stephanopoulos, G., The oxidative pentose

393

phosphate pathway is the primary source of NADPH for lipid overproduction from

394

glucose in Yarrowia lipolytica. Metabolic engineering 2015, 30, 27-39.

395

24. Zhu, Z.; Zhang, S.; Liu, H.; Shen, H.; Lin, X.; Yang, F.; Zhou, Y. J.; Jin, G.; Ye,

396

M.; Zou, H., A multi-omic map of the lipid-producing yeast Rhodosporidium

397

toruloides. Nature communications 2012, 3, 1112.

398

25. Chen, Y.; Daviet, L.; Schalk, M.; Siewers, V.; Nielsen, J., Establishing a platform

399

cell factory through engineering of yeast acetyl-coa metabolism. Metabolic

400

Engineering 2013, 15, 48-54.

401

26.

Lian, J.; Si, T.; Nair, N.U.; Zhao, H., Design and construction of acetyl-CoA

402

overproducing Saccharomyces cerevisiae strains. Metabolic Engineering 2014, 24,

403

139-149.

404

27. Cao, X.; Lv, Y.-B.; Chen, J.; Imanaka, T.; Wei, L.-J.; Hua, Q., Metabolic

405

engineering of oleaginous yeast Yarrowia lipolytica for limonene overproduction.

406

Biotechnology for biofuels 2016, 9, 214.

407

28. Yang, X.; Nambou, K.; Wei, L.; Hua, Q., Heterologous production of α-farnesene

20

ACS Paragon Plus Environment

Page 20 of 37

Page 21 of 37

Journal of Agricultural and Food Chemistry

408

in metabolically engineered strains of Yarrowia lipolytica. Bioresource technology

409

2016, 216, 1040-1048.

410

29. Matthäus, F.; Ketelhot, M.; Gatter, M.; Barth, G., Production of lycopene in the

411

non-carotenoid-producing yeast Yarrowia lipolytica. Appl. Environ. Microbiol.

412

2014, 80, 1660-1669.

413

30. Liu, H.-H.; Ji, X.-J.; Huang, H., Biotechnological applications of Yarrowia

414

lipolytica: Past, present and future. Biotechnology Advances 2015, 33, 1522-1546.

415

31. Zhao, C.; Gu, D.; Nambou, K.; Wei, L.; Chen, J.; Imanaka, T.; Hua, Q.,

416

Metabolome analysis and pathway abundance profiling of Yarrowia lipolytica

417

cultivated on different carbon sources. Journal of biotechnology 2015, 206, 42-51.

418

32. Papanikolaou, S.; Muniglia, L.; Chevalot, I.; Aggelis, G.; Marc, I., Yarrowia

419

lipolytica as a potential producer of citric acid from raw glycerol. Journal of

420

applied microbiology 2002, 92, 737-744.

421

33. Nicaud, J.-M.; Madzak, C.; van den Broek, P.; Gysler, C.; Duboc, P.; Niederberger,

422

P.; Gaillardin, C., Protein expression and secretion in the yeast Yarrowia lipolytica.

423

FEMS yeast research 2002, 2, 371-379.

424

34. Guo, X.; Sun, J.; Li, D.; Lu, W., Heterologous biosynthesis of (+)-nootkatone in

425

unconventional yeast Yarrowia lipolytica. Biochemical Engineering Journal 2018.

426

35. Liu, H.-H.; Madzak, C.; Sun, M.-L.; Ren, L.-J.; Song, P., Huang, H.; Ji, X.-J.,

427

Engineering Yarrowia lipolytica for arachidonic acid production through rapid

428

assembly of metabolic pathway. Biochemical Engineering Journal 2017, 119, 52-

429

58.

21

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

430

Page 22 of 37

36. Sun, J.; Zhang, C.; Nan, W.; Li, D.; Ke, D.; Lu, W., Glycerol improves

431

heterologous

biosynthesis

of

betulinic

acid

432

lipolytica. Chemical Engineering Science. 2018.

in

engineered

Yarrowia

433

37. Abad, S.; Kitz, K.; Hörmann, A.; Schreiner, U.; Hartner, F. S.; Glieder, A.,

434

Real‐time PCR‐based determination of gene copy numbers in Pichia pastoris.

435

Biotechnology Journal: Healthcare Nutrition Technology 2010, 5, 413-420.

436

38. Zhang, X.; Xue, C.; Zhao, F.; Li, D.; Yin, J.; Zhang, C.; Caiyin, Q.; Lu, W.,

437

Suitable extracellular oxidoreduction potential inhibit rex regulation and effect

438

central carbon and energy metabolism in Saccharopolyspora spinosa. Microbial

439

cell factories 2014, 13, 98.

440

39. Zhao, F.; Bai, P.; Liu, T.; Li, D.; Zhang, X.; Lu, W.; Yuan, Y., Optimization of a

441

cytochrome P450 oxidation system for enhancing protopanaxadiol production in

442

Saccharomyces cerevisiae. Biotechnology and bioengineering 2016, 113, 1787-

443

1795.

444 445

40. Keasling, J. D., Synthetic biology and the development of tools for metabolic engineering. Metabolic engineering 2012, 14, 189-195.

446

41. Gao, S.; Tong, Y.; Zhu, L.; Ge, M.; Zhang, Y.; Chen, D.; Jiang, Y.; Yang, S.,

447

Iterative integration of multiple-copy pathway genes in Yarrowia lipolytica for

448

heterologous β-carotene production. Metabolic engineering 2017, 41, 192-201.

449

42. Verwaal, R.; Wang, J.; Meijnen, J.-P.; Visser, H.; Sandmann, G.; van den Berg, J.

450

A.; van Ooyen, A. J., High-level production of beta-carotene in Saccharomyces

451

cerevisiae by successive transformation with carotenogenic genes from

22

ACS Paragon Plus Environment

Page 23 of 37

Journal of Agricultural and Food Chemistry

452

Xanthophyllomyces dendrorhous. Applied and environmental microbiology 2007,

453

73, 4342-4350.

454

36. Sun, J.; Zhang, C.; Nan, W.; Li, D.; Ke, D.; Lu, W., Glycerol improves

455

heterologous biosynthesis of betulinic acid in engineered Yarrowia lipolytica.

456

Chemical Engineering Science 2018.

457

43. Dai, Z.; Liu, Y.; Zhang, X.; Shi, M.; Wang, B.; Wang, D.; Huang, L.; Zhang, X.,

458

Metabolic engineering of Saccharomyces cerevisiae for production of

459

ginsenosides. Metabolic engineering 2013, 20, 146-156.

460

44. Yuan, J.; Ching, C. B., Combinatorial engineering of mevalonate pathway for

461

improved amorpha‐4, 11‐diene production in budding yeast. Biotechnology and

462

bioengineering 2014, 111, 608-617.

463

45. Westfall, P. J.; Pitera, D. J.; Lenihan, J. R.; Eng, D.; Woolard, F. X.; Regentin, R.;

464

Horning, T.; Tsuruta, H.; Melis, D. J.; Owens, A., Production of amorphadiene in

465

yeast, and its conversion to dihydroartemisinic acid, precursor to the antimalarial

466

agent artemisinin. Proceedings of the National Academy of Sciences 2012, 109,

467

E111-E118.

468

46. Polakowski, T.; Stahl, U.; Lang, C., Overexpression of a cytosolic

469

hydroxymethylglutaryl-CoA reductase leads to squalene accumulation in yeast.

470

Applied microbiology and biotechnology 1998, 49, 66-71.

471

47. Donald, K.; Hampton, R. Y.; Fritz, I. B., Effects of overproduction of the catalytic

472

domain of 3-hydroxy-3-methylglutaryl coenzyme A reductase on squalene

473

synthesis in Saccharomyces cerevisiae. Applied and environmental microbiology

23

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

474

1997, 63, 3341-3344.

475

48. Albertsen, L.; Chen, Y.; Bach, L. S.; Rattleff, S.; Maury, J.; Brix, S.; Nielsen, J.;

476

Mortensen, U. H., Diversion of flux toward sesquiterpene production in

477

Saccharomyces cerevisiae by fusion of host and heterologous enzymes. Applied

478

and environmental microbiology 2011, 77, 1033-1040.

479

24

ACS Paragon Plus Environment

Page 24 of 37

Page 25 of 37

Journal of Agricultural and Food Chemistry

480

Figure Captions

481

Figure 1. Biosynthesis pathway of CK in Yarrowia lipolytica.

482

ERG10,

483

synthase; HMG1/HMG2, HMG-CoA reductase, ERG12, mevalonate kinase; ERG8,

484

phosphomevalonate kinase; ERG19, diphosphomevalonate decarboxylase; IDI,

485

Isopentenyl pyrophosphate isomerase; ERG20, FPP synthase; ERG9, squalene

486

synthase; ERG1, squalene monooxygenase; DS, DMD synthase; PPDS, PPD synthase;

487

ATR1, NADPH-P450 reductase; UGT1, UDP-glucose glucosyltransferase. HMG-CoA,

488

hydroxymethylglutaryl-CoA; MVA, Mevalonate acid; DMD, Dammarenediol-II; PPD,

489

Protopanaxadiol; CK, Compound K. Green represents the overexpressed genes in the

490

final strain. Blue represents the overexpressed native genes and red represents the

491

optimized heterologous genes. Three dotted arrows represent multi-step reactions.

492

Figure 2. Metabolites analysis of different engineered strains. The production of

493

DMD, squalene, and 2, 3-oxidosqualene in strains YL-OpDS0 and YL-OpDS1. YL-

494

OpDS0 represents the strain in which DS was overexpressed with single copy in Y.

495

lipolytica ATCC 201249, and YL-OpDS1 represents the strain in which DS was

496

overexpressed with multiple copies in Y. lipolytica ATCC 201249. Each data represents

497

an average ± 1 standard deviation of three parallel fermentations.

498

Figure 3. Production of PPD by metabolic engineering. (A) Comparison of PPD

499

yields between YL-PPDS0 and the strains YL-PPDS1-11. YL-PPDS1-11 represent the

500

Y. lipolytica strains that overexpressed ERG1, ERG8, ERG9, ERG10, ERG12, ERG13,

acetyl-CoA

C-acetyltransferase;

ERG13,

25

ACS Paragon Plus Environment

hydroxymethylglutaryl-CoA

Journal of Agricultural and Food Chemistry

501

ERG19, ERG20, IDI, HMG1 and tHMG1, respectively. (B) Production of PPD by

502

multiple genes overexpression. YL-PPDS12-18 represents the Y. lipolytica strains YL-

503

PPDS0 were inserted into the overexpression plasmids pINA1269-HE1, pINA1269-

504

HE9, pINA1269-HE12, pINA1269-HE19, pINA1269-HE20, pINA1269-HE9-20 and

505

pINA1269-HE1-9, respectively. (C) The PPD and DMD production of fusion P450

506

oxidation system in Y. lipolytica. Each data represents an average ± 1 standard deviation

507

of three parallel fermentations.

508

Figure 4. The biosynthesis of CK in engineered Y. lipolytica. YL-CK0 represents the

509

Y. lipolytica YL-PPDS0 was inserted into the plasmid pINA1269-UGT1. YL-CK1

510

represents the Y. lipolytica YL-CK0 was inserted in pINA1269-HE9-20. YL-MVA-CK

511

represents Y. lipolytica ATCC 201249 overexpressed tHMG1, ERG9 and ERG20 in zeta

512

sites, and overexpressed OpDS, PPDS-linker2-ATR1 and UGT1 in rDNA sites. Each

513

data represents an average ± 1 standard deviation of three parallel fermentations.

514

Figure 5. Fed-batch fermentation of the engineered strain YL-MVA-CK in 5-L

515

bioreactor. The feed solution was added at 24 h. Each data point represents an average

516

± 1 standard deviation of three individual bioreactor fermentations.

517

26

ACS Paragon Plus Environment

Page 26 of 37

Page 27 of 37

Journal of Agricultural and Food Chemistry

Table 1. Plasmids used in this study. Plasmids

Description

Source

Y. lipolytica integrative plasmid, hp4d promoter, (Nicaud et pINA1269 XPR2 terminator, leu2 selection marker pINA1269-DS

al. 2002)

pINA1269 vector containing DS from Panax This study ginseng

pINA1269-OpDS

pINA1269 vector containing codon-optimized DS

This study

pINA1269-UGT1

pINA1269 vector containing codon-optimized This study UGT1

pINA1269-ERG1

pINA1269

vector

containing

ERG1

from This study

vector

containing

ERG8

from This study

vector

containing

ERG9

from This study

vector

containing

ERG10

from This study

vector

containing

ERG12

from This study

vector

containing

ERG13

from This study

Y.lipolytica pINA1269-ERG8

pINA1269 Y.lipolytica

pINA1269-ERG9

pINA1269 Y.lipolytica

pINA1269-ERG10

pINA1269 Y.lipolytica

pINA1269-ERG12

pINA1269 Y.lipolytica

pINA1269-ERG13

pINA1269 Y.lipolytica

27

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

pINA1269-ERG19

pINA1269

Page 28 of 37

vector

containing

ERG19

from This study

vector

containing

ERG20

from This study

Y.lipolytica pINA1269-ERG20

pINA1269 Y.lipolytica

pINA1269-IDI

pINA1269 vector containing IDI from Y.lipolytica

pINA1269-HMG1

pINA1269

This study

vector

containing

HMG1

from This study

vector

containing

tHMG1

from This study

Y.lipolytica pINA1269-tHMG1

pINA1269 Y.lipolytica

pINA1269-HE1

pINA1269 vector containing tHMG1 and ERG1

This study

pINA1269-HE9

pINA1269 vector containing tHMG1 and ERG9

This study

pINA1269-HE12

pINA1269 vector containing tHMG1 and ERG12

This study

pINA1269-HE19

pINA1269 vector containing tHMG1 and ERG19

This study

pINA1269-HE20

pINA1269 vector containing tHMG1 and ERG20

This study

pINA1269-HE9-20

pINA1269 vector containing tHMG1, ERG9 and This study ERG20

pINA1269-HE1-9

pINA1269 vector containing tHMG1, ERG1 and This study ERG9

pINA1269-HE9-

pINA1269 vector containing tHMG1, ERG9 and This study

UGT1

UGT1

28

ACS Paragon Plus Environment

Page 29 of 37

Journal of Agricultural and Food Chemistry

Table 2.

Strains used in this study.

Strains Y.

Description

Source

lipolytica MATA,ura3-302,leu2-270,lys8-11,pex17-ha

ATCC 201249

(Papanikola ou S 2002)

YL-DS

ATCC 201249-pINA1269-DS

This study

YL-OpDS0

ATCC 201249-pINA1269-OpDS

This study

YL-OpDS1

ATCC 201249 rDNA :: TEF1p-OpDS-TXPR2t

This study

YL-PPDS0

ATCC 201249 rDNA :: TEF1p-OpDS-TXPR2t, This study EXP1p-PPDS-MIG1t, GPD1p-ATR1-LIP2t

YL-PPDS1

YL-PPDS0-pINA1269-ERG1

This study

YL-PPDS2

YL-PPDS0-pINA1269-ERG8

This study

YL-PPDS3

YL-PPDS0-pINA1269-ERG9

This study

YL-PPDS4

YL-PPDS0-pINA1269-ERG10

This study

YL-PPDS5

YL-PPDS0-pINA1269-ERG12

This study

YL-PPDS6

YL-PPDS0-pINA1269-ERG13

This study

YL-PPDS7

YL-PPDS0-pINA1269-ERG19

This study

YL-PPDS8

YL-PPDS0-pINA1269-ERG20

This study

YL-PPDS9

YL-PPDS0-pINA1269-IDI

This study

YL-PPDS10

YL-PPDS0-pINA1269-HMG1

This study

YL-PPDS11

YL-PPDS0-pINA1269-tHMG1

This study

YL-PPDS12

YL-PPDS0-pINA1269-HE1

This study

29

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 30 of 37

YL-PPDS13

YL-PPDS0-pINA1269-HE9

This study

YL-PPDS14

YL-PPDS0-pINA1269-HE12

This study

YL-PPDS15

YL-PPDS0-pINA1269-HE19

This study

YL-PPDS16

YL-PPDS0-pINA1269-HE20

This study

YL-PPDS17

YL-PPDS0-pINA1269-HE9-20

This study

YL-PPDS18

YL-PPDS0-pINA1269-HE1-9

This study

YL-fuPPDS1

ATCC 201249 rDNA :: TEF1p-OpDS-TXPR2t, This study EXP1p-PPDS-linker1-ATR1-MIG1t

YL-fuPPDS2

ATCC 201249 rDNA :: TEF1p-OpDS-TXPR2t, This study EXP1p-PPDS-linker2-ATR1-MIG1t

YL-CK0

YL-PPDS0-pINA1269-UGT1

This study

YL-CK1

YL-CK0-pINA1269-HE9-20

This study

YL-MVA

ATCC 201249 zeta :: FBAINp-tHMG1-XPR2t , This study EXP1p-ERG9-MIG2t, GPD1p-ERG20-CYC1t

YL-MVA-CK

YL-MVA rDNA :: TEF1p-OpDS-XPR2t, EXP1p- This study PPDS-linker2-ATR1-MIG1t, GPD1p-UGT1-LIP2t

30

ACS Paragon Plus Environment

Page 31 of 37

Journal of Agricultural and Food Chemistry

Table 3. Metabolite analysis in different strains Strain

DMD

PPD

Squalene

2,3-oxidosqualene

Ergosterol

(mg/L)

(mg/L)

(mg/L)

(mg/L)

(mg/L)

YL-DS1

17.9±2.3

-

6.5±1.8

5.6±1.3

29.4±1.3

YL-PPDS0

4.8±0.6

12.9±1.3

5.9±1.4

4.8±0.7

27.1±2.7

YL-PPDS17 32.1±4.4

82.5±7.5

9.1±0.5

13.6±1.2

38.5±2.2

31

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 1. Biosynthesis pathway of CK in Yarrowia lipolytica.

32

ACS Paragon Plus Environment

Page 32 of 37

Page 33 of 37

Journal of Agricultural and Food Chemistry

Figure 2. Metabolites analysis of different engineered strains.

33

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 3. Production of PPD by metabolic engineering.

34

ACS Paragon Plus Environment

Page 34 of 37

Page 35 of 37

Journal of Agricultural and Food Chemistry

Figure 4. The biosynthesis of CK in engineered Y. lipolytica.

35

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 5. Fed-batch fermentation of the engineering strain YL-MVA-CK in 5-L bioreactor.

36

ACS Paragon Plus Environment

Page 36 of 37

Page 37 of 37

Journal of Agricultural and Food Chemistry

For Table of Contents Only

37

ACS Paragon Plus Environment