Heterologous biosynthesis and manipulation of crocetin in Saccharomyces cerevisiae
Microbial Cell Factoriesvolume 16, Article number: 54 (2017)
Due to excellent performance in antitumor, antioxidation, antihypertension, antiatherosclerotic and antidepressant activities, crocetin, naturally exists in Crocus sativus L., has great potential applications in medical and food fields. Microbial production of crocetin has received increasing concern in recent years. However, only a patent from EVOVA Inc. and a report from Lou et al. have illustrated the feasibility of microbial biosynthesis of crocetin, but there was no specific titer data reported so far. Saccharomyces cerevisiae is generally regarded as food safety and productive host, and manipulation of key enzymes is critical to balance metabolic flux, consequently improve output. Therefore, to promote crocetin production in S. cerevisiae, all the key enzymes, such as CrtZ, CCD and ALD should be engineered combinatorially.
By introduction of heterologous CrtZ and CCD in existing β-carotene producing strain, crocetin biosynthesis was achieved successfully in S. cerevisiae. Compared to culturing at 30 °C, the crocetin production was improved to 223 μg/L at 20 °C. Moreover, an optimal CrtZ/CCD combination and a titer of 351 μg/L crocetin were obtained by combinatorial screening of CrtZs from nine species and four CCDs from Crocus. Then through screening of heterologous ALDs from Bixa orellana (Bix_ALD) and Synechocystis sp. PCC6803 (Syn_ALD) as well as endogenous ALD6, the crocetin titer was further enhanced by 1.8-folds after incorporating Syn_ALD. Finally a highest reported titer of 1219 μg/L at shake flask level was achieved by overexpression of CCD2 and Syn_ALD. Eventually, through fed-batch fermentation, the production of crocetin in 5-L bioreactor reached to 6278 μg/L, which is the highest crocetin titer reported in eukaryotic cell.
Saccharomyces cerevisiae was engineered to achieve crocetin production in this study. Through combinatorial manipulation of three key enzymes CrtZ, CCD and ALD in terms of screening enzymes sources and regulating protein expression level (reaction temperature and copy number), crocetin titer was stepwise improved by 129.4-fold (from 9.42 to 1219 μg/L) as compared to the starting strain. The highest crocetin titer (6278 μg/L) reported in microbes was achieved in 5-L bioreactors. This study provides a good insight into key enzyme manipulation involved in serial reactions for microbial overproduction of desired compounds with complex structure.
Crocetin, a kind of carotenoid existing in Crocus sativus L. , has great potential medical applications due to various pharmacological activities, such as antitumor [2, 3], antioxidation , antihypertension , antiatherosclerotic  and antidepressant . Additionally, crocetin can be also used as edible pigment. Currently, since crocetin manufacture mainly relied on extraction and purification from Crocus stigmas, deficient resource and low extraction rate restricted the large-scale application for commercialization. De novo synthesis of crocetin from simple carbon (glucose etc.) in engineered heterologous hosts would be an important complement to traditional sources. For crocetin biosynthesis, the conversion of β-carotene to crocetin required three steps catalyzed by β-carotene hydroxylase (CrtZ), carotenoid cleavage dioxygenase (CCD) and aldehyde dehydrogenase (ALD), respectively (Fig. 1a) . It is speculated that balancing metabolic flux mediated by the above three enzymes is a big challenge for high output. To date, only a patent from EVOVA Inc.  and a report from Lou et al.  have just illustrated the feasibility of heterologous biosynthesis of crocetin, and there was no specific titer data reported yet. For promoting crocetin production, combinatorial manipulation of the CrtZ, CCD and ALD would be a promising solution to overcome this challenge.
Screening enzymes sources and regulating protein expression level have been proved to be efficient strategies for manipulating the key enzymes for balancing metabolic flux, consequently improving production [11–13]. Cao et al.  once improved odd-chain fatty alcohols production in Escherichia coli through balancing the expression level of TesA, αDOX, AHRs and the genes involved in fatty acids metabolism pathway. Meanwhile, through combinatorially screening the carotenogenic enzymes (CrtE, CrtB and CrtI) from diverse organisms and fine-tuning the expression level of CrtI, an optimal enzymes combination with the highest lycopene yield was obtained in Saccharomyces cerevisiae . In crocetin biosynthesis fields, CrtZ, CCD and ALD have been characterized separately in the last decades. Li et al.  once achieved zeaxanthin titer as 43.46 mg/L in a recombinant E. coli strain by integrating Pantoea ananatis CrtZ into a β-carotene producing strain. Meanwhile, Crocus ZCD was firstly annotated as 7, 8 (7′, 8′)-zeaxanthin cleavage dioxygenase in 2003 . However, Frusciante et al.  demonstrated this enzyme could not achieve crocetin synthesis in E. coli and corn. Another two Crocus CCDs, CCD2  and ZCD1 , have been proved to cleavage of zeaxanthin at the 7, 8- and 7′, 8′-positions for forming crocetin dialdehyde in E. coli and Chlorella vulgaris, respectively. Moreover, even though EVOVA Inc.  and Lou et al.  realized crocetin synthesis by using endogenous ALD in yeast and algae, respectively, there was no titer data uncovered yet. It could guess that the complexity of fine-tuning CrtZ, CCD and ALD was the main obstacle. Therefore, it is urgent to explore CrtZ, CCD and ALD systematically for crocetin higher production.
Saccharomyces cerevisiae has been reported as a safe (Generally Recognized as Safe, GRAS) and robust host cell to produce heterologous carotenoids, including lycopene , β-carotene  and astaxanthin . Thus, in our study, crocetin was successfully synthesized in S. cerevisiae through incorporating heterologous CrtZ and CCD in an existing β-carotene producing strain SyBE_Sc0014CY06 (with β-carotene titer of 220 mg/L) (Table 1). A higher crocetin titer was achieved by adjusting the culture temperature from 30 to 20 °C. The production of crocetin was further enhanced by 2.8-fold via screening of CrtZ/CCD combination and ALD sources. Moreover, the crocetin titer was reached to 1219 μg/L by increasing the copy numbers of ccd and ald. Finally, the highest reported crocetin titer as 6278 μg/L was archived in 5-L bioreactors. This study sets a good example of fine-tuning multiple enzymes systematically for heterologous biosynthesis of desired pharmaceuticals and chemicals.
Construction of plasmids and strains
Primers and plasmids used in this study were listed in Additional file 1: Table S1; Table 1, respectively. All the heterologous genes including crtZ, ccd, and ald were codon optimized (Additional file 1: Table S2) and synthesized by GENEWIZ (Suzhou, China). All these genes were delivered as pUC57-simple serious plasmids (Table 1). Promoters (P GAL1 , P GAL7 and P GAL10 ), terminators (T HIS5 , T TEF2 , and T PGI1 ) and integration homologous arms (ho_L and ho_R) were amplified from the genomic DNA of S. cerevisiae CEN.PK2-1C, as well as the auxotroph marker URA3 was amplified from the plasmid pRS416. Cassette ho_L-P GAL1 -T HIS5 -URA3-ho_R was assembled by overlap extension PCR (OE-PCR) and cloned into pJET1.2, obtaining the plasmid pJET1.2-Z-01 (Table 1; Additional file 1: Figure S1). Genes crtZ were recovered by BsaI digestion from pUC57-Simple-01–09 and inserted into the same site of pJET1.2-Z-01, generating pJET1.2-Z series plasmids (CrtZ expression cassette plasmids pJET1.2-Z-02–10, Table 1; Additional file 1: Figure S1). Then the CrtZ expression cassette ho_L-P GAL1 -CrtZ-T HIS5 -URA3-ho_R were cut from pJET1.2-Z series plasmids by PmeI and transformed into S. cerevisiae SyBE_SC0014CY06 for genomic integration (Fig. 1b) via the lithium acetate method . Marker URA3 was deleted according to Boeke et al. , obtaining zeaxanthin producing strains SyBE_Sc0123Cz10-18 (Table 1) as the host cell in our study.
For constructing the initial crocetin producing strain and screening CrtZ/CCD combination, only heterologous CCDs were carried by single copy plasmid pRS416 and introduced into zeaxanthin producing strains (Fig. 1b). Genes ccd were amplified from the plasmid pUC57-Simple-10–13 and assembled together with promoter P GAL10 , terminators T HIS5 and T TEF2 into CCD expression cassette T HIS5 -P GAL10 -CCD-T TEF2 by OE-PCR. The products were inserted into the NotI site of plasmid pRS416, obtaining pRS416-C serious plasmids (CCD expression plasmids pRS416-C-01-04, Table 1; Additional file 1: Figure S2). These plasmids were transferred into zeaxanthin producing strains according to Table 1, producing crocetin producing strains (Table 1).
For screening ALD sources, heterologous CCD and ALD were carried by centromeric plasmid pRS416 and introduced into zeaxanthin producing strain (Fig. 1c). Cassette T TEF2 -P GAL7 -T PGI1 was also assembled by OE-PCR and cloned into pRS425 K, obtaining the plasmid pRS425 K-A-01 at first (Table 1; Additional file 1: Figure S3). Genes ald were recovered by BsaI digestion from pUC57-Simple-14–16 and inserted into the same site of pRS425 K-A-01, generating pRS425 K-A series plasmids (pRS425 K-A-02–04, Table 1; Additional file 1: Figure S3). Meanwhile, cassette T TEF2 -T PGI1 was assembled by OE-PCR. The product was incubated with XhoI/SacI and inserted into the same sites of pRS416, producing pRS416-A-01. Then cassettes T HIS5 -P GAL10 -ccd2-T TEF2 (digested from pRS416-C-01 by NotI), T TEF2 -P GAL7 -ald-T PGI1 (digested from pRS425 K-A-02–04 by PstI/BamHI) and linearized vector pRS416-A-01 (digested by BamHI) were assembled based on RADOM method in the particular zeaxanthin producing strain (producing strains SyBE_Sc0123C048–50 harboring plasmids pRS416-A-02–04 respectively, Table 1; Additional file 1: Figure S3) . For adjusting the expression level of CCD and ALD, multiple plasmid pRS426, instead of pRS416, was employed to carry CCD and ALD expression cassettes. Similar procedures were taken as motioned above, which were presented in Additional file 1: Figure S3.
Strains and culture conditions
Escherichia coli DH5α or TransT1 was used for plasmid construction, which was cultured at 37 °C in Luria–Bertani medium  supplemented with 50 μg/mL kanamycin or 100 μg/mL ampicillin for selection. Meanwhile, all the engineered yeast strains summarized in Table 1 were based on an existing β-carotene producing strain, S. cerevisiae SyBE_SC0014CY06. Engineered yeast strains were cultured on YPD medium or synthetic complete (SC) medium lacking appropriate nutrient component for selection . When needed, 1% (w/v) d-(+)-galactose were used as the inducer in fermentations and supplied into YPD medium (generating YPDG medium).
For shake-flask cultivation, colonies on solid plates were picked up and cultured in 3 mL SC medium for overnight growth at 30 °C. Then the preculture was transferred into 25 mL fresh SC medium and grew until reaching to mid-log phase. After that, the seed culture was inoculated into 50 mL YPD medium with an initial OD600 of 0.1 and cultivated at 30 °C for 72 h or 20 °C for 96 h. All the fermentation experiments were performed in triplicate.
The strain SyBE_Sc0123C053 was used for fed-batch fermentation. 100 µL glycerol-stock was inoculated into 25 mL SC medium and cultured at 30 °C, 250 rpm for overnight growth. Then the preculture was transferred to 200 mL fresh SC medium and grew until entering mid-exponential phase. Seed cultures were transferred to 1.8 L YPD medium (20 g/L glucose) in a 5 L bioreactor (BLBIO-5GJG-2, Shanghai, China) at a 10% (v/v) inoculum. The pH was automatically controlled at 5.5 with ammonia hydroxide (6 M). And the dissolved oxygen was kept at 40% by agitation cascade from 400 to 600 rpm, while the air flow was set at 2.5 vvm.
As the crocetin production modules were controlled by employed galactose-inducible system, the fed-batch fermentation should be divided into two stages: cell growth stage and crocetin accumulation stage. During the period of the cell growth stage, fermentation was carried out at 30 °C. The glucose concentration was monitored every 2 h and the glucose consumption rate was obtained accordingly. Based on this data, the glucose concentration was maintained less than 1 g/L by adding an appropriate volume of concentrated glucose solution (500 g/L) continuously. And 5 g yeast extract was added into the bioreactor every 12 h by feeding 400 g/L yeast extract stock solution. When the cell growth fell into stable phase, fermentation entered the second stage: crocetin accumulation stage. Then after fermentation temperature reduced to 20 °C, 10 g/L of d-(+)-galactose was fed to induce crocetin biosynthesis. As glucose was exhausted, cells begun to use ethanol as carbon source. The ethanol concentration was controlled below 5 g/L through adjusting the feeding rate of ethanol until harvest. Duplicate samples were collected to determine the cell density, glucose concentration, ethanol concentration and crocetin production. To avoid the spontaneous degradation from light, bioreactor should be covered with foils.
Extraction and analysis of carotenoids
To determine carotenoids accumulation, standards of lycopene, β-carotene and zeaxanthin were purchased from Sigma (Sigma-Aldrich, MO, USA), and standard of crocetin was purchased from Yuanye Bio-Technology (Shanghai, China). The procedures for extracting and analyzing carotenoids were modified according to Xie et al. . To be specific, after harvested cells were washed with distilled water, the cell pellet was re-suspended in 3 N HCl and boiled for 2 min, and then immediately cooled in ice for 3 min. Then cells debris were harvested and resuspended in acetone containing 1% (w/v) butylated hydroxytoluene. The above mixture was vortexed until colorless. After centrifugation, the acetone phase containing the extracted carotenoid was collected and evaporated by nitrogen blow. The products were analyzed by high-performance liquid chromatography system (HPLC, Waterse2695, Waters Corp, USA) equipped with a BDS HYPERSIL C18 column (150 mm × 4.6 mm, 5 μm, Thermo Scientific) and a UV/VIS detector (Waters 2489). To characterize lycopene, β-carotene and zeaxanthin, the product was dissolved in acetone and the signals were detected at 450 nm. The mobile phase consisting of acetonitrile-methanol (65:35 v/v) was chosen with a flow rate of 0.8 mL/min and the column temperature was set at 25 °C. In the meanwhile, for crocetin analysis, sample was dissolved in methanol-dimethylformamide (7:1 v/v) and crocetin was detected at 430 nm. 70% (v/v) methanol–water (containing 2% formic acid) was utilized as the mobile phase with a flow rate of 1 mL/min at 40 °C. Notably, considering that carotenoids are extremely unstable and susceptible to light, brown centrifugal tubes were used in the above procedures to avoid exposure to light.
Bioinformatics and structural analysis of CCD
The protein identified sequences of the target CCD from different taxa were queried from protein knowledgebase (UniProtKB) available at http://www.uniprot.org/, using the key term “carotenoid cleavage dioxygenase”, and subjected to a brief bioinformatics analysis to guarantee suitable diversity. Initially the CCD protein sequences were aligned by means of clustal W with default settings . Phylogenetic tree of CCD gene family was conducted in MEGA7  and inferred by Neighbor-Joining method . The bootstrap consensus tree deduced from 1000 replicates was taken to represent the evolutionary history of the taxa analyzed .
The structures of the CCD2 and CCD3 were both modeled based on the target-template (PDB ID: 2biw) alignment using SWISS-MODEL [30, 31]. And the Coordinates which are conserved between the targets and the template are copied from the template to the model. Insertions and deletions are remodeled using a fragment library. Side chains are then rebuilt. Finally, the geometry of the resulting model is regularized by using a force field. The modeled structures of target proteins were resolved with PyMol software .
Results and discussion
Construction of inducible crocetin biosynthesis pathway
To realize crocetin biosynthesis, heterologous crtZ and ccd were codon optimized and introduced into an existing β-carotene producer (S. cerevisiae SyBE_SC0014CY06), which processed endogenous ALDs to catalyze the final step in crocetin synthesis pathway (Fig. 1a) . At first, crtZ was integrated into the ho locus of the chromosome, while ccd was carried by centromeric plasmid pRS416. The expression of CrtZ and CCD were under the control of galactose-regulated GAL promoters GAL1 and GAL10, respectively (Fig. 1b). Because a highest zeaxanthin production was once achieved in yeast strain harboring CrtZ from Erwinia uredovora (Eu_CrtZ) among nine selected CrtZ species , Eu_CrtZ were also selected and intergraded into the chromosome of strain SyBE_SC0014CY06, generating strain SyBE_Sc0123Cz12 as a host cell in our study. In the meanwhile, CCD2 from Crocus was also selected to convert zeaxanthin to crocetin dialdehyde, obtaining strain SyBE_Sc0123C009. Strains SyBE_Sc0123C009 and SyBE_Sc0123Cz12 together with the parent strain SyBE_SC0014CY06 were cultured in shake-flask with YPDG medium at 30 °C and their products were analyzed by HPLC after 72 h incubation. As shown in Fig. 1d, crocetin (peak III) was successfully detected with a titer as 9.42 μg/L in strain SyBE_Sc0123C009, indicating that a functional crocetin biosynthesis pathway succeeded here. To be notably, there was no distinct β-carotene accumulation in zeaxanthin producing strain SyBE_Sc0123Cz12, while an amount of β-carotene (peak II), zeaxanthin (peak I), as well as other unidentified byproducts or intermediates were observed in crocetin producing strain SyBE_Sc0123C009 (Fig. 1e), suggesting that the step catalyzed by CCD was rate-limiting here and the selected CrtZ/CCD combination did not match well, which needed to be optimized further.
Optimization of cultivation temperature
It is reported by Shi et al.  that low temperature was benefit for carotenoids accumulation in Phaffia rhodozyma. In our study, by cultivation of series zeaxanthin producing strains at 20 and 30 °C respectively, it was also found that the production of zeaxanthin was higher at 20 °C than that at 30 °C (Additional file 1: Figure S4), indicating lower temperature benefited much for zeaxanthin production, which would provide more sufficient precursor supplies for higher crocetin production. Moreover, concerning that root development and flower emergence occur at low temperature for Crocus plants, and the expression of CCD were induced by low temperature in Crocus [35–37], the effect of culture temperature was also investigated here. Thus, for higher crocetin titer, the culture temperature for strain SyBE_Sc0123C009 was decreased from 30 °C, via 25 to 20 °C. The cell density, zeaxanthin accumulation and crocetin production were measured during the time course. As a result, in case of cell growth, there was a longer lag phase under lower temperature, compared to cultivating at 30 °C (Fig. 2a). Meanwhile, a dramatical increase on crocetin production along with a decrease on zeaxanthin accumulation was achieved by reducing cultivation temperature (Fig. 2b, c), suggesting 20 °C was the optimal temperature for converting zeaxanthin to crocetin. Javiera López et al.  once reported that β-ionone producing yeast strain processing CCD1, the homologue of CCD2, worked much better at low temperature, which showed similar results as our study. Finally, the crocetin titer reached 223 μg/L at 20 °C after 96 h fermentation in shake-flask (Fig. 2b). And 20 °C was used as the culture temperature in further study.
Optimal CrtZ/CCD combination by screening enzymes from diverse sources
As mentioned above, combinatorially screening enzymes from diverse sources has been proved to be a promising method to obtain the best combination in terms of substrate selectivity, catalytic activity and host cell compatibility, which would lead to higher productivity of the target compound [39–42]. Through blastp searching through NCBI database (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastp&PAGE_TYPE=Blast) Search&LINK_LOC = blasthome based on the sequence of CCD2, CCD3 showed a 97% identity with CCD2 (Fig. 4a). Hence besides three crocetin synthesis related CCDs (ZCD, ZCD1 and CCD2) described before [10, 17, 18], CCD3 was selected as potential candidate in our study. Here, these four CCDs together with nine CrtZs from E. uredovora (Eu_CrtZ), Pantoea agglomerans (Pa_CrtZ), Sulfolobus solfataricus P2 (Ss_CrtZ), Pantoea stewartii (Ps_CrtZ), Brevundimonas sp. SD212 (B.SD_CrtZ), Brevundimonas vesicularis DC263 (B.DC_CrtZ), Haematococcus pluvialis (Hp_CrtZ), Agrobacterium aurantiacum (Aa_CrtZ), Alcaligenes sp. PC-1 (As_CrtZ) were expressed in strain SyBE_SC0014CY06, generating 36 strains with diverse CrtZ/CCD combinations (Fig. 3a; Table 1). Nine strains carrying different CrtZs without CCDs introduced were used as the blank control (Fig. 3b, c; Table 1). All the above strains were cultured in YPDG medium to analyze the accumulation of zeaxanthin and crocetin. As illustrated in Fig. 3b, only the strain harboring CCD2 instead of other three CCDs could achieve crocetin accumulation in yeast, furtherly demonstrating CCD was a rate-limiting enzyme in crocetin synthesis pathway. Rather than CrtZ, CCD seemed to be more crucial for crocetin production. Finally, the combination as Ps_CrtZ/CCD2 achieved the highest crocetin titer as 351 μg/L in strain SyBE_Sc0123C017. This optimal combination would be a promising candidate for further optimization.
In this study, ZCD, ZCD1 and CCD3 could not achieve crocetin production in yeast, which required sequential cleavage at C7–C8 and C7′–C8′ double bonds adjacent to the 3-OH-β-ionone ring . Even though there was no crocetin detected in strains carrying these three enzymes separately, zeaxanthin accumulations were consumed at varying degrees in these strains, suggesting their cleave activities in yeast might at other position or only at one side of the molecules. Among the five subfamilies of plant CCDs, the CCD1 and CCD4 families were the only two involved in the cleave activities at 7, 8/7′, 8′ positions [44, 45]. A phylogenetic analysis of CCD sequences from diversity sources belonging to CCD1 and CCD4 families, illustrated that CCD2 and CCD3 belonged to CCD1 subfamily, while ZCD and ZCD1 were members of CCD4 subfamily (Additional file 1: Figure S5).
For CCD2 and CCD3, they shared 97% identities and exhibited dramatically diversity on enzyme activities. Through alignment of their protein sequences, seven dissimilar short fragments were detected (Fig. 4a). In order to further characterize these differences, the structural models of CCD2 and CCD3 were generated based on the crystal structure of their homological protein apocarotenoid cleavage oxygenase from Synechocystis (PDB accession ID: 2biw). As shown in Fig. 4b, CCD comprised seven bladed β-propellers, which is highly conserved among all CCDs and covered by a less rigid dome formed by a series of loops . To be notably, there is a tunnel perpendicular to the propeller axis of CCD. As reported, the tunnel acted as a channel for the passage of their hydrophobic substrates to the active site, and was consisted of hydrophobic residues (mainly Phe, Val, Leu) interacting with their lipophilic substrates via hydrophobic forces to guarantee both the specificity and correct orientation of substrate for the cleavage reactions . Thus, when the high hydrophobicity of the tunnel was subsided by the alteration in Fragment 5 (which located at the tunnel) as the residues of K320-F321 from CCD2 and E321-I322 from CCD3 (Fig. 4b), the substrate entrance to CCD3 was impacted for change on substrate specificity consequently. Meanwhile, the entrance of the tunnel located in a large hydrophobic patch for membrane insertion, which provided an appropriate environment for lipophilic substrates accommodation and enzyme contraction. The function of this hydrophobic patch mainly depended on the stable α-helices region, which was involved in the Fragment 3 and Fragment 4. The structure of CCD3 in Fragment 3 and Fragment 4 showed the longer and more unstable loops than that of CCD2 (Fig. 4b). As illustrated in Fig. 3c, besides substrates selectivity, CCD3 exhibited lower cleavage activity on zeaxanthin than CCD2, no matter cooperated with what kind of CrtZ sources. These results could be explained by above descriptions. Moreover, there were still some variances between CCD3 and CCD2 which could not support above results by current protein model. Therefore, a more delicate phylogenetic analysis of CCD sequences only from CCD1 family members were performed and showed that those unexplained different residues were highly conserved among all the tested CCD1 subfamily members except CCD3 (Additional file 1: Figure S7), suggesting the alternation on these conserved regions which might be essential to CCD function would reduce enzyme activities.
For ZCD and ZCD1, they share 96% identities (Additional file 1: Figure S6), and both truncated at the N-terminal as lacking a blade of β-propeller and part of the dome in classic CCD4 subfamily members. The truncation was once proved to lead to loss on any cleavage activity for ZCD in E. coli . ZCD1 was reported to once achieve crocetin production in C. vulgaris . However, in our study, both these two enzymes could not sequentially cleave zeaxanthin on 7, 8/(7′, 8′) positions in yeast (Fig. 3b). These conflicting dates highlight the importance of host cell compatibility on the performance of heterologous enzymes, which were also corroborated by the reports from Greene et al. .
Screening ALD sources and fine-tuning of CCD/ALD
As so far, there is no ALD has been identified in Crocus for crocetin synthesis. Meanwhile, except endogenous ALDs in yeast (such as ALD6) and algae, none heterologous ALD has been reported yet to realize crocetin producing. Since the current crocetin titer, which was achieved by yeast endogenous ALDs, was still low, it is urgent to search and screen ALD isozymes from other organisms for higher crocetin production. Here, besides yeast endogenous ALD6 , two heterologous ALD originated from Bixa orellana (Bix_ALD)  and Synechocystis sp. PCC6803 (Syn_ALD) , whose substrates share the similar structure with crocetin dialdehyde, were selected and introduced together with CCD2 into the strain with Ps_CrtZ integrated in its chromosome (Fig. 5a). CCD2 and ALD were carried by single copy plasmid pRS416 and placed under the control of promoters GAL10 and GAL7, respectively (Fig. 1c). After growing in YPDG medium for 96 h, strain SyBE_Sc0123C050 harboring Syn_ALD achieved higher crocetin titer as 633 μg/L (Fig. 5b). Moreover, by increasing copy numbers of CCD2 and Syn_ALD via interchange of vector pRS416 into multicopy plasmid pRS426, the crocetin titer further improved to 1219 μg/L (Fig. 5b), obtaining strain SyBE_Sc0123C053 for bioreactor experiment.
Optimization of crocetin production in bioreactor
To evaluate the production performance of the engineered strain SyBE_Sc0123C053, fed-batch fermentation was performed at a 2 L scale using YPD as the medium (Fig. 6). During cell growth stage, based on carbon restriction strategy, glucose concentration was strictly restricted. Cell density reached 96 for 35 h cultivation at 30 °C. There was also no acetate observed in this stage (data was not shown). When the culture temperature reduced to 20 °C at 36 h, d-(+)-galactose was added to induce crocetin production. After the initial ethanol generated by glucose was consumed below 5 g/L, additional 100 mL ethanol was fed into the medium periodically to maintain ethanol concentration around at 5 g/L until harvest. Eventually, a crocetin titer of 6278 μg/L was obtained after 124 h cultivation (Fig. 6), which was the highest reported titer in eukaryotic cell to date. However, the absolute titer (6278 μg/L) and the production yield based on ethanol consumption (YP/S = 0.012%) were far away from commercialization, strain engineering by metabolic engineering as well as synthetic biology and process innovation would be two basic but efficient aspects to promote crocetin output. In terms of strain engineering, increasing the catalysis activity of CCD and strain tolerance to product were the main challenges. Combinatorial engineering of S. cerevisiae and crocetin biosynthesis pathway in parallel would probably meet the demand . As recent study in process optimization has demonstrated great potential in isoprene overproduction (up to 24 g/L) , we believe that crocetin production by our engineered strain would be further improved by continuous efforts in metabolic engineering, synthetic biology and fermentation optimization.
In our study, crocetin biosynthesis pathway was successfully established in S. cerevisiae through incorporating heterologous CrtZ and CCD in an existing β-carotene producing strain. Then the effects of culture temperature, combination of CrtZ/CCD, ALD from different species, as well as the expression level of CCD and ALD on crocetin were investigated respectively. Compared to culturing at 30 °C, the crocetin accumulation performed much better at 20 °C. The accumulation of crocetin was further promoted by 2.8-fold by screening of CrtZ/CCD combination and ALD sources. Moreover, the crocetin titer was reached to 1219 μg/L by overexpression of ccd and ald. Consequently, the highest reported crocetin titer of 6278 μg/L was obtained in 5-L bioreactors. This study promotes the opportunities for industrialization of crocetin and crocin. This study also sets a good reference for microbial production of pharmaceuticals and chemicals in complex structure by fine-tuning multiple enzymes systematically.
carotenoid cleavage dioxygenase
- Crocus :
Crocus sativus L.
Sulfolobus solfataricus P2
Brevundimonas sp. SD212
Brevundimonas vesicularis DC263
Alcaligenes sp. PC-1
Synechocystis sp. PCC6803
Winterhalter P. Carotenoid-derived aroma compounds: an introduction. 2001.
Dhar A, Mehta S, Dhar G, et al. Crocetin inhibits pancreatic cancer cell proliferation and tumor progression in a xenograft mouse model. Mol Cancer Ther. 2009;8(2):315–23.
Bathaie SZ, Hoshyar R, Miri H, et al. Anticancer effects of crocetin in both human adenocarcinoma gastric cancer cells and rat model of gastric cancer. Biochem Cell Biol. 2013;91(6):397–403.
Hsu JD, Chou FP, Lee MJ, et al. Suppression of the TPA-induced expression of nuclear-protooncogenes in mouse epidermis by crocetin via antioxidant activity. Anticancer Res. 1999;19(5B):4221–7.
Higashino S, Sasaki Y, Giddings JC, et al. Crocetin, a carotenoid from Gardenia jasminoides Ellis, protects against hypertension and cerebral thrombogenesis in stroke-prone spontaneously hypertensive rats. Phytother Res. 2014;28(9):1315–9.
Zheng S, Qian Z, Tang F, et al. Suppression of vascular cell adhesion molecule-1 expression by crocetin contributes to attenuation of atherosclerosis in hypercholesterolemic rabbits. Biochem Pharmacol. 2005;70(8):1192–9.
Amin B, Nakhsaz A, Hosseinzadeh H. Evaluation of the antidepressant-like effects of acute and sub-acute administration of crocin and crocetin in mice. Avicenna J Phytomedicine. 1900;5(5):458–68.
Pfander H, Schurtenberger H. Biosynthesis of C20-carotenoids in Crocus sativus. Phytochemistry. 1982;21(5):1039–42.
Raghavan S, Hansen J, Sonkar S, et al. Methods and materials for recombinant production of saffron compounds: WO, WO/2013/021261[P]. 2013.
Lou S, Wang L, He L, et al. Production of crocetin in transgenic Chlorella vulgaris expressing genes crtRB and ZCD1. J Appl Phycol. 2016;28(3):1657–65.
Ajikumar PK, Xiao WH, Tyo KE, et al. Isoprenoid pathway optimization for taxol precursor overproduction in Escherichia coli. Science. 2010;330(6000):70–4.
Song MC, Kim EJ, Kim E, et al. Microbial biosynthesis of medicinally important plant secondary metabolites. Nat Product Rep. 2014;31(11):1497–509.
Yao YF, Wang CS, Qiao J, et al. Metabolic engineering of Escherichia coli for production of salvianic acid A via an artificial biosynthetic pathway. Metab Eng. 2013;19(5):79–87.
Cao YX, Xiao WH, Liu D, et al. Biosynthesis of odd-chain fatty alcohols in Escherichia coli. Metab Eng. 2015;29:113–23.
Chen Y, Xiao W, Wang Y, et al. Lycopene overproduction in Saccharomyces cerevisiae through combining pathway engineering with host engineering. Microb Cell Fact. 2016;15(1):1–13.
Li XR, Tian GQ, Shen HJ, et al. Metabolic engineering of Escherichia coli to produce zeaxanthin. J Ind Microbiol Biotechnol. 2015;42(4):627–36.
Bouvier F, Suire C, Mutterer J, et al. Oxidative remodeling of chromoplast carotenoids identification of the carotenoid dioxygenase CsCCD and CsZCD genes involved in crocus secondary metabolite biogenesis. Plant Cell. 2003;15(1):47–62.
Frusciante S, Diretto G, Bruno M, et al. Novel carotenoid cleavage dioxygenase catalyzes the first dedicated step in saffron crocin biosynthesis. Proc Natl Acad Sci USA. 2014;111(33):12246–51.
Xie W, Lv X, Ye L, et al. Construction of lycopene-overproducing Saccharomyces cerevisiae by combining directed evolution and metabolic engineering. Metab Eng. 2015;30:69–78.
Xie W, Liu M, Lv X, et al. Construction of a controllable β-carotene biosynthetic pathway by decentralized assembly strategy in Saccharomyces cerevisiae. Biotechnol Bioeng. 2014;111(1):125–33.
Zhou P, Ye L, Xie W, et al. Highly efficient biosynthesis of astaxanthin in Saccharomyces cerevisiae by integration and tuning of algal crtZ and bkt. Appl Microbiol Biotechnol. 2015;99(20):8419–28.
Gietz RD, Schiestl RH. High-efficiency yeast transformation using the LiAc/SS carrier DNA/PEG method. Nat Protoc. 2007;2:31–4.
Boeke JD, Lacroute F, Fink GR. A positive selection for mutants lacking orotidine-5′-phosphate decarboxylase activity in yeast: 5-fluoro-orotic acid resistance. Mol Genet Genom. 1984;197(2):345–6.
Lin Q, Jia B, Mitchell LA, et al. RADOM, an efficient in vivo method for assembling designed DNA fragments up to 10 kb long in Saccharomyces cerevisiae. Acs Synth Biol. 2015;4(3):213–20.
Gietz RD, Woods RA. Yeast transformation by the LiAc/SS carrier DNA/PEG method. Methods Mol Biol. 2014;313:107–20.
Thompson JD, Higgins DG, Gibson TJ. CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, positions-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994;22:4673–80.
Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Sanderson MJ, Wojciechowski MF. Improved bootstrap confidence limits in large-scale phylogenies, with an example from Neo-Astragalus (Leguminosae). Syst Biol. 2000;49(4):671–85.
Biasini M, Bienert S, Waterhouse A, et al. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014;42(w1):252–8.
Arnold K, Bordoli L, Kopp J, et al. The SWISS-MODEL workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006;22(2):195–201.
Delano WL. The PyMOL molecular graphics system. My publications, 2010.
Mei XA, Chen Y, Wang R, et al. Construction of zeaxanthin pathway in Saccharomyces cerevisiae. China Biotechnol. 2016;36(8):64–72.
Shi F, Zhan W, Li Y, et al. Temperature influences β-carotene production in recombinant Saccharomyces cerevisiae expressing carotenogenic genes from Phaffia rhodozyma. World J Microbiol Biotechnol. 2014;30(1):125–33.
Ahrazem O, Rubiomoraga A, Nebauer SG, et al. Saffron: its phytochemistry, developmental processes, and biotechnological prospects. J Agric Food Chem. 2015;63(40):8751–64.
Molina RV, Valero M, Navarro Y, et al. Low temperature storage of corms extends the flowering season of saffron (Crocus sativus L.). J Hortic Sci Biotechnol. 2005;80(3):319–26.
Molina RV, Valero M, Navarro Y, et al. Temperature effects on flower formation in saffron (Crocus sativus L.). Sci Hortic. 2005;103(3):361–79.
López J, Essus K, Kim IK, et al. Production of β-ionone by combined expression of carotenogenic and plant CCD1 genes in Saccharomyces cerevisiae. Microb Cell Fact. 2015;14(1):1–13.
Kim E, Moore BS, Yoon YJ. Reinvigorating natural product combinatorial biosynthesis with synthetic biology. Nat Chem Biol. 2015;11(9):649–59.
Ding MZ, Yan HF, Li LF, et al. Biosynthesis of taxadiene in Saccharomyces cerevisiae: selection of geranylgeranyl diphosphate synthase directed by a computer-aided docking strategy. PLoS ONE. 2014;9(10):e109348.
Chang JJ, Thia C, Lin HY, et al. Integrating an algal β-carotene hydroxylase gene into a designed carotenoid-biosynthesis pathway increases carotenoid production in yeast. Bioresour Technol. 2015;184:2–8.
Sarria S, Wong B, Martín HG, et al. Microbial synthesis of pinene. Acs Synth Biol. 2014;3(7):466–75.
Moraga AR, Nohales PF, Pérez JAF, et al. Glucosylation of the saffron apocarotenoid crocetin by a glucosyltransferase isolated from Crocus sativus stigmas. Planta. 2004;219(6):955–66.
Ahrazem O, Rubio-Moraga A, Berman J, et al. The carotenoid cleavage dioxygenase CCD2 catalysing the synthesis of crocetin in spring crocuses and saffron is a plastidial enzyme. New Phytol. 2016;209(2):650–63.
Priya R, Siva R. Phylogenetic analysis and evolutionary studies of plant carotenoid cleavage dioxygenase gene. Gene. 2014;548(2):223–33.
Oussama Ahrazem, Gómez-Gómez Lourdes, Rodrigo María J, et al. Carotenoid cleavage oxygenases from microbes and photosynthetic organisms: features and functions. Int J Mol Sci. 2016;17:1781.
Sui X, Kiser PD, Lintig JV, et al. Structural basis of carotenoid cleavage: from bacteria to mammals. Arch Biochem Biophys. 2011;539(2):203–13.
Greene JJ. Host cell compatibility in protein expression. Methods Mol Biol. 2004;267:3–14.
Saint-Prix F, Bönquist L, Dequin S. Functional analysis of the ALD gene family of Saccharomyces cerevisiae during anaerobic growth on glucose: the NADP+-dependent Ald6p and Ald5p isoforms play a major role in acetate formation. Microbiology. 2004;150(7):2209–20.
Bouvier F, Dogbo O, Camara B. Biosynthesis of the food and cosmetic plant pigment bixin (annatto). Science. 2003;300(5628):2089–91.
Trautmann D, Beyer P, Al-Babili S. The ORF slr0091 of Synechocystis sp. PCC6803 encodes a high-light induced aldehyde dehydrogenase converting apocarotenals and alkanals. FEBS J. 2013;280(15):3685–96.
Chen Y, Xiang G, Yu J, et al. Synergy between methylerythritol phosphate pathway and mevalonate pathway for isoprene production in Escherichia coli. Metab Eng. 2016;37:79–91.
FC, WX and YY conceived of the study and drafted the manuscript. FC and XM carried out the molecular genetic studies. FC and WX carried out the fed-batch fermentation experiments. YW participated in design and coordination of the study and helped to draft the manuscript. MY carried out the protein analysis. YC and HL participated in strain construction and HPLC analysis respectively. WX supervised the whole research and revised the manuscript. All authors read and approved the final manuscript.
The authors are grateful for the financial support from the International S&T Cooperation Program of China (2015DFA00960), the National Natural Science Foundation of China (31600052 and 21676192) and Innovative Talents and Platform Program of Tianjin (16PTSYJC00050).
Availability of data and materials
The material and data supporting their findings can be found in the main paper and the additional file.
The authors declare that they have no competing interests.
The International S&T Cooperation Program of China (2015DFA00960), the National Natural Science Foundation of China (31600052 and 21676192) and Innovative Talents and Platform Program of Tianjin (16PTSYJC00050) supported this work.
Fenghua Chai, Ying Wang and Xueang Mei contributed equally to this work