The photosynthetic bacteria Rhodobacter capsulatus and Synechocystis sp. PCC 6803 as new hosts for cyclic plant triterpene biosynthesis

Abstract

Cyclic triterpenes constitute one of the most diverse groups of plant natural products. Besides the intriguing biochemistry of their biosynthetic pathways, plant triterpenes exhibit versatile bioactivities, including antimicrobial effects against plant and human pathogens. While prokaryotes have been extensively used for the heterologous production of other classes of terpenes, the synthesis of cyclic triterpenes, which inherently includes the two-step catalytic formation of the universal linear precursor 2,3-oxidosqualene, is still a major challenge. We thus explored the suitability of the metabolically versatile photosynthetic α-proteobacterium Rhodobacter capsulatus SB1003 and cyanobacterium Synechocystis sp. PCC 6803 as alternative hosts for biosynthesis of cyclic plant triterpenes. Therefore, 2,3-oxidosqualene production was implemented and subsequently combined with different cyclization reactions catalyzed by the representative oxidosqualene cyclases CAS1 (cycloartenol synthase), LUP1 (lupeol synthase), THAS1 (thalianol synthase) and MRN1 (marneral synthase) derived from model plant Arabidopsis thaliana. While successful accumulation of 2,3-oxidosqualene could be detected by LC-MS analysis in both hosts, cyclase expression resulted in differential production profiles. CAS1 catalyzed conversion to only cycloartenol, but expression of LUP1 yielded lupeol and a triterpenoid matching an oxidation product of lupeol, in both hosts. In contrast, THAS1 expression did not lead to cyclic product formation in either host, whereas MRN1-dependent production of marnerol and hydroxymarnerol was observed in Synechocystis but not in R. capsulatus. Our findings thus indicate that 2,3-oxidosqualene cyclization in heterologous phototrophic bacteria is basically feasible but efficient conversion depends on both the respective cyclase enzyme and individual host properties. Therefore, photosynthetic α-proteo- and cyanobacteria are promising alternative candidates for providing new bacterial access to the broad class of triterpenes for biotechnological applications.

About the Authors:

Anita Loeschcke

Contributed equally to this work with: Anita Loeschcke, Dennis Dienst, Vera Wewer

Roles Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing

Affiliations Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, Forschungszentrum Jülich, Jülich, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Dennis Dienst

Contributed equally to this work with: Anita Loeschcke, Dennis Dienst, Vera Wewer

Roles Formal analysis, Methodology, Visualization, Writing – original draft, Writing – review & editing

Affiliations Institute for Synthetic Microbiology, Department of Biology, Heinrich Heine University Düsseldorf, Düsseldorf, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Vera Wewer

Contributed equally to this work with: Anita Loeschcke, Dennis Dienst, Vera Wewer

Roles Formal analysis, Investigation, Visualization, Writing – review & editing

Affiliations MS Platform, Department of Biology, University of Cologne, Cologne, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Jennifer Hage-Hülsmann

Roles Investigation, Methodology, Writing – review & editing

Affiliations Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, Forschungszentrum Jülich, Jülich, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Maximilian Dietsch

Roles Investigation, Writing – review & editing

Affiliations Institute for Synthetic Microbiology, Department of Biology, Heinrich Heine University Düsseldorf, Düsseldorf, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Sarah Kranz-Finger

Roles Methodology, Writing – review & editing

Affiliations Institute of Biochemistry II, Department of Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Vanessa Hüren

Roles Investigation, Writing – review & editing

Affiliations Institute for Synthetic Microbiology, Department of Biology, Heinrich Heine University Düsseldorf, Düsseldorf, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Sabine Metzger

Roles Formal analysis, Validation, Writing – review & editing

Affiliations MS Platform, Department of Biology, University of Cologne, Cologne, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Vlada B. Urlacher

Roles Conceptualization, Writing – review & editing

Affiliations Institute of Biochemistry II, Department of Chemistry, Heinrich Heine University Düsseldorf, Düsseldorf, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Tamara Gigolashvili

Roles Conceptualization, Writing – review & editing

Affiliations Botanical Institute, University of Cologne, Cologne, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Stanislav Kopriva

Roles Conceptualization, Supervision, Writing – review & editing

Affiliations Botanical Institute, University of Cologne, Cologne, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Ilka M. Axmann

Roles Conceptualization, Supervision, Validation, Writing – review & editing

* E-mail: [email protected] (IMA); [email protected] (TD)

Affiliations Institute for Synthetic Microbiology, Department of Biology, Heinrich Heine University Düsseldorf, Düsseldorf, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

Thomas Drepper

Roles Conceptualization, Methodology, Supervision, Validation, Writing – review & editing

* E-mail: [email protected] (IMA); [email protected] (TD)

Affiliations Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, Forschungszentrum Jülich, Jülich, Germany, Cluster of Excellence on Plant Sciences (CEPLAS)

ORCID http://orcid.org/0000-0002-0096-8084

Karl-Erich Jaeger

Roles Conceptualization, Funding acquisition, Supervision, Writing – review & editing

Affiliations Institute of Molecular Enzyme Technology, Heinrich Heine University Düsseldorf, Forschungszentrum Jülich, Jülich, Germany, Cluster of Excellence on Plant Sciences (CEPLAS), Institute of Bio- and Geosciences (IBG-1), Forschungszentrum Jülich, Jülich, Germany

Introduction

Plant secondary metabolites comprise a large variety of structurally divergent compounds that can serve as signaling molecules or as protecting agents against microbial pathogens, herbivorous attacks or plant competitors (overview in [1]). Among these secondary metabolites, terpenoids including the cyclic plant-type triterpenes constitute one of the largest and most diverse groups exhibiting important functions in plant physiology and development. For example, the linear triterpene precursor squalene plays a role in defense elicitation against herbivore attacks [2], the sterol cycloartenol is important for membrane functionality, plastid biogenesis and cell viability [3], β-amyrin is involved in root development [4], marneral affects general plant growth and development [5], and lupeol is essential for root nodule formation [6]. A group of glycosylated triterpenes, referred to as saponins, can act, for instance, as natural barriers against pathogenic microbial penetration due to their detergent-like characteristics (reviewed in [7]). Furthermore, triterpenes represent commercial potential in the pharmaceutical sector due to their bioactivities such as antiviral, anticancer, anti-inflammatory or wound-healing properties [8–13]. Antifertility effects of lupeol were recently demonstrated [14] and contraceptive application is discussed.

Like all terpenes, plant triterpenes are synthesized from C5 isoprene units. As a first step in the formation of triterpenes, the NAD(P)H-dependent enzyme squalene synthase (SQS) catalyzes the formation of the linear C30 isoprenoid squalene from C15 farnesyl pyrophosphate (FPP) by condensation of two FPPs to presqualene pyrophosphate and a subsequent reductive rearrangement to form squalene [15]. Subsequent squalene epoxidase (SQE)-catalyzed oxygenation of squalene yields 2,3-oxidosqualene, which is the linear key precursor of cyclic plant type triterpenes. SQE is a flavin-binding monooxygenase which requires NADPH and molecular oxygen to oxidize the substrate, and is reported to obtain electrons from cytochrome P450 reductase [16]. 2,3-oxidosqualene, in turn, can be further converted by diverse oxidosqualene cyclases (OSC) to a multitude of cyclic products. These plant triterpene products can be divided into sterol type molecules, which are generated via substrate folding into a transition state with chair-boat-chair (CBC) conformation, and the large group of triterpenes whose multiple skeletons are generated after substrate folding into a chair-chair-chair (CCC) conformation. A portion of the cyclic plant triterpene scaffolds is subsequently oxidized and further decorated thereby generating a larger structural diversity [17].

Consistent with the wide array of triterpene bioactivities, there is considerable effort in developing strategies for pathway identification and characterization of the corresponding enzymatic activities [18]. However, heterologous biosynthesis of cyclic plant-type triterpenes in microbial cells was thus far largely restricted to yeast systems, most notably Saccharomyces cerevisiae [17], and could only recently also be demonstrated for the triterpene dammarenediol-II in Escherichia coli [19]. To broaden the range of microbial systems for the functional heterologous biosynthesis of plant triterpenes, we chose to employ the physiologically versatile photosynthetic bacteria Rhodobacter capsulatus SB1003 and Synechocystis sp. PCC 6803 (hereafter: Synechocystis) as host strains. Both are well-established model organisms for anaerobic and aerobic bacterial photosynthesis, respectively. They have both gained attention as biotechnological production hosts for a variety of biochemicals [20, 21] and their plant growth promoting properties are of interest in the agricultural sector [22, 23]. Engineered strains synthesizing bioactive compounds in an agricultural context may contribute to novel farming strategies. Both bacteria represent promising platform hosts for triterpene biosynthesis due to unique phototrophy-associated physiological properties, i.e. a large intracellular membrane system, which can incorporate hydrophobic compounds and membrane-bound or -associated enzymes. In addition, both strains intrinsically harbor an effective isoprene metabolism for biosynthesis of carotenoids that serve as photopigments [24–26]. Isoprenoid pathway engineering has been demonstrated for different cyanobacterial and Rhodobacter strains, resulting in significant accumulation of hemi- (C5), mono- (C10), sesqui- (C15), di- (C20), tri- (C30) and tetraterpenes (C40) [27–35]. Recently, implementation of squalene biosynthesis has been demonstrated in both R. capsulatus [36] and Synechocystis [27, 37]. Notably, to enable C30 isoprenoid biosynthesis in R. capsulatus, heterologous expression of an SQS enzyme is essential [36]. In contrast, the natural hopanoid pathway of Synechocystis can be modified by disrupting the squalene hopene cyclase gene shc to accumulate squalene [37].

Here, we present for the first time the heterologous biosynthesis of the key plant triterpene precursor 2,3-oxidosqualene in the phototrophic prokaryotes R. capsulatus and Synechocystis by additional expression of SQE [38] in squalene producing cells. Further, the ability of both strains to accumulate cyclic plant triterpene compounds was evaluated by co-expression of selected OSC enzymes cycloartenol synthase CAS1 [39], lupeol synthase LUP1 [40], thalianol synthase THAS1 [41], and marneral synthase MRN1 [42], all derived from the model plant Arabidopsis thaliana.

Materials and methods

Bacterial strains and cultivation conditions

Escherichia coli strains DH5α [43], NEB5α (New England Biolabs), J53 [44] and S17-1 [45], used for cloning and conjugation, respectively, were cultivated on LB-agar plates or in liquid LB medium (Luria/Miller, Carl Roth) at 37 °C. Antibiotics were added to E. coli culture medium to the following final concentrations [μg mL-1]: 100 (ampicillin), 50 (kanamycin), 20 (spectinomycin). The photosynthetic purple non-sulfur α-proteobacterium R. capsulatus SB1003 [46] was cultivated on 2% (w/v) agar (Bacto agar; Difco) containing PY plates [47] or in RCV liquid medium [48] at 30 °C. The non-motile, glucose-tolerant (GT) strain of Synechocystis sp. PCC 6803 as well as the corresponding mutant strain Δshc (slr2089) used in this study were kindly provided by Pia Lindberg (Uppsala University, Sweden) and cultivated on 0.75% (w/v) agar (Bacto agar; Difco) plates containing BG-11 mineral medium or in BG-11 liquid medium [49] at 30 °C. Saccharomyces cerevisiae strain GIL77 (gal2 hem3-6 erg7 ura3-167) [50], carrying derivatives of vector pYES/DEST-52 (Invitrogen) with MRN1 or THAS1, or as empty vector control [51], were cultivated at 30 °C using synthetic complete medium without uracil, supplemented with glucose, ergosterol, hemin, and Tween, as previously described [52].

Construction of expression vectors

A. thaliana derived genes SQS1 (At4g34640), SQE1 (At1g58440), CAS1 (At2g07050), LUP1 (At1g78970), THAS1 (At5g48010), and MRN1 (At5g42600) were obtained as synthetic DNA by Eurofins Genomics with adapted codon-usage for expression in both bacterial hosts, and with RBS and restriction endonuclease recognition sequences where appropriate for initial cloning steps. An overview of constructs and genetic features, as well as primers used in this study is given in Table 1 and S1 Table, respectively. All constructs were verified by sequencing of relevant sequences (performed as commercial service by Eurofins Genomics, Seqlab-Microsynth or Sequiserve), as are specified in S2 Table. Plasmid maps are shown in S1 Fig.

[Figure omitted. See PDF.]

Table 1. Strain constructs used in this study and their genetic features.

https://doi.org/10.1371/journal.pone.0189816.t001

Vectors for expression in R. capsulatus.

Genes were cloned into expression vector pRhon5Hi-2, which is based on the expression vector pRhotHi-2 [53]. It harbors the native R. capsulatus promotor of the nitrogenase nif genes (Pnif), which was inserted as an NheI/XbaI-fragment (NCBI Genbank Accession MG208548; Özgür, Drepper et al, unpublished) into the respective sites of vector pRhotHi-2 [53]. Use of this promoter enables tight regulation via ammonia provision or limitation.

Squalene synthase gene SQS1, which was designed as a synthetic gene with an NdeI-site at the 5’-end as well as an XhoI- and HindIII-site at the 3’-end, was cloned into the vector as an NdeI/HindIII fragment, creating pRhon5Hi-2-SQS1. To generate vector pRhon5Hi-2-SQS1-SQE1, both genes were first cloned in pUC18 successively, using NdeI/HindIII for SQS1, and XhoI/HindIII for SQE1 with the pET vector RBS. The SQS1-SQE1 cassette was excised and cloned into the expression vector as NdeI/HindIII fragment. For construction of expression cassettes with the structure OSC-SQS1-SQE1, OSC genes were PCR-amplified for introduction into the expression vector pRhon5Hi-2-SQS1-SQE1 upstream of SQS1. For cloning of CAS1, LUP1 and MRN1, the synthetic DNA fragments were used as templates, and PCR primers (#510 - #515) were adapted to introduce NdeI sites on both ends to enable cloning at the first position of the operon, as well as the RBS of R. capsulatus nifK at the 3’ end of the fragment for the downstream gene SQS1. For cloning of THAS1, vector pRhon5Hi-2-H6-THAS1-SQS1-SQE1 was used as template (find sequence in S2 Table) for amplification of the synthetic gene together with the nifK RBS (using primers #516/ #517), which already carried the synthetic OSC gene with the following nifK RBS in front of the SQS1 gene. Therefore, each gene was equipped with its own RBS, as specified in Table 1 and S2 Table: The first operon gene was preceded by the RBS of the pRhotHi-2 vector, SQE1 was cloned with the pET vector RBS, and in OSC-containing constructs, SQS1 was equipped with the nifK RBS from R. capsulatus.

Vectors for expression in Synechocystis.

Expression plasmids for Synechocystis were all derivatives of the conjugative shuttle vector pVZ-spec [54, 55]. As basic cloning backbone, a pJET1.2 (Thermo Scientific) -based cloning plasmid pJPVCS (NCBI Genbank Accession MG191280; Schmelling, Dienst et al., unpublished) was used that harbored the native cobalt-inducible promoter PcoaT from Synechocystis including the upstream-located repressor gene coaR. The promoter-5’UTR sequence of coaT was fused to the mVenus CDS via an NheI site constituting the 2nd and 3rd codon (Ala-Ser). The whole expression cassette was flanked by the Biobrick prefix and suffix, which were further flanked by NsiI- and SalI sites, respectively, on this plasmid.

For construction of plasmids pVZ-PcoaT-SQE1-MRN1 and pVZ-PcoaT-SQE1-THAS1, the sequential cassettes SQE1-MRN1 and SQE1-THAS1 were PCR-amplified from plasmids pRhofHi2-SQS1-SQE1-MRN1 and pRhofHi2-SQS1-SQE1-THAS1 (Loeschcke, unpublished), carrying the synthetic plant genes with bacterial codon usage (as specified in S2 Table), respectively. The used primers were #252 (introducing an NheI site to the 2nd and 3rd codon of SQE1) and #253 (adding transcriptional terminator T7 106 bp downstream of the stop codon and an SpeI site to the 3’terminus of the amplicon). Both PCR products were digested with SpeI/ NheI and ligated into the correspondingly cut plasmid pJPVCS, thereby placing SQE1 downstream of the coaR-PcoaT-5’UTR-construct. Translation of the OSC genes was mediated by the Biobrick RBS BBa_0034. The expression cassette was excised by NsiI/ SalI digestion and transferred to the PstI/ SalI-cut pVZ-spec plasmid. Cloning of plasmid pVZ-PcoaT-SQE1 was based on a gBlock fragment (IDT) harboring the 3’-terminal 408 bp of SQE1 (starting with the CDS´s AgeI site at pos. 2794), the T7 terminator with a 29 bp spacer to the stop codon and a 3’-terminal SpeI site. The fragment was digested with AgeI and SpeI and ligated in the correspondingly cut vector pJET-PcoaT-SQE1-MRN1, thereby excising the MRN1 cassette including the contiguous T7 terminator. The coaR-PcoaT-5’UTR-SQE1 construct was excised from the pJET backbone by NsiI/ SalI digestion and transferred to the PstI/ SalI-cut pVZ-spec plasmid. For construction of plasmids pVZ-PcoaT-SQE1-CAS1 and pVZ-PcoaT-SQE1-LUP1, we used the AQUA (advanced quick assembly) technique [56] using plasmid pJET-PcoaT-SQE1-MRN1 as vector backbone. In order to exchange MRN1 with CAS1, AQUA cloning was performed with each three overlapping PCR fragments, generated with primer pairs #531/ #532, #528/ #529 and #526/ #527. For corresponding LUP1 cloning, primer pairs #533/ #532, #528/ #530 and #524/ #525 were used to amplify the AQUA fragments, before the assembly reaction was conducted in E. coli strain NEB5α. From the resulting plasmids, pJET-PcoaT-SQE1-CAS1 and pJET-PcoaT-SQE1-LUP1, inserts SQE1-CAS1 and SQE1-LUP1 were subsequently excised by NheI/ SalI digestion and transferred to NheI/ SalI digested plasmid pVZ-PcoaT-SQE1-MRN1.

Heterologous expression of plant triterpene biosynthesis genes in R. capsulatus

For introduction of A. thaliana triterpene biosynthesis genes in R. capsulatus, respective pRhon5Hi-2-based plasmids were transferred to the host via conjunctional transfer as previously described [47]. Thereafter, exconjugants were selected and further cultivated on PY agar, containing 25 μg mL-1 kanamycin and 25 μg mL-1 rifampicin.

For heterologous expression of triterpene biosynthesis genes, cultivation was conducted in liquid RCV medium with 4g/L malate containing 25 μg mL-1 kanamycin and 25 μg mL-1 rifampicin under non-phototrophic, micro-aerobic conditions (which is sufficient for the induction of carotenoid synthesis and formation of the intracellular membrane system in this organism) in 100 mL Erlenmeyer flasks at 30°C under permanent agitation at 130 rpm in a Multitron Standard incubation shaker (Infors HT) in the dark. Pre-cultures of 35 mL RCV medium containing 0.1% (NH4)2SO4 were inoculated and incubated for 48 h. Expression cultures were inoculated from pre-cultures to an initial OD660nm of 0.05 in 60 mL RCV medium containing 0.1% serine as sole nitrogen source. The omission of ammonium together with implementation of micro-aerobic conditions ensured derepression of the Pnif promoter and target protein expression. After incubation of cultures for 2 days (52 h), when an OD660nm of ~1.5 was reached, cell samples corresponding to OD660nm = 15 were harvested by centrifugation (10’, 1780 g, 4°C) and pellets stored at -20 °C.

Heterologous expression of plant triterpene biosynthesis genes in Synechocystis

For expression of plant triterpene biosynthesis genes in Synechocystis, pVZ-based expression plasmids were transferred to strain Δshc by conjugation [55], whereupon exconjugants were selected and further cultivated on BG11 agar containing, 20 μg mL-1 spectinomycin and 40 μg mL-1 kanamycin.

For triterpene biosynthesis experiments, cells were grown under continuous illumination with white light (50 μmol photons m–2 s–1) at 30 °C and 150 rpm in a New Brunswick Innova42 incubator shaker (Eppendorf). The culture volume was 80 mL within 250 mL Erlenmeyer flasks. Precultures were inoculated in BG-11 medium containing 20 μg mL-1 spectinomycin and 40 μg mL-1 kanamycin, and 10 mM TES buffer (pH 8.0) without cobalt supply, in order to repress gene expression from promoter PcoaT. For induction, cultures were supplemented with 20 μM CoCl2 at an optical density (OD750nm) of ~0.4 (Specord 200 PLUS UV/Vis spectrophotometer, Analytik Jena). After cultivation for 4 days in cobalt-repleted medium, when cultures reached an OD750nm of ~1.0, cells were harvested from 20 mL culture portions by centrifugation (10 min, 4 °C and 4800 g). Resulting pellets were frozen in liquid nitrogen and stored at -80 °C.

Analysis of R. capsulatus growth, pigmentation, and dry cell weight

Growth of R. capsulatus in expression cultures was measured as optical cell densities at 660 nm, and pigmentation by whole cell absorbance spectra determined semi-daily during the cultivation (GENESYS 10S UV-Vis-Spectrophotometer, Thermo Scientific). For measurement of absorbance spectra, cell samples were washed with water to remove rifampicin in order to avoid interference by the red colored antibiotic, and spectra of whole cells in water were recorded from 300 to 900 nm and cell-density normalized. Carotenoid-specific absorption at 500 nm was extracted from the data.

To determine the correlation of dry cell weight (DCW) of R. capsulatus grown micro-aerobically in RCV medium with antibiotics and serine to measured optical cell density at 660 nm, six differently concentrated samples of cells were harvested from expression cultures by centrifugation. Supernatants were discarded and pellets were dried at 30 °C for one hour under vacuum in a Concentrator 5301 (Eppendorf). Differential weighing of sample tubes before and after filling and drying was employed to determine cell weight. An influence of the drying procedure on the weight of sample tubes was excluded by subjecting empty control tubes to the same treatment. Linear regression enabled determining a correlation of 1 mL OD660nm = 1 with 0.58 mgDCW.

Analysis of Synechocystis growth, pigmentation, and dry cell weight

Growth of Synechocystis cultures was monitored by semi-daily measurement of the OD750nm until the end of the cultivation period. Whole cell absorbance spectra of Synechocystis cells in the clear medium were recorded from 400 to 750 nm on a Specord 200 PLUS spectrophotometer (Analytik Jena) and were normalized to OD750nm. Carotenoid-specific absorption at 500 nm was extracted.

For determination of DCW of Synechocystis, cells of 25 mL of liquid cultures were collected on Supor 0.45-μm membrane filters (Pall) by vacuum filtration and incubated at 60 °C. The filters were weighed before and after cell collection. To exclude effects from desiccation of the filter material, empty control filters were treated the same way. Linear regression enabled determining a correlation of 1 mL OD750nm = 1 with 0.44 mg DCW.

Extraction of triterpenes from R. capsulatus and Synechocystis

For extraction of triterpenes from bacterial cells, all organic solvents (acetone, n-hexane, and chloroform/methanol (2:1)) were supplemented with 0.05% (w/v) butylated hydroxytoluene (BHT). Frozen cell pellets were extracted with acetone that was afore supplemented with a total of 10 nmol β-sitosterol as internal standard. R. capsulatus cell samples (equivalent to OD660nm = 15) were extracted twice with 0.7 mL acetone, Synechocystis cell samples (equivalent to OD750nm = 20) were extracted twice with 1 mL acetone under gentle agitation at 50 °C for 15 min. After centrifugation (3 min at 1780 g and RT), the supernatants were transferred to a fresh tube, supplemented with 1.5 mL 1M NaCl (R. capsulatus) or 1 mL 1M NaCl (Synechocystis) and mixed. Samples were subsequently extracted twice with a total volume of 1.8 ml n-hexane by vigorous mixing for 30 sec. To collect the (upper) hexane phase, samples were centrifuged for phase separation (1 min at 1780 g and RT), and n-hexane extracts were transferred to a fresh tube before drying in a vacuum centrifuge for 20 min at 30 °C. The dried extracts were resuspended in 150 μL chloroform/methanol (2:1), transferred into HPLC vials and stored at -20 °C until analysis.

Extraction of triterpenes from S. cerevisiae

To obtain analytic references for marnerol, hydroxymarnerol and thalianol, lanosterol synthase deficient S. cerevisiae GIL77, carrying pYES/DEST-52 with MRN1, THAS1 or as empty vector control [51], was cultivated in 400 mL batches for plant triterpene production as previously described [52], utilizing galactose supplemented medium for gene expression. After cultivation, 15 mL portions were harvested and pellets stored at -80 °C until extraction. Pellets were incubated with 2 mL 20% (w/v) KOH in 50% (v/v) EtOH for 30 min at 70 °C, prior to extraction of triterpenes with 2 mL n-hexane. 100 μL of the n-hexane phase were dried under a stream of nitrogen gas and re-dissolved in methanol for LC-MS (liquid chromatography–mass spectrometry) analysis.

LC-MS analysis of triterpenes

Extracts from Synechocystis and R. capsulatus, dissolved in chloroform/methanol (2:1) with 0.05% BHT, were diluted 1:100 with methanol (+ 0.05% BHT) prior to injection of 10 μl sample volumes. Triterpenes were separated on a Dionex HPG 3200 HPLC system (Thermo Scientific) equipped with a 150 x 2.1 mm, 2.7 μm, Cortecs C8 column (Waters) with a binary gradient system. Mobile phase A consisted of water + 0.1% formic acid (FA) and mobile phase B consisted of methanol + 0.1% FA. The mobile phase gradient was as follows: Starting conditions were 75% mobile phase B, increased to 85% B within 4 min and then further increased to 100% B in 13.5 min, the plateau was held for another 4.5 min and the system was returned to starting conditions within 3 min. The flow rate was 0.5 mL/min. Triterpenes were analyzed by Q-TOF (quadrupole-time-of-flight) MS and MS/MS on a maXis 4G instrument (Bruker Daltonics) equipped with an ESI (electrospray ionization) source. The instrument was operated in positive-ion mode and the operating conditions were as follows: dry gas (nitrogen): 8.0 L/min, dry heater: 220 °C, nebulizer pressure: 1.8 bar, capillary voltage: 4500 V.

Quantification of triterpenes

Squalene, 2,3-oxidosqualene, cycloartenol and lupeol were quantified using calibration curves obtained with commercial references (Sigma-aldrich). β-sitosterol (Sigma-aldrich) was used as an internal standard. Triterpene production was calculated based on LC-MS data as product titer in the culture (mg L-1), as cell-density normalized specific titer (mg L-1 OD-1), and as specific yield per dry cell weight (mg gDCW-1), all determined as end point values after the cultivation duration of 52 h (R. capsulatus) and 168 h (i.e. 116 h after cobalt-induction; Synechocystis), respectively. In addition, the volumetric and specific productivities (μg L-1 h-1 and μg gDCW-1 h-1) were calculated as averaged values per h cultivation. Squalene productivity in Synechocystis was calculated for the complete cultivation duration (168 h), since squalene production was not induced. Productivity for further triterpenes was calculated in reference to the 116 h period after induction.

Results

Pathway design for triterpene biosynthesis in R. capsulatus and Synechocystis

To implement different plant triterpene biosynthetic pathways in R. capsulatus and Synechocystis, we aimed to link the intrinsic isoprene metabolism of the phototrophic microbes to heterologous pathway modules, introduced as codon optimized genes from A. thaliana (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Triterpene pathway design.

(A) Pathways for targeted biosynthesis of triterpenes in R. capsulatus and Synechocystis by implementation of A. thaliana biosynthesis modules. In both organisms, the common triterpene precursor FPP is provided by the native 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway. To establish the first step of the triterpene-specific precursor module in Synechocystis, a knock-out mutant of the shc gene encoding the squalene-hopene cyclase of the intrinsic hopanoid biosynthesis was employed [37]. In R. capsulatus, heterologous expression of SQS1 was implemented to form squalene. For monooxygenation of squalene to the central precursor 2,3-oxidosqualene, SQE1 was heterologously expressed in both host strains. Subsequently, the cyclization module was designed to convert the linear precursor 2,3-oxidosqualene into plant sterols and further cyclic triterpenes. To cover different triterpene scaffolds, the OSC enzymes CAS1, LUP1, THAS1, and MRN1 from A. thaliana were expressed in each host to synthesize cycloartenol (sterol), lupeol (exhibiting the often occurring pentacyclic scaffold) as well as thalianol and marneral (representing more unusual tri- or monocyclic structures). Respective substrate folding is indicated (CBC, chair-boat-chair; CCC, chair-chair-chair; CB, chair-boat). (B) Schematic representation of expression constructs used for triterpenoid biosynthesis in R. capsulatus and Synechocystis. The respective genes were arranged as synthetic operons, using native ammonium (Pnif) and cobalt (PcoaT) regulated promoters for transcription activation. FPP, farnesyl pyrophosphate; SQS, squalene synthase; SQE, squalene epoxidase; OSC, oxidosqualene cyclase.

https://doi.org/10.1371/journal.pone.0189816.g001

R. capsulatus harbors the 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway leading to the formation of isoprene units and thereof derived FPP, which is a natural intermediate in biosynthesis of the tetraterpenes spheroidene and spheroidenone [20]. We first introduced SQS1 (original cDNA sequence: At4g34640 [57]) encoding squalene synthase 1 in order to redirect the flux of FPP toward heterologous squalene production. In Synechocystis, the MEP pathway provides precursors for the biosynthesis of several carotenoids including β-carotene, echinenone and a set of xanthophylls [58]. Moreover, the pathway is directly linked to native squalene biosynthesis, leading to the formation of hopanoids. For this work, we therefore used the recently described squalene hopene cyclase deficient Δshc mutant, which is characterized by squalene accumulation [37]. For subsequent squalene oxygenation and production of the triterpene precursor 2,3-oxidosqualene, squalene epoxidase 1 encoding SQE1 (At1g58440 [38]) was introduced into both hosts.

The ability of the chosen host organisms to synthesize different triterpene scaffolds was evaluated by introducing four OSC enzymes, namely cycloartenol synthase 1, lupeol synthase 1, thalianol synthase 1 or marneral synthase 1, encoded by CAS1 (At2g07050 [39]), LUP1 (At1g78970 [40]), THAS1 (At5g48010 [41]), and MRN1 (At5g42600 [42]), respectively. While cycloartenol and lupeol represent typical plant triterpenes with tetracyclic plant sterol and pentacyclic plant triterpene scaffolds, respectively, thalianol and marneral exhibit more unusual tri- and monocyclic structures (Fig 1A).

For heterologous expression, synthetic operons were successively assembled encompassing the pathway genes as well as host-specific promoters and regulatory elements (Fig 1B). In R. capsulatus, we used the expression plasmid pRhon5Hi-2, which harbors the native NH4+-depletion inducible Pnif of the nitrogenase-encoding nifHDK operon. In Synechocystis, joint expression of SQE1 and the OSC genes was mediated by the native Co2+-inducible promoter PcoaT of the Co2+ efflux system-encoding coaT gene [59] on plasmid pVZ-spec [54, 60]. The sequences of the synthetic genes and relevant operon elements are given in S2 Table and a graphic representation of the plasmids is shown in S1 Fig.

After plasmid transfer, the R. capsulatus and Synechocystis strains were cultivated in liquid medium as appropriate for respective expression cultures. Cells were harvested in the early stationary (after 52 h, R. capsulatus), or the late logarithmic growth phase (after 168 h, i.e. 116 h after induction, Synechocystis) for LC-MS (liquid chromatography mass-spectrometry) based triterpene analysis.

Implementation of triterpene precursor supply: squalene and 2,3-oxidosqualene

In order to assess the applicability of R. capsulatus and Synechocystis for the biosynthesis of cyclic plant triterpenes, their potentials for accumulating the linear precursors squalene and 2,3-oxidosqualene were determined. To this end, cell extracts of strains Rc_SQS1 and Syn_Δshc were analyzed for squalene accumulation by LC-MS. Specific accumulation of squalene could clearly be determined in specific extracted ion chromatograms (EICs) in both production strains (Fig 2A, two upper panels), and the identification was further validated by comparison of MS fragment spectra (MS/MS spectra) with a commercial standard (S2 Fig). Quantitative analysis revealed comparable specific squalene yields of 9.44 mg gDCW-1 and 9.70 mg gDCW-1 in Rc_SQS1 and Syn_Δshc, respectively. Besides specific yields, also product titers, specific titers and productivities were calculated to enable consideration of respective cell masses, culture volumes, and times required for production (Table 2).

[Figure omitted. See PDF.]

Fig 2. LC-MS detection of triterpenoids produced in engineered R. capsulatus SB1003 (Rc) and Synechocystis sp. PCC 6803 (Syn) strains.

Cells of the respective strains (indicated in chromatograms and color-coded with red for Rc-, green for Syn-strains) were subjected to extraction and analysis after expression of A. thaliana triterpenoid biosynthesis genes (SQS1, squalene synthase; SQE1, squalene epoxidase; CAS1, cycloartenol synthase; LUP1, lupeol synthase; MRN1, marneral synthase). Triterpenoid identity was verified by comparison of retention times, as well as MS and MS/MS spectra to commercial or biological references (S2 Fig). (A) Signals in EICs of m/z 411.399 at RT 14.9 min correspond to squalene (C30H50), signals in EICs of m/z 427.393 at RT 11.8 min correspond to 2,3-oxidosqualene (C30H50O). (B) Signals in EICs of m/z 427.393 at RT 10.9 min correspond to cycloartenol (C30H50O). (C) Signals in EICs of m/z 409.383 at RT 9.9 min correspond to lupeol (C30H50O-H2O). Signals in EICs of m/z 427.393 (and 409.383) at RT 6.7 min correspond to lupX (tentatively identified as lupanediol (C30H52O2-H2O)). In both EICs, signals at RT 11.7 min correspond to 2,3-oxidosqualene (C30H50O). (D) Signals in EICs of m/z 429.409 at RT 10.7 min correspond to marnerol (C30H52O), signals in EICs of m/z 445.404 at RT 7.1 min correspond to hydroxymarnerol (C30H52O2). Signals in EICs of m/z 429.409 at RT 11.8 min correspond to the (M+2) isotopic peak of 2,3-oxidosqualene (C30H52O). Peaks in EICs of m/z 445.404 also detected in Syn_Δshc-SQE1 (control) also correspond to 2,3-oxidosqualene-derived signals. As a reference, chromatograms of samples from S. cerevisiae GIL77 (Sc), carrying pYES/DEST-52 with MRN1 or as empty vector control (EVC), are shown. Shown chromatograms are representative for replicate measurements from at least three independent cultivations. The corresponding quantitative data are summarized in Table 2.

https://doi.org/10.1371/journal.pone.0189816.g002

[Figure omitted. See PDF.]

Table 2. Triterpene levels in engineered R. capsulatus SB1003 (Rc) and Synechocystis sp. PCC 6803 (Syn) strains.

https://doi.org/10.1371/journal.pone.0189816.t002

Moreover, strains Rc_SQS1-SQE1 and Syn_Δshc-SQE1, additionally carrying the epoxidase, were assayed for 2,3-oxidosqualene accumulation which was detected (Fig 2A, two lower panels; S2 Fig) with final yields of 0.32 mg gDCW-1 and 2.09 mg gDCW-1, respectively (Table 2). Although these quantities correspond to only ~3% and ~20% of squalene yields previously observed without SQE1 expression, amounts of the non-oxygenated precursor were reduced in these extracts below the limit of quantification, indicating effective epoxidation of squalene in both strains.

Biosynthesis of the tetracyclic sterol cycloartenol catalyzed by CAS1

Cycloartenol is the common precursor for the biosynthesis of most sterol compounds in plants, which is synthesized via the CBC-type cyclization of 2,3-oxidosqualene resulting in the characteristic tetracyclic arrangement. For biosynthesis of cycloartenol in R. capsulatus and Synechocystis, strains Rc_CAS1-SQS1-SQE1 and Syn_Δshc-SQE1-CAS1 were generated, expressing the precursor module genes together with the CAS1 gene from A. thaliana.

The specific accumulation of cycloartenol in both strains was verified by LC-MS. Cycloartenol was identified by comparison of signal retention times in characteristic EICs (Fig 2B) and MS/MS spectra with a commercial reference (S2 Fig). Quantitative analysis revealed yields of 0.4 mg gDCW-1 and 2.06 mg gDCW-1 for Rc_CAS1-SQS1-SQE1 and Syn_Δshc-SQE1-CAS1, respectively, basically reflecting 2,3-oxidosqualene levels in the preceding analysis of strains not expressing CAS1 (Table 2). At the same time, no significant concomitant accumulation of the precursors squalene and 2,3-oxidosqualene could be detected (compare 2,3-oxidosqualene signals in strains not expressing CAS1 in lower panels of Fig 2A, and absent signals in CAS1 expressing strains in Fig 2B). Therefore, engineering the triterpenoid pathway toward cycloartenol apparently has led to a virtually quantitative cyclization of the linear precursor in both host strains.

Biosynthesis of the pentacyclic triterpene lupeol and a putative derivative by LUP1

As a representative of canonical pentacyclic plant triterpenes, lupeol was chosen as target candidate. Lupeol is synthesized from 2,3-oxidosqualene via the CCC-type cyclization. However, the selected OSC from A. thaliana, LUP1, was shown previously to catalyze the formation of multiple scaffolds besides lupeol, including β-amyrin as well as triterpene alcohols and diols [40, 61]. Thus, the detection of a variety of further minor cyclization products was expected in strains Rc_LUP1-SQS1-SQE1 and Syn_Δshc-SQE1-LUP1. LC-MS analyses of extracts were therefore expanded to the corresponding masses for alternative alcohol and diol compounds. As deducible from the EICs (Fig 2C, two upper panels) and respective MS as well as MS/MS spectra (S2 Fig) in comparison to the lupeol reference, both host strains accumulated lupeol. Here, the EIC of m/z 409.383 was chosen for representation, since this signal, corresponding to [M+H-H2O]+ resulting from loss of H2O (S3 Table) due to in-source decay during electrospray ionization was predominant. While in Rc_LUP1-SQS1-SQE1, amounts were below the limits of quantification, in Syn_Δshc-SQE1-LUP1, the specific yield was 0.29 mg gDCW-1 (Table 2). In addition, a stronger signal for m/z 427.393, accompanied by lower signals for m/z 409.383 and m/z 467.386 was detected deviating in retention time (Fig 2C, lower panels) from the lupeol standard, likely corresponding to another product of LUP1-mediated cyclization. The emerging m/z were calculated to correspond to the [M+H-H2O]+ (427.393), [M+H-(H2O)2]+ (409.383) and [M+Na]+ (467.386) of a compound with the formula C30H52O2. The [M+H]+ was not detected, presumably due to in-source decay of the compound (S2 Fig). The retention time of the compound and deduced formula may point to lupanediol, which has previously been described as a product of LUP1 [61, 62]. However, since the structure of the proposed hydroxylated lupeol derivative cannot be unambiguously assigned from the data, it is designated as lupX at this point. While only trace amounts of both squalene and 2,3-oxidosqualene were detected in Synechocystis extracts, suggesting quantitative precursor cyclization by LUP1, residual 2,3-oxidosqualene in R. capsulatus (0.05 mg gDCW-1) indicates poorer conversion than observed with CAS1.

Expression of THAS1 and MRN1 for biosynthesis of unusual tricyclic and seco-triterpenes

The OSC THAS1 from A. thaliana catalyzes the formation of thalianol via the CCC-type cyclization of 2,3-oxidosqualene, constituting an atypical triterpenoid scaffold, characterized by a tricyclic arrangement [41]. Extracts of strains Rc_THAS1-SQS1-SQE1 and Syn_Δshc-SQE1-THAS1 were analyzed in comparison to samples from a reference yeast strain, i.e. S. cerevisiae GIL77 carrying pYES/DEST-52 with THAS1, which was previously shown to accumulate the product [51]. However, in neither of the engineered strains, the conversion of the substrate to thalianol or any other cyclic product was detected (S2 Fig). In Synechocystis, this is further evident from the observation that the 2,3-oxidosqualene levels were comparable to those the control strains expressing only SQE1. The substrate levels in R. capsulatus co-expressing the OSC were significantly lower than in the control expressing SQS1-SQE1. However, a direct comparison of both strains is not possible here, since in R. capsulatus, the OSC gene was placed upstream of the precursor module genes SQS1 and SQE1 in order to ensure high expression levels of the cyclase, probably at the expense of the downstream genes. The absence of a cyclic triterpene product in THAS1-expressing strains can be caused by diverse problems on multiple levels, including transcription, mRNA stability, translation, protein stability, and activity of the heterologous enzyme. However, this observation most probably points to difficulties at the level of enzyme folding or assembly, since limitations at DNA and transcription levels could be excluded by DNA sequencing as well as RT-qPCR (data not shown) or precursor accumulation (indicating concerted expression of all genes within the synthetic operon THAS1-SQS-SQE).

The OSC MRN1 from A. thaliana catalyzes the formation of the seco-triterpene marneral via an unusual CB substrate conformation, resulting in a further atypical, monocyclic triterpenoid structure [42]. LC-MS analysis revealed that the primary cyclization product marneral was not detected in either the extracts from Rc_MRN1-SQS1-SQE1 or Syn_Δshc-SQE1-MRN1 (data not shown). While in the R. capsulatus extracts none of the expected products were detectable, strain Syn_SQE1-MRN1 apparently accumulated both marnerol (m/z = 429.409) and an additional product with m/z signals (m/z = 427.394, 445.404, 467.386) corresponding to C30H52O2, presumably hydroxymarnerol (Fig 2D, two lower panels and S2 Fig). Both marnerol and the putative hydroxymarnerol were also detected in a reference yeast strain expressing MRN1 [51] (Fig 2D, two upper panels and S2 Fig). However, still significant amounts of 2,3-oxidosqualene were present in the Synechocystis extracts, suggesting low efficiency of substrate cyclization by MRN1.

Cell growth and pigmentation during triterpene biosynthesis

Bacterial growth was monitored during expression experiments to investigate potential effects of triterpene formation on cell viability (Fig 3A). Since all strains reached comparable cell densities during the respective cultivations, cell growth was apparently unimpaired by heterologous gene expression and triterpene biosynthesis.

[Figure omitted. See PDF.]

Fig 3. Comparative physiological characterization of triterpene producing strains of R. capsulatus and Synechocystis.

(A) Growth of R. capsulatus (left) and Synechocystis (right) expressing A. thaliana triterpene biosynthesis genes. Cell densities were determined by turbidity measurements at 660 nm (R. capsulatus) or 750 nm (Synechocystis) to record growth curves. Strains carrying the precursor biosynthetic modules for the generation of squalene and 2,3-oxidosqualene were compared to strains that additionally harbored different cyclization modules or an empty vector as control. Data represent mean values from three independent cultivations, as well as the respective standard deviations. The time point of cobalt (Co2+)-induction in Synechocystis is indicated. The y-axes are scaled logarithmically (base e). (B) Whole cell absorbance spectra of R. capsulatus (upper panel) and Synechocystis (lower panel) expressing the precursor synthesis modules. Pigment profiles were recorded at time point 52 h (R. capsulatus) and 164 h (Synechocystis) of the respective cultivation period (cf. A), and are representative for replicate data obtained from at least three independent cultivations. Typical pigment absorption is indicated. All values were normalized to OD660nm (R. capsulatus) and OD750nm (Synechocystis), respectively. BChl, bacteriochlorophyll a; Car, carotenoids; Chl, chlorophyll a; PC, phycocyanin. (C) Carotenoid content of R. capsulatus (upper panel) and Synechocystis (lower panel). Normalized values at 500 nm were plotted as a measure of carotenoid content in the host cells. R. capsulatus was analyzed at time points 24 h and 52 h during the cultivation. Pigments of Synechocystis were measured before and after cobalt-induction at time points 48 h and 164 h, respectively. Data represent mean results from three independent experiments, as well as the respective standard deviations. EVC, empty vector control; WT, wild type and Δshc both harbor the empty vector.

https://doi.org/10.1371/journal.pone.0189816.g003

In addition, photopigment accumulation was analyzed in strains with precursor modules by recording whole cell absorbance spectra. Both hosts exhibited typical absorbance spectra after cultivation with expected maxima of cellular pigments including bacteriochlorophyll a and spheroidenone for R. capsulatus and chlorophyll a, carotenoids and phycobilins for Synechocystis, indicating that overall pigment synthesis is not affected by transferring the chosen plant triterpene pathways (Fig 3B). To analyze whether triterpene precursor formation results in a decrease of carotenoid (i.e. tetraterpene) formation via FPP depletion, carotenoid-related absorption at λ = 500 was specifically investigated (Fig 3C). Carotenoid-specific absorption remained unchanged in comparison to the control strains suggesting that the introduction of heterologous plant triterpene pathways were supplied by the cellular isoprene pool without interference with native tetraterpene biosynthesis.

Discussion

In this work, we have demonstrated the applicability of the two phototrophic bacteria R. capsulatus SB1003 and Synechocystis sp. PCC 6803 as promising new prokaryotic hosts for the heterologous biosynthesis of diverse cyclic plant triterpene scaffolds.

Basically, the functional heterologous expression of triterpene biosynthesis enzymes requires suitable expression hosts that offer: i) a compatible metabolic background providing the isoprene-derived substrate molecules in high amounts; ii) an appropriate cellular environment including sufficient membrane space, supporting the activity of enzymes, such as OSCs, that favor hydrophobic environments [63–65] as well as storage of hydrophobic intermediates and products; and iii) a broad tolerance toward heterologous biosynthetic products. The heterologous expression of plant triterpene pathway genes in microbial host systems is hitherto almost exclusively analyzed in yeast cells. Only recently, also E. coli has been used as the first bacterial strain for functional expression of an OSC from plants leading to the formation of dammarenediol-II [19]. In this case, however, 2,3-oxidosqualene synthesis further required heterologous co-expression of an NADPH-cytochrome P450 reductase (CPR) in order to efficiently transfer electrons to the employed squalene epoxidase from Methylococcus capsulatus. In contrast, Banta and co-workers recently described an approach using the SQE from Methylomicrobium alcaliphilum to produce 2,3-oxidosqualene as precursor for cyclic triterpenoid synthesis likewise in E. coli, but without CPR co-expression [66]. In this context, our data clearly demonstrate that R. capsulatus and Synechocystis strains expressing plant SQE effectively convert squalene without coexpressing auxiliary reductases. Therefore, the phototrophic host cells likely provide native reductases that efficiently support SQE activity. However, while squalene quantitatively dissipated upon expression of SQE1, 2,3-oxidosqualene did not accumulate to equal amounts, with an especially pronounced decline in product levels in Rhodobacter (~30-fold) as compared to Synechocystis (~5-fold), suggesting either a further conversion by so far unknown enzymes or its export from the cells or a lower stability compared to squalene.

The expression of four different OSC enzymes yielded very different results regarding observed products in the two employed hosts. As for the expression of CAS1, both strains showed effective conversion of 2,3-oxidosqualene to cycloartenol as sole cyclization product, supporting previous reports on yeast strains expressing cycloartenol synthase homologs from different plant species [39,67]. In contrast, LUP1 expression led to formation of lupeol as well as the triterpenoid product lupX with observed m/z argueing for a further oxidized cyclization product in both strains. Previous studies on LUP1 activity in yeast cells suggested the synthesis of a variety of different cyclic products, including 3,20-lupanediol [40,61,62], which would match the here proposed structure of lupX. Furthermore, Krokida and co-workers reported that a closely related OSC from Lotus japonicus, AMY2, converts 2,3-oxidosqualene to both β-amyrin and dihydro-lupeol, without lupeol as intermediate when expressed in Nicotiana benthamiana [68]. In this host, the dihydro-lupeol product was further oxidized to 20-hydroxy-lupeol and 20-hydroxy-betulinic acid by the intrinsic cytochrome P450 (CYP) enzyme CYP71D353. However, when expressed in yeast cells, AMY2 was shown to produce β-amyrin and lupeol [69]. The discrepancies of catalytic products between diverging host systems thus indicates a striking influence of the cellular environment on the nature of the cyclization product(s). In R. capsulatus and Synechocystis, LUP1-mediated formation of oxidized lupeol-derivatives appears to be favored, potentially supported by native, yet undefined CYP acitvities. Likewise, MRN1 expression in Synechocystis leads to the formation of marnerol and–intriguingly–putative hydroxymarnerol. Marnerol accumulation was previously likewise reported for S. cerevisiae strains expressing MRN1, which has been attributed to a spontaneous conversion of the primary product marneral [42,51]. Accumulation of hydroxymarnerol, as here also observed in MRN1 expressing yeast strains (cf. Fig 2D), was previously shown in roots of A. thaliana as a product of the cytochrome P450 type marneral oxidase (MRO) activity [51]. Therefore, the presence of promiscuous oxidative enzyme activity in Synechocystis cells is likely to explain this observation. Indeed, the high potential for native CYP activity in cyanobacterial cells has previously been reported and was further supported by functional studies [70–72]. Among them are CYPs of Nostoc sp. capable of oxygenating sesquiterpenes when heterologously expressed in E. coli [73]. For Synechocystis, so far, oxygenation of C20 retinoids by CYP120A1 (slr0574) was demonstrated in vitro [74]. Notably, no MRN1 cyclization product could be detected in R. capsulatus, further argueing for the reasonability of establishing a variety of microbial platforms for triterpenoid biosynthesis. Moreover, lacking activity of THAS1 in both strains corroborates that successful biosynthesis of cyclic plant triterpenes depends not only on the choice of bacterial host, but also on the specific catalytic properties of the chosen OSC enzyme. In this context, it is worth mentioning that prokaryotic systems appear as yet underrepresented in triterpene production—despite their otherwise thorough exploitation with regard to other terpene classes (recently reviewed in [75–77]). This factum, together with our findings may serve to define plant OSC enzymes as generally rather ‘difficult-to-express’ proteins in bacteria, again underlining the demand for increasing the diversity of expression systems.

As indicated by unimpaired cell growth during heterologous biosynthesis, our host systems represent promising candidates for robust triterpene production, despite previous reports on an antibacterial activity of compounds from this group [78]. Possibly, the enlarged membrane space especially enables photosynthetic bacteria to store certain amounts of terpene precursors and products without significant interference with cellular functions. The two hosts showed rather different production profiles (Table 2). While R. capsulatus features faster growth and thus shorter production times, Synechocystis was capable of accumulating more triterpene products. However, the product levels of 2,3-oxidosqualene and cycloartenol obtained in both host strains as yet are lower than those previously described for 2,3-oxidosqualene and dammarenediol-II in E. coli [19]. Notably, in that study, genes of the FPP building MEP pathway were heterologously overexpressed. Employment of this strategy, typically co-expressing rate-limiting enzymes of the MEP pathway, and/or introducing the entire mevalonate (MVA) pathway has proven highly effective in enhancing bacterial terpene production. In Rhodobacter, such optimization has led to the successful production of terpene compounds to substantial levels [20] as for example Coenzyme Q10, which comprises a C10 isoprene side chain (138 mg L-1 [79]), the C15 sesquiterpene valencene (352 mg L-1 [33]), and non-cyclic C30+ triterpene botryococcene (110 mg L-1 [36]). As for cyanobacteria, similar efforts were successful for establishing significant yields [72] of e.g. C5 isoprene in Synechococcus elongatus (1.26 g L-1 [28]) and C10 β-phellandrene in Synechocystis (10 mg gDCW-1 [80]). In line with the here employed strategy for R. capsulatus, overexpression of squalene synthase might further increase triterpene yields in Synechocystis, where the precursor’s accumulation currently relies on disruption of native hopanoid biosynthesis. Accordingly, Choi and colleagues presented elevated yields of squalene in S. elongatus by combining heterologous SQS expression with MEP pathway optimization [27]. In both hosts, down-regulation of intrinsic carotenoid biosynthesis may represent an attractive strategy to increase the available pool of isoprene building blocks for heterologous terpene production. Recently, Wang and co-workers moreover demonstrated the importance of the terpene synthase step, reporting a 100-fold increase in monoterpene productivity in S. elongatus by enhancing the product sink, which was contrasted with rather marginal improvements by upstream pathway engineering [81]. In addition to metabolite flow into triterpene biosynthetic routes, we identified 2,3-oxidosqualene availability for the pathway as one limiting aspect for product yields–especially in R. capsulatus. Potentially, implementation of downstream OSC enzymes with improved activity for fast 2,3-oxidosqualene conversion, together with strain and process optimization might further increase yields of cyclic triterpene products in the future. In summary, the unique characteristics of photosynthetic bacteria, together with metabolic engineering approaches, may lead to production strains which may be able to produce high levels of triterpenoids from CO2.

Besides providing the chassis for efficient biotechnological production of secondary metabolites, microorganisms can support agricultural processes [82], particularly targeting crop yields and the nutritional value of staple foods [83]. In this context, plant growth promoting bacteria can be applied in order to optimize nutrient supply, phytohormone levels or resistance against biotic as well as abiotic stresses [23,84]. Using Rhodobacter or cyanobacterial species for co-cultivations, e.g. as part of artificial biofilms, has been shown to have the potential to improve e.g. yields and nutritional value of rice crops [22,85–87]. Furthermore, a basic experimental setup for in vitro co-cultivation of rice roots with our model chassis has been established in our lab (data not shown). Engineered biosynthesis of plant-protecting agents like glycosylated triterpenes (saponins) in these bacteria might further contribute to novel strategies in desease management, particularly targeting fungal pathogens. Such strategies further require the heterologous expression of additional product-decorating enzymes [17]. Importantly, the basic suitability of photosynthetic bacteria for functional expression of plant cytochrome P450 and glycosyltransferase enzymes was recently demonstrated using the example of dhurrin pathway reconstitution in Synechocystis 6803 [88]. As an overall future perspective, a vast chemical space of triterpene structures has come into reach, not only by availability of increasing numbers of plant OSC genes, but also by rational enzyme engineering approaches allowing for qualitatively altering product specificities. In this sense, it was recently demonstrated that substrate and product specificities of OSC enzymes can be altered by point mutations of a single amino acid residue close to the active site of the enzyme [89]. In summary, the variety of plant triterpenes harbor a huge potential for different biotechnological applications, particularly addressing the biomedical and agricultural sectors. The here presented host systems shall contribute to realize these applications.

Supporting information

[Figure omitted. See PDF.]

S1 Table. List of primers used in this study.

https://doi.org/10.1371/journal.pone.0189816.s001

(PDF)

S2 Table. Sequence table.

Word document containing DNA sequences of synthetic operons in expression vectors.

https://doi.org/10.1371/journal.pone.0189816.s002

(DOCX)

S3 Table. Sum formulas and calculated m/z of selected triterpenes.

https://doi.org/10.1371/journal.pone.0189816.s003

(PDF)

S1 Fig. Plasmid maps.

Document containing plasmid maps of vectors for expression in R. capsulatus (pRhon5Hi-2, pRhon5Hi-2-SQS1, pRhon5Hi-2-SQS1-SQE1, and pRhon5Hi-2-LUP1-SQS1-SQE1) and of vectors for expression in Synechocystis (pVZ-spec, pVZ-PcoaT-SQE1, and pVZ-PcoaT-SQE1-LUP1).

https://doi.org/10.1371/journal.pone.0189816.s004

(DOCX)

S2 Fig. LC-MS data.

Document containing MS spectra of all analyzed triterpenoids as detected in bacteria; MS/MS spectra of triterpenoids squalene, 2,3-oxidosqualene, cycloartenol, lupeol and lupX as detected in bacteria; LC-MS chromatrograms of extracts from THAS1-expressing bacteria; MS/MS spectra of marnerol/hydroxymarnerol as detected in Synechocystis.

https://doi.org/10.1371/journal.pone.0189816.s005

(DOCX)

Acknowledgments

We thank Nicolas Schmelling and Anna Behle for technical support. Thanks to Pia Lindberg (Uppsala Universitet, Sweden) for kindly providing the Δshc mutant of Synechocystis. We thankfully acknowledge Anne Osbourn (John Innes Centre, Norwich Research Park, UK) for kindly providing S. cerevisiae reference strains carrying plasmids for MRN1 and THAS1 expression.

Citation: Loeschcke A, Dienst D, Wewer V, Hage-Hülsmann J, Dietsch M, Kranz-Finger S, et al. (2017) The photosynthetic bacteria Rhodobacter capsulatus and Synechocystis sp. PCC 6803 as new hosts for cyclic plant triterpene biosynthesis. PLoS ONE 12(12): e0189816. https://doi.org/10.1371/journal.pone.0189816

References

1. Wink M. Modes of Action of Herbal Medicines and Plant Secondary Metabolites. Medicines. 2015;2:251. pmid:28930211

2. Dutton A, Mattiacci L, Amadò R, Dorn S. A novel function of the triterpene squalene in a tritrophic system. J Chem Ecol. 2002;28:103–116. pmid:11868668

3. Babiychuk E, Bouvier-Navé P, Compagnon V, Suzuki M, Muranaka T, Van Montagu M, et al. Allelic mutant series reveal distinct functions for Arabidopsis cycloartenol synthase 1 in cell viability and plastid biogenesis. Proc Natl Acad Sci USA. 2008;105:3163–3168. pmid:18287026

4. Kemen AC, Honkanen S, Melton RE, Findlay KC, Mugford ST, Hayashi K, et al. Investigation of triterpene synthesis and regulation in oats reveals a role for β-amyrin in determining root epidermal cell patterning. Proc Natl Acad Sci USA. 2014;111:8679–8684. pmid:24912185

5. Go YS, Lee SB, Kim HJ, Kim J, Park H-Y, Kim J-K, et al. Identification of marneral synthase, which is critical for growth and development in Arabidopsis. Plant J. 2012;72:791–804. pmid:22882494

6. Delis C, Krokida A, Georgiou S, Peña-Rodríguez LM, Kavroulakis N, Ioannou E, et al. Role of lupeol synthase in Lotus japonicus nodule formation. New Phytol. 2011;189:335–346. pmid:20868395

7. Augustin JM, Kuzina V, Andersen SB, Bak S. Molecular activities, biosynthesis and evolution of triterpenoid saponins. Phytochemistry. 2011;72:435–457. pmid:21333312

8. Chen M, Zhong Z, Tan W, Wang S, Wang Y. Recent advances in nanoparticle formulation of oleanolic acid. Chin Med. 2011;6:20. pmid:21619582

9. Xi J, Tang HY, Zheng Y. Oral dosage forms of oleanolic acid and their pharmacokinetics. Chin J New Drugs. 2009;18:507–515.

10. Yamaguchi H, Noshita T, Kidachi Y, Umetsu H, Hayashi M, Komiyama K, et al. Isolation of ursolic acid from apple peels and its specific efficacy as a potent antitumor agent. J Health Sci. 2008;54:654–660.

11. Řezanka T, Siristova L, Sigler K. Sterols and triterpenoids with antiviral activity. Antiinfect Agents Med Chem. 2009;8:193–210.

12. Díaz AM, Abad MJ, Fernández L, Recuero C, Villaescusa L, Silván AM, et al. In vitro anti-inflammatory activity of iridoids and triterpenoid compounds isolated from Phillyrea latifolia L. Biol Pharm Bull 2000;23:1307–1313. pmid:11085357

13. Szakiel A, Pączkowski C, Pensec F, Bertsch C. Fruit cuticular waxes as a source of biologically active triterpenoids. Phytochem Rev. 2012;11:263–284. pmid:23519009

14. Mannowetz N, Miller MR, Lishko PV. Regulation of the sperm calcium channel CatSper by endogenous steroids and plant triterpenoids. Proc Natl Acad Sci USA. 2017;114(22):5743–5748. pmid:28507119

15. Nakashima T, Inoue T, Oka A, Nishino T, Osumi T, Hata S. Cloning, expression, and characterization of cDNAs encoding Arabidopsis thaliana squalene synthase. Proc Natl Acad Sci USA. 1995;92:2328–2332. pmid:7892265

16. Laden BP, Tang Y, Porter TD. Cloning, heterologous expression, and enzymological characterization of human squalene monooxygenase. Arch Biochem Biophys. 2000;374:381–388. pmid:10666321

17. Thimmappa R, Geisler K, Louveau T, O'Maille P, Osbourn A. Triterpene biosynthesis in plants. Annu Rev Plant Biol. 2014;65:225–257. pmid:24498976

18. Dziggel C, Schäfer H, Wink M. Tools of pathway reconstruction and production of economically relevant plant secondary metabolites in recombinant microorganisms. Biotechnol J. 2017;12:1600145.

19. Li D, Zhang Q, Zhou Z, Zhao F, Lu W. Heterologous biosynthesis of triterpenoid dammarenediol-II in engineered Escherichia coli. Biotechnol Lett. 2016;38:603–609. pmid:26739962

20. Heck A, Drepper T. Engineering Photosynthetic α-Proteobacteria for Microbial Production of Recombinant Proteins and Terpenes. In Hallenbeck P (ed), Modern Topics in the Phototrophic Prokaryotes (Springer International Publishing); 2017. pp. 395–425.

21. Al-Haj L, Lui YT, Abed RM, Gomaa MA, Purton S. Cyanobacteria as Chassis for Industrial Biotechnology: Progress and Prospects. Life (Basel). 2016;6(4). pii:E42.

22. Gamal-Eldin H, Elbanna K. Field evidence for the potential of Rhodobacter capsulatus as biofertilizer for flooded rice. Curr Microbiol. 2011;62:391–395. pmid:20697715

23. Singh S. A review on possible elicitor molecules of cyanobacteria: their role in improving plant growth and providing tolerance against biotic or abiotic stress. J Appl Microbiol. 2014;117:1221–1244. pmid:25069397

24. Armstrong GA, Schmidt A, Sandmann G, Hearst JE. Genetic and biochemical characterization of carotenoid biosynthesis mutants of Rhodobacter capsulatus. J Biol Chem. 1990;265:8329–8338. pmid:2159477

25. Armstrong GA. Eubacteria show their true colors: genetics of carotenoid pigment biosynthesis from microbes to plants. J Bacteriol. 1994;176:4795–4802. pmid:8050991

26. Zhang L, Selão TT, Selstam E, Norling B. Subcellular localization of carotenoid biosynthesis in Synechocystis sp. PCC 6803. PLoS One. 2015;10:e0130904. pmid:26083372

27. Choi SY, Lee HJ, Choi J, Kim J, Sim SJ, Um Y, et al. Photosynthetic conversion of CO2 to farnesyl diphosphate-derived phytochemicals (amorpha-4,11-diene and squalene) by engineered cyanobacteria. Biotechnol Biofuels. 2016;9:202. pmid:27688805

28. Gao X, Gao F, Liu D, Zhang H, Nie X, Yang C. Engineering the methylerythritol phosphate pathway in cyanobacteria for photosynthetic isoprene production from CO2. Energy & Environ Sci. 2016;9:1400–1411.

29. Englund E, Andersen-Ranberg J, Miao R, Hamberger B, Lindberg P. Metabolic engineering of Synechocystis sp. PCC 6803 for production of the plant diterpenoid manoyl oxide. ACS Synth Biol. 2015;4:1270–1278. pmid:26133196

30. 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. pmid:25152894

31. 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:71–86. pmid:24157609

32. Bentley FK, García-Cerdán JG, Chen H-C, Melis A. Paradigm of monoterpene (β-phellandrene) hydrocarbons production via photosynthesis in cyanobacteria. Bioenerg Res. 2013;6:917–929.

33. Beekwilder J, van Houwelingen A, Cankar K, van Dijk AD, de Jong RM, Stoopen G, et al. Valencene synthase from the heartwood of Nootka cypress (Callitropsis nootkatensis) for biotechnological production of valencene. Plant Biotechnol J. 2014;12:174–182. pmid:24112147

34. Loeschcke A, Markert A, Wilhelm S, Wirtz A, Rosenau F, Jaeger KE, et al. TREX: a universal tool for the transfer and expression of biosynthetic pathways in bacteria. ACS Synth Biol. 2013;2:22–33. pmid:23656323

35. Pattanaik B, Lindberg P. Terpenoids and their biosynthesis in cyanobacteria. Life (Basel). 2015;5:269–293.

36. Khan NE, Nybo SE, Chappell J, Curtis WR. Triterpene hydrocarbon production engineered into a metabolically versatile host-Rhodobacter capsulatus. Biotechnol Bioeng. 2015;112:1523–1532. pmid:25728701

37. Englund E, Pattanaik B, Ubhayasekera SJ, Stensjö K, Bergquist J, Lindberg P. Production of squalene in Synechocystis sp. PCC 6803. PLoS One. 2014;9:e90270. pmid:24625633

38. Rasbery JM, Shan H, LeClair RJ, Norman M, Matsuda SP, Bartel B. Arabidopsis thaliana squalene epoxidase 1 is essential for root and seed development. J Biol Chem. 2007;282:17002–17013. pmid:17426032

39. Corey EJ, Matsuda SP, Bartel B. Isolation of an Arabidopsis thaliana gene encoding cycloartenol synthase by functional expression in a yeast mutant lacking lanosterol synthase by the use of a chromatographic screen. Proc Natl Acad Sci USA. 1993;90:11628–11632. pmid:7505443

40. Herrera JB, Bartel B, Wilson WK, Matsuda SP. Cloning and characterization of the Arabidopsis thaliana lupeol synthase gene. Phytochemistry. 1998;49:1905–1911. pmid:9883589

41. Fazio GC, Xu R, Matsuda SP. Genome mining to identify new plant triterpenoids. J Am Chem Soc. 2004;126:5678–5679. pmid:15125655

42. Xiong Q, Wilson WK, Matsuda SP. An Arabidopsis oxidosqualene cyclase catalyzes iridal skeleton formation by Grob fragmentation. J Am Chem Soc. 2006;45:1285–1288.

43. Hanahan D. Studies on transformation of Escherichia coli with plasmids. J Mol Biol. 1983;166:557–580. pmid:6345791

44. Elhai J, Wolk CP. Conjugal transfer of DNA to cyanobacteria. Methods Enzymol. 1988;167:747–754. pmid:3148842

45. Simon R, Priefer U, Pühler A. A broad host range mobilization system for in vivo genetic engineering: transposon mutagenesis in gram negative bacteria. Nat Biotechnol. 1983;1: 784–791.

46. Strnad H, Lapidus A, Paces J, Ulbrich P, Vlcek C, Paces V, et al. Complete genome sequence of the photosynthetic purple nonsulfur bacterium Rhodobacter capsulatus SB 1003. J Bacteriol. 2010;192:3545–3546. pmid:20418398

47. Klipp W, Masepohl B, Pühler A. Identification and mapping of nitrogen fixation genes of Rhodobacter capsulatus: duplication of a nifA-nifB region. J Bacteriol. 1988;170:693–699. pmid:2828320

48. Weaver PF, Wall JD, Gest H. Characterization of Rhodopseudomonas capsulata. Arch Microbiol. 1975;105:207–216. pmid:1103769

49. Rippka R, Deruelles J, Waterbury JB, Herdman M, Stanier RY. Generic assignments, strain histories and properties of pure cultures of cyanobacteria. J Gen Microbiol. 1979;111:1–61.

50. Gollub EG, Liu K-P, Dayan J, Adlersberg M, Sprinson DB. Yeast mutants deficient in heme biosynthesis and a heme mutant additionally blocked in cyclization of 2,3-oxidosqualene. J Biol Chem. 1977;252:2846–2854. pmid:323256

51. Field B, Fiston-Lavier AS, Kemen A, Geisler K, Quesneville H, Osbourn AE. Formation of plant metabolic gene clusters within dynamic chromosomal regions. Proc Natl Acad Sci USA. 2011;108:16116–16121. pmid:21876149

52. Kushiro T, Shibuya M, Ebizuka Y. β-amyrin synthase-cloning of oxidosqualene cyclase that catalyzes the formation of the most popular triterpene among higher plants. Eur J Biochem. 1998;256:238–244. pmid:9746369

53. Katzke N, Arvani S, Bergmann R, Circolone F, Markert A, Svensson V, et al. A novel T7 RNA polymerase dependent expression system for high-level protein production in the phototrophic bacterium Rhodobacter capsulatus. Protein Expr Purif. 2010;69:137–146. pmid:19706327

54. Mitschke J, Georg J, Scholz I, Sharma CM, Dienst D, Bantscheff J, et al. An experimentally anchored map of transcriptional start sites in the model cyanobacterium Synechocystis sp. PCC6803. Proc Natl Acad Sci USA. 2011;108:2124–2129. pmid:21245330

55. Zinchenko VV, Piven IV, Melnik VA, Shestakov SV. Vectors for the complementation analysis of cyanobacterial mutants. Russ J Genet. 1999;35:228–232.

56. Beyer HM, Gonschorek P, Samodelov SL, Meier M, Weber W, Zurbriggen MD. AQUA cloning: a versatile and simple enzyme-free cloning approach. PLoS One. 2015;10:e0137652. pmid:26360249

57. Busquets A, Keim V, Closa M, del Arco A, Boronat A, Arró M, et al. Arabidopsis thaliana contains a single gene encoding squalene synthase. Plant Mol Biol. 2008;67:25–36. pmid:18236008

58. Takaichi S, Mochimaru M. Carotenoids and carotenogenesis in cyanobacteria: unique ketocarotenoids and carotenoid glycosides. Cell Mol Life Sci. 2007;64:2607–2619. pmid:17643187

59. Peca L, Kós PB, Vass I. Characterization of the activity of heavy metal-responsive promoters in the cyanobacterium Synechocystis PCC 6803. Acta Biol Hung. 2007;58 Suppl:11–22.

60. Georg J, Dienst D, Schürgers N, Wallner T, Kopp D, Stazic D, et al. The small regulatory RNA SyR1/PsrR1 controls photosynthetic functions in cyanobacteria. Plant Cell. 2014;26:3661–3679. pmid:25248550

61. Segura MJ, Meyer MM, Matsuda SP. Arabidopsis thaliana LUP1 converts oxidosqualene to multiple triterpene alcohols and a triterpene diol. Org Lett. 2000;2:2257–2259. pmid:10930257

62. Kushiro T, Hoshino M, Tsutsumi T, Kawai K, Shiro M, Shibuya M, et al. Stereochemical course in water addition during LUP1-catalyzed triterpene cyclization. Org Lett. 2006;8:5589–5592. pmid:17107079

63. Liang Y-L, Zhao S-J, Xu L-X, Zhang X-Y. Heterologous expression of dammarenediol synthase gene in an engineered Saccharomyces cerevisiae. Lett Appl Microbiol. 2012;55:323–329. pmid:22897704

64. Milla P, Athenstaedt K, Viola F, Oliaro-Bosso S, Kohlwein SD, Daum G, et al. Yeast oxidosqualene cyclase (Erg7p) is a major component of lipid particles. J Biol Chem. 2002;277:2406–2412. pmid:11706015

65. Milla P, Viola F, Oliaro-Bosso S, Rocco F, Cattel L, Joubert BM, et al. Subcellular localization of oxidosqualene cyclases from Arabidopsis thaliana, Trypanosoma cruzi, and Pneumocystis carinii expressed in yeast. Lipids. 2002;37:1171–1176. pmid:12617471

66. Banta AB, Wei JH, Gill CC, Giner J-L, Welander PV. Synthesis of arborane triterpenols by a bacterial oxidosqualene cyclase. Proc Natl Acad Sci USA. 2016;114:245–250. pmid:28028245

67. Hayashi H, Hiraoka N, Ikeshiro Y, Kushiro T, Morita M, Shibuya M, et al. Molecular cloning and characterization of a cDNA for Glycyrrhiza glabra cycloartenol synthase. Biol Pharm Bull. 2000;23:231–234. pmid:10706391

68. Krokida A, Delis C, Geisler K, Garagounis C, Tsikou D, Peña-Rodríguez LM, et al. A metabolic gene cluster in Lotus japonicus discloses novel enzyme functions and products in triterpene biosynthesis. New Phytol. 2013;200:675–690. pmid:23909862

69. Iturbe-Ormaetxe I, Haralampidis K, Papadopoulou K, Osbourn AE. Molecular cloning and characterization of triterpene synthases from Medicago truncatula and Lotus japonicus. Plant Mol Biol. 2003;51:731–743. pmid:12683345

70. Robert FO, Pandhal J, Wright PC. Exploiting cyanobacterial P450 pathways. Curr Opin Microbiol. 2010;13:301–306. pmid:20299274

71. Ke N, Baudry J, Makris TM, Schuler MA, Sligar SG. A retinoic acid binding cytochrome P450: CYP120A1 from Synechocystis sp. PCC 6803. Arch Biochem Biophys. 2005;436:110–120. pmid:15752715

72. Xue Y, He Q. Cyanobacteria as cell factories to produce plant secondary metabolites. Front Bioeng Biotechnol. 2015;3:57. pmid:25973419

73. Agger SA, Lopez-Gallego F, Hoye TR, Schmidt-Dannert C. Identification of sesquiterpene synthases from Nostoc punctiforme PCC 73102 and Nostoc sp. strain PCC 7120. J Bacteriol. 2008;190:6084–6096. pmid:18658271

74. Alder A, Bigler P, Werck-Reichhart D, Al-Babili S. In vitro characterization of Synechocystis CYP120A1 revealed the first nonanimal retinoic acid hydroxylase. FEBS J. 2009;276:5416–5431. pmid:19703230

75. Jongedijk E, Cankar K, Buchhaupt M, Schrader J, Bouwmeester H, Beekwilder J. Biotechnological production of limonene in microorganisms. Appl Microbiol Biotechnol. 2016;100:2927–2938. pmid:26915992

76. Xie M, Wang W, Zhang W, Chen L, Lu X. Versatility of hydrocarbon production in cyanobacteria. Appl Microbiol Biotechnol. 2017;101:905–919. pmid:28032195

77. Majdi M, Ashengroph M, Abdollahi MR. Sesquiterpene lactone engineering in microbial and plant platforms: parthenolide and artemisinin as case studies. Appl Microbiol Biotechnol. 2016;100:1041–1059. pmid:26567019

78. Chung PY, Navaratnam P, Chung LY. Synergistic antimicrobial activity between pentacyclic triterpenoids and antibiotics against Staphylococcus aureus strains. Ann Clin Microbiol Antimicrob. 2011;10:25. pmid:21658242

79. Lu W, Ye L, Xu H, Xie W, Gu J, Yu H. Enhanced production of coenzyme Q10 by self-regulating the engineered MEP pathway in Rhodobacter sphaeroides. Biotechnol Bioeng. 2014;111:761–769. pmid:24122603

80. Formighieri C, Melis A. Sustainable heterologous production of terpene hydrocarbons in cyanobacteria. Photosynth Res. 2016;130: 123–135. pmid:26895437

81. Wang X, Liu W, Xin C, Zheng Y, Cheng Y, Sun S, et al. Enhanced limonene production in cyanobacteria reveals photosynthesis limitations. Proc Natl Acad Sci USA. 2016;113:14225–14230. pmid:27911807

82. Moshelion M, Altman A. Current challenges and future perspectives of plant and agricultural biotechnology. Trends Biotechnol. 2015;33:337–342. pmid:25842169

83. McKersie B. Planning for food security in a changing climate. J Exp Bot. 2015;66:3435–3450. pmid:25614663

84. Lugtenberg B, Kamilova F. Plant-growth-promoting rhizobacteria. Annu Rev Microbiol. 2009;63:541–556. pmid:19575558

85. Priya H, Prasanna R, Ramakrishnan B, Bidyarani N, Babu S, Thapa S, et al. Influence of cyanobacterial inoculation on the culturable microbiome and growth of rice. Microbiol Res. 2015;171:78–89. pmid:25644956

86. Adak A, Prasanna R, Babu S, Bidyarani N, Verma S, Pal M, et al. Micronutrient enrichment mediated by plant-microbe interactions and rice cultivation practices. J Plant Nutr. 2016;39:1216–1232.

87. Prasanna R, Adak A, Verma S, Bidyarani N, Babu S, Pal M, et al. Cyanobacterial inoculation in rice grown under flooded and SRI modes of cultivation elicits differential effects on plant growth and nutrient dynamics. Ecol Eng. 2015;84:532–541.

88. Wlodarczyk A, Gnanasekaran T, Nielsen AZ, Zulu NN, Mellor SB, Luckner M, et al. Metabolic engineering of light-driven cytochrome P450 dependent pathways into Synechocystis sp. PCC 6803. Metab Eng. 2016;33:1–11. pmid:26548317

89. Salmon M, Thimmappa RB, Minto RE, Melton RE, Hughes RK, O'Maille PE, et al. A conserved amino acid residue critical for product and substrate specificity in plant triterpene synthases. Proc Natl Acad Sci USA. 2016;113:E4407–4414. pmid:27412861