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
© Rodriguez et al.; licensee BioMed Central Ltd. 2013
Received: 21 June 2013
Accepted: 28 August 2013
Published: 30 September 2013
During the last two decades many efforts have been directed towards obtaining efficient microbial processes for the production of shikimic acid (SA); however, feeding high amounts of substrate to increase the titer of this compound has invariably rendered low conversion yields, leaving room for improvement of the producing strains. In this work we report an alternative platform to overproduce SA in a laboratory-evolved Escherichia coli strain, based on plasmid-driven constitutive expression of six genes selected from the pentose phosphate and aromatic amino acid pathways, artificially arranged as an operon. Production strains also carried inactivated genes coding for phosphotransferase system components (ptsHIcrr), shikimate kinases I and II (aroK and aroL), pyruvate kinase I (pykF) and the lactose operon repressor (lacI).
The strong and constitutive expression of the constructed operon permitted SA production from the beginning of the cultures, as evidenced in 1 L batch-mode fermentors starting with high concentrations of glucose and yeast extract. Inactivation of the pykF gene improved SA production under the evaluated conditions by increasing the titer, yield and productivity of this metabolite compared to the isogenic pykF+ strain. The best producing strain accumulated up to 43 g/L of SA in 30 h and relatively low concentrations of acetate and aromatic byproducts were detected, with SA accounting for 80% of the produced aromatic compounds. These results were consistent with high expression levels of the glycolytic pathway and synthetic operon genes from the beginning of fermentations, as revealed by transcriptomic analysis. Despite the consumption of 100 g/L of glucose, the yields on glucose of SA and of total aromatic compounds were about 50% and 60% of the theoretical maximum, respectively. The obtained yields and specific production and consumption rates proved to be constant with three different substrate concentrations.
The developed production system allowed continuous SA accumulation until glucose exhaustion and eliminated the requirement for culture inducers. The obtained SA titers and yields represent the highest reported values for a high-substrate batch process, postulating the strategy described in this report as an interesting alternative to the traditionally employed fed-batch processes for SA production.
Although the implementation of these modifications along with bioengineering strategies has led to diverse E. coli strains capable of accumulating SA, the yields obtained to date are still far from the theoretical maximum [10, 11, 14, 15]. This can be partially attributed to the fact that most expression systems used involve genes controlled by a mixture of inducible and native promoters of variable strengths, contained in more than one type of plasmid. These imbalances often cause a metabolic burden and heterogeneities on the intensity and temporality of gene expression, which may translate into suboptimal production capabilities of the recombinant strains, resulting in low productivity and yield of SA [16–19]. Consequently, optimized DNA expression systems and genetic backgrounds are needed for promoting a more efficient carbon channeling towards SA formation.
With the goal of producing aromatic compounds, our group has constructed and characterized strains lacking PTS, the major glucose transport system . One of such strains is PB11, which grows poorly on glucose due to the inactivation of PTS [21, 22]. Strain PB12, a derivative of PB11 with a 400% increased growth rate, was isolated in a short laboratory adaptive evolution process to foster derivatives growing in glucose [21, 22]. This strain can simultaneously utilize glucose and other carbon sources (acetate, glycerol and various carbohydrates) in minimal medium due to the lack of catabolite repression exerted by PTS [21, 23]. Whole genome analysis allowed the identification of the genetic changes that occurred in PB12, suggesting that the deletion of 12 genes, including rppH, galR and mutH, is the main reason for its rapid growth on glucose .
It was reported that PB12, which assimilates glucose by the non-PTS symporter GalP , can be engineered to accumulate SA in culture media containing glucose (Glc) and yeast extract (YE). For instance, when PB12 was transformed with two plasmids encoding four biosynthetic genes, the variant with both functional pyruvate kinases accumulated the highest SA concentration (up to 7 g/L), but the highest yield of aromatic compounds was achieved by a derivative with an inactivated pykF gene . This result may be related to other reported effects caused by the inactivation of pykF, such as an increase in plasmid copy number per cell , low acetate production due to less glycolytic overflux [13, 27, 28], or higher concentrations of the AAA pathway enzymes . Interestingly, in spite of the aforesaid features that can be beneficial for SA production, the metabolic engineering efforts to overproduce this compound have been mainly applied to strains with a pykF+ background, probably because of their typically higher glucose consumption rates compared to the pykF- counterparts [11, 30].
Results and discussion
Construction of strains derived from PB12 aroK-aroL- containing a plasmid designed for the constitutive expression of a synthetic operon used in the production of shikimic acid
Unpublished evidence from our laboratory indicates that the production of aromatic compounds in the laboratory-evolved strain PB12 can attain higher levels when the transcriptional induction of the genes involved in canalizing carbon flux into the AAA pathway occurs at the beginning of fermentations. Taking into account this observation, a new strategy was developed for optimizing the production of SA in PB12 carrying inactive aroK and aroL genes (Figure 1). This strategy included the design and construction of a plasmid for the strong and stable expression of six key genes arranged in the form of a synthetic operon, controlled exclusively by a single Trc promoter. In order to reduce metabolic burden, a single plasmid derived from pBR327 carrying the par locus for increased plasmid stability was utilized as the vector , after incorporating a fragment containing the promoter, polylinker, and transcriptional terminators from pTrc99A (Figure 2).
Escherichia coli strains and plasmids utilized in this report
F´ tra D36 proA + proB + lacI q lac ZΔM15/supE thi Δ(lac-proAB) rpoS(33 am)
JM101 Δ(ptsH, ptsI, crr)::kan
PB11, PTS- Glc+; laboratory-evolved strain
AR2 + pTrcAro6 (Trc/aroB + tktA + aroGfbr + aroE + aroD + zwf + )
AR3 + pTrcAro6 (Trc/aroB + tktA + aroGfbr + aroE + aroD + zwf + )
AR2 + pTrc327par (plasmid vector without synthetic operon)
AR3 + pTrc327par (plasmid vector without synthetic operon)
PCR template for amplification of chloramphenicol resistance gene flanked by homologous sequences
Datsenko and Wanner 
Plasmid expressing λ-Red recombinase system with thermosensitive origin of replication
Datsenko and Wanner 
FLP recombinase expression plasmid
Cherepanov and Wackernagel 
Derivative of pBR322 exhibiting increased copy number and segregational stability
Zurita et al. 
Multipurpose expression plasmid bearing lacI gene and a polylinker in front of Trc promoter
Amann et al. 
Contains the promoter, polylinker, and terminators of pTrc99A, and par and ori regions of pBR327par
This work (Figure 2)
pTrc327par containing a 6-gene synthetic operon to enhance the production of shikimate
This work (Figure 2)
The spatial arrangement of the coding sequences that constitute the synthetic operon in pTrcAro6, flanked by the Trc promoter and transcriptional terminators, is shown in Figure 2. aroB is the first gene in the operon since several evidences indicate that its low expression is one of the limiting steps in the production of aromatic compounds [33–35]. Plasmid pTrcAro6 also carries the tktA and aroGfbr genes, whose products are involved in E4P synthesis and its condensation with PEP to form DAHP, the first aromatic compound (Figure 1). aroD and aroE genes were also included to promote an efficient conversion of DHQ to SA. Additionally, this plasmid carries the zwf gene, coding for the first enzyme of the PPP (Figure 1). The decision to include this gene was based on the following observations: 1) the overexpression of zwf substantially recovered the growth rate loss due to plasmid metabolic load in strain JM101 growing on glucose as only carbon source ; 2) it has been reported that strain PB12 displays a particularly low carbon flux partition at the glucose 6-phosphate (G6P) node towards the PPP (5% of the consumed G6P compared to 22% in the parental strain JM101) . Therefore, an overexpression of this gene should increase NADPH availability, required in catalytic amounts by the enzyme shikimate dehydrogenase (AroE), and may alleviate potential growth affectations by redirecting more G6P towards nucleotide and amino acid biosynthesis in strains derived from PB12 . However, the experiments presented in this report did not aim to dissect the specific effect of any utilized gene but instead sought to characterize the consequences of expressing all of them as an operon.
In order to promote an efficient translation of every gene, each coding sequence was amplified using designated primers that introduced a consensus Shine-Dalgarno sequence located 8 bp upstream of the translation start site. The nucleotide sequence of the constructed operon is presented in Additional file 1.
Assessment of the effects caused by pykF inactivation in strains expressing the Aro6 operon
To evaluate the effects caused by pykF inactivation on the production of SA, the performance of production strains AR26 (pykF+) and AR36 (pykF-) was compared using shake flasks containing 15 g/L of Glc and 5 g/L of YE. As a control, the same strains containing an empty pTrc327par plasmid (without the Aro6 operon), AR2e and AR3e, were also included.
To determine if the higher acetate production and lower SA production in AR26 compared to AR36 is a consequence of the inherently low oxygen availability and acidification of the medium in shake flask cultures, both strains were cultured in 1 L batch fermentors under controlled conditions of pH and dissolved oxygen tension (DOT). As an approach to increase the SA titer, the initial concentration of Glc in these experiments was raised to 100 g/L, and the YE concentration was concomitantly increased to 15 g/L to allow higher biomass generation.
Comparative data from 1 L batch fermentations of strains AR26 and AR36, using 100 g/L of Glc and 15 g/L of YE as substrates
SA titer (g/L)
12.95 ± 0.64
41.80 ± 2.83
Glc consumed (g/L)
82.65 ± 4.88
103.70 ± 6.79
Duration of culture (h)
Y SA/Glc (mol/mol)
0.16 ± 0.02
0.42 ± 0.00
Acetate titer (g/L)
29.35 ± 0.21
11.90 ± 0.14
X max (g/L)
6.18 ± 0.10
6.54 ± 0.09
μ (h -1 )
0.45 ± 0.01
0.45 ± 0.02
Qp global (gSA/L*h)
0.27 ± 0.02
0.75 ± 0.07
Qs global (gGlc/L*h)
-1.76 ± 0.10
-1.73 ± 0.11
It is remarkable that such large differences in acetate and SA production were observed by disrupting only one gene, which demonstrates the advantages of the combined inactivation of PTS and pykF when using a constitutive expression system in an evolved E. coli strain. To account for the observed improvements in SA production, we suggest that the early and constant expression of enzymes encoded in the operon could maintain a steady consumption of glycolytic intermediates throughout the cultures, preventing high fluctuations in their intracellular concentrations. We hypothesize that the combination of this steady metabolic state with a reduced flux from PEP to pyruvate caused by the inactivation of the pykF gene may increase the availability of PEP and other glycolytic precursors for SA production without decreasing the Glc consumption rate. However, we acknowledge that in the absence of measured intracellular metabolite concentrations, these remarks are speculative.
Fermentation profiles of AR36 in batch cultures
Taking into account the previous results, AR36 was selected for further characterization of its kinetic and stoichiometric performance in 1 L fermentors. To accomplish such purpose, the production of SA was tested with three different high-substrate culture conditions. Growth, Glc and byproducts were measured for each case, which in turn allowed a comparison of the productivities and yields.
Comparison of measured metabolites and calculated kinetic and stoichiometric parameters between three fermentations of strain AR36 with different substrate concentrations
Batch 50 g/L Glc + 15 g/L YE
Batch 100 g/L Glc + 15 g/L YE
Batch 100 g/L Glc + 30 g/L YE
SA titer (g/L)
23.80 ± 0.00
41.80 ± 2.83
43.30 ± 0.57
Glc consumed (g/L)
52.65 ± 1.20
103.70 ± 6.79
105.55 ± 4.45
Duration of culture (h)
Y SA/Glc (mol/mol)
0.47 ± 0.01
0.42 ± 0.00
0.42 ± 0.01
Acetate titer (g/L)
1.45 ± 0.00
11.90 ± 0.14
8.65 ± 0.92
X max (g/L)
6.30 ± 0.09
6.54 ± 0.09
12.54 ± 0.06
μ (h -1 )
0.53 ± 0.03
0.45 ± 0.02
0.45 ± 0.00
Qp global (gSA/L*h)
0.74 ± 0.00
0.75 ± 0.07
1.44 ± 0.02
Qs global (gGlc/L*h)
-1.65 ± 0.04
-1.74 ± 0.11
-3.52 ± 0.15
qp exp (gSA/gDCW*h)
0.38 ± 0.00
0.34 ± 0.03
0.46 ± 0.06
qs exp (gGlc/gDCW*h)
-1.14 ± 0.15
-1.11 ± 0.00
-1.18 ± 0.03
qp sta (gSA/gDCW*h)
0.16 ± 0.01
0.19 ± 0.02
0.20 ± 0.01
qs sta (gGlc/gDCW*h)
-0.36 ± 0.01
-0.42 ± 0.01
-0.42 ± 0.03
The effect of increasing the YE on SA productivity was investigated with a third set of experiments, using 100 g/L of Glc and 30 g/L of YE. Although the biomass was doubled when using twice the concentration of YE, the SA titer, μ and YSA/Glc were very similar to those obtained in the culture with 100 g/L of Glc and 15 g/L of YE (Figure 4b and Figure 4c). In conjunction with data obtained from the other two conditions, these findings suggest that the amount of YE primarily determines the maximum biomass that can be achieved. Additionally, an increment in the initial YE concentration did not alter the SA titer, and supports the observation that SA is mainly being produced from glucose. The direct relation between the initial YE concentration and the maximum biomass generated, regardless of the initial Glc concentration tested in these growth conditions, suggests that one or more limiting nutrients are being supplied by the YE. It would also appear that such nutrients cannot be synthesized from Glc, hence their depletion from YE limits growth long before Glc is exhausted. It is expected that the aromatic amino acids and vitamins present in the YE that are needed to counteract AR36 auxotrophy will become limiting; however, other compounds in this complex media may also play a role in growth limitation over time.
For a starting YE concentration of 30 g/L, a total of 106 g/L of Glc and 43 g/L of SA were consumed and produced, respectively, in approximately half the time than the fermentation with 15 g/L of YE. With 30 g/L of YE, the Qsglobal and Qpglobal increased twofold, in comparison with the fermentations with 15 g/L of YE, even though the SA titer remained unchanged (Table 3). Since the biomass also increased twofold, the calculated qp and qs were similar between the three experiments, both in exponential and stationary phases, exhibiting the metabolic robustness of the engineered strain under the tested conditions.
In addition, the results showed that an increase in YE concentration did not increase considerably the concentration of SA pathway intermediates (Figure 5c). In this respect, it has been acknowledged that the presence of high quantities of pathway intermediates can negatively impact the recovery of SA from the fermentation broth [7, 38]. This concern has directed some efforts into the subject, leading to the testing of culture conditions, genetic backgrounds, and the use of non-metabolizable glucose analogs, as attempts to minimize byproduct generation .
In these experiments, a high proportion of SA relative to byproducts was detected without applying any further modification to the strain or process. The concentration of each pathway intermediate was compared against the sum of all aromatic intermediates, and their percentages were used to calculate the molar ratio of SA to each byproduct at the end of the fermentations (Figure 6). The ratio of SA turned out to be higher than 10 for DHS, QA, or DAHP, and higher than 40 for GA or DHQ for all the substrate concentrations tested. Remarkably, in all the conditions the obtained SA yields were close to 50% of the theoretical maximum and the yields of total aromatic compounds (TAC) were above 60% of the theoretical maximum, estimated as 0.86 molTAC/molGlc (see Methods and Figure 6). This reflects the efficient redirection of glucose towards the AAA pathway in strain AR36, even when using high-glucose batch cultures. The ratio of SA to byproducts, as well as the obtained SA and TAC yields are fairly constant for all the conditions evaluated, and represent to our knowledge the highest reported values for a SA production fermentation process. These improvements can be justified by taking into account that the platform present in the engineered strain allows a more homogeneous expression of the necessary enzymes on an efficient genetic background. This, in contrast with other expression systems where the required genes are expressed from separate plasmids, under different promoters, or in strains not optimized for efficient use of high levels of Glc. In addition to the advantages concerning the dynamics of gene expression, the fact that IPTG is not needed to induce the Aro6 operon represents an important economical benefit for the production process, since the high price of IPTG restricts its use in large-scale fermentations.
Insights on the glycolytic and acetate metabolisms of strain AR36 by RT-qPCR
To gain a deeper insight of the metabolic changes induced by the constitutive expression of the Aro6 synthetic operon in strain AR36, transcript levels of several genes were measured at three different growth stages in cultures with 50 g/L of Glc and 15 g/L of YE. As detailed in Methods, data obtained from early exponential phase (EE), late exponential phase (LE), and stationary phase (ST) were normalized against the values measured from strain AR3e at EE, grown under the same culture conditions.
On the other side, the transcriptional levels of genes coding for enzymes involved in acetate biosynthesis (poxB, ackA and pta) were not modified by the presence of the synthetic operon, while actP and acs, coding for enzymes involved in acetate assimilation, were strongly upregulated in the EE and LE phases (Figure 7b). Upregulation of actP and acs genes has also been detected in the exponential growth phase in the parental strain PB12 that is capable of co-utilizing Glc and acetate in minimal medium . These findings correlate with the low levels of acetate in the assayed growth condition (Figure 4a). Importantly, the transcriptional values of these genes involved in acetate assimilation were low in ST phase (Figure 7b). If this response is representative of the other growth conditions used, it could partially explain the acetate accumulation observed in fermentations with 100 g/L of Glc, which consume higher amounts of Glc during stationary phase (Figure 4b and Figure 4c). These results highlight actP and acs as potential gene targets to artificially increase their expression in late culture stages, taking advantage of the expected capabilities of strain AR36 to utilize simultaneously Glc and acetate, present in its parental strain PB12 [21, 40].
The genes present in the synthetic operon showed very strong expression levels (even in stationary phase), reflecting the constitutive nature of the promoter and high copy number of the plasmid (Figure 7c). These results correlate with the uninterrupted Glc consumption and SA production observed during the entire fermentation (Figure 4a), suggesting that the enzymes coded by the genes in the operon are present throughout the cultivation time. It can be seen in Figure 7c that the transcript levels of aroD and zwf are comparatively higher and lower, respectively, than the other four genes in the operon. This observation should be taken with caution because the six genes in the operon are being compared to the ones present in the chromosome of reference strain AR3e. Since the values obtained for the six genes are not normalized between them, variations amongst their chromosomal expression in strain AR3e can alter the relative comparisons with strain AR36. Nevertheless, the transcriptomic data is consistent with the high ratio of SA to aromatic intermediates obtained in the tested conditions, which is to be expected if all the genes in the operon were adequately expressed. Together with kinetic and stoichiometric data, these results highlight the benefits of employing a constitutively-expressed synthetic operon as an alternate strategy to increase the yield of SA from Glc in an evolved strain that lacks PTS and pykF.
E. coli is the microorganism that has given the best results for SA production, with engineered strains that can accumulate up to 85 g/L using 10 L fermentors in fed-batch processes . In this report, we showed that the constitutive and synchronous expression of a six-gene synthetic operon, in a laboratory-evolved strain bearing simultaneous PTS and pykF inactivations, resulted in a competitive process for the production of SA. The expression of Aro6 operon in the PTS-pykF- derivative resulted in higher Qpglobal and similar Qsglobal than its PTS-pykF+ counterpart. In addition, the Glc consumption and SA production profiles of strain AR36 are consistent with the observed increase in transcription levels of glycolytic genes as a response to the constitutive expression of the operon in this strain. These features translated into significant improvements in growth and production parameters in strain AR36 (producing 0.74 gSA/L*h with 54% yield, using 50 g/L of Glc + 15 g/L of YE), compared to the PTS-pykF- strain reported by Escalante et al. in 2010 (producing 0.11 gSA/L*h with 26% yield, using 25 g/L of Glc + 15 g/L YE).
Albeit fed-batch fermentations in the past have given the best results with respect to SA production, here we report that under the appropriate conditions, batch cultures of strain AR36 with an initially high concentration of Glc can also be efficiently used to produce SA. However, the production profiles obtained suggest that fed-batch fermentations could also yield good results with this strain. The fact that the SA production ceases only when the glucose is exhausted suggests that higher titers could be achieved by adequate Glc feeding strategies and improvements in the acetate uptake capabilities of this strain, considering that it lacks PTS and could co-utilize Glc and acetate in these growing conditions. The fermentations reported here yielded an elevated ratio of SA to other byproducts of the pathway (10 times more SA than the main byproducts generated). Besides increasing the SA yield, this behavior is relevant for purifying SA from the culture broth. In fact, preliminary experiments concerning SA purification from the broth obtained from these cultures, resulted in an almost quantitative purification process (unpublished results). Furthermore, the highest yield of total aromatic compounds obtained represents 67% of the theoretical maximum, demonstrating the efficient redirection of carbon to the AAA pathway by strain AR36. Nevertheless, the relatively low cellular concentration present in the cultures, even when administering high concentrations of YE, represents a significant problem to this production system because it restricts the productivity of the process. Other strategies need to be utilized in order to increase the biomass concentration without increasing the supplemented YE, which will constitute an important improvement for scaling-up the process. Minimizing the metabolic load imposed by a high-copy plasmid while maintaining a sufficient gene dosage of the operon, should improve the distribution of resources that are directed towards biomass generation and SA production.
Construction of Escherichia coli derivatives and plasmids
The laboratory-evolved strain PB12, a derivative of PB11 (obtained by the inactivation of PTS in strain JM101 ), was the receptor of the genetic modifications described in this work [21, 22, 24]. The strains and plasmids used in this report are listed in Table 1. The chromosomal inactivations of aroK, aroL, and lacI genes were performed sequentially by homologous recombination of PCR products . In all cases, plasmid pKD3 was used as a PCR template in conjunction with tailored oligonucleotides containing 45 bp homology with the target chromosomal sequence (Additional file 2). Plasmid pKD46 expressed the Red recombinase system of bacteriophage lambda, and plasmid pCP20 allowed removal of the chloramphenicol resistance cassette after each event . Every step was verified by PCR, identifying the clones that presented the expected amplicon sizes when using different sets of oligonucleotides and chromosomal DNA as a template (Additional file 2). The pykF inactivated gene was transduced to PB12 aroK-aroL- using a P1 phage lysate obtained from strain PB28 (pykF::gen) . Transductants were selected on gentamycin plates (10 μg/ml), and the inactivation was confirmed by PCR.
The construction of the Aro6 operon and expression vector was accomplished in several steps (Figure 2 and Additional file 3). First, aroB, aroGfbr, tktA and aroE genes were amplified by PCR using Pfu DNA polymerase and ligated sequentially into the polylinker of plasmid pBRINT-Ts Cm . Chromosomal DNA from strain JM101 was used as a template for amplification of the required genes, with the exception of aroGfbr, which was amplified from plasmid pJLBaroGfbrtktA. Different sets of oligonucleotides were employed for the amplification of each gene (Additional file 2), which also generated flanking restriction sites and consensus Shine-Dalgarno sequences (AGGAGG) situated 8 bp upstream of the start of each coding sequence. The PCR products were inserted into the polylinker in the following order: aroB in the SmaI site, aroGfbr in the XhoI site, tktA in the EcoRV site, and aroE in the ApaI site. Simultaneously, plasmid pTrc327parlacI+ (Additional file 3) was built by ligating a PCR-amplified fragment containing the lacI gene, Trc promoter, polylinker, and transcriptional terminators of pTrc99A , into the ScaI and PvuI sites of pBR327par . The 4-gene operon present in plasmid pBRINT-Ts Cm was amplified by PCR with a unique set of oligonucleotides (Additional file 1) and ligated into pTrc327parlacI+ after digesting both with SacI and XbaI. Later, aroD was amplified by PCR (flanked by NheI sites) and ligated into compatible XbaI site of pTrc327parlacI+. Because of our interest in expressing the operon in a constitutive manner, a lacI- derivative of the initial pTrc327parlacI+ plasmid (without synthetic operon) was generated, and called pTrc327par (Additional file 3). The 5-gene operon was then transferred into SacI and NcoI sites of pTrc327par, giving rise to pTrcAro5. Finally, the zwf gene was inserted into the XbaI site, creating a 6-gene operon in pTrc327par. The resulting plasmid was named pTrcAro6 and transformed into AR2 and AR3 strains, generating AR26 and AR36, respectively (Table 1). The transformed strains were selected in LB plates supplemented with tetracycline (30 μg/ml).
Each step in the gene cloning and plasmid construction schemes was screened by endonuclease digestion and PCR, visualized with gel electrophoresis, and verified by DNA sequencing (3730, Perkin-Elmer/Applied Biosystems, USA). All the enzymes and reagents used in the molecular biology procedures were purchased from Fermentas (USA) and New England Biolabs (USA). When required, kits for the purification of PCR, plasmid, and agarose-embedded DNA were utilized (Roche, Switzerland). TOP10 cells (Invitrogen, USA) were used as a host for screening of DNA ligations during intermediate steps in vector construction.
Cultivation media and growth conditions
Composition of production medium
SA production medium (adjusted to pH 7.0 with 10 N NaOH) contained per liter: K2HPO4 (7.5 g), KH2PO4 (7.5 g), citric acid monohydrate (2.1 g), ammonium iron (III) citrate (0.3 g), concentrated H2SO4 (1.2 ml), MgSO4 (0.64 g), CaCl2 (0.06 g), (NH4)6(Mo7O24) (0.0037 g), ZnSO4 (0.0029 g), H3BO3 (0.0247 g), CuSO4 (0.0025 g), MnCl2 (0.0158 g), CoCl2 (0.00129 g), thiamine (0.001 g), and betaine (0.234 g) as an osmoprotectant. Tetracycline (30 μg/ml) was added to inocula and cultures whenever needed for plasmid maintenance. Glucose (filter-sterilized) and yeast extract (added before autoclaving) were supplied at the concentrations indicated for each experiment. The glucose was purchased from JT Baker (USA) and the autolysed yeast extract from BD Difco (USA).
Shake flask cultures
The inoculum preparation for the shake flask cultures started by the addition of 1 ml frozen aliquots to 250 ml shake flasks containing 25 ml of production medium supplemented with glucose (25 g/L) and yeast extract (15 g/L). The inoculum was grown at 37°C and 300 rpm until mid-exponential phase, and approximately 5% of the final volume was transferred to the test shake flasks and incubated under the same controlled conditions with media containing 15 g/L of glucose and 5 g/L of yeast extract. Cell growth was measured by monitoring the optical density at 600 nm (OD600) in a DU700 spectrophotometer (Beckman, USA), and samples were taken periodically, centrifuged, and the supernatant was stored at -20°C for metabolite analysis. These experiments were performed at least in triplicate. All cultures started at approximately 0.3 OD600.
Batch cultures were performed at least in duplicate using 1 L autoclavable glass bioreactors (Applikon, The Netherlands) with 500 ml of working volume. Bioreactors were connected to an Applikon ADI 1010 BioController and ADI 1025 controllers to monitor temperature, pH, impeller speed, and dissolved oxygen tension (DOT). The pH was kept at 7.0 by the addition of H3PO3 (3.3%) and NH4OH (10%). DOT in the culture medium was maintained by a continuous supply of filtered air (1 vvm), and by manually controlling the impeller speed (ranging from 500 to 1000 rpm) to ensure that DOT was kept above 20% at all times. The inoculum preparation for the fermentors started by the addition of 1 ml frozen aliquots to 500 ml shake flasks containing 50 ml of production medium supplemented with glucose (25 g/L) and yeast extract (15 g/L). The strains were grown at 37°C and 300 rpm until mid-exponential phase and approximately 5% of the final volume was transferred from each inoculum to previously prepared bioreactors containing the production medium. All fermentations were performed in presence of tetracycline (30 μg/ml). Cell growth was measured by monitoring optical density at 600 nm (OD600) in a spectrophotometer (DU700, Beckman, USA), and samples were taken periodically, centrifuged, and the supernatant was stored at -20°C for metabolite analysis. All the fermentations started at approximately 0.3 OD600.
The supernatant from each sample was properly diluted and filtered through 0.45 μM nylon membranes. Shikimic acid (SA), 3-dehydroshikimic acid (DHS), 3-dehydroquinic acid (DHQ), quinic acid (QA), gallic acid (GA), acetic acid, and glucose (Glc) concentrations were determined by HPLC using a Waters system (600E quaternary pump, 717 automatic injector, 2410 refraction index, and 996 photodiode array detectors; USA) equipped with an Aminex HPX-87H column (300 × 7.8 mm; 9 μm; Bio-Rad, USA). The mobile phase was 5 mM H2SO4, with a flow rate of 0.5 ml/min, maintained at 50°C. 3-deoxy-D-arabinoheptulosonate 7-phosphate (DAHP) concentrations were determined colorimetrically by the thiobarbituric acid assay . This method does not distinguish between DAHP and its unphosphorylated form, DAH, therefore in this work DAHP levels correspond to the sum of both compounds.
Data analysis and calculations
The measured concentrations of metabolites and biomass were normalized to the starting volume conditions to account for changes derived from pH control in fermentors. Data from independent experiments were averaged and presented in the corresponding graphs, where the error bars indicate the standard deviation for each point. Biomass concentration (X) was determined with a calibration curve between dry cellular weight and OD600, resulting in the equation X = 0.3587*OD600. Specific growth rate (μ) was determined by linearly fitting the biomass concentration to time during exponential phase with the following equation: lnX = lnXo + μ*t (where t is time, and X0 is the biomass concentration at initial time), displaying R2values >0.97. The yield of SA from Glc (YSA/Glc) was calculated with the average molar concentrations of SA and Glc, produced and consumed, respectively, at the point of highest SA concentration. The yield of total aromatic compounds from glucose (YTAC/Glc) was calculated with the combined molar yields of DAHP, DHQ, DHS, SA, QA, and GA at the point of highest SA concentration. The maximum theoretical yield of aromatic compounds was previously estimated as 0.86 molTAC/molGlc for a PTS- strain growing on glucose as only carbon source . The global volumetric SA productivity (Qpglobal) and the global volumetric Glc consumption rate (Qsglobal) were calculated taking into account the time needed to reach the maximum SA concentration. Besides the previous calculations, linearizations were made to obtain apparent biomass on substrate (YX/S) and product on biomass (YP/X) yields. Although these apparent yields do not take into account the yeast extract consumption, correlation values for linearizations in all experiments were found to be >0.95, allowing comparisons between them. These yields were used to calculate the specific productivity and specific consumption rate on the exponential phase (qpexp and qsexp, respectively) with the following equations: qpexp = YP/X*μ; qsexp = μ/YX/S.
The volumetric productivity and volumetric Glc consumption rate in stationary phase were determined by linearization of the first concentration data points at this stage versus time. The volumetric rates were utilized for calculation of specific production and consumption rates at stationary phase (qpsta and qssta, respectively) by dividing them by the average biomass concentration.
RNA extraction, cDNA synthesis and RT-qPCR analysis
Samples from batch fermentations of strain AR36 with 50 g/L of Glc and 15 g/L of YE were collected for RNA extraction at early exponential phase, EE (2 h ~ 1 OD600), late exponential phase, LE (8 h ~ 12 OD600), and stationary phase, ST (24 h ~ 17 OD600), to determine gene expression levels. For comparison of the data, samples from early exponential phase (3.5 h ~ 1 OD600) of strain AR3e (bearing an empty pTrc327par plasmid carrying only the tetracycline-resistance gene, Figure 2) cultured under the same conditions were also collected and processed. RNA was extracted using hot phenol equilibrated with water, and cDNA synthesis was performed using RevertAid H First Strand cDNA Synthesis kit (Fermentas, USA) and a mixture of specific DNA primers, as reported previously [21, 24]. qPCR experiments were performed with the ABI Prism 7300 Real Time PCR System (Applied Biosystems, USA) using Maxima SYBRGreen PCR Master Mix (Fermentas, USA) and reaction conditions previously described [21, 24]. The quantification technique used to compare data was the 2-∆∆CT method  and the results were normalized using the ihfB gene as an internal control. The same reproducible expression level for this gene was detected in all the strains and conditions analyzed . All qPCR experiments complied with the MIQE guidelines for publication of quantitative real-time PCR experiments . Using cells from two separate fermentations, RNA extraction and cDNA synthesis reactions were performed for each biological replicate at the indicated times and the gene expression values were measured by triplicate for each sample. Average values were graphed, with error bars representing standard deviation. Standard deviation was less than 30% in all cases.
The authors thank Ramón de Anda for technical support, José González for assistance with strain construction, Enrique Morett for helpful discussions, Paul Gaytán, Jorge Yáñez and Eugenio López for the synthesis of oligonucleotides and sequencing of DNA (Instituto de Biotecnología-UNAM, México), as well as Mercedes Enzaldo and Aurelia González for technical assistance. A. Rodriguez was supported by grants 215044 and 290604 from Consejo Nacional de Ciencia y Tecnología (CONACyT, México). This work was supported by grants 105782 and 177568 (CONACyT), 44126 and 126793 (FONSEC/SSA/ISSSTE/CONACyT), and IN224709 and IN205811 (DGAPA-PAPIIT-UNAM). The funding sources had no involvement in study design, collection and analysis of data, writing of the report or decision to publish.
- Farina V, Brown JD: Tamiflu: the supply problem. Angew Chem Int Ed Engl. 2006, 45: 7330-7334. 10.1002/anie.200602623View ArticleGoogle Scholar
- Kim HK, Park KJJ: A new efficient synthesis of oseltamivir phosphate (Tamiflu) from (-)-shikimic acid. Tetrahedron Lett. 2012, 53: 1561-1563. 10.1016/j.tetlet.2012.01.017.View ArticleGoogle Scholar
- Nie LD, Ding W, Shi XX, Quan N, Lu X: A novel and high-yielding asymmetric synthesis of oseltamivir phosphate (Tamiflu) starting from (-)-shikimic acid. Tetrahedron: Asymmetry. 2012, 23: 742-747. 10.1016/j.tetasy.2012.05.014.View ArticleGoogle Scholar
- Patnaik R, Spitzer RG, Liao JC: Pathway engineering for production of aromatics in Escherichia coli: confirmation of stoichiometric analysis by independent modulation of AroG, TktA, and Pps activities. Biotechnol Bioeng. 1995, 46: 361-370. 10.1002/bit.260460409View ArticleGoogle Scholar
- Gosset G, Yong-Xiao J, Berry A: A direct comparison of approaches for increasing carbon flow to aromatic biosynthesis in Escherichia coli. J Ind Microbiol. 1996, 17: 47-52. 10.1007/BF01570148View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Johansson L, Lindskog A, Silfversparre G, Cimander C, Nielsen KF, Lidén G: Shikimic acid production by a modified strain of E. coli (W3110.shik1) under phosphate-limited and carbon-limited conditions. Biotechnol Bioeng. 2005, 92: 541-552. 10.1002/bit.20546View ArticleGoogle Scholar
- Chen K, Dou J, Tang S, Yang Y, Wang H, Fang H, Zhou C: Deletion of the aroK gene is essential for high shikimic acid accumulation through the shikimate pathway in E. coli. Bioresour Technol. 2012, 119: 141-147.View ArticleGoogle Scholar
- Flores N, Xiao J, Berry A, Bolivar F, Valle F: Pathway engineering for the production of aromatic compounds in Escherichia coli. Nat Biotechnol. 1996, 14: 620-623. 10.1038/nbt0596-620View ArticleGoogle Scholar
- 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/bp025769pView ArticleGoogle Scholar
- Escalante A, Calderón R, Valdivia A, De Anda R, Hernández G, Ramírez 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-21View ArticleGoogle Scholar
- Juminaga D, Baidoo EEK, Redding-Johanson AM, Batth TS, Burd H, Mukhopadhyay A, Petzold CJ, Keasling JD: Modular engineering of L-tyrosine production in Escherichia coli. Appl Environ Microbiol. 2012, 78: 89-98. 10.1128/AEM.06017-11View ArticleGoogle Scholar
- Meza E, Becker J, Bolivar F, Gosset G, Wittmann C: Consequences of phosphoenolpyruvate:sugar phosphotranferase system and pyruvate kinase isozymes inactivation in central carbon metabolism flux distribution in Escherichia coli. Microb Cell Fact. 2012, 11: 127- 10.1186/1475-2859-11-127View ArticleGoogle Scholar
- Ahn JO, Lee HW, Saha R, Park MS, Jung JK, Lee DY: Exploring the effects of carbon sources on the metabolic capacity for shikimic acid production in Escherichia coli using in silico metabolic predictions. J Microbiol Biotechnol. 2008, 18: 1773-1784.Google Scholar
- Ghosh S, Chisti Y, Banerjee UC: Production of shikimic acid. Biotechnol Adv. 2012, 30: 1425-1431. 10.1016/j.biotechadv.2012.03.001View ArticleGoogle Scholar
- Glick BR: Metabolic load and heterologous gene expression. Biotechnol Adv. 1995, 13: 247-261. 10.1016/0734-9750(95)00004-AView ArticleGoogle Scholar
- Keasling JD: Gene-expression tools for the metabolic engineering of bacteria. Trends Biotechnol. 1999, 17: 452-460. 10.1016/S0167-7799(99)01376-1View ArticleGoogle Scholar
- Silva F, Queiroz JA, Domingues FC: Evaluating metabolic stress and plasmid stability in plasmid DNA production by Escherichia coli. Biotechnol Adv. 2012, 30: 691-708. 10.1016/j.biotechadv.2011.12.005View ArticleGoogle Scholar
- Carneiro S, Ferreira EC, Rocha I: Metabolic responses to recombinant bioprocesses in Escherichia coli. J Biotechnol. 2013, 164: 396-408. 10.1016/j.jbiotec.2012.08.026View ArticleGoogle Scholar
- Postma PW, Lengeler JW, Jacobson GR: Phosphoenolpyruvate: carbohydrate phosphotransferase systems. Escherichia coli and Salmonella: Cellular and Molecular Biology. Edited by: Neidhart FC. 1996, 1149-1174. USA: ASM Press, 2Google Scholar
- Flores N, Flores S, Escalante A, De Anda R, Leal L, Malpica R, Georgellis D, Gosset G, Bolivar F: Adaptation for fast growth on glucose by differential expression of central carbon metabolism and gal regulon genes in an Escherichia coli strain lacking the phosphoenolpyruvate:carbohydrate phosphotransferase system. Metab Eng. 2005, 7: 70-87. 10.1016/j.ymben.2004.10.002View ArticleGoogle Scholar
- Flores N, Leal L, Sigala JC, De Anda R, Escalante A, Martínez A, Ramírez OT, Gosset G, Bolivar F: Growth recovery on glucose under aerobic conditions of an Escherichia coli strain carrying a phosphoenolpyruvate: carbohydrate phosphotransferase system deletion by inactivating arcA and overexpressing the genes coding for glucokinase and galactose permease. J Mol Microbiol Biotechnol. 2007, 13: 105-116. 10.1159/000103602View ArticleGoogle Scholar
- Martínez K, De Anda R, Hernández G, Escalante A, Gosset G, Ramírez OT, Bolivar F: Coutilization of glucose and glycerol enhances the production of aromatic compounds in an Escherichia coli strain lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system. Microb Cell Fact. 2008, 7: 1- 10.1186/1475-2859-7-1View ArticleGoogle Scholar
- Aguilar C, Escalante A, Flores N, De Anda R, Riveros-McKay F, Gosset G, Morett E, Bolivar F: Genetic changes during a laboratory adaptive evolution process that allowed fast growth in glucose to an Escherichia coli strain lacking the major glucose transport system. BMC Genomics. 2012, 13: 385- 10.1186/1471-2164-13-385View ArticleGoogle Scholar
- Flores S, Gosset G, Flores N, De Graaf AA, Bolivar F: Analysis of carbon metabolism in Escherichia coli strains with an inactive phosphotransferase system by (13)C labeling and NMR spectroscopy. Metab Eng. 2002, 4: 124-137. 10.1006/mben.2001.0209View ArticleGoogle Scholar
- Cunningham DS, Liu Z, Domagalski N, Koepsel RR, Ataai MM, Domach MM: Pyruvate kinase-deficient Escherichia coli exhibits increased plasmid copy number and cyclic AMP levels. J Bacteriol. 2009, 191: 3041-3049. 10.1128/JB.01422-08View ArticleGoogle Scholar
- Ponce E: Effect of growth rate reduction and genetic modifications on acetate accumulation and biomass yields in Escherichia coli. J Biosci Bioeng. 1999, 87: 775-780. 10.1016/S1389-1723(99)80152-2View ArticleGoogle Scholar
- Zhu T, Phalakornkule C, Koepsel RR, Domach MM, Ataai MM: Cell growth and by-product formation in a pyruvate kinase mutant of E. coli. Biotechnol Prog. 2001, 17: 624-628. 10.1021/bp0100575View ArticleGoogle Scholar
- Kedar P, Colah R, Shimizu K: Proteomic investigation on the pyk-F gene knockout Escherichia coli for aromatic amino acid production. Enzyme Microb Technol. 2007, 41: 455-465. 10.1016/j.enzmictec.2007.03.018.View ArticleGoogle Scholar
- Al Zaid Siddiquee K, Arauzo-Bravo MJ, Shimizu K: Metabolic flux analysis of pykF gene knockout Escherichia coli based on 13 C-labeling experiments together with measurements of enzyme activities and intracellular metabolite concentrations. Appl Microbiol Biotechnol. 2004, 63: 407-417. 10.1007/s00253-003-1357-9View ArticleGoogle Scholar
- Amann E, Ochs B, Abel K: Tightly regulated tac promoter vectors useful for the expression of unfused and fused proteins in Escherichia coli. Gene. 1988, 69: 301-315. 10.1016/0378-1119(88)90440-4View ArticleGoogle Scholar
- Zurita M, Bolivar F, Soberon X: Construction and characterization of new cloning vehicles VII. Construction of plasmid pBR327par, a completely sequenced, stable derivative of pBR327 containing the par locus of pSC101. Gene. 1984, 28: 119-122. 10.1016/0378-1119(84)90094-5View ArticleGoogle Scholar
- Dell KA, Frost JW: Identification and removal of impediments to biocatalytic synthesis of aromatics from D-Glucose: rate-limiting enzymes in the common pathway of aromatic amino acid biosynthesis. J Am Chem Soc. 1993, 115: 11581-11589. 10.1021/ja00077a065.View ArticleGoogle Scholar
- Oldiges M, Kunze M, Degenring D, Sprenger GA, Takors R: Stimulation, monitoring, and analysis of pathway dynamics by metabolic profiling in the aromatic amino acid pathway. Biotechnol Prog. 2004, 20: 1623-1633. 10.1021/bp0498746View ArticleGoogle Scholar
- Lütke-Eversloh T, Stephanopoulos G: Combinatorial pathway analysis for improved L-tyrosine production in Escherichia coli: identification of enzymatic bottlenecks by systematic gene overexpression. Metab Eng. 2008, 10: 69-77. 10.1016/j.ymben.2007.12.001View ArticleGoogle Scholar
- Flores S, De Anda-Herrera R, Gosset G, Bolívar F: Growth-rate recovery of Escherichia coli cultures carrying a multicopy plasmid, by engineering of the pentose-phosphate pathway. Biotechnol Bioeng. 2004, 87: 485-494. 10.1002/bit.20137View ArticleGoogle Scholar
- Sprenger GA: Genetics of pentose-phosphate pathway enzymes of Escherichia coli K-12. Arch Microbiol. 1995, 164: 324-330. 10.1007/BF02529978View ArticleGoogle Scholar
- Knop DR, Draths KM, Chandran SS, Barker JL, Von Daeniken R, Weber W, Frost JW: Hydroaromatic equilibration during biosynthesis of shikimic acid. J Am Chem Soc. 2001, 123: 10173-10182. 10.1021/ja0109444View ArticleGoogle Scholar
- 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.001View ArticleGoogle Scholar
- Gimenez R, Nuñez MF, Badia J, Aguilar J, Baldoma L: The gene yjcG, cotranscribed with the gene acs, encodes an acetate permease in Escherichia coli. J Bacteriol. 2003, 185: 6448-6455. 10.1128/JB.185.21.6448-6455.2003View ArticleGoogle Scholar
- Messing J: Recombinant DNA technical bulletin. NIH Publ. 1979, 2: 43-48.Google Scholar
- Datsenko K, Wanner B: 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.120163297View ArticleGoogle Scholar
- Cherepanov PP, Wackernagel W: Gene disruption in Escherichia coli: TcR and KmR cassettes with the option of Flp-catalyzed excision of the antibiotic-resistance determinant. Gene. 1995, 158: 9-14. 10.1016/0378-1119(95)00193-AView ArticleGoogle Scholar
- Ponce E, Flores N, Martinez A, Valle F, Bolivar F: Cloning of the two pyruvate kinase isoenzyme structural genes from Escherichia coli: the relative roles of these enzymes in pyruvate biosynthesis. J Bacteriol. 1995, 177: 5719-5722.Google Scholar
- Le Borgne S, Palmeros B, Valle F, Bolivar F, Gosset G: pBRINT-Ts: a plasmid family with a temperature-sensitive replicon, designed for chromosomal integration into the lacZ gene of Escherichia coli. Gene. 1998, 223: 213-219. 10.1016/S0378-1119(98)00168-1View ArticleGoogle Scholar
- Balderas-Hernández VE, Sabido-Ramos A, Silva P, Cabrera-Valladares N, Hernández-Chávez G, Baez-Viveros JL, Martínez A, Bolivar F, Gosset G: Metabolic engineering for improving anthranilate synthesis from glucose in Escherichia coli. Microb Cell Fact. 2009, 8: 19- 10.1186/1475-2859-8-19View ArticleGoogle Scholar
- Weissbach A, Hurwitz J: The Formation of 2-Keto-3-deoxyheptonic Acid in Extracts of Escherichia coli B. J Biol Chem. 1958, 234: 705-709.Google Scholar
- 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.1262View ArticleGoogle Scholar
- Bustin SA, Benes V, Garson JA, Hellemans J, Hugget J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.