BioMed Research International
Volume 2016 (2016), Article ID 9845927, 8 pages
http://dx.doi.org/10.1155/2016/9845927

Research Article

Expression of Codon-Optimized Plant Glycosyltransferase UGT72B14 in Escherichia coli Enhances Salidroside Production

Feiyan Xue,1,2 Huili Guo,1,2 Yingying Hu,1,2 Ran Liu,1,2 Lina Huang,1,2 Heshu Lv,1,2 Chunmei Liu,1,2 Mingfeng Yang,1,2 and Lanqing Ma1,2,3

¹Key Laboratory of Urban Agriculture (North), Ministry of Agriculture, Beijing 102206, China
²College of Biological Science and Engineering, Beijing University of Agriculture, Beijing 102206, China
³Beijing Collaborative Innovation Center for Eco-Environmental Improvement with Forestry and Fruit Trees, Beijing 102206, China

Received 6 May 2016; Revised 18 July 2016; Accepted 21 July 2016

Academic Editor: Yeo J. Yoon

Copyright © 2016 Feiyan Xue et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Salidroside, a plant secondary metabolite in Rhodiola, has been demonstrated to have several adaptogenic properties as a medicinal herb. Due to the limitation of plant source, microbial production of salidroside by expression of plant uridine diphosphate glycosyltransferase (UGT) is promising. However, glycoside production usually remains hampered by poor expression of plant UGTs in microorganisms. Herein, we achieved salidroside production by expression of Rhodiola UGT72B14 in Escherichia coli (E. coli) and codon optimization was accordingly applied. UGT72B14 expression was optimized by changing 278 nucleotides and decreasing the G+C content to 51.05% without altering the amino acid sequence. The effect of codon optimization on UGT72B14 catalysis for salidroside production was assessed both in vitro and in vivo. In vitro, salidroside production by codon-optimized UGT72B14 is enhanced because of a significantly improved protein yield (increased by 4.8-fold) and an equivalently high activity as demonstrated by similar kinetic parameters ( and ), compared to that by wild-type protein. In vivo, both batch and fed-batch cultivation using the codon-optimized gene resulted in a significant increase in salidroside production, which was up to 6.7 mg/L increasing 3.2-fold over the wild-type UGT72B14.

1. Introduction

Glycosylation of small molecules has a profound impact on their biological properties including bioactivity, solubility, and stability. The majority of glycosylation reactions in nature are mediated by glycosyltransferases, which comprise a large and divergent polyphyletic multigene superfamily [1, 2]. In plants, glycosylation catalysed by uridine diphosphate (UDP) glycosyltransferases (UGTs) [2–4] is involved in the synthesis of many secondary metabolites. Studies on the recombinant expression of UGTs have resulted in the production of various plant secondary metabolites as glycosides [2, 5, 6]. However, many valuable plant glycosides remain difficult to produce in high yield and in a timely manner.

Microbes, especially Escherichia coli and Saccharomyces cerevisiae, can be engineered to produce plant natural products by introducing key enzymes [7, 8]. E. coli is the most widely used, well-characterized heterologous host and proved to result in more products yield than S. cerevisiae [4, 9, 10]. Notably, poor expression of wild-type plant genes in E. coli has restricted the yield of products [11, 12]. One strategy for overcoming this limitation is codon optimization [12–14].

Salidroside (8-O-β-D-glucoside of tyrosol) is one of the main plant secondary metabolites in Rhodiola and has been demonstrated to have adaptogenic and ergogenic capacity [15]. Consequently, demand for salidroside compounds has grown, but Rhodiola has a long growth period and low salidroside content, and an alternative method of salidroside production is thus needed. Glycosylation typically occurs as one of the last steps during natural product synthesis [5], and this is the case for salidroside [16, 17]; the last step involves glycosylation of tyrosol, catalysed by UGT (Figure 1). In our previous work, three different UGT genes, UGT72B14, UGT73B6, and UGT74R1, were isolated from Rhodiola sachalinensis and characterized [16, 17]. Among the three UGTs, UGT72B14 demonstrated the highest activity for salidroside production, with a catalytic efficiency (/) of 12.3, which was 620% and 180% higher than UGT74R1 and UGT73B6, respectively [17]. Accordingly, we investigated the effect of codon optimization on the heterologous expression of UGT72B14 in E. coli and assessed the impact on salidroside production in the present study.

Figure 1: Reaction catalysed by UGT during salidroside production.

2. Methods

2.1. Strains and Plasmids

The host E. coli (DE3) strain was from a frozen stock stored at the Key Laboratory of Urban Agriculture (North) of the Ministry of Agriculture of China, Beijing, China. This strain was used for salidroside production following transformation with the appropriate plasmid. The pET-28a (+) plasmid was purchased from Novagen (Darmstadt, Germany) and used for cloning and expression of UGT72B14 in the E. coli host strain.

2.2. Codon Optimization

The wild-type UGT72B14 gene was isolated from Rhodiola sachalinensis (GenBank Accession Number EU567325) and optimized by replacing codons predicted to be less frequently used in E. coli with more favored codons, according to the Codon Usage Database and graphical codon usage analyzer. The optimized gene with introduction of restriction sites (Nco I and Xho I) was synthesized by Inovogen Tech. Co. (Beijing, China).

2.3. Protein Expression

The synthetic gene was digested and ligated into the prokaryotic expression vector pET-28a (+) at the same restriction sites (Nco I and Xho I), under the control of T7 promoter. The resulting plasmid was transformed into the E. coli host. Following identification of a positive transformant, cells were grown at 37°C with shaking at 180 rpm in LB (Luria-Bertani) medium containing 50 μg/mL kanamycin. When the desired OD₆₀₀ was reached, the temperature was reduced to 20°C, shaking was slowed to 120 rpm, and expression was induced with isopropyl β-D-1-thiogalactopyranoside (IPTG) at a final concentration of 0.8 mM. Cultures were harvested at various time points to measure OD₆₀₀ and do sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE) analysis. Cell pellets, obtained by centrifugation at 12,857 ×g for 5 min, were resuspended in sterile H₂O and subjected to SDS-PAGE analysis.

2.4. Protein Purification and Measurement of Enzyme Activity

Recombinant E. coli cultures were incubated at 37°C with shaking at 180 rpm in LB medium containing 50 μg/mL kanamycin. When the OD₆₀₀ reached , E. coli harbouring the optimized UGT72B14 plasmid was induced at 20°C with shaking at 120 rpm as described above for a further 6–8 h. Cell pellets were harvested by centrifugation (8,000 ×g; 10 min; 4°C) and resuspended in lysis buffer (20 mM Tris, pH7.5, 300 mM NaCl, 20 mM imidazole, 1 mg/mL leupeptin, 2 mg/mL aprotinin, and 100 mg/mL lysozyme). Cells were disrupted by sonication (30 × 3 s) with ice bath. Supernatant was separated from the sonication homogenate by centrifugation (10,000 ×g; 10 min; 4°C) and passed through Ni-NTA His-Bind™ Resin (Novagen, Madison, WI, USA) column containing Ni²⁺ as an affinity ligand. After washing with elution buffer (20 mM Tris, pH7.5, 300 mM NaCl, and 20 mM imidazole), unbound contaminant proteins were removed. His-tagged proteins were eluted with washing buffer (20 mM Tris, pH7.5, 300 mM NaCl, and 300 mM imidazole) and purified using a PD-10 column (Amersham Pharmacia Biotech, Uppsala, Sweden) and then concentrated by a centrifugal filter unit (YM-50 kD, Millipore). The efficiency of purification was monitored by SDS-PAGE analysis. The concentration of purified UGT72B14 was determined using BCA protein assay kit (Sigma-Aldrich, USA), and then the purified protein yield was assessed by calculating the purified protein content per litre of culture broth. For recombinant E. coli harbouring the wild-type UGT72B14 plasmid, induction was initiated at OD₆₀₀ = and growth continued for a further 9 h. Cells were then harvested and protein was purified as described above. To measure the kinetic parameters, tyrosol (Sigma-Aldrich, Saint Louis, Missouri, USA) and UDP-glucose substrates were incubated with purified protein at 30°C for 30 min, and the salidroside product was identified using High Performance Liquid Chromatography (HPLC). Enzyme assays were carried out as previously described [17].

2.5. Measuring Salidroside Production in the Recombinant E. coli Strains

Cultures were grown and induced as described above, and tyrosol, glucose, and sodium citrate (all analytical grades) were added in equimolar concentrations and culturing continued at 30°C with shaking at 150 rpm. Samples were harvested at various time-points during fermentation by centrifugation at 12,857 ×g for 10 min. Supernatants were analyzed directly by HPLC or freeze-dried and extracted with an equal volume of ethanol, evaporated, and dissolved in methanol to prepare Liquid Chromatography and Mass Spectrometry (LC-MS) samples as described previously [17].

3. Results

3.1. Codon Optimization and Its Effect on Protein Expression

UGT72B14 was found to have the highest activity for salidroside production and was thus considered the ideal candidate for testing the effect of codon optimization. The full-length gene sequence (1671 bp) has a G+C content of 56.75% and includes an open reading frame (ORF) of 1422 bp that encodes a polypeptide of 473 amino acids [17]. Codons such as AGA (Arg), AGG (Arg), CGA (Arg), CGG (Arg), AUA (Ile), CUA (Leu), CUC (Leu), and CUU (Leu), present in the wild-type plant sequence, are used at less than 10% in E. coli. Bioinformatics analysis was performed to identify codons that could be modified to better match the host codon usage preferences without altering the resulting amino acid sequence. We replaced most of them with those that are more frequently used in E. coli BL21 (DE3). Codons encoding Arg, Gly, and Leu were altered in 93%, 67%, and 81% of cases, respectively. The termination codon was also modified as UAA. It is well documented that ORFs containing high G+C content are often poorly expressed in A+T rich hosts [14, 18]. After altering 278 nucleotides, the G+C content was decreased to 51.05% in the optimized sequence. The codon-optimized UGT72B14 was submitted to NCBI (GenBank Accession Number KU523897).

To study the effect of codon optimization on UGT72B14 expression, the optimized gene was synthesized and used to construct the recombinant plasmid that was subsequently transformed into E. coli, and protein expression was assessed by SDS-PAGE. The cell density at the time of IPTG induction is believed to be important for protein expression, and an optical density of 600 nm (OD₆₀₀) of 0.6–0.8 was previously demonstrated to be optimum for expression of both UGT73B6 and UGT72B14 [16, 17, 19]. This cell density was thus used for induction of expression of the optimized UGT72B14 gene in the present study. Considerable cell density (Figure 2(a)) and maximum accumulation of the target protein 9 h after induction (Figure 2(b)) were observed, and the variation of cell growth and protein expression was consistent with our previous work [17]. However, expression levels of the codon-optimized UGT72B14 were even lower than the wild-type gene when induction was performed at OD₆₀₀ ~ 0.6 (Figure 2(b)). Induction was thus tested at OD₆₀₀ ~ 0.1, 0.3, 0.6, and 0.9, and OD₆₀₀ ~ 0.3 was found to result in the highest expression of the codon-optimized UGT72B14 gene (Figure 3(a)). SDS-PAGE analysis showed that induction at OD₆₀₀ ~ 0.3 and growth for a further 6–8 h resulted in peak target protein levels (Figure 3(d)), which was overall significantly superior to that induction at OD₆₀₀ ~ 0.6 (Figure 3(c)), though lower in cell growth (Figure 3(b)). To our best knowledge, this interesting result is firstly observed in Rhodiola UGT expression. However, it is consistent with a previous report of other proteins [20] that IPTG addition decreased cell growth but enhanced expression of the target gene when added at a relatively low cell density (OD₆₀₀ ~ 0.3).

(a)

(b)

Figure 2: Growth profile and protein expression of the recombinant strains harbouring codon-optimized and wild-type UGT72B14 gene. (a) Cell growth of recombinant strains. Data are presented as mean of three biological replicates, and error bar represents the standard deviation (SD). (b) SDS-PAGE analysis of protein expression. Sample of 20 μL was used for SDS-PAGE analysis. UGT72B14-opt = codon-optimized UGT72B14; UGT72B14-wt = wild-type UGT72B14. Lane M: protein molecular weight markers. Lanes 1–3: protein expression at 4, 6, and 9 h after induction. Target protein is indicated by arrows.

(a)

(b)

(c)

(d)

Figure 3: Optimization of UGT72B14-opt expression. (a) Protein expression at 10 h after induction with IPTG at varying OD₆₀₀. (b) Growth profile of recombinant strain during induction. (c) Time courses of protein expression with IPTG induction at OD₆₀₀ = . (d) Time courses of protein expression with IPTG induction at OD₆₀₀ = . Lane M: protein molecular weight markers. Lanes I–IV: UGT72B14-opt expression at OD₆₀₀ ~ 0.1, 0.3, 0.6, and 0.9 with 10 h after induction. Lanes 1–7: UGT72B14-opt expression at 0, 2, 4, 6, 8, 9, and 10 h after induction. Cell pellets were resuspended, the OD₆₀₀ was adjusted to , and 15 μL was used for SDS-PAGE analysis. Target protein is indicated by arrows. Data are presented as mean of three biological replicates, and error bar represents the SD.

3.2. Effect of Codon Optimization on Salidroside Production In Vitro

To examine whether codon-optimized UGT72B14 are capable of converting tyrosol to salidroside, the His-tagged UGT72B14 recombinant proteins were purified using Ni affinity chromatography. Result of SDS-PAGE analysis indicated that the codon-optimized gene resulted in not only more target protein expression but also more efficient purification (Figure 4(a)). It was 27.55 mg/L of purified protein obtained from the culture of strain harbouring the codon-optimized UGT72B14, which was an increase of 4.8-fold over that from the wild-type gene (Figure 4(b)). To investigate the effect of codon optimization on salidroside production in vitro, the activity of codon-optimized and wild-type UGT72B14 protein was tested using tyrosol and UDP-glucose as substrates, and production of salidroside was monitored by HPLC and LC-MS. There was no significant difference in specific activity between codon-optimized and wild-type enzymes, as demonstrated by comparable kinetic parameters ( and ; Figure 4(b)). Wild-type UGT72B14 has been already proved to exhibit highest level of activity for salidroside production, compared to the other two UGTs we have cloned [17]. Here, an equivalently high enzyme activity and a significant improvement in protein yield of codon-optimized UGT72B14 will offer great promise for salidroside production in vitro.

(a)

(b)

Figure 4: Result of the recombinant protein yield and salidroside production in vitro. (a) SDS-PAGE analysis of the recombinant protein expression and purification. Lane M: protein molecular weight markers. Lane 1: total protein in uninduced strain. Lane 2: total protein in induced strain. Lane 3: concentrated protein after being purified by Ni affinity chromatography and PD-10 column. Sample of 15 μL was used for SDS-PAGE analysis. Target protein is indicated by arrows. (b) Result of salidroside production catalysed by UGT72B14 in vitro. Protein yield data of the optimized UGT72B14 was obtained by purified protein after induction of 7 h initiated with OD₆₀₀ = while that of the wild-type UGT72B14 was obtained by purified protein after induction of 9 h initiated with OD₆₀₀ = . The recombinant enzyme catalysis reaction system (100 μL) contained 50 mM Tris-HCl (pH 7.5), 2 mM UDP-glucose, 250 μM tyrosol, and the enzyme protein of 0.2 mg, proceeded for 30 min at 30°C, and terminated with 200 μL MeOH. Data were presented as mean ± SD.

3.3. Effect of Codon Optimization on Salidroside Production In Vivo

Utilizing recombinant enzymes to synthesize plant secondary metabolites in E. coli has many advantages, including simple nutritional requirements, easy cultivation, and the ability to prevent enzyme inactivation [6, 7, 21, 22]. The plant secondary metabolite salidroside is synthesized via a one-step glycosylation reaction from tyrosol in Rhodiola [16, 17, 23]. Following confirmation of codon-optimized UGT72B14 expression and catalytic activity in vitro, substrates were added to growing cultures at 6 h after induction with IPTG and conversion of tyrosol to salidroside was investigated. Tyrosol was tested at 0, 5, 50, and 500 mg/L, and most salidroside was observed at 50 mg/L tyrosol (Figure 5). This tyrosol concentration was subsequently used for a time-course experiment on cultures harbouring codon-optimized or wild-type UGT72B14 plasmids. As shown in Figure 6(a), salidroside was first detected at 2 h after addition of tyrosol, and the concentration gradually increased, peaked at 8 h, and then declined over the following 2 h. A similar pattern was observed for both recombinant strains, but the codon-optimized culture exhibited a much wider range. Codon optimization increased the maximum salidroside production by 3-fold. Time-course experiments revealed a tendency for the target protein yield to decline at the end of batch fermentation, which was consistent with previous reports [18, 24]. We tentatively conclude that batch cultivation resulted in substrate inadequate utilization at the early stages of fermentation, and the low product yield encouraged us to investigate other ways of improving salidroside production.

Figure 5: Effect of tyrosol concentration on salidroside production in E. coli expressing codon-optimized UGT72B14. Data are presented as mean of three biological replicates, and error bar represents the SD. Different letters above bars indicate a significant difference at according to Duncan’s test.

(a)

(b)

Figure 6: Result of salidroside production in vivo. (a) Batch cultivation. (b) Fed-batch cultivation. Sal-opt and Sal-wt refer to salidroside synthesized by strains harbouring codon-optimized and wild-type UGT72B14 plasmids, respectively. Tyr-opt and Tyr-wt refer to tyrosol concentration in the culturing broth of strains harbouring codon-optimized and wild-type UGT72B14 plasmids, respectively. Black arrows indicate substrate additions, and red arrows indicate substrate feeding during cultivation. Samples were tested three times. Data are presented as mean, and error bar represents the SD.

Fed-batch cultivation has been shown to be an effective way to enhance the expression of target proteins and the yield of target products in E. coli [25–28]. In the present study, we performed fed-batch experiments by feeding cultures with four additions of 50 mg/L tyrosol (an initial addition of 20 mg/L, and three more additions of 10 mg/L at 2 h intervals). In these experiments, salidroside was first detected in the culture broth at 2 h after the first tyrosol addition and was then maintained at a high level (with slight fluctuations) between 4 and 11 h after the first addition (Figure 6(b)). Up to 6.7 mg/L of salidroside was accumulated at 9 h in the codon-optimized culture, which was an increase of 3.2-fold over the culture harbouring the wild-type UGT72B14 plasmid. Compared with salidroside production in the batch cultivation method, fed-batch fermentation achieved a higher titre, and salidroside production was stable over a longer period. A gradual supplement of tyrosol is thus beneficial for salidroside production, and optimizing the amount of tyrosol added could further increase salidroside yield. This may explain the higher yield of salidroside that was achieved following reconstitution of the tyrosol biosynthetic pathway in E. coli [19].

4. Discussion

Glycosyltransferases have proved to be very important for plant secondary metabolism synthesis and have become an important research field in improving catalytic efficiency [3, 29]. In terms of UGTs, they have emerged as promising catalyst for UDP-sugar based glycosylation focusing on small molecules like salidroside [4]. Although many UGTs have been cloned and expressed in E. coli glycoside production remains hampered by poor expression of plant genes in this prokaryote host [11]. Codon optimization, as an effective approach to improving expression of heterologous protein, has attracted considerable attention and was applied successfully in some cases [12–14, 17]. However, no systematic studies of codon optimization have been reported for plant UGTs.

The UGTs of Rhodiola sachalinensis as one of the most important plants for salidroside supply [30, 31] thus aroused our great interest. Three Rhodiola UGTs were firstly cloned and characterized in our previous work [17]. With the sequence of UGT73B6, one of the three UGTs we published, salidroside has been synthesized in vivo by Bai et al. [19]. Here with the sequence of UGT72B14, codon optimization was systematically studied and proved favorable for enhancing salidroside production both in vivo and in vitro.

5. Conclusions

The optimized UGT72B14 expression at a high level was benefited from induction with IPTG at a reasonably low cell density (OD₆₀₀ ~ 0.3). A fed-batch cultivation method proved optimal for salidroside production in terms of overall yield and the duration over which maximum production occurred. As demonstrated favorable for enhancing UGT72B14 expression and salidroside production in a timely manner in E. coli, codon optimization therefore deserved to be used to improve plant UGTs expression in heterologous hosts to enhance their corresponding glycoside production.

Competing Interests

The authors declare that they have no competing interests.

Authors’ Contributions

Feiyan Xue and Huili Guo contributed equally to this work and performed most of the experiments and drafted the paper. Yingying Hu, Ran Liu, and Lina Huang participated in selection of positive clones, protein expression, and cultivation runs. Heshu Lv and Mingfeng Yang contributed to some detailed experiments designing. Chunmei Liu helped in doing HPLC-MS analysis. Lanqing Ma coordinated the whole project and completed the paper. All authors read and approved the final version of the paper.

Acknowledgments

This research was financially supported by “Beijing Natural Science Foundation (2164059),” “Beijing Outstanding Talent Project (2013D005021000003),” “Funding Project for Scientific Research Quality Improvement in Beijing University of Agriculture (GJB2013006, GZL2015012),” “Beijing Municipal Education Commission (CEFF-PXM2016_014207_000038),” “National Natural Science Foundation of China (31300620, 31370674),” and “Foundation of Beijing Municipal Education Committee (KM20140020001).” The authors also thank Mr. Mingcheng Zhang and Jianguo Yang for technical support.

References

J. Ross, Y. Li, E.-K. Lim, and D. J. Bowles, “Higher plant glycosyltransferases,” Genome Biology, vol. 2, article 3004, 2001.

C. M. M. Gachon, M. Langlois-Meurinne, and P. Saindrenan, “Plant secondary metabolism glycosyltransferases: the emerging functional analysis,” Trends in Plant Science, vol. 10, no. 11, pp. 542–549, 2005.

T. Vogt and P. R. Jones, “Glycosyltransferases in plant natural product synthesis: characterization of a supergene family,” Trends in Plant Science, vol. 5, no. 9, pp. 380–386, 2000.

F. De Bruyn, J. Maertens, J. Beauprez, W. Soetaert, and M. De Mey, “Biotechnological advances in UDP-sugar based glycosylation of small molecules,” Biotechnology Advances, vol. 33, no. 2, pp. 288–302, 2015.

K. S. Hansen, C. Kristensen, D. B. Tattersall et al., “The in vitro substrate regiospecificity of recombinant UGT85B1, the cyanohydrin glucosyltransferase from Sorghum bicolor,” Phytochemistry, vol. 64, no. 1, pp. 143–151, 2003.

F. Bourgaud, A. Gravot, S. Milesi, and E. Gontier, “Production of plant secondary metabolites: a historical perspective,” Plant Science, vol. 161, no. 5, pp. 839–851, 2001.

L. B. Pickens, Y. Tang, and Y.-H. Chooi, “Metabolic engineering for the production of natural products,” Annual Review of Chemical and Biomolecular Engineering, vol. 2, pp. 211–236, 2011.

J. D. Keasling, “Manufacturing molecules through metabolic engineering,” Science, vol. 330, no. 6009, pp. 1355–1358, 2010.

J. Beekwilder, R. Wolswinkel, H. Jonker, R. Hall, C. H. De Rie Vos, and A. Bovy, “Production of resveratrol in recombinant microorganisms,” Applied and Environmental Microbiology, vol. 72, no. 8, pp. 5670–5672, 2006.

E. Gruchattka, O. Hädicke, S. Klamt, V. Schütz, and O. Kayser, “In silico profiling of Escherichia coli and Saccharomyces cerevisiae as terpenoid factories,” Microbial Cell Factories, vol. 12, article 84, 2013.

V. J. J. Martin, D. J. Piteral, S. T. Withers, J. D. Newman, and J. D. Keasling, “Engineering a mevalonate pathway in Escherichia coli for production of terpenoids,” Nature Biotechnology, vol. 21, no. 7, pp. 796–802, 2003.

H. Liu, J. Li, G. Du, J. Zhou, and J. Chen, “Enhanced production of α-cyclodextrin glycosyltransferase in Escherichia coli by systematic codon usage optimization,” Journal of Industrial Microbiology and Biotechnology, vol. 39, no. 12, pp. 1841–1849, 2012.

H. G. Menzella, “Comparison of two codon optimization strategies to enhance recombinant protein production in Escherichia coli,” Microbial Cell Factories, vol. 10, article 15, 2011.

B. Maertens, A. Spriestersbach, U. Von Groll et al., “Gene optimization mechanisms: a multi-gene study reveals a high success rate of full-length human proteins expressed in Escherichia coli,” Protein Science, vol. 19, no. 7, pp. 1312–1326, 2010.

H.-M. Chiang, H.-C. Chen, C.-S. Wu, P.-Y. Wu, and K.-C. Wen, “Rhodiola plants: chemistry and biological activity,” Journal of Food and Drug Analysis, vol. 23, no. 3, pp. 359–369, 2015.

L.-Q. Ma, B.-Y. Liu, D.-Y. Gao et al., “Molecular cloning and overexpression of a novel UDP-glucosyltransferase elevating salidroside levels in Rhodiola sachalinensis,” Plant Cell Reports, vol. 26, no. 7, pp. 989–999, 2007.

H.-S. Yu, L.-Q. Ma, J.-X. Zhang, G.-L. Shi, Y.-H. Hu, and Y.-N. Wang, “Characterization of glycosyltransferases responsible for salidroside biosynthesis in Rhodiola sachalinensis,” Phytochemistry, vol. 72, no. 9, pp. 862–870, 2011.

Y. J. Jeong, S. G. Woo, C. H. An et al., “Metabolic engineering for resveratrol derivative biosynthesis in Escherichia coli,” Molecules and Cells, vol. 38, no. 4, pp. 318–326, 2015.

Y. Bai, H. Bi, Y. Zhuang et al., “Production of salidroside in metabolically engineered Escherichia coli,” Scientific Reports, vol. 4, article 6640, 2014.

B. A. Fong and D. W. Wood, “Expression and purification of ELP-intein-tagged target proteins in high cell density E. coli fermentation,” Microbial Cell Factories, vol. 9, article 77, 2010.

A. Schmid, J. S. Dordick, B. Hauer, A. Kiener, M. Wubbolts, and B. Witholt, “Industrial biocatalysis today and tomorrow,” Nature, vol. 409, no. 6817, pp. 258–268, 2001.

M. C. Song, E. J. Kim, E. Kim, K. Rathwell, S.-J. Nam, and Y. J. Yoon, “Microbial biosynthesis of medicinally important plant secondary metabolites,” Natural Product Reports, vol. 31, no. 11, pp. 1497–1509, 2014.

X. Lan, K. Chang, L. Zeng et al., “Engineering salidroside biosynthetic pathway in hairy root cultures of Rhodiola crenulata based on metabolic characterization of tyrosine decarboxylase,” PLoS ONE, vol. 8, no. 10, Article ID e75459, 2013.

S. Liu, Q. Qi, N. Chao et al., “Overexpression of artificially fused bifunctional enzyme 4CL1-CCR: a method for production of secreted 4-hydroxycinnamaldehydes in Escherichia coli,” Microbial Cell Factories, vol. 14, no. 1, article 118, 2015.

A. Bahrami, S. A. Shojaosadati, R. Khalilzadeh, and E. V. Farahani, “Two-stage glycerol feeding for enhancement of recombinant hG-CSF production in a fed-batch culture of Pichia pastoris,” Biotechnology Letters, vol. 30, no. 6, pp. 1081–1085, 2008.

M. Krause, K. Ukkonen, T. Haataja et al., “A novel fed-batch based cultivation method provides high cell-density and improves yield of soluble recombinant proteins in shaken cultures,” Microbial Cell Factories, vol. 9, article no. 11, 2010.

J. Li, J. Jaitzig, P. Lu, R. D. Süssmuth, and P. Neubauer, “Scale-up bioprocess development for production of the antibiotic valinomycin in Escherichia coli based on consistent fed-batch cultivations,” Microbial Cell Factories, vol. 14, article 83, 2015.

C. Wang, J. Zhang, H. Wu, Z. Li, and Q. Ye, “Heterologous gshF gene expression in various vector systems in Escherichia coli for enhanced glutathione production,” Journal of Biotechnology, vol. 214, pp. 63–68, 2015.

J. Härle, S. Günther, B. Lauinger et al., “Rational design of an Aryl-C-glycoside catalyst from a natural product O-glycosyltransferase,” Chemistry and Biology, vol. 18, no. 4, pp. 520–530, 2011.

S. Nakamura, X. Li, H. Matsuda et al., “Bioactive constituents from Chinese natural medicines. XXVI. Chemical structures and hepatoprotective effects of constituents from roots of Rhodiola sachalinensis,” Chemical and Pharmaceutical Bulletin, vol. 55, no. 10, pp. 1505–1511, 2007.

S. Wu, Y. Zu, and M. Wu, “High yield production of salidroside in the suspension culture of Rhodiola sachalinensis,” Journal of Biotechnology, vol. 106, no. 1, pp. 33–43, 2003.