Metabolic engineering of Synechococcus elongatus PCC 7942 for improvement of 1,3-propanediol and glycerol production based on in silico simulation of metabolic flux distribution
© The Author(s) 2017
Received: 11 October 2017
Accepted: 13 November 2017
Published: 25 November 2017
Production directly from carbon dioxide by engineered cyanobacteria is one of the promising technologies for sustainable future. Previously, we have successfully achieved 1,3-propanediol (1,3-PDO) production using Synechococcus elongatus PCC 7942 with a synthetic metabolic pathway. The strain into which the synthetic metabolic pathway was introduced produced 3.48 mM (0.265 g/L) 1,3-PDO and 14.3 mM (1.32 g/L) glycerol during 20 days of incubation. In this study, the productivities of 1,3-PDO were improved by gene disruption selected by screening with in silico simulation.
First, a stoichiometric metabolic model was applied to prediction of cellular metabolic flux distribution in a 1,3-PDO-producing strain of S. elongatus PCC 7942. A genome-scale model of S. elongatus PCC 7942 constructed by Knoop was modified by the addition of a synthetic metabolic pathway for 1,3-PDO production. Next, the metabolic flux distribution predicted by metabolic flux balance analysis (FBA) was used for in silico simulation of gene disruption. As a result of gene disruption simulation, NADPH dehydrogenase 1 (NDH-1) complexes were found by screening to be the most promising candidates for disruption to improve 1,3-PDO production. The effect of disruption of the gene encoding a subunit of the NDH-1 complex was evaluated in the 1,3-PDO-producing strain.
Results and Conclusions
During 20 days of incubation, the ndhF1-null 1,3-PDO-producing strain showed the highest titers: 4.44 mM (0.338 g/L) 1,3-PDO and 30.3 mM (2.79 g/L) glycerol. In this study, we successfully improved 1,3-PDO productivity on the basis of in silico simulation of gene disruption.
In silico simulation is an powerful methodology to improve productivity with low experimental costs. The stoichiometric metabolic model not requiring kinetic parameters of metabolic reactions is comparatively applicable to analysis of whole-cell metabolism. Flux balance analysis (FBA) can predict metabolic flux distribution under constraint conditions . FBA using the genome-scale metabolic model (GSM) has been applied to improvement, for example, of succinic acid , threonine , valine , and lycopene  production by Escherichia coli. For screening of effective candidates for gene manipulation, various strategies have been established for prediction of the effects of a gene knockout, knockdown, and overexpression . Shastri and Morgan  were the first to construct a stoichiometric metabolic model of cyanobacteria composed of the central metabolism and photosystem of Synechocystis sp. PCC 6803. So far, the metabolic models of some cyanobacteria used for bioproduction (Synechocystis sp. PCC 6803, Synechococcus sp. PCC 7002, Synechococcus elongatus PCC 7942, Arthrospira platensis, Cyanothece sp. ATCC 51142) have been constructed on the basis of an annotated genome sequence [30–32]. Anfelt et al.  reported that introduction of the phosphoketolase pathway increases cellular acetyl-CoA level and butanol productivity in Synechocystis sp. PCC 6803. FBA was used for validation of this genetic modification to increase theoretical productivity of butanol. Furthermore, some simulations based on FBA have predicted various metabolic modifications (knockout and overexpression) to increase theoretical production by cyanobacteria, for example, of biofuels (1-butanol, 1-octanol, limonene, ethanol) and succinic acid by Synechocystis sp. PCC 6803 [33–35] and biofuels (butanol, ethanol) by Synechococcus sp. PCC 7002 . Nonetheless, there are only a few reports that these predicted strategies for productivity improvement applied experimentally .
In this study, we applied FBA to prediction of metabolic flux distribution in a 1,3-PDO-producing strain of S. elongatus PCC 7942 constructed in our previous work . The candidates for the gene knockout to increase 1,3-PDO productivity were estimated by FBA with the minimization of metabolic adjustment (MOMA) algorithm . The positive effect of the gene knockout on 1,3-PDO productivity after screening by in silico simulation was experimentally validated in the 1,3-PDO-producing strain.
Chemicals and reagents
All the chemicals were purchased from Wako Pure Chemical Industries, Ltd. (Osaka, Japan) unless specified otherwise. Restriction enzymes, phosphatase (New England Biolabs; Ipswich, MA, USA), ligase (Rapid DNA Ligation Kit, Roche; Mannheim, Germany), and DNA polymerase (KOD Plus Neo DNA polymerase, TOYOBO Co., Ltd.; Osaka, Japan) were used for cloning. Oligonucleotides were synthesized by Life Technologies Japan, Ltd. (Tokyo, Japan).
A modified BG11 medium supplemented with 20 mM HEPES–NaOH (pH 7.5) for pH stabilization was used for cultivation of S. elongatus PCC 7942 (Life Technologies Corporation; Carlsbad, CA, USA). Hereafter, this medium is referred to as the BG11 medium. The composition of the BG11 medium used in this study was the same as previously reported . Antibiotics (10 μg/mL kanamycin, 20 µg/mL spectinomycin, and 5 µg/mL chloramphenicol) were added to the BG11 medium when appropriate.
The incubation conditions for 1,3-PDO production were the same as those described previously . Cell density (OD730) was measured using Infinite 200 PRO (TECAN; Männedorf, Switzerland). For determination of dry cell weight (DCW), cells in culture were harvested on the Ominipore™ membrane filter (JHWP04700; Merck Millipore; Darmstadt, Germany) by filtration under reduced pressure. Cell harvest filters were dried at 65 °C for 12 h. The weight differences before and after cell harvesting were determined as DCW in culture.
Plasmids and primers used in this study are listed in Additional file 1: Tables S1 and S2, respectively. To construct a plasmid for gene disruption, the upper and lower regions of a target gene were amplified for homologous recombination by PCR. The primer sets for the upstream region of ndhD1, ndhD2, ndhD3, ndhD4, and ndhF1 were T2711–T2712, T3189–T3190, T2723–T2724, T2719–T2720, and T3208–3209, respectively. The primer sets for the downstream region of ndhD1, ndhD2, ndhD3, ndhD4, and ndhF1 were T3187–T3188, T3191–T3192, T2787–T2726, T2721–2722, and T2713–T2714, respectively. The amplified upstream and downstream regions except for ndhF1 were digested by KpnI–HindIII and HindIII–BamHI, respectively, and ligated into KpnI–BamHI sites of pTA703. Only for ndhF1, the upstream and downstream regions were digested with SalI–HindIII and HindIII–BamHI, respectively, and ligated into SalI–BamHI sites of pTA703. A chloramphenicol resistance cassette was amplified from pZA31-luc (Expressys; Ruelzheim, Germany) by PCR using primers T2128–T2129. The amplified fragment digested with HindIII was ligated into the HindIII site of each plasmid harboring upstream and downstream regions of a target gene, thus creating plasmids pTA1811 (ndhD1 disruption), pTA1812 (ndhD2 disruption), pTA1641 (ndhD3 disruption), pTA1640 (ndhD4 disruption), and pTA1828 (ndhF1 disruption).
The constructed plasmids (pTA1811, 1812, 1641, 1640, and 1828) were used for gene disruption in TA1297 (wild-type strain), TA2984 (1,3-PDO-producing strain), and TA3800 (glycerol-producing strain). TA2984 and TA3800 were constructed in our previous studies [21, 38]. The strains used in this study are listed in Additional file 1: Table S1.
Culture supernatants were centrifuged (20,000×g, 10 min) and filtered using Minisart RC4 (Sartorius; Goettingen, Germany). Concentrations of glycerol and 1,3-PDO in filtered samples were determined using a high-performance liquid chromatograph (LC-20AD, Shimadzu; Kyoto, Japan) equipped with an autosampler (AOC-20; Shimadzu), a SUGAR SP-G guard column (Shodex; Tokyo, Japan), a SUGAR SP0810 column (Shodex), and a refractive index detector (RID-10A; Shimadzu). Operating conditions were as follows: injection volume, 5 µL; mobile phase, MilliQ water; flow rate, 1.0 mL/min in isocratic mode; column temperature, 80 °C.
In silico simulation of metabolic flux distribution
FBA for TA2984 (1,3-PDO-producing strain)
The stoichiometric metabolic model based on the annotated genome sequence of S. elongatus PCC 7942 published by Knoop  (hereafter called GSM_7942 in this study) served as a basic metabolic model of the wild-type strain in this study. GSM_7942 is composed of 595 genes, 666 reactions, and 533 metabolites. The information on the introduced synthetic metabolic pathway for 1,3-PDO production (8 genes, 6 reactions, and 3 metabolites) was added to GSM_7942, resulting in a metabolic model of TA2984 (hereafter called GSM_7942_PDO in this study). GSM_7942_PDO is composed of 603 genes, 672 reactions, and 536 metabolites. GSM_7942_PDO was used to perform FBA using the COBRA toolbox 2.05  on Matlab (Mathworks Inc., Natick, MA).
Gene disruption simulation in TA2984 with the MOMA algorithm
This simulation was performed on GSM_7942_PDO with scripts provided in COBRA toolbox 2.05 using Gurobi 6.03 (Gurobi optimization Inc., USA) as a solver.
Respiratory and photosynthetic activities
Chlorophyll concentrations were measured by extraction in methanol . Oxygen consumption (respiration) and evolution (photosynthesis) from whole cells were measured with a Clark-type electrode (OXYT-1, Hansatech Instruments, Norfolk, UK) at a chlorophyll concentration of 10 µg/mL. The temperature of the electrode chamber was maintained at 30°C, and the chamber contents were continuously stirred. Saturated light (1000 μE/m2/s) and growth light (100 μE/m2/s) provided by a digital LED source (PD3-5024-4-EI, CCS, Kyoto, Japan) were used for measurement of photosynthetic activity. Photon flux density was measured with IKS-27 (Koito; Tokyo, Japan). Respiratory activity was similarly measured in complete darkness.
Measurement of ATP content
ATP extraction and measurement were basically performed by the method of Wang et al. . The appropriate number of cells at growth phase (corresponding to 1-mL culture at OD730 of 1.0) was collected, and ATP was extracted from pelleted cells by resuspension in a 1% (w/v) trichloroacetic acid solution. After centrifugation at 20,000×g for 10 min at 4 °C, the supernatant was neutralized with 1 M Tris–acetate buffer (pH 7.8). The concentration of ATP in the neutralized sample was determined using the ATP determination kit (A22066, Molecular Probes Inc., Oregon, USA).
Measurement of intracellular metabolites
Intracellular metabolites of TA2984 (1,3-PDO-producing strain) and TA4058 (ndhF1-null 1,3-PDO-producing strain) were measured by LC-QqQ-MS and GC-Q-MS analyses. Sample preparation and analysis were performed by the method described by Soma and Hanai .
Measurement of glycogen content
Glycogen extraction was basically performed by the method of Suzuki et al. . Cells at early stationary phase (days 10) were used for glycogen measurement. The extracted glycogen was mixed with glucoamylase (Oriental Yeast Co., Ltd.; Tokyo, Japan) in 2.5 mM Na-acetate (pH 5.0) and incubated at 40 °C for 1 h. After centrifugation, glucose generated from glycogen was measured using BF-7 (Oji Scientific Instruments, Hyogo, Japan) equipping glucose electrode (ED07-0003, Oji Scientific Instruments).
Results and discussion
FBA of the 1,3-PDO-producing strain
We previously constructed a 1,3-PDO-producing strain, designated as TA2984, by introduction of a synthetic metabolic pathway. A four-step reaction of the introduced pathway converts DHAP, one of the metabolites in the Calvin cycle, into 1,3-PDO via glycerol (Fig. 1). Construction of the stoichiometric metabolic model based on the annotated genome sequence of S. elongatus PCC 7942 was first reported by Triana et al. . In the same year, the comparable scale of a metabolic model of this strain was published by Knoop et al. . We used Knoop’s model composed of 595 genes, 666 reactions, and 533 metabolites as a basic metabolic model of the wild-type strain of S. elongatus PCC 7942 (hereafter called GSM_7942 in this study). The additional information on the introduced synthetic metabolic pathway for 1,3-PDO production (8 genes, 6 reactions, and 3 metabolites) was added to GSM_7942, resulting in a metabolic model of the 1,3-PDO-producing strain composed of 603 genes, 672 reactions, and 536 metabolites (hereafter called GSM_7942_PDO in this study).
Simulation of gene disruption in TA2984 (1,3-PDO-producing strain)
The metabolic flux distribution of TA2984 using FBA was applied to the gene disruption simulation based on the MOMA algorithm . This is one of the algorithms used for simulation of sequential gene deletion assuming that the cell adapts to a mutation to minimize the changes of metabolic flux distribution by the gene deletion. The estimated candidates for gene deletion based on this algorithm successfully increased lycopene production in E. coli . Because the objective of our simulation was to screen promising candidates for gene disruption for improvement of 1,3-PDO productivity, we set the specific production rate of 1,3-PDO obtained experimentally [0.0250 mmol/(g DCW)/h] as the lower limit. The metabolic flux distribution corresponding to each of the 595 gene disruptions in GSM_7942_PDO except for the 8 genes of 1,3-PDO producing synthetic metabolic pathway was simulated and the simulated results on 581 strains with single-gene disruption were plotted according to the estimated specific growth rate and estimated specific production rate of 1,3-PDO (Fig. 3b). The simulated results on the 581 strains with single-gene disruption were roughly subdivided into four groups. These four groups are presented as different symbols in Fig. 3b.
The first group (crosses in the figure) was composed of 74 strains that could not satisfy the constraint conditions and lower limit of 1,3-PDO production. The second group (triangles) included 254 strains showing a higher specific production rate of 1,3-PDO, but a null value of the specific growth rate. This may also mean that each gene disrupted in the 254 strains was essential for growth. It was likely that the disruption of any of the 328 genes classified into these two groups was not applicable to improvement of 1,3-PDO productivity. The metabolic flux distributions of the other 267 strains with single-gene disruption could show a certain specific growth rate and higher specific production rate of 1,3-PDO as compared with the experimental data. We screened the 267 candidates for more desirable ones using the product of the specific growth rate and specific production rate as an indicator. Out of the 267 candidates, 231 showed a lower product of the specific growth rate and specific production rate and were classified into the third group (squares).
Finally, 36 genes were classified into the fourth group (circles) as promising candidates to be disrupted for increasing 1,3-PDO productivity (Additional file 1: Table S3). Only 10 patterns of metabolic flux distribution were obtained from 36 strains in the fourth group. First, in GSM_7942, the annotated-gene information was defined by a gene–protein reaction (GPR) association . Under the definition of the GPR association, each disruption of a subunit of a protein complex indicates the same flux distribution. In this study, 22 genes involved in NADPH dehydrogenase-1 (NDH-1) complexes and four genes involved in cytochrome c oxidase were classified into the fourth group. Second, two cases indicating the same metabolic flux distribution were observed in different two-gene disruptions (Synpcc7942_0191 and 1501, and Synpcc7942_2079 and 2080). These genes were annotated as genes encoding an enzyme catalyzing a serial reaction without any blanching pathway. The other six metabolic flux distributions were examined by single-gene disruption. Based on the value of the product of the specific growth rate and specific production, the disruption of 22 genes involved in the NDH-1 complex was found by screening as the most promising candidate for improvement of 1,3-PDO production among the 10 patterns of the fourth group.
Disruption of genes involved in NDH-1 complexes
NDH-1 complexes in cyanobacteria are involved in several cellular functions (respiration, cyclic electron flow, and carbon dioxide acquisition) and are subdivided into four variants (NDH-1L, L′, MS, and MS′) by function and subunit composition according to reverse-genetic studies . The phenotypic analyses of gene knockout mutants revealed the functional differences among NDH-1 complexes, indicating that NDH-1L and L′ are involved in respiration and cyclic electron flow, whereas NDH-1MS and MS′ participate in carbon dioxide acquisition [48, 49]. Because NdhB is one of the common subunits for all four variants, it was expected that ndhB disruption would be desirable to represent the metabolic flux distribution observed in the gene disruption simulation (Fig. 4b, Additional file 1: Table S3). Nonetheless, all functions of NDH-1 complexes were drastically decreased in the ndhB-null mutant, in particular, this mutant could not grow under ambient air conditions [50–52]. Indeed, ndhB could not be completely disrupted in the wild-type strain and TA2984 (1,3-PDO-producing strain) in our experiment. On the other hand, it was reported that four mutants in which each gene encoding the unique subunit in four NDH-1 complexes is disrupted can grow under photoautotrophic and ambient air conditions . Therefore, we tried to disrupt the genes encoding these unique subunits in TA2984.
Output fluxes observed by experiments (days 2–4) and the simulated output fluxes
Experimental or simulated data
Sum of products
1,3-PDO + glycerol
1,3-PDO/sum of products
TA2984 (parental strain for 1,3-PDO production)
Simulated data of NDH-1 elimination
TA4021 (ndhD1-null strain of TA2984)
TA4022 (ndhD2-null strain of TA2984)
TA3672 (ndhD3-null strain of TA2984)
TA3671 (ndhD4-null strain of TA2984)
TA4058 (ndhF1-null strain of TA2984)
Previously, we have constructed a glycerol-producing strain (TA3800) into which only the pathway converting DHAP to glycerol was introduced . The effect of ndhF1 disruption was also evaluated in TA3800. As in the 1,3-PDO producing strain, the ndhF1 disruption resulted in a growth defect and improvement of the production titer (Additional file 2: Figure S2). The growth defect caused by the ndhF1 disruption in the 1,3-PDO-producing strain was more significant than that in the wild-type strain and glycerol-producing strain (Additional file 2: Figure S2A–C). The highest titer of glycerol produced by the ndhF1-null glycerol-producing strain designated as TA4059 reached 36.7 mM (3.38 g/L; Additional file 2: Figure S2D). This value was almost equal to the sum of the titers of 1,3-PDO and glycerol (34.7 mM) shown by TA4058 (Fig. 5) and about twofold higher than that of TA3800 (Additional file 2: Figure S2E). Although the reason for the significant growth defect caused by ndhF1 disruption in the 1,3-PDO-producing strain was not clear, glycerol production was successfully stimulated by the ndhF1 disruption, as was 1,3-PDO production. The titer of glycerol produced by TA4059 was the highest among other chemicals (glycerol, 1,3-PDO, and 1,2-PDO) produced from DHAP by engineered cyanobacteria [18–20, 38] and this titer corresponded to high titer compared to other chemicals produced not from DHAP by engineered cyanobacteria, such as ethanol , sucrose , 2,3-butanediol , and lactate . Although various chemicals have been produced by engineered cyanobacteria as summarized in some reviews [1, 2], the research resulting from the titer of product over 3 g/L directly from carbon dioxide seems to be still rare.
Metabolism changes caused by ndhF1 disruption
The disruption of ndhF1 successfully increased metabolic flux of DHAP into the 1,3-PDO-producing pathway as expected on the basis of the in silico simulation (Table 1) and improved the titers of 1,3-PDO and glycerol production (Fig. 5 and Additional file 2: Figure S2). This result suggests that in silico simulation is an effective method for metabolic improvement of bioproduction in cyanobacteria. To identify the causes of improved productivity, the metabolic changes caused by the ndhF1 disruption were analyzed. It was expected that the ndhF1 disruption affected to respiration and cyclic electron flow . Therefore, at first, the respiration activities of TA4058 (ndhF1-null 1,3-PDO-producing strain) and TA2984 (parental strain) were measured. As the result, TA4058 showed lower respiratory activity in the dark and almost the same photosynthetic activity under saturated light conditions (1000 μE/m2/s) as compared with TA2984 (Additional file 2: Figure S3). This finding was consistent with the phenotype of a ndhF1-null strain of Synechococcus sp. PCC 7002 observed in another study . Under growth light conditions (100 μE/m2/s), the photosynthetic activity of TA4058 was lower than that of TA2984, being consistent with the growth defect of TA4058 (Fig. 5).
The physiological role of cyclic electron flow is known to increase the ATP/NADPH ratio produced by photosynthetic electron flow . So, next, ATP contents in TA4058 and TA2984 were measured. TA4058 showed lower cellular ATP content as compared with TA2984 (Additional file 2: Figure S4). In the respiratory chain, the reducing power (NADH and FADH2) mainly produced by glycolysis and the TCA cycle is used for generating a proton gradient and H2O from O2. The generated proton gradient functions as a driving force of ATP synthase for ATP production. Whole or partial loss of cyclic electron flow and respiratory activity under the influence of the ndhF1 disruption may decrease cellular ATP content.
The intracellular metabolites involved in a reaction with ATP seemed to reflect the ATP deficiency in TA4058. For example, ATP is consumed in reactions producing RuBP and ADP-glucose, and is produced in the reaction producing pyruvate. In ATP-consuming reactions, substrates (ribulose-5-phosphate: Ru5P and glucose-1-phosphate: G1P) and products (RuBP and ADP-glucose) were detected at high and low concentrations, respectively. Conversely, lower concentration of a substrate (phosphoenolpyruvate: PEP) and higher concentration of a product (pyruvate) were observed in an ATP-producing reaction. In particular, the metabolites around the glycogen synthesis pathway seemed to be affected by the ATP deficiency. Fructose-6-phosphate (F6P) is a metabolite of the branching point for glycogen synthesis from the Calvin cycle and is sequentially converted to glucose-6-phosphate (G6P), glucose-1-phosphate (G1P), and ADP-glucose. In this glycogen synthesis pathway, significant accumulation (more than fivefold) was observed for metabolites except for ADP-glucose produced with ATP consumption. This result suggests that glycogen synthesis is repressed in TA4058 caused by ATP deficiency. Because other studies showed that higher accumulation of glycogen is observed in stationary phase [44, 59], the glycogen concentrations of TA2984 and TA4058 were measured at 10 days (Additional file 2: Figure S6). As expected, the intracellular glycogen concentration of TA4058 was about half of that of TA2984 (Additional file 2: Figure S6A).
It is well known that the respiratory activity is strongly related to catabolic metabolism (glycolysis and TCA cycle) in heterotrophic organisms. We hypothesized that the repressed respiratory activity can decrease the metabolic flux in catabolic metabolism generating reducing power in photosynthetic organisms, too. The decreasing of 3-PG-consuming flux into glycolysis by the ndhF1 disruption in the simulation result (Fig. 4) was consistent with the hypothesis. Quantification of intracellular metabolites also supported this hypothesis (Additional file 2: Figure S5). Most of the metabolites in the Calvin cycle showed higher intracellular concentration in TA4058 than in TA2984. On the other hand, the metabolites in the glycolysis and TCA cycle were mostly lower in TA4058 compared with TA2984.
Not only greater flux of 3-PG into the Calvin cycle but also repression of glycogen synthesis may increase the productivity values of 1,3-PDO and glycerol. These results indicate that the productivity improvement was caused by the metabolism adapted to energy imbalance induced by the ndhF1 disruption. Hasunuma et al.  successfully improved the carbon fixation rate in Synechocystis sp. PCC 6803 by overexpression of the flv3 gene involved in NADPH consumption. The strategy of acceleration of NADPH consumption also showed the energy balance to change cellular metabolism of cyanobacteria. In the present study, the metabolism changes intended to adapt to the energy imbalance (ATP deficiency) were applied to chemical production (1,3-PDO and glycerol). Taking into consideration of the points mentioned above, ndhF1 disruption is a good strategy for enhancement of production of other chemicals (for instance, 1,2-PDO and sucrose) from metabolites in the Calvin cycle in cyanobacteria [18, 38, 55].
The simulation result of the metabolic flux distribution of the ndhF1-null and wild type strain (Fig. 4) showed the 3-PG-consuming flux to 2-PG or TCA cycle was decreased by the ndhF1 disruption. This result suggested that the ndhF1-null strain was not suitable for the chemical productions from pyruvate or acetyl-CoA. Yoshikawa et al. , however, showed the ndhF1-null PCC6803 strain with ethanol producing pathway was produced 145% more ethanol (0.132 g/L) from pyruvate than the wild type strain with the same pathway (0.191 g/L). It is difficult to explain clearly the reason why the suggestion based on our simulation result is different from their result, but one possible reason would be the driving force by ATP generation by the reaction from PEP to pyruvate to compensate the ATP deficiency in the ndhF1-null strain.
Recently, Broddrick et al.  constructed GEM of metabolism in PCC7942 composed of 785 genes. This model was constructed by using the information about essential genes revealed by the random barcode transposon site sequencing and physiological data for photoautotrophic metabolism. Using this model or the larger models than that we used, there are possibilities to screen out the other candidate genes for knockout to improve 1,3-PDO productivity.
In this study, the productivity of 1,3-PDO and glycerol in S. elongatus PCC 7942 were successfully improved by gene manipulation based on in silico simulation of gene disruptions. The best candidate for improvement of 1,3-PDO productivity was found to be the elimination of NDH-1 complexes involved in several cellular functions: respiration, cyclic electron flow, and carbon dioxide uptake (Additional file 1: Table S3). Genes encoding NDH-1 complexes were actually disrupted in TA2984 (parental strain for 1,3-PDO production). The effective improvement of productivity was achieved by disruption of the NDH-1 complexes that are related to respiration and cyclic electron flow (Additional file 2: Figure S1). Finally, TA4058 (ndhF1-null 1,3-PDO-producing strain) was found to produce 1.28-fold more 1,3-PDO (4.44 mM, 0.338 g/L) and 2.13-fold more glycerol (30.3 mM, 2.79 g/L) than TA2984 did (Fig. 5). The metabolism adapted to the energy imbalance caused by the ndhF1 disruption may enhance the 3-PG-consuming flux into the Calvin cycle and can be suitable for 1,3-PDO and glycerol production (Additional file 2: Figure S5). This result opens up exciting opportunities for applying in silico simulation to metabolic improvement of cyanobacteria and is suggestive of the feasibility of metabolic changes focused on the energy balance.
YH, SM, and TH designed the experiments. YH and SM carried out all experiments and SM, HH, and FM carried out theoretical analysis and YH drafted the figures and manuscript. TH edited the manuscript. All authors read and approved the final manuscript.
This study was financially supported by the Core Research of Evolutional Science and Technology program (CREST) from the Japan Science and Technology Agency (JST).
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used in the current study are available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
The Core Research of Evolutional Science and Technology program (CREST) from the Japan Science and Technology Agency (JST).
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.
- Angermayr SA, Rovira AG, Hellingwerf KJ. Metabolic engineering of cyanobacteria for the synthesis of commodity products. Trends Microbiol. 2015;33:352–61.Google Scholar
- Oliver NJ, Rabinovitch-Deere CA, Carroll AL, Nozzi NE, Case AE, Atsumi S. Cyanobacterial metabolic engineering for biofuel and chemical production. Curr Opin Chem Biol. 2016;35:43–50.View ArticleGoogle Scholar
- Atsumi S, Higashide W, Liao JC. Direct photosynthetic recycling of carbon dioxide to isobutyraldehyde. Nat Biotechnol. 2009;27:1177–80.View ArticleGoogle Scholar
- Li X, Shen CR, Liao JC. Isobutanol production as an alternative metabolic sink to rescue the growth deficiency of the glycogen mutant of Synechococcus elongatus PCC 7942. Photosynth Res. 2014;120:301–10.View ArticleGoogle Scholar
- Gao Z, Zhao H, Li Z, Tan X, Lu X. Photosynthetic production of ethanol from carbon dioxide in genetically engineered cyanobacteria. Energy Environ Sci. 2012;5:9857–65.View ArticleGoogle Scholar
- Oliver JWK, Machado IM, Yoneda H, Atsumi S. Cyanobacterial conversion of carbon dioxide to 2,3-butanediol. Proc Natl Acad Sci USA. 2013;110:1249–54.View ArticleGoogle Scholar
- Oliver JWK, Machado IM, Yoneda H, Atsumi S. Combinatorial optimization of cyanobacterial 2,3-butanediol production. Metab Eng. 2014;22:76–82.View ArticleGoogle Scholar
- Shen CR, Liao JC. Photosynthetic production of 2-methyl-1-butanol from CO2 in cyanobacterium Synechococcus elongatus PCC 7942 and characterization of the native acetohydroxyacid synthase. Energy Environ Sci. 2012;5:9574–83.View ArticleGoogle Scholar
- Anfelt J, Kaczmarzyk D, Shabestary K, Renberg B, Rockberg J, Nielsen J, Uhlén M, Hudson EP. Genetic and nutrient modulation of acetyl-CoA levels in Synechocystis for n-butanol production. Microb Cell Fact. 2015;14:167–78.View ArticleGoogle Scholar
- Lan EI, Ro SY, Liao JC. Oxygen-tolerant coenzyme A-acylating aldehyde dehydrogenase facilitates efficient photosynthetic n-butanol biosynthesis in cyanobacteria. Energy Environ Sci. 2013;6:2672–81.View ArticleGoogle Scholar
- Chwa JW, Kim WJ, Sim SJ, Um Y, Woo HM. Engineering of a modular and synthetic phosphoketolase pathway for photosynthetic production of acetone from CO2 in Synechococcus elongatus PCC 7942 under light and aerobic condition. Plant Biotechnol J. 2016;14:1768–76.View ArticleGoogle Scholar
- Zhou J, Zhang H, Zhang Y, Li Y, Ma Y. Designing and creating a modularized synthetic pathway in cyanobacterium Synechocystis enables production of acetone from carbon dioxide. Metab Eng. 2012;14:394–400.View ArticleGoogle Scholar
- Wang B, Pugh S, Nielsen DR, Zhang W, Meldrum DR. Engineering cyanobacteria for photosynthetic production of 3-hydroxybutyrate directly from CO2. Metab Eng. 2013;16:68–77.View ArticleGoogle Scholar
- Lee T-C, Xiong W, Paddock T, Carrieri D, Chang I-F, Chiu H-F, Ungerer J, Juo S-HH, Maness P-C. Engineered xylose utilization enhances bio-products productivity in the cyanobacterium Synechocystis sp. PCC 6803. Metab Eng. 2015;30:179–89.View ArticleGoogle Scholar
- Ungerer J, Tao L, Davis M, Ghirardi M, Maness PC, Yu J. Sustained photosynthetic conversion of CO2 to ethylene in recombinant cyanobacterium Synechocystis 6803. Energy Environ Sci. 2012;5:8998–9006.View ArticleGoogle Scholar
- Bentley FC, Melis Anastasios. Diffusion-based process for carbon dioxide uptake and isoprene emission in gaseous/aqueous two-phase photobioreactors by photosynthetic microorganisms. Biotechnol Bioeng. 2012;109:100–9.View ArticleGoogle Scholar
- Gao X, Gao F, Liu D, Zhang H, Nie X, Yang C. Engineering the methylerythritol phosphate pathway in cyanobacteria for photosynthetic isoprene production from CO2. Energy Environ Sci. 2016;9:1400.View ArticleGoogle Scholar
- Li H, Liao JC. Engineering a cyanobacterium as the catalyst for the photosynthetic conversion of CO2 to 1,2-propanediol. Microb Cell Fact. 2013;12:4.View ArticleGoogle Scholar
- Savakis P, Tan X, Du W, Branco dos Santos F, Lu X, Hellingwerf KJ. Photosynthetic production of glycerol by a recombinant cyanobacterium. J Bacteriol. 2015;195:46–51.Google Scholar
- Wang Y, Tao F, Ni J, Li C, Xu P. Production of C3 platform chemicals from CO2 by genetically engineered cyanobacteria. Green Chem. 2015;17:3100–10.View ArticleGoogle Scholar
- Hirokawa Y, Maki Y, Hanai T. Cyanobacterial production of 1,3-propanediol directly from carbon dioxide using a synthetic metabolic pathway. Metab Eng. 2016;34:97–103.View ArticleGoogle Scholar
- Saxena RK, Anand P, Saran S, Isar J. Microbial production of 1,3-propanediol: recent developments and emerging opportunities. Biotechnol Adv. 2009;27:895–913.View ArticleGoogle Scholar
- Orth JD, Thiele I, Palsson BØ. What is flux balance analysis? Nat Biotechnol. 2010;28:245–8.View ArticleGoogle Scholar
- Hong SH, Lee SY. Metabolic flux analysis for succinic acid production by recombinant Escherichia coli with amplified malic enzyme activity. Biotechnol Bioeng. 2001;74:89–95.View ArticleGoogle Scholar
- Lee KH, Park JH, Kim TY, Kim HU, Lee SY. Systems metabolic engineering of Escherichia coli for l-threonine production. Biol: Mol Syst; 2007. p. 3.Google Scholar
- Park JH, Lee KH, Kim TY, Lee SY. Metabolic engineering of Escherichia coli for the production of l-valine based on transcriptome analysis and in silico gene knockout simulation. Proc Natl Acad Sci USA. 2007;104:7797–802.View ArticleGoogle Scholar
- Choi HS, Lee SY, Kim TY, Woo HM. In silico identification of gene amplification targets for improvement of lycopene production. Appl Environ Microbiol. 2010;76:3097–105.View ArticleGoogle Scholar
- Yen JY, Tanniche I, Fischer AK, Gillaspy GE, Bevan DR, Senger RS. Designing metabolic engineering strategies with genome-scale metabolic flux modeling. Adv Genomics Genet. 2015;5:93–105.Google Scholar
- Shastri AA, Morgan JA. Flux balance analysis of photoautotrophic metabolism. Biotechnol Prog. 2005;21:1617–26.View ArticleGoogle Scholar
- Baroukh C, Muñoz-Tamayo R, Bernard O, Steyer J-P. Mathematical modeling of unicellular microalgae and cyanobacteria metabolism for biofuel production. Curr Opin Biotechnol. 2015;33:198–205.View ArticleGoogle Scholar
- Triana J, Montagud A, Siurana M, Fuente D, Urchueguía A, Gamermann D, Torres J, Tena J, Fernández de Córdoba P, Urchueguía JF. Generation and evaluation of a genome-scale metabolic network model of Synechococcus elongatus PCC 7942. Metabolites. 2014;4:680–98.View ArticleGoogle Scholar
- Yoshikawa K, Aikawa S, Kojima Y, Toya Y, Furusawa C, Kondo A, Shimizu H. Construction of a genome-scale metabolic model of Arthrospira platensis NIES-39 and metabolic design for cyanobacterial bioproduction. PLoS ONE. 2015;7:e0144430.View ArticleGoogle Scholar
- Shabestary K, Hudson EP. Computational metabolic engineering strategies for growth-coupled biofuel production by Synechocystis. Metab Eng Commun. 2016;3:216–26.View ArticleGoogle Scholar
- Shirai T, Osanai T, Kondo A. Designing intracellular metabolism for production of target compounds by introducing a heterologous metabolic reaction based on a Synechocystis sp. PCC 6803 genome-scale model. Microb Cell Fact. 2016;15:13–9.View ArticleGoogle Scholar
- Yoshikawa K, Toya Y, Shimizu H. Metabolic engineering of Synechocystis sp. PCC6803 for enhanced ethanol production based on flux balance analysis. Bioprocess Biosyst Eng. 2017. https://doi.org/10.1007/s00449-017-1744-8.Google Scholar
- Hendry JI, Prasannan CB, Joshi A, Dasgupta S, Wangikar PP. Metabolic model of Synechococcus sp. PCC 7002: prediction of flux distribution and network modification for enhanced biofuel production. Bioresour Technol. 2016;213:190–7.View ArticleGoogle Scholar
- Segrè D, Vitkup D, Church GM. Analysis of optimality in natural and perturbed metabolic networks. Proc Natl Acad Sci USA. 2002;99:15112–7.View ArticleGoogle Scholar
- Hirokawa Y, Maki Y, Hanai T. Improvement of 1,3-propanediol production using an engineered cyanobacterium, Synechococcus elongatus by optimization of the gene expression level of a synthetic metabolic pathway and production conditions. Metab Eng. 2017;39:192–9.View ArticleGoogle Scholar
- Knoop H. Integrative analyse des cyanobakteriellen stoffwechsels. PhD thesis. 2014. https://doi.org/10.18452/17055.
- Schellenberger J, Que R, Fleming RM, Thiele I, Orth JD, Feist AM, Zielinski DC, Bordbar A, Lewis NE, Rahmanian S, Kang J, Hyduke DR, Palsson BØ. Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc. 2011;6:1290–307.View ArticleGoogle Scholar
- Grimme LH, Boardman NK. Photochemical activities of a particle fraction P1 obtained from the green algae Chlorella fucsa. Biochem Biophys Res Commun. 1972;49:1617–23.View ArticleGoogle Scholar
- Wang X, Liu W, Xin C, Zheng Y, Cheng Y, Sun S, Li R, Zhu X-G, Dai SY, Rentzepis PM, Yuan JS. Enhanced limonene production in cyanobacteria reveals photosynthesis limitations. Proc Natl Acad Sci USA. 2016;113:14225–30.View ArticleGoogle Scholar
- Soma Y, Hanai T. Self-induced metabolic state switching by a tunable cell density sensor for microbial isopropanol production. Metab Eng. 2015;30:7–15.View ArticleGoogle Scholar
- Suzuki E, Ohkawa H, Moriya K, Matsubara T, Nagaike Y, Iwasaki I, Fujiwara S, Tsuzuki M, Nakamura Y. Carbohydrate metabolism in mutants of the cyanobacterium Synechococcus elongatus PCC 7942 defective in glycogen synthesis. Appl Environ Microbiol. 2010;76:3153–9.View ArticleGoogle Scholar
- Alper H, Jin Y-S, Moxley JF, Stephanopoulos G. Identifying gene tragets for the metabolic engineering of lycopene biosynthesis in Escherichia coli. Metab Eng. 2005;7:155–64.View ArticleGoogle Scholar
- Thiele I, Palsson BØ. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc. 2010;5:93–121.View ArticleGoogle Scholar
- Battchikova N, Eisenhut M, Aro E-M. Cyanobacterial NDH-1 complexes: novel insights and remaining puzzles. Biochim Biophys Acta. 2011;1807:935–44.View ArticleGoogle Scholar
- Ohkawa H, Pakrasi HB, Ogawa T. Two types of functionally distinct NAD(P)H dehydrogenases in Synechocystis sp. strain PCC 6803. J Biol Chem. 2000;275:31630–4.View ArticleGoogle Scholar
- Shibata M, Ohkawa H, Kaneko T, Fukuzawa H, Tabata S, Kaplan A, Ogawa T. Distinct constitutive and low-CO2-induced CO2 uptake systems in cyanobacteria: genes involved and their phylogenetic relationship with homologous genes in other organisms. Proc Natl Acad Sci USA. 2001;98:11789–94.View ArticleGoogle Scholar
- Marco E, Ohad N, Schwarz R, Lieman-Hurwitz J, Gabay C, Kaplan A. High CO2 concentration alleviates the block in photosynthetic electron transport in an ndhB-inactivated mutant of Synechococcus sp. PCC 7942. Plant Physiol. 1993;101:1047–53.View ArticleGoogle Scholar
- Mi H, Endo T, Schreiber U, Ogawa T, Asada K. Electron donation from cyclic and respiratory flows to the photosynthetic intersystem chain is mediated by pyridine nucleotide dehydrogenase in the cyanobacterium Synechocystis PCC 6803. Plant Cell Physiol. 1992;33:1233–7.Google Scholar
- Ogawa T. A gene homologous to the subunit-2 gene of NADH dehydrogenase is essential to inorganic carbon transport of Synechocystis PCC 6803. Proc Natl Acad Sci USA. 1991;88:4275–9.View ArticleGoogle Scholar
- Zhang P, Battchikova N, Jansen T, Appel J, Ogawa T, Aro E-M. Expression and functional roles of the two distinct NDH-1 complexes and the carbon acquisition complex NdhD3/F3/CupA/Sll1735 in Synechocystis sp. PCC 6803. Plant Cell. 2004;16:3326–40.View ArticleGoogle Scholar
- Wulfhorst H, Franken LE, Wessinghage T, Boekema EJ, Nowaczyk MM. The 5 kDa protein NdhP is essential for stable NDH-1L assembly in Thermosynechococcus elongatus. PLoS ONE. 2014;9:e103584.View ArticleGoogle Scholar
- Ducat DC, Avelar-Rivas JA, Way JC, Silver PA. Rerouting carbon flux to enhance photosynthetic productivity. Appl Environ Microbiol. 2012;78:2660–8.View ArticleGoogle Scholar
- Angermayr SA, Hellingwerf KJ. On the use of metabolic control analysis in the optimization of cyanobacterial biosolar cell factories. J Phys Chem B. 2013;117:11169–75.View ArticleGoogle Scholar
- Schluchter WM, Zhao J, Bryant DA. Isolation and characterization of the ndhF gene of Synechococcus sp. strain PCC 7002 and initial characterization of an interposon mutant. J Bacteriol. 1993;175:3343–52.View ArticleGoogle Scholar
- Rochaix J-D. Reprint of: regulation of photosynthetic electron transport. Biochem Biophys Acta. 2011;1807:876–86.Google Scholar
- Song K, Tan X, Liang Y, Lu X. The potential of Synechococcus elongatus UTEX 2973 for sugar feedstock production. Appl Microbiol Biotechnol. 2016;100:7865–75.View ArticleGoogle Scholar
- Hasunuma T, Matsuda M, Senga Y, Aikawa S, Toyoshima M, Shimakawa G, Miyake C, Kondo A. Overexpression of flv3 improves photosynthesis in the cyanobacterium Synechocystis sp. PCC 6803 by enhancement of alternative electron flow. Biotechnol Biofuels. 2014;7:493–503.View ArticleGoogle Scholar
- Broddrick JT, Rubin RE, Welkie DG, Du N, Mih N, Diamond S, Lee JJ, Golden SS, Palsson BO. Unique attributes of cyanobacterial metabolism revealed by improved genome-scale metabolic modeling and essential gene analysis. Proc Natl Acad Sci USA. 2016;113:E8344–53.View ArticleGoogle Scholar