Effects of glutamate addition on cell growth and γ-PGA synthesis
To assess the effect of glutamate addition on fermentation, B. subtilis SCP017-03 was cultured under two conditions: one with glutamate addition (experimental group: GA) and another without glutamate addition (control group: CK). The differences in glucose and glutamate consumption, cell growth, and γ-PGA yields were compared. As depicted in Fig. 1a, the results revealed no significant differences in cell growth or carbon source consumption between both groups during the initial 0–4 h of fermentation. However, after 10 h, cell growth in the GA group lagged significantly behind that in the CK group. Between 18 and 48 h of fermentation, the differences in cell growth and glucose consumption between the experimental and control groups were more pronounced. At the 36-hour point, the OD600 value difference between the two groups reached the maximum, with CK in the control group at 6.4, while GA in the experimental group measured 4.8. After 12 h, glucose consumption in the GA group slowed, with the glucose content stabilizing at 25 g/L after 72-hour fermentation. However, the control group exhibited faster sugar consumption, exhausting the glucose supply after 60th hour of fermentation. Sugar consumption correlated precisely with differences in cell growth. Furthermore, the GA group initiated γ-PGA synthesis after 8 h of fermentation and continued to accumulate. After 18 h, the γ-PGA content reached 16.7 g/L, peaking at 36.5 g/L after 48 h of fermentation, remaining relatively stable until the end of fermentation. The residual concentration of glutamate in the fermentation broth of the GA group gradually decreased from 50 g/L to 35 g/L after 60 h of fermentation, indicating a modest reduction of 15 g/L over the course of fermentation. This observation suggested that the bacterial cells of this strain synthesized a certain amount of glutamate.
The γ-PGA generated in the fermentation broth was assessed using SDS-PAGE (Fig. 1b). The figure illustrates that the molecular weight of γ-PGA produced by this strain significantly exceeded 200 kD. However, as fermentation progressed, there was a slight reduction in molecular weight.
Transcriptome analysis of strain scp017-03 in response to glutamate addition
The RNA samples from each group were subjected to sequencing using an Illumina NovaSeq 6000 high-throughput sequencing platform. The sequencing output data exhibited base quality values of Q30 and Q20 exceeding 95% and 98%, respectively, with a sample sequencing error rate of less than 0.03%. These metrics complied with stringent quality standards for Illumina sequencing. Table 1 shows the alignment results for each sample against the reference genome. All six samples achieved a total mapping rate, were successfully aligned to the genome via sequence mapping, and each exhibited unique alignment positions on the reference sequence.
Table 1
Mapping results of strain SCP017-0 transcriptome
Sample.name | CK1 | CK2 | CK3 | GA1 | GA2 | GA3 |
Total reads | 8754356 | 7523564 | 7800772 | 7724204 | 7734902 | 7438034 |
Total mapped reads | 8621421 | 7418727 | 7684771 | 7636601 | 7638796 | 7325040 |
Uniquely mapped reads | 8425126 | 7275289 | 7526589 | 7500556 | 7495239 | 7168957 |
Multiple mapped reads | 196295 | 143438 | 158182 | 136045 | 143557 | 156083 |
Total mapping rate | 98.48% | 98.61% | 98.51% | 98.87% | 98.76% | 98.48% |
Uniquely mapping rate | 96.24% | 96.7% | 96.49% | 97.1% | 96.9% | 96.38% |
Multiple mapping rate | 2.24% | 1.91% | 2.03% | 1.76% | 1.86% | 2.1% |
Statistical analysis of differentially expressed genes was conducted using the following screening criteria: log2 (FoldChange) > 0 or padj < 0.05. A total of 1867 differentially expressed genes were identified, comprising 862 up-regulated and 1005 down-regulated genes. Subsequently, ClusterProfiler software was utilized to perform GO function enrichment analysis and KEGG pathway enrichment analysis on the sets of differential genes, with a significance threshold set at padj < 0.05. GO function enrichment analysis indicated the presence of both up-regulated and down-regulated genes within three categories: biological process (BP), cellular component (CC), and molecular function (MF) (Fig. 2). These findings suggested that the addition of exogenous glutamic acid significantly affected various metabolic pathways in this strain. Furthermore, the enrichment results of the KEGG pathway analysis encompassed the differentially expressed genes and certain genes related to γ-PGA synthesis (Fig. 3).
The phosphotransferase system (PTS) is the primary mechanism employed by bacteria to uptake hexose, hexitol, disaccharides, and other carbohydrates. As shown in Fig. 3, the addition of exogenous glutamate led to the differential expression of ybbF and sacX genes within the sucrose metabolism system (EIIBCA or EIIBC component). Their log2FC values were − 2.03 and + 2.65, respectively. The carbohydrates transported into the membrane entered not only the central metabolic pathway but also the polysaccharide synthesis (EPS). Within the EMP pathway consisting of energy consumption and productivity stages, gene expression in the energy consumption stage, including pgi (glucose-6-phosphate isomerase), glpX (fructose 1,6-diphosphatase), and pfkA (ATP-dependent phosphofructokinase), was slightly down-regulated. Conversely, during the production stage of Pyr synthesis from G-3-P, the expression of these five genes was significantly up-regulated. These genes include fbaA (fructose-diphosphate aldolase), pgm (glycerophosphate mutase), eno (enolase), gapA, and gapB (glyceraldehyde-3-phosphate dehydrogenase). Consequently, the addition of exogenous glutamate enhanced the EMP pathway. In the PPP pathway, the expression of zwf (glucose-6-phosphate dehydrogenase) and ykgB (glucose-6-phosphate esterase) was down-regulated by factors of 1.07 and 2.03, respectively, while ganZ (6-glucose-6-phosphate dehydrogenase) witnessed an up-regulation in expression.
The TCA cycle is initiated by the condensation of acetyl-CoA and oxaloacetic acid to produce citric acid. Acetyl-CoA is formed through the oxidative decarboxylation of pyruvate under aerobic conditions. Figure 4 illustrates that most genes within the TCA cycle pathway were down-regulated in the GA group. Notably, the log2FC values for citZ and icd were − 1.76 and − 0.72, respectively, with only the sdhABC gene (encoding succinate dehydrogenase) displaying up-regulation. In addition, the coding genes responsible for the pyruvate dehydrogenase complex, which catalyzed the decarboxylation of pyruvate to synthesize acetyl coenzyme A, were down-regulated. For example, the acoC gene (encoding dihydrooctylamide transacetylase) exhibited a down-regulation of 4.69 times, while the idh gene (encoding lactate dehydrogenase), responsible for transforming pyruvate into lactic acid, displayed an up-regulation of 2.74 times. These findings suggested that the synthesis of high-viscosity γ-PGA by cells resulted in reduced dissolved oxygen levels in the culture solution, thereby reducing the flow of substances into the tricarboxylic acid pathway.
L-glutamic acid was converted into D-glutamic acid through the action of RacE (glutamate racemase), and these compounds were subsequently utilized in the synthesis of γ-PGA by γ-PGA synthetase. Figure 3 demonstrates an up-regulation in the expression of the glnA gene (encoding glutamine synthetase) and racE, with log2FC values of 0.86 and 1.77, respectively. This upregulation enhanced the production of the precursor D-glutamic acid for γ-PGA, thereby strengthening γ-PGA synthesis. Notably, the pgsBCA gene cluster (encoding γ-PGA synthetase) was significantly up-regulated, with log2FC values of 2.63, 2.83, and 2.78, respectively. The pgsBCA operon was regulated by DegS-DegU and ComA-ComP two-component systems. In the CK group without γ-PGA synthesis, the FPKM value of the degU gene was 2396, and the FPKM values of the pgsBCA gene cluster were 859, 744, and 1243, indicating baseline expression of the γ-PGA synthesis system in the absence of exogenous glutamate. However, the pgsBCA operon was further activated upon addition of exogenous glutamate, with the FPKM value of the pgsA gene reaching 8533. These two-component systems can interact with phosphorylated ComA to induce degQ expression. Nevertheless, the expression levels of comA, comX, comQ, and degQ genes in this strain remained low. For example, the FPKM value of degQ in the GA group was only 69. After γ-PGA synthesis, the cells required the action of γ-PGA hydrolase to release γ-PGA. This process may involve three genes, pgdS, cwlO, and ggt. Transcriptome analysis revealed that these three genes exhibited low expression in the GA group, with FPKM values of 27, 282, and 266, respectively. Furthermore, the expression of ggt in the CK group was significantly higher than that in the GA group.
qPCR analysis of genes related to γ-PGA biosynthesis
Selected genes related to γ-PGA biosynthesis were subjected to qPCR analysis to validate the findings of transcriptome analysis. To investigate dynamic changes in key genes and regulatory factors during fermentation, samples were collected for qPCR analysis at three time points: 10 h (immediately following γ-PGA synthesis initiation), 20 h, and 48 h (when γ-PGA accumulation peaked). The results are presented in Table 2. After 10-hour fermentation, most genes in the TCA cycle, PPP pathway, and glutamic acid synthesis pathway exhibited slight down-regulation, except for a few genes such as citA, gltP, and rocA. In the γ-PGA synthesis module, only the expression of the regulator DegQ showed a 0.67-fold decrease, and the genes encoding γ-PGA synthetase were significantly up-regulated. Specifically, pgsB expression increased by 4.94 times, and pgsA expression increased by 4.63 times. Moreover, the genes pgdS and cwlO, associated with γ-PGA hydrolysis, were up-regulated, with fold changes of 2.73 and 1.47, respectively. At the 20-hour point, genes in the central metabolic pathway generally showed down-regulation, except for the eno gene in the EMP pathway and the sdhA gene in the TCA cycle, consistent with the transcriptome analysis results. Notably, pgdS exhibited slight downregulation. In addition, the γ-PGA synthesis module consistently displayed up-regulated expression. In the late stage of fermentation (48 h), owing to significant nutrient depletion in the culture medium, gene expression in the central pathway decreased. Concurrently, as γ-PGA accumulation peaked, some γ-PGA was hydrolyzed by bacterial cells to maintain its state, as evidenced by the substantial upregulation of cwlO gene expression, with a fold change of 21.93. Throughout the fermentation process, the regulatory factor degQ was consistently down-regulated, whereas degS and degU were up-regulated, indicating complex interactions within the DegS-DegU two-component system.
Table 2
Comparison of expression of key genes in γ-PGA synthesis by qPCR of strain SCP017-03 grown with glutamate and without glutamate at different fermentation times
Gene | 10h | 20h | 48h |
pfkA | 4.30 | 0.60 | 2.18 |
pgi | 2.93 | 0.75 | 1.63 |
eno | 6.30 | 2.97 | 2.18 |
pdhA | 2.52 | 0.05 | 0.11 |
citB | 1.87 | 1.99 | 0.54 |
icd | 1.50 | 0.40 | 0.53 |
fumC | 2.69 | 0.44 | 0.43 |
ykgB | 1.06 | 0.23 | 0.22 |
zwf | 1.97 | 0.88 | 3.34 |
ganZ | 2.80 | 0.47 | 0.20 |
gltP | 0.67 | 0.16 | 0.04 |
rocA | 0.20 | 2.86 | 5.96 |
racE | 4.42 | 0.35 | 0.01 |
degQ | 0.67 | 0.32 | 0.05 |
degU | 1.36 | 1.02 | 1.41 |
degS | 7.59 | 2.03 | 1.13 |
pgsA | 4.63 | 2.45 | 18.83 |
pgsB | 4.94 | 4.82 | 21.26 |
pgsC | 3.23 | 9.38 | 12.21 |
glnA | 0.92 | 1.64 | 0.94 |
putM | 1.55 | 0.03 | 0.06 |
pgdS | 2.73 | 0.76 | 1.71 |
comA | 2.45 | 0.58 | 6.70 |
cwlO | 1.47 | 4.89 | 21.93 |
Notes: The numbers in the table represent gene expression levels, which have been normalized to the 16S rRNA internal reference gene. |
Proteomic analysis of strain SCP017-03 in response to glutamic acid addition
Strain SCP017-03 was subjected to fermentation in culture media with or without exogenous glutamate (GA group). Cell samples were collected at 6 and 20 h after fermentation. Protein analysis was conducted using 4D-label-free technology, leading to the identification of 2,161 proteins and 13,243 peptides. Each sample (CK6h, GA6h, CK20h, and GA20h) was then quantified. Differentially expressed proteins were identified based on criteria where the FC > 2.0 or < 0.5 and Pvalue < 0.05, with a significance level of P < 0.05. After 6 h of fermentation, a comparison between the GA and CK groups (Fig. 4a) revealed 328 differentially expressed proteins. These included 162 up-regulated and 166 down-regulated proteins. After 20 h of fermentation (Fig. 4b), 279 differentially expressed proteins were identified, comprising 142 up-regulated upregulated and 137 down-regulated proteins. A comparison within the GA group, specifically between 20 and 6 h of fermentation (Fig. 4c), yielded 297 differentially expressed proteins. Among them, 158 were up-regulated and 139 were down-regulated. Similarly, 296 differentially expressed proteins were observed in the CK group (Fig. 4d), with 162 up-regulated and 136 down-regulated proteins. The volcano plot of the differential proteins visually illustrated that the majority of proteins (represented in black) exhibited no significant differences. Proteins with differences in expression are denoted in red (up-regulated) and green(down-regulated) (Fig. 4).
The enrichment category of differential proteins GO in the GA group exhibited significant changes between the 6-hour and 20-hour fermentation periods, whereas the corresponding changes in the CK group were relatively minor (Fig. 5a and 5b). Compared to differentially expressed proteins associated with BP in the CK group, those in the GA group at 6 h of fermentation were primarily involved in processes related to localization (25 proteins), transportation (21 proteins), and coenzyme biosynthesis (6 proteins). Furthermore, differentially expressed proteins linked to “cell composition (CC)” were predominantly associated with specific cellular structures. As the fermentation duration reached 20 h, the BP processes in the GA group emphasized transmembrane transport (13 proteins), transport (23 proteins), and localization (24 proteins), while the CC processes centered on membrane-related components (28 proteins) and the ABC transport complex (2 proteins). In the MF category, the focus shifted towards hydrolase activity (four proteins). Similarly, the analysis of GO enrichment for differential proteins in both the GA and CK groups (Fig. 5c and 5d) revealed that the differentially expressed proteins in the GA group were primarily concentrated in the BP (22 localized proteins) and CC (30 membrane-related proteins) categories, whereas in the CK group, the emphasis was mainly on the CC (28 membrane-related proteins) and MF (38 proteins associated with hydrolase activities) categories.
Table 2 illustrates that after 6 h of fermentation, the protein expression within the central metabolic pathways of the GA group exhibited slight changes compared to the CK group. For example, in the EMP pathway, with the exception of GlpX (fructose 1,6-diphosphatase), which showed a modest 0.31-fold up-regulation, other proteins displayed down-regulated expression. Notably, Pgi, Fbp, and Eno had log2FC values of -0.58, -0.93, and − 0.82, respectively. In the TCA cycle, only PdhB (β subunit of pyruvate dehydrogenase E1 component) and PDHC (α subunit of DNA polymerase III) exhibited slight increases in protein expression, with log2FC values of 0.47 and 0.59, respectively. Conversely, the expression of other proteins, including SdhA (protein subunit of succinate dehydrogenase), notably decreased, with a log2FC value of -1.25. Changes in protein expression within the PPP pathway were relatively minimal. After 20-hour fermentation, the log2FC values for certain proteins in these metabolic pathways underwent substantial changes, with most proteins in the EMP pathway showing increased expression. In the TCA cycle, aside from PdhABC, which was up-regulated by 0.71, 1.38, and 1.16, the expression of other proteins aligned with the transcriptome results, indicating down-regulation. Similarly, the expression of PPP pathway proteins was consistent with the results of the transcriptome analysis. Furthermore, this study analyzed protein expression at different stages of fermentation (20 and 6 h) for different treatments. The results revealed that in the GA group, protein expression in the EMP pathway was predominantly up-regulated, whereas the CK group exhibited the opposite trend. However, enzyme expression in the TCA cycle and PPP pathway remained largely consistent between the two groups, with PPP pathway enzymes being down-regulated.
After a 6-hour fermentation period, a comparison between the GA and CK groups revealed notable changes in the expression of proteins involved in the glutamate synthesis pathway (Table 2). Specifically, the expressions of RacE and RocG in the glutamate synthesis pathway were up-regulated in the GA group, with log2FC values of 1.78 and 0.15, respectively. In contrast, the expression of other proteins within this pathway decreased, with GltB (glutamate synthase domain 3) showing the most significant down-regulation, marked by a log2FC value of -3.28. However, proteins and regulatory factors related to γ-PGA synthesis exhibited only slight changes, such as the log2FC values of -0.39 for DegQ, -0.78, DegU, and − 0.19 for PgsA. However, after a 20-hour fermentation period, substantial changes were observed in the expression of histones in both the GA and CK groups. In the glutamic acid synthesis pathway, the log2FC value of PurQ increased from − 0.29 to 1.89, while RacE exhibited an increase in log2FC from 1.76 to 4.12. Concurrently, the expression of other proteins in this pathway was also down-regulated. Within the γ-PGA synthesis system, both PgsA and PgsB were up-regulated, with log2FC values of 1.68 and 3.05, respectively. Notably, the protein expression of DegQ and DegS in the DegS-DegU two-component system diverged from the transcriptome analysis. Specifically, the log2FC values for regulatory factors DegQ and DegS were 1.9 and − 0.44, respectively. However, the protein expression of the ComP-ComA system aligned with aforementioned transcriptional analysis, with the ComA regulatory factor exhibiting down-regulation. Additionally, there was a general increase in enzymatic activity within the glutamate synthesis pathway. For example, the activity of glutamate racemase (RacE), which was responsible for the transformation of D-glutamic acid and L-glutamic acid, increased by 4.12 times, while the expression of γ-PGA hydrolase (Pgds) decreased by 2.24 times.
Table 2 presents a significant increase in γ-PGA polymerase activity within the GA group after 20 h compared to the 6-hour mark. Notably, the log2FC values for PgsA and PgsB were 2.04 and 3.83, respectively. Furthermore, it was observed that γ-PGA synthetase activity in the CK group after 20 h of fermentation was lower than that at the 6-hour mark. γ-PGA synthetase activity notably increased by a factor of 2.05. This observation suggested that the expression of the γ-PGA synthetase gene exogenous glutamic acid. However, the translation of this mRNA into protein may have been influenced by other genes or regulatory factors, resulting in the reduced activity of γ-PGA synthetase.
Optimization of fermentation conditions
Preliminary optimizations were conducted for carbon and nitrogen sources in an effort to enhance γ-PGA yield during fermentation. The results of carbon source optimization are presented in Fig. 6a. Among the tested options, glucose as the carbon source yielded the highest γ-PGA production at 38.35 g/L, followed by sucrose (36.5 g/L), starch (33.9 g/L), and xylose (32.6 g/L). Conversely, when glycerol was employed as the carbon source, the γ-PGA yield was the lowest at 22.8 g/L. Figure 6b illustrates that yeast powder was a superior nitrogen source, yielding a γ-PGA production of 40.12 g/L. Other nitrogen sources, including pancreatic peptone, soybean peptone, and bean cake powder, yielded approximately 37.5 g/L γ-PGA. In contrast, beef paste and corn steep liquor, when applied as nitrogen sources, resulted in lower yields of 14.32 g/L and 16.2 g/L, respectively. Consequently, the optimal carbon source was glucose, and the preferred nitrogen source was yeast powder.
The results revealed that during the synthesis of γ-PGA, the transcriptional activity of related genes, including those associated with glycolysis, tricarboxylic acid cycle, and PPP pathway, decreased. This decrease could be related to the increased fermentation broth viscosity and the impact of dissolved oxygen. To optimize fermentation, four organic carriers (w/v: 0.5%) were selected. The organic carriers included n-hexane, n-heptane, n-dodecane, and n-hexadecane. Figure 7a illustrates that when n-hexadecane was introduced at the onset of fermentation, the highest γ-PGA yield was achieved, reaching 42.9 g/L. This was followed by n-hexane and n-heptane, with γ-PGA yields of 38.2 g/L and 36.2 g/L, respectively. In contrast, n-dodecane afforded the lowest γ-PGA yield, only reaching 9.24 g/L. Considering the potential impact of organic matter added at the early stages of fermentation on bacterial growth, the organic carrier was introduced into the fermentation for 24 h. Figure 7b shows that, with the addition of the organic carrier after 24-hour fermentation, the γ-PGA production increased by 9.4 g/L in the n-heptane group and by 17.06 g/L in the n-dodecane group. Among these, the n-heptane group achieved the highest γ-PGA output of 45.6 g/L. Notably, there was no significant difference observed between the n-hexadecane group and the n-hexane group, with respective outputs of 43.3 g/L and 35.9 g/L. Therefore, n-heptane was selected for liquid fermentation for 24 h.
Scale-up experiment of batch feeding in 5L fermentor
In the context of shaking flask conditions, the fermentation duration was insufficient, resulting in a low dissolved oxygen concentration and inadequate carbon source concentration in the later stages. These limitations restricted the successful synthesis of γ-PGA. Therefore, fermentation processes were optimized using a 5 L fermentor. The fermentation parameters were set to 500r/min for agitation, 2 L/min for aeration, and a pH of 7. Glucose was added when the total sugar content decreased below 50 g/L. Figure 8 illustrates the dynamic changes in various fermentation parameters. After 24 h of fermentation, the cell growth (OD600) reached 3.195, the total sugar residue was 24.09 g/L, and the γ-PGA yield reached 36.29 g/L. However, glutamate content in the fermentation broth remained relatively stable. At this point, an additional 100 g of glucose was introduced, increasing the total sugar content to 55.82 g/L. By the 32-hour point, the glutamate in the fermentation broth decreased to 33.72 g/L, indicating a consumption of 16.27 g/L. Bacterial cell growth increased to 7.68, prompting the second addition of 100 g glucose. At the 48-hour point, bacterial growth (OD600) peaked at 10.82, while the total sugar residue reached 43.67 g/L. To sustain fermentation, 50 g of glucose was added. During the fermentation period from 54 to 96 h, as bacterial cells entered the decay phase and nutrient consumption decreased, the last feeding occurred at 54 h, with 50 g of glucose added. At 82 h of fermentation, the total glucose consumption reached 84.18 g/L, and γ-PGA production peaked at 92.5 g/L, subsequently stabilizing. Towards the end of fermentation, bacterial growth (OD600) decreased to approximately 8.0, with total sugar and glutamate levels in the tank measuring 18.3 g/L and 31.6 g/L, respectively.