- Open Access
Deciphering cyanobacterial phenotypes for fast photoautotrophic growth via isotopically nonstationary metabolic flux analysis
Biotechnology for Biofuelsvolume 10, Article number: 273 (2017)
Synechococcus elongatus UTEX 2973 is the fastest growing cyanobacterium characterized to date. Its genome was found to be 99.8% identical to S. elongatus 7942 yet it grows twice as fast. Current genome-to-phenome mapping is still poorly performed for non-model organisms. Even for species with identical genomes, cell phenotypes can be strikingly different. To understand Synechococcus 2973’s fast-growth phenotype and its metabolic features advantageous to photo-biorefineries, 13C isotopically nonstationary metabolic flux analysis (INST-MFA), biomass compositional analysis, gene knockouts, and metabolite profiling were performed on both strains under various growth conditions.
The Synechococcus 2973 flux maps show substantial carbon flow through the Calvin cycle, glycolysis, photorespiration and pyruvate kinase, but minimal flux through the malic enzyme and oxidative pentose phosphate pathways under high light/CO2 conditions. During fast growth, its pool sizes of key metabolites in central pathways were lower than suboptimal growth. Synechococcus 2973 demonstrated similar flux ratios to Synechococcus 7942 (under fast growth conditions), but exhibited greater carbon assimilation, higher NADPH concentrations, higher energy charge (relative ATP ratio over ADP and AMP), less accumulation of glycogen, and potentially metabolite channeling. Furthermore, Synechococcus 2973 has very limited flux through the TCA pathway with small pool sizes of acetyl-CoA/TCA intermediates under all growth conditions.
This study employed flux analysis to investigate phenotypic heterogeneity among two cyanobacterial strains with near-identical genome background. The flux/metabolite profiling, biomass composition analysis, and genetic modification results elucidate a highly effective metabolic topology for CO2 assimilatory and biosynthesis in Synechococcus 2973. Comparisons across multiple Synechococcus strains indicate faster metabolism is also driven by proportional increases in both photosynthesis and key central pathway fluxes. Moreover, the flux distribution in Synechococcus 2973 supports the use of its strong sugar phosphate pathways for optimal bio-productions. The integrated methodologies in this study can be applied for characterizing non-model microbial metabolism.
Efforts toward a sustainable bio-based economy have focused on phototrophic hosts for the production of chemicals from CO2 and light. Cyanobacteria are of special biotechnological interest due to their metabolic flexibility, unicellular anatomy, and high photosynthetic efficiency. Additionally, cyanobacteria can use flue gas from power plants to mitigate CO2 emissions that contribute to climate change . To achieve metrics required for commercialization, cyanobacterial photo-biorefineries must have comparable biosynthesis capability to commonly used heterotrophic organisms. Recently, Synechococcus elongatus UTEX 2973 was isolated, whose growth rate reaches a doubling time of 2 h under high light and high CO2 conditions . In comparison, under optimal growth conditions, S. elongatus PCC 7942 exhibits a doubling time of ~ 5 h although its genome sequence is 99.8% identical to Synechococcus 2973 (55 single nucleotide polymorphisms and a large 188 kb inversion between Synechococcus 2973 and 7942) . To understand how cyanobacteria achieve maximal growth rates, this study describes Synechococcus 2973 flux topology under both optimal and suboptimal growth conditions.
Metabolic flux analysis (MFA) can provide a quantitative description of the metabolic network, link genome profiling to phenome analysis, and reveal pathway regulations through comparative studies. Currently, the cyanobacterial strain, Synechocystis sp. PCC 6803 (doubling time ~ 8 h), is considered to be the model cyanobacterium whose metabolism has been extensively profiled by flux analysis tools [3,4,5,6]. Synechocystis 6803 has significant flux through malic enzyme and oxidative pentose phosphate pathways (OPPP) under the photoautotrophic and photomixtrophic conditions. It also operates a cyclic TCA cycle via the γ-aminobutyric acid (GABA) shunt, which forms succinate through the intermediates glutamate and succinate semialdehyde, despite the missing enzyme from 2-oxoglutarate to succinyl-CoA . To profile Synechocystis 6803 photoautotrophic metabolism, 13C-bicarbonate pulse experiments and isotopically nonstationary metabolic flux analysis (INST-MFA) were developed . Using the software package, INCA, mass isotopomer data from dynamic labeling experiments can be used to quantify fluxes without the need to precisely determine metabolite pool sizes (which are fitted as parameters to account for transient labeling data) , and is therefore more convenient than other flux profiling methods . In the current study, INST-MFA, gene knockouts, and metabolite analysis were performed to obtain insights into the physiology and metabolic regulations of Synechococcus 2973 under different bioreactor conditions. Meanwhile, aspects of biomass composition were measured to reveal changes in macromolecule tradeoffs that correlate to cell growth and bioreactor conditions . The outcome highlights the advantages and hurdles of establishing Synechococcus 2973 as a new platform organism for bioproduction. The comparative studies among metabolisms of different cyanobacterial strains may also offer new insights into flux dependency on adaptive evolution.
Synechococcus 2973 growth and biomass compositions
In optimal photobioreactor (PBR, 500 μmol photons/m2s continuous light) conditions , Synechococcus 2973 exhibits very rapid growth (0.33 ± 0.05 h−1; Fig. 1a). For comparison, maximal growth rate of Synechococcus 7942 was 0.14 ± 0.02 h−1 at 300 μmol photons/m2s continuous light due to photo-inhibition at greater light intensities. Synechococcus strains were also grown under suboptimal light conditions in shake flasks (100 μmol photons/m2s continuous light or diurnal light irradiation; Fig. 1b). Suboptimal cultivations greatly impaired cyanobacterial photosynthesis, while Synechococcus 2973 still demonstrated moderately faster growth rate than Synechococcus 7942. Based on growth and biomass composition analyses, Synechococcus 2973 assimilated ~ 12.2 mmol-C/gDCW/h in the PBR but only ~ 6.7 mmol-C/gDCW/h under the shaking flask (SF) conditions. The estimated values of carbon fixation were further confirmed by the measured net CO2 uptake rate using gas chromatography methods (Additional file 1: Figure S1) . This study additionally considers two Synechococcus 2973 mutants (∆zwf and ∆pgl) whose genes from the oxidative pentose phosphate pathway (OPPP) were removed. The two mutants showed unimpaired growth rates compared to the wild-type strain in our tested conditions (Fig. 1).
The composition of Synechococcus 2973 biomass was measured in both PBR and SF conditions under continuous light (Additional file 1: Part 1). The total protein was ~ 50% of cyanobacterial biomass, while protein from the PBR had more glutamate/glutamine than that from the SF (p value 0.054 using two-tailed equal variance Student’s t Test). Synechococcus 2973 contained 9 ± 1% lipids and 1.5 ± 0.5% glycogen in the PBR and 13 ± 2% lipid and 6.0 ± 1.0% glycogen in SF conditions. For comparison, Synechococcus 7942 biomass composition was quantified under optimal PBR conditions and comprised 11 ± 1% lipids, 13 ± 4% glycogen, and 41 ± 0.4% protein. Proteinogenic glutamate/glutamine content in Synechococcus 7942 was significantly lower (p < 0.05) than that in Synechococcus 2973 biomass regardless of conditions. Fatty acid compositions of lipids in both strains were predominantly C16:1 and C16:0 fatty acids (85%), followed by C18:1 and C18:0. Relevant to photosynthetic light-harvesting, less chlorophyll a was produced in PBR conditions (4.4 ± 0.5 µg/mL/OD730 for 2973; 5.2 ± 0.6 µg/mL/OD730 for 7942) compared to SF (7.4 ± 0.3 µg/mL/OD730 for 2973, p < 0.01). Increased chlorophyll levels in SF cells may be a cellular adjustment to compensate for insufficient light conditions.
INST-MFA of photoautotrophic metabolism
INST-MFA relies on measurement of transient labeling data of central metabolites after a 13C-pulse. To quench metabolism, we used ice-cold media and liquid N2 bath to quickly freeze time-course samples. The approach avoided a traditional cold methanol quenching strategy that causes significant metabolite leakage for gram-negative bacteria (i.e., cyanobacteria)  and resulted in improved LC–MS peak quality (Additional file 1: Figure S2a, b) . The INCA software package (isotopomer network compartmental analysis)  was used to calculate flux values (Fig. 2) based on labeling data and the metabolic network (Additional file 1: Tables S1, S2). The fitted model for the Synechococcus 2973 in the PBR was statistically acceptable [Sum of square residuals SSR = 713, with an accepted range of (597 740)], whereas the best fit of the SF model could not be completely explained by measurement errors (SSR = 981, with an accepted range of [382 471]). The fitted model for Synechococcus 7942 in the PBR was statistically acceptable (SSR = 616, range of [543 679]). Biomass composition information and CO2 uptake rate were used to improve the accuracy of the flux maps. Flux-partitioning across conditions and strains was compared through relative flux values that were normalized to CO2 uptake rate and are plotted in Fig. 2 (Confidence intervals and dilution factors are described in Additional file 1: Tables S3, S4, S5; Additional file 1: Figures S3, S4 show isotopomer fitting quality for individual metabolites). Additionally, labeling dynamics for citrate and malate are compared between Synechococcus 2973 and 7942 in Fig. 5. 13C-enrichment for key metabolites is presented in Fig. 5; Additional file 1: Figures S5, S6. The slower labeling pattern of certain metabolites may be indicative of the presence of metabolically inactive pools in bulk cytoplasm. To account for this problem, INST-MFA has employed dilution factors . The isotopic dilutions in certain metabolites during non-stationary labeling suggest a non-homogeneous intracellular environment.
Synechococcus flux maps in PBR and SF cultures
Figure 2 provides a comparison of the flux distributions between Synechococcus 2973 PBR and SF cultures and Synechococcus 7942 PBR cultures. Although the cells have significantly different CO2 fixation rates and biosynthesis fluxes (i.e., Synechococcus 2973 PBR cultures contained more protein synthesis fluxes, while Synechococcus 2973 SF and Synechococcus 7942 PBR cultures had higher carbon allocation to glycogen), the flux ratios in central metabolism (the Calvin cycle and the TCA pathways) were similar after normalization to CO2 uptake rate. Regardless of the strain and growth rate, Synechococcus did not exhibit flux through OPPP, and approximately 12–14% of the fixed carbon was directed to the TCA cycle similar to another recent study of Synechococcus 7942 . Here, Synechococcus had no flux from aKG to FUM, resulting in a linear TCA pathway that is consistent with those in leaves or other photosynthetic tissues [3, 4, 15]. The flux distribution explained the absence of a number of key genes in Synechococcus 2973 genome, including succinyl-CoA ligase, malate dehydrogenase, α-ketoglutarate dehydrogenase, and GABA shunt genes. Although cyanobacteria are thought to employ a cyclic TCA cycle via GABA shunt for respiratory energy production , a complete TCA cycle seems non-essential for the Synechococcus strains possibly due to evolutionary non-necessity or evolved energy balancing . Based on genome annotation, Synechococcus synthesizes fumarate from aspartate via purine metabolism, which can then be converted into malate and succinate (dynamic labeling data shown in Additional file 1: Figure S5) providing a possible source of several TCA organic acids. Synechococcus 2973 had significant anaplerotic fluxes through PEP carboxylase for the synthesis of TCA cycle intermediates, but reduced flux through the malic enzyme (less than ~ 1% of total fixed CO2 in both the PBR and SF conditions). Among tested culture conditions, there were noticeable differences in respiratory fluxes. Specifically, photorespiration in Synechococcus 2973 was elevated in the PBR condition possibly to rebalance carbon and dissipate excessive reducing equivalents generated by high light. In contrast, Synechococcus 7942 demonstrated lower photorespiration activity under PBR conditions.
Isotopic metabolite profiling for relative pool size estimation
We analyzed the change of metabolite concentrations in different cyanobacterial cultures under continuous light PBR conditions . Metabolite leakage for gram negative cyanobacterial cells makes precise determination of cyanobacterial pool sizes difficult; therefore, we used isotopic ratio method to benchmark the change of cyanobacterial metabolite pools against intracellular metabolites from a standard generated from E. coli. Figure 3 shows the relative ratio of metabolite pool sizes normalized by the amount of biomass. Compared to E. coli, cyanobacterial R5P/Ru5P and other sugar phosphates were abundant, whereas the pool sizes of TCA metabolites were significantly low (Fig. 3c). Synechococcus 2973 had generally smaller metabolite pools than Synechocystis 6803 (with the exception of R5P/Ru5P) and similar metabolite pools compared to Synechococcus 7942 (Fig. 3c) with the exception of UDP glucose (i.e., Synechococcus 7942 contained more UDP glucose for carbohydrate synthesis) (Fig. 3b). Notably, Synechococcus 2973 had the highest NADPH concentration among the three cyanobacterial species (Fig. 3a). In addition, Synechococcus 2973 PBR cultivations had decreased central metabolite concentrations (mainly sugar phosphates) relative to SF cultures (Fig. 3d). Comparing two growth conditions of Synechococcus 2973 (Fig. 3d), the SF cells contained more ADP and AMP than PBR cells, which may suggest a change in energy charges (relative ATP concentration to overall adenosine phosphate concentrations) under suboptimal light conditions.
Effect of exogenous organic acids on Synechococcus 2973’s growth
Synechococcus 2973 shows significant sugar phosphate interconversion, but the TCA cycle is incomplete with limited flux and reduced metabolite pools. We hypothesized that the addition of exogenous organic substrates may alleviate such biosynthesis bottlenecks; however, cell growth was not enhanced after supplying OAA, malate, αKG or organic acids (pyruvate and acetate) in either PBR or SF conditions (Fig. 4a), even though the cells demonstrated the ability to incorporate exogenous nutrients into proteinogenic amino acids (especially for OAA; Fig. 4b). Unlabeled malate from the culture medium could be significantly incorporated into aspartate, suggesting an uncharacterized malate dehydrogenase to synthesize oxaloacetate (the precursor of aspartate). In summary, Synechococcus 2973 shows limited capability to efficiently utilize organic acids and convert these nutrients into cyanobacteria biomass .
Flux comparisons among cyanobacterial species
Synechococcus 2973 and 7942 have a number of phenotypic differences. The CO2 uptake rate in Synechococcus 2973 was twofold greater than Synechococcus 7942. Synechococcus 2973 showed up to a 3.5-fold greater 13C-enrichment over Synechococcus 7942 for metabolites such as FBP and S7P (Additional file 1: Figure S6) at 40 s, indicating greater Calvin cycle fluxes. From the normalized flux map, there are differences in the carbon partitioning as presented in Fig. 2 and Additional file 1: Figure S7. Synechococcus 2973 fluxes through fructose-biphosphate aldolase and phosphatase were increased while the fluxes towards glycogen and the malic enzyme were reduced. Flux from F6P to G6P, which is the key step towards OPPP and sugar storage metabolism, was significantly decreased in Synechococcus 2973 under both conditions. Malate and citrate were labeled rapidly in Synechococcus 2973 due to the enhanced PEPC (phosphoenolpyruvate carboxylase) reaction. Unlike Synechocystis 6803 that maintains a functional GABA shunt , Synechococcus elongatus apparently does not contain the genes to operate a complete TCA cycle (as annotated by the KEGG Genome database), but based on labeling data presented here may utilize a linear pathway from OAA to make citrate and aKG. The GABA shunt contains one decarboxylating reaction and the absence of this shunt reduces CO2 loss. A succinate dehydrogenase deletion mutant was created, and the mutant showed no growth defects under continuous light conditions, which supports the absence of fluxes from aKG to fumarate. Unlike Synechocystis 6803, Synechococcus 2973 pyruvate is predominately made from flux through pyruvate kinase with little flux through the malic enzyme (MAL → PYR + CO2). Pyruvate kinase is a key step for producing pyruvate, but its activity in some cyanobacterial species can be inhibited by high ATP/ADP ratios during photosynthesis . For example, the majority of pyruvate synthesis in Synechocystis 6803  and Synechococcus 7942  during slow growth with suboptimal conditions was a result of malate degradation. Overexpression of pyruvate kinase for pyruvate synthesis improved cyanobacterial lactate and alcohol productions [20, 21]. Synechococcus 2973 has naturally evolved to overcome this pyruvate synthesis bottleneck, de-coupling pyruvate kinase inhibition from high photosynthesis rates. Furthermore, as shown in Fig. 2 and genome scale modeling , strong photorespiration is essential for high light photosynthesis and maximal Synechococcus 2973 growth despite its net carbon loss. Synechococcus 2973 exhibits twice the O2 evolution rate of Synechococcus 7942 . Water splitting results in heightened local oxygen concentration and increased linear electron flow, both of which may increase the need for photorespiration or cyclic electron flow .
The OPPP can generate NADPH and oxidize sugar phosphates as a means to derive more reduced compounds but its activation at the same time as the Calvin Cycle would generate a futile cycle that reduces the efficient use of carbon and energy, which is why in higher plants the two are coordinately regulated by multiple mechanisms including redox state and pH [24, 25]. Due to lack of sufficient regulation or Calvin cycle bottlenecks, slow growing Synechocystis 6803 demonstrated ~ 12% CO2 loss through OPPP under continuous light , while Synechococcus does not exhibit significant flux through this pathway. To further test this observation, we constructed Synechococcus 2973 zwf and pgl deletion mutants by disrupting the gene encoding G6P dehydrogenase and 6-phosphogluconolactonase, respectively, so that G6P could not be converted to 6PG (Additional file 1: Figure S2e, f). The mutants showed similar growth rate to the wild-type strain under continuous light and diurnal bioreactor conditions (Fig. 1), indicating the OPPP did not provide a benefit to the organism and was futile for photoautotrophic metabolism in Synechococcus 2973 . These mutants and metabolic insights from the quantification of central carbon fluxes suggest new strategies for the redirection of carbon flux for the bioengineering of photosynthetic organisms. The observation that Synechococcus 2973 exhibits a similar growth phenotype to Synechococcus 7942 under diurnal conditions emphasizes the differences that can occur as a result of changes in light provided by optimal bioreactor conditions.
The optimal photoautotrophic metabolism (reduced pool sizes, enhanced energy levels, and repartitioning of biomass composition)
Metabolite concentrations are tied to input/output fluxes and can provide insight into biochemical level regulation [27, 28]. The lack of correlation between pool sizes and absolute fluxes has been shown in literature [15, 29], as well as here between Synechococcus 2973, Synechococcus 7942, and Synechocystis 6803. Compared to SF cultures and Synechocystis 6803, the concentrations of a number of central metabolites from glycolysis and TCA in Synechococcus 2973 were reduced during the PBR conditions. Less metabolite accumulation in central metabolism under PBR conditions could indicate reduced feedback or higher rates of metabolite turnover. From literature, it has been suggested that as cell growth rates and carbon fluxes increase, the fraction of anabolic enzymes is increased to pull central metabolites towards biomass synthesis and thus decreases their pool sizes . Conversely, the heterotroph E. coli had intracellular acetyl-CoA/organic acid concentrations much higher than cyanobacteria (Fig. 3c). Such observations, together with cultivation experiments using organic acids, indicated that the TCA cycle and its associated anabolic pathways in Synechococcus 2973 operate at lower rates than in E. coli, which would limit production of chemicals from the cyanobacterial TCA cycle . In Synechococcus 2973, the photosynthetic capacity generates sugar phosphates readily (FBP phosphatase/aldolase, fructose-bisphosphatase, transketolase, and RuBisCO reactions) that can enhance its growth [31,32,33] and could potentially be applied for biotechnological productions .
Additionally, the higher energy charge in the PBR than SF for Synechococcus 2973 may benefit anabolic metabolism and cell tolerance to stress conditions . Cyanobacterial biosynthesis has been closely tied to photosynthesis and changing bioreactor conditions . In SF cultures, there was an increase in lower energy molecules (NADP+, ADP, and AMP) inhibiting anabolism, possibly resulting in higher levels of some central metabolites. The PBR grown Synechococcus 2973 showed the highest NADPH level that may reflect enhanced light-harvesting that can facilitate organic carbon assimilation and biomass growth [20, 37]. Synechococcus 2973 adjusts its PSII/PSI ratio throughout growth although it was found to have a slightly lower PSII/PSI ratio than Synechococcus 7942 while exhibiting greater O2 evolution . Additionally, it may have a more optimal ATP/NADPH ratio to sustain metabolism while not exceeding ATP demand , which may be important to further unlocking cyanobacteria growth constraints. While CO2 fixation rate is a key factor for differences in phenotypes, CO2 fixation rate is influenced by optimal energy balance, biomass distribution/central carbon fluxes, and compartmentalization of photosynthetic reactions.
Synechococcus 2973 PBR cells produce little glycogen (8 ~ 10 times less than PBR cells of Synechococcus 7942 and Synechocystis 6803 ). When less carbon is allocated to make glycogen, more can be partitioned toward synthesis of biomass, resulting in greater protein levels in Synechococcus 2973 than 7942. Consequently, Synechococcus 2973 may have more photosystems/RuBisCO proteins for effective photosynthesis. In contrast, Synechococcus 2973 produced significantly more glycogen under suboptimal SF conditions, leading to longer cell doubling times.
Metabolic features of Synechococcus 2973 for bio-production applications
Understanding and exploring native pathways with high metabolic strengths is a promising direction for future microbial cell factories . Our flux results indicate that Synechococcus 2973 is advantageous for the synthesis of targeted products from sugar phosphate pathways under optimal bioreactor conditions. Recently, researchers have engineered Synechococcus 2973 to produce ~ 9 g L−1 of sucrose under potassium chloride stress . This engineered strain also demonstrated potential to produce valuable polysaccharide products. It has been observed that while cyanobacteria have a flexible photosynthetic metabolism, there are a high number of essential metabolic genes involved in photoautotroph growth. This presents additional challenges for mutant strain generation or extensive genetic modification, and the knowledge gained from this study will help surmount the challenges facing cyanobacterial photo-biorefineries .
Subpopulations, metabolic inactive pools, and substrate channeling
Modeling results indicated the presence of unlabeled pools that required additional dilution factors for isotopic data fitting. For SF cultures, an increase in dilution factors was necessary to improve fitting, still resulting in an unaccepted SF fit. The dilution of the metabolite active pool by unlabeled metabolites infers phenotypic heterogeneity, i.e., the presence of inactive/non-growth cells or the presence of pre-existing sources of carbon that are more slowly turned over after 13C-bicarbonate pulses . Inactive cells would not be involved in significant photosynthesis within the duration of the labeling period and may induce reflux of unlabeled carbon into the metabolic network . Inactive subpopulations can be increased by poor light transmittance due to suboptimal culture conditions, cell self-shading and/or insufficient mixing conditions . This may explain why Synechococcus 2973 does not exhibit significant growth advantages in SF conditions.
Metabolite channeling could also contribute to the observed phenotype. Channels pass metabolite intermediates between enzymes without intracellular diffusion. Channeling increases pathway efficiency  and is often associated with microcompartments (i.e., carboxysomes) or enzyme proximity [45,46,47]. For fast growing E. coli species, glycolysis channeling has been evident . Substrate channeling can be inferred from transient labeling experiments . When tracing 13C-enrichment from 3PGA to downstream metabolites in Synechocystis 6803, it was previously found that certain downstream metabolites could be labeled faster than precursors (e.g., PEP was labeled faster than 3PGA) [3, 29]. In this study, we observed 13C-enrichment might not completely follow expected precursor-product relationships (Additional file 1: Figure S6). Comparing to Synechococcus 7942, Synechococcus 2973 showed rapid 13C-enrichment in S7P over precursors. This can be explained by possible channels in the Calvin cycle that may influence the labeling of certain metabolites and provide advantages for CO2 fixation. Synechococcus 2973 also shows evidence for channeling of glycolysis intermediates towards the TCA pathway with faster citrate M+4 distribution over time that exceeds the 13C-enrichment levels of 3PGA (Fig. 5). However, whether differences are a result of traditional channeling definitions or represent the presence of multiple pools separated through additional spatial compartments cannot be conclusively determined from our studies. These observations suggest a heterogenetic distribution of intracellular metabolites as well as possible spatial organization of enzymes to confer growth benefits.
Synechococcus 2973 demonstrates high photosynthesis, efficient carbon fixation and small fluxes towards carbon loss and transient storage pathways under optimal growth conditions. These metabolic traits suggest this strain may be a promising platform to produce high-value compounds in well-controlled PBR conditions. Under suboptimal conditions, Synechococcus 2973 has inferior photosynthesis and decreased biomass synthesis, causing an accumulation of central carbon metabolites. Moreover, Synechococcus strains with different growth rates maintain similar flux distributions in central pathways. These observations offer insight into flux response to adaptive evolution. From a technological perspective, INST-MFA of cyanobacteria is not only important for genome-to-phenome mapping but also crucial for the rational application of platform photo-biorefineries. Currently, INST-MFA is still challenging, and flux results may be influenced by substrate channeling or subpopulations. Sorting out these factors will require advanced labeling experiments [49,50,51,52] to complement computer modeling and labeling experiments. Therefore, mutant creations and biomass composition analysis were performed to offer information resources to aid this study.
Cultivation conditions, and transient labeling experiments
Synechococcus 2973 was grown in BG-11 medium (pH 8.0–8.5) at 38 °C with 3% CO2 aeration (2000 mL min−1) and continuous 500 μmol photons/m2s under photobioreactors (Additional file 1: photo 1). Synechococcus 2973 was also grown in shaking flasks (atmospheric CO2, 38 °C and continuous and diurnal (12 h intervals) 100 μmol photons/m2s light, 250 rpm). Synechococcus 7942 was grown in identical conditions, except in its optimal PBR conditions, where the light condition was 300 μmol photons/m2s. Higher light irradiance inhibits Synechococcus 7942 growth . For the isotopic pool-size experiment, Synechocystis 6803 was grown in identical photobioreactors at 30 °C, 3% CO2 aeration and continuous 300 μmol photons/m2s light. For each labeling experiment, cultures were diluted to an OD730 0.05 in BG-11 medium without citrate (so citrate could not be used a carbon source or affect bulk citrate labeling dynamics) and grown to OD730 0.6 (exponential phase) before the experiment was initiated. Prior to the experiment, the 3% CO2 aeration was replaced with air (0.04% CO2) and a 2-mL aliquot of saturated NaH13CO3 (> 98% purity, Sigma Aldrich, St Louis) was injected into each PBR or shaking flask (SF) for a final concentration of 4 g L−1 of NaH13CO3 to saturate the cell with CO2. The use of 13C-bicarbonate rather than 13CO2 is sufficient because 13C-bicarbonate quickly equilibrates within the media without gas–liquid mass transfer limitations.
After the 13C-pulse, each identical culture in the PBR and SF was quenched at different time intervals (20, 30, 40, 60, 90, 150, 300, 600, 1200, 3600, 7200 s). These time-courses samples were analyzed to capture labeling dynamics of metabolites. Timed biological duplicates were used to generated standard deviations for experimentally measured MID values; however, a minimum standard deviation of 3% was used in INCA model.
Metabolite quenching, extraction, and analysis
Fast-cooling via cold solvents followed by centrifugation is reliable for recovering intracellular labeled metabolites . After a 13C-pulse, samples were quenched by mixing with ice-cold minimal medium (< −5 °C, without a carbon source) in a falcon tube, which was immediately bathed in liquid N2 to sustain cold temperature (To avoid ice formation, the culture sample was rapidly swirled for several seconds). Then the quenched biomass samples (< 0 °C) were harvested via refrigerator centrifuge (3 min, 6000 rpm). For LC–MS analysis, the pelleted samples were extracted with a chloroform–methanol method . Ion-pairing LC–MS/MS was performed at the Proteomics and Mass Spectrometry Facility, Donald Danforth Plant Science Center, St. Louis (details in Additional file 1: Part 1). The labeling of energy molecules (e.g., NADH and ATP), organic acids and acetyl-CoA were quantified using hydrophilic interaction liquid chromatography (HILIC) coupled to electrospray time-of-flight MS  at Joint Bio-Energy Institute, CA.
Estimation of relative pool sizes of metabolites
Changes of metabolite pool sizes in cyanobacteria were analyzed by an MS isotopomer ratio approach  (using fully labeled cell extracts as internal standards for semi-quantitative metabolomics ). Specifically, E. coli K-12 was cultured with uniformly labeled 13C-glucose and 13C-sodium bicarbonate in an M9 minimal media. The estimation of metabolite pool sizes was performed by mixing a known amount of labeled E. coli biomass with unlabeled cyanobacteria cultures (three biological replicates). Then, the mixtures were quenched by a liquid N2 bath and extracted for HILIC-MS analysis (Additional file 1: Figure S2c). The isotopic ratio of each metabolite (labeled vs. unlabeled) was normalized by the amount of E. coli and cyanobacterial biomass, respectively.
Biomass composition analysis
The cells were harvested during exponential growth by centrifugation, and cell pellets were washed with 0.9% NaCl and ddH2O and then freeze-dried. Protein and amino acid compositional analysis was performed by the Molecular Structure Facility, University of California (Davis, CA). Carbohydrates, lipids, chlorophyll a, and ash were measured by previously reported methods (details in Additional file 1: Part 1).
GC–MS analysis of proteinogenic amino acids
To confirm the capability of Synechococcus 2973 to uptake organic acids, unlabeled seed cultures were inoculated (4% inoculation ratio) into the medium with fully labeled NaH13CO3 and unlabeled substrates (6 mM) including acetate, TCA cycle intermediates (e.g., malate and citrate), and pyruvate in shaking flasks or PBRs for cultivation of 48 h. The labeling in proteinogenic amino acids was then analyzed using a TBDMS method .
Isotopically nonstationary MFA
MFA is estimated based on isotope labeling dynamics of the free metabolites. The MFA model included the Calvin cycle, pentose phosphate pathway, TCA cycle, the glyoxylate shunt, anaplerotic pathways, and photorespiration pathway. A list of reactions in the network model with their atom transitions is provided in Additional file 1: Table S1. A lumped biomass equation based on biomass composition analysis and biochemical equations was used to fit the model to biomass production (Additional file 1: Table S2). The INCA platform  used a custom MATLAB ODE solver to fit over 890, 600, and 800 individual mass fragments for Synechococcus 2973 PBR, SF, and Synechococcus 7942 PBR to obtain fluxes in the assumed cyanobacterial network, respectively. INCA provided a goodness of fit via a Chi square statistical test as well as confidence intervals of all estimated parameters (Model formulation details are included in the Additional file 1: Part 1).
Construction of the zwf, pgl, and the succinate dehydrogenase deletion mutants
Three mutants were built, which served to independently validate the related aspects of metabolism described by the modeling process. The ∆ zwf mutant was constructed by insertional mutagenesis of the first enzyme of the OPPP, glucose-6-phosphate 1 dehydrogenase. The ∆ zwf mutant was characterized via growth rate (Fig. 1). The mutant activity was confirmed via IP-LC–MS/MS by the absence of the downstream product of its disrupted gene, 6-phosphogluconate (Additional file 1: Figure S2d–f). The ∆pgl and ∆sdh mutant was constructed by a cpf1 based CRISPR system was used to delete succinate dehydrogenase from start codon to stop codon from the chromosome (Additional file 1: Part 1).
- ADP glucose:
mass isotopomer distribution
- NAD+ :
nitcotinamide adenine dinucleotide
nitcotinamide adenine dinucleotide (reduced)
- NADP+ :
nitcotinamide adenine dinucleotide phosphate
nitcotinamide adenine dinucleotide phosphate (reduced)
oxidative pentose phosphate pathway
tricarboxylic acid cycle
glucose uridine diphosphoglucose
Microalgae Sayre R. The potential for carbon capture. Bioscience. 2010;60:722–7.
Yu J, Liberton M, Cliften PF, Head RD, Jacobs JM, Smith RD, Koppenaal DW, Brand JJ, Pakrasi HB. Synechococcus elongatus UTEX 2973, a fast growing cyanobacterial chassis for biosynthesis using light and CO2. Sci Rep. 2015;5:8132.
Young JD, Shastri AA, Stephanopoulos G, Morgan JA. Mapping photoautotrophic metabolism with isotopically nonstationary 13C flux analysis. Metab Eng. 2011;13:656–65.
Xiong W, Morgan JA, Ungerer J, Wang B, Maness PC, Yu J. The plasticity of cyanobacterial metabolism supports direct CO2 conversion to ethylene. Nat Plants. 2015;1:15053.
Yang C, Hua Q, Shimizu K. Energetics and carbon metabolism during growth of microalgal cells under photoautotrophic, mixotrophic and cyclic light-autotrophic/dark-heterotrophic conditions. Biochem Eng J. 2000;6:87–102.
You L, Berla B, He L, Pakrasi HB, Tang YJ. 13C-MFA delineates the photomixotrophic metabolism of Synechocystis sp. PCC 6803 under light-and carbon-sufficient conditions. Biotechnol J. 2014;9:684–92.
Xiong W, Brune D, Vermaas WFJ. The γ-aminobutyric acid shunt contributes to closing the tricarboxylic acid cycle in S ynechocystis sp. PCC 6803: the γ-aminobutyric acid shunt in Synechocystis. Mol Microbiol. 2014;93:786–96.
Young JD. INCA: a computational platform for isotopically non-stationary metabolic flux analysis. Bioinformatics. 2014;30:1333–5.
Yuan J, Bennett BD, Rabinowitz JD. Kinetic flux profiling for quantitation of cellular metabolic fluxes. Nat Protoc. 2008;3:1328–40.
Pramanik J, Keasling JD. Stoichiometric model of Escherichia coli metabolism: incorporation of growth-rate dependent biomass composition and mechanistic energy requirements. Biotechnol Bioeng. 1997;56:398–421.
Mueller TJ, Ungerer J, Pakrasi HB, Maranas CD. Identifying the metabolic differeneces of a fast-growth phenotype in Synechococcus UTEX 2973. Sci Rep. 2017;7:41569.
Bolten CJ, Kiefer P, Letisse F, Portais J-C, Wittmann C. Sampling for metabolome analysis of microorganisms. Anal Chem. 2007;79:3843–9.
Millard P, Massou S, Wittmann C, Portais J-C, Létisse F. Sampling of intracellular metabolites for stationary and non-stationary 13C metabolic flux analysis in Escherichia coli. Anal Biochem. 2014;465:38–49.
Jazmin LJ, Xu Y, Cheah YE, Adebiyi AO, Johnson CH, Young JD. Isotopically nonstationary 13C flux analysis of cyanobacterial isobutyraldehyde production. Metab Eng. 2017;42:9–18.
Ma F, Jazmin LJ, Young JD, Allen DK. Isotopically nonstationary 13C flux analysis of changes in Arabidopsis thaliana leaf metabolism due to high light acclimation. Proc Natl Acad Sci. 2014;111:16967–72.
Rubin BE, Wetmore KM, Price MN, Diamond S, Shultzaberger RK, Lowe LC, Curtin G, Arkin AP, Deutschbauer A, Golden SS. The essential gene set of a photosynthetic organism. Proc Natl Acad Sci. 2015;112:E6634–43.
Bennett BD, Kimball EH, Gao M, Osterhout R, van Dien SJ, Rabinowitz JD. Absolute metabolite concentrations and implied enzyme active site occupancy in Escherichia coli. Nat Chem Biol. 2009;5:593–9.
Yan R, Zhu D, Zhang Z, Zeng Q, Chu J. Carbon metabolism and energy conversion of Synechococcus sp. PCC 7942 under mixotrophic conditions: comparison with photoautotrophic condition. J Appl Phycol. 2011;24:657–68.
Bricker TM, Zhang S, Laborde SM, Mayer PR, Frankel LK, Moroney JV. The Malic enzyme is required for optimal photoautotrophic growth of Synechocystis sp. strain PCC 6803 under continuous light but not under a diurnal light regimen. J Bacteriol. 2004;186:8144–8.
Oliver JWK, Atsumi S. A carbon sink pathway increases carbon productivity in cyanobacteria. Metab Eng. 2015;29:106–12.
Angermayr SA, van der Woude AD, Correddu D, Vreugdenhil A, Verrone V, Hellingwerf KJ. Exploring metabolic engineering design principles for the photosynthetic production of lactic acid by Synechocystis sp. PCC 6803. Biotechnol Biofuels. 2014;7:99.
Nogales J, Gudmundsson S, Knight EM, Palsson BO, Thiele I. Detailing the optimality of photosynthesis in cyanobacteria through systems biology analysis. Proc Natl Acad Sci. 2012;109:2678–83.
Kramer DM, Evans JR. The importance of energy balance in improving photosynthetic productivity. Plant Physiol. 2011;155:70–8.
Buchanan BB. Role of light in the regulation of chloroplast enzymes. Annu Rev Plant Physiol. 1980;31:341–74.
Werdan K, Heldt HW, Milovancev M. The role of pH in the regulation of carbon fixation in the chloroplast stroma. Studies on CO2 fixation in the light and dark. Biochim Biophys Acta. 1975;396:276–92.
Scanlan DJ, Sundaram S, Newman J, Mann NH, Carr NG. Characterization of a zwf mutant of Synechococcus sp. strain PCC 7942. J Bacteriol. 1995;177:2550–3.
Schwender J, Hebbelmann I, Heinzel N, Hildebrandt T, Rogers A, Naik D, Klapperstuck M, Braun HP, Schreiber F, Denolf P, et al. Quantitative multilevel analysis of central metabolism in developing oilseeds of oilseed rape during in vitro culture. Plant Physiol. 2015;168:828–48.
Buescher JM, Antoniewicz MR, Boros LG, Burgess SC, Brunengraber H, Clish CB, DeBerardinis RJ, Feron O, Frezza C, Ghesquiere B, et al. A roadmap for interpreting 13C metabolite labeling patterns from cells. Curr Opin Biotechnol. 2015;34:189–201.
Huege J, Goetze J, Schwarz D, Bauwe H, Hagemann M, Kopka J. Modulation of the major paths of carbon in photorespiratory mutants of Synechocystis. PLoS ONE. 2011;6:e16278.
You C, Okano H, Hui S, Zhang Z, Kim M, Gunderson CW, Wang YP, Lenz P, Yan D, Hwa T. Coordination of bacterial proteome with metabolism by cyclic AMP signalling. Nature. 2013;500:301–6.
Liang F, Lindblad P. Effects of overexpressing photosynthetic carbon flux control enzymes in the cyanobacterium Synechocystis PCC 6803. Metab Eng. 2016;38:56–64.
Uematsu K, Suzuki N, Iwamae T, Inui M, Yukawa H. Increased fructose 1,6-bisphosphate aldolase in plastids enhances growth and photosynthesis of tobacco plants. J Exp Bot. 2012;63:3001–9.
Bernstein HC, McClure RS, Hill EA, Markillie LM, Chrisler WB, Romine MF, McDermott JE, Posewitz MC, Bryant DA, Konopka AE, et al. Unlocking the constraints of cyanobacterial productivity: acclimations enabling ultrafast growth. mBio. 2016;7:e00949–1016.
Markou G, Nerantzis E. Microalgae for high-value compounds and biofuels production: a review with focus on cultivation under stress conditions. Biotechnol Adv. 2013;31:1532–42.
Srivastava AK, Rai AN, Neilan BA. Stress biology of cyanobacteria: molecular mechanisms to cellular responses. Boca Raton: CRC Press; 2013.
Carrieri D, Broadbent C, Carruth D, Paddock T, Ungerer J, Maness PC, Ghirardi M, Yu J. Enhancing photo-catalytic production of organic acids in the cyanobacterium Synechocystis sp. PCC 6803 ΔglgC, a strain incapable of glycogen storage. Microb Biotechnol. 2015;8:275–80.
Zhou J, Zhang F, Meng H, Zhang Y, Li Y. Introducing extra NADPH consumption ability significantly increases the photosynthetic efficiency and biomass production of cyanobacteria. Metab Eng. 2016;38:217–27.
Iijima H, Nakaya Y, Kuwahara A, Hirai MY, Osanai T. Seawater cultivation of freshwater cyanobacterium Synechocystis sp. PCC 6803 drastically alters amino acid composition and glycogen metabolism. Front Microbiol. 2015;6:326.
Nielsen J, Keasling JD. Engineering cellular metabolism. Cell. 2016;164:1185–97.
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.
Abed RM, Dobretsov S, Sudesh K. Applications of cyanobacteria in biotechnology. J Appl Microbiol. 2009;106:1–12.
Mohr W, Vagner T, Kuypers MMM, Ackermann M, LaRoche J. Resolution of conflicting signals at the single-cell level in the regulation of cyanobacterial photosynthesis and nitrogen fixation. PLoS ONE. 2013;8:e66060.
Nargund S, Misra A, Zhang X, Coleman GD, Sriram G. Flux and reflux: metabolite reflux in plant suspension cells and its implications for isotope-assisted metabolic flux analysis. Mol BioSyst. 2014;10:1496–508.
He L, Wu SG, Wan N, Reding AC, Tang YJ. Simulating cyanobacterial phenotypes by integrating flux balance analysis, kinetics, and a light distribution function. Microb Cell Factories. 2015;14:206.
Wheeldon I, Minteer SD, Banta S, Barton SC, Atanassov P, Sigman M. Substrate channelling as an approach to cascade reactions. Nat Chem. 2016;8:299–309.
Conrado RJ, Varner JD, DeLisa MP. Engineering the spatial organization of metabolic enzymes: mimicking nature’s synergy. Curr Opin Biotechnol. 2008;19:492–9.
Jandt U, You C, Zhang YHP, Zeng AP. Compartmentalization and metabolic channeling for multienzymatic biosynthesis: practical strategies and modeling approaches. In: Zeng AP, editor. Fundamentals and application of new bioproduction systems, advances in biochemical engineering/biotechnology. Berlin: Springer; 2013. p. 41–65.
Shearer G, Lee JC, Koo J, Kohl DH. Quantitative estimation of channeling from early glycolytic intermediates to CO2 in intact Escherichia coli. FEBS J. 2005;272:3260–9.
Allen DK, Laclair RW, Ohlrogge JB, Shachar-Hill Y. Isotope labelling of Rubisco subunits provides in vivo information on subcellular biosynthesis and exchange of amino acids between compartments. Plant Cell Environ. 2012;35:1232–44.
Mandy DE, Goldford JE, Yang H, Allen DK, Libourel IGL. Metabolic flux analysis using 13C peptide label measurements. Plant J. 2014;77:476–86.
Allen DK, Goldford J, Gierse JK, Mandy D, Diepenbrock C, Libourel IG. Quantification of peptide m/z distributions from 13C-labeled cultures with high-resolution mass spectrometry. Anal Chem. 2014;86:1894–901.
Allen DK, Evans BS, Libourel IG. Analysis of isotopic labeling in peptide fragments by tandem mass spectrometry. PLoS ONE. 2014;9:e91537.
Hollinshead WD, Rodriguez S, Martin HG, Wang G, Baidoo EEK, Sale KL, Keasling JD, Mukhopadhyay A, Tang YJ. Examining Escherichia coli glycolytic pathways, catabolite repression, and metabolite channeling using ∆pfk mutants. Biotechnol Biofuels. 2016;9:212.
Bennett BD, Yuan J, Kimball EH, Rabinowitz JD. Absolute quantitation of intracellular metabolite concentrations by an isotope ratio-based approach. Nat Protoc. 2008;38:1299–311.
Wu L, Mashego MR, van Dam JC, Proell AM, Vinke JL, Ras C, van Winden WA, van Gulik WM, Heijnen JJ. Quantitative analysis of the microbial metabolome by isotope dilution mass spectrometry using uniformly 13C-labeled cell extracts as internal standards. Anal Biochem. 2005;336:164–71.
You L, Page L, Feng X, Berla B, Pakrasi HB, Tang YJ. Metabolic pathway confirmation and discovery through 13C-labeling of proteinogenic amino acids. J Vis Exp. 2012;59:e3583.
HBP, CDM, and YJT initiated the project. MHA, FM, YJT and DKA designed and performed 13C-MFA and LC–MS analysis. JY, HBP, MHA and WH performed growth experiments and physiological analysis. HBP, ML, JU, MHA constructed mutants and helped with strain characterizations. CDM, SG and LH assisted the flux modeling and statistics. All authors read and approved the final manuscript.
We are thankful for Dr. Edward Baidoo from Joint Bio-Energy Institute for advices on LC–MS analysis and on Dr. Wei Liao at Michigan State University for biomass analysis.
The authors declare that they have no competing interests.
Availability of data and materials
All datasets used and/or analyzed during the current study are available in this published article, its supplementary information files, or is available from the corresponding author on reasonable request.
Consent for publication
Ethical approval and consent to participate
This work was supported by DOE-BER (DESC0012722) to CDM, Y.J.T. and H.B.P. and NSF (EF-1105249), DOE (AE-AR0000202) and USDA-ARS funding to support F.M. and D.K.A. Analyses were in part enabled by major research instrumentation grants from the NSF (DBI-1427621 and DBI-0521250) that supported the acquisition of mass spectrometers within the Proteomics and Mass Spectrometry Core Facility at the Donald Danforth Plant Science Center. The core facility is additionally acknowledged for help with metabolite evaluations.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.