Production of shikimic acid from Escherichia coli through chemically inducible chromosomal evolution and cofactor metabolic engineering
Microbial Cell Factories volume 13, Article number: 21 (2014)
Shikimic acid (SA) produced from the seeds of Chinese star anise (Illicium verum) is a key intermediate for the synthesis of neuraminidase inhibitors such as oseltamivir (Tamiflu®), an anti-influenza drug. However, plants cannot deliver a stable supply of SA. To avoid the resulting shortages and price fluctuations, a stable source of affordable SA is required. Although recent achievements in metabolic engineering of Escherichia coli strains have significantly increased SA productivity, commonly-used plasmid-based expression systems are prone to genetic instability and require constant selective pressure to ensure plasmid maintenance. Cofactors also play an important role in the biosynthesis of different fermentation products. In this study, we first constructed an E. coli SA production strain that carries no plasmid or antibiotic marker. We then investigated the effect of endogenous NADPH availability on SA production.
The pps and csrB genes were first overexpressed by replacing their native promoter and integrating an additional copy of the genes in a double gene knockout (aroK and aroL) of E. coli. The aroG fbr, aroB, aroE and tktA gene cluster was integrated into the above E. coli chromosome by direct transformation. The gene copy number was then evolved to the desired value by triclosan induction. The resulting strain, E. coli SA110, produced 8.9-fold more SA than did the parental strain E. coli (ΔaroK ΔaroL). Following qRT-PCR analysis, another copy of the tktA gene under the control of the 5Ptac promoter was inserted into the chromosome of E. coli SA110 to obtain the more productive strain E. coli SA110. Next, the NADPH availability was increased by overexpressing the pntAB or nadK genes, which further enhanced SA production. The final strain, E. coli SA116, produced 3.12 g/L of SA with a yield on glucose substrate of 0.33 mol/mol.
An SA-producing E. coli strain that carries neither a plasmid nor an antibiotic marker was constructed by triclosan-induced chromosomal evolution. We present the first demonstration that increasing NADPH availability by overexpressing the pntAB or nadK genes significantly enhances SA production.
Shikimic acid is a key intermediate for the synthesis of the neuraminidase inhibitor oseltamivir (Tamiflu®), an anti-influenza treatment . The main commercial source of SA is seeds of the Illicium plant, such as I. verum or I. anistatum. However, the conventional method of producing SA from I. verum is typically low-yield and costly. Therefore, researchers have developed several metabolic engineering approaches to overproduce SA in E. coli[2, 3]. These approaches are based on genetic modifications of central carbon metabolism and SA pathways. Specifically, the aromatic amino acid pathway is blocked after the SA production stage by deleting the aroK and aroL genes encoding shikimate kinase I and II. To increase the carbon flux from central carbon metabolism entering the aromatic amino acid pathway, researchers have amplified the feedback resistant 3-deoxy-D-arabino heptulosonate 7-phosphate (DAHP) synthase genes aroF fbr or aroG fbr[4–8]. These modifications are commonly complemented with over-expression of the aroB and aroE genes [4–8]. Over-expression of the tktA gene (encoding transketolase) enhances the availability of erythrose 4-phosphate (E4P), and consequently increases the SA titer from 38 g/L to 52 g/L . Over-expression of the pps gene (encoding phosphoenolpyruvate (PEP) synthase) has elevated the SA titer to 66 g/L . SA production can also be increased by inactivating the phosphoenolpyruvate:carbohydrate phosphotransferase system (PTS) operon, combined with expressing the ATP dependent uptake and phosphorylation system comprising the glucose facilitator and the glucokinase from Zymomonas mobili[7, 8]. E. coli engineered in this way yielded SA concentrations as high as 87 g/L in minimum medium supplemented with yeast extract . The constitutive and synchronous expression of a six-gene synthetic operon (aroG fbr, aroB, aroD, aroE, tktA and zwf), in a laboratory-evolved strain bearing simultaneous PTS and pykF inactivations, was recently reported to increase the SA yield on glucose to 42% mol/mol, which represents the highest reported yield .
However, all of the above studies use plasmids for gene expression. Among the drawbacks of plasmid-based expression systems are structural and segregational instability, and allele segregation [10–12]. These plasmid instabilities cause genetic instability, with decreased productivity of the desired compound. Tyo et al.  reported that plasmid-carrying strains lose poly-3-hydroxybutyrate productivity after 40 rounds of subculturing. Moreover, the markers used for selecting and maintaining plasmids in hosts during cultivation are usually antibiotic resistance genes. However, antibiotics are both costly and banned from food and pharmaceutical production processes. In addition, the potential spread of antibiotic-resistant marker to natural microbes requires serious consideration. The likely outcome of this scenario is rapid emergence of multidrug-resistant organisms (e.g., superbacteria) [14, 15]. To overcome these drawbacks of plasmid constructs, Tyo et al.  developed a plasmid-free method that achieves high copy numbers of the desired genes, termed chemically induced chromosomal evolution (CIChE). However, the λInCh genomic integration protocol used in CIChE is complicated and time-consuming, and Tyo et al.’s CIChE strains still carry an antibiotic resistance marker (chloramphenicol resistance) . In our previous paper , we reported a modified CIChE that overcomes the drawbacks of Tyo et al’s original CIChE protocol. The resulting strain, developed by triclosan-induced chromosomal evolution, carries neither a plasmid nor an antibiotic marker. Thus, our first task in this study was to construct an SA-producing E. coli strain by triclosan-induced chromosomal evolution.
Cofactors are also known to play an important role in the biosynthesis of different fermentation products. Once the enzyme levels are no longer limiting, cofactor availability can become limiting and productivity must be boosted by cofactor manipulation. In silico flux analysis has identified a potentially important role for intracellular NADPH concentration in SA biosynthesis . Thus, we also investigated the effect of endogenous NADPH availability on SA production by over-expressing the genes involved in NADPH synthesis.
Results and discussion
Inactivation of genes encoding shikimate kinase I and II, and chromosomal promoter replacement
To construct an SA-accumulating host strain, we initially focused on preventing the conversion of SA to chorismic acid in E. coli. To this end, we deleted the aroK and aroL genes to obtain E. coli BW25113 (ΔaroKΔaroL). PEP synthase (encoded by the pps gene) converts pyruvic acid to PEP. Carbon storage regulator CsrA controls glycogen synthesis and modulates the levels of three enzymes that directly participate in PEP metabolism: pyruvate kinase, PEP carboxykinase and PEP synthase (which synthesize PEP from oxaloacetate and pyruvate, respectively). The first of these is positively regulated, while both PEP enzymes are negatively regulated, by CsrA [18, 19]. Because CsrA activity is antagonized by CsrB , PEP availability was increased by replacing the native promoters of the pps and csrB genes with the P lacQ1 promoter in E. coli BW25113 (ΔaroKΔaroL), yielding E. coli BW25113 (ΔaroKΔaroL, P pps ::P lacQ1 , P csrB ::P lacQ1 ). SA production in this strain was 53.3% higher than in E. coli BW25113 (ΔaroKΔaroL) ((0.23 ± 0.02) g/L vs (0.15 ± 0.04) g/L; P < 0.01; Table 1). To further improve PEP availability, an additional copy of the pps and csrB genes controlled by the T5 promoter was inserted into E. coli BW25113 (ΔaroKΔaroL, P pps ::P lacQ1 , P csrB ::P lacQ1 ) to obtain E. coli BW25113 (ΔaroKΔaroL, P pps ::P lacQ1 , P csrB ::P lacQ1 , P T5 -pps, P T5 -csrB). SA production in this strain was 73.3% higher than in E. coli BW25113 (ΔaroKΔaroL) ((0.26 ± 0.02) g/L vs (0.15 ± 0.04) g/L; P < 0.01; Table 1). These modifications also resulted in an increment in yield (mol SA/mol glucose) of about 80%. Over-expressing the pps gene increased DAHP production almost twofold, to approach the theoretical maximum . DAHP is the precursor of SA and the condensation product of PEP and E4P, mediated by DAHP synthase. They thought that over-expressing the pps gene drives pyruvate to be recycled back to PEP . Yakandawala et al.  achieved an approximately twofold increase in phenylalanine production by over-expressing the csrB gene. SA is the precursor of phenylalanine. Chandran et al. reported that over-expressing the pps gene raised the SA titer of 26% (from 52 to 66 g/L) . Thus, the increments in SA titer and yield could be the consequence of overexpression of both pps and csrB genes.
Quantitative real-time PCR analysis demonstrates that the above modifications indeed enhanced the transcription levels of the pps and csrB genes (Figure 1). The replacement of the native promoter resulted in the upregulations of the transcription levels of the pps and csrB genes about 6.9- and 2.8-fold (P < 0.01), respectively. Integrating an additional copy of the two genes further upregulated the transcription levels of the pps and csrB genes about 10.4- and 20.9-fold (P < 0.01), respectively.
Chemically induced chromosomal evolution
Many papers reported that plasmid-based over-expression of feedback-resistant DAHP synthases (coded by aroF fbr or aroG fbr), shikimate dehydrogenase (coded by aroE), transketolase (ecoded by tktA), and DHQ synthase enzymes (ecoded by aroB) [2–9]. To overcome the drawbacks of plasmid expression systems, the aroG fbr, tktA, aroB and aroE gene cluster was inserted into the chromosome of E. coli BW25113 (ΔaroKΔaroL, P pps ::P lacQ1 , P csrB ::P lacQ1 , P T5 -pps, P T5 -csrB) by direct transformation. The transformed strain was evolved to higher gene copy number by exposure to increasing triclosan concentrations. Figure 2 plots the SA production in CIChE strains resistant to different triclosan concentrations. The maximum SA concentration (1.34 ± 0.15) g/L was obtained by the CIChE strains resistant to 2 μM triclosan. When the triclosan concentration exceeded 2 μM, SA production of the CIChE strains did not further increase. Thus, the recA gene of the CIChE strain resistant to 2 μM triclosan was deleted to obtain E. coli SA110. Homologous recombination, which can potentially reduce the copy number, is inhibited in this strain. Figure 3 shows that the copy numbers and transcription levels of the aroE gene of the CIChE strains increase with triclosan concentration during chromosomal evolution. At a triclosan concentration of 2 μM, the copy number reached about 8 in the CIChE strain. When the triclosan concentration was above 2 μM, the gene copy numbers and transcription levels still increased; however, the SA production of the CIChE strains did not increase. The transcription levels of the other three genes show the similar trend (Additional file 1: Figure S1). The results indicated that there is an optimal copy number and transcription level of the gene cluster for efficient production of SA.
We also analyzed the transcription levels of the E. coli SA110 genes by qRT-PCR, and compared the data with those of E. coli BW25113 (ΔaroKΔaroL) (Figure 4). The transcription levels of the tktA, aroG fbr , aroB and aroE genes in E. coli SA110 were about 1.6, 14.7, 35.2 and 8.1 times higher (p < 0.05) than in E. coli BW25113 (ΔaroKΔaroL). The transcription level of the tktA gene was the lowest. It may be because the distance of the tktA gene from the promoter is the longest . They thought that the abundance of mRNA decreased monotonically with the increasing distance of the gene from the promoter in E. coli. The gene ranked in front will be translated primarily, followed by the translation of the subsequent genes along with the mRNA transcription for polycistronic operons in E. coli. However, the transcription level of the second gene (the aroB gene) from the promoter was the highest. This could be the consequence of codon usage biases.
Because the tktA gene was least upregulated among these genes, we examined whether tktA is the bottleneck for SA production in E. coli SA110. An additional copy of the tktA gene under the control of the 5Ptac promoter was inserted into the chromosome of E. coli SA110 to obtain E. coli SA112. This modification raised the SA titer to (1.70 ± 0.01) g/L and the SA yield to 0.25 mol/mol (Table 1). Figure 4 shows that this integration upregurated the transcription level of the tktA gene. The transcription level of the tktA gene in E. coli SA112 was much higher than that in E. coli SA110. This may be because the additional integrated tktA gene was directly controlled by the stronger 5Ptac promoter.
Effect of NADPH availability
To investigate the effect of NADPH availability on SA production, transhydrogenase (encoded by sthA or pntAB) and NAD kinase (encoded by nadK) were first amplified by plasmid-based over-expression in E. coli SA112. The results are presented in Table 2. Overexpression of the pntAB or nadK genes increased both SA production and the intracellular NADPH concentration. E. coli generates NADPH in one of three ways: the pentose phosphate (PP) pathway, isocitrate dehydrogenase and transhydrogenase. Two isoforms of transhydrogenase exist in E. coli, membrane-bound transhydrogenase (PntAB) and soluble transhydrogenase (SthA). During aerobic batch growth with glucose, 35–45% of the NADPH required for biosynthesis is produced via PntAB, while the PP pathway and isocitrate dehydrogenase generate 35–45% and 20–25%, respectively . ATP-dependent NAD kinase encoded by nadK catalyzes the phosphorylation of NAD to NADP. Overexpression of nadK increases the size of the NADP pool, which potentially increases the abundance of NADPH. It can be seen in Table 2 that overexpression of the sthA gene increased SA titer and reduced the intracellular NADPH concentration. Moreover, the SA titer and intracellular NADPH concentration of the strain overexpressing sthA was lower than in strains overexpressing the pntAB or nadK genes. PntAB catalyzes the transfer of reducing power from NADH to NADP+ in an energy-dependent manner at low intracellular NADPH levels, while SthA mainly catalyzes the reoxidation of NADPH in an energy-independent manner when NADPH is abundant . Thus, if the pntAB or sthA gene is overexpressed, the intracellular NADPH concentration is increased or decreased, respectively. As shown in Table 3, overexpressing the sthA gene indeed caused a reduction in intracellular NADPH concentration, with consequent reduction of the SA titer.
Because plasmid-based overexpression of sthA showed lower the SA titer compared to overexpression of the pntAB or nadK genes, to alleviate the metabolic burden caused by the plasmid, the pntAB and nadK genes were integrated into the chromosome of E. coli SA112 to obtain E. coli SA114 and SA116, respectively. Chromosomal overexpression enhanced the intracellular NADPH concentration and SA titer relative to plasmid overexpression (Table 2), and also improved the cell growth. It can be seen in Figure 5 that the transcription levels of the pntAB and nadK genes were enhanced about 3.3- and 3.6-fold (P < 0.01) after chromosomal integration, respectively.
This study demonstrated a strong correlation between NADPH availability and SA production. In the SA biosynthesis pathway, shikimate dehydrogenase (encoded by the aroE gene) catalyzes the reduction of 3-dehydroshikimate to shikimate, which requires NADPH as a cofactor. These CIChE strains (e.g. SA110, 112, 114, 116), in which the aroE gene was overexpressed and therefore required more NADPH for SA biosynthesis. Although NADPH can be provided by the oxidative part of the PP pathway, this branch is unfavorable for SA production. Stoichiometric analysis shows that the oxidative part of the PP pathway is not required in the reaction scheme. Glucose enters the nonoxidative part of the PP pathway, where it provides E4P for maximum theoretical SA yield . Accordingly, the overall reaction for SA biosynthesis from glucose is: 7 glucose + 6 ATP + 6 NAD + 6 NADPH → 6 SA + 6 ADP + 6 NADH + 6 NADH + 6 NADP + 12H+ + 6 HPO4 2+. Increased intracellular NADPH via overexpression of the pntAB or nadK genes favored SA production (Table 2), suggesting that more carbon flux was channeled into the nonoxidative part of the PP pathway for E4P biosynthesis. In silico flux analysis has demonstrated a potentially important role for intracellular NADPH concentration in SA biosynthesis. Under different genetic and environmental conditions, E. coli cells consume the same quantities of NADPH to maximize their SA production . According to some studies, maintaining NADPH availability improves metabolite production. For instance, overexpressing E. coli pntAB in Corynebacterium glutamicum enhances NADPH availability, in turn increasing the intracellular levels of L-lysine  and L-ornithine . Overexpression of nadk in E. coli increases the NADPH/NADP ratio, thereby enhancing thymidine biosynthesis . In C. glutamicum, overexpressing the NAD kinase gene improves L-lysine  and L-ornithine  production. Simultaneous chromosomal overexpression of transhydrogenase (pntAB) and NAD kinase (yfjB) genes in E. coli increases the NADPH supply and improves anaerobic isobutanol production .
The SA yield obtained in the present study is below the highest value (0.42 mol/mol) reported in the literature , achieved by the constitutive and synchronous expression of a six-gene synthetic operon, in a laboratory-evolved strain bearing simultaneous PTS and pykF inactivations in a 1 L bioreactor containing 100 g/L of glucose and 30 g/L of yeast extract. However, the SA yield obtained in this research is on part with the 33% yield reported by Chandran et al. . They reported the highest SA titer (87 g/L). It remains far from the theoretical maximum (86% mol/mol) , implying a need for further research to improve this value. Currently, we are exploring whether inactivation of PTS, pykF and the oxidative part of the PP pathway ensures a high SA yield. The results of these modifications will be reported in a subsequent paper.
To overcome the drawbacks of plasmid-based expression systems, the aroG fbr, aroB, aroE and tktA gene cluster was integrated into the E. coli chromosome by direct transformation. The gene copy number was then evolved to the desired value by triclosan induction. Following qRT-PCR analysis, SA production was further enhanced by inserting a second copy of the tktA gene under the control of the 5Ptac promoter into the chromosome of the CIChE strain. The effect of NADPH availability on SA production was also investigated. NADPH availability and SA production were found to be strongly correlated. Plasmid-based or chromosomal overexpression of the pntAB or nadK genes enhanced the intracellular NADPH concentration and consequently the SA titer. This is the first report of an engineered SA-producing strain of E. coli that lacks both a plasmid and an antibiotic marker. Using this strain, no resistance-conferring compound was required during the fermentation process.
Strains, plasmids and primers
The strains and plasmids used in this study are listed in Table 4. E. coli DH5α was used for plasmid construction. The parent strain for SA production was E. coli BW2511 . The primers used in this study are listed in Table 5.
Gene knockout and integration
Gene knockouts and replacement of the native promoter were carried out by PCR product recombination  using the pSIM6 plasmid  expressing the lambda red recombination system and pKD3  as the template for PCR. Gene knockouts were verified by colony PCR using appropriate primers (listed in Table 5).
Chemically induced chromosomal evolution
The aroE, aroB and aroG genes were amplified from K12 genomic DNA using the corresponding primer pairs (Table 5) and cloned into pMD18-T to obtain pMD-aroE, pMD-aroB and pMD-aroG, respectively. The feedback-inhibition-resistant (fbr) aroG gene, which contained an Asp-146-Asn substitution, was obtained by site-directed mutagenesis using the MutanBEST kit (Takara, Dalian, China), following the manufacturer’s instructions. The aroE and aroB fragments were digested by their corresponding restriction enzymes and inserted into the Eco RI/Nco I sites of pBAD24  to form pBEB. The aroG fbr fragment was digested from pMD-aroGfbr and inserted into the Pst I/Sph I sites of pBEB to obtain pBEBG. The aroG fbr , aroB and aroE gene cluster was excised from pBEBG using restriction enzymes Eco RI and Sph I. The resulting gene cluster was cloned into the Eco RI/Sph I sites of pHKKF3T5b  to obtain pHKEBG. The tktA gene was amplified from E. coli using the primers TP1 and TP2. The fragment was digested along with pHKEBG using Sph I and Sac I, and ligated together to form pHKEBGT. Assisted by a helper plasmid, pAH69, which expresses the phage integrase , the resulting integration vector pHKEBGT was inserted into the bacterial attachment (attB) site of E. coli by direct transformation, as described by Chen et al.  and Chiang et al. .
CIChE of the above construct was carried out by subculturing the resulting strains in 5 mL Super Optimal Broth (SOB) medium with increasing concentrations of triclosan in 15 mL culture tubes, as described by Chen et al.  and Tyo et al. . The strains were grown to stationary phase in 1 μM triclosan. Fifty milliliters of the culture was subcultured into a new culture tube, in which the triclosan concentration was doubled from 1 to 2 μM and allowed to grow to stationary phase. The process was repeated until the desired concentration (as high as 32 μM) was reached. The recA gene of the CIChE strain was then deleted.
The multiple cloning site (ATGCATGACGTCGGGCCCGCATGCCACGTGGAGCTCGGTACCATAAAACGAAAGGCTCAGTCGAAAGACTGGGCCTTTCGTTTTATCAATTGCTGCAGCCCGGGCTCGAGTCTAGAGTCGACCCGCGG) was synthesized by Takara Biotechnology. This fragment was digested along with p5TG  using Nsi I/Sac II, and ligated with p5GT to form pMP5. The sthA, pntAB and nadK genes were amplified from E. coli using the corresponding primer pairs (Table 5) and then inserted into the corresponding sites of pMP5 to form pMPsthA, pMPpntAB and pMPnadK, respectively.
To enhance the expression of the genes on the integrate vector, the T5 promoter of the integrate expression plasmid was replaced with the 5Ptac promoter (five tandem repeats of the core-tac-promoter) . The 5Ptac promoter was amplified using primers 5TacF and 5TacR from p5TG . The resulting fragment was digested along with pHKKT5b or pP21KT5b [33, 34] using Mlu I and Sph I, and was ligated into these plasmids to form pHKK5Tacb or pP21K5Tacb, respectively. The tktA gene was excised from pHKEBGT by Sph I and Sac I and ligated into the Sph I/Sac I sites of pP21K5Tacb to obtain the integrate vector pP21K5Tacb-tktA. The pntAB and nadK genes were excised from pMPpntAB and pMPnadk by their corresponding restriction enzymes and then inserted into the corresponding sites of pP21K5Tacb to form the integrate vectors pP21K5Tacb-pntAB and pP21K5Tacb-nadK, respectively.
Assisted by a helper plasmid, pAH121, expressing the phage integrase , the above integrate vectors (pP21KT5b-csrB-pps, pP21K5Tacb-tktA, pP21K5Tacb-pntAB and pP21K5Tacb-nadK) were inserted into the bacterial attachment (attB) site of E. coli by direct transformation, as described by Chen et al.  and Chiang et al. .
Quantitative real-time PCR (qRT-PCR)
Total RNA from E. coli cells grown for 60 h in shake flasks was isolated using an RNA extraction kit (Dongsheng Biotech, Guangzhou, China), following the manufacturer’s instructions. The first-strand cDNA was synthesized using an All-in-One™ First-Strand cDNA Synthesis Kit (GeneCopoeia, Guangzhou, China). The qRT-PCR was performed with the All-in-One™ qPCR Mix kit (GeneCopoeia) on an iCycler iQ5 Real Time PCR system (Bio-Rad Laboratories, California, USA). The template was 100 ng of cDNA. The PCR conditions were: 95°C for 10 min, followed by 45 cycles of denaturation at 95°C for 10 s, annealing at 60°C for 20 s, and extension at 72°C for 15 s. The primers for qRT-PCR are presented in Table 4. Data were analyzed by the 2-ΔΔCt method described by Livak and Schmittgen , and normalized by cysG gene expression.
Gene copy numbers were measured by qPCR on genomic DNA isolated from the appropriate CIChE strains. qPCR was performed as above. The primers QEF and QER (Table 5) were used to measure the copy number of aroE.
For SA production, E. coli cells were precultured overnight in a falcon tube containing 5 mL Luria Broth (LB) at 37°C. The main cultures were incubated in the fermentation medium (pH 7.0) containing (g/L): glucose 10, peptone 1, proline 1.24, KH2PO4 3, Na2HPO4 · 7H2O 13, MgCl2·6H2O 0.24, NaCl 0.5, NH4Cl 1, CaCl2 0.1. The main cultures were inoculated with a starting OD600 of 0.1 and incubated at 37°C for 72 h in a rotary shaking incubator at 150 rpm.
Cell growth was measured by optical density at 600 nm and converted to dry cell weight (DCW; g/L) using a standard curve. SA concentration was determined by HPLC using a Shimadzu system (LC-20A, Shimadzu, Japan) equipped with an Inertsil ODS-SP column (5 μm, 4.6 × 150 mm, GL Sciences Inc, Tokyo, Japan). The mobile phase was 5 mmol/L H2SO4, with a flow rate of 0.5 mL/min, at 50°C. SA was detected by a photodiode array detector (SPD-M20A) operating at 210 nm, and quantified by a standard curve constructed from serial dilutions of a SA standard stock solution (1 mg/ml, J&K Scientific Ltd, Beijing, China). Glucose concentration was determined by glucose oxidase using a glucose assay kit (Shanghai Rongsheng Biotech Corporation, Shanghai, China).
Following aerobic cultivation of E. coli on a rotary shaker (150 rpm) at 37°C for 54 h, the cells were harvested by centrifugation and washed twice with water. Intracellular NADPH was extracted and quantified using the Enzychrom™ NADP+/NADPH Assay kit (BioAssay Systems, Hayward, CA), following the manufacturer’s instructions.
All experiments were conducted in triplicate, and data were averaged and presented as the mean ± standard deviation. Significant differences were determined by one-way analysis of variance followed by Tukey’s test, using the OriginPro (version 7.5) package. Statistical significance was defined as p < 0.05.
Kim CU, Lew W, Williams MA, Liu H, Zhang L, Swaminathan S, Bischofberger N, Chen MS, Mendel DB, Tai CY, Laver WG, Stevens RC: Influenza neuraminidase inhibitors possessing a novel hydrophobic interaction in the enzyme active site: design, synthesis, and structural analysis of carbocyclic sialic acid analogs with potent anti-influenza activity. J Am Chem Soc. 1997, 119: 681-690. 10.1021/ja963036t.
Krämer M, Bongaerts J, Bovenberg R, Kremer S, Müller U, Orf S, Wubbolts M, Raeven L: Metabolic engineering for microbial production of shikimic acid. Metab Eng. 2003, 5: 277-283. 10.1016/j.ymben.2003.09.001.
Ghosh S, Chisti Y, Banerjee Uttam C: Production of shikimic acid. Biotechnol Adv. 2012, 30: 1425-1431. 10.1016/j.biotechadv.2012.03.001.
Escalante A, Calderon R, Valdivia A, de Anda R, Hernández G, Hernández OT, Gosset G, Bolivar F: Metabolic engineering for the production of shikimic acid in an evolved Escherichia coli strain lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system. Microb Cell Fact. 2010, 9: 21- 10.1186/1475-2859-9-21.
Draths KM, Knop DR, Frost JW: Shikimic acid and quinic acid: replacing isolation from plant sources with recombinant microbial biocatalysis. J Am Chem Soc. 1999, 121: 1603-1604. 10.1021/ja9830243.
Knop DR, Draths KM, Chandran SS, Barker JL, Frost JW: Hydroaromatic equilibrium during biosynthesis of shikimic acid. J Am Chem Soc. 2001, 123: 10173-10182. 10.1021/ja0109444.
Chandran SS, Yi J, Draths KM, Von Daeniken R, Weber W, Frost JW: Phosphoenolpyruvate availability and the biosynthesis of shikimic acid. Biotechnol Prog. 2003, 19: 808-814. 10.1021/bp025769p.
Yi J, Draths KM, Li K, Frost JW: Altered glucose transport and shikimate pathway product yields in E. coli. Biotechnol Prog. 2003, 19: 1450-1459. 10.1021/bp0340584.
Rodriguez A, Martínez JA, Báez-Viveros JL, Flores N, Hernández-Chávez G, Ramírez OT, Gosset G, Bolivar F: Constitutive expression of selected genes from the pentose phosphate and aromatic pathways increases the shikimic acid yield in high-glucose batch cultures of an Escherichia coli strain lacking PTS and pykF. Microb Cell Fact. 2013, 12: 86- 10.1186/1475-2859-12-86.
Bentley WE, Mirjalili N, Andersen DC, Davis RH, Kompala DS: Plasmid-encoded protein: the principal factor in the metabolic burden associated with recombinant bacteria. Biotechnol Bioeng. 1990, 35: 668-681. 10.1002/bit.260350704.
Noack D, Roth M, Geuther R, Muller G, Undisz K, Hoffmeier C, Gaspar S: Maintenance and genetic stability of vector plasmids pBR322 and pBR325 in Escherichia coli K12 strains grown in a chemostat. Mol Gen Genet. 1981, 184: 121-124. 10.1007/BF00271207.
O’Connor M, Peifer M, Bender W: Construction of large DNA segments in Escherichia coli. Science. 1989, 244: 1307-1312. 10.1126/science.2660262.
Tyo KEJ, Ajikumar PK, Stephanopoulos G: Stabilized gene duplication enables long-term selection-free heterologous pathway expression. Nat Biotechnol. 2009, 27: 760-765. 10.1038/nbt.1555.
Chiang C-J, Chen PT, Chao YP: Replicon-free and markerless methods for genomic insertion of DNAs in phage attachment sites and controlled expression of chromosomal genes in Escherichia coli. Biotechnol Bioeng. 2008, 101: 985-995. 10.1002/bit.21976.
Goh S, Good L: Plasmid selection in Escherichia coli using an endogenous essential gene marker. BMC Biotechnol. 2008, 8: 61- 10.1186/1472-6750-8-61.
Chen Y-Y, Shen H-J, Cui Y-Y, Chen S-G, Weng Z-M, Zhao M, Liu J-Z: Chromosomal evolution of Escherichia coli for the efficient production of lycopene. BMC Biotechnol. 2013, 13: 6- 10.1186/1472-6750-13-6.
Ahn J, Chung BKS, Lee D-Y, Park M, Karimi IA, Jung J-K, Lee H: NADPH-dependent pgi-gene knockout Escherichia coli metabolism producing shikimate on different carbon sources. FEMS Microbiol Lett. 2011, 324: 10-16. 10.1111/j.1574-6968.2011.02378.x.
Yakandawala N, Romeo T, Friesen AD, Madhyastha S: Metabolic engineering of Escherichia coli to enhance phenylalanine production. Appl Microbiol Biotechnol. 2008, 78: 283-291. 10.1007/s00253-007-1307-z.
Suzuki K, Babitzke P, Kushner SR, Romeo T: Identification of a novel regulatory protein (CsrD) that targets the global regulatory RNAs CsrB and CsrC for degradation by RNase E. Genes Dev. 2006, 20: 2605-2617. 10.1101/gad.1461606.
Patnaik R, Liao JC: Engineering of Escherichia coli central metabolism for aromatic metabolite production with near theoretical yield. Appl Environ Microbiol. 1994, 60: 3903-3908.
Nishizaki T, Tsuge K, Itaya M, Doi N, Yanagawa H: Metabolic engineering of carotenoid biosynthesis in Escherichia coli by ordered gene assembly in Bacillus subtilis. Appl Environ Microbiol. 2007, 73: 1355-1361. 10.1128/AEM.02268-06.
Sauer U, Canonaco F, Heri S, Perrenoud A, Fischer E: The soluble and membrane-bound transhydrogenases udhA and pntAB have divergent functions in NADPH metabolism of Escherichia coli. J Biol Chem. 2004, 279: 6613-6619.
Kabus A, Georgi T, Wendisch VF, Michael B: Expression of the Escherichia coli pntAB genes encoding a membrane-bound transhydrogenase in Corynebacterium glutamicum improves L-lysine formation. Appl Microbiol Biotechnol. 2007, 75: 47-53. 10.1007/s00253-006-0804-9.
Jiang LY, Chen SG, Zhang YY, Liu JZ: Metabolic evolution of Corynebacterium glutamicum for increased production of L-ornithine. BMC Biotechnol. 2013, 13: 47- 10.1186/1472-6750-13-47.
Lee HC, Kim JS, Jang W, Kim SY: High NADPH/NADP ratio improves thymidine production by a metabolically engineered Escherichia coli strain. J Biotechnol. 2010, 149: 24-32. 10.1016/j.jbiotec.2010.06.011.
Lindner SN, Niederholtmeyer H, Schmitz K, Schoberth SM, Wendisch VF: Polyphosphate/ATP-dependent NAD kinase of Corynebacterium glutamicum: biochemical properties and impact of ppnK overexpression on lysine production. Appl Microbiol Biotechnol. 2010, 87: 583-593. 10.1007/s00253-010-2481-y.
Shi A, Zhu X, Lu J, Zhang X, Ma Y: Activating transhydrogenase and NAD kinase in combination for improving isobutanol production. Metab Eng. 2013, 16: 1-10.
Datsenko KA, Wanner BL: One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci U S A. 2000, 97: 6640-6645. 10.1073/pnas.120163297.
Sharan SK, Thomason LC, Kuznetsov SG, Court DL: Recombineering: a homologous recombination-based method of genetic engineering. Nat Protoc. 2009, 4: 206-223. 10.1038/nprot.2008.227.
Guzman LM, Belin D, Carson MJ, Beckwith J: Tight regulation, modulation, and high-level expression by vectors containing the arabinose PBAD promoter. J Bacteriol. 1995, 177: 4121-4130.
Haldimann A, Wanner BL: Conditional-replication, integration, excision, and retrieval plasmid-host systems for gene structure-function studies of bacteria. J Bacteriol. 2001, 183: 6384-6393. 10.1128/JB.183.21.6384-6393.2001.
Li M, Wang J, Geng Y, Li Y, Wang Q, Liang Q, Qi Q: A strategy of gene overexpression based on tandem repetitive promoters in Escherichia coli. Microb Cell Fact. 2012, 11: 19- 10.1186/1475-2859-11-19.
Liu JZ, Huang MT, Cui YY, Chen YY: A series of expression plasmids for chromosomal integration and evolution. China patent. 2012, 201210060042-
Huang MT, Chen YY, Liu JZ: Chromosomal engineering of Escherichia coli for the efficient production of coenzyme Q10. Chin J Chem Eng. In press
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCt Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
We are grateful to both the National Natural Science Foundation of China (Grant No. 30970089, 20876181, 21276289), the Natural Science Foundation of Guangdong Province (No. S2011010001396) and the Key Project of the Scientific and Technical Innovation of Higher Education of Guangdong Province (2012CXZD0002) for their financial support.
The authors declare that they have no competing interests.
YYC carried out most of the experiments. CL amplified some of the genes. YYZ constructed some of the plasmids. JH deleted the aroK and aroL genes. JZL developed the concept and designed the method, led the project and drafted the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Cui, YY., Ling, C., Zhang, YY. et al. Production of shikimic acid from Escherichia coli through chemically inducible chromosomal evolution and cofactor metabolic engineering. Microb Cell Fact 13, 21 (2014). https://doi.org/10.1186/1475-2859-13-21