Skip to main content

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

Abstract

Background

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).

Results

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.

Conclusions

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.

Background

Shikimic acid (SA) is an intermediate compound in the aromatic amino acid (AAA) biosynthetic pathway in plants and bacteria (Figure 1). This metabolite is utilized as starting material in the chemical synthesis of oseltamivir phosphate (Tamiflu), used for influenza treatment [13]. Several genetic strategies have been reported for improving SA productivity and yield in Escherichia coli. These strategies aim to increase the availability of the direct precursors of the AAA pathway, erythrose 4-phosphate (E4P) and phosphoenolpyruvate (PEP), by genetic alterations that promote a convenient redistribution of carbon fluxes in the central metabolism [4, 5]. Complementary approaches include the interruption of the AAA pathway after SA formation by inactivation of the genes coding for shikimate kinases (aroK and aroL), as well as enhancements in carbon channeling towards SA by overexpression of feedback-resistant DAHP synthases, shikimate dehydrogenase, transketolase, and DHQ synthase enzymes (coded by aroFGHfbr, aroE, tktA, and aroB, respectively) (Figure 1) [69]. In an attempt to further increase the intracellular availability of PEP, strains overexpressing PEP synthase (coded by ppsA), or lacking the PEP:carbohydrate phosphotransferase system (PTS) and the pyruvate kinase isozymes (coded by pykF and pykA), have also been evaluated (Figure 1) [1013].

Figure 1
figure 1

Genetic modifications applied in this report to enhance the production of shikimic acid from glucose in the laboratory-evolved E. coli strain PB12. Inactivated genes are indicated with a cross and plasmid-expressed genes are circled (see Methods for details). Dashed arrows indicate more than one catalytic step. G6P = glucose 6-phosphate; F6P = fructose 6-phosphate; GAP = glyceraldehyde 3-phosphate; 6PGNL = 6-phosphogluconolactone; Ru5P = ribulose 5-phosphate; R5P = ribose 5-phosphate; Xu5P = xylulose 5-phosphate; S7P = sedoheptulose 7-phosphate; E4P = erythrose 4-phosphate; PEP = phosphoenolpyruvate; PYR = pyruvate; ACoA = acetyl-coenzyme A; Ace-P = acetyl phosphate; CIT = citrate; OAA = oxaloacetate; DAHP = 3-deoxy-D-arabinoheptulosonate 7-phosphate; DHQ = 3-dehydroquinic acid; DHS = 3-dehydroshikimic acid; QA = quinic acid; GA = gallic acid; SA = shikimic acid; S3P = shikimate 3-phosphate; CHO = chorismate; IICBGlc = membrane component of glucose-specific PTS permease; E1 = PTS enzyme 1; Hpr = PTS histidine protein; IIAGlc = cytosolic component of glucose-specific PTS permease. Genes coding for enzymes not named in the figure: galP, galactose permease; glk, glucokinase; pgi, phosphoglucose isomerase; pfkA, 6-phosphofructokinase I; fbaA, fructose bisphosphate aldolase class II; gapA, glyceraldehyde 3-phosphate dehydrogenase; eno, enolase; actP, acetate permease; acs, acetyl-coenzyme A synthetase; pta, phosphate acetyltransferase; ackA, acetate kinase; poxB, pyruvate oxidase.

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 [1619]. 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 [20]. 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 [24].

It was reported that PB12, which assimilates glucose by the non-PTS symporter GalP [25], 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 [11]. 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 [26], low acetate production due to less glycolytic overflux [13, 27, 28], or higher concentrations of the AAA pathway enzymes [29]. 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].

Here, we propose that a PTS-pykF- strain has the potential to increase the yield and titer of SA when compared to an isogenic pykF+ strain, provided that the gene expression system permits an appropriate temporal coordination in the synthesis of the enzymes required to channel the carbon towards SA, while reducing the accumulation of acetate and intermediate compounds in the AAA pathway. In order to accomplish this goal a synthetic operon was constructed containing the coding sequences of six genes selected from the pentose phosphate (PPP) and AAA pathways (Figure 1), controlled by a single constitutive Trc promoter [31] and inserted it into a high-copy plasmid containing a region that confers segregational stability (Figure 2) [32]. The resulting plasmid was transformed into a modified PB12 strain with inactive aroK, aroL, pykF, and lacI genes, and was cultured in fermentors using mineral media containing Glc and supplemented with YE. Overall, the strategy proposed in this report allowed the overproduction of SA from the beginning of the culture, resulting in a high titer and yield of SA with relatively low accumulation of acetate and aromatic byproducts. It was also found that, under the high-substrate conditions tested, the SA titer was independent of the YE concentration and the maximum biomass produced depended exclusively on the initial YE concentration but not on the amount of glucose.

Figure 2
figure 2

Simplified scheme of the steps required in the construction of plasmid pTrcAro6, carrying 6 synthetic genes under the control of the Trc promoter. First, a segment of pTrc99A was amplified and ligated into pBR327par, creating plasmid pTrc327par. A synthetic operon comprising the required genes (aroB, tktA, aroGfbr, aroE, aroD, and zwf) was assembled separately and transferred to pTrc327par, generating pTrcAro6. The dotted lines indicate the site and orientation of some of the performed ligation reactions. Only the relevant restriction sites are displayed. A more detailed scheme of the constructions is presented in Additional file 3.

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 [32], after incorporating a fragment containing the promoter, polylinker, and transcriptional terminators from pTrc99A (Figure 2).

The initial part of the operon was constructed by sequential amplification and ligation of the first 4 coding sequences (aroB, tktA, aroGfbr, and aroE) into the polylinker of plasmid pBRINT-Ts Cm, used as a cloning scaffold (see Methods). Later, the 4-gene construction was transferred to the hybrid plasmid pTrc327par in conjunction with 2 more genes (aroD and zwf), leading to an 8Kb operon contained in a 12Kb plasmid (Figure 2). The resulting plasmid, termed pTrcAro6, was transformed into the PB12 aroK-aroL- strain devoid of the lacIq gene, allowing constitutive expression of the genes of interest (Table 1). For simplicity, the generated PB12 aroK-aroL-lacI- strain was termed AR2. After the pykF gene was inactivated in AR2, the resulting strain was named AR3. Strains derived from AR2 and AR3 carrying plasmid pTrcAro6 were named AR26 and AR36, respectively (Table 1).

Table 1 Escherichia coli strains and plasmids utilized in this report

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 [3335]. 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 [36]; 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) [25]. 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 [37]. 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.

Even though SA accumulated in all cases, as expected for mutants in aroK and aroL, the strains containing pTrcAro6 reached higher SA concentrations than the ones with an empty plasmid (Figure 3b). Moreover, the SA titer was almost two times higher in AR36 than in AR26 (6.1 g/L vs. 3.3 g/L). A decrease in Glc consumption was observed in strain AR26 after approximately 18 h of culture, correlating with high acetate concentration and an arrest in the production of SA. In contrast, strain AR36 exhibited constant Glc consumption and negligible amounts of acetate were produced (Figure 3c, 3d). These results demonstrate that the genes present in the artificial operon are functional and promote the production of SA since the beginning of the culture. Their constitutive expression diminished the specific growth rate (μ) by 25% in the pykF+ background, and marginally increased it in the pykF- variant, but did not cause significant changes to the maximum biomass produced (Xmax) compared to strains with an empty plasmid (Figure 3a). Remarkably, in the operon-expressing strains under these growth conditions, the inactivation of the pykF gene increased the production of SA, eliminated the accumulation of acetate, and allowed steady Glc consumption.

Figure 3
figure 3

Behavior of strains AR26, AR36, and their empty-plasmid derivatives, AR2e ( pykF+) and AR3e ( pykF-), using shake flasks containing 15 g/L of Glc and 5 g/L of YE (a,b,c,d), and 1 L fermentors containing 100 g/L of Glc and 15 g/L of YE (e). a) Growth; b) SA production; c) Glc consumption; d) acetate production; e) Glc consumption and SA production of AR26 and AR36 in fermentors. Error bars represent standard deviation.

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.

Under these conditions strain AR36 produced 42 g/L of SA in 60 h, consuming all the Glc, and accumulating 12 g/L of acetate. In contrast, after 47 h strain AR26 produced a maximum of 13 g/L of SA, did not exhaust the Glc, and accumulated 29 g/L of acetate (Figure 3e and Table 2). Regardless of the controlled conditions in the fermentors, where the pH was kept at 7 and the DOT was higher than 20% at all times, the production profiles of both strains resembled the behavior observed in shake flasks, with AR26 producing more acetate and less SA. Even when the global volumetric Glc consumption rate (Qsglobal), μ and Xmax attained by both strains were similar, the productivity, yield, and titer were more than twofold higher in AR36 than in AR26 (Figure 3e and Table 2).

Table 2 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

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.

First, 50 g/L of Glc and 15 g/L of YE were utilized (Figure 4a). Growth occurred during the first 10 h, generating 6.3 g/L of dry cell weight with a μ of 0.53 h-1. Under this condition, 24 g/L of SA were produced in 32 h. Glc consumption and SA production occurred since the beginning of the fermentation and lasted until Glc exhaustion, although the specific Glc consumption rate (qs) and specific SA productivity (qp) were higher in exponential phase (Table 3). The resulting yield of SA on Glc (YSA/Glc) was 0.47 mol/mol and the global volumetric SA productivity (Qpglobal) was 0.74 gSA/L*h (Table 3). With respect to the accumulation of byproducts in the SA pathway, concentrations of 2.4 g/L of DAHP, 2.1 g/L of DHS, 1.4 g/L of QA, 0.4 g/L of GA, and 0.3 g/L of DHQ, were present in the supernatant at the end of the fermentation (Figure 5a). Under these conditions, virtually no acetate was produced during the course of the fermentation, reaching a maximum concentration of 1.5 g/L after 32 h (Figure 4a).

Figure 4
figure 4

Fermentation profile of strain AR36 cultivated in 1 L bioreactors with three different substrate concentrations. a) 50 g/L of Glc and 15 g/L of YE; b) 100 g/L of Glc and 15 g/L of YE; c) 100 g/L of Glc and 30 g/L of YE. Glc: circles; SA: squares; acetate: open triangles; biomass concentration: inverted triangles. Error bars represent standard deviation.

Table 3 Comparison of measured metabolites and calculated kinetic and stoichiometric parameters between three fermentations of strain AR36 with different substrate concentrations
Figure 5
figure 5

Aromatic byproducts of the SA pathway detected in 1 L fermentor cultures of strain AR36 using three different substrate concentrations. a) 50 g/L of Glc and 15 g/L of YE; b) 100 g/L of Glc and 15 g/L of YE; c) 100 g/L of Glc and 30 g/L of YE. Diamonds: DAHP (3-deoxy-D-arabinoheptulosonate 7-phosphate); squares: DHQ (3-dehydroquinic acid); circles: DHS (3-dehydroshikimic acid); triangles: QA (quinic acid); inverted triangles: GA (gallic acid). Error bars represent standard deviation.

Considering that 50 g/L of Glc were consumed completely, a second batch experiment was initiated with 100 g/L of Glc and 15 g/L of YE. As stated in the comparison with AR26 in the previous section, AR36 grown under these conditions produced approximately 42 g/L of SA in 60 h (Figure 4b). In this case, after consuming about 100 g/L of glucose and attaining the maximum concentration of SA, the strain produced 12 g/L of acetate. The values obtained for YSA/Glc, Qpglobal, Qsglobal, Xmax, and μ, were similar to those obtained with 50 g/L of Glc and 15 g/L of YE (Table 3). These experiments show that when using the same YE concentration, twice the amount of Glc is consumed in almost twice the time, indicating that the average glucose consumption rate is maintained between both culture conditions. Concentrations of 4.8 g/L of DAHP, 2.8 g/L of DHS, 3.4 g/L of QA, 0.7 g/L of GA, and 0.9 g/L of DHQ, were present in the supernatant after 60 h (Figure 5b). Interestingly, when doubling the Glc concentration the intermediate products of the AAA pathway increased in a fairly proportional manner with the SA, indicating that the consumption of 100 g/L of Glc did not apparently generate new carbon flux bottlenecks. As a result, the amount of SA formed with respect to the total aromatic compounds produced was close to 80% in both experiments (Figure 6).

Figure 6
figure 6

Molar percentage of each aromatic compound produced in strain AR36 with respect to the total in batch cultures starting with: a) 50 g/L of Glc and 15 g/L of YE; b) 100 g/L of Glc and 15 g/L of YE; c) 100 g/L of Glc and 30 g/L of YE. Calculated yields and molar ratios of the produced aromatic compounds are shown below each bar. The comparisons were made with the concentrations measured in the supernatant at the end of the fermentations. YSA/Glc = yield of SA from Glc; YTAC/Glc = yield of total aromatic compounds from Glc; Ymax = maximum theoretical yield of aromatic compounds.

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 [39].

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.

The results indicate that the presence and expression of the operon in strain AR36 increases the transcriptional levels of several genes coding for glycolytic enzymes during the EE and LE phases (Figure 1 and Figure 7a). The rise in expression of genes galP and glk is particularly interesting because it has been reported that their products control the import and phosphorylation of glucose in PB12, the parental strain of AR36 [21, 25]. Furthermore, there is a significant increase in the transcriptional levels of pgi and eno, but not pykA. These changes may translate into higher availability of PEP and fructose 6-P (which can be directly converted into E4P by plasmid-encoded transketolase), increasing the yield of aromatic compounds. We theorize that the observed upregulation of glycolytic genes in strain AR36 could be one of the consequences to low levels of some glycolytic intermediates (glucose 6-phosphate, fructose 6-phosphate and PEP), caused by the strong and constitutive expression of the operon-encoded enzymes that consume these metabolites.

Figure 7
figure 7

Transcriptional changes resulting from the expression of Aro6 synthetic operon in strain AR36. For comparison, the transcription level of each gene was determined at three different points in the growth curve of strain AR36, grown with 50 g/L of Glc and 15 g/L of YE (see Figure 4a). All data were normalized against the values obtained from strain AR3e at early exponential growth phase. a) Genes coding for glycolytic enzymes; b) genes involved in acetate assimilation and biosynthesis; c) genes comprised in the synthetic Aro6 operon (see Figure 1 and Figure 2). Black bars: early exponential phase; gray bars: late exponential phase; white bars: stationary phase. Error bars represent standard deviation.

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 [21]. 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.

Conclusions

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 [10]. 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.

Methods

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 [41]), 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 [42]. 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 [43]. 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) [44]. 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 [45]. 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[46]. 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 [31], into the ScaI and PvuI sites of pBR327par [32]. 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.

Fermentor cultures

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.

Metabolite quantification

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 [47]. 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 [10]. 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 [48] 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 [24]. All qPCR experiments complied with the MIQE guidelines for publication of quantitative real-time PCR experiments [49]. 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.

References

  1. Farina V, Brown JD: Tamiflu: the supply problem. Angew Chem Int Ed Engl. 2006, 45: 7330-7334. 10.1002/anie.200602623

    Article  CAS  Google Scholar 

  2. 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.

    Article  CAS  Google Scholar 

  3. 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.

    Article  CAS  Google Scholar 

  4. 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.260460409

    Article  CAS  Google Scholar 

  5. 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/BF01570148

    Article  CAS  Google Scholar 

  6. 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.

    Article  CAS  Google Scholar 

  7. 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.20546

    Article  CAS  Google Scholar 

  8. 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.

    Article  CAS  Google Scholar 

  9. 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-620

    Article  CAS  Google Scholar 

  10. 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

    Article  CAS  Google Scholar 

  11. 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-21

    Article  Google Scholar 

  12. 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-11

    Article  CAS  Google Scholar 

  13. 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-127

    Article  CAS  Google Scholar 

  14. 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.

    CAS  Google Scholar 

  15. Ghosh S, Chisti Y, Banerjee UC: Production of shikimic acid. Biotechnol Adv. 2012, 30: 1425-1431. 10.1016/j.biotechadv.2012.03.001

    Article  CAS  Google Scholar 

  16. Glick BR: Metabolic load and heterologous gene expression. Biotechnol Adv. 1995, 13: 247-261. 10.1016/0734-9750(95)00004-A

    Article  CAS  Google Scholar 

  17. Keasling JD: Gene-expression tools for the metabolic engineering of bacteria. Trends Biotechnol. 1999, 17: 452-460. 10.1016/S0167-7799(99)01376-1

    Article  CAS  Google Scholar 

  18. 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.005

    Article  CAS  Google Scholar 

  19. 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.026

    Article  CAS  Google Scholar 

  20. 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, 2

    Google Scholar 

  21. 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.002

    Article  CAS  Google Scholar 

  22. 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/000103602

    Article  CAS  Google Scholar 

  23. 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-1

    Article  Google Scholar 

  24. 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-385

    Article  CAS  Google Scholar 

  25. 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.0209

    Article  CAS  Google Scholar 

  26. 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-08

    Article  CAS  Google Scholar 

  27. 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-2

    Article  CAS  Google Scholar 

  28. 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/bp0100575

    Article  CAS  Google Scholar 

  29. 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.

    Article  CAS  Google Scholar 

  30. 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-9

    Article  CAS  Google Scholar 

  31. 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-4

    Article  CAS  Google Scholar 

  32. 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-5

    Article  CAS  Google Scholar 

  33. 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.

    Article  CAS  Google Scholar 

  34. 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/bp0498746

    Article  CAS  Google Scholar 

  35. 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.001

    Article  Google Scholar 

  36. 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.20137

    Article  CAS  Google Scholar 

  37. Sprenger GA: Genetics of pentose-phosphate pathway enzymes of Escherichia coli K-12. Arch Microbiol. 1995, 164: 324-330. 10.1007/BF02529978

    Article  CAS  Google Scholar 

  38. 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/ja0109444

    Article  CAS  Google Scholar 

  39. 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

    Article  Google Scholar 

  40. 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.2003

    Article  CAS  Google Scholar 

  41. Messing J: Recombinant DNA technical bulletin. NIH Publ. 1979, 2: 43-48.

    CAS  Google Scholar 

  42. 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.120163297

    Article  CAS  Google Scholar 

  43. 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-A

    Article  CAS  Google Scholar 

  44. 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.

    CAS  Google Scholar 

  45. 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-1

    Article  CAS  Google Scholar 

  46. 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-19

    Article  Google Scholar 

  47. 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 

  48. 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

    Article  CAS  Google Scholar 

  49. 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.112797

    Article  CAS  Google Scholar 

Download references

Acknowledgements

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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Guillermo Gosset.

Additional information

Competing interests

The authors have declared no competing interests.

Authors’ contributions

AR, JLB and FB conceived the study and designed the experiments. AR and JAM carried out all fermentation experiments and data analysis. JAM performed the calculation of fermentation parameters. NF carried out RNA extraction and RT-qPCR experiments. GH was involved in metabolite quantification and analysis. AR performed the molecular biology and strain construction procedures and wrote the manuscript. JAM, OTR, GG and FB critically revised the results and the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

12934_2013_881_MOESM1_ESM.pdf

Additional file 1:Nucleotide sequence of the synthetic operon constructed in this work and present in plasmid pTrcAro6. The aroGfbr gene included in this construction was a gift from DuPont™-Genencor®, therefore its coding sequence cannot be disclosed. Each nucleotide of aroGfbr is indicated with an “n” except for the ones corresponding to its start and stop codons. (PDF 14 KB)

Additional file 2:Oligonucleotides utilized in this work.(PDF 53 KB)

Additional file 3:Detailed scheme of the construction of plasmid pTrcAro6.(PDF 896 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Rodriguez, A., Martínez, J.A., Báez-Viveros, J.L. et al. 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 12, 86 (2013). https://doi.org/10.1186/1475-2859-12-86

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2859-12-86

Keywords