Metabolic flux analysis of Cyanothece sp. ATCC 51142 under mixotrophic conditions

Precision fermentation with mass spectrometry-based spent media analysis

Optimization and monitoring of bioprocesses requires the measurement of several process parameters and quality attributes. Mass spectrometry (MS)-based techniques such as those coupled to gas chromatography (GCMS) and liquid Chromatography (LCMS) enable the simultaneous measurement of hundreds of metabolites with high sensitivity. When applied to spent media, such metabolome analysis can help determine the sequence of substrate uptake and metabolite secretion, consequently facilitating better design of initial media and feeding strategy. Furthermore, the analysis of metabolite diversity and abundance from spent media will aid the determination of metabolic phases of the culture and the identification of metabolites as surrogate markers for product titer and quality. This review covers the recent advances in metabolomics analysis applied to the development and monitoring of bioprocesses. In this regard, we recommend a stepwise workflow and guidelines that a bioprocesses engineer can adopt to develop and optimize a fermentation process using spent media analysis. Finally, we show examples of how the use of MS can revolutionize the design and monitoring of bioprocesses.

Quantitative insight into the metabolism of isoprene-producing Synechocystis sp. PCC 6803 using steady state 13C-MFA

Cyanobacteria are photosynthetic bacteria, widely studied for the conversion of atmospheric carbon dioxide to useful platform chemicals. Isoprene is one such industrially important chemical, primarily used for production of synthetic rubber and biofuels. Synechocystis sp. PCC 6803, a genetically amenable cyanobacterium, produces isoprene on heterologous expression of isoprene synthase gene, albeit in very low quantities. Rationalized metabolic engineering to re-route the carbon flux for enhanced isoprene production requires in-dept knowledge of the metabolic flux distribution in the cell. Hence, in the present study, we undertook steady state ¹³C-metabolic flux analysis of glucose-tolerant wild-type (GTN) and isoprene-producing recombinant (ISP) Synechocystis sp. to understand and compare the carbon flux distribution in the two strains. The R-values for amino acids, flux analysis predictions and gene expression profiles emphasized predominance of Calvin cycle and glycogen metabolism in GTN. Alternatively, flux analysis predicted higher activity of the anaplerotic pathway through phosphoenolpyruvate carboxylase and malic enzyme in ISP. The striking difference in the Calvin cycle, glycogen metabolism and anaplerotic pathway activity in GTN and ISP suggested a possible role of energy molecules (ATP and NADPH) in regulating the carbon flux distribution in GTN and ISP. This claim was further supported by the transcript level of selected genes of the electron transport chain. This study provides the first quantitative insight into the carbon flux distribution of isoprene-producing cyanobacterium, information critical for developing Synechocystis sp. as a single cell factory for isoprenoid/terpenoid production.

Multi-Objective Optimization of Microalgae Metabolism: An Evolutive Algorithm Based on FBA

Studies enabled by metabolic models of different species of microalgae have become significant since they allow us to understand changes in their metabolism and physiological stages. The most used method to study cell metabolism is FBA, which commonly focuses on optimizing a single objective function. However, recent studies have brought attention to the exploration of simultaneous optimization of multiple objectives. Such strategies have found application in optimizing biomass and several other bioproducts of interest; they usually use approaches such as multi-level models or enumerations schemes. This work proposes an alternative in silico multiobjective model based on an evolutionary algorithm that offers a broader approximation of the Pareto frontier, allowing a better angle for decision making in metabolic engineering. The proposed strategy is validated on a reduced metabolic network of the microalgae Chlamydomonas reinhardtii while optimizing for the production of protein, carbohydrates, and CO2 uptake. The results from the conducted experimental design show a favorable difference in the number of solutions achieved compared to a classic tool solving FBA.

Metabolic systems biology and multi-omics of cyanobacteria: Perspectives and future directions

Cyanobacteria are oxygenic photoautotrophs whose metabolism contains key biochemical pathways to fix atmospheric CO₂ and synthesize various metabolites. The development of bioengineering tools has enabled the manipulation of cyanobacterial chassis to produce various valuable bioproducts photosynthetically. However, effective utilization of cyanobacteria as photosynthetic cell factories needs a detailed understanding of their metabolism and its interaction with other cellular processes. Implementing systems and synthetic biology tools has generated a wealth of information on various metabolic pathways. However, to design effective engineering strategies for further improvement in growth, photosynthetic efficiency, and enhanced production of target biochemicals, in-depth knowledge of their carbon/nitrogen metabolism, pathway fluxe distribution, genetic regulation and integrative analyses are necessary. In this review, we discuss the recent advances in the development of genome-scale metabolic models (GSMMs), omics analyses (metabolomics, transcriptomics, proteomics, fluxomics), and integrative modeling approaches to showcase the current understanding of cyanobacterial metabolism.

Robustifying Experimental Tracer Design for13C-Metabolic Flux Analysis

¹³C metabolic flux analysis (MFA) has become an indispensable tool to measure metabolic reaction rates (fluxes) in living organisms, having an increasingly diverse range of applications. Here, the choice of the¹³C labeled tracer composition makes the difference between an information-rich experiment and an experiment with only limited insights. To improve the chances for an informative labeling experiment, optimal experimental design approaches have been devised for¹³C-MFA, all relying on some a priori knowledge about the actual fluxes. If such prior knowledge is unavailable, e.g., for research organisms and producer strains, existing methods are left with a chicken-and-egg problem. In this work, we present a general computational method, termed robustified experimental design (R-ED), to guide the decision making about suitable tracer choices when prior knowledge about the fluxes is lacking. Instead of focusing on one mixture, optimal for specific flux values, we pursue a sampling based approach and introduce a new design criterion, which characterizes the extent to which mixtures are informative in view of all possible flux values. The R-ED workflow enables the exploration of suitable tracer mixtures and provides full flexibility to trade off information and cost metrics. The potential of the R-ED workflow is showcased by applying the approach to the industrially relevant antibiotic producer Streptomyces clavuligerus, where we suggest informative, yet economic labeling strategies.

The Effect of Promoter and RBS Combination on the Growth and Glycogen Productivity of Sodium-Dependent Bicarbonate Transporter (SbtA) Overexpressing Synechococcus sp. PCC 7002 Cells

Sodium dependent bicarbonate transporter, SbtA is a high-affinity, inducible bicarbonate transporter in cyanobacterial cells. Our previous work has shown that overexpression of this transporter can significantly increase growth and glycogen accumulation in Synechococcus sp. PCC 7002 cells. In this work, we have tested the effect of two different RBS sequences (RBS1: GGAGGA and RBS2: AGGAGA) and three different promoters (P_cpcB, P_cpcB560, and P_rbcL2) on the growth and glycogen production in SbtA-overexpressing Synechococcus sp. PCC 7002 cells. Our results show that P_cpcB or P_cpcB560 were more effective than P_rbcL2 in increasing the growth and glycogen content. The choice of RBS sequence had relatively minor effect, though RBS2 was more effective than RBS1. The transformant E, with P_cpcB560 and RBS2, showed the highest growth. The biomass after 5 days of growth on air or 1% CO₂ was increased by about 90% in the strain E compared to PCC 7002 cells. All transformants overexpressing SbtA had higher glycogen content. However, growing the cells with bubbling of 1% CO₂ did not increase cellular glycogen content any further. The strain E had about 80% higher glycogen content compared to WT PCC 7002 cells. Therefore, the glycogen productivity of the strain E grown with air-bubbling was about 2.5-fold that of the WT PCC 7002 cells grown similarly. Additionally, some of the transformants had higher chlorophyll content while all the transformants had higher carotenoid content compared to the PCC 7002 cells, suggesting interaction between carbon transport and pigment levels. Thus, this work shows that the choice of photosynthetic promoters and RBSs sequences can impact growth and glycogen accumulation in SbtA-overexpressing cells.

Liquid Chromatography Methods for Separation of Polar and Charged Intracellular Metabolites for 13C Metabolic Flux Analysis

Accurate quantification of mass isotopolog distribution (MID) of intracellular metabolites is a key requirement for ¹³C metabolic flux analysis (¹³C-MFA). Liquid chromatography coupled with mass spectrometry (LC/MS) has emerged as a frontrunner technique that combines two orthogonal separation strategies. While metabolomics requires separation of monoisotopic peaks, ¹³C-MFA imposes additional demands for chromatographic separation as isotopologs of metabolites significantly add to the number of analytes. In this protocol chapter, we discuss two liquid chromatography methods, namely, reverse phase ion-pairing and hydrophilic interaction chromatography (HILIC) that together can separate a wide variety of metabolites that are typically used for ¹³C metabolic flux analysis.

Glycolytic Shunts Replenish the Calvin-Benson-Bassham Cycle as Anaplerotic Reactions in Cyanobacteria

The recent discovery of the Entner-Doudoroff (ED) pathway as a third glycolytic route beside Embden-Meyerhof-Parnas (EMP) and oxidative pentose phosphate (OPP) pathway in oxygenic photoautotrophs requires a revision of their central carbohydrate metabolism. In this study, unexpectedly, we observed that deletion of the ED pathway alone, and even more pronounced in combination with other glycolytic routes, diminished photoautotrophic growth in continuous light in the cyanobacterium Synechocystis sp. PCC 6803. Furthermore, we found that the ED pathway is required for optimal glycogen catabolism in parallel to an operating Calvin-Benson-Bassham (CBB) cycle. It is counter-intuitive that glycolytic routes, which are a reverse to the CBB cycle and do not provide any additional biosynthetic intermediates, are important under photoautotrophic conditions. However, observations on the ability to reactivate an arrested CBB cycle revealed that they form glycolytic shunts that tap the cellular carbohydrate reservoir to replenish the cycle. Taken together, our results suggest that the classical view of the CBB cycle as an autocatalytic, completely autonomous cycle that exclusively relies on its own enzymes and CO2 fixation to regenerate ribulose-1,5-bisphosphate for Rubisco is an oversimplification. We propose that in common with other known autocatalytic cycles, the CBB cycle likewise relies on anaplerotic reactions to compensate for the depletion of intermediates, particularly in transition states and under fluctuating light conditions that are common in nature.

Metabolic engineering of a fast-growing cyanobacterium Synechococcus elongatus PCC 11801 for photoautotrophic production of succinic acid

BACKGROUND

Cyanobacteria, a group of photosynthetic prokaryotes, are being increasingly explored for direct conversion of carbon dioxide to useful chemicals. However, efforts to engineer these photoautotrophs have resulted in low product titers. This may be ascribed to the bottlenecks in metabolic pathways, which need to be identified for rational engineering. We engineered the recently reported, fast-growing and robust cyanobacterium, Synechococcus elongatus PCC 11801 to produce succinate, an important platform chemical. Previously, engineering of the model cyanobacterium S. elongatus PCC 7942 has resulted in succinate titer of 0.43 g l^-1 in 8 days.

RESULTS

Building on the previous report, expression of α-ketoglutarate decarboxylase, succinate semialdehyde dehydrogenase and phosphoenolpyruvate carboxylase yielded a succinate titer of 0.6 g l^-1 in 5 days suggesting that PCC 11801 is better suited as host for production. Profiling of the engineered strains for 57 intermediate metabolites, a number of enzymes and qualitative analysis of key transcripts revealed potential flux control points. Based on this, we evaluated the effects of overexpression of sedoheptulose-1,7-bisphosphatase, citrate synthase and succinate transporters and knockout of succinate dehydrogenase and glycogen synthase A. The final construct with seven genes overexpressed and two genes knocked out resulted in photoautotrophic production of 0.93 g l^-1 succinate in 5 days.

CONCLUSION

While the fast-growing strain PCC 11801 yielded a much higher titer than the model strain, the efficient photoautotrophy of this novel isolate needs to be harnessed further for the production of desired chemicals. Engineered strains of S. elongatus PCC 11801 showed dramatic alterations in the levels of several metabolites suggesting far reaching effects of pathway engineering. Attempts to overexpress enzymes deemed to be flux controlling led to the emergence of other potential rate-limiting steps. Thus, this process of debottlenecking of the pathway needs to be repeated several times to obtain a significantly superior succinate titer.

Evaluation of freely available software tools for untargeted quantification of 13C isotopic enrichment in cellular metabolome from HR-LC/MS data

¹³C Metabolic Flux Analysis (¹³C-MFA) involves the quantification of isotopic enrichment in cellular metabolites and fitting the resultant data to the metabolic network model of the organism. Coverage and resolution of the resultant flux map depends on the total number of metabolites and fragments in which ¹³C enrichment can be quantified accurately. Experimental techniques for tracking ¹³C enrichment are evolving rapidly and large volumes of data are now routinely generated through the use of Liquid Chromatography coupled with High-Resolution Mass Spectrometry (HR-LC/MS). Therefore, the current manuscript is focused on the challenges in high-throughput analyses of such large datasets. Current ¹³C-MFA studies often have to rely on the targeted quantification of a small subset of metabolites, thereby leaving a large fraction of the data unexplored. A number of public domain software tools have been reported in recent years for the untargeted quantitation of isotopic enrichment. However, the suitability of their application across diverse datasets has not been investigated. Here, we test the software tools X¹³CMS, DynaMet, geoRge, and HiResTEC with three diverse datasets. The tools provided a global, untargeted view of ¹³C enrichment in metabolites in all three datasets and a much-needed automation in data analysis. Some inconsistencies were observed in results obtained from the different tools, which could be partially ascribed to the lack of baseline separation and potential mass conflicts. After removing the false positives manually, isotopic enrichment could be quantified reliably in a large repertoire of metabolites. Of the software tools explored, geoRge and HiResTEC consistently performed well for the untargeted analysis of all datasets tested.

Mass Isotopologue Distribution of dimer ion adducts of intracellular metabolites for potential applications in 13C Metabolic Flux Analysis

¹³C Metabolic Flux Analysis (¹³C-MFA) is a powerful tool for quantification of carbon flux distribution in metabolic pathways. However, the requirement to obtain accurate labeling patterns, especially for compounds with low abundance, poses a challenge. Chromatographic separation and high sensitivity of the modern mass spectrometers (MS) alleviate this problem to a certain extent. However, the presence of derivatives such as in-source fragments, multimer ion adducts, and multiply charged ions result in reduced intensity of the molecular ion. While multimer ion adducts have been reported in the field of metabolomics, their presence is considered undesirable in quantitative studies. Here, we demonstrate a novel application of dimer ion adducts in calculating the mass isotopologue distribution (MIDs) of the corresponding monomer ions for public domain and in-house generated datasets comprising of ¹³C-labeling time-course experiments. Out of the 100 standard compounds analyzed, we could detect multimer ion adducts in 24 of the intermediate metabolites. Further, a subset of these multimer ions were detected in all the biological samples analyzed. Majority of these ion adducts were either not detected in the original study or labeled as a putative features. Regression analysis was performed to estimate the monomer MIDs from those of the dimer. This resulted in accurate estimation regardless of the biological system, chromatographic method, the MS hardware, or the relative abundance of the dimer ion. We argue that this analysis may be useful in cases where satisfactory data cannot be extracted from the chromatographic peaks of the monomer ions.

The effect of CO2 in enhancing photosynthetic cofactor recycling for alcohol dehydrogenase mediated chiral synthesis in cyanobacteria

The light harvesting photosystem in cyanobacteria offers a potential pathway for the regeneration of the nicotinamide cofactor NADPH, thereby facilitating the application of cyanobacteria as excellent whole cell biocatalysts in oxidoreductase-mediated biotransformation. The use of cyanobacterial metabolism for cofactor recycling improves the atom economy of the process compared to the commonly employed enzyme-coupled cofactor recycling using enzymes such as glucose dehydrogenase. Here we report the asymmetric conversion of acetophenone to chiral 1-phenylethanol by recombinant Synechococcus elongatus PCC 7942 whole cell biocatalyst that expresses the NADPH dependent L. kefir alcohol dehydrogenase. Besides light, it was observed that carbon dioxide levels play a critical role in improving the bioconversion efficiency possibly due to the enhanced growth rate and improved cofactor availability at elevated CO₂ levels. Complete reduction of acetophenone to optically pure (R)-1-phenylethanol at 99% enantiomeric excess was achieved within 6 h with a relatively low cell density of 0.66 g/l by coupling optimum light and CO₂ levels and without the need for a co-substrate.

An improved method for extraction of polar and charged metabolites from cyanobacteria

A key requirement for 13C Metabolic flux analysis (13C-MFA), a widely used technique to estimate intracellular metabolic fluxes, is an efficient method for the extraction of intermediate metabolites for analysis via liquid chromatography mass spectrometry (LC/MS). The 13C isotopic labeling results in further distribution of an already sparse pool of intermediate metabolites into isotopologues, each appearing as a separate chromatographic feature. We examined some of the reported solvent systems for the extraction of polar intracellular metabolites from three strains of cyanobacteria of the genus Synechococcus, viz., Synechococcus sp. PCC 7002, Synechococcus elongatus PCC 7942, and a newly isolated Synechococcus elongatus PCC 11801 (manuscript under review). High resolution-LC/MS was used to assess the relative abundance of the extracted metabolites. The different solvent systems used for extraction led to statistically significant changes in the extraction efficiency for a large number of metabolites. While a few hundred m/z features or potential metabolites were detected with different solvent systems, the abundance of over a quarter of all metabolites varied significantly from one solvent system to another. Further, the extraction methods were evaluated for a targeted set of metabolites that are important in 13C-MFA studies of photosynthetic organisms. While for the strain PCC 7002, the reported method using methanol-chloroform-water system gave satisfactory results, a mild base in the form of NH4OH had to be used in place of water to achieve adequate levels of extraction for PCC 7942 and PCC 11801. While minor changes in extraction solvent resulted in dramatic changes in the extraction efficiency of a number of compounds, certain metabolites such as amino acids and organic acids were adequately extracted in all the solvent systems tested. Overall, we present a new improved method for extraction using a methanol-chloroform-NH4OH system. Our method improves the extraction of polar compounds such as sugar phosphates, bisphosphates, that are central to 13C-MFA studies.

Rerouting of Metabolism into Desired Cellular Products by Nutrient Stress: Fluxes Reveal the Selected Pathways in Cyanobacterial Photosynthesis

Boosting cellular growth rates while redirecting metabolism to make desired products are the preeminent goals of gene engineering of photoautotrophs, yet so far these goals have been hardly achieved owing to lack of understanding of the functional pathways and their choke points. Here we apply a ¹³C mass isotopic method (INST-MFA) to quantify instantaneous fluxes of metabolites during photoautotrophic growth. INST-MFA determines the globally most accurate set of absolute fluxes for each metabolite from a finite set of measured ¹³C-isotopomer fluxes by minimizing the sum of squared residuals between experimental and predicted mass isotopomers. We show that the widely observed shift in biomass composition in cyanobacteria, demonstrated here with Synechococcus sp. PCC 7002, favoring glycogen synthesis during nitrogen starvation is caused by (1) increased flux through a bottleneck step in gluconeogenesis (3PG → GAP/DHAP), and (2) flux overflow through a previously unrecognized hybrid gluconeogenesis-pentose phosphate (hGPP) pathway. Our data suggest the slower growth rate and biomass accumulation under N starvation is due to a reduced carbon fixation rate and a reduced flux of carbon into amino acid precursors. Additionally, ¹³C flux from α-ketoglutarate to succinate is demonstrated to occur via succinic semialdehyde, an alternative to the conventional TCA cycle, in Synechococcus 7002 under photoautotrophic conditions. We found that pyruvate and oxaloacetate are synthesized mainly by malate dehydrogenase with minimal flux into acetyl coenzyme-A via pyruvate dehydrogenase. Nutrient stress induces major shifts in fluxes into new pathways that deviate from historical metabolic pathways derived from model bacteria.

SWATH Tandem Mass Spectrometry Workflow for Quantification of Mass Isotopologue Distribution of Intracellular Metabolites and Fragments Labeled with Isotopic 13C Carbon

Accurate quantification of mass isotopologue distribution (MID) of metabolites is a prerequisite for ¹³C-metabolic flux analysis. Currently used mass spectrometric (MS) techniques based on multiple reaction monitoring (MRM) place limitations on the number of MIDs that can be analyzed in a single run. Moreover, the deconvolution step results in amplification of error. Here, we demonstrate that SWATH MS/MS, a data independent acquisition (DIA) technique allows quantification of a large number of precursor and product MIDs in a single run. SWATH sequentially fragments all precursor ions in stacked mass isolation windows. Co-fragmentation of all precursor isotopologues in a single SWATH window yields higher sensitivity enabling quantification of MIDs of fragments with low abundance and lower systematic and random errors. We quantify the MIDs of 53 precursor and product ions corresponding to 19 intracellular metabolites from a dynamic ¹³C-labeling of a model cyanobacterium, Synechococcus sp. PCC 7002. The use of product MIDs resulted in an improved precision of many measured fluxes compared to when only precursor MIDs were used for flux analysis. The approach is truly untargeted and allows additional metabolites to be quantified from the same data.

Recent advances in synthetic biology of cyanobacteria

Cyanobacteria are attractive hosts that can be engineered for the photosynthetic production of fuels, fine chemicals, and proteins from CO₂. Moreover, the responsiveness of these photoautotrophs towards different environmental signals, such as light, CO₂, diurnal cycle, and metals make them potential hosts for the development of biosensors. However, engineering these hosts proves to be a challenging and lengthy process. Synthetic biology can make the process of biological engineering more predictable through the use of standardized biological parts that are well characterized and tools to assemble them. While significant progress has been made with model heterotrophic organisms, many of the parts and tools are not portable in cyanobacteria. Therefore, efforts are underway to develop and characterize parts derived from cyanobacteria. In this review, we discuss the reported parts and tools with the objective to develop cyanobacteria as cell factories or biosensors. We also discuss the issues related to characterization, tunability, portability, and the need to develop enabling technologies to engineer this "green" chassis.

Rewiring of Cyanobacterial Metabolism for Hydrogen Production: Synthetic Biology Approaches and Challenges

With the demand for renewable energy growing, hydrogen (H₂) is becoming an attractive energy carrier. Developing H₂ production technologies with near-net zero carbon emissions is a major challenge for the "H₂ economy." Certain cyanobacteria inherently possess enzymes, nitrogenases, and bidirectional hydrogenases that are capable of H₂ evolution using sunlight, making them ideal cell factories for photocatalytic conversion of water to H₂. With the advances in synthetic biology, cyanobacteria are currently being developed as a "plug and play" chassis to produce H₂. This chapter describes the metabolic pathways involved and the theoretical limits to cyanobacterial H₂ production and summarizes the metabolic engineering technologies pursued.

13C-assisted metabolic flux analysis to investigate heterotrophic and mixotrophic metabolism in Cupriavidus necator H16

INTRODUCTION

Cupriavidus necator H16 is a gram-negative bacterium, capable of lithoautotrophic growth by utilizing hydrogen as an energy source and fixing carbon dioxide (CO2) through Calvin-Benson-Bassham (CBB) cycle. The potential to utilize synthesis gas (Syngas) and the prospects of rerouting carbon from polyhydroxybutyrate synthesis to value-added compounds makes C. necator an excellent chassis for industrial application.

OBJECTIVES

In the context of lack of sufficient quantitative information of the metabolic pathways and to advance in rational metabolic engineering for optimized product synthesis in C. necator H16, we carried out a metabolic flux analysis based on steady-state 13C-labelling.

METHODS

In this study, steady-state carbon labelling experiments, using either d-[1-13C]fructose or [1,2-13C]glycerol, were undertaken to investigate the carbon flux through the central carbon metabolism in C. necator H16 under heterotrophic and mixotrophic growth conditions, respectively.

RESULTS

We found that the CBB cycle is active even under heterotrophic condition, and growth is indeed mixotrophic. While Entner-Doudoroff (ED) pathway is shown to be the major route for sugar degradation, tricarboxylic acid (TCA) cycle is highly active in mixotrophic condition. Enhanced flux is observed in reductive pentose phosphate pathway (redPPP) under the mixotrophic condition to supplement the precursor requirement for CBB cycle. The flux distribution was compared to the mRNA abundance of genes encoding enzymes involved in key enzymatic reactions of the central carbon metabolism.

CONCLUSION

This study leads the way to establishing 13C-based quantitative fluxomics for rational pathway engineering in C. necator H16.

Rerouting of carbon flux in a glycogen mutant of cyanobacteria assessed via isotopically non-stationary 13 C metabolic flux analysis

Cyanobacteria, which constitute a quantitatively dominant phylum, have attracted attention in biofuel applications due to favorable physiological characteristics, high photosynthetic efficiency and amenability to genetic manipulations. However, quantitative aspects of cyanobacterial metabolism have received limited attention. In the present study, we have performed isotopically non-stationary ¹³ C metabolic flux analysis (INST-¹³ C-MFA) to analyze rerouting of carbon in a glycogen synthase deficient mutant strain (glgA-I glgA-II) of the model cyanobacterium Synechococcus sp. PCC 7002. During balanced photoautotrophic growth, 10-20% of the fixed carbon is stored in the form of glycogen via a pathway that is conserved across the cyanobacterial phylum. Our results show that deletion of glycogen synthase gene orchestrates cascading effects on carbon distribution in various parts of the metabolic network. Carbon that was originally destined to be incorporated into glycogen gets partially diverted toward alternate storage molecules such as glucosylglycerol and sucrose. The rest is partitioned within the metabolic network, primarily via glycolysis and tricarboxylic acid cycle. A lowered flux toward carbohydrate synthesis and an altered distribution at the glucose-1-phosphate node indicate flexibility in the network. Further, reversibility of glycogen biosynthesis reactions points toward the presence of futile cycles. Similar redistribution of carbon was also predicted by Flux Balance Analysis. The results are significant to metabolic engineering efforts with cyanobacteria where fixed carbon needs to be re-routed to products of interest. Biotechnol. Bioeng. 2017;114: 2298-2308. © 2017 Wiley Periodicals, Inc.

Gene essentiality, conservation index and co-evolution of genes in cyanobacteria

Cyanobacteria, a group of photosynthetic prokaryotes, dominate the earth with ~ 1015 g wet biomass. Despite diversity in habitats and an ancient origin, cyanobacterial phylum has retained a significant core genome. Cyanobacteria are being explored for direct conversion of solar energy and carbon dioxide into biofuels. For this, efficient cyanobacterial strains will need to be designed via metabolic engineering. This will require identification of target knockouts to channelize the flow of carbon toward the product of interest while minimizing deletions of essential genes. We propose "Gene Conservation Index" (GCI) as a quick measure to predict gene essentiality in cyanobacteria. GCI is based on phylogenetic profile of a gene constructed with a reduced dataset of cyanobacterial genomes. GCI is the percentage of organism clusters in which the query gene is present in the reduced dataset. Of the 750 genes deemed to be essential in the experimental study on S. elongatus PCC 7942, we found 494 to be conserved across the phylum which largely comprise of the essential metabolic pathways. On the contrary, the conserved but non-essential genes broadly comprise of genes required under stress conditions. Exceptions to this rule include genes such as the glycogen synthesis and degradation enzymes, deoxyribose-phosphate aldolase (DERA), glucose-6-phosphate 1-dehydrogenase (zwf) and fructose-1,6-bisphosphatase class1, which are conserved but non-essential. While the essential genes are to be avoided during gene knockout studies as potentially lethal deletions, the non-essential but conserved set of genes could be interesting targets for metabolic engineering. Further, we identify clusters of co-evolving genes (CCG), which provide insights that may be useful in annotation. Principal component analysis (PCA) plots of the CCGs are demonstrated as data visualization tools that are complementary to the conventional heatmaps. Our dataset consists of phylogenetic profiles for 23,643 non-redundant cyanobacterial genes. We believe that the data and the analysis presented here will be a great resource to the scientific community interested in cyanobacteria.

Metabolic model of Synechococcus sp. PCC 7002: Prediction of flux distribution and network modification for enhanced biofuel production

Flux Balance Analysis was performed with the Genome Scale Metabolic Model of a fast growing cyanobacterium Synechococcus sp. PCC 7002 to gain insights that would help in engineering the organism as a production host. Gene essentiality and synthetic lethality analysis revealed a reduced metabolic robustness under genetic perturbation compared to the heterotrophic bacteria Escherichia coli. Under glycerol heterotrophy the reducing equivalents were generated from tricarboxylic acid cycle rather than the oxidative pentose phosphate pathway. During mixotrophic growth in glycerol the photosynthetic electron transport chain was predominantly used for ATP synthesis with a photosystem I/photosystem II flux ratio higher than that observed under autotrophy. An exhaustive analysis of all possible double reaction knock outs was performed to reroute fixed carbon towards ethanol and butanol production. It was predicted that only ∼10% of fixed carbon could be diverted for ethanol and butanol production.

13 C flux analysis of cyanobacterial metabolism

(13)C metabolic flux analysis (MFA) has made important contributions to our understanding of the physiology of model strains of E. coli and yeast, and it has been widely used to guide metabolic engineering efforts in these microorganisms. Recent advancements in (13)C MFA methodology combined with publicly available software tools are creating new opportunities to extend this approach to examine less characterized microbes. In particular, growing interest in the use of cyanobacteria as industrial hosts for photosynthetic production of biofuels and biochemicals has led to a critical need to better understand how cyanobacterial metabolic fluxes are regulated in response to changes in growth conditions or introduction of heterologous pathways. In this contribution, we review several prior studies that have applied isotopic steady-state (13)C MFA to examine heterotrophic or mixotrophic growth of cyanobacteria, as well as recent studies that have pioneered the use of isotopically nonstationary MFA (INST-MFA) to study autotrophic cultures. We also provide recommendations for the design and analysis of INST-MFA experiments in cyanobacteria, based on our previous experience and a series of simulation studies used to assess the selection of measurements and sample time points. We anticipate that this emerging knowledgebase of prior (13)C MFA studies, optimized experimental protocols, and public software tools will catalyze increasing use of (13)C MFA techniques by the cyanobacteria research community.

Challenges and opportunities for microalgae-mediated CO2 capture and biorefinery

Aquacultures of microalgae are frontrunners for photosynthetic capture of CO2 from flue gases. Expedient implementation mandates coupling of microalgal CO2 capture with synthesis of fuels and organic products, so as to derive value from biomass. An integrated biorefinery complex houses a biomass growth and harvesting area and a refining zone for conversion to product(s) and separation to desired purity levels. As growth and downstream options require energy and incur loss of carbon, put together, the loop must be energy positive, carbon negative, or add substantial value. Feasibility studies can, thus, aid the choice from among the rapidly evolving technological options, many of which are still in the early phases of development. We summarize basic engineering calculations for the key steps of a biorefining loop where flue gases from a thermal power station are captured using microalgal biomass along with subsequent options for conversion to fuel or value added products. An assimilation of findings from recent laboratory and pilot-scale experiments and life cycle analysis (LCA) studies is presented as carbon and energy yields for growth and harvesting of microalgal biomass and downstream options. Of the biorefining options, conversion to the widely studied biofuel, ethanol, and manufacture of the platform chemical, succinic acid are presented. Both processes yield specific products and do not demand high-energy input but entail 60-70% carbon loss through fermentative respiration. Thermochemical conversions, on the other hand, have smaller carbon and energy losses but yield a mixture of products.

Coupling of Cellular Processes and Their Coordinated Oscillations under Continuous Light in Cyanothece sp. ATCC 51142, a Diazotrophic Unicellular Cyanobacterium

Unicellular diazotrophic cyanobacteria such as Cyanothece sp. ATCC 51142 (henceforth Cyanothece), temporally separate the oxygen sensitive nitrogen fixation from oxygen evolving photosynthesis not only under diurnal cycles (LD) but also in continuous light (LL). However, recent reports demonstrate that the oscillations in LL occur with a shorter cycle time of ~11 h. We find that indeed, majority of the genes oscillate in LL with this cycle time. Genes that are upregulated at a particular time of day under diurnal cycle also get upregulated at an equivalent metabolic phase under LL suggesting tight coupling of various cellular events with each other and with the cell's metabolic status. A number of metabolic processes get upregulated in a coordinated fashion during the respiratory phase under LL including glycogen degradation, glycolysis, oxidative pentose phosphate pathway, and tricarboxylic acid cycle. These precede nitrogen fixation apparently to ensure sufficient energy and anoxic environment needed for the nitrogenase enzyme. Photosynthetic phase sees upregulation of photosystem II, carbonate transport, carbon concentrating mechanism, RuBisCO, glycogen synthesis and light harvesting antenna pigment biosynthesis. In Synechococcus elongates PCC 7942, a non-nitrogen fixing cyanobacteria, expression of a relatively smaller fraction of genes oscillates under LL condition with the major periodicity being 24 h. In contrast, the entire cellular machinery of Cyanothece orchestrates coordinated oscillation in anticipation of the ensuing metabolic phase in both LD and LL. These results may have important implications in understanding the timing of various cellular events and in engineering cyanobacteria for biofuel production.

Characterization of three putative xylulose 5-phosphate/fructose 6-phosphate phosphoketolases in the cyanobacterium Anabaena sp. PCC 7120

Xylulose 5-phosphate/fructose 6-phosphate phosphoketolase (Xfp) is a key enzyme in the central carbohydrate metabolism in heterofermentative bacteria, in which enzymatic property of Xfps is well characterized. This is not the case in other microbes. The cyanobacterium Anabaena sp. PCC 7120 possesses three putative genes encoding Xfp, all1483, all2567, and alr1850. We purified three putative Xfps as recombinant proteins. The results of gel filtration indicated that these proteins form homomultimer complex. All1483 and All2567 showed phosphoketolase activity, whereas Alr1850 did not show the activity. Kinetic analyses demonstrated that substrates, fructose 6-phosphate and inorganic phosphate, are cooperatively bound to enzymes positively and negatively, respectively.

Extremophilic micro-algae and their potential contribution in biotechnology

Micro-algae have potential as sustainable sources of energy and products and alternative mode of agriculture. However, their mass cultivation is challenging due to low survival under harsh outdoor conditions and competition from other, undesired, species. Extremophilic micro-algae have a role to play by virtue of their ability to grow under acidic or alkaline pH, high temperature, light, CO2 level and metal concentration. In this review, we provide several examples of potential biotechnological applications of extremophilic micro-algae and the ranges of tolerated extremes. We also discuss the adaptive mechanisms of tolerance to these extremes. Analysis of phylogenetic relationship of the reported extremophiles suggests certain groups of the Kingdom Protista to be more tolerant to extremophilic conditions than other taxa. While extremophilic microalgae are beginning to be explored, much needs to be done in terms of the physiology, molecular biology, metabolic engineering and outdoor cultivation trials before their true potential is realized.

High cell density lipid rich cultivation of a novel microalgal isolate Chlorella sorokiniana FC6 IITG in a single-stage fed-batch mode under mixotrophic condition

A single-stage mixotrophic cultivation strategy was developed with a novel microalgal isolate Chlorella sorokiniana FC6 IITG for high cell density lipid-rich biomass generation. The strain was evaluated for growth and lipid content under different physico-chemical parameters, nutritional conditions and trophic modes. Finally, a single-stage mixotrophic fed-batch cultivation strategy was demonstrated with intermittent feeding of key nutrients along with dynamic increase in light intensity for high cell density biomass and sodium acetate as elicitor for lipid enrichment. The key findings: (i) glucose and sodium acetate was identified as growth supporting and lipid inducing nutrients, respectively; (ii) mixotrophic batch cultivation resulted in maximum biomass and lipid productivity (mgL(-1)day(-1)) of 455.5 and 111.85, respectively; (iii) single-stage mixotrophic fed-batch cultivation showed maximum biomass productivity of 1.93gL(-1)day(-1) (biomass titer 15.81gL(-1)) and lipid productivity of 550mgL(-1)day(-1); (iv) biodiesel properties were in accordance with international standards.

Integrated metabolic flux and omics analysis of Synechocystis sp. PCC 6803 under mixotrophic and photoheterotrophic conditions

Cyanobacteria have flexible metabolic capability that enables them to adapt to various environments. To investigate their underlying metabolic regulation mechanisms, we performed an integrated analysis of metabolic flux using transcriptomic and metabolomic data of a cyanobacterium Synechocystis sp. PCC 6803, under mixotrophic and photoheterotrophic conditions. The integrated analysis indicated drastic metabolic flux changes, with much smaller changes in gene expression levels and metabolite concentrations between the conditions, suggesting that the flux change was not caused mainly by the expression levels of the corresponding genes. Under photoheterotrophic conditions, created by the addition of the photosynthesis inhibitor atrazine in mixotrophic conditions, the result of metabolic flux analysis indicated the significant repression of carbon fixation and the activation of the oxidative pentose phosphate pathway (PPP). Moreover, we observed gluconeogenic activity of upstream of glycolysis, which enhanced the flux of the oxidative PPP to compensate for NADPH depletion due to the inhibition of the light reaction of photosynthesis. 'Omics' data suggested that these changes were probably caused by the repression of the gap1 gene, which functions as a control valve in the metabolic network. Since metabolic flux is the outcome of a complicated interplay of cellular components, integrating metabolic flux with other 'omics' layers can identify metabolic changes and narrow down these regulatory mechanisms more effectively.

Toward a systems-level understanding of gene regulatory, protein interaction, and metabolic networks in cyanobacteria

Cyanobacteria are essential primary producers in marine ecosystems, playing an important role in both carbon and nitrogen cycles. In the last decade, various genome sequencing and metagenomic projects have generated large amounts of genetic data for cyanobacteria. This wealth of data provides researchers with a new basis for the study of molecular adaptation, ecology and evolution of cyanobacteria, as well as for developing biotechnological applications. It also facilitates the use of multiplex techniques, i.e., expression profiling by high-throughput technologies such as microarrays, RNA-seq, and proteomics. However, exploration and analysis of these data is challenging, and often requires advanced computational methods. Also, they need to be integrated into our existing framework of knowledge to use them to draw reliable biological conclusions. Here, systems biology provides important tools. Especially, the construction and analysis of molecular networks has emerged as a powerful systems-level framework, with which to integrate such data, and to better understand biological relevant processes in these organisms. In this review, we provide an overview of the advances and experimental approaches undertaken using multiplex data from genomic, transcriptomic, proteomic, and metabolomic studies in cyanobacteria. Furthermore, we summarize currently available web-based tools dedicated to cyanobacteria, i.e., CyanoBase, CyanoEXpress, ProPortal, Cyanorak, CyanoBIKE, and CINPER. Finally, we present a case study for the freshwater model cyanobacteria, Synechocystis sp. PCC6803, to show the power of meta-analysis, and the potential to extrapolate acquired knowledge to the ecologically important marine cyanobacteria genus, Prochlorococcus.

Rhythmic and sustained oscillations in metabolism and gene expression of Cyanothece sp. ATCC 51142 under constant light

Cyanobacteria, a group of photosynthetic prokaryotes, oscillate between day and night time metabolisms with concomitant oscillations in gene expression in response to light/dark cycles (LD). The oscillations in gene expression have been shown to sustain in constant light (LL) with a free running period of 24 h in a model cyanobacterium Synechococcus elongatus PCC 7942. However, equivalent oscillations in metabolism are not reported under LL in this non-nitrogen fixing cyanobacterium. Here we focus on Cyanothece sp. ATCC 51142, a unicellular, nitrogen-fixing cyanobacterium known to temporally separate the processes of oxygenic photosynthesis and oxygen-sensitive nitrogen fixation. In a recent report, metabolism of Cyanothece 51142 has been shown to oscillate between photosynthetic and respiratory phases under LL with free running periods that are temperature dependent but significantly shorter than the circadian period. Further, the oscillations shift to circadian pattern at moderate cell densities that are concomitant with slower growth rates. Here we take this understanding forward and demonstrate that the ultradian rhythm under LL sustains at much higher cell densities when grown under turbulent regimes that simulate flashing light effect. Our results suggest that the ultradian rhythm in metabolism may be needed to support higher carbon and nitrogen requirements of rapidly growing cells under LL. With a comprehensive Real time PCR based gene expression analysis we account for key regulatory interactions and demonstrate the interplay between clock genes and the genes of key metabolic pathways. Further, we observe that several genes that peak at dusk in Synechococcus peak at dawn in Cyanothece and vice versa. The circadian rhythm of this organism appears to be more robust with peaking of genes in anticipation of the ensuing photosynthetic and respiratory metabolic phases.

Highlights from the Indo-US workshop "Cyanobacteria: molecular networks to biofuels" held at Lonavala, India during December 16-20, 2012

An Indo-US workshop on "Cyanobacteria: molecular networks to biofuels" was held December 16-20, 2012 at Lagoona Resort, Lonavala, India. The workshop was jointly organized by two of the authors, PPW, a chemical engineer and LAS, a biologist, thereby ensuring a broad and cross-disciplinary participation. The main objective of the workshop was to bring researchers from academia and industry of the two countries together with common interests in cyanobacteria or microalgae and derived biofuels. An exchange of ideas resulted from a series of oral and poster presentations and, importantly, through one-on-one discussions during tea breaks and meals. Another key objective was to introduce young researchers of India to the exciting field of cyanobacterial physiology, modeling, and biofuels. PhD students and early stage researchers were especially encouraged to participate and about half of the 75 participants belonged to this category. The rest were comprised of senior researchers, including 13-15 invited speakers from each country. Overall, twenty-four institutes from 12 states of India were represented. The deliberations, which are being compiled in the present special issue, revolved mainly around molecular aspects of cyanobacterial biofuels including metabolic engineering, networks, genetic regulation, circadian rhythms, and stress responses. Representatives of some key funding agencies and industry provided a perspective and opportunities in the field and for bilateral collaboration. This article summarizes deliberations that took place at the meeting and provides a bird's eye view of the ongoing research in the field in the two countries.