- Open Access
Rewiring glycerol metabolism for enhanced production of poly-γ-glutamic acid in Bacillus licheniformis
- Yangyang Zhan†1,
- Bojie Sheng†2,
- Huan Wang1,
- Jiao Shi1,
- Dongbo Cai1,
- Li Yi1,
- Shihui Yang1,
- Zhiyou Wen3, 4,
- Xin Ma1 and
- Shouwen Chen1, 2Email authorView ORCID ID profile
© The Author(s) 2018
- Received: 22 January 2018
- Accepted: 1 November 2018
- Published: 9 November 2018
Poly-γ-glutamic acid (γ-PGA) is a natural polymer with great potential applications in areas of agriculture, industry, and pharmaceutical. The biodiesel-derived glycerol can be used as an attractive feedstock for γ-PGA production due to its availability and low price; however, insufficient production of γ-PGA from glycerol is limitation.
The metabolic pathway of Bacillus licheniformis WX-02 was rewired to improve the efficiency of glycerol assimilation and the supply of NADPH for γ-PGA synthesis. GlpK, GlpX, Zwf, and Tkt1 were found to be the key enzymes for γ-PGA synthesis using glycerol as a feedstock. Through combinational expression of these key enzymes, the γ-PGA titer increased to 19.20 ± 1.57 g/L, which was 1.50-fold of that of the wild-type strain. Then, we studied the flux distributions, gene expression, and intracellular metabolites in WX-02 and the recombinant strain BC4 (over-expression of the above quadruple enzymes). Our results indicated that over-expression of the quadruple enzymes redistributed metabolic flux to γ-PGA synthesis. Furthermore, using crude glycerol as carbon source, the BC4 strain showed a high productivity of 0.38 g/L/h, and produced 18.41 g/L γ-PGA, with a high yield of 0.46 g γ-PGA/g glycerol.
The approach to rewiring of metabolic pathways enables B. licheniformis to efficiently synthesize γ-PGA from glycerol. The γ-PGA productivity reported in this work is the highest obtained in glutamate-free medium. The present study demonstrates that the recombinant B. licheniformis strain shows significant potential to produce valuable compounds from crude glycerol.
- Bacillus licheniformis
- Poly-γ-glutamic acid
- Glycerol metabolism
- Gluconeogenesis pathway
- Combinatorial optimization
Poly-γ-glutamic acid (γ-PGA) is a group of natural polymers consisting of l- and/or d-glutamic acid monomers polymerized through γ-glutamyl bonds . Due to its features of cation chelating, hygroscopicity, water-solubility, biodegradability, and non-toxic towards humans and the environment, γ-PGA has been widely used in foods, medicine, cosmetics, and agriculture industries . The bacteria Bacillus has been used as an efficient γ-PGA producer through various metabolically engineering manipulations [3–5]. In those studies, however, glucose and glutamate are commonly used as substrates, resulting in high production cost . It is therefore necessary to use cost-effective substrates for the production of γ-PGA.
The biodiesel-derived glycerol is an ideal substrate for γ-PGA production due to its abundance and low price. About 10 kg of crude glycerol will be generated from every 100 kg of biodiesel production . The rapid development of biodiesel industry has produced a large amount of crude glycerol, which makes this chemical as under-valued material . There are some chemical applications utilizing pure glycerol as feedstock, but it is not economical to refine the crude glycerol into pure glycerol . Crude glycerol is biologically converted into the value-added products such as acetoin , succinate , n-butanol , 3-hydroxypropionic acid , 1,3-propanediol , and poly-3-hydroxybutyrate . The bioconversion of glycerol will provide an opportunity to recycle the waste disposal and produce valuable chemicals. Moreover, glycerol can generate more reducing agents than other carbon sources, and thus is able to produce high amounts of biochemical products .
NADPH is proven to be an important cofactor for γ-PGA biosynthesis  and is mainly generated from PPP in B. licheniformis. When glycerol is used as the substrate, the PPP may be weakened in several ways to cause a shortage of NADPH . For instance, the activities of isocitrate dehydrogenase (ICDH, EC 184.108.40.206), glucose-6-phosphate dehydrogenase (Zwf, EC 220.127.116.11), and 6-phosphogluconate dehydrogenase (Gnd, EC 18.104.22.168) are decreased compared to the glucose as the substrate . The shortage of NADPH may be another reason for the insufficient production of γ-PGA from glycerol. Moreover, NADPH can prevent from cell damage resulting from reactive oxygen species [20, 23]. When cells are cultured in crude glycerol, the insufficient NADPH generation may aggravate the inhibitory effects of impurities of crude glycerol such as methanol, salts, fatty acids, and heavy metals on cell proliferation and metabolism . Collectively, improvement of NADPH supply is required for γ-PGA synthesis regarding crude glycerol utilization as feedstock.
Bacillus licheniformis WX-02 strain has been proven as an efficient γ-PGA producer . However, insufficient production of γ-PGA from glycerol was still unsolved using glycerol as a substrate . The aim of this work is to improve γ-PGA production by genetically modifying the glycerol assimilation pathway, gluconeogenesis steps, and the availability of NADPH in B. licheniformis when using glycerol as substrate. To gain an insight into mechanism for the γ-PGA overproduction, the gene expression, intracellular metabolites, and flux distributions in the wild-type strain and the engineered strain were studied.
Improve glycerol kinase expression to increase glycerol assimilation efficiency
Previous research demonstrated that the activity of glycerol kinase was the limiting step for glycerol assimilation in many microorganisms [18–20]. Because of the insufficient coordination of the glycerol kinase, the low glycerol assimilation efficiency limits cell growth and metabolic production [18–20]. To improve the expression of glycerol kinase, the expression levels of glpK, glpP, and tglpFK were modulated. As a terminator-dependent carbon catabolite repression (CCR) mechanism is operative for the glpFK operon by the terminator tglpFK, whose formation is prevented by the Gly3P-activated antiterminator GlpP , the terminator tglpFK was knocked out by homologous recombination, and the glpK and glpP were overexpressed, respectively, by introducing an additional copy of each gene under P43 promoter.
Identification of the limiting enzyme in the gluconeogenesis pathway
Insufficient G6P generation in the absence of glucose leads to scant flux towards PPP. To identify limiting enzymes in upper gluconeogenesis pathway for γ-PGA synthesis, each of genes (pgi, glpX and fba) was overexpressed individually, and the resultant strains were denoted as WX02-pgi, WX02-glpX, and WX02-fba, respectively. As shown in Fig. 2b, over-expression of glpX significantly increased γ-PGA titer (27.67%) and glycerol consumption (26.81%) compared to the wild-type strain. However, the glycerol utilization efficiency and γ-PGA titer were decreased obviously in WX02-fba strain. Strengthening of fba could increase the amount of FBP, which further stimulated the kinase function of HprK/P and triggered the CcpA-dependent CCR mechanism . As a result, the glycerol metabolism was reduced in fba overexpressed strain. Additionally, over-expression of pgi had no significant effect on γ-PGA production. Collectively, this experiment identified GlpX, the second enzyme in the upper gluconeogenesis pathway, as a bottleneck in glycerol metabolism, and this gene was considered as the metabolic engineering target for γ-PGA overproduction.
Strengthening pentose phosphate pathway flux to enhance γ-PGA production
Strengthening PPP flux provides more NADPH supply, which may be beneficial for product biosynthesis . In this study, the effects of over-expression of these five genes (zwf, gnd, rpe, tkt1, and tkt2) (Fig. 1) on γ-PGA synthesis were investigated. As shown in Fig. 2c, the highest γ-PGA titer of 15.03 g/L was obtained by overexpressing zwf, which was 17.14% enhancement over the wild type. Also, the glycerol consumption of WX02-zwf was increased by 14.19%. Compared to the wild type, the γ-PGA titer was increased by 14.42% in WX02-tkt1 strain; while over-expression of gnd, tkt2, and rpe presented almost no influence on the γ-PGA synthesis and glycerol consumption (Fig. 2c). Thus, the Zwf and Tkt1 were chosen to direct our subsequent metabolic engineering efforts to boost γ-PGA production.
Combinatorial over-expression of glpK, glpX, zwf, and tkt1 genes for enhancement of γ-PGA synthesis
Comparison of γ-PGA fermentation between WX-02 and mutant strains
B. licheniformis strains
Glycerol consumption (g/L)
γ-PGA titer (g/L)
10.99 ± 0.34
30.51 ± 0.98
12.83 ± 0.27
12.39 ± 0.65
37.17 ± 0.85
17.42 ± 1.17
13.35 ± 0.32
38.95 ± 0.92
18.85 ± 0.53
13.33 ± 0.69
40.19 ± 0.73
19.20 ± 1.57
Combinatorial over-expression of glpK, glpX, zwf, and tkt1 triggers metabolic flux redistribution
Cofactors such as NADH and NADPH play important roles in regulating redox relevant reactions and generating energy. NADH is involved in ATP generation of oxidative phosphorylation and NADH-dependent products synthesis, while NADPH is essential for γ-PGA biosynthesis [26, 27]. The intracellular NADPH and NADH pools of BC4 and WX-02 strains were also compared (Additional file 1). The strain BC4 produced 20.05 μmol/gDCW NADPH, a 30.45% increase than that of wild-type strain (15.37 μmol/gDCW), while the NADH content was decreased by 20%. The ratio of NADPH/NADH was 2.05, increased by 62.70%. These results suggested that combinatorial over-expression of quadruple genes might affect the overall metabolic network, and further led to the metabolic flux redistribution.
Transcriptional levels and intracellular metabolites analysis
To gain a further insight into mechanisms for the γ-PGA overproduction, the intracellular metabolite concentrations in WX-02 and BC4 were measured by GC–MS-based metabolome approach [28, 29]. Based on our results, the metabolites (containing intermediates of glycolysis, PPP, TCA cycle, and amino acids pools) with significant differences among the two strains were listed in Fig. 5b. In terms of the glycerol pathway, the pool sizes of Gly3P in engineered strains were higher than that of the WX-02. The metabolites involved in the PPP, such as Ru5P, uracil, and tyrosine, were increased in engineered strains (increased by 40.60–328.66%). TCA cycle metabolites, such as succinate, fumarate, and citrate, showed a higher abundance in the mutants compared to the wild-type strain WX-02. The glycolysis metabolite 3-phosphoglyceric acid (3PG) was increased by 9.35-fold. These results were consistent with the above-performed transcription analysis, explaining γ-PGA production enhancement. Due to the difficulty in analyzing glycolytic intermediates, some metabolites have not been detected like G6P, F6P, and FBP in cell extracts.
γ-PGA production by recombinant BC4 strain with crude glycerol as substrate
Comparison of γ-PGA production by Bacillus strains using glycerol as substrate
Main nutrients (g/L)
Culture time (h)
γ-PGA titer (g/L)
γ-PGA productivity (g/L/h)
B. licheniformis ATCC 9945a
l-Glutamic acid, glycerol, citric acid, NH4Cl
B. licheniformis CCRC 12826
B. licheniformis WBL-3
B. subtilis RKY3
l-Glutamic acid, glycerol, citric acid, NH4Cl, yeast extract
l-Glutamic acid, glycerol, citric acid, NH4Cl
B. subtilis IFO 3335
l-Glutamic acid, glycerol, citric acid, (NH4)2SO4
B. subtilis BL53
l-Glutamic acid, glycerol, citric acid, NH4Cl
B. licheniformis WX02-P43glpFK
Sodium glutamate, crude glycerol, citric acid, NH4Cl
B. methylotrophicus SK19.001
Glycerol, sodium citrate, peptone
B. subtilis C1
Glycerol, citric acid, NH4Cl
B. licheniformis WX-02
Glycerol, sodium citrate, NaNO3, NH4Cl
B. licheniformis BC4
B. licheniformis BC4
CGbase, sodium citrate, NaNO3, NH4Cl
B. licheniformis BC4
CGacid, sodium citrate, NaNO3, NH4Cl
Glycerol has been considered as a promising feedstock for biochemicals production due to its abundance, low price, and high degree of reduction . Much efforts have been made to engineer glycerol catabolism and downstream metabolic pathways for production of various valuable products in E. coli, Corynebacterium glutamicum, Klebsiella species, and Clostridium species [8, 13, 30–32]. Nevertheless, few related studies have been performed in B. licheniformis [17, 33]. As a consequence, the metabolic pathway of B. licheniformis was rationally engineered for efficient production of γ-PGA from crude glycerol (Fig. 1).
In this study, the gene glpK in B. licheniformis was overexpressed by using three different methods to increase glycerol assimilation efficiency. According to our results, introducing an additional copy of glpK was an efficient approach to improve glycerol consumption and γ-PGA yield. With one added copy of glpK on the chromosome, the glpK transcription level and GlpK enzyme activity were increased by 12.06-fold and 19%, respectively, compared to the wild-type strain. The expression levels of glpF and glpD genes were also found to be up-regulated in glpK overexpressing strain due to an increase in Gly3P generated by GlpK (Fig. 5b). These results were consistent with the previous reports [19, 20], and indicated that the GlpK was a rate-limiting step in glycerol metabolism in B. licheniformis. Another limiting step in glycerol metabolism is the necessity for formation of G6P through gluconeogenic pathway. In this study, the limiting enzymes in upper gluconeogenesis pathway were identified in B. licheniformis cultured with glycerol. Over-expression of glpX obviously improved glycerol assimilation efficiency and γ-PGA accumulation, which might be due to the fact that the increasing G6P pool was beneficial for enhancing the PPP flux.
The supply of NADPH, which has been proven to be critical for γ-PGA synthesis, may be another limited factor. Several attempts have been conducted to increase the availability of NADPH by overexpressing glucose 6-phosphate dehydrogenase (encoded by zwf) and 6-phosphogluconate dehydrogenase (encoded by gnd) [4, 21]. The roles of Zwf and Gnd in γ-PGA synthesis from glycerol were evaluated in this study. Our results implied that over-expression of zwf had the positive effect on production. Similarly, over-expression of zwf led to a 35% improvement of γ-PGA titer using glucose as carbon source . Apart from the Zwf and Gnd, the ribulose-5-phosphate 3-epimerase and transketolase were also studied. Although two transketolase isoenzymes are preserved in B. licheniformis, the transketolase encoded by tkt1 was firstly identified as the dominant isoenzyme and overexpressed for improving γ-PGA biosynthesis in B. licheniformis. Also, over-expression of transketolase in E. coli strain could improve the phenylalanine production from glycerol . However, over-expression of rpe had almost no influence on the γ-PGA production and glycerol consumption. These results indicated that more precursors and NADPH might be generated after over-expression of zwf and tkt1.
In this study, the combinatorial effects of over-expression of glpK, glpX, zwf, and tkt1 on γ-PGA titer were evaluated. In the presence of pure glycerol, the recombinant strain BC4 displayed a 49.65% increase in γ-PGA titer and 31.76% increase in glycerol consumption, compared to the wild-type strain. The circuit of the γ-PGA overproduction in the BC4 strain was further investigated. In the present study, the transcriptional levels of genes involved in glycerol catabolic pathway were enhanced in BC4 strain, as well as the GlpK activity and Gly3P accumulation (Fig. 5b), and these results were consistent with the improvement of glycerol consumption. Additionally, TCA cycle is an essential part of γ-PGA synthetic pathways and cell growth . In BC4 strain, more carbon flux was channeled to TCA cycle and glutamate for γ-PGA biosynthesis (Figs. 4 and 5b). The transcription of genes and concentrations of metabolites involved in TCA cycle and glutamate biosynthesis were increased obviously (Fig. 5a, b), further demonstrating that more carbon flux was channeled towards TCA cycle to generate sufficient ATP and glutamate for γ-PGA synthesis.
In NADH-dependent production systems, NADH availability and the ratio of NADH/NAD+ may play a pivotal role in metabolic regulation . Over-expression of glpK, glpX, zwf, and tkt1 was accompanied with the cofactor conversion from NADH to NADPH. The decreasing of NADH might result in reduction of 2,3-BD synthesis. The decreased overflow metabolism may help reducing the carbon waste and boosting carbon flux towards γ-PGA. Moreover, over-expression of glpX, zwf, and tkt1 caused an obvious flux redirection at the G6P node, and enhancement of PPP flux in comparison to the wild-type strain (Fig. 5a, b). These results suggested that combinatorial over-expression of quadruple genes could change the metabolic network, and redistribute metabolic flux to γ-PGA production.
The use of crude glycerol for microbial γ-PGA production not only reduces the cost of γ-PGA production, but also provides a solution for the excess crude glycerol from biodiesel production. This study successfully developed a recombinant strain BC4 for efficient production of γ-PGA from crude glycerol. As shown in Table 2, the final γ-PGA yield was 0.46 g/g glycerol in the glutamic acid-free medium containing crude glycerol, sodium citrate, NaNO3, NH4Cl, and the γ-PGA yield reached 68.2% of theoretical yield, which exhibits a significant increase compared to the wild-type strain. The γ-PGA titer (18.41 g/L) produced by BC4 was slightly lower than those produced by many γ-PGA producers [15, 16, 33, 35–37]; however, the cost of feedstocks was the lowest value compared to other strains. Moreover, compared with the pure glycerol, by using the industry byproduct as feedstock, nearly a half of material cost could be saved (Additional file 3). To the best of our knowledge, the productivity of γ-PGA achieved in this study is the highest value reported using crude glycerol as the substrate [17, 38]. Our results implied that the crude glycerol might be an ideal substitute for glucose for commercial production of γ-PGA using the metabolically engineered B. licheniformis.
The metabolic pathway of B. licheniformis was rewired to improve γ-PGA production from crude glycerol. After systematically investigating the three module pathways genes glpP, tglpFK, glpK, fba, glpX, pgi, zwf, gnd, rpe, tkt1, and tkt2, we confirmed that GlpK, GlpX, Zwf, and Tkt1 were the key enzymes for γ-PGA synthesis in glycerol fermentation. Through combinational expression of these key genes, the titer and productivity of γ-PGA was increased to 19.20 g/L and 0.4 g/L/h, respectively, which were about 1.50-fold of those of the wild-type strain. We also achieved the high γ-PGA titer of 18.4 g/L using crude glycerol as substrate, comparing with other studies reported, with the highest γ-PGA productivity of 0.38 g/L/h. Our studies successfully showed that the crude glycerol can be efficiently converted to γ-PGA by metabolically engineered B. licheniformis, and the developed strategy might be applied to high-level production of other valuable compounds from crude glycerol.
Bacterial strains, media, and culture conditions
Strains and plasmids in this study
B. subtilis 168
Wild-type strain, CCTCC M208065
WX-02 derivative, over-expression of glpK
WX-02 derivative, over-expression of glpP
WX-02 derivative, defective in tglpFK
WX-02 derivative, over-expression of zwf
WX-02 derivative, over-expression of tkt1
WX-02 derivative, over-expression of tkt2
WX-02 derivative, over-expression of rpe
WX-02 derivative, over-expression of gndA
WX-02 derivative, over-expression of pgi
WX-02 derivative, over-expression of fba
WX-02 derivative, over-expression of glpX
WX02-glpK derivative, over-expression of glpX
BC2 derivative, over-expression of zwf
BC3 derivative, over-expression of tkt1
E. coli DH5a
F−Φ80d/lacZΔM15, Δ(lacZYA-argF)U169, recA1, endA1, hsdR17 (r K − , m K + ), phoA, supE44, λ−, thi-1, gyrA96, relA1
Bacillus knockout vector; Kan r
T2(ori)-glpK(A+B+glpK); to over-express glpK
T2(ori)-glpP(A+B+glpP); to over-express glpP
T2(ori)-tglpFK(A+B); to knock out tglpFK
T2(ori)-zwf(A+B+glpK); to over-express glpK
T2(ori)-tkt1(A+B+tkt1); to over-express tkt1
T2(ori)-tkt2 (A+B+tkt2); to over-express tkt2
T2(ori)-rpe(A+B+rpe); to over-express rpe
T2(ori)-gndA(A+B+gndA); to over-express gndA
T2(ori)-pgi(A+B+pgi); to over-express pgi
T2(ori)-glpX(A+B+glpX); to over-express glpX
T2(ori)-fba(A+B+fba); to over-express fba
The fermentation for γ-PGA production was carried out as described previously . The medium for γ-PGA production contained (per liter) 60 g glycerol, 12 g sodium citrate, 15 g NaNO3, 8 g ammonium chloride, 1 g K2HPO4·3H2O, 1 g MgSO4·7H2O, 1 g ZnSO4·7H2O, 1 g CaCl2, and 0.15 g MnSO4·H2O with pH 7.2. B. licheniformis recombinants were cultivated in 250 mL Erlenmeyer’s flasks containing 50 mL medium at 37 °C and 230 rpm for 48 h. The alkaline crude glycerol (CGbase) with 82% purity (w/w) or acid crude glycerol (CGacid) with 72% purity (w/w) was added to reach 60 g/L glycerol concentration. Crude glycerol purification by treatment with activated carbon was performed according to the previous report .
Gene knockout in B. licheniformis
The gene tglpFK was deleted by homologous recombination based on the previous described method . Briefly, the homologous arms of tglpFK were amplified by the corresponding primers and fused by splicing overlap extension PCR (SOE-PCR). The fused fragments were then inserted into T2(2)-ori at the restriction sites XbaI/SacI, generating the recombinant plasmid T2-∆tglpFK. Furthermore, the T2-∆tglpFK was transformed into B. licheniformis WX-02 strain via electroporation. The positive transformants were incubated in LB medium containing kanamycin at 45 °C, and then streaked onto LB agar containing kanamycin incubation to obtain the single-crossover recombinants. The selected colonies were grown in LB medium at 37 °C with serial subcultures to obtain double crossover colonies. The gene-deleted strains were further verified by PCR and DNA-sequencing.
Gene integration in B. licheniformis
The prophage regions were chosen for the integrated expression of genes in B. licheniformis . Briefly, the homologous arms of xkdG, glpX, and amyL terminator were, respectively, amplified from WX-02 genome, and P43 promoter was amplified from B. subtilis 168 genome. The homologous arms of xkdG, P43 promoter, glpX, and amyL terminator were fused by SOE-PCR, then inserted into T2(2)-ori, resulting in the plasmid T2-GglpX. The integration mutants were obtained according to the gene knockout method described in above section (Additional file 5). Other gene integrated expression strains were constructed using the same method.
Transcriptional level analysis
The total RNA was extracted by according to the previous report . The DNA was digested by DNase I enzyme (TaKaRa, Japan), and RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher, USA) was used for first-stand cDNA synthesis following the manufacturer’s instructions. The real-time PCR was performed by using SYBR® Select Master Mix (ABI, USA). The primers listed in Additional file 4 were used for amplifying the corresponding genes, and the 16s rRNA gene was served as the internal reference gene and real-time-PCR was conducted in triplicate for each sample. The expression levels for the recombinant strains were compared to the wild-type strain after normalization.
Determination of enzymatic activities
Cells were harvested in mid-log-phase (cultivated for 20 h), washed with cold potassium phosphate buffer, and then suspended in the same buffer. Cell disruption was achieved by sonication (150 W, 20 kHz, pulse: 2 s on; 2 s off; total: 8 min;), and disruption solution was centrifuged at 12,000×g for 20 min at 4 °C to remove cell debris. The supernatant was stored on 4 °C for subsequent experiments.
Activities of Zwf in crude cell extracts were measured by determining the increase in A340 of NADPH at 30 °C as described before . Glycerol kinase activity was measured using a pyruvate kinase and lactate dehydrogenase coupled assay as described previously . The transketolase activity was measured according to the previous described method . The activity of FBPase was measured following the increase in phosphate concentration liberated using a malachite green/ammonium molybdate solution .
Intracellular metabolites analysis
The cells were collected and analyzed as described previously with minor modifications [43, 44]. Briefly, an appropriate volume of culture broth was harvested at 20 h, and mixed with an equivalent volume of 0.25 M perchloric acid by vortexing transitorily. The volume of transferred cells was adjusted based on OD600 at each sample to satisfy the formula: sampling volume (mL) × OD600 = 10.0 and washed 2 times with 1.0 mL of cold 0.80% NaCl solution. 1.0 mL of pre-cooling 75% ethanol was used for cell metabolites extraction. The samples were treated at 90 °C for 10 min and subsequent at − 40 °C for 5 min, and centrifuged at 4 °C, 12,000×g for 5 min. Supernatants were collected and dried in a rotary vacuum centrifuge device (LABCONCO) for GC–MS analysis. In each experiment, sampling was carried out sevenfold in parallel.
Samples were derivatized before GC–MS analyzing . Firstly, the extracts were oximated with 50 μL of 20 mg/mL methoxyamine hydrochloride (Sigma-Aldrich) in pyridine at 37 °C for 90 min. Before derivatization, 10 μL of 0.5% phenylethyl acetate was added as an analytical internal standard. Subsequently, the acidic protons in samples were derivatized by 37 °C reaction for 30 min with addition of 50 μL of N-methyl-N-(trimethylsilyl)-trifluoroacetamide (MSTFA, Sigma-Aldrich). 1 μL aliquot of derivatized samples was injected into a DB-5MS column (250 μm ID × 30 m, 0.25 μm, Thermo Fisher Scientific) using splitless model. The temperature-programmed profile was set as follows: 50 °C for initial temperature, increase with 10 °C/min to 110 °C and keep for 5 min, then ramp to 165 °C at a rate of 2 °C/min, lastly raise to 220 °C at a rate of 3 °C/min, and hold at 220 °C for 10 min. Helium was used as the carrier gas and the speed of the flow was kept at 1.2 mL/min constantly. For the MS aspect, electron impact ionization (EI) was selected and the ionization was of 70 eV energy with 8000 V acceleration voltage. The MS transfer line temperature and ion source temperature were 280 °C and 300 °C, respectively. The scanned range was 50–650 m/z.
The compounds were identified by comparing mass spectra to National Institute of Standards and Technology (NIST) MS search [version 2.2 (2014)]. TranceFinder 4.1 (Thermo Fisher Scientific) was used for automatic peak picking and calculation of the peak area. The selected peaks were also confirmed manually, and the peak area was normalized by internal standard. The metabolite levels were further compared based on their peak areas of GC–MS chromatogram . Principal component analysis (PCA) and orthogonal partial least squares discrimination analysis (OPLS-DA) were performed by SIMCA 13 (UMETRICS, Sweden) with the metabolome data set (Additional file 6).
Quantification of NADPH and NADH
The intracellular concentrations of NADPH and NADH were determined in the exponential growth phase according to the previous method . For NADH, the reaction mixture (200 μL) contained 50 mM HEPES (pH 7.5), 2 mM EDTA, 120 μM DCPIP, 1 mM PMS, 10 μL alcohol dehydrogenase, and 20 μL disrupted cell solution. Ethanol (15 μL) was then added to the wells to start the reaction. The decrease of A600 was monitored for 5 min. The contents were calculated by reference to the standards run concurrently (0–40 pmol NADH). The same reaction mixture was applied to determine NADPH concentration; however, the reaction was started by the addition of 10 μL glucose-6-phosphate dehydrogenase. The decrease in A600 was monitored for 5 min. The rates were calculated with the standards, similar to the NADH analysis.
Cell growth was measured at 600 nm with a spectrophotometer (Bio-Rad, USA). A standard curve relating OD to cell dry weight (CDW) was developed (1 OD600 = 0.363 g CDW/L). The concentrations of γ-PGA were determined using a gel permeation chromatography (GPC) system equipped with TSK Gel G6000 PWXL column (7.8 mm × 300 mm, Tosoh, Tokyo, Japan) . The mobile phase was a mixture of 25 mM sodium sulfate solution: acetonitrile (8:1) at a flow rate of 0.5 mL/min and detected at 220 nm. The calibration curve for the determination of the γ-PGA concentration was established using purified γ-PGA. Samples were appropriately diluted with ethanol, and then mixed by vortex. After centrifugation, the supernatant was used for subsequent analysis. Citric acid content was determined by a high-performance liquid chromatography (HPLC) using Zorbax SB-Aq (4.6 mm ID × 250 mm, 5 μm) column by means of a UV light detector. A mobile phase of 99% 20 mM Na2HPO4 and 1% acetonitrile (pH 2.0) at a flow rate of 0.5 mL/min was used. The acetoin, 2,3-butanediol (2,3-BD), and acetate concentrations were determined by gas chromatography method as described previously . Glycerol was quantified by HPLC (Agilent, USA) using Inertsil NH2 column (250 mm × 4.6 mm, 5 μm; GL Science Tokyo, Japan) and an evaporative light-scattering detector (ELSD, Agilent, USA). A mobile phase of 80% acetonitrile at 1.0 mL/min flow rate was used .
All experiments were performed in three replicates. The data were represented as the mean value ± standard deviation. All data were conducted to analyze the variance at P < 0.05 and P < 0.01, and a t test was applied to compare the mean values using the Statistica 6.0 software package.
Gly3P: glycerol-3-phosphate; DHAP: dihydroxyacetone phosphate; GAP: glyceraldehyde 3-phosphate; FBP: fructose 1.6-bisphosphate; F6P: fructose 6-phosphate; G6P: glucose-6-phosphate; 6PG: 6-phosphogluconate; R5P: ribulose 5-phosphate; Ru5P: ribose 5-phosphate; X5P: xylulose 5-phosphate; S7P: sedoheptulose 7-phosphate; E4P: erythrose 4-phosphate; 1,3 BPG: 1,3-bisphosphoglyceric acid; 3PG: 3-phosphoglyceric acid; PEP: phosphoenolpyruvate; PYR: pyruvate; AceCoA: acetyl-CoA; CIT: citrate; ICIT: isocitrate; AKG: 2-oxoglutarate; SUC-CoA: succinyl-CoA; SUC: succinate; FUM: fumarate; MAL: malate; OAA: oxaloacetate; Glu: glutamate; γ-PGA: poly-γ-glutamic acid; AC: acetoin; 2,3-BD: 2,3-butanediol.
Enzymes and others
GlpP: glycerol uptake operon antiterminator; GlpF: glycerol transport facilitator; GlpK: glycerol kinase; GlpD: glycerol-3-phosphate dehydrogenase; Fba: fructose-bisphosphate aldolase; GlpX: fructose-1,6-bisphosphatase; Pgi: glucose-6-phosphate isomerase; Zwf: glucose-6-phosphate dehydrogenase; GndA: 6-phosphogluconate dehydrogenase; Rpe: d-ribulose-5-phosphate 3-epimerase; Tkt: transketolase; PYC: pyruvate carboxylase; CitZ: citrate synthase; CitB: aconitate hydratase; ICDH: isocitrate dehydrogenase; OGDC: α-oxoglutarate dehydrogenase; SucCD: succinyl-CoA ligase; SDH: succinate dehydrogenase; FH: fumarate hydratase; Mdh: malate dehydrogenase; GDH: glutamate dehydrogenase; GOGAT: glutamate synthase; PTA: phosphate acetyltransferase; AK: acetate kinase; AlsS: acetolactate synthase; AlsD: alpha-acetolactate decarboxylase; 2,3-BDH: 2,3-butanediol dehydrogenase; GDH: glycerol dehydrogenase; PPP: pentose phosphate pathway; TCA: tricarboxylic acid; CGacid: acid crude glycerol; CGbase: alkaline crude glycerol; CCR: carbon catabolite repression; MSTFA: N-methyl-N-(trimethylsilyl)-trifluoroacetamide; PCA: principal component analysis; OPLS-DA: orthogonal partial least squares discrimination analysis; CDW: cell dry weight; GPC: gel permeation chromatography.
YZ and SC designed and supervised the study. YZ, BS, HW, and JS performed the experiments. YZ, BS, DC, LY, SY, ZW, XM, and SC analyzed the data and wrote the manuscript. All authors read and approved the final manuscript.
This study was funded by the National Program on Key Basic Research Project (973 Program, No. 2015CB150505), the China National Tobacco Corporation, Grant (110201502018), and the Technical Innovation Special Fund of Hubei Province (2018ACA149).
The authors declare that they have no competing interests.
Availability of supporting data
Data supporting the results of the article are included within this manuscript.
Consent for publication
All authors have approved the manuscript for publication.
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Luo Z, Guo Y, Liu J, Qiu H, Zhao M, Zou W, Li S. Microbial synthesis of poly-gamma-glutamic acid: current progress, challenges, and future perspectives. Biotechnol Biofuels. 2016;9:134.View ArticleGoogle Scholar
- Buescher JM, Margaritis A. Microbial biosynthesis of polyglutamic acid biopolymer and applications in the biopharmaceutical, biomedical and food industries. Crit Rev Biotechnol. 2007;27:1–19.View ArticleGoogle Scholar
- Bajaj I, Singhal R. Poly (glutamic acid)—an emerging biopolymer of commercial interest. Bioresour Technol. 2011;102:5551–61.View ArticleGoogle Scholar
- Cai D, He P, Lu X, Zhu C, Zhu J, Zhan Y, Wang Q, Wen Z, Chen S. A novel approach to improve poly-γ-glutamic acid production by NADPH regeneration in Bacillus licheniformis WX-02. Sci Rep. 2017;7:43404.View ArticleGoogle Scholar
- Feng J, Gu Y, Quan Y, Cao M, Gao W, Zhang W, Wang S, Yang C, Song C. Improved poly-gamma-glutamic acid production in Bacillus amyloliquefaciens by modular pathway engineering. Metab Eng. 2015;32:106–15.View ArticleGoogle Scholar
- Chen Z, Liu D. Toward glycerol biorefinery: metabolic engineering for the production of biofuels and chemicals from glycerol. Biotechnol Biofuels. 2016;9:205.View ArticleGoogle Scholar
- Dobson R, Gray V, Rumbold K. Microbial utilization of crude glycerol for the production of value-added products. J Ind Microbiol Biotechnol. 2012;39:217–26.View ArticleGoogle Scholar
- Chen X, Tian K, Niu D, Shen W, Algasan G, Singh S, Wang Z. Efficient bioconversion of crude glycerol from biodiesel to optically pure d-lactate by metabolically engineered Escherichia coli. Green Chem. 2014;16:342–50.View ArticleGoogle Scholar
- Wang Y, Tao F, Xin B, Liu H, Gao Y, Zhou NY, Xu P. Switch of metabolic status: redirecting metabolic flux for acetoin production from glycerol by activating a silent glycerol catabolism pathway. Metab Eng. 2017;39:90–101.View ArticleGoogle Scholar
- Blankschien MD, Clomburg JM, Gonzalez R. Metabolic engineering of Escherichia coli for the production of succinate from glycerol. Metab Eng. 2010;12:409–19.View ArticleGoogle Scholar
- Saini M, Wang ZW, Chiang CJ, Chao YP. Metabolic engineering of Escherichia coli for production of n-butanol from crude glycerol. Biotechnol Biofuels. 2017;10:173.View ArticleGoogle Scholar
- Kim K, Kim SK, Park YC, Seo JH. Enhanced production of 3-hydroxypropionic acid from glycerol by modulation of glycerol metabolism in recombinant Escherichia coli. Bioresour Technol. 2014;156:170–5.View ArticleGoogle Scholar
- Gonzalez-Pajuelo M, Meynial-Salles I, Mendes F, Andrade JC, Vasconcelos I, Soucaille P. Metabolic engineering of Clostridium acetobutylicum for the industrial production of 1,3-propanediol from glycerol. Metab Eng. 2005;7:329–36.View ArticleGoogle Scholar
- Zhu C, Nomura CT, Perrotta JA, Stipanovic AJ, Nakas JP. Production and characterization of poly-3-hydroxybutyrate from biodiesel-glycerol by Burkholderia cepacia ATCC 17759. Biotechnol Prog. 2010;26:424–30.PubMedGoogle Scholar
- Cromwick AM, Birrer GA, Gross RA. Effects of pH and aeration on gamma-poly(glutamic acid) formation by Bacillus licheniformis in controlled batch fermentor cultures. Biotechnol Bioeng. 1996;50:222–7.View ArticleGoogle Scholar
- Shih IL, Van YT, Chang YN. Application of statistical experimental methods to optimize production of poly (γ-glutamic acid) by Bacillus licheniformis CCRC 12826. Enzyme Microb Technol. 2002;31:213–20.View ArticleGoogle Scholar
- Zhan Y, Zhu C, Sheng B, Cai D, Wang Q, Wen Z, Chen S. Improvement of glycerol catabolism in Bacillus licheniformis for production of poly-gamma-glutamic acid. Appl Microbiol Biotechnol. 2017;101:7155–64.View ArticleGoogle Scholar
- Darbon E, Servant P, Poncet S, Deutscher J. Antitermination by GlpP, catabolite repression via CcpA and inducer exclusion triggered by P-GlpK dephosphorylation control Bacillus subtilis glpFK expression. Mol Microbiol. 2002;43:1039–52.View ArticleGoogle Scholar
- Gottlieb K, Albermann C, Sprenger GA. Improvement of l-phenylalanine production from glycerol by recombinant Escherichia coli strains: the role of extra copies of glpK, glpX, and tktA genes. Microb Cell Fact. 2014;13:96.View ArticleGoogle Scholar
- Hao G, Chen H, Gu Z, Zhang H, Chen W, Chen YQ. Metabolic engineering of Mortierella alpina for arachidonic acid production with glycerol as carbon source. Microb Cell Fact. 2015;14:205.View ArticleGoogle Scholar
- Kim YM, Cho HS, Jung GY, Park JM. Engineering the pentose phosphate pathway to improve hydrogen yield in recombinant Escherichia coli. Biotechnol Bioeng. 2011;108:2941–6.View ArticleGoogle Scholar
- Tian G, Wang Q, Wei X, Ma X, Chen S. Glutamate dehydrogenase (RocG) in Bacillus licheniformis WX-02: enzymatic properties and specific functions in glutamic acid synthesis for poly-gamma-glutamic acid production. Enzyme Microb Technol. 2017;99:9–15.View ArticleGoogle Scholar
- Kumar C, Igbaria A, D’Autreaux B, Planson AG, Junot C, Godat E, Bachhawat AK, Delaunay-Moisan A, Toledano MB. Glutathione revisited: a vital function in iron metabolism and ancillary role in thiol-redox control. EMBO J. 2011;30:2044–56.View ArticleGoogle Scholar
- Khanna S, Goyal A, Moholkar VS. Microbial conversion of glycerol: present status and future prospects. Crit Rev Biotechnol. 2012;32:235–62.View ArticleGoogle Scholar
- Wang J, Yuan H, Wei X, Chen J, Chen S. Enhancement of poly-γ-glutamic acid production by alkaline pH stress treatment in Bacillus licheniformis WX-02. J Chem Technol Biotechnol. 2016;91:2399–403.View ArticleGoogle Scholar
- Spaans SK, Weusthuis RA, van der Oost J, Kengen SW. NADPH-generating systems in bacteria and archaea. Front Microbiol. 2015;6:742.View ArticleGoogle Scholar
- Yao R, Xiong D, Hu H, Wakayama M, Yu W, Zhang X, Shimizu K. Elucidation of the co-metabolism of glycerol and glucose in Escherichia coli by genetic engineering, transcription profiling, and 13C metabolic flux analysis. Biotechnol Biofuels. 2016;9:175.View ArticleGoogle Scholar
- Wang Z, Chen T, Ma X, Shen Z, Zhao X. Enhancement of riboflavin production with Bacillus subtilis by expression and site-directed mutagenesis of zwf and gnd gene from Corynebacterium glutamicum. Bioresour Technol. 2011;102:3934–40.View ArticleGoogle Scholar
- Yao R, Pan K, Peng H, Feng L, Hu H, Zhang X. Engineering and systems-level analysis of Pseudomonas chlororaphis for production of phenazine-1-carboxamide using glycerol as the cost-effective carbon source. Biotechnol Biofuels. 2018;11:130.View ArticleGoogle Scholar
- Oh BR, Seo JW, Heo SY, Hong WK, Luo LH, Joe MH, Park DH, Kim CH. Efficient production of ethanol from crude glycerol by a Klebsiella pneumoniae mutant strain. Bioresour Technol. 2011;102:3918–22.View ArticleGoogle Scholar
- Meiswinkel TM, Rittmann D, Lindner SN, Wendisch VF. Crude glycerol-based production of amino acids and putrescine by Corynebacterium glutamicum. Bioresour Technol. 2013;145:254–8.View ArticleGoogle Scholar
- Gonzalez-Pajuelo M, Meynial-Salles I, Mendes F, Soucaille P, Vasconcelos I. Microbial conversion of glycerol to 1,3-propanediol: physiological comparison of a natural producer, Clostridium butyricum VPI 3266, and an engineered strain, Clostridium acetobutylicum DG1(pSPD5). Appl Environ Microbiol. 2006;72:96–101.View ArticleGoogle Scholar
- Jeong JH, Kim JN, Wee YJ, Ryu HW. The statistically optimized production of poly(γ-glutamic acid) by batch fermentation of a newly isolated Bacillus subtilis RKY3. Bioresour Technol. 2010;101:4533–9.View ArticleGoogle Scholar
- Yu W, Chen Z, Ye H, Liu P, Li Z, Wang Y, Li Q, Yan S, Zhong C, He N. Effect of glucose on poly-gamma-glutamic acid metabolism in Bacillus licheniformis. Microb Cell Fact. 2017;16:22.View ArticleGoogle Scholar
- Jung DY, Jung S, Yun JS, Kim JN, Wee YJ, Jang HG, Ryu HW. Influences of cultural medium component on the production of poly (γ-glutamic acid) by Bacillus sp. RKY3. Biotechnol Bioprocess Eng. 2005;10:289–95.View ArticleGoogle Scholar
- Peng Y, Jiang B, Zhang T, Mu W, Miao M, Hua Y. High-level production of poly(γ-glutamic acid) by a newly isolated glutamate-independent strain, Bacillus methylotrophicus. Process Biochem. 2015;50:329–35.View ArticleGoogle Scholar
- Du G, Yang G, Qu Y, Chen J, Lun S. Effects of glycerol on the production of poly(γ-glutamic acid) by Bacillus licheniformis. Process Biochem. 2005;40:2143–7.View ArticleGoogle Scholar
- Shih IL, Wu PJ, Shieh CJ. Microbial production of a poly(γ-glutamic acid) derivative by Bacillus subtilis. Process Biochem. 2005;40:2827–32.View ArticleGoogle Scholar
- Yangtse W, Zhou Y, Lei Y, Qiu Y, Wei X, Ji Z, Qi G, Yong Y, Chen L, Chen S. Genome sequence of Bacillus licheniformis WX-02. J Bacteriol. 2012;194:3561–2.View ArticleGoogle Scholar
- Qiu Y, Zhang J, Li L, Wen Z, Nomura CT, Wu S, Chen S. Engineering Bacillus licheniformis for the production of meso-2,3-butanediol. Biotechnol Biofuels. 2016;9:117.View ArticleGoogle Scholar
- Cai D, Wei X, Qiu Y, Chen Y, Chen J, Wen Z, Chen S. High-level expression of nattokinase in Bacillus licheniformis by manipulating signal peptide and signal peptidase. J Appl Microbiol. 2016;121:704–12.View ArticleGoogle Scholar
- Klaus A, Pfirrmann T, Glomb MA. Transketolase A from E. coli significantly suppresses protein glycation by glycolaldehyde and glyoxal in vitro. J Agric Food Chem. 2017;65:8196–202.View ArticleGoogle Scholar
- Li X, Long D, Ji J, Yang W, Zeng Z, Guo S, Ji Z, Qi G, Chen S. Sample preparation for the metabolomics investigation of poly-gamma-glutamate-producing Bacillus licheniformis by GC-MS. J Microbiol Methods. 2013;94:61–7.View ArticleGoogle Scholar
- Guan Y, Yin D, Du X, Ye X. Functional metabolomics approach reveals the reduced biosynthesis of fatty acids and TCA cycle is required for pectinase activity in Bacillus licheniformis. J Ind Microbiol Biotechnol. 2018. https://doi.org/10.1007/s10295-018-2071-z.View ArticlePubMedGoogle Scholar
- Goto A, Kunioka M. Biosynthesis and hydrolysis of poly(gamma-glutamic acid) from Bacillus subtilis IF03335. Biosci Biotechnol Biochem. 1992;56:1031–5.View ArticleGoogle Scholar
- da Silva GP, Mack M, Contiero J. Glycerol: a promising and abundant carbon source for industrial microbiology. Biotechnol Adv. 2009;27:30–9.View ArticleGoogle Scholar
- da Silva SB, Cantarelli VV, Ayub MAZ. Production and optimization of poly-γ-glutamic acid by Bacillus subtilis BL53 isolated from the Amazonian environment. Bioprocess Biosyst Eng. 2014;37:469–79.View ArticleGoogle Scholar