Metabolome- and genome-scale model analyses for engineering of Aureobasidium pullulans to enhance polymalic acid and malic acid production from sugarcane molasses
Biotechnology for Biofuels volume 11, Article number: 94 (2018)
Polymalic acid (PMA) is a water-soluble biopolymer with many attractive properties for food and pharmaceutical applications mainly produced by the yeast-like fungus Aureobasidium pullulans. Acid hydrolysis of PMA, resulting in release of the monomer l-malic acid (MA), which is widely used in the food and chemical industry, is a competitive process for producing bio-based platform chemicals.
In this study, the production of PMA and MA from sucrose and sugarcane molasses by A. pullulans was studied in shake flasks and bioreactors. Comparative metabolome analysis of sucrose- and glucose-based fermentation identified 81 intracellular metabolites and demonstrated that pyruvate from the glycolysis pathway may be a key metabolite affecting PMA synthesis. In silico simulation of a genome-scale metabolic model (iZX637) further verified that pyruvate carboxylase (pyc) via the reductive tricarboxylic acid cycle strengthened carbon flux for PMA synthesis. Therefore, an engineered strain, FJ-PYC, was constructed by overexpressing the pyc gene, which increased the PMA titer by 15.1% compared with that from the wild-type strain in a 5-L stirred-tank fermentor. Sugarcane molasses can be used as an economical substrate without any pretreatment or nutrient supplementation. Using fed-batch fermentation of FJ-PYC, we obtained the highest PMA titers (81.5, 94.2 g/L of MA after hydrolysis) in 140 h with a corresponding MA yield of 0.62 g/g and productivity of 0.67 g/L h.
We showed that integrated metabolome- and genome-scale model analyses were an effective approach for engineering the metabolic node for PMA synthesis, and also developed an economical and green process for PMA and MA production from renewable biomass feedstocks.
Polymalic acid (PMA) is a water-soluble biopolymer composed of l-malic acid (MA) monomers and is mainly produced by the yeast-like fungus Aureobasidium pullulans . PMA has free carboxyl groups, making it simple to perform chemical modifications and create various derivatives or carrier-linked pro-drugs. Due to its unique properties, including high water solubility, biocompatibility, and biodegradability, PMA has attracted an increasing attention as a drug carrier or biomaterial in the past few years and is expected to have applications in the preparation of various polymeric micelles, microparticles, nanoconjugates, and nanoparticles for drug delivery systems [2,3,4,5]. PMA also has the potential to be used as a nano-imaging agent for safe and noninvasive diagnosis in the clinical setting . In addition, its monomer MA can be easily generated from PMA via acid hydrolysis. Like other dicarboxylic acids, MA is also an important organic acid and is regarded as a C4 platform chemical in the food and pharmaceutical industries, for which there is a strong global market of over 600,000 tons/year, with an annual growth rate of 4% [7, 8].
Inexpensive renewable feedstocks and media for fermentation, as well as high PMA titers with good productivity are also important constraints in developing industrial PMA fermentation processes. The previous studies have demonstrated fermentative production of PMA from low-cost substrates, such as sweet potato , corncob [10, 11], soy molasses , and wheat straw . However, starchy and lignocellulosic feedstocks must first be subjected to thermochemical, dilute acid, or enzymatic treatments to release fermentable sugars; these limitations must be overcome to achieve an economically competitive process, including enhanced strain tolerance and fermentation performance in the presence of inhibitors released from dilute acid-pretreated lignocellulose .
In this study, we aimed to develop an economical and environmentally friendly bioprocess for PMA and MA production from sugarcane molasses as a feedstock. Sugarcane molasses is the main by-product of sugar production and contains ~ 50% reducing sugars, mainly sucrose and some glucose and fructose; in addition, sugarcane molasses can be used directly in fermentation without requiring any pretreatment [14, 15], and its price is ~ 40% lower than corn (starch and dextrose) based on the fermentable sugar content according to the current market prices. Therefore, sugarcane molasses is considered a relatively inexpensive and renewable biomass feedstocks for biorefinery applications due to its rich sugar content and cost-effectiveness. As omics analysis developing rapidly over the last decade, metabolome analysis has been employed to characterize intracellular metabolic states for understanding cell metabolism and improving the production of target metabolites . Metabolomics has been a powerful tool for understanding the intracellular metabolism with broad range of applications in various fields, including medical science , synthetic biology , medicine , plant systems , and microbial systems . In this study, we evaluated the performance of the strain Aureobasidium pullulans for the utilization of different sugars, including sucrose, fructose, and glucose contained in sugarcane molasses, and compared the metabolic processes of sucrose- and glucose-based fermentation through metabolome analysis. In silico simulation of a genome-scale model of A. pullulans was further verified, and the metabolic node containing pyruvate carboxylase via the reductive tricarboxylic acid (TCA) cycle was engineered to enhance PMA production. Finally, direct sugarcane molasses fermentation, without any pretreatment or nutrient supplementation, was developed for PMA production. This work will be beneficial for the development of an economical and green process for PMA and bio-based MA production from renewable biomass feedstocks.
Strains, media, and culture conditions
The strain A. pullulans CCTCC M2012223 was isolated by our laboratory and can be obtained from the China Center for Type Culture Collection (Wuhan, China). This strain was maintained on potato dextrose agar slants. The seed culture medium contained 60-g/L glucose, 2-g/L NH4NO3, 0.1-g/L KH2PO4, 0.1-g/L MgSO4, 0.1-g/L ZnSO4, 0.5-g/L KCl, and 20-g/L CaCO3. The seed culture was grown in a 500-mL shake flask containing 50-mL liquid medium and incubated at 25 °C on a rotary shaker (180 rpm) for 2 days. The fermentation medium contained 90-g/L sugar (glucose, sucrose, xylose, arabinose, or fructose), 2-g/L NH4NO3, 0.1-g/L KH2PO4, 0.1-g/L MgSO4, 0.1-g/L ZnSO4, 0.5-g/L KCl, and 20-g/L CaCO3. Escherichia coli DH5α cells were employed for routine DNA manipulations. E. coli was grown in LB medium (5 g/L yeast extract, 10 g/L tryptone, 10 g/L NaCl, pH 7.0) at 37 °C, and kanamycin or ampicillin (50 mg/L) was added when required. Agrobacterium tumefaciens AGL1 was grown on YEB medium (10-g/L yeast extract, 5-g/L tryptone, 5-g/L sucrose, 0.5-g/L MgSO4·7H2O, pH 7.0) at 28 °C, and kanamycin or carbenicillin (50 mg/L) was added when required.
Sampling, quenching, and extraction of intracellular metabolites
The cells of strain A. pullulans CCTCC M2012223 were collected after culturing for 48 and 72 h with glucose- and sucrose-based fermentation, respectively, and intracellular metabolism was quenched by immediately adding five volumes of prechilled 60% (v/v) methanol. After quenching at − 40 °C, cells were pelleted in a centrifuge (4000×g, 4 °C, 10 min). The pellets were washed with phosphate-buffered saline (pH 7.4), frozen in liquid nitrogen, and stored at − 80 °C. Next, 60 mg per sample with six biological replicates was ball-milled into a fine powder under frozen conditions with 1 mL of extraction mix consisting of deionized water (50%), methanol (50%), and the internal standards nonadecanoic acid (0.2 mg/mL) and isotope alanine (10 mM), and then lyophilized. The samples were derivatized by vortexing with 60 μL of 15-mg/mL methoxypyridine hydrochloride for 30 s and then incubating for 2 h at 37 °C. The samples were then trimethylsilylated by adding 60-μL bis(trimethylsilyl) trifluoroacetamide (BSTFA) with 1% trimethylchlorosilane and incubating for 90 min at 37 °C. Subsequently, an extract containing all intracellular metabolites was collected via centrifugation at 12,000×g for 10 min at 4 °C .
Metabolite detection with gas chromatography–mass spectrometry (GC–MS)
GC–MS analysis was performed to detect the metabolites in the samples. Briefly, 1 μL sample was injected into an Agilent 7890A/5975C GC–MS system (Agilent Technologies, CA, USA) with a 20:1 split injection ratio. The system was equipped with an HP-5MS capillary column (5% phenyl methyl silox: 30 m × 250 μm i.d., 0.25 μm; Agilent J&W Scientific, Folsom, CA, USA). The inlet temperature was maintained at 280 °C, the interface temperature was 150 °C, and the ion source temperature was 250 °C. The GC temperature program was as follows: hold at 70 °C for 2 min, ramp at 10 °C/min to 300 °C, and a final hold at 300 °C for 5 min. Helium was used as carrier gas at a constant flow rate of 1 mL/min, with a total run time of 30 min. MS conditions were as follows: electron ionization (EI); full-scan mode (m/z 35–780); electron energy, 70 eV . The chromatographic differences between sample groups were obvious, and the retention time was reproducible and stable, demonstrating the reliability of metabolomics analysis.
Data processing and statistical analysis
XCMS software was used to identify and quantify mass spectral peaks. Data preprocessing, including raw data filtering, peak detection, alignment, normalization, and identification, was performed automatically. Metabolites were identified by comparison of mass spectra with the National Institution of Standards and Technology (NIST), Wiley metabolome, Golm metabolome (GM), and in-house databases (Additional file 1: Table S1, red color) . Six independent experiments were performed for each conclusion. Pattern recognition methods based on principal component analysis (PCA) and orthogonal-partial least squares-discriminant analysis (OPLS-DA) were then carried out using SIMCA-P V13.0 (Umetrics AB, Umea, Sweden). All variables were unit variance scaled before PCA and OPLS-DA. Hierarchical clustering analysis was carried out using R software (version 3.1.3) to visualize and group metabolite profiles. A list of variance metabolites that contributed mostly to the model was generated from the P values obtained using Student’s t tests (P < 0.05) and variable importance in the projection values (> 1). In addition, the false discovery rate (FDR) significance criterion was performed with an FDR limit of 0.05 to avoid false-positive results in analyzing differential metabolite correlations. Metabolic pathways were constructed according to pathway analysis of the Metaboanalyst and KEGG metabolic databases [25, 26].
In silico analysis of a genome-scale metabolic model
The reconstructed genome-scale metabolic network (iZX637) of A. pullulans, containing 1347 reactions and 1133 metabolites, was used for constraint-based flux analyses  to determine the most effective route for maximizing PMA production and minimizing the influence of biomass accumulation. The synthesis of PMA is closely related to MA synthesis in A. pullulans. Therefore, three routes of MA synthesis, the oxidative branches of the TCA cycle, the reductive branches of the TCA cycle, and the glyoxylate shunt, were examined. In silico, constrained reactions of the oxidative branches of the TCA cycle contained oxoglutarate dehydrogenase, aconitate hydratase, fumarase, succinate dehydrogenase, and succinate-CoA ligase. Constrained reactions of the reductive branches of the TCA cycle contained malate dehydrogenase, and constrained reactions of the glyoxylate shunt contained fumarase, isocitrate lyase, and glyoxylate lyase. For in silico simulation, the reaction fluxes corresponding to one route of MA synthesis (the oxidative branches of the TCA cycle, the reductive branches of the TCA cycle, or the glyoxylate shunt) were constrained from 0 to 140 mmol/gDW/h (the reaction fluxes of other two routes of MA synthesis were set to 0 gDW/h), and the cell growth rate was constrained from 0 to 1 gDW/h. The specific PMA synthesis rate was then maximized as an objective function and the substrate uptake rates for glucose in the simulation were constrained by 10 mmol/gDW/h. The specific PMA synthesis rate at different flux values of the cell growth rate and route of MA synthesis were plotted on the Z-axis of the flux profile graph. In addition, because the reductive branches of the TCA cycle were assumed to be located only in the cytoplasm in silico, the flux of malate dehydrogenase located in the mitochondria was set to 0 mmol/gDW/h in simulation of three routes. The basic tools used for the model analysis were flux balance analysis (FBA) . GLPK (freeware) (http://www.gnu.org/software/glpk/) was used as a linear programming solver. The COBRA Toolbox 2.0  was used for calculation and analysis in the MATLAB environment.
Gene cloning and expression in the host strain
The DNA of A. pullulans was extracted using a TIANamp Yeast DNA Kit (TIANGEN, China) and was used as a template to amplify the pyruvate carboxylase (pyc, EC: 22.214.171.124) gene. First, one 3.748-kb pyc fragment was amplified using SmaI-pyc-S and MunI-pyc-A, as listed in Additional file 1: Table S1. Second, the pyc fragment obtained was double digested with SmaI and MunI and then ligated into plasmid pBARGPE1, which was also treated with SmaI and EcoRI (MunI isocaudarner) to generate the plasmid pBARGPE1-pyc, which contained a gpdA promoter (PgpdA) and a trpC terminator (TtrpC) flanking the pyc gene. The PgpdA-pyc-TtrpC cassette was amplified by primers PgpdA.pyc.TtrpC-S and PgpdA.pyc.TtrpC-A. The PgpdA-pyc-TtrpC cassette was cloned into the EcoRI site of the pK2-PtrpC-hyg-TtrpC plasmid using a One Step Cloning Kit (Vazyme, USA), forming the plasmid pK2-hyg-pyc for A. tumefaciens-mediated transformation (ATMT), as described in our previous study . T-DNA insertion sites were analyzed using a Genome Walking Kit (TaKaRa, Japan) in ATMT-derived clones. Polymerase chain reaction (PCR) primers were used to amplify the T-DNA flanking sequence (Additional file 1: Table S1).
Batch and fed-batch fermentation in a 5-L stirred-tank fermentor
Batch fermentation kinetics were studied in a 5-L stirred-tank fermentor (Shanghai Baoxing Co., Ltd., China) containing 3.5 L of fermentation media (with glucose as carbon resource) or only diluted sugarcane molasses. The fermentation medium was inoculated with 350 mL of the seed culture grown in a shake flask for 48 h and operated at 25 °C with agitation and aeration at 400–700 rpm and 0.8 vvm, respectively. During fermentation, the agitation was controlled to maintain the dissolved oxygen level at greater than 20%. For sugarcane molasses fermentation, the crude sugarcane molasses (containing ~ 420 g/kg of total sugar and 14.29 ± 0.27 g/kg of total nitrogen) were diluted by water to ~ 90 g/L total sugar. For fed-batch fermentation, when the total sugar concentration was lower than 30 g/L, sugarcane molasses (approximately 270-g/L reducing sugars) was continuously added to maintain the total sugar concentration at approximately 10–20 g/L. Broth samples were collected periodically for the analysis of residual sugar, biomass, and PMA titers. All trials were performed in triplicate.
Quantitative reverse transcription (RT)-PCR
Total RNA was extracted using a Fungal RNA Kit (Omega, USA) and reverse transcribed into cDNA using reverse transcriptase (TIANGEN). The quantitative PCR assay was performed using the SYBR Green method (FastQuant RT Super Mix; TIANGEN). Primers for the pyc gene are shown in Additional file 1: Table S1. The gene encoding β-actin was used as the reference gene.
Biomass was determined by the dry cell weight method. Prior to measurement, excessive CaCO3 was eliminated from the broth by adding 1-M HCl. The cell suspension was centrifuged at 4000×g and then dried overnight at 105 °C . PMA was analyzed by centrifuging the fermentation broth and the resulting supernatant was mixed with an equal volume of 2-M H2SO4 in an incubator at 85 °C for 8 h. The hydrolyzed PMA sample was analyzed by high-performance liquid chromatography (Agilent 1260, USA) using a Spursil C18-EP organic acid column at 40 °C and eluted with 5-mM H2SO4 at a rate of 0.6 mL/min to determine its malic acid content . The glucose concentration was measured using a fully automatic residual sugar analyzer (Biology Institute of Shandong Academy of Sciences, China).
Results and discussion
Evaluation of carbon sources in shake-flask culture
Sugarcane molasses is the main by-product of sugar production, and mainly comprises sucrose and some glucose and fructose. In this study, fermentation with different carbon sources in shake-flask culture was investigated (Fig. 1). Of the carbon sources used, sucrose and fructose showed superior PMA production of 28.30 ± 2.58 and 22.58 ± 1.08 g/L, respectively, which was higher than that (20.38 ± 2.32 g/L) from glucose as the carbon source. Moreover, five-carbon sugar (arabinose and xylose) fermentation was also superior to that with glucose for PMA production. In our previous study, a novel PMA-producing strain, A. pullulans YJ6-11, was also found to utilize xylose as a superior carbon source compared with glucose . The fermentation kinetics of sucrose and glucose in shake-flask culture are shown in Fig. 2. A total of 33.91 ± 3.70 g/L PMA (~ 38.98 g/L MA after hydrolysis) were produced from ~ 90 g/L sucrose in 96 h via batch fermentation; this was increased by 40.7% compared with that of glucose (24.10 ± 3.30 g/L). These results showed that sucrose present in sugarcane molasses was the most suitable sugar for PMA biosynthesis. However, the underlying factors responsible for the differences between sucrose and glucose were still unclear.
Analysis of key metabolites and pathways associated with PMA biosynthesis
To explain the differences in metabolites of A. pullulans, the metabolic states of sucrose- and glucose-based fermentation were compared by metabolomics analysis. As shown in Additional file 2: Table S2, in total, 81 intracellular metabolites were identified and divided into eight categories. These groups mainly included organic acids, amino acids, and sugar metabolism that contained 21, 15, and 14 metabolites, respectively. Among these metabolites, compared with culture at 48 h, 13 and 3 differential metabolites were increased at 72 h in glucose and sucrose, respectively. 27 and 30 differentially abundant metabolites from glucose versus sucrose after culture for 48 and 72 h, respectively, were validated. Of these metabolites, both at 48 and 72 h, 16 in sucrose fermentation, including fructose, ribose, glycolic acid, oleic acid, and heptanoic acid, were higher than that of glucose (Additional files 3 and 4: Figure S1 and Table S3). These results indicated that the differences of metabolites were varied in the different sugar fermentation. Moreover, compared with glucose fermentation at 48 h, 13 differential metabolites were increased at 72 h. Notably, at 72 h, citric acid was increased by 1.84-fold, whereas other intermediate metabolites associated with the TCA cycle (such as pyruvic acid, succinic acid, and fumaric acid) did not differ significantly at 48 h (Additional file 4: Table S3). However, at 72 h, MA was decreased by 1.04-fold compared with that at 48 h, which was negative correlated with citric acid (Fig. 3). This may induce a decrease in PMA biosynthesis, because accumulation of citric acid resulted in carbon flux that was not efficiently channeled to the objective product. Pyruvate is a key node in central carbon metabolism and plays a major role in pathways related to organic acid synthesis . The carbon flux can be channeled into different products through pyruvate, including succinate , lactate , and malate . The reductive TCA cycle has been identified as the most efficient pathway for many products, such as succinate  and malate . Thus, the accumulation of citric acid also indicated that metabolic flux could not effectively direct pyruvate to PMA synthesis or the high-efficiency PMA synthesis route in native metabolic regulation. Therefore, the enhanced carboxylation of pyruvate to oxaloacetate, catalyzed by pyruvate carboxylase, may be an important target for improving PMA synthesis.
In addition, pyruvate is also a major component in pathways related to amino acid metabolism. Amino acids are crucial to microorganism metabolism and have roles in synthesis of proteins (e.g., leucine, threonine, and phenylalanine), regulation of gene expression (e.g., arginine and leucine), osmoregulation (e.g., citrulline and proline), control of enzyme activity (e.g., alanine and leucine), transamination (e.g., aspartate, alanine, and glutamate), and ammonia detoxification (e.g., citrulline, arginine, and glutamate) . As shown in Fig. 3, amino acids present at higher levels, particularly leucine and proline, could be directly or indirectly transformed from pyruvate and α-oxoglutarate . Most amino acids were observed at higher levels in sucrose than in glucose because of the abundant precursors in the glycolytic pathway and TCA cycle. The redox balance is maintained by modifying the distribution of the metabolic flux of pyruvate .
In addition, when sucrose was used as the carbon source, relative concentration of trehalose was significantly decreased 2.37- and 2.70-fold at 48 and 72 h, respectively, compared with that when glucose was used as the carbon source (Additional file 4: Table S3). Accumulation of trehalose in yeasts has been shown to be an important mechanism mediating tolerance against adverse stress conditions, and trehalose functions as a reserve carbohydrate that can be synthesized when the exogenous energy exceeds the cellular needs for growth and biosynthesis [42, 43]. These results suggested that glucose as a rapidly metabolizable carbon source can induce a regulatory response (e.g., osmotic stress) to restore excess exogenous energy, which would cause some of the carbon source to be used for other processes, not PMA synthesis, thereby decreasing PMA titers.
In silico analysis of a genome-scale metabolic model
PMA is synthesized from MA, which is an intermediate in the TCA cycle in aerobic metabolism. Genome-scale flux sensitivity analyses were performed to determine the most effective route for maximizing PMA synthesis. Three metabolic routes involved in MA synthesis, including oxidative branches of the TCA cycle, reductive branches of the TCA cycle, and glyoxylate shunt, were simulated using the genome-scale metabolic model (iZX637). As shown in Fig. 4, the optimal PMA synthesis rate, using glucose or sucrose as the carbon source, reached 15 mmol/gDW/h (indicated by an arrow), increasing the reductive TCA malate synthesis rate to 15 mmol/gDW/h and the biomass growth rate to 0 mmol/gDW/h. When the carbon flux turned to oxidative branches of the TCA cycle and the glyoxylate shunt, the optimal PMA synthesis rates reached only as high as 12.89 mmol/gDW/h (sucrose as a carbon source: 12.58 mmol/gDW/h) and 14.67 mmol/gDW/h (sucrose as a carbon source: 14.33 mmol/gDW/h) with corresponding malate synthesis rates of 13 and 8 mmol/gD/h, respectively. In addition, the malate synthesis rate of the corresponding route constantly increased, resulting in a negative influence on the optimal PMA synthesis rate. When the malate synthesis rates of the reductive TCA route, oxidative TCA route, and glyoxylate shunt reached 120, 37, and 37 mmol/gDW/h, respectively, the PMA synthesis rate and biomass growth rate were reduced to 0 mmol/gDW/h (glucose or sucrose as a carbon source). Furthermore, the PMA synthesis rate decreased gradually as the biomass synthesis rate increased for all simulation results. These results showed that PMA and biomass synthesis showed a competitive relationship. However, in the reductive TCA route, biomass synthesis could be maintained with a relatively high PMA synthesis rate compared with the other two routes, regardless of whether glucose or sucrose was used as the carbon source. In our previous study, we found that high levels of the nitrogen source favored cell growth but decreased PMA synthesis . Therefore, the efficiency of the PMA synthesis route should balance cell growth and PMA formation. In our previous study, in silico analysis of the carbon flux distribution and changes in PMA synthesis rates showed that, in the context of a high PMA synthesis rate, a large amount of carbon flux was transformed into pyruvate and channeled into the reductive TCA route by pyruvate carboxylase . These in silico results showed that the reductive TCA route, as the optimal PMA synthesis route, showed a relatively high optimal PMA synthesis rate and mitigated the competitive relationship between PMA and biomass synthesis. In addition, the simulation results also showed that the reductive TCA route required a relatively high route flux (the reductive TCA malate synthesis rate to 15 mmol/gDW/h), yielding a high PMA synthesis rate (the optimal PMA synthesis rate to 15 mmol/gDW/h). Therefore, pyruvate carboxylase (pyc) in the reductive TCA route was selected as the target to be overexpressed to improve PMA yield and productivity.
Moreover, according to the simulation result from sucrose and glucose as carbon sources, sucrose was not dominant compared with glucose as a carbon source for the optimal PMA synthesis rate and cell growth by A. pullulans strains in the conversion efficiency of metabolic pathways. This could be explained by the observations that the reaction of sucrose hydrolysis exhibited no energy consumption or carbon loss in silico.
Overexpression of the pyc gene in A. pullulans
To evaluate the effects of the pyc gene on PMA biosynthesis, the endogenous pyc gene was amplified using the genomic DNA of A. pullulans as a template (Additional file 3: Figure S2). The PgpdA promoter of Aspergillus nidulans, which directs the constitutively high expression of the glyceraldehyde-3-phosphate dehydrogenase (gpd) gene, can be used to induce high levels of constitutive expression of genes of interest or marker genes and to drive the overexpression of various genes in fungi [45,46,47]. Therefore, the PgpdA promoter of A. nidulans was used to constitute the pyc cassette by cloning into pBARGPE1 (Additional file 3: Figure S3). In our previous study, we developed a simple and efficient system for genetic transformation of A. pullulans using A. tumefaciens . Therefore, A. tumefaciens carrying the pK2-hyg-pyc binary plasmid was used to transform A. pullulans. Because ATMT-mediated integration is random, five different ATMT-derived clones were tested for PMA production in shake-flask culture. The PMA titers of five clones compared with the wild-type strain were improved from 3.5 to 7.1% (Additional file 3: Figure S4). The differences among mutants were due to random integration, causing interruption of gene function by insertion of T-DNA, and were not caused by intrinsic factors . This result indicated that the increased PMA titer was mainly dependent on pyc gene overexpression. Thus, we selected the highest PMA producer (E10 strain; as called FJ-PYC) for further analysis of the sequences flanking the T-DNA in the genome. As shown in Additional file 3: Figure S5, the hygromycin B resistance cassette and overexpression pyc cassette were inserted into the noncoding regions of g4940.t1 and g4941.t1 on the A. pullulans genome. The insertions did not disrupt the g4940.t1 and g4941.t1 open reading frames. Furthermore, we found that the expression level of the pyc gene in strain E-10 was increased by 9.5-fold compared with that in the wild-type strain (Fig. 5). These results also revealed that the promoter PgpdA was an effective tool for regulating gene expression in A. pullulans.
Subsequently, the engineered strain E10 was tested with glucose as the carbon source in a 5-L stirred-tank fermentor. As shown in Fig. 6, about 36.2 g/L of PMA was produced via batch fermentation, with a PMA productivity of 0.53 g/L h, which was increased by 15.1 and 12.7% compared with the control, respectively. This result indicated that overexpression of the pyc gene could improve the production of PMA, which was consistent with the results of in silico simulation showing that strengthening of the reductive TCA route resulted in increased carbon flux toward PMA synthesis.
Batch and fed-batch fermentation with sugarcane molasses
Sugarcane molasses is a low-cost by-product of sugar production and mainly contains sucrose, with some fructose and glucose. To develop an economical fermentation process, sugarcane molasses without any pretreatment or nutrient supplementation was employed for PMA production by FJ-PYC strain in a 5-L stirred-tank fermentor. As shown in Fig. 7a, during the early fermentation stage, sucrose was rapidly consumed and transformed into glucose and fructose, and glucose was then used for cell growth and PMA synthesis compared with fructose. This phenomenon, called glucose-mediated carbon catabolite repression (CCR), is widely observed in bacteria and yeasts [48,49,50]. About 31.5-g/L PMA (36.5-g/L MA after hydrolysis) was produced from ~ 95-g/L mixed sugar (61.8-g/L sucrose, 18.5-g/L glucose, and 15.3-g/L fructose from sugarcane molasses) in 60 h via batch fermentation, with an MA yield of 0.44 g/g. Moreover, without adding any other media, the productivity of PMA (0.53 g/L h) was the same as that from glucose fermentation. In fed-batch fermentation (Fig. 7b), a high final PMA titer of 81.5 g/L (94.2-g/L MA after hydrolysis) was achieved in 140 h with a corresponding MA yield of 0.62 g/g and productivity of 0.67 g/L h. As shown in Table 1, compared with other renewable feedstocks, PMA fermentation from sugarcane molasses by strain FJ-PYC resulted in the highest titer and productivity. Nevertheless, the titer and productivity could be further increased by continuous fed operation, because the PMA synthesis rate maintained a high speed. In addition, the fed-batch fermentation with FJ-PYC resulted in an overall MA yield of 0.62 g/g from sugarcane molasses without nitrogen supplementation, which was the same as the highest yield from sugarcane juice . However, sugarcane juice contained only ~ 15% reducing sugars and relatively fewer inhibitors, requiring it to be condensed for high-sugar PMA fermentation and dramatically decreasing its economic efficiency. These results indicated that the strain FJ-PYC exhibited extremely high adaptation and tolerance against sugarcane molasses.
In this study, different carbon sources were evaluated, and sugarcane molasses was assessed as a potential feedstock for economic PMA and MA production. Among the different carbon sources examined in this study, sucrose was found to be the optimal carbon source for PMA biosynthesis. Metabolomics analysis of sucrose- and glucose-based fermentation identified a total of 81 intracellular metabolites and showed that pyruvate from the glycolysis pathway may be a key metabolite affecting PMA synthesis. In silico analysis of a genome-scale metabolic model (iZX637) verified that the pyc gene via the reductive TCA cycle was a target affecting PMA synthesis. Consistent with this, overexpression of the pyc gene in A. pullulans strain FJ-PYC increased the PMA titer by 15.1% compared with the control.
Moreover, sugarcane molasses could be directly utilized by this recombinant strain without any pretreatment or nutrient supplementation, producing 31.5-g/L PMA (36.5-g/L MA) with a high-level PMA productivity of 0.53 g/L h in batch fermentation. In fed-batch fermentation, compared with the other renewable feedstock, the highest final PMA titer of 81.5 g/L (94.2-g/L MA after hydrolysis) was achieved in 140 h, with a corresponding MA yield of 0.62 g/g and productivity of 0.67 g/L h. This study demonstrated the great potential of sugarcane molasses for the economical production of PMA and MA on the industrial scale.
- pyc :
pyruvate carboxylase gene
genome scale metabolic model
high-performance liquid chromatography
gas chromatography–mass spectrometry
National Institution of Standards and Technology
principal component analysis
orthogonal-partial least squares-discriminant analysis
false discovery rate
- PgpdA :
- TtrpC :
A. tumefaciens-mediated transformation
polymerase chain reaction
- gpd :
glyceraldehyde-3-phosphate dehydrogenase gene
carbon catabolite repression
Chi Z, Liu GL, Liu CG, Chi ZM. Poly(beta-l-malic acid) (PMLA) from Aureobasidium spp. and its current proceedings. Appl Microbiol Biotechnol. 2016;100:3841–51.
Yang YT, Lee SJ, Nai YS, Kim S, Kim JS. Up-regulation of carbon metabolism-related glyoxylate cycle and toxin production in Beauveria bassiana JEF-007 during infection of bean bug, Riptortus pedestris (Hemiptera: Alydidae). Fungal Biol. 2016;120:1236–48.
Zhou Q, Yang T, Qiao Y, Guo S, Zhu L, Wu H. Preparation of poly(beta-l-malic acid)-based charge-conversional nanoconjugates for tumor-specific uptake and cellular delivery. Int J Nanomed. 2015;10:1941–52.
Lanz-Landazuri A, Portilla-Arias J, Martinez de Ilarduya A, Garcia-Alvarez M, Holler E, Ljubimova J, Munoz-Guerra S. Nanoparticles of esterified polymalic acid for controlled anticancer drug release. Macromol Biosci. 2014;14:1325–36.
Loyer P, Cammas-Marion S. Natural and synthetic poly(malic acid)-based derivates: a family of versatile biopolymers for the design of drug nanocarriers. J Drug Target. 2014;22:556–75.
Patil R, Gangalum PR, Wagner S, Portilla-Arias J, Ding H, Rekechenetskiy A, Konda B, Inoue S, Black KL, Ljubimova JY, Holler E. Curcumin targeted, polymalic acid-based mri contrast agent for the detection of a plaques in alzheimer’s disease. Macromol Biosci. 2015;15:1212–7.
Goldberg I, Rokem JS, Pines O. Organic acids: old metabolites, new themes. Chem Technol Biotechnol. 2006;81:1601–11.
Alonso S, Rendueles M, Diaz M. Microbial production of specialty organic acids from renewable and waste materials. Crit Rev Biotechnol. 2015;35:497–513.
Zan Z, Zou X. Efficient production of polymalic acid from raw sweet potato hydrolysate with immobilized cells of Aureobasidium pullulans CCTCC M2012223 in aerobic fibrous bed bioreactor. J Chem Technol Biotechnol. 2013;88:1822–7.
Wang J, Li L, Niu X, Zou D. Phosphine-induced phosphorus mobilization in the rhizosphere of rice seedlings. J Soils Sediments. 2016;16:1735–44.
Zou X, Wang Y, Tu G, Zan Z, Wu X. Adaptation and transcriptome analysis of Aureobasidium pullulans in corncob hydrolysate for increased inhibitor tolerance to malic acid production. PLoS ONE. 2015;10:e0121416.
Cheng C, Zhou Y, Lin M, Wei P, Yang S-T. Polymalic acid fermentation by Aureobasidium pullulans for malic acid production from soybean hull and soy molasses: fermentation kinetics and economic analysis. Bioresour Technol. 2017;223:166–74.
Leathers TD, Manitchotpisit P. Production of poly(beta-l-malic acid) (PMA) from agricultural biomass substrates by Aureobasidium pullulans. Biotechnol Lett. 2013;35:83–9.
Chan S, Kanchanatawee S, Jantama K. Production of succinic acid from sucrose and sugarcane molasses by metabolically engineered Escherichia coli. Bioresour Technol. 2012;103:329–36.
Jung MY, Park BS, Lee J, Oh MK. Engineered Enterobacter aerogenes for efficient utilization of sugarcane molasses in 2,3-butanediol production. Bioresour Technol. 2013;139:21–7.
Putri SP, Nakayama Y, Matsuda F, Uchikata T, Kobayashi S, Matsubara A, Fukusaki E. Current metabolomics: practical applications. J Biosci Bioeng. 2013;115:579–89.
Teul J, Ruperez FJ, Garcia A, Vaysse J, Balayssac S, Gilard V, Malet-Martino M, Martin-Ventura JL, Blanco-Colio LM, Tunon J, et al. Improving metabolite knowledge in stable atherosclerosis patients by association and correlation of GC–MS and H-1 NMR fingerprints. J Proteome Res. 2009;8:5580–9.
Villas-Boas SG, Akesson M, Nielsen J. Biosynthesis of glyoxylate from glycine in Saccharomyces cerevisiae. FEMS Yeast Res. 2005;5:703–9.
Bugrim A, Nikolskaya T, Nikolsky Y. Early prediction of drug metabolism and toxicity: systems biology approach and modeling. Drug Discov Today. 2004;9:127–35.
Kim JK, Bamba T, Harada K, Fukusaki E, Kobayashi A. Time-course metabolic profiling in Arabidopsis thaliana cell cultures after salt stress treatment. J Exp Bot. 2007;58:415–24.
Mitsunaga H, Meissner L, Palmen T, Bamba T, Buechs J, Fukusaki E. Metabolome analysis reveals the effect of carbon catabolite control on the poly(gamma-glutamic acid) biosynthesis of Bacillus licheniformis ATCC 9945. J Biosci Bioeng. 2016;121:413–9.
Vanholme R, Storme V, Vanholme B, Sundin L, Christensen JH, Goeminne G, Halpin C, Rohde A, Morreel K, Boerjan W. A systems biology view of responses to lignin biosynthesis perturbations in Arabidopsis. Plant Cell. 2012;24:3506–29.
Wu Q, Zhang L, Xia H, Yu C, Dou K, Li Y, Chen J. Omics for understanding synergistic action of validamycin A and Trichoderma asperellum GDFS1009 against maize sheath blight pathogen. Sci Rep. 2017;7:40140.
Smith CA, Want EJ, O’Maille G, Abagyan R, Siuzdak G. XCMS: processing mass spectrometry data for metabolite profiling using Nonlinear peak alignment, matching, and identification. Anal Chem. 2006;78:779–87.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.
Kopka J, Schauer N, Krueger S, Birkemeyer C, Usadel B, Bergmuller E, Dormann P, Weckwerth W, Gibon Y, Stitt M, et al. GMD@CSB.DB: the Golm metabolome database. Bioinformatics. 2005;21:1635–8.
Feng J, Yang J, Li X, Guo M, Wang B, Yang ST, Zou X. Reconstruction of a genome-scale metabolic model and in silico analysis of the polymalic acid producer Aureobasidium pullulans CCTCC M2012223. Gene. 2017;607:1–8.
Orth JD, Thiele I, Palsson BO. What is flux balance analysis? Nat Biotechnol. 2010;28:245–8.
Schellenberger J, Que R, Fleming RM, Thiele I, Orth JD, Feist AM, Zielinski DC, Bordbar A, Lewis NE, Rahmanian S, et al. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA toolbox v2.0. Nat Protoc. 2011;6:1290–307.
Tu G, Wang Y, Feng J, Li X, Guo M, Zou X. Agrobacterium tumefaciens-mediated transformation of Aureobasidium pullulans and high-efficient screening for polymalic acid producing strain. Chin J Biotechnol. 2015;31:1063–72.
Zou X, Zhou YP, Yang ST. Production of polymalic acid and malic acid by Aureobasidium pullulans fermentation and acid hydrolysis. Biotechnol Bioeng. 2013;110:2105–13.
Zou X, Tu GW, Zan ZQ. Cofactor and CO2 donor regulation involved in reductive routes for polymalic acid production by Aureobasidium pullulans CCTCC M2012223. Bioprocess Biosyst Eng. 2014;37:2131–6.
Zou X, Yang J, Tian X, Guo M, Li Z, Li Y. Production of polymalic acid and malic acid from xylose and corncob hydrolysate by a novel Aureobasidium pullulans YJ 6–11 strain. Process Biochem. 2016;51:16–23.
Yin X, Li J, Shin H-d DuG, Liu L, Chen J. Metabolic engineering in the biotechnological production of organic acids in the tricarboxylic acid cycle of microorganisms: advances and prospects. Biotechnol Adv. 2015;33:830–41.
Gokarn RR, Evans JD, Walker JR, Martin SA, Eitemann MA, Altman E. The physiological effects and metabolic alterations caused by the expression of Rhizobium etli pyruvate carboxylase in Escherichia coli. Appl Microbiol Biotechnol. 2001;56:188–95.
Su Y, Rhee MS, Ingram LO, Shanmugam KT. Physiological and fermentation properties of Bacillus coagulans and a mutant lacking fermentative lactate dehydrogenase activity. J Ind Microbiol Biotechnol. 2011;38:441–50.
Zelle RM, de Hulster E, van Winden WA, de Waard P, Dijkema C, Winkler AA, Geertman JMA, van Dijken JP, Pronk JT, van Maris AJA. Malic acid production by Saccharomyces cerevisiae: engineering of pyruvate carboxylation, oxaloacetate reduction, and malate export. Appl Environ Microbiol. 2008;74:2766–77.
Sanchez AM, Bennett GN, San KY. Efficient succinic acid production from glucose through overexpression of pyruvate carboxylase in an Escherichia coli alcohol dehydrogenase and lactate dehydrogenase mutant. Biotechnol Prog. 2005;21:358–65.
Wu G. Amino acids: metabolism, functions, and nutrition. Amino Acids. 2009;37:1–17.
Xia M, Huang D, Li S, Wen J, Jia X, Chen Y. Enhanced FK506 production in Streptomyces tsukubaensis by rational feeding strategies based on comparative metabolic profiling analysis. Biotechnol Bioeng. 2013;110:2717–30.
Thierry A, Deutsch S-M, Falentin H, Dalmasso M, Cousin FJ, Jan G. New insights into physiology and metabolism of Propionibacterium freudenreichii. Int J Food Microbiol. 2011;149:19–27.
Jiang H, Liu G-L, Chi Z, Hu Z, Chi Z-M. Genetics of trehalose biosynthesis in desert-derived Aureobasidium melanogenum and role of trehalose in the adaptation of the yeast to extreme environments. Curr Genet. 2018;64:479–91.
Tapia H, Young L, Fox D, Bertozzi CR, Koshland D. Increasing intracellular trehalose is sufficient to confer desiccation tolerance to Saccharomyces cerevisiae. PNAS. 2015;112:6122–7.
Wang Y, Song X, Zhang Y, Wang B, Zou X. Effects of nitrogen availability on polymalic acid biosynthesis in the yeast-like fungus Aureobasidium pullulans. Microb Cell Fact. 2016;15:1–12.
Cao Y, Peng G, He Z, Wang Z, Yin Y, Xia Y. Transformation of Metarhizium anisopliae with benomyl resistance and green fluorescent protein genes provides a tag for genetically engineered strains. Biotechnol Lett. 2007;29:907–11.
Jin K, Zhang Y, Fang W, Luo Z, Zhou Y, Pei Y. Carboxylate transporter gene JEN1 from the entomopathogenic fungus Beauveria bassiana is involved in conidiation and virulence. Appl Environ Microbiol. 2010;76:254–63.
Zhu ZW, Zhang SF, Liu HW, Shen HW, Lin XP, Yang F, Zhou YJJ, Jin GJ, Ye ML, Zou HF, Zhao ZBK. A multi-omic map of the lipid-producing yeast Rhodosporidium toruloides. Nat Commun. 2012;3:11.
Wei P, Lin M, Wang Z, Fu H, Yang H, Jiang W, Yang ST. Metabolic engineering of Propionibacterium freudenreichii subsp. shermanii for xylose fermentation. Bioresour Technol. 2016;219:91–7.
Wei PL, Cheng C, Lin M, Zhou YP, Yang ST. Production of poly(malic acid) from sugarcane juice in fermentation by Aureobasidium pullulans: kinetics and process economics. Bioresour Technol. 2017;224:581–9.
Yao R, Shimizu K. Recent progress in metabolic engineering for the production of biofuels and biochemicals from renewable sources with particular emphasis on catabolite regulation and its modulation. Process Biochem. 2013;48:1409–17.
FJ and ZX designed and carried out the model simulations and fermentations. YJ and JC evaluated effects of carbon sources on cell growth and PMA production different carbon sources and analyzed the data of metabolome. FJ and YJ drafted the manuscript. YWW assisted in shake-flask fermentation. ZX and MJ supervised the project and revised the manuscript. All authors read and approved the final manuscript.
This study was supported in part by Grants from the National Natural Science Foundation of China (Grant No. 31571816), National High Technology Research and Development Program of China (863 Program) (2015AA021005), Chongqing Social and People’s Livelihood Guarantee Special Program (cstc2016shmszx80075), and State Key Laboratory of Materials-Oriented Chemical Engineering (Nanjing Tech University), China (KL15-10).
The authors declare that they have no competing interests.
Availability of supporting data
The GEM iZX637 during the current study is available in the PubMed repository https://doi.org/10.1016/j.gene.2016.12.034.The authors declare that all other data supporting the findings of this study are available within the article its Additional files.
Consent for publication
The authors consent for publication.
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
List of strains, plasmids, and primers.
Identified metabolites by comparison of mass spectra.
The statistical numbers of differential metabolites between groups. Red bars represent the numbers of relative concentration increased metabolites and blue bars represent the numbers of relative concentration decreased metabolites. Figure S2. PCR amplification of the endogenous pyc gene using the genomic DNA of A. pullulans. Figure S3. Construction of the pyc cassette by cloning into the plasmid pBARGPE1 with the PgpdA promoter. Figure S4. Batch fermentation of different ATMT-derived clones with glucose as carbon source in shake flask. Figure S5. Analysis of sequences flanking to T-DNA in genome of the strain E10.
The changes of differential metabolites in glucose- and sucrose-based fermentation.
About this article
Cite this article
Feng, J., Yang, J., Yang, W. et al. Metabolome- and genome-scale model analyses for engineering of Aureobasidium pullulans to enhance polymalic acid and malic acid production from sugarcane molasses. Biotechnol Biofuels 11, 94 (2018). https://doi.org/10.1186/s13068-018-1099-7
- Aureobasidium pullulans
- Metabolic engineering
- Polymalic acid
- Pyruvate carboxylase
- Sugarcane molasses