- Open Access
Photosynthetic conversion of CO2 to farnesyl diphosphate-derived phytochemicals (amorpha-4,11-diene and squalene) by engineered cyanobacteria
Biotechnology for Biofuels volume 9, Article number: 202 (2016)
Metabolic engineering of cyanobacteria has enabled photosynthetic conversion of CO2 to value-added chemicals as bio-solar cell factories. However, the production levels of isoprenoids in engineered cyanobacteria were quite low, compared to other microbial hosts. Therefore, modular optimization of multiple gene expressions for metabolic engineering of cyanobacteria is required for the production of farnesyl diphosphate-derived isoprenoids from CO2.
Here, we engineered Synechococcus elongatus PCC 7942 with modular metabolic pathways consisting of the methylerythritol phosphate pathway enzymes and the amorphadiene synthase for production of amorpha-4,11-diene, resulting in significantly increased levels (23-fold) of amorpha-4,11-diene (19.8 mg/L) in the best strain relative to a parental strain. Replacing amorphadiene synthase with squalene synthase led to the synthesis of a high amount of squalene (4.98 mg/L/OD730). Overexpression of farnesyl diphosphate synthase is the most critical factor for the significant production, whereas overexpression of 1-deoxy-d-xylulose 5-phosphate reductase is detrimental to the cell growth and the production. Additionally, the cyanobacterial growth inhibition was alleviated by expressing a terpene synthase in S. elongatus PCC 7942 strain with the optimized MEP pathway only (SeHL33).
This is the first demonstration of photosynthetic production of amorpha-4,11-diene from CO2 in cyanobacteria and production of squalene in S. elongatus PCC 7942. Our optimized modular OverMEP strain (SeHL33) with either co-expression of ADS or SQS demonstrated the highest production levels of amorpha-4,11-diene and squalene, which could expand the list of farnesyl diphosphate-derived isoprenoids from CO2 as bio-solar cell factories.
Terpenoids extracted from plants are widely used as flavors, fragrances, and medicines. Because these compounds are naturally produced in small quantities, plant metabolic engineering has been applied in different compartments for the overproduction of target terpenoids: mono- and diterpenes and carotenoids in the plastids and sesqui- and triterpenes in the cytosol [1, 2]. The methylerythritol phosphate (MEP) pathway has been engineered to enhance the production of volatilized terpenoids [i.e., geranyl diphosphate (GPP)-derived monoterpene] such as linalool , geraniol , and limonene  by overexpressing the corresponding terpene synthases with engineered promoters in the plastids of carnation, tomato, and tobacco, respectively. On the other hand, metabolic engineering studies for higher production of sesquiterpene (i.e., artemisinin) has been extensively performed in the cytosol of Artemisia annua L. by overexpressing the mevalonate pathway enzymes , amorpha-4,11-diene synthase (ADS) , transcriptional factors to increase the precursor supply and repressing enzymes that compete for the precursor using the RNA interference-based suppression . Also, synthetic biology has been employed to produce artemisinin by reconstructing the heterologous metabolic pathway and optimizing the gene expression of ADS and its cognitive enzymes in a non-native fast-growing crop plant (i.e., tobacco)  or to produce various farnesyl diphosphate (FPP)-derived chemicals by altering the catalytic functions of terpene cyclases in A. annua L. .
Heterologous terpenoids including mono- and sesquiterpenes and carotenoids have also been produced in engineered microbes without compartments. For microbial production of carotenoids via the MEP pathway, the strains have been engineered by increasing the accumulation of precursors  or optimizing phenotype through high-throughput screening methods [12, 13] with expressing heterologous carotenoid biosynthetic enzymes. Besides, the heterologous mevalonate pathway has been introduced to produce high amounts of heterologous monoterpenes such as limonene and pinene in Escherichia coli [14, 15] and sesquiterpenes such as amorpha-4,11-diene and farnesene in E. coli [16, 17]. The recently, engineered Saccharomyces cerevisiae has produced heterologous artemisinic acid [18, 19], and the E. coli-S. cerevisiae co-culture system has been engineered to produce oxygenated taxanes and functionalized sesquiterpenes (i.e., nootkatone) .
As the MEP pathway exists in the plastid of plants, cyanobacteria, CO2-fixing photoautotrophs, have the sole MEP pathway to form terpenoids and precursors of chlorophyll, of which engineering has taken advantage of higher rates of photosynthesis in cyanobacteria than land plants . Most recent metabolic engineering of cyanobacteria has been carried out to directly convert CO2 to various biochemicals  including isobutyrlaldehyde (1.1 g/L) , n-butanol (404 mg/L) , and 2,3-butanediol (2.3 g/L) . However, the production levels of isoprenoids in engineered cyanobacteria showed significantly low titers when heterologous terpene synthase was expressed only in cyanobacteria: isoprene (approximately, 0.3 mg/L) , the monoterpenes limonene (4 mg/L)  and β-phellandrene (0.07 mg/L) , the sesquiterpene bisabolene (0.6 mg/L)  and triterpene squalene (0.6 mg/L/OD730) , and carotenoids (1.1 mg/L/OD730) . Because the enzymatic activity of terpene synthase is quite low, an optimized cyanobacterial chassis for isoprenoid production has not been constructed. Also, a majority of cyanobacteria do not produce sesquiterpene and triterpene, although cyanobacterial sesquiterpene synthases have been recently identified in N2-fixing filamentous cyanobacteria (Nostoc punctiforme PCC 73102 and Nostoc sp. Strain PCC 7120) . Thus, a modular metabolic engineering of cyanobacteria is required to provide a synthetic cyanobacterial chassis and to produce high levels of GPP-derived or FPP-derived isoprenoids production, respectively.
In this work, we describe the metabolic engineering of Synechococcus elongatus PCC 7942 as a model system with modular expression of the genes encoding the MEP pathway enzymes and terpene phosphate synthases (TPSs). This engineering allowed for the exogenous production of amorpha-4,11-diene and endogenous production of squalene directly from CO2.
Results and discussion
Designing modularized pathways for the production of amorpha-4,11-diene and squalene
S. elongatus PCC 7942 possesses the MEP pathway that is used to synthesize isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) from CO2. In plants and eukaryotes, biosynthesis of amorpha-4,11-diene and squalene are performed from a substrate FPP or two molecules of FPP, respectively. Thus, we modularized the biosynthesis pathways into two modules (OverMEP and TPS) for optimization [16, 32] (Fig. 1a). First, to increase the intracellular concentration of the FPP substrate, we overexpressed the heterologous genes encoding key enzymes in the MEP pathway and an IspA enzyme of E. coli using chromosomal integration vectors (pSe1B1s ) targeting neutral site I (NSI) of S. elongatus PCC 7942. Four E. coli genes (dxs, dxr, idi, ispA) documented to increase lycopene [13, 34, 35] were selected and a combinatorial approach of metabolic engineering was applied to identify the bottleneck of the cyanobacterial MEP pathway toward FPP (Fig. 1b).
For the OverMEP modules, overexpressions of the single gene of either the dxs (encoding for 1-deoxy-d-xylulose-5-phosphate synthase) or idi (encoding for isopentenyl diphosphate isomerase) were tested to confirm if they are good engineering targets as shown previously, resulting in the strains SeHL11A and SeHL12A. In addition, the dxr gene (encoding for 1-deoxy-d-xylulose 5-phosphate reductase) and/or ispA (encoding for farnesyl diphosphate synthase) genes were chosen for the overexpression in a combination with the dxs or/and idi genes to investigate whether or not Dxr or IspA is the rate-limiting enzyme in the MEP pathway. As a result, nine recombinant S. elongatus strains were generated for the OverMEP modules (Table 1).
Next, a heterologous TPS gene was introduced in S. elongatus PCC 7942 at neutral site II (NSII) using the second chromosomal integration vector (pSe2B1 k) as a TPS module. For the production of amorpha-4,11-diene and squalene, the ADS gene  encoding for amorpha-4,11-diene synthase from A. annua and the SQS gene  encoding for squalene synthase from S. cerevisiae, respectively, were also inserted into the chromosome. Based on the strains with OverMEP modules, nine and three S. elongatus strains were constructed for amorpha-4,11-diene and squalene, respectively (Table 1).
Cyanobacterial growth for the production of amorpha-4,11-diene from CO2 with in situ hexadecane overlay
To reduce the cellular toxicity of mono- or sesquiterpenes in microbial cell culture, solvent overlays are applied for in situ extraction [16, 37]. The dodecane overlay has been used to extract limonene and bisabolene from Synechococcus sp. PCC 7002 . The relative activity of the cells using a biocompatible organic solvent was correlated with the log Poctanol, an indicator of the hydrophobicity of the solvent . Here, solvents with log Poctanol > 5.5 (dodecane, hexadecane, tetradecane) were selected as a biocompatible organic solvent for S. elongatus PCC 7942. Then, the cellular toxicity of each solvent in S. elongatus PCC 7942 was investigated by adding 20 % solvent into the culture medium. The growth of S. elongatus was slightly inhibited by the addition of dodecane (81.3 % growth of the wild type at 5 days), whereas hexadecane or tetradecane overlay did not inhibit the cell growth (103 or 93.6 %) (Fig. 2). Thus, we selected 20 % (v/v) hexadecane overlay as in situ extraction of for amorpha-4,11-diene.
With the hexadecane overlay, the growth of the recombinant S. elongatus strains (SeHL11A, SeHL12A, SeHL31A, SeHL33A) for amorpha-4,11-diene production showed the 18, 8, 17, 18 % growth reductions, respectively, compared to the wild type (Fig. 3). The strain SeHL22A (Dxs-Idi overexpressed) showed severe growth inhibition (less than 48 % growth of the wild type), although a strain with a single overexpression of either dxs or idi (SeHL11A and SeHL12A) grew better than SeHL22A. Interestingly, strain SeHL21A, SeHL23A SeHL32A, and SeHL41A, where overexpression of Dxr was included, showed severe growth inhibition as well (<54, 66, 55, 54 % growth of the wild type, respectively).
Photosynthetic production of amorpha-4,11-diene from CO2
Using in situ extraction and GC–MS analysis, nine recombinant S. elongatus strains with the OverMEP and ADS produced amorpha-4,11-diene that was identical to the amorpha-4,11-diene produced by recombinant E. coli strain 2p  and secreted into the hexadecane overlay (Fig. 3a, b), whereas the wild type did not produce amorpha-4,11-diene at all. Strain SeHL11A (Dxs overexpressed) and SeHL12A (Idi overexpressed) produced 0.87 mg/L and 0.84 mg/L of amorpha-4,11-diene, respectively. For a combinatorial approach, the co-expression of Dxs and Idi (strain SeHL22A) resulted in an increased level (1.34 mg/L) of amorpha-4,11-diene by 1.54- and 1.59-fold over those produced by strain SeHL11A and SeHL12A, respectively. The next combinatorial engineering was co-expression of Dxr in previously engineered strains. However, co-expression of Dxr with either Dxs or Idi did not increase the levels of amorpha-4,11-diene (SeHL21A and SeHL23A) (Fig. 3c). Strain SeHL31A (Dxs-Dxr-Idi overexpressed) showed even lower levels (0.435 mg/L) than strain SeHL22A by 2.95-fold. Besides Dxr, co-overexpression of IspA was applied to strain SeHL22A and SeHL23A to improve the production levels. As a result, strain SeHL32A (Dxr-Idi-IspA overexpressed) showed 12-fold increased levels of amorpha-4,11-diene (3.59 mg/L) relative to strain SeHL23A (Dxr-Idi overexpressed; 0.3 mg/L of amorpha-4,11-diene). Moreover, strain SeHL33A (Dxs-Idi-IspA overexpressed) showed the significantly increased levels of amorpha-4,11-diene (19.8 mg/L) than SeHL22A (Dxs-Idi overexpressed) by 14.8-fold. As a full construct, SeHL41A (Dxs-Dxr-Idi-IspA overexpressed) showed only 2.54 mg/L of amorpha-4,11-diene, 7.8-fold lower than SeHL33A. Finally, we engineered S. elongatus PCC 7942 to produce 19.8 mg/L amorpha-4,11-diene, which is the highest levels among cyanobacterial terpenoids production reported so far.
Based on the results of the production of amorpha-4,11-diene and combinatorial engineering of the MEP pathway, we identified three key factors to produce high levels of amorpha-4,11-diene in S. elongatus PCC 7942 using the MEP pathway. First, co-overexpression of Dxs and Idi from E. coli in S. elongatus PCC 7942 was a good engineering choice to enhance the production of amorpha-4,11-diene from IPP and DMAPP as shown in other studies of isoprenoids [30, 40, 41]. Limonene production  (0.27 mg/L for 14 days) was increased by 2.3-fold in Anavaena sp. PCC 7120 by co-expressing the limonene synthase (LS) with synthetic Dxs-Idi from Haematococcus pluvialis and geranyl diphosphate synthase (GPPS) from Mycoplasma tuberculosis, compared to the strain with limonene synthase only. Also, co-expression of the native Dxs-Idi with CrtE (or, GPPS) has led 1.4-fold higher production of limonene (0.05 mg/L/d) in Synechocystis sp. PCC 6803 than the parental strain with limonene synthase only.
However, co-expression of Dxr negatively impacted the levels of amorpha-4,11-diene in S. elongatus PCC 7942. Engineered strains expressing with Dxr (SeHL21A, SeHL23A, SeHL41A) showed lower levels of both amorpha-4,11-diene production and cell growth than the parental strains (SeHL11A, SeHL12A, SeHL33A), respectively. Strain SeHL31A showed lower levels of amorpha-4,11-diene than the parental strains (SeHL22A). This negative effect has been also observed in Synechococcus leopoliensis  in which the DMAPP content was not influenced by additional overexpression of Dxr. In addition, overexpression of Dxr has not shown much impact on metabolic balances of the MEP pathway in Arabidopsis [43, 44]. As little is known about the multiple levels of the MEP pathway regulation in S. elongatus PCC 7942, it is not clear why overexpression of Dxr negatively affected production. Fine-tuning gene expression, such as with the RBS calculator  to fine-tune production or by adding protein degradation tags  to Dxr to fine-tune degradation, may be required to balance the gene expression of Dxr for the MEP pathway flux, because E. coli Dxr has lower K m value (30 µM) and faster k cat (100 s−1) than the cyanobacterial enzyme (K m of 134 µM and k cat of 5 s−1) .
Overexpression of E. coli’s IspA was the most crucial factor to increase amorpha-4,11-diene production in S. elongatus PCC 7942, where the best producer, strain SeHL33A, showed 22.8-fold higher production than the first engineered strain (SeHL11A) without growth retardation. Whereas E. coli has been engineered with IspA for isoprenoid production [17, 35], overexpression of Idi and Dxs was regarded as a key factor to improve production [35, 48, 49]. There have been no FPP synthase activities experimentally verified in S. elongatus PCC 7942. Thus, heterologous expression of sesquiterpene synthases without co-expression of FPP synthases in cyanobacteria could produce relatively low amounts of sesquiterpenes [27, 50] compared to the amounts of monoterpenes produced when monoterpene synthases are expressed (i.e., limonene [40, 41]). Overall, multi-omics technologies that integrate metabolomics, proteomics and network-modeling will allow understanding the complex regulation and metabolic aspects including identification of the rate-limiting step in metabolic engineering of cyanobacteria for isoprenoid productions . Furthermore, directed evolution and protein fusions of IspA to GPPS or TPS [17, 52] of E. coli could provide a chance to improve production of sesquiterpenes with chain length specificity. Besides IspA of E. coli, overexpression of IspA from different sources could be also useful to engineer cyanobacteria.
Photosynthetic production of squalene from CO2
To determine whether the OverMEP module can be used as a basis for isoprenoid production, we engineered the strain SeHL33 (Dxs-Idi-IspA overexpressed) with co-expression with SQS as another TPS module, since SQS catalyzes the formation of squalene, a triterpene synthesized from two molecules of FPP. Three recombinant S. elongatus strains were constructed for the production of squalene (Fig. 4). The wild-type S. elongatus PCC 7942 strain did not produce squalene, and a very low amount of squalene (0.0001 mg/L/OD730) was produced in a cell by expressing SQS alone (SeSC00S). Similar to SeHL33A, the SeSC33S strain (Dxs-Idi-IspA overexpressed) produced a high amount of squalene in S. elongatus PCC 7942 (4.98 mg/L/OD730 ±0.90), which was 50,000-fold higher than the strain SeSC00S and the highest levels reported so far in cyanobacteria. In the strain SeSC41S (Dxs, Dxr, Idi, and IspA overexpressed), the squalene content (0.13 mg/L/OD730 ±0.01) was reduced by 38.3-fold compared to SeSC33S. The increased or decreased levels were consistent with the previous results of amorpha-4,11-diene production using the SeHL33 or SeHL41 module, respectively. Thus, squalene production studies also support that overexpression of either Dxr was a negative factor or IspA was the most crucial factor in cyanobacteria to increase FPP-derived isoprenoids such as sesquiterpenes and triterpenes along with overexpression of Dxs and Idi, respectively. With respect to the enhanced production of amorpha-4,11-diene and squalene, the SeHL33 with modularized OverMEP proved to be orthogonally functional and expandable to produce other sesquiterpenes or triterpenes with the corresponding TPS’s.
Alleviating the cyanobacterial growth inhibition by expressing a terpene synthase in S. elongatus PCC 7942 strain with the optimized MEP pathway only (SeHL33)
To expand the FPP-derived isoprenoids production using cyanobacterial strain expressing Dxs, Idi, and IspA (SeHL33), we compared the strain SeHL33 as a basis with amorphadiene-producing SeHL33A and squalene-producing SeSC33S. Previously, the toxicity of FPP, which can accumulate in engineered E. coli, has been alleviated by increasing the activity of the terpene cyclase . Thus, we speculated that the strain SeHL33 could also accumulate FPP. In absence of cyanobacterial analysis of proteomics and intracellular FPP accumulation patterns , we examined whether or not the strain SeHL33 showed a growth defect due to possible FPP toxicity. As results, the cultures of three recombinant S. elongatus strains (strains SeHL33, SeHL33A and SeSC33S) were monitored, showing that the growth of SeHL33 was significantly lower than that of SeHL33A and SeSC33S (Fig. 5a). Moreover, cyanobacterial growth can be measured indirectly using the changes of chlorophyll content . Also, the biosynthesis of Chl a is connected to the engineered MEP pathway via geranylgeranyl diphosphate. Chlorophyll a contents of the wild-type, SeHL33A, SeSC33S, or SeHL33 strains were determined at 2 and 8 days, resulting in that the Chl a content of the strain SeHL33 at 8 days was significantly lower than the Chl a contents of other strains (Fig. 5b). Thus, low Chl a content (mg/L/OD730) of the SeHL33 strain might be related to the growth defect of of SeHL33.
However, our engineered strain SeHL33 with either an ADS or SQS did not show the negative effect on the growth. Thus, alleviation of the cellular toxicity in the engineered MEP pathway with co-expression of terpene synthase in this study was consistent with the previous result in engineered E. coli . Furthermore, genome-wide transcriptomics approach could be employed using the strain SeHL33 to increase the levels of isoprenoid production by identifying the stress-response promoters for engineering dynamic pathway regulation . In-depth metabolite analysis (or metabolomics) for the strains could be also needed to understand the bottlenecks and to improve isoprenoid biosynthesis. Since our engineered strains with a modular OverMEP and either an ADS or SQS did not show a negative effect on the growth, we could expand engineering of SeHL33 with different TPSs converting FPP to various sesqui- or triterpenes.
Beyond the proof of concept of isoprenoid productions in cyanobacteria [27, 29], the effort of metabolic engineering of cyanobacteria expanded the list of products from CO2 as bio-solar cell factories in this study. Although isoprenoid titers may be influenced by the type of strains, the characteristics of TPS and culture conditions and the high levels of amounts of amorpha-4,11-diene and squalene were produced in this work when pathway optimization was a main strategy of metabolic engineering for isoprenoid production (Fig. 6; Additional file 1: Table S1). Combined with the pathway optimization, increasing the low activity (slow k cat) of terpene synthase has been shown to be very crucial for the enhanced production of isoprenoid in cyanobacteria [55, 56]. Recently, CpcB as fusion partner or a key MEP enzyme (Idi) was fused to a terpene synthase that resulted in significantly high levels of production of β-phellandrene [55, 57] and isoprene , respectively. The engineered strains for isoprene production showed the 40 % of photosynthetically fixed carbon fluxes toward isoprenoid biosynthesis pathways. However, there is still much room left for engineering cyanobacteria to improve the titers, compared to other microbial hosts (293 mg/L/OD600 of amorphadiene in E. coli  and 242 mg/L/OD600 of amorphadiene in S. cerevisiae ; 150 mg/L and 55 mg/gDW of squalene in E. coli ). For example, installing efflux pump protein could improve the amorphadiene in cyanobacteria, of which screening strategy has been successful in E. coli .
To our knowledge, our optimized modular OverMEP strain (SeHL33) with either co-expression of ADS or SQS demonstrated the highest production levels of amorpha-4,11-diene and squalene from CO2.
Chemicals and reagents
All chemicals were purchased from Sigma-Aldrich (St. Louis, MO, USA) unless otherwise specified. Restriction enzymes, Phusion DNA polymerase, and ligases were purchased from Fermentas (Thermo Fisher Scientific Inc., MA, USA).
E. coli strain DH5α  was used for gene cloning and grown in Luria–Bertani medium (containing per liter: 10 g tryptone, 5 g yeast extract, and 10 g NaCl) at 37 °C, and when appropriate the medium was supplemented with 50 μg/mL kanamycin or 100 μg/mL spectinomycin. All plasmids were derived from SyneBrick expression plasmids (standard vectors for chromosomal integration at the neutral site I (NSI) or II (NSII) using a pSyn_1; pSe1Bb1 s-gfp  and pSe2Bb1 k-gfp ) using the BglBrick standard cloning method  (Table 1). To overexpress key enzymes of the MEP pathway, the dxs gene encoding for 1-deoxy-d-xylulose-5-phosphate synthase (Dxs), the dxr gene encoding for 1-deoxy-d-xylulose 5- phosphate reductase (Dxr), the idi gene encoding for isopentenyl diphosphate isomerase (Idi), and the ispA gene encoding for farnesyl diphosphate synthase (IspA) from E. coli were codon-optimized to S. elongatus PCC 7942 using Gene Designer 2.0 software (DNA2.0, Menlo Park, CA, USA) with a codon usage matrix of S. elongatus PCC 7942 and synthesized and cloned into a pSe1Bb1 s vector (targeting at NSI). For production of amorpha-4,11-diene and squalene, the ADS gene  encoding for amorpha-4,11-diene synthase from A. annua and the truncated SQS gene  encoding for squalene synthase from S. cerevisiae were also synthesized and cloned into a pSe2Bb1 k vector (targeting at NSII), respectively.
Cyanobacterial strain construction and transformation
Transformation of S. elongatus PCC 7942 was performed as described previously . The cyanobacterial strains were transformed by incubating cells at a mid-log phase (OD730 of 1–2) with 100 ng of plasmid DNA for 24 h in the dark. The mixed culture was then spread on BG-11 plates supplemented with appropriate antibiotics for selection of successful recombination. For selection and culture maintenance, 10 μg/mL spectinomycin or/and 10 μg/mL kanamycin were added into BG-11 agar plates and the BG-11 medium where appropriate. Sub-culture of a single colony was performed to prevent chromosomal segregation. The strains were confirmed by PCR to verify chromosomal integration of targets into either the NSI or NSII (Fig. 1) and the DNA sequences were also correctly verified using a pair of oligonucleotides (Additional file 1: Table S1). S. elongatus PCC 7942 strains were constructed by transforming the overexpressed MEP strains with plasmids pSe2Bb1 k-ADS and pSe2Bb1 k-SQS, respectively. Genotypes of recombinant S. elongatus strains are listed in Table 1.
Growth condition for amorphadiene and squalene from engineered cyanobacteria
S. elongatus PCC 7942 derivatives for production of amorpha-4,11-diene and squalene were cultivated at 30 °C in the 100 mL culture under continuous fluorescent light (100 µmol photons/m2/s) in BG-11 medium supplanted with 10 mM MOPS (pH 8.0). 5 % (v/v) CO2 gas was supplied at a constant flow rate of 10 mL/min into the medium. 10 μg/mL spectinomycin or/and 10 μg/mL kanamycin were supplemented for selection pressure. 0.5 mM isopropyl-β-d-1-thiogalactopyranoside (IPTG) was supplemented into the culture medium at 24 h after inoculation for inductions. Especially, 20 % (v/v) hexadecane was added for in situ extraction of amorphadiene.
Quantification of amorpha-4,11-diene and squalene
Quantification of amorpha-4,11-diene (mg of amorpha-4,11-diene per liter of culture volume) was performed as described previously [16, 64]. The hexadecane layer in the culture was collected and diluted 1:10 with ethyl acetate containing 5 μg/mL of (-)-trans-caryophyllene as an internal standard. Hexadecane/ethyl acetates extracts were analyzed using a GC–MS (Agilent 6890 gas chromatograph interfaced with an Agilent 5975 MSD). The extracted sample (1 μL) was injected in split mode (10:1) at an injector temperature of 220 °C and separated using an Ultra-2 capillary column (33 m × 0.2 mm i.d., 0.11 μm film thickness; Agilent Technologies). The initial oven temperature was 100 °C for 3 min and ramped up to 200 °C at 7 °C/min for a total run time of 17.29 min. Helium (99.9999 %) was used as the carrier gas (1.0 mL/min. constant flow at an oven temperature of 150 °C). The ion source temperature was 220 °C. The mass spectrometer was operated at 70 eV in the electron ionization mode with selected ion monitoring (SIM) and the mass range was 35–400 m/z. Biosynthetic amorphadiene was obtained from the culture of E. coli DH1 harboring the pBbA5c-MevT-MBIS plasmid and pADS (Strain 2p) .
For quantification of squalene, the extraction method was performed as described previously  with slight modification. 50 mL of the culture was used for extraction. After centrifugation at 3000×g for 10 min, the cell pellets were re-suspended with a 2 mL mixture of chloroform and methanol (1:2 ratio) and liquid–liquid extraction occurred for 30 min at room temperature. After additional centrifugation (16,000×g for 3 min), the supernatants (200 μL) were collected and the samples supplemented with 20 μg/mL of 1-phenyloctadecane as internal standard were analyzed using a GC–MS. The extracted sample was analyzed using a gas chromatography equipped with a gas chromatography–mass spectrometry (GC-MSD; Agilent Technologies, Santa Clara, CA) equipped with a U2 capillary column (33 m × 0.25 mm, film thickness 0.25 mm, Agilent Technologies, Santa Clara, CA), carrier gas: He (1 mL/min), oven temperature: 90–290 °C (increase rate 6 °C/min). The mass spectrometer was operated at 70 eV in the electron ionization mode with selected ion monitoring (SIM; 69 m/z) and mass range was 35–400 m/z.
Measurement of chlorophyll a
Cyanobacterial cell cultures (1.5 mL) were harvested and the pigments of cell pellets were extracted in 100 % acetone at 50 °C for 15 min. After centrifugation (16,000×g for 3 min), the pigment extract was analyzed for chlorophyll a (Chl a) by following the previous spectrophotometric assay . The supernatant was obtained and the absorbance was measured at 400–700 nm on Cary 60 UV–Vis spectrophotometer (Agilent technologies, CA, USA).
gas chromatography–mass spectrometry
geranyl diphosphate (GPP)
terpene phosphate synthase
Kirby J, Keasling JD. Biosynthesis of plant isoprenoids: perspectives for microbial engineering. Annu Rev Plant Biol. 2009;60:335–55.
Aharoni A, Jongsma MA, Bouwmeester HJ. Volatile science? Metabolic engineering of terpenoids in plants. Trends Plant Sci. 2005;10(12):594–602.
Lavy M, Zuker A, Lewinsohn E, Larkov O, Ravid U, Vainstein A, et al. Linalool and linalool oxide production in transgenic carnation flowers expressing the Clarkia breweri linalool synthase gene. Mol Breed. 2002;9(2):103–11.
Davidovich-Rikanati R, Lewinsohn E, Bar E, Iijima Y, Pichersky E, Sitrit Y. Overexpression of the lemon basil alpha-zingiberene synthase gene increases both mono- and sesquiterpene contents in tomato fruit. Plant J. 2008;56(2):228–38.
Ohara K, Ujihara T, Endo T, Sato F, Yazaki K. Limonene production in tobacco with Perilla limonene synthase cDNA. J Exp Bot. 2003;54(393):2635–42.
Ram M, Khan MA, Jha P, Khan S, Kiran U, Ahmad MM, et al. HMG-CoA reductase limits artemisinin biosynthesis and accumulation in Artemisia annua L. plants. Acta Physiol Plant. 2010;32(5):859–66.
Ma CF, Wang HH, Lu X, Wang H, Xu GW, Liu BY. Terpenoid metabolic profiling analysis of transgenic Artemisia annua L. by comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry. Metabolomics. 2009;5(4):497–506.
Zhang L, Jing FY, Li FP, Li MY, Wang YL, Wang GF, et al. Development of transgenic Artemisia annua (Chinese wormwood) plants with an enhanced content of artemisinin, an effective anti-malarial drug, by hairpin-RNA-mediated gene silencing. Biotechnol Appl Biochem. 2009;52:199–207.
Farhi M, Marhevka E, Ben-Ari J, Algamas-Dimantov A, Liang Z, Zeevi V, et al. Generation of the potent anti-malarial drug artemisinin in tobacco. Nat Biotechnol. 2011;29(12):1072–4.
Salmon M, Laurendon C, Vardakou M, Cheema J, Defernez M, Green S, et al. Emergence of terpene cyclization in Artemisia annua. Nat Commun. 2015;6:6143.
Kim SW, Keasling JD. Metabolic engineering of the nonmevalonate isopentenyl diphosphate synthesis pathway in Escherichia coli enhances lycopene production. Biotechnol Bioeng. 2001;72(4):408–15.
Alper H, Miyaoku K, Stephanopoulos G. Construction of lycopene-overproducing E. coli strains by combining systematic and combinatorial gene knockout targets. Nat Biotechnol. 2005;23(5):612–6.
Choi HS, Lee SY, Kim TY, Woo HM. In silico identification of gene amplification targets for improvement of lycopene production. Appl Environ Microb. 2010;76(10):3097–105.
Alonso-Gutierrez J, Chan R, Batth TS, Adams PD, Keasling JD, Petzold CJ, et al. Metabolic engineering of Escherichia coli for limonene and perillyl alcohol production. Metab Eng. 2013;19:33–41.
Sarria S, Wong B, Martin HG, Keasling JD, Peralta-Yahya P. Microbial synthesis of pinene. Acs Synth Biol. 2014;3(7):466–75.
Martin VJ, Pitera DJ, Withers ST, Newman JD, Keasling JD. Engineering a mevalonate pathway in Escherichia coli for production of terpenoids. Nat Biotechnol. 2003;21(7):796–802.
Wang C, Zhou J, Jang HJ, Yoon SH, Kim JY, Lee SG, et al. Engineered heterologous FPP synthases-mediated Z, E-FPP synthesis in E. coli. Metab Eng. 2013;18:53–9.
Westfall PJ, Pitera DJ, Lenihan JR, Eng D, Woolard FX, Regentin R, et al. Production of amorphadiene in yeast, and its conversion to dihydroartemisinic acid, precursor to the antimalarial agent artemisinin. Proc Natl Acad Sci USA. 2012;109(3):111–8.
Ro D-K, Paradise EM, Ouellet M, Fisher KJ, Newman KL, Ndungu JM, et al. Production of the antimalarial drug precursor artemisinic acid in engineered yeast. Nature. 2006;440(7086):940–3.
Zhou K, Qiao KJ, Edgar S, Stephanopoulos G. Distributing a metabolic pathway among a microbial consortium enhances production of natural products. Nat Biotechnol. 2015;33(4):377–83.
Melis A. Solar energy conversion efficiencies in photosynthesis: minimizing the chlorophyll antennae to maximize efficiency. Plant Sci. 2009;177(4):272–80.
Savakis P, Hellingwerf KJ. Engineering cyanobacteria for direct biofuel production from CO2. Curr Opin Biotechnol. 2015;33:8–14.
Atsumi S, Higashide W, Liao JC. Direct photosynthetic recycling of carbon dioxide to isobutyraldehyde. Nat Biotechnol. 2009;27(12):1177–80.
Lan EI, Ro SY, Liao JC. Oxygen-tolerant coenzyme A-acylating aldehyde dehydrogenase facilitates efficient photosynthetic n-butanol biosynthesis in cyanobacteria. Energy Environ Sci. 2013;6(9):2672–81.
Oliver JW, Machado IM, Yoneda H, Atsumi S. Cyanobacterial conversion of carbon dioxide to 2,3-butanediol. Proc Natl Acad Sci USA. 2013;110(4):1249–54.
Bentley FK, Zurbriggen A, Melis A. Heterologous expression of the mevalonic acid pathway in cyanobacteria enhances endogenous carbon partitioning to isoprene. Mol Plant. 2014;7(1):71–86.
Davies FK, Work VH, Beliaev AS, Posewitz MC. Engineering limonene and bisabolene production in wild type and a glycogen-deficient mutant of synechococcus sp. pcc 7002. Front Bioeng Biotechnol. 2014;2:21.
Formighieri C, Melis A. Regulation of β-phellandrene synthase gene expression, recombinant protein accumulation, and monoterpene hydrocarbons production in Synechocystis transformants. Planta. 2014;240(2):309–24.
Englund E, Pattanaik B, Ubhayasekera SJ, Stensjo K, Bergquist J, Lindberg P. Production of squalene in Synechocystis sp. PCC 6803. PLoS ONE. 2014;9(3):e90270.
Kudoh K, Kawano Y, Hotta S, Sekine M, Watanabe T, Ihara M. Prerequisite for highly efficient isoprenoid production by cyanobacteria discovered through the over-expression of 1-deoxy-d-xylulose 5-phosphate synthase and carbon allocation analysis. J Biosci Bioeng. 2014;118(1):20–8.
Agger SA, Lopez-Gallego F, Hoye TR, Schmidt-Dannert C. Identification of sesquiterpene synthases from Nostoc punctiforme PCC 73102 and Nostoc sp. strain PCC 7120. J Bacteriol. 2008;190(18):6084–96.
Xu P, Gu Q, Wang W, Wong L, Bower AG, Collins CH, et al. Modular optimization of multi-gene pathways for fatty acids production in E. coli. Nat Commun. 2013;4:1409.
Chwa JW, Kim WJ, Sim SJ, Um Y, Woo HM. Engineering of a modular and synthetic phosphoketolase pathway for photosynthetic production of acetone from CO2 in Synechococcus elongatus PCC 7942 under light and aerobic condition. Plant Biotechnol J. 2016;14:1768–76.
Jin YS, Stephanopoulos G. Multi-dimensional gene target search for improving lycopene biosynthesis in Escherichia coli. Metab Eng. 2007;9(4):337–47.
Wang HH, Isaacs FJ, Carr PA, Sun ZZ, Xu G, Forest CR, et al. Programming cells by multiplex genome engineering and accelerated evolution. Nature. 2009;460(7257):894–8.
Zhang D, Jennings SM, Robinson GW, Poulter CD. Yeast squalene synthase: expression, purification, and characterization of soluble recombinant enzyme. Arch Biochem Biophys. 1993;304(1):133–43.
Peralta-Yahya PP, Ouellet M, Chan R, Mukhopadhyay A, Keasling JD, Lee TS. Identification and microbial production of a terpene-based advanced biofuel. Nat Commun. 2011;2:483.
Leon R, Martin M, Vigara J, Vilchez C, Vega JM. Microalgae mediated photoproduction of beta-carotene in aqueous-organic two phase systems. Biomol Eng. 2003;20(4–6):177–82.
Redding-Johanson AM, Batth TS, Chan R, Krupa R, Szmidt HL, Adams PD, et al. Targeted proteomics for metabolic pathway optimization: application to terpene production. Metab Eng. 2011;13(2):194–203.
Halfmann C, Gu L, Zhou R. Engineering cyanobacteria for the production of a cyclic hydrocarbon fuel from CO2 and H2O. Green Chem. 2014;16(6):3175–85.
Kiyota H, Okudac Y, Ito M, Hirai MY, Ikeuchi M. Engineering of cyanobacteria for the photosynthetic production of limonene from CO2. J Biotechnol. 2014;185:1–7.
Miller B, Heuser T, Zimmer W. Functional involvement of a deoxy-d-xylulose 5-phosphate reductoisomerase gene harboring locus of Synechococcus leopoliensis in isoprenoid biosynthesis. FEBS Lett. 2000;481(3):221–6.
Xing S, Miao J, Li S, Qin G, Tang S, Li H, et al. Disruption of the 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR) gene results in albino, dwarf and defects in trichome initiation and stomata closure in Arabidopsis. Cell Res. 2010;20(6):688–700.
Mahmoud SS, Croteau RB. Metabolic engineering of essential oil yield and composition in mint by altering expression of deoxyxylulose phosphate reductoisomerase and menthofuran synthase. Proc Natl Acad Sci USA. 2001;98(15):8915–20.
Borujeni AE, Channarasappa AS, Salis HM. Translation rate is controlled by coupled trade-offs between site accessibility, selective RNA unfolding and sliding at upstream standby sites. Nucleic Acids Res. 2014;42(4):2646–59.
Torella JP, Ford TJ, Kim SN, Chen AM, Way JC, Silver PA. Tailored fatty acid synthesis via dynamic control of fatty acid elongation. Proc Natl Acad Sci USA. 2013;110(28):11290–5.
Murkin AS, Manning KA, Kholodar SA. Mechanism and inhibition of 1-deoxy-d-xylulose-5-phosphate reductoisomerase. Bioorg Chem. 2014;57:171–85.
Wang CW, Oh MK, Liao JC. Engineered isoprenoid pathway enhances astaxanthin production in Escherichia coli. Biotechnol Bioeng. 1999;62(2):235–41.
Yuan LZ, Rouviere PE, Larossa RA, Suh W. Chromosomal promoter replacement of the isoprenoid pathway for enhancing carotenoid production in E. coli. Metab Eng. 2006;8(1):79–90.
Reinsvold RE, Jinkerson RE, Radakovits R, Posewitz MC, Basu C. The production of the sesquiterpene β-caryophyllene in a transgenic strain of the cyanobacterium Synechocystis. J Plant Physiol. 2011;168(8):848–52.
Brunk E, George Kevin W, Alonso-Gutierrez J, Thompson M, Baidoo E, Wang G, et al. Characterizing strain variation in engineered E. coli using a multi-omics-based workflow. Cell Syst. 2016;2(5):335–46.
Lee PC, Petri R, Mijts BN, Watts KT, Schmidt-Dannert C. Directed evolution of Escherichia coli farnesyl diphosphate synthase (IspA) reveals novel structural determinants of chain length specificity. Metab Eng. 2005;7(1):18–26.
Dahl RH, Zhang F, Alonso-Gutierrez J, Baidoo E, Batth TS, Redding-Johanson AM, et al. Engineering dynamic pathway regulation using stress-response promoters. Nat Biotechnol. 2013;31(11):1039–46.
Li Y, Lin Y, Loughlin PC, Chen M. Optimization and effects of different culture conditions on growth of Halomicronema hongdechloris—a filamentous cyanobacterium containing chlorophyll f. Front Plant Sci. 2014;5:67.
Formighieri C, Melis A. Sustainable heterologous production of terpene hydrocarbons in cyanobacteria. Photosynth Res. 2016. doi:10.1007/s11120-016-0233-2.
Gao X, Gao F, Liu D, Zhang H, Nie X, Yang C. Engineering the methylerythritol phosphate pathway in cyanobacteria for photosynthetic isoprene production from CO2. Energy Environ Sci. 2016;9(4):1400–11.
Formighieri C, Melis A. A phycocyanin phellandrene synthase fusion enhances recombinant protein expression and beta-phellandrene (monoterpene) hydrocarbons production in Synechocystis (cyanobacteria). Metab Eng. 2015;32:116–24.
Anthony JR, Anthony LC, Nowroozi F, Kwon G, Newman JD, Keasling JD. Optimization of the mevalonate-based isoprenoid biosynthetic pathway in Escherichia coli for production of the anti-malarial drug precursor amorpha-4,11-diene. Metab Eng. 2009;11(1):13–9.
Katabami A, Li L, Iwasaki M, Furubayashi M, Saito K, Umeno D. Production of squalene by squalene synthases and their truncated mutants in Escherichia coli. J Biosci Bioeng. 2015;119(2):165–71.
Zhang C, Chen X, Stephanopoulos G, Too HP. Efflux transporter engineering markedly improves amorphadiene production in Escherichia coli. Biotechnol Bioeng. 2016;113:1755–63.
Hanahan D. Studies on transformation of Escherichia coli with plasmids. J Mol Biol. 1983;166(4):557–80.
Lee TS, Krupa RA, Zhang F, Hajimorad M, Holtz WJ, Prasad N, et al. BglBrick vectors and datasheets: a synthetic biology platform for gene expression. J Biol Eng. 2011;5:12.
Golden SS, Brusslan J, Haselkorn R. Genetic engineering of the cyanobacterial chromosome. Methods Enzymol. 1987;153:215–31.
Woo HM, Murray GW, Batth TS, Prasad N, Adams PD, Keasling JD, et al. Application of targeted proteomics and biological parts assembly in E. coli to optimize the biosynthesis of an anti-malarial drug precursor, amorpha-4,11-diene. Chem Eng Sci. 2013;103:21–8.
Jeffrey SW, Humphrey GF. New spectrophotometric equations for determining chlorophyll a, b, c1 and c2 in higher plants, algae and natural phytoplankton. Biochem Physiol Pfl. 1975;167(2):191–4.
SYC carried out the molecular genetic studies, performed the biological experiments, and drafted the manuscript. HJL carried out the molecular genetic studies, performed the biological experiments, and drafted the manuscript. JC conducted the GC–MS analyses, performed the biological experiments, and drafted the manuscript. JK conducted the GC–MS analyses, performed the biological experiments, and drafted the manuscript. SJS participated in the design of the study, helped to analyze the data and helped to draft the manuscript. YU participated in the design of the study and helped to analyze the data to draft the manuscript. YK participated in the design of the study and helped to analyze the data and to draft the manuscript. TSL participated in the design of the study and helped to analyze the data and to draft the manuscript. JDK participated in the design of the study and helped to analyze the data and to draft the manuscript. HMW conceived the study, participated in its design and coordination, and helped to analyze the data and to draft the manuscript. All authors read and approved the final manuscript.
The authors thank Dr. Min-Kyung Kang for the technical support.
JDK has financial conflicts of interest in Amyris Biotechnologies, LS9 Inc. and Lygos Inc. The authors declare that they have no competing interests.
Availability of data and materials
The online version of this article contains an Additional file.
Consent for publication
The authors have consented for publication.
This work was supported by Korea CCS R&D Center (KCRC) (Grant No. 2014M1A8A1049277) and SJS was supported by Grant No. NRF-2013R1A2A1A01015644, funded by the Korean Government (Ministry of Science, Information and Communications Technology (ICT) & Future Planning) (2016, University-Institute Cooperation program by NRF).
Sun Young Choi and Hyun Jeong Lee contributed equally to this work