Skip to main content

Engineering cyanobacteria as a new platform for producing taxol precursors directly from carbon dioxide

Abstract

Taxol serves as an efficient natural anticancer agent with extensive applications in the treatment of diverse malignancies. Although advances in synthetic biology have enabled the de novo synthesis of taxol precursors in various microbial chassis, the total biosynthesis of taxol remains challengable owing to the restricted oxidation efficiency in heterotrophic microbes. Here, we engineered Synechocystis sp. PCC 6803 with modular metabolic pathways consisting of the methylerythritol phosphate pathway enzymes and taxol biosynthetic enzymes for production of taxadiene-5α-ol (T5α-ol), the key oxygenated intermediate of taxol. The best strain DIGT-P560 produced up to 17.43 mg/L of oxygenated taxanes and 4.32 mg/L of T5α-ol. Moreover, transcriptomic analysis of DIGT-P560 revealed that establishing a oxygenated taxane flux may enhance photosynthetic electron transfer efficiency and central metabolism in the engineered strain to ameliorate the metabolic disturbances triggered by the incorporation of exogenous genes. This is the first demonstration of photosynthetic production of taxadiene-5α-ol from CO2 in cyanobacteria, highlighting the broad prospects of engineered cyanobacteria as bio-solar cell factories for valuable terpenoids production and expanding the ideas for further rational engineering and optimization.

Introduction

Paclitaxel (Taxol) is a highly potent natural antineoplastic agent that exhibits remarkable efficacy against breast, ovarian, and pulmonary malignancies. It impedes the proliferation of neoplastic cells by arresting the cell cycle at the G2/M phase (the cell phase before mitosis) by stabilizing microtubule formation [1]. Presently, taxol production predominantly involves three methodologies: botanical extraction, chemical/semi-chemical synthesis, and biosynthesis. Taxol was initially derived from the bark of Taxus brevifolia [2]. Subsequently, researchers developed structural analogs with more extensive taxane frameworks from needles and employed chemically modified approaches for large-scale production [3]. Nonetheless, the current processes of plant extraction and semi-chemical synthesis depend heavily on natural resources and the cultivation industry. The utilization of ecologically sustainable and green biosynthetic strategies has the potential to comprehensively resolve the supply–demand incongruity of such bioactive agents. Furthermore, it has substantial strategic significance for research on natural-product production.

Numerous structural and regulatory genes involved in taxol biosynthesis have been identified [4]. The initial step in taxol biosynthesis is catalyzed by taxadiene synthase (TASY), which involves cyclization of the common diterpenoid precursor geranylgeranyl pyrophosphate (GGPP) to yield taxadiene and iso-taxadiene. Subsequently, the tricyclic structure of taxadiene undergoes oxidation and acylation at various sites through a series of cytochrome P450 oxidases (referred to hereafter as P450s) and acetyltransferases. Taxadiene-5α-ol (referred to hereafter as T5α-ol), the primary oxygenated taxane, is catalyzed by taxadiene-5α-hydroxylase (CYP725A4, referred to as T5αH hereafter), which introduces a hydroxyl group at the C5 position of taxadiene and transfers the double bond from C4(5) to C4(20) [5]. Meanwhile, a cytochrome P450 reductase (CPR) is required to donate two NADPH-derived electrons for the catalytic reaction [6].

With the rapid advancement of synthetic biology, metabolic engineering strategies, protein engineering, and multigroup techniques have found extensive applications in taxol biosynthesis research. To date, specific taxol precursors have been de novo synthesized in various microorganisms, including Escherichia coli [7,8,9] and Sacharomyces cerevisiae [10,11,12]. The highest reported yield of taxadiene in E. coli has reached to 1 g/L through a multivariate modular metabolic engineering strategy, although only 58 mg/L of T5α-ol was produced [7]. By balancing P450 module expression, engineering CPR interactions and modifying the N-terminus, Biggs et al. have increased the yield of total oxygenated taxanes to 570 mg/L in E. coli [8]. Another ambitious endeavor was made to introduce eight taxol biosynthetic genes simultaneously into S. cerevisiae for the biosynthesis of Baccatin III; however, only trace amounts of T5α-ol were detected [10]. Coupled with advanced biological process strategies, the engineered S. cerevisiae yields 34 mg/L of T5α-ol and 11 mg/L of taxadien-5α-yl-acetate [12], which represents the highest titers in yeast cell factories. However, it is still challenging to assemble the complex taxol biosynthesis pathway in E. coli because of the absence of internal membranes, which are essential for P450s integration [13, 14], and the inefficient regeneration of nicotinamide adenine dinucleotide phosphate (NADPH) might jeopardize the high activity of P450s in microbial chassis such as E. coli and yeast [14].

As prokaryotic and photoautotrophic microorganisms, cyanobacteria are considered as non-food-based feedstock resources that can convert solar energy and carbon dioxide for the synthesis of various carbon-based compounds [15]. Since cyanobacteria naturally synthesize numerous photosynthetic pigments such as carotenoids and chlorophylls [16], their native GGPP pools are supposed to be ample for stable supply of diterpene precusor, making them well suited for terpenoids production. On the other hand, P450s generally need to receive two electrons from NADPH or NADH to promote the catalytic reactions [14], the presence of photosynthetic electron transport chain in cyanobacteria ensures the regeneration of reducing power, thereby directly support the high activity of P450s. Most importantly, the successful integration of eukaryotic P450s into the photosynthetic thylakoid membranes has recasted cyanobacteria as promising hosts for P450s expression [17,18,19]. Given all that, cyanobacteria are competent candidates for the direct conversion of CO2 into valuable secondary metabolites such as terpenoids [20,21,22].

In this study, as a proof of concept, multiple endeavors were undertaken to achieve the efficient biosynthesis of various taxol precursors in the cyanobacterial model strain Synechocystis sp. PCC 6803, hereafter referred to as Syn6803. The taxadiene biosynthetic pathway was initially established in Syn6803 by introducing heterologous GGPPS and TASY. By optimizing the upstream modules and increasing the biomass, the yield of taxadiene augmented approximately 11-fold. Subsequently, codon-optimized T5αH and its reductase partner CPR were integrated into the taxadiene-producing chassis with modular modifications implemented to enhance the efficiency of downstream oxygenated taxanes synthesis, including T5α-ol. Complemented by transcriptomic analysis, our aim was to elucidate the effects of introducing a heterologous taxol synthesis pathway on the intracellular energy and material metabolism of cyanobacterial chassis. Together, we sought to provide insights for establishing a more efficient photosynthetic platform for the direct conversion of CO2 into specialty compounds such as taxol and other terpenoids.

Materials and methods

Strain and cultivation conditions

All strains developed in this study were derived from Synechocystis sp. PCC 6803 wild type (hereafter described as WT). For flask shaking cultivation, the strains were grown at 30 °C in BG-11 medium [23] under continuous illumination of 25 µmol m−2 s−1, with shaking at 130 rpm for normal growth in a photoincubator shaker (ZQZY-AGS8, Zhichu, China). Genetically engineered strains were cultivated in BG11 medium containing the appropriate antibiotics (20 μg mL−1 spectinomycin, 20 μg mL−1 chloramphenicol and 50 μg mL−1 kanamycin). The initial inoculation concentration of the strains was OD730 = 0.2. Cell density was assessed by measuring the optiamal density at 730 nm using a Cary 60 UV spectrophotometer (Agilent, Germany). All strains were stored in BG-11 medium with 25% glycerol at − 80 °C.

For high density cultivation, the HDC 6.10B Starter Kit (CellDEG, Germany) was used to cultivate Syn6803 cells according to the manufacturer’s instructions. The exact composition of the nutrient-enriched media, referred to as CD medium, is available in protocols.io (dx.doi.org/https://doi.org/10.17504/protocols.io.2bxgapn) [24]. For taxanes production tests, HDC cultures were shaken at 200 rpm under continuous illumination of 100 µmol m−2 s−1 at 30 °C for the first two days. Afterwards, the light intensity was increased to 200 µmol m−2 s−1. Simultaneously, 1 mL overlay of dodecane (D221104, Sigma, USA) was added to 9 mL of cell cultures, and 1 mL of CD medium containing appropriate antibiotics was supplemented every 4 days to compensate the water evaporation.

Plasmids construction

The host for plasmid construction was E. coli DH5α strain. All vectors, strains and primers used in this study are listed in Supplementary Tables S1 and S2. pBluescript II KS( +) purchased from Beijing qualityard biotechnology Co., Ltd. was used as a backbone. The target genes except for T5αH and CPR were amplified with standard PCR reactions using KOD-Plus-Neo high-fidelity DNA polymerase (KOD-401, Toyobo, Japan). The resulting fragments were assembled using Hieff Clone Plus Muiti One Step Cloning Kit (10912ES10, Yeasen, China) or via standard T4 DNA ligation reaction (15224025, Thermo Fisher Scientific, USA) according to the manufacturer’s instructions. All primers and genes were commercially synthesized by Suzhou Genewiz Inc. All plasmids constructed were verified by sanger sequencing (Genewiz Inc., China) and subsequently used for natural transformation.

Chromosomal integration

For transformation of constructed vectors, 2 mL of fresh WT cells with an optical density (OD730) ranging from 0.4 to 0.6 were harvested by centrifugation at 2000g for 5 min and washed twice with ddH2O. The resulting cell pellet was resuspended in 400 μL fresh BG-11 medium. Subsequently, 5 μg plasmid was added to the cell mixture and incubated at 30 ℃ under constant light at an intensity of 25 µmol m−2 s−1 for 6–8 h. 200 μL of the cells was initially spread onto blank BG-11 agar plates for overnight incubation, then was transferred to the BG-11 agar plates supplemented with appropriate antibiotics (10 μg mL−1 spectinomycin, 10 μg·mL−1 chloramphenicol and 20 μg mL−1 kanamycin). After about 2 weeks, single clones appeared on the agars were screened by PCR verification to get positive clones. Several rounds of streaking were performed with increasing antibiotic concentrations to achieve complete segregation.

Pigments extraction and measurement

After 72 h cultivation, 1 mL of cells from HDC cultivator were centrifuged (15,000×g, 7 min, 4 °C) and then resuspended in 1 mL of pre-cooled methanol. After shaking and mixing, the cell mixtures stand at 4 °C for 30 min, then were centrifuged (15,000g, 10 min, 4 °C), the organic phase containing chlorophyll a and phycocyanin was aliquoted to a 96-well plate. The OD value was detected at different wavelengths measured with Synergy2 Microplate Reader (BioTeK, USA), including 625 nm, 678 nm, 720 nm, and 750 nm. The following formulae were used for quantifications [25]:

For chlorophyll content: 14.94 (OD678–OD750)–0.616 (OD720–OD750).

For phycocyanin content: 0.138445 (OD625–OD750)–0.0354 (OD678–OD750).

All values were corrected using the optical density of the medium or methanol before calculation, and three replicates of the test were performed for each sample.

Taxanes extraction and measurement

The productions of taxadiene and oxygenated taxanes were monitored by gas chromatography/mass spectrometry (GC/MS), using β-caryophyllene (75541, Sigma, USA) as an internal standard. The dodecane phase from different cell cultures was collected on the 10th day of cultivation. After centrifugation at 12,000 rpm at room temperature, the supernatant was diluted ten times in hexane containing 2.5 µg/mL β-Caryophyllene. Internal standards were prepared by diluting with hexane to concentrations between 0.625 and 50 μg/mL.

GC–MS analysis was conducted by using an Agilent Technologies 7890B GC System equipped with an Agilent 5977B mass-selective detector, a DB-5 ms capillary column (30 m × 250 μm × 0.25 μm, Agilent) was used for separation. Injections were made in pulsed splitless mode with an inlet temperature of 260 °C and pressure of 7.7758 psi. The GC method was following the method described by Biggs et al.[8]. The MS data were collected in the range of 45–550 m/z for each sample. Helium was used as the carrier gas at a flow rate of 1.0111 mL/min. The titer of taxanes were counted based on the peak area in comparison with the standard. GC-MS profiles depicting taxadiene and the four major oxygenated taxanes can be found in the Supplementary Figs. S1–S3.

RNA-sequencing (RNA-Seq)

Two engineered Syn6803 strains (DIGT, DIGT-P560) were grown up to the middle exponential phase and then collected by centrifugation. For each sample, three biological replicates were conducted to ensure reproducibility. Total RNA from 5 mL of HDC culture had been isolated using RNA isolater Total RNA Extraction Reagent (R401-01, Vazyme, China) following the modified method described by Schwarzkopf et al. [26]. Library construction and high-throughput sequencing were conducted by Shanghai Xu Ran Genomics Biotechnology Co. according to standard protocols. After sequencing, the raw data underwent initial processing, which involved the removal of adaptors, unknown bases, and low-quality reads. Subsequently, the filtered data were mapped to the Synechocystis reference genome in the NCBI Database (https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_000009725.1/).

For gene expression analysis, raw sequence counts of known genes were quantified using the StringTie software. The expression levels of known genes were calculated using FPKM (Fragments Per Kilobase of transcript per million fragments mapped) metrics. Differential gene expression (DEG) analysis between different sample groups was performed using DESeq2 software, with the absolute value of log2FC ≥ 1 and p ≤ 0.05 to define DEGs. To gain full-scale insights into gene function, TopGO software was utilized for Gene Ontology functional analysis, which involved functional annotation and categorization of the entire gene set and the target gene set. In addition, pathway functional annotation and categorization of these genes were performed using KEGG functions.

Results

Introducing upstream pathway of taxol biosynthesis into Synechocystis chassis for taxadiene production

In cyanobacteria, geranylgeranyl diphosphate (GGPP), a diterpene backbone comprising 20 carbon atoms, is originating from the condensation of three isopentenyl pyrophosphate (IPP) molecules and one dimethylallyl pyrophosphate (DMAPP) molecule via the methyl-D-erythritol 4-phosphate (MEP) pathway (Fig. 1). This pathway is prevalent in bacteria, cyanobacteria, green microalgae, and plant plastids [27]. As GGPP synthase (GGPPS) is recognized as the rate-limiting enzyme for GGPP production [28, 29], our initial step involved co-expression of heterologous GGPPS and TASY in the Syn6803 wild-type strain to divert the metabolic flux from the MEP pathway toward taxol synthesis. The ggpps and tasy genes were cloned from plasmid pGB259 (Supplementary Table S1) and co-overexpressed as an operon under the strong promoter Pcpc560 [30] and integrated into the psbA2 locus (slr1311) of Syn6803 genome [31, 32] through homologous recombination, creating the foundational engineered Strain TAXA-GT (Fig. 2A). To circumvent the solubility challenge of diterpenes in microbial cell cultures, engineered Syn6803 cells were cultivated in liquid medium overlaid with dodecane, a biocompatible organic solvent that mitigates cellular toxicity [33]. This two-phase culture system is commonly employed to recover secreted hydrophobic products from cells [34, 35]. According to the taxadiene profiling results, Strain TAXA-GT produced approximately 0.24 mg/L of taxadiene after 7-day cultivation in shake flasks, while no trace of taxadiene was detected in the WT strain (Fig. 2B, Supplementary Fig. S1).

Fig. 1
figure 1

Engineered biosynthetic pathway in Synechocystis sp. PCC 6803. In our model cyanobacteria, the taxol heterologous synthesis pathway starts from the isoprenoid universal precursors IPP and DMAPP (provided by the native MEP pathway) to generate GGPP, then taxadiene and small amounts of its isomer iso-taxadiene. Under the catalysis of P450 oxidase (T5αH) and its cognate reductase (CPR), taxadiene is oxidized into a series of oxygenated taxanes (as shown in orange square), of which the specific product T5α-ol will undergo multiple rounds of stereospecific oxidations, acylations, and benzoylation to form Taxol (as shown in black dashed square). Genes highlighted in red represent heterologous genes that have been overexpressed through genome integration. G3P: glyceraldehyde-3 phosphate, DXP: deoxyxylulose 5-phosphate, MEP: methylerythritol 4-phosphate, HMBPP: (E)-4-Hydroxy-3-methyl-but-2-enyl pyrophosphate, IPP: isopentenyl diphosphate, DMAPP: dimethylallyl diphosphate, DXS: 1-deoxyxylulose-5-phosphate synthase, IspA: farnesyl phosphate synthase, GGPPS: geranylgeranyl diphosphate synthase, TASY: taxadiene synthase, T5αH: taxadiene-5α-hydroxylase, CPR: cytochrome P450 reductase

Fig. 2
figure 2

Metabolic engineering of Synechocystis sp. PCC 6803 for photosynthetic taxadiene production from CO2. A Schematic diagrams of the construction of Syn6803 strains that produce taxadiene. The heterologous genes dxs and ispA derived from E. coli were introduced to a neutral site slr3031, while ggpps and tasy with truncation of plastid transport peptides were introduced to another site slr1311. B Taxadiene titers of Strains TAXA-GT and DIGT in shake flasks. C Biomass measured as absorbance at 730 nm of WT and engineered strains. All data are the mean ± standard deviation (SD) from triplicated shake flasks

While the introduction of the mevalonate (MVA) pathway is often preferred over engineering the native MEP pathway in E. coli to enhance terpenoid production [34], the latter is more suitable for manipulation to benifit terpenoid biosynthesis in cyanobacteria. This is due to the relatively larger pools of glyceraldehyde-3-phosphate (GAP) and pyruvate compared to acetyl-CoA in cyanobacteria [36]. 1-deoxyxylulose-5-phosphate synthase (DXS) is generally considered as the key rate-limiting enzyme of the MEP pathway in various hosts [37, 38], while farnesyl phosphate synthase (IspA) is responsible for catalyzing the condensation of IPP and DMAPP to generate the terpenoid backbone molecule farnesyl diphosphate (FPP), plays a pivotal role in directing precursors toward the terpenoid synthesis pathway [21, 39] (Fig. 1). Therefore, the dxs and ispA genes cloned from E. coli were inserted between slr2030 and slr2031 (refered to as slr3031) [40] of Strain TAXA-GT to generate Strain DIGT (Fig. 2A). After 7-day cultivation, Strain DIGT accumulated up to 0.96 mg/L of taxadiene (Fig. 2B), which was four times the yield of Strain TAXA-GT. This result demonstrated that co-overexpression of DXS and IspA was beneficial for taxadiene biosynthesis by increasing the metabolic flux through the MEP pathway.

To assess whether the introduction of heterologous genes hindered cell growth, the cell densities of the WT, TAXA-GT, and DIGT strains were monitored for 12 d under normal conditions. The growth pattern of engineered strains exhibited a slghtly sluggish phenotype compared with WT (Fig. 2C). Given that alleviation of cytotoxicity caused by engineering the MEP pathway by overexpression of terpene synthase has been observed in both E. coli [34] and cyanobacteria [21], the growth defect of the taxadiene-producing strains suggests that maintenance of the metabolic flux balance between the MEP and terpenoid pathway is critical to preserve the robustness of the microbial host.

Promoting the titer of taxadiene by high density cultivation

High-density cultivation (HDC) represents an efficient cultivation system designed specifically for cyanobacteria [24, 41]. By leveraging membrane-mediated CO2 supply technology, optimized nutrient-enriched CD medium, and intense illumination, HDC enables rapid and sustainable biomass accumulation. We harnessed this system to enhance biomass production and unlock its potential for generating ample taxol precursors in our engineered Syn6803 strains.

Strains WT, TAXA-GT, and DIGT were cultured in nutrient-rich CD medium within the HDC 6.10B device (CellDEG, Fig. 3A). In this setup, a concentrated KHCO3/K₂CO₃ buffer (the molar ratio of KHCO3 and K₂CO₃ was 9:1) in the container base continuously delivered 90 mbar of CO2-partial pressure, where the gas diffused through the hydrophobic membranes into turbulent cell suspensions. To maintain a homogenous culture, 9 mL of each culture was supplemented with 10% [v/v] dodecane overlay and agitated at 200 rpm under constant illumination with increasing light intensities.

Fig. 3
figure 3

Growth characteristics and taxadiene production in High Density Cultivation (HDC) cultures. A Schematic representation of the HDC 6.10B system (CellDEG). B Taxadiene titers of Strains TAXA-GT and DIGT in HDC cultures. C Biomass measured as absorbance at 730 nm of WT and engineered strains. D and E Levels of chlorophyll a and phycocyanin contents from engineered cyanobacteria. All data are the mean ± standard deviation (SD) from triplicated HDC vessels

By in situ extration and GC–MS analysis, taxadiene titers of 1.48 mg/L for the TAXA-GT strain and 2.94 mg/L for the DIGT strain were discovered in HDC media (Fig. 3B), which were increased by 5 times and two times, respectively, compared with shake flask culture. Cell density (OD730) and photosynthetic pigments content were recorded daily over a 10-day period in three independent cultivation runs. As Fig. 3C shown, the growth curve exhibited a similar trend to that of shake flask culture, with growth gradually slowing down after 10-day cultivation. The difference was that the growth rates of the engineered strains were close to that of WT before the sixth day, but then gradually lagged behind. Normally, the photosynthetic activity of cyanobacteria can be measured sideways by the content of photosynthetic pigment. There was a distinct difference in the chlorophyll a content between the WT and the engineered strain, especially for Strain DIGT which was almost half that of the WT (Fig. 3D). Phycocyanin is an important photosynthetic accessory pigment that transfers absorbed light energy from phycobilisomes to the chlorophyll-based photosynthetic apparatus [42]. The phycocyanin content in WT strain showed a steady upward trend with culture time, while that of the engineered strains fluctuated greatly and suddenly dropped after the 8th day (Fig. 3E), indicating that the photosynthetic efficiency of taxadiene producing strains may be disrupted in the later phase. Despite this, compared to shake flask, the growth rate of the engineered strains under HDC condition was not significantly different from that of the WT, suggesting that high light may be able to compensate to a certain extent for the adverse effects of insufficient photosynthetic pigments on photosynthesis. Given that the elevated biomass concentration will improve the efficiency of in situ extraction of terpenes, the HDC strategy was applied in subsequent experiments to introduce the downstream oxidation pathway of taxadiene.

Combinatorial optimization of downstream pathway to produce oxygenated taxanes in Synechocysis

A pivotal step in engineering microbes for taxol production is the establishment of cytochrome P450s-based oxidation within the cellular environment. In previous studies conducted using E. coli chassis [7, 8], considering the absence of compartmentalization structures such as the endoplasmic reticulum in prokaryotes, the N-terminal domains of T5αH and CPR were engineered to ensure their proper function in cytoplasm. In this context, we introduced the intact coding sequences of T5αH and CPR (Supplementary Table S3) directly into the taxadiene-producing strain DIGT by employing codon optimization to enhance the expression efficiency of these heterologous proteins (Fig. 4A). The chromosomal loci slr0168 [43] was used as the neutral site for expressing the target genes. It is well known that gene expression can be fine-tuned at the transcriptional level through promoter optimization or at the translational level by modifying ribosome-binding sites (RBSs). To facilitate genetic manipulation of Syn6803, a previous study established an RBS library consisting of 12 modified RBSs [44]. Notably, RBSv33 with a disrupted AT-rich region exhibited the highest expression level of the reporter gene, reaching 254% that of RBSv4, a relative strong RBS that was used to regulate the expression of ethylene forming enzyme in Syn6803 [45]. Accordingly, we recruited RBSv4 and RBSv33 into our P450 expression module to investigate whether T5αH or CPR plays a more critical role in taxane oxygenation under the genetic background of Syn6803. Strains DIGT-TC, DIGT-mT and DIGT-mC harbored different combinations of RBSs under the regulation of a chimeric promoter PpsbA* [44], while the expression of P450 module in Strain DIGT-P560 was under the control of Pcpc560 (Fig. 4A). In addition, in the purpose of steering the metabolic flux toward taxane downstream pathway, T5αH was physically linked to its native CPR via a short flexible linker to create a chimera, resulting in Strain DIGT-TLC (Fig. 4A).

Fig. 4
figure 4

Combinatorial optimization of T5αH and CPR expression to produce oxygenated taxanes in Syn6803. A Schematic diagrams of the construction of Syn6803 strains that produce oxygenated taxanes. The heterologous genes t5αh and cpr were codon-optimized and introduced together to the neutral site slr0168 of taxadiene-producing strain DIGT. Five different plasmids were constructed with different combination of promoters, RBSs and linker to regulate the expression levels between t5αh and cpr. B The distribution of oxygenated taxanes produced in downstream strains. T5α-ol: Taxadiene-5α-ol. C Total oxygenated taxanes production in downstream strains. D Growth rate of DIGT and all downstream strains under high density cultivation, OD730 values were measured daily. E and F Levels of chlorophyll a and phycocyanin contents from engineered cyanobacteria. All data are the mean ± standard deviation (SD) from the triplicated culture, and  in C represents statistical significance as indicated by Student’s t-test. p-value of < 0.05; p-value of < 0.01; p-value of < 0.001; p-value of < 0.0001

The product distributions of all downstream strains were determined by GC–MS analysis of the dodecane extracts. As shown in Supplementary Figs. S1 and S2, multiple taxane products were identified, including taxadiene and four major monohydroxylated compounds (parent ion at 288 m/z). Through comparison with gas chromatographs and mass spectra found in the literature [11], the peaks at 21.095, 21.378, 21.805, and 22.082 min were identified as iso-OCT, Diterpenoid 1, OCT, and T5α-ol, respectively (Supplementary Fig. S3). OCT and iso-OCT are considered as the major products in most T5αH-expressing hosts, such as E. coli, S. cerevisiae, and Yarrowia lipolytica [9, 46]. However, Diterpenoid 1, an isomer of T5α-ol, was emerging as the predominant oxygenated product in all downstream Syn6803 strains (Fig. 4B), which was consistent with a recent work conducted in S. cerevisiae [11]. This result highlighted the significant influence of host organisms and cultivation conditions on the selectivity of T5αH. Simultaneously, T5α-ol was the second most abundant compound in all downstream strains, accounting for approximately 20–27% of the total products (Supplementary Table S4).

Considering the overall productivity of oxygenated taxanes, Strain DIGT-P560, featuring the “super strong” Pcpc560 promoter, produced the highest titer of total oxygenated taxanes as 17.43 mg/L (Fig. 4C). Conversely, Strain DIGT-TC harboring the regulatory module of relatively low strength exhibited the lowest yield of 6.98 mg/L. It is noteworthy that a substantial amount of taxadiene was left over in the organic extracts of DIGT-TC, DIGT-mC, and DIGT-TLC (Supplementary Figs. S1, S2), suggesting the inefficiency of P450-mediated oxygenation toward taxadiene in those strains. In addition, Strain DIGT-mT, which might possess relatively higher T5αH expression than CPR resulted from the regulation of RBSv33 on T5αH, had produced more oxygenated products than Strain DIGT-mC that exhibited higher CPR expression (11.49 mg/L vs 8.71 mg/L). Taken together, it is plausible that a higher T5αH expression level could have advantageous effect on taxadiene oxygenation in Syn6803.

Compared with the basal Strain DIGT, the growth rates of some downstream strains transformed with T5αH/CPR were slightly lower (Fig. 4D), indicating that the expression of heterologous P450 enzyme may increase the metabolic burden of the cells. Interestingly, in terms of the synthesis of photosynthetic pigments, most downstream strains performed better than the upstream strain DIGT (Fig. 4E, F). It seems that the introduction of P450 enzyme improved the photosynthetic efficiency of the engineered strains to a certain extent. In conclusion, by manipulating the taxadiene downstream pathway via T5αH/CPR modular modification, we successfully developed an optimally engineered strain, DIGT-P560, which produced 17.43 mg/L of oxygenated taxanes and 4.32 mg/L of T5α-ol. While the rate of biomass and taxane accumulation may not currently rival that of heterotrophic microorganisms, efficient expression of the heterologous P450 enzyme led to a relatively high proportion of T5α-ol in our recombinant strains. This characteristic underscores the advantages of cyanobacterial chassis for further in vivo modification of the paclitaxel biosynthetic pathway when compared to other microbial hosts.

Transcriptomic analysis for exploring the comprehensive effects of heterologous P450 expression on a whole-genome dimension of Synechocystis

To elucidate the potential mechanisms of unnatural diterponoid biosynthesis in cyanobacteria, a comprehensive whole-genome transcriptomic analysis was conducted to assess gene expression changes in DIGT and DIGT-P560 strains. Strains for RNA-seq were collected during the exponential growth phase. The Pearson Correlation Coefficient heat map shows the differences between samples within each group, allowing us to visually evaluate the similarities between biological replicates in the same group, as well as the differences between samples in different groups. As shown in Fig. 5A, there were significant differences between samples from DIGT and DIGT-P560 groups. The samples within the DIGT group were highly similar, while there were certain differences between the three replicates within the DIGT-P560 group. In comparison to the DIGT group, 1367 differentially expressed genes (DEGs) with Log2 of fold change (FC) ≥ 1 or ≤ -1 and p ≤ 0.05 were found in DIGT-P560 (Fig. 5B). With assignment to KEGG orthology database, apart from hypothetical and other functional genes, the rest DEGs were mainly categorized into 11 enriched pathways, including energy metabolism, metabolism of cofactors and vitamins, translation, amino acid metabolism, et al. (Fig. 5C). As expression of the majority of genes related to translation was significantly downregulated (Supplementary Data Set S1), the decreased levels of ribosomal proteins could potentially lead to disruption of ribosome assembly, which in turn may further diminish the translation efficiency of mature mRNAs in DIGT-P560 strains. Therefore, there were more downregulated genes than upregulated genes in DIGT-P560 (732 vs 635).

Fig. 5
figure 5

Differentially expressed genes (DEGs) analysis in downstream Strain DIGT-P560 relative to upstream Strain DIGT. A Pearson correlation heat map of all samples. B Volcano plots of differentially expressed genes in DIGT and DIGT-P560. Grey plot: Not DEGs; Blue plot: Up-regulated DEGs; Purple plot: Down-regulated DEGs. C DEGs number of up-regulated and down-regulated genes in annotated KEGG metabolic pathways. Light brown column: Up-regulated DEGs; Navy blue column: Down-regulated DEGs

For DIGT-P560, energy metabolism pathway had the largest number of DEGs, the expression of 80 related genes was significantly changed, of which nearly two-thirds were downregulated (Fig. 5C). The down-regulated genes were mainly concentrated in the photosynthesis, antenna proteins, and oxidative phosphorylation pathways (Fig. 6, Supplementary Data Set S1). This may be related to the fact that heterologous P450 enzymes were located on the thylakoid membrane, thereby affecting the expression of other photosystem proteins and F-type ATP synthase on the membrane. Meanwhile, the expression levels of genes related to biological processes such as nitrogen metabolism and sulfur metabolism showed an upward trend (Supplementary Data Set S1). It is noteworthy that the most significantly up-regulated genes in energy metabolism were nifJ (log2FC = 4.86) and ppsA (log2FC = 6.00), where nifJ is an oxidoreductase required for the transfer of electrons from pyruvate to flavodoxin, and ppsA is responsible for catalyzing the phosphorylation reaction of pyruvate [47]. The substantial up-regulation of these two genes suggested that the intracellular pyruvate content of DIGT-P560 could be higher than DIGT.

Fig. 6
figure 6

Comparative transcriptomic analysis of Strain DIGT-P560 and DIGT. Transcript abundance changes of the unigene involved in the TCA cycle (yellow panel), Glycolysis (navy blue panel), Gluconeogenesis (light blue panel), ED pathway (flesh panel), CBB cycle (green panel) and the photosynthetic chain (bottom panel). Up-regulated and down-regulated genes were labeled in red and green. TCA: tricarboxylic acid; ED: Entner-Doudoroff; CBB: Calvin–Benson–Bassham

In cyanobacteria, carbohydrate metabolism typically includes central metabolic pathways, including glycolysis, the pentose phosphate pathway, the Calvin–Benson–Bassham (CBB) cycle, the Entner–Doudoroff (ED) pathway, and the tricarboxylic acid (TCA) cycle [47]. For Strain DIGT-P560, there were 44 DEGs belonging to carbohydrate metabolism, half of which were up-regulated (Fig. 5C) and mainly concentrated in the TCA cycle, pentose phosphate pathway and pyruvate metabolism (Supplementary Data Set S1). The TCA cycle, which is responsible for providing building blocks for growth, typically maintains a low flux rate under photoautotrophic conditions [48]. However, in DIGT-P560, five related genes were significantly upregulated, including gltA, 2-OGDC, SSADH, sdhA, and citH (Fig. 6, Supplementary Data Set S1), suggesting the metabolic flux might be redirected to the TCA cycle to compensate the NADPH consumed by P450s, thereby maintaining the overall energy balance. In addition, the enhanced flux through the TCA cycle may lead to the increased synthesis of other essential compounds, such as amino acids and fatty acids, as evidenced by the higher number of upregulated genes related to amino acid and lipid metabolism in DIGT-P560 (Fig. 5C, Supplementary Data Set S1). Notably, lipids serve multiple crucial functions in cyanobacteria, including energy storage and structural support for cell membranes and thylakoid membranes [49]. In response to the challenges posed by the introduction of membrane-anchored proteins T5αH and CPR, DIGT-P560 may enhance its lipid synthesis capacity to maintain the normal function of thylakoids and provide a lipid-rich environment conducive to full activity of the P450 enzyme [6].

When exploring the expression level of genes in signal transduction and membrane transport pathways, most genes were significantly upregulated in DIGT-P560 (Fig. 5C). Specifically, 22 out of 29 genes encoding two-component system proteins and 17 out of 26 genes encoding ATP-binding cassette (ABC) transporters showed upregulation (Supplementary Data Set S1). The two-component system is recognized for its tight regulation of ABC transporters, which play crucial roles in nutrient transport and cytoplasmic pH regulation [50, 51]. Both systems have the capacity to sense biotic and abiotic stresses as well as substrates such as peptides, amino acids, sugars, and antibiotics. Therefore, the activation of such genes in DIGT-P560 suggested that there were increased transmembrane material transport activities aimed at alleviating cellular stress in the taxane-producing strains. In addition, the expression levels of most genes related to DNA replication and repair, protein folding, sorting, and degradation were significantly higher in DIGT-P560 than in DIGT (Fig. 5C, Supplementary Data Set S1). This implied that cells under greater stress would have increased demands for restoration through processes such as DNA repair, homologous recombination, and RNA degradation.

Discussion

As one of the most widely administered anticancer agents, the economical biosynthesis of taxol has remained a significant challenge for decades. While certain achievements have been made through heterologous expression of complex biosynthetic pathways in fast-growing microbial hosts, such as E. coli, B. subtilis, and S. cerevisiae, there has been no research on taxol biosynthesis by cyanobacteria despite their advantages of low culture cost, rapid growth, and ability to utilize CO2 as a carbon source. Therefore, we conducted a series of attempts in cyanobacteria model strain Synechocystis sp. PCC 6803 to explore a new platform for taxol biosynthesis.

The final conversion of taxadiene to taxol requires effective hydroxylation at eight specific sites: C1, C2, C4, C5, C7, C9, C10, and C13. The initial step involved introducing a hydroxyl group at the C5 position of taxadiene and relocating the double bond from C4(5) to C4(20) through the action of T5αH to yield T5α-ol [5]. This reaction is considered as the rate-limiting step in the entire oxygenation process. The native T5αH of T. cuspidata possesses a putative N-terminal sequence of 137 residues that is cleaved during maturation in plastids [52]. Truncating and modifying the N-terminal regions of T5αH have been found to promote the production of more active and soluble proteins in E. coli, addressing the challenges posed by the absence of internal membranes [7, 8]. In this study, the successful expression of intact T5αH and its reductase partner CPR clearly demonstrates the adaptability of the cyanobacterial chassis for the expression of N-terminally anchored cytochrome P450s, which will greatly facilitate further reconstruction of the complex taxol metabolic pathways in cyanobacterial cell factories.

In addition, it has been hypothesized that maintaining a balance between P450 and CPR is crucial for efficient catalysis of taxane substrates. Previous reports have shown that physically fusing P450/CPR does not provide advantages for protein expression or product synthesis [8], whereas the ratio of P450 to CPR in plant membranes is approximately 15:1 [53]. The oxygenated taxane profiles observed in our engineered strains were consistent with these findings. Strain DIGT-TLC exhibited the lowest oxygenated taxane productivity owing to the unsatisfied activity of the T5αH-CPR fusion protein. The concentration of total oxygenated taxanes, from highest to lowest, was observed in DIGT-P560, DIGT-mT, DIGT-mC, DIGT-TC, and DIGT-TLC, while the proportion of T5α-ol from highest to lowest was observed in DIGT-TLC, DIGT-mT, DIGT-mC, DIGT-P560, and DIGT-TC. These results indicate that manipulating the expression pattern of the T5αH/CPR combination at the translational or transcriptional level can significantly influence the productivity and distribution of downstream oxygenated products.

The activity of cytochrome P450s in the taxol pathway appears to be limited in prokaryotic hosts, as the production of taxadiene oxides significantly decreases when T5αH is introduced [7]. Although the delicate P450 module optimizations contribute to the effective enhancement of oxygenated taxanes in engineered E. coli, the titer of downstream products can not surpass their substrate, taxadiene [8]. Notably, in our study with Syn6803, even if the basal strain DIGT produced only 3 mg/L of taxadiene, the downstream strains transformed with the T5αH/CPR module produced up to 17 mg/L of oxygenated taxanes, significantly surpassing the substrate concentration. The transcriptomic analysis toward DIGT and DIGT-P560 indicated that, majority of downregulated DEGs were categorized into energy metabolism and translation, suggesting that the introduction of heterologous P450s could interfere with the photosystems and normal biological process of Syn6803. However, at the same time, certain genes were significantly upreglulated in DIGT-P560, particularly those related to the TCA cycle (Supplementary Data Set S1). Activation of such pathways may benefit the production of glyceraldehyde-3-P and pyruvate [47], providing more substrates for the MEP pathway to increase the GGPP pool for more efficient diterpene synthesis.

Although the photosystems of our engineered Syn6803 strains were impaired owing to a significant reduction in the expression of related genes, the differences in growth rates between the wild-type and recombinant strains were relatively minor. How did these recombinant strains survive with their compromised photosystems? Phylloquinone is a prenylated naphthoquinone that plays a crucial role as an electron carrier in photosystem I and an electron acceptor for the formation of protein disulfide bonds [54]. Cyanobacterial strains lacking genes involved in phylloquinone biosynthesis, such as menA, menB, menD, and menE, have been shown to have slower photoautotrophic growth rates than the wild type [55, 56]. Plastoquinone (PQ), another important quinol, is located on the thylakoid membrane and forms a PQ pool with the reduced plastoquinone (PQH2) within the thylakoids. In cyanobacteria, PQ primarily mediates electron transfer between photosystem II and cytochrome b6f complex [57]. Our RNA-seq results revealed that genes related to phylloquinone (menB, menC, menE) and plastoquinone biosynthesis (ubiC) were significantly upregulated in DIGT-P560 (Supplementary Data Set S1). This discrepancy indicated a potential enhancement in the electron transport rate through the linear and cyclic photosynthetic electron transport chains. Meanwhile, several NAD biosynthetic genes were upregulated, such as pntA/pntB encoding pyridine nucleotide transhydrogenase (Supplementary Data Set S1), which catalyses electron transfer between NAD+ and NADP+ and regulates the NAD+/NADP+ ratio [58], thereby could provide sufficient reducing power for light-powered P450 activity. Intriguingly, previous studies have illustrated that P450s can serve as artificial electron sinks for excess electrons derived from photosynthesis and can increase the maximum rate of photosynthetic electron flow in cyanobacteria [13, 17, 59, 60]. It has been proved that cellular ATP and NADPH concentrations are directly proportional to P450 acitivity [13]. Our transcriptomic data may provide evidence of a potential mechanism for how the heterologous P450 electron sink influences photosynthetic electron transport and energy metabolism in cyanobacteria.

Conclusion

In summary, our integrated approach led to the successful construction of a taxane-producing platform in photoautotrophic cyanobacteria. Given the flexibility of the metabolic network in this chassis, it is reasonable to expect that further improvements in the productivity of valuable oxygenated taxanes can be achieved through extensive research in the near future. For example, inducible promoters may be more suitable for heterologous cytochrome P450s enzymes to mitigate the metabolic burdens on the photosystems in cyanobacteria. Although the total biosynthesis of taxol in microbial cell factories remains a significant challenge, our work lays a solid foundation for the direct transformation of CO2 into taxol precursors in cyanobacteria in vivo. We hope that this will promote the widespread application of biosynthesis of various complex and valuable natural products in photosynthetic microorganisms.

Data availability

No data sets were generated or analysed during the current study.

References

  1. Schiff PB, Fant J, Horwitz SB. Promotion of microtubule assembly in vitro by taxol. Nature. 1979;277(5698):665–7.

    Article  CAS  PubMed  Google Scholar 

  2. Wani MC, Taylor HL, Wall ME, Coggon P, McPhail AT. Plant antitumor agents. VI. The isolation and structure of taxol, a novel antileukemic and antitumor agent from Taxus brevifolia. J Am Chem Soc. 1971;93(9):2325–7.

    Article  CAS  PubMed  Google Scholar 

  3. Frense D. Taxanes: perspectives for biotechnological production. Appl Microbiol Biotechnol. 2007;73(6):1233–40.

    Article  CAS  PubMed  Google Scholar 

  4. Wang T, Li L, Zhuang W, Zhang F, Shu X, Wang N, Wang Z. Recent research progress in taxol biosynthetic pathway and acylation reactions mediated by taxus acyltransferases. Molecules. 2021;26(10):2855.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Jennewein S, Long RM, Williams RM, Croteau R. Cytochrome p450 taxadiene 5alpha-hydroxylase, a mechanistically unusual monooxygenase catalyzing the first oxygenation step of taxol biosynthesis. Chem Biol. 2004;11(3):379–87.

    Article  CAS  PubMed  Google Scholar 

  6. Munro AW, Girvan HM, Mason AE, Dunford AJ, McLean KJ. What makes a P450 tick? Trends Biochem Sci. 2013;38(3):140–50.

    Article  CAS  PubMed  Google Scholar 

  7. Ajikumar PK, Xiao WH, Tyo KE, Wang Y, Simeon F, Leonard E, Mucha O, Phon TH, Pfeifer B, Stephanopoulos G. Isoprenoid pathway optimization for Taxol precursor overproduction in Escherichia coli. Science. 2010;330(6000):70–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Biggs BW, Lim CG, Sagliani K, Shankar S, Stephanopoulos G, De Mey M, Ajikumar PK. Overcoming heterologous protein interdependency to optimize P450-mediated Taxol precursor synthesis in Escherichia coli. Proc Natl Acad Sci USA. 2016;113(12):3209–14.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Sagwan-Barkdoll L, Anterola AM. Taxadiene-5-ol is a minor product of CYP725A4 when expressed in Escherichia coli. Biotechnol Appl Bioc. 2018;65(3):294–305.

    Article  CAS  Google Scholar 

  10. Dejong JM, Liu Y, Bollon AP, Long RM, Jennewein S, Williams D, Croteau RB. Genetic engineering of taxol biosynthetic genes in Saccharomyces cerevisiae. Biotechnol Bioeng. 2006;93(2):212–24.

    Article  CAS  PubMed  Google Scholar 

  11. Walls LE, Malci K, Nowrouzi B, Li RA, d’Espaux L, Wong J, Dennis JA, Semiao AJC, Wallace S, Martinez JL, et al. Optimizing the biosynthesis of oxygenated and acetylated Taxol precursors in Saccharomyces cerevisiae using advanced bioprocessing strategies. Biotechnol Bioeng. 2021;118(1):279–93.

    Article  CAS  PubMed  Google Scholar 

  12. Walls LE, Martinez JL, Del Rio Chanona EA, Rios-Solis L. Definitive screening accelerates Taxol biosynthetic pathway optimization and scale up in Saccharomyces cerevisiae cell factories. Biotechnol J. 2022;17(1): e2100414.

    Article  PubMed  Google Scholar 

  13. Berepiki A, Gittins JR, Moore CM, Bibby TS. Rational engineering of photosynthetic electron flux enhances light-powered cytochrome P450 activity. Synth Biol. 2018;3(1): yay009.

    Article  Google Scholar 

  14. Lassen LM, Nielsen AZ, Ziersen B, Gnanasekaran T, Moller BL, Jensen PE. Redirecting photosynthetic electron flow into light-driven synthesis of alternative products including high-value bioactive natural compounds. ACS Synth Biol. 2014;3(1):1–12.

    Article  CAS  PubMed  Google Scholar 

  15. Parmar A, Singh NK, Pandey A, Gnansounou E, Madamwar D. Cyanobacteria and microalgae: a positive prospect for biofuels. Bioresour Technol. 2011;102(22):10163–72.

    Article  CAS  PubMed  Google Scholar 

  16. Ruiz-Sola MA, Coman D, Beck G, Barja MV, Colinas M, Graf A, Welsch R, Rutimann P, Buhlmann P, Bigler L, et al. Arabidopsis GERANYLGERANYL DIPHOSPHATE SYNTHASE 11 is a hub isozyme required for the production of most photosynthesis-related isoprenoids. New Phytol. 2016;209(1):252–64.

    Article  CAS  PubMed  Google Scholar 

  17. Berepiki A, Hitchcock A, Moore CM, Bibby TS. Tapping the unused potential of photosynthesis with a heterologous electron sink. ACS Synth Biol. 2016;5(12):1369–75.

    Article  CAS  PubMed  Google Scholar 

  18. Jensen K, Jensen PE, Moller BL. Light-driven cytochrome p450 hydroxylations. ACS Chem Biol. 2011;6(6):533–9.

    Article  CAS  PubMed  Google Scholar 

  19. Nielsen AZ, Mellor SB, Vavitsas K, Wlodarczyk AJ, Gnanasekaran T, de Perestrello Ramos HJM, King BC, Bakowski K, Jensen PE. Extending the biosynthetic repertoires of cyanobacteria and chloroplasts. Plant J. 2016;87(1):87–102.

    Article  CAS  PubMed  Google Scholar 

  20. Vickers CE, Bongers M, Liu Q, Delatte T, Bouwmeester H. Metabolic engineering of volatile isoprenoids in plants and microbes. Plant Cell Environ. 2014;37(8):1753–75.

    Article  CAS  PubMed  Google Scholar 

  21. Choi SY, Lee HJ, Choi J, Kim J, Sim SJ, Um Y, Kim Y, Lee TS, Keasling JD, Woo HM. Photosynthetic conversion of CO2 to farnesyl diphosphate-derived phytochemicals (amorpha-4,11-diene and squalene) by engineered cyanobacteria. Biotechnol Biofuels. 2016;9:202.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Diao J, Song X, Zhang L, Cui J, Chen L, Zhang W. Tailoring cyanobacteria as a new platform for highly efficient synthesis of astaxanthin. Metab Eng. 2020;61:275–87.

    Article  CAS  PubMed  Google Scholar 

  23. Rosmarie Rippka JD, John BW, Herdman M, Stanier RY. Generic assignments, strain histories and properties of pure cultures of cyanobacteria. Microbiology. 1979;111(1):61.

    Google Scholar 

  24. Dienst D, Wichmann J, Mantovani O, Rodrigues JS, Lindberg P. High density cultivation for efficient sesquiterpenoid biosynthesis in Synechocystis sp. PCC 6803. Sci Rep. 2020;10(1):5932.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Khan RI, Wang Y, Afrin S, Wang B, Liu Y, Zhang X, Chen L, Zhang W, He L, Ma G. Transcriptional regulator PrqR plays a negative role in glucose metabolism and oxidative stress acclimation in Synechocystis sp. PCC 6803. Sci Rep. 2016;6:32507.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Schwarzkopf M, Yoo YC, Huckelhoven R, Park YM, Proels RK. Cyanobacterial phytochrome2 regulates the heterotrophic metabolism and has a function in the heat and high-light stress response. Plant Physiol. 2014;164(4):2157–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Lichtenthaler HK. Non-mevalonate isoprenoid biosynthesis: enzymes, genes and inhibitors. Biochem Soc Trans. 2000;28(6):785–9.

    Article  CAS  PubMed  Google Scholar 

  28. Jiang M, Stephanopoulos G, Pfeifer BA. Toward biosynthetic design and implementation of Escherichia coli-derived paclitaxel and other heterologous polyisoprene compounds. Appl Environ Microbiol. 2012;78(8):2497–504.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Bai C, Capell T, Berman J, Medina V, Sandmann G, Christou P, Zhu C. Bottlenecks in carotenoid biosynthesis and accumulation in rice endosperm are influenced by the precursor-product balance. Plant Biotechnol J. 2016;14(1):195–205.

    Article  CAS  PubMed  Google Scholar 

  30. Zhou J, Zhang H, Meng H, Zhu Y, Bao G, Zhang Y, Li Y, Ma Y. Discovery of a super-strong promoter enables efficient production of heterologous proteins in cyanobacteria. Sci Rep. 2014;4:4500.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Kuchmina E, Wallner T, Kryazhov S, Zinchenko VV, Wilde A. An expression system for regulated protein production in Synechocystis sp. PCC 6803 and its application for construction of a conditional knockout of the ferrochelatase enzyme. J Biotechnol. 2012;162(1):75–80.

    Article  CAS  PubMed  Google Scholar 

  32. 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.

    Article  CAS  PubMed  Google Scholar 

  33. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 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.

    Article  CAS  PubMed  Google Scholar 

  35. 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.

    Article  PubMed  Google Scholar 

  36. 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;4:12.

    Google Scholar 

  37. 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.

    Article  CAS  PubMed  Google Scholar 

  38. 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.

    Article  CAS  PubMed  Google Scholar 

  39. Wang C, Zhou J, Jang HJ, Yoon SH, Kim JY, Lee SG, Choi ES, Kim SW. Engineered heterologous FPP synthases-mediated Z, E-FPP synthesis in E. coli. Metab Eng. 2013;18:53–9.

    Article  CAS  PubMed  Google Scholar 

  40. Sun T, Li S, Song X, Pei G, Diao J, Cui J, Shi M, Chen L, Zhang W. Re-direction of carbon flux to key precursor malonyl-CoA via artificial small RNAs in photosynthetic Synechocystis sp. PCC 6803. Biotechnol Biofuels. 2018;11:26.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Guljamow A, Kreische M, Ishida K, Liaimer A, Altermark B, Bahr L, Hertweck C, Ehwald R, Dittmann E. High-density cultivation of terrestrial nostoc strains leads to reprogramming of secondary metabolome. Appl Environ Microbiol. 2017;83(23): e01510.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Rockwell NC, Martin SS, Lagarias JC. Elucidating the origins of phycocyanobilin biosynthesis and phycobiliproteins. Proc Natl Acad Sci USA. 2023;120(17): e2300770120.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Kunert A, Hagemann M, Erdmann N. Construction of promoter probe vectors for Synechocystis sp. PCC 6803 using the light-emitting reporter systems Gfp and LuxAB. J Microbiol Methods. 2000;41(3):185–94.

    Article  CAS  PubMed  Google Scholar 

  44. Wang B, Eckert C, Maness PC, Yu J. A genetic toolbox for modulating the expression of heterologous genes in the cyanobacterium Synechocystis sp. PCC 6803. ACS Synth Biol. 2018;7(1):276–86.

    Article  CAS  PubMed  Google Scholar 

  45. Xiong W, Morgan JA, Ungerer J, Wang B, Maness PC, Yu JP. The plasticity of cyanobacterial metabolism supports direct CO2 conversion to ethylene. Nat Plants. 2015;1(5):5053.

    Article  Google Scholar 

  46. Edgar S, Zhou K, Qiao K, King JR, Simpson JH, Stephanopoulos G. Mechanistic insights into taxadiene epoxidation by taxadiene-5alpha-hydroxylase. ACS Chem Biol. 2016;11(2):460–9.

    Article  CAS  PubMed  Google Scholar 

  47. Mills LA, McCormick AJ, Lea-Smith DJ. Current knowledge and recent advances in understanding metabolism of the model cyanobacterium Synechocystis sp. PCC 6803. Biosci Rep. 2020;40(4): BSR20193325.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  48. Wan N, DeLorenzo DM, He L, You L, Immethun CM, Wang G, Baidoo EEK, Hollinshead W, Keasling JD, Moon TS, et al. Cyanobacterial carbon metabolism: fluxome plasticity and oxygen dependence. Biotechnol Bioeng. 2017;114(7):1593–602.

    Article  CAS  PubMed  Google Scholar 

  49. Tasaka Y, Gombos Z, Nishiyama Y, Mohanty P, Ohba T, Ohki K, Murata N. Targeted mutagenesis of acyl-lipid desaturases in Synechocystis: evidence for the important roles of polyunsaturated membrane lipids in growth, respiration and photosynthesis. EMBO J. 1996;15(23):6416–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Uchiyama J, Itagaki A, Ishikawa H, Tanaka Y, Kohga H, Nakahara A, Imaida A, Tahara H, Ohta H. Characterization of ABC transporter genes, sll1180, sll1181, and slr1270, involved in acid stress tolerance of Synechocystis sp. PCC 6803. Photosynth Res. 2019;139(1–3):325–35.

    Article  CAS  PubMed  Google Scholar 

  51. Kurdrid P, Yutthanasirikul R, Saree S, Senachak J, Saelee M, Hongsthong A. Hik28-dependent and Hik28-independent ABC transporters were revealed by proteome-wide analysis of DeltaHik28 under combined stress. BMC Mol Cell Biol. 2022;23(1):27.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Kitaoka N, Wu Y, Xu M, Peters RJ. Optimization of recombinant expression enables discovery of novel cytochrome P450 activity in rice diterpenoid biosynthesis. Appl Microbiol Biotechnol. 2015;99(18):7549–58.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Jensen K, Moller BL. Plant NADPH-cytochrome P450 oxidoreductases. Phytochemistry. 2010;71(2–3):132–41.

    Article  CAS  PubMed  Google Scholar 

  54. Basset GJ, Latimer S, Fatihi A, Soubeyrand E, Block A. Phylloquinone (Vitamin K1): occurrence, biosynthesis and functions. Mini Rev Med Chem. 2017;17(12):1028–38.

    Article  CAS  PubMed  Google Scholar 

  55. Wade Johnson T, Naithani S, Stewart C Jr, Zybailov B, Daniel Jones A, Golbeck JH, Chitnis PR. The menD and menE homologs code for 2-succinyl-6-hydroxyl-2,4-cyclohexadiene-1-carboxylate synthase and o-succinylbenzoic acid-CoA synthase in the phylloquinone biosynthetic pathway of Synechocystis sp. PCC 6803. Biochim Biophys Acta. 2003;1557(1–3):67–76.

    Article  CAS  PubMed  Google Scholar 

  56. Johnson TW, Shen G, Zybailov B, Kolling D, Reategui R, Beauparlant S, Vassiliev IR, Bryant DA, Jones AD, Golbeck JH, et al. Recruitment of a foreign quinone into the A(1) site of photosystem I: I. Genetic and physiological characterization of phylloquinone biosynthetic pathway mutants in Synechocystis sp. pcc 6803. J Biol Chem. 2000;275(12):8523–30.

    Article  CAS  PubMed  Google Scholar 

  57. Lea-Smith DJ, Bombelli P, Vasudevan R, Howe CJ. Photosynthetic, respiratory and extracellular electron transport pathways in cyanobacteria. Biochim Biophys Acta. 2016;1857(3):247–55.

    Article  CAS  PubMed  Google Scholar 

  58. Kamarainen J, Huokko T, Kreula S, Jones PR, Aro EM, Kallio P. Pyridine nucleotide transhydrogenase PntAB is essential for optimal growth and photosynthetic integrity under low-light mixotrophic conditions in Synechocystis sp. PCC 6803. New Phytol. 2017;214(1):194–204.

    Article  PubMed  Google Scholar 

  59. Torrado A, Connabeer HM, Rottig A, Pratt N, Baylay AJ, Terry MJ, Moore CM, Bibby TS. Directing cyanobacterial photosynthesis in a cytochrome c oxidase mutant using a heterologous electron sink. Plant Physiol. 2022;189(4):2554–66.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Santos-Merino M, Torrado A, Davis GA, Rottig A, Bibby TS, Kramer DM, Ducat DC. Improved photosynthetic capacity and photosystem I oxidation via heterologous metabolism engineering in cyanobacteria. Proc Natl Acad Sci USA. 2021;118(11): e2021523118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors thank Dr. Guangkai Bian for providing the plasmid pGB259.

Funding

This study was funded by National Natural Science Foundation of China, 31901022, National Key Research and Development Program of China, 2019YFA0904603, Natural Science Foundation of Shanghai, China, 21ZR1431700.

Author information

Authors and Affiliations

Authors

Contributions

Y.S.W., and G.M. conceived and designed the experiments. J.L.Z., Y.S.W., and Z.Y.C. performed the experiments. J.L.Z., Y.S.W., Z.Y.C., and G.M. analyzed the data. Y.S.W., Y.Y., L.H., T.G.L., and G.M. contributed reagents/materials/analysis tools. J.L.Z., Y.S.W. and G.M. wrote the paper. All authors reviewed the manuscript.

Corresponding author

Correspondence to Gang Ma.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhong, J., Wang, Y., Chen, Z. et al. Engineering cyanobacteria as a new platform for producing taxol precursors directly from carbon dioxide. Biotechnol Biofuels 17, 99 (2024). https://doi.org/10.1186/s13068-024-02555-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13068-024-02555-9

Keywords