A new synthetic biology approach allows transfer of an entire metabolic pathway from a medicinal plant to a biomass crop

  1. Paulina Fuentes
  2. Fei Zhou
  3. Alexander Erban
  4. Daniel Karcher
  5. Joachim Kopka
  6. Ralph Bock  Is a corresponding author
  1. Max-Planck-Institut für Molekulare Pflanzenphysiologie, Germany
7 figures and 3 tables

Figures

Metabolic pathway of artemisinin biosynthesis.

The canonical pathway of artemisinin synthesis starts with the conversion of IPP/DMAPP (C5 isoprenoids produced by the MVA pathway in the cytosol or the MEP pathway in the chloroplast) into farnesyl pyrophosphate (FPP), catalyzed by farnesyl pyrophosphate synthase (FPS). Amorpha-4,11-diene synthase (ADS) converts FPP into amorpha-4,11-diene in the first committed step of the pathway. Amorpha-4,11-diene is then successively oxidized to artemisinic alcohol, artemisinic aldehyde and artemisinic acid by the cytochrome P450 monooxygenase CYP71AV1 (CYP) and its redox partner, the cytochrome P450 reductase (CPR). In A. annua, artemisinic aldehyde is converted to dihydroartemisinic aldehyde by DBR2, and then to dihydroartemisinic acid by ALDH1. Artemisinin is generated by the spontaneous oxidation of dihydroartemisinic acid in planta, and can be produced by chemical conversion of artemisinic acid in vitro. Enzymes depicted in red improve the efficiency of different oxidation steps in yeast (Paddon et al., 2013; Paddon and Keasling, 2014). See text for details.

https://doi.org/10.7554/eLife.13664.003
Implementation of the canonical pathway of artemisinic acid biosynthesis in chloroplasts.

Synthetic codon-optimized genes for the four enzymes required to produce artemisinic acid (Figure 1) were introduced into the tobacco plastid genome by stable genetic transformation with four different synthetic operon constructs (pAO1-4). The constructs differ in gene arrangement and in the translation signals that drive synthesis of the key pathway enzyme (ADS) catalyzing the first committed step. (A) Physical map of the plastid genome region (ptDNA) used for integration of the synthetic artemisinic acid operons and maps of the transgenic loci in the generated transplastomic tobacco lines (Nt-AO1-4). The artemisinic acid operon genes are depicted as light blue boxes. Chloroplast promoters and terminators are shown in green, the aadA selectable marker gene for chloroplast transformation is represented as a white box, and genes in flanking plastid sequences used for transgene targeting via homologous recombination are in black. Genes above the line are transcribed from left to right, genes below the line are transcribed in the opposite direction. The four transgenes are arranged in two dicistronic operons. FPS and CYP are driven by the Chlamydomonas reinhardtii plastid ribosomal RNA operon promoter (Cr Prrn) and the g10 leader sequence from phage T7 (T7 Lg10). The second operon containing ADS and CPR is driven by the C. reinhardtii psbA promoter (Cr PpsbA) and either the T7 Lg10 or the psbA leader sequence from C. reinhardtii (Cr LpsbA). This operon is arranged either in sense and downstream of the first operon (AO1, 3) or in antisense, downstream of the aadA cassette (AO2, 4). The genes in each operon are separated by an intercistronic expression element (IEE) conferring intercistronic RNA processing and, in this way, enhancing expression of downstream cistrons of the operon (Zhou et al., 2007; Drechsel and Bock, 2010). The BamHI restriction sites used in RFLP analyses and the expected fragment sizes are indicated. The location of the hybridization probe is shown as a black bar. Cr: C. reinhardtii; Nt: N. tabacum; T7: bacteriophage T7; P: promoter; L: leader sequence; T: terminator; SD: Shine-Dalgarno sequence. (B) RFLP analysis of transplastomic plants. Two independently isolated transplastomic lines are shown for constructs pAO1-3 and one for pAO4. (C) Seed assays confirming the homoplasmic state of the transplastomic plants. Seeds were germinated on medium containing 500 mg/L spectinomycin (Nt-AO2-1, Nt-AO3-1, Nt-wt) or antibiotic-free medium (Nt-wt).

https://doi.org/10.7554/eLife.13664.004
Figure 3 with 2 supplements
Phenotype of transplastomic tobacco plants and accumulation of artemisinic compounds.

(A) Transplastomic lines Nt-AO1-1 and Nt-AO3-1 display a slightly pale and growth-delayed phenotype at the juvenile stage. WAT: weeks after transfer from tissue culture to soil; scale bars: 10 cm. (B) Amorpha-4,11-diene is synthesized in all transplastomic lines, but accumulates to lower levels in the lines displaying an altered phenotype (purple bars). (C) Artemisinic alcohol is detected in similar amounts in all transplastomic plants. (D) Accumulation of artemisinic acid correlates with the altered phenotype of Nt-AO1-1 and Nt-AO3-1. Relative accumulation of amorpha-4,11-diene was profiled by GC-MS analysis of volatile organic compounds (VOCs). Relative accumulation of the sum of free and conjugated artemisinic alcohol and artemisinic acid were determined by GC-MS analysis of the soluble metabolite fraction after saponification (see Materials and methods; Figures 6 and 7). In agreement with previous reports (van Herpen et al., 2010), these compounds were found to be present mainly as conjugates. Expanding leaves of 5–6 plants per line were used for each measurement. Error bars represent the SD. Different letters above the bars indicate significant differences as determined by One-way ANOVA (p<0.001) and the Holm-Sidak post-hoc test.

https://doi.org/10.7554/eLife.13664.005
Figure 3—figure supplement 1
Phenotypes of Nt-AO2-1 and Nt-AO3-1 plants throughout development.

Six plants per line (Nt-wt, Nt-AO2-1 and Nt-AO3-1) were grown under standard greenhouse conditions and photographs were taken of one representative plant per line at different time points: young plants (before flowering, stage 1), flowering plants (stage 2) and old plants (seed capsules formed, stage 3). Light-green leaves and slightly delayed growth of line Nt-AO3-1 are more evident at the young stage. At later stages, all Nt-AO lines display a wild type-like phenotype and produce viable seeds in indistinguishable amounts. y: young leaf; i: expanding (intermediate) leaf; m: fully expanded (mature) leaf.

https://doi.org/10.7554/eLife.13664.006
Figure 3—figure supplement 2
Isoprenoids levels throughout development in wild-type Nicotiana tabacum plants (Nt-wt) and the transplastomic lines Nt-AO2-1 and Nt-AO3-1.

Plants were grown under standard greenhouse conditions and samples were taken from young (y), expanding (i), and fully expanded (m) leaves at three developmental stages (1–3; cf. Figure 3—figure supplement 1). Metabolite levels were determined by UPLC analysis. The values represent the peak height for each compound divided by 10.000 and normalized to the fresh weight (FW), resulting in the normalized response/FW. Error bars represent the SD (n = 3 plants per line). Different letters above the bars indicate significant differences as determined by One-way ANOVA (p<0.05) and the Holm-Sidak post-test.

https://doi.org/10.7554/eLife.13664.007
Figure 4 with 1 supplement
Production of artemisinic acid is maintained throughout plant development and correlates with a high CYP/CPR expression ratio.

Artemisinic compounds and expression levels of the transgenes were measured in young (stage 1), flowering (stage 2) and old plants (stage 3; see Figure 3—figure supplement 1). (A) Amorpha-4,11-diene accumulates to higher levels in line Nt-AO2-1 than in Nt-AO3-1. (B) Artemisinic alcohol is present at similar levels in early (1) and late stages (3) of development in lines Nt-AO2-1 and Nt-AO3-1, but it is slightly higher in the flowering stage (2) of line Nt-AO3-1. (C) Artemisinic acid accumulates to high levels during all developmental stages of line Nt-AO3-1, whereas in line Nt-AO2-1, it is detectable only in mature leaves of young and flowering plants. Relative accumulation of amorpha-4,11-diene and artemisinic alcohol was profiled, the tissue content of artemisinic acid was quantified using an authenticated reference standard (n = 5–6 plants per line; Figures 6 and 7). The sum of free and conjugated artemisinic alcohol and artemisinic acid were determined. y: young leaf; i: expanding (intermediate) leaf; m: fully expanded (mature) leaf. Error bars represent the SD. (D–G) Northern blot analysis of the expression of the four transgenes. Total RNA samples from N. tabacum wild-type (Nt-wt) plants and the transplastomic lines Nt-AO2-1 and Nt-AO3-1 (at the developmental stages 1–3) were separated in denaturing 1.5% agarose gels, blotted and hybridized to strand-specific RNA probes. Below each blot, the rRNA-containing region of the ethidium bromide-stained gel prior to blotting is shown as a control for RNA integrity and equal loading. The Nt-wt sample corresponds to RNA extracted from a fully expanded leaf of a N. tabacum wild-type plant at developmental stage 2. The smallest labeled band in each blot corresponds to the monocistronic mRNA. Larger bands represent unprocessed polycistronic precursor transcripts and read-through transcripts (which are common in plastids; e.g., Elghabi et al., 2011; Lu et al., 2013). CYP transcripts accumulate to higher levels in line Nt-AO3-1, while CPR transcripts accumulate to higher levels in line Nt-AO2-1, resulting in a higher CYP/CPR expression ratio in line Nt-AO3-1.

https://doi.org/10.7554/eLife.13664.008
Figure 4—figure supplement 1
Quantitation of the expression of CYP and CPR by qRT-PCR analysis.

Samples from young (y), expanding (intermediate; i) and fully expanded (mature; m) leaves were measured in three technical replicates for each line, in early, flowering and late developmental stages (1–3; cf. Figure 3—figure supplement 1). The wild-type sample (wt) corresponds to a sample from a fully expanded leaf of N. tabacum cv. Petit Havana at the flowering stage. The qRT-PCR data confirm the northern blot analyses that had revealed a higher CYP to CPR expression ratio in transplastomic line Nt-AO3-1 than in line Nt-AO2-1.

https://doi.org/10.7554/eLife.13664.009
Figure 5 with 2 supplements
Isolation of combinatorially supertransformed transplastomic lines with a strong increase in artemisinic acid accumulation.

Transplastomic lines Nt-AO2-1 and Nt-AO3-1 were combinatorially supertransformed with genes for additional enzymes of the pathway (Figure 1) to facilitate large-scale screening for increased artemisinic acid production. (A) Relative artemisinic acid levels (given in response/FW; see Materials and methods) in the T1 generation of two combinatorially supertransformed lines obtained with transplastomic recipient line Nt-AO2-1 (Nt-AO2-CS) and eight lines obtained with transplastomic recipient line Nt-AO3-1 (Nt-AO3-CS). An up to 77-fold increase in artemisinic acid was achieved (line Nt-AO3-CS180) in comparison to recipient line Nt-AO3-1. For a complete list of screened supertransformed lines, see Figure 5—source data 1. (B) Inverse relationship between artemisinic acid accumulation and artemisinic alcohol accumulation in supertransformed lines. Fully expanded leaves of 5–6 plants per line (at the flowering stage) were used for metabolite profiling. The sum of free and conjugated artemisinic alcohol and artemisinic acid were determined. (C) qRT-PCR analysis of transgene expression suggests a predominant role of ALDH1 in boosting artemisinic acid synthesis. 2–3 plants per line were measured, and the expression levels were ranked after One-way ANOVA comparison (p<0.05). Brown color indicates absence of gene expression. (D) Combinatorially supertransformed lines with a high increase in artemisinic acid (Nt-AO2-CS132 and Nt-AO3-CS180) display a similar phenotype as the corresponding transplastomic recipient line. WAT: weeks after transfer from tissue culture to soil; scale bars: 10 cm.

https://doi.org/10.7554/eLife.13664.010
Figure 5—source data 1

Metabolic screening (phenotyping) and genotyping of the T0 generation of combinatorially supertransformed Nt-AO-CS lines.

The supertransformed lines Nt-AO2-CS and Nt-AO3-CS are arranged according to their artemisinic acid content, from low to high. Fresh weight corrected response values (R/FW) for amorpha-4,11-diene were multiplied by 1,000 and expressed as Rx1000/FW. Dihydroartemisinic acid was only detectable in line Nt-AO3-CS180 at a low level of 0.03 R/FW. Asterisks mark the selected candidate lines further analyzed in the T1 generation. Lines were clustered according to artemisinic acid content using hierarchical cluster analysis based on Ward’s method. Cnd: cluster nd (artemisinic acid not detected); C1: cluster 1; C2: cluster 2; C3: cluster 3; C4: cluster 4; C5: cluster 5. Genes detected in the genomic PCR assays are numbered as follows: 1: dxr; 2: CYB5; 3: ADH1; 4: ALDH1; 5: DBR2; 0: no gene detected; ?: unclear result. nd: not detected; nm: not measured; -: not determined. Amorpha-4,11-diene values are from one measurement per line, values for artemisinic alcohol, dihydroartemisinic alcohol and artemisinic acid represent averages of three technical replicates per line. SD: standard deviation.

https://doi.org/10.7554/eLife.13664.011
Figure 5—figure supplement 1
Amorpha-4,11-diene, dihydroartemisinic alcohol and dihydroartemisinic acid accumulation in the T1 generation of combinatorially supertransformed plants.

Levels of amorpha-4,11-diene in Nt-AO-CS lines do not correlate with the strong increase in artemisinic acid (Figure 5). Dihydroartemisinic alcohol accumulates to higher levels in lines Nt-AO2-CS95 and Nt-AO3-CS145. Dihydroartemisinic acid is detectable only in two of the best-performing COSTREL lines. Amorpha-4,11-diene levels were determined by GC-MS of VOCs. Dihydroartemisinic alcohol and dihydroartemisinic acid were measured by GC-MS analysis of soluble saponified metabolites. The sum of free and conjugated artemisinic compounds was determined. Error bars represent the SD (n = 5–6 plants per line).

https://doi.org/10.7554/eLife.13664.012
Figure 5—figure supplement 2
Measurements of plant height and total leaf biomass of COSTREL line Nt-AO3-CS180 (progeny of four different T1 lines), its transplastomic progenitor line Nt-AO3-1 and the wild type (wt).

(A) Comparison of plant height at the same age. Plants were measured when the wild type started to flower. (B) Comparison of plant height at the same developmental stage. Transplastomic line Nt-AO3-1 and COSTREL line Nt-AO3-CS180 were measured five days later than the wild type to compensate for their slightly delayed onset of flowering. (C) Comparison of total leaf biomass (fresh weight, FW) at the same plant age. Plants were measured when the wild type started to flower. (D) Comparison of total leaf biomass at the same developmental stage. Transplastomic line Nt-AO3-1 and COSTREL line Nt-AO3-CS180 were measured five days later than the wild type to compensate for their delayed flowering. Note that, at the same plant age, transplastomic line Nt-AO3-1 and COSTREL line Nt-AO3-CS180 are slightly shorter and produce less total leaf biomass than the wild type. Once the transplastomic plants and the COSTREL lines flower (same developmental stage, five days later), the COSTREL plants reach a height and a total leaf biomass that is close to the values measured for the wild type. On average, the COSTREL plants are 7% shorter and produce 13% less leaf biomass than wild-type plants. No significant difference in either height or total leaf FW was observed between COSTREL plants and their transplastomic progenitor line Nt-AO3-1. Error bars represent the SD (n = 6). Different letters above the bars indicate significant differences as determined by One-way ANOVA (p<0.05) and the Holm-Sidak post-hoc test.

https://doi.org/10.7554/eLife.13664.013
Chromatograms and mass spectra of amorpha-4,11-diene and artemisinic alcohol.

Characteristic peaks for one specific fragment at the expected retention time or index are displayed for each compound. (A) Amorpha-4,11-diene-specific mass feature 119 at a retention time of 1564 s. This metabolite is present in all Nt-AO lines, and at slightly higher levels in lines Nt-AO2-1 and Nt-AO4-1. (B) Artemisinic alcohol-specific mass feature 202 at a retention index of 1784. The compound is present at similar levels in all Nt-AO lines. Both compounds are absent from the wild-type sample. In addition to the chromatograms, the characteristic mass spectrum (m/z) of each compound is shown for the standard and for one of the artemisinic acid operon lines. EPY224: yeast strain that produces amorpha-4,11-diene (Ro et al., 2006). One representative plant per line is depicted. Mass spectra and mass features of trimethylsilylated artemisinic alcohol are shown.

https://doi.org/10.7554/eLife.13664.016
Chromatograms and mass spectra of the artemisinic compounds artemisinic acid, dihydroartemisinic alcohol and dihydroartemisinic acid.

Characteristic peaks for one specific fragment at the expected retention index are shown for each compound. (A) Artemisinic acid-specific mass feature 216 is shown at a retention index of 1850. This compound accumulates to higher levels in lines Nt-AO1-1 and Nt-AO3-1. (B) Dihydroartemisinic alcohol-specific mass feature 162 at a retention index of 1789. The compound is present at high levels in COSTREL line Nt-AO2-CS95, but is absent from transplastomic line Nt-AO2-1. (C) Dihydroartemisinic acid-specific mass feature 163 at a retention index of 1859. This compound accumulates in COSTREL line Nt-AO3-CS180, but is absent from transplastomic line Nt-AO3-1. All compounds are absent from the wild-type sample. In addition to the chromatograms, the characteristic mass spectrum of each compound is shown for the standard and for one of the artemisinic acid operon lines. Mass spectra and mass features of trimethylsilylated artemisinic compounds are shown.

https://doi.org/10.7554/eLife.13664.017

Tables

Table 1

Correlation analysis of artemisinic compounds and transgenes introduced into transplastomic line Nt-AO2-1 by combinatorial supertransformation. The levels of the artemisinic compounds amorpha-4,11-diene, artemisinic alcohol, dihydroartemisinic alcohol and artemisinic acid, and the presence of the transgenes dxr, CYB5, ADH1, ALDH1 and DBR2 were correlated using Spearman’s method in the 39 Nt-AO2-CS lines analyzed by genomic PCR in the T0 generation (see Figure 5—source data 1) using the SPSS software. Dihydroartemisinic acid was excluded from this analysis, because it was not detectable in any of the Nt-AO2-CS lines. CC: correlation coefficient. Positive values indicate positive correlations and negative values indicate negative correlations. *: p<0.05; **: p<0.01; N: number of samples where both variables are present.

https://doi.org/10.7554/eLife.13664.014
Nt-AO2-CSAmorpha-4,11-dieneArtemisinic alcoholDihydroartemi-sinic alcoholArtemisinic aciddxrCYB5ADH1ALDH1DBR2
Amorpha-4,11-dieneCC1.0000.602**0.147-0.1120.067-0.020-0.009-0.0150.156
N352913243535353535
Art. alcoholCC0.602**1.000-0.557*-0.141-0.159-0.165-0.1700.013-0.217
N292913242929292929
Dihydroart. alcoholCC0.147-0.557*1.0000.642*0.4620.2480.3830.1800.496
N131313131313131313
Art. acidCC-0.112-0.1410.642*1.0000.4040.1590.1790.2070.317
N242413242424242424
dxrCC0.067-0.1590.4620.4041.0000.643**0.402*0.546**0.578**
N352913243939393939
CYB5CC-0.020-0.1650.2480.1590.643**1.0000.1920.507**0.793**
N352913243939393939
ADH1CC-0.009-0.1700.3830.1790.402*0.1921.0000.372*0.270
N352913243939393939
ALDH1CC-0.0150.0130.1800.2070.546**0.507**0.372*1.0000.420**
N352913243939393939
DBR2CC0.156-0.2170.4960.3170.578**0.793**0.2700.420**1.000
N352913243939393939
Table 2

Correlation analysis of artemisinic compounds and transgenes introduced into transplastomic line Nt-AO3-1 by combinatorial supertransformation. The levels of the artemisinic compounds amorpha-4,11-diene, artemisinic alcohol, dihydroartemisinic alcohol and artemisinic acid, and the presence of the transgenes dxr, CYB5, ADH1, ALDH1 and DBR2 were correlated using Spearman’s method in the 61 Nt-AO3-CS lines analyzed by genomic PCR in the T0 generation (see Figure 5—source data 1) using the SPSS software. Dihydroartemisinic acid had to be excluded from this analysis, because it was detectable only in one of the Nt-AO3-CS lines in the T0 generation. Note that the negative correlation between artemisinic alcohol and artemisinic acid (Figure 5A,B) is restricted to those lines that display increased artemisinic acid contents, and therefore is not statistically significant over all COSTREL lines analyzed (cf. Figure 5—source data 1). CC: correlation coefficient. Positive values indicate positive correlations and negative values indicate negative correlations. *: p<0.05; **: p<0.01; N: number of samples where both variables are present.

https://doi.org/10.7554/eLife.13664.015
Nt-AO3-CSAmorpha-4,11-dieneArtemisinic alcoholDihydroartemi-sinic alcoholArtemisinic aciddxrCYB5ADH1ALDH1DBR2
Amorpha-4,11-dieneCC1.0000.363**0.0790.2370.0060.0820.439**-0.0320.049
N605839596060606060
Art. alcoholCC0.363**1.0000.124-0.081-0.088-0.1190.0860.026-0.009
N585939595959595959
Dihydroart. alcoholCC0.0790.1241.0000.029-0.1690.1530.3060.1950.219
N393939393939393939
Art. acidCC0.237-0.0810.0291.0000.2410.1830.2170.386**0.141
N595939606060606060
dxrCC0.006-0.088-0.1690.2411.0000.1530.2240.415**0.170
N605939606161616161
CYB5CC0.082-0.1190.1530.1830.1531.0000.1850.322*0.503**
N605939606161616161
ADH1CC0.439**0.0860.3060.2170.2240.1851.0000.339**0.252*
N605939606161616161
ALDH1CC-0.0320.0260.1950.386**0.415**0.322*0.339**1.0000.444**
N605939606161616161
DBR2CC0.049-0.0090.2190.1410.1700.503**0.252*0.444**1.000
N605939606161616161
Table 3

List of oligonucleotides used in this study. The reverse primers (_R) for amplification of the genes FPS, ADS, CYP and CPR contain the sequence of the T7 promoter (bold) to facilitate in vitro transcription.

https://doi.org/10.7554/eLife.13664.018
GenePrimerSequence (5'→3')Purpose
dxrDXR_FCAACCTATGTACGTTGTTGGAGAAGAGGGqPCR
DXR_RCTGGAGCACCAGCAATCAATGTCTC
CYB5CYB5_FCCAGGAGGAGATGAAGTTCTTTTGGCTGqPCR
CYB5_RGCTGGAGGAACGTAAGCTCTCTTCTTTG
ADH1qADH1_F2TCCAGGTCATGAAGGTGTTGqPCR
qADH1_R2ATTGTCCACACTCACCAAGG
ALDH1ALDH1_FCCTGTTTCTTTGGAATTGGGTGGTAAGTCqPCR
ALDH1_RCAGCAACACACATCTCACCTTTGTTAGTG
DBR2DBR2_FGAGCAAGTTGAGGGTTGGAAGAAAGTTGqPCR
DBR2_RTAGAAGAGATAGGAGCAGCTCCACCTG
CYPCYP_qFCCTGAACCTTGGAGATTACCqPCR
CYP_qRGCCCATTTAGGAGAAGATACAAC
CPRCPR_qFCCTGTTGGAATGGGTGATGqPCR
CPR_qRCCTACAGCAGCAGTATAAGGAG
ACTINqTac9actin fCCTGAGGTCCTTTTCCAACCAqPCR
qTac9actin rGGATTCCGGCAGCTTCCATT
FPSFPS_probe_FCCTGCTTTTGAATTTGATGATGRNA probe
FPS_probe_RTAATACGACTCACTATAGGGCGAAACCAACAAGGTTGTCC
ADSADS_probe_FCTGAAGCTGTTGAAAGATGGTCRNA probe
ADS_probe_RTAATACGACTCACTATAGGGGGAGCAGATACAGCCCATTC
CYPCYP_probe_FTCCTCATCGAGGAGTACGAGRNA probe
CYP_probe_RTAATACGACTCACTATAGGGTTACAGGTCGTCCAGATCCAG
CPRCPR_probe_FATGATTGGTCCTGGAACTGGRNA probe
CPR_probe_RTAATACGACTCACTATAGGGGCCATTCCTTTAGCATCTCC

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Paulina Fuentes
  2. Fei Zhou
  3. Alexander Erban
  4. Daniel Karcher
  5. Joachim Kopka
  6. Ralph Bock
(2016)
A new synthetic biology approach allows transfer of an entire metabolic pathway from a medicinal plant to a biomass crop
eLife 5:e13664.
https://doi.org/10.7554/eLife.13664