- Research article
- Open Access
Relative expression of genes of terpene metabolism in different tissues of Artemisia annuaL
BMC Plant Biology volume 11, Article number: 45 (2011)
Recently, Artemisia annua L. (annual or sweet wormwood) has received increasing attention due to the fact that the plant produces the sesquiterpenoid endoperoxide artemisinin, which today is widely used for treatment of malaria. The plant produces relatively small amounts of artemisinin and a worldwide shortage of the drug has led to intense research in order to increase the yield of artemisinin. In order to improve our understanding of terpene metabolism in the plant and to evaluate the competition for precursors, which may influence the yield of artemisinin, we have used qPCR to estimate the expression of 14 genes of terpene metabolism in different tissues.
The four genes of the artemisinin biosynthetic pathway (amorpha-4,11-diene synthase, amorphadiene-12-hydroxylase, artemisinic aldehyde ∆11(13) reductase and aldehyde dehydrogenase 1) showed remarkably higher expression (between ~40- to ~500-fold) in flower buds and young leaves compared to other tissues (old leaves, stems, roots, hairy root cultures). Further, dihydroartemisinic aldehyde reductase showed a very high expression only in hairy root cultures. Germacrene A and caryophyllene synthase were mostly expressed in young leaves and flower buds while epi-cedrol synthase was highly expressed in old leaves. 3-Hydroxy-3-methyl-glutaryl coenzyme A reductase exhibited lower expression in old leaves compared to other tissues. Farnesyldiphosphate synthase, squalene synthase, and 1-deoxy-D-xylulose-5-phosphate reductoisomerase showed only modest variation in expression in the different tissues, while expression of 1-deoxy-D-xylulose-5-phosphate synthase was 7-8-fold higher in flower buds and young leaves compared to old leaves.
Four genes of artemisinin biosynthesis were highly expressed in flower buds and young leaves (tissues showing a high density of glandular trichomes). The expression of dihydroartemisinic aldehyde reductase has been suggested to have a negative effect on artemisinin production through reduction of dihydroartemisinic aldehyde to dihydroartemisinic alcohol. However, our results show that this enzyme is expressed only at low levels in tissues producing artemisinin and consequently its effect on artemisinin production may be limited. Finally, squalene synthase but not other sesquiterpene synthases appears to be a significant competitor for farnesyl diphosphate in artemisinin-producing tissues.
The genus Artemisia, belonging to the Asteraceae family, contains a large number of aromatic plants. During the latest decades, A. annua L. (annual or sweet wormwood) has received increasing attention due to the fact that the plant produces the sesquiterpenoid endoperoxide artemisinin, which today is widely used for the treatment of malaria . The plant produce relatively small amounts of artemisinin and a worldwide shortage of the drug has led to intense research in order to increase the yield of artemisinin in the plant or to develop alternative methods of artemisinin production . Artemisinin and a number of other terpenoids are produced in glandular secretory trichomes present on aerial surfaces of the plant. The glandular trichomes are 10-cell structures with three pairs of secretory cells . The apical cells are transparent proplastid- or leucoplast-containing cells, and the two cell pairs below the apical cells, the sub-apical cells, are green chloroplast-containing cells. The substances that are produced in these cells are excreted into, and stored in the subcuticular sac, which is covering the secretory cells . This has led to the focus on these structures when studying genes involved in terpene biosynthesis in general and for artemisinin biosynthesis in particular.
The major components of the essential oil from A. annua are mono- and sesquiterpenes produced in trichomes . The large number of structurally different terpenoids indicates that several mono- and sesquiterpene synthases are expressed in A. annua. In fact, so far five different sesquiterpene synthases (epi-cedrol synthase (ECS) , amorpha-4,11-diene synthase (ADS) ; β-caryophyllene synthase (CPS) , germacrene A synthase (GAS)  and β-farnesene synthase (BFS)  have been cloned from the plant. However, it is likely that additional sesquiterpene synthases are expressed in A. annua since other classes of sesquiterpenes such as eudesmanes (e.g. β-selinene), guaianes (e.g. α-guaiene), longipinanes (e.g. β-longipinene) and santalanes (e.g. α-santalol) have been isolated from the plant . All these sesquiterpene synthases may compete for the same substrate, farnesyl diphosphate (FDP), which is a precursor of artemisinin. Consequently, the activity of various sesquiterpene synthases may influence the yield of artemisinin in the plant. In addition, FDP is used for the synthesis of squalene, which is the precursor of sterols and triterpenes.
Figure 1 shows the biosynthetic pathway leading to artemisinin as it is understood today along with other pathways of terpene metabolism in A. annua. The first committed step in artemisinin biosynthesis is the conversion of FDP to amorpha-4,11-diene by ADS [6, 11]. In the following step, amorpha-4,11-diene is hydroxylated to yield artemisinic alcohol. This reaction is catalyzed by a cytochrome P450 dependent hydroxylase (CYP71AV1) . This enzyme can also oxidize the alcohol to artemisinic aldehyde and then further on to artemisinic acid . It has long been assumed that artemisinic acid is a direct precursor of artemisinin. However, recent feeding experiments with artemisinic acid  and dihydroartemisinic acid  have shown that the latter substance is the precursor of artemisinin. Dihydroartemisinic acid is formed from artemisinic aldehyde in two steps via dihydroartemisinic aldehyde. The reduction is catalyzed by artemisinic aldehyde Δ11(13) reductase (DBR2)  and the oxidation to the acid by aldehyde dehydrogenase 1 (ALDH1) . It has not been fully evaluated if the CYP71AV1 enzyme can catalyze the oxidation of dihydroartemisinic aldehyde to the corresponding acid [16, 17].
The conversion of dihydroartemisinic acid to artemisinin is believed to be a non-enzymatic spontaneous reaction . In a similar way, artemisinic acid is converted to arteannuin B as indicated in Figure 1.
Recently, a dihydroartemisinic aldehyde reductase (RED1) has been cloned from A. annua . This enzyme can potentially convert dihydroartemisinic aldehyde into dihydroartemisinic alcohol, a substance that appears to be a "dead end product", thereby affecting the yield of artemisinin in a negative way.
The genetic variation within A. annua appears to be high. At least two chemotypes with different compositions of the essential oil during the vegetative period have been described . One chemotype shows high content of dihydroartemisinic acid and artemisinin, while the second chemotype shows high content of artemisinic acid and arteannuin B but low amounts of artemisinin. According to previous investigations, there is no conversion in planta of artemisinic acid to dihydroartemisinic acid or vice versa . It has been suggested that arteannuin B can be converted to artemisinin in planta .
In order to increase our understanding of terpene metabolism in the plant A. annua and to evaluate the competition for precursors, which may influence the yield of artemisinin in the plant, we have used qPCR to estimate the expression of genes of terpene metabolism in different plant tissues. cDNA made from flower buds, young leaves, old leaves, stems, roots and hairy roots has been used as templates for the qPCR studies. In total, 14 different transcripts have been monitored using three reference genes. The different tissues have also been studied by fluorescence microscopy to determine trichome densities.
Results and Discussion
The glandular trichomes of A. annua have been extensively studied and it is well established that the production of specific terpenoids, including artemisinin, takes place within the secretory cells of such trichomes [21, 22]. Furthermore the yield of terpenoids in various plants is highly dependent on trichome abundance . In order to investigate the differences in trichome density on different plant tissues, glandular trichomes of A. annua were visualized by fluorescence microscopy (Figure 2) and found in all aerial tissues of the plant such as flower buds (panels A and B), leaves (panels C-F), stems (panels G and H) at different densities, but as expected not on roots (panels I and J) or hairy root cultures of A. annua (panels K and L).
From Figure 2, it is evident that young leaves (panels C and D) carry considerably higher density of trichomes than older leaves (panels E and F). All these tissues were taken from 5-6 months old plants with flower buds, which had been induced by a 21 days short day treatment (8 h day). We have estimated, by counting trichomes, that on young leaves there are around 300±40 trichomes/cm2 while the corresponding number for old leaves was 130 ± 25 trichomes/cm2. The original number of trichomes seem to be spread out on a larger surface as the leaf expands with age.
In order to examine the gene expression of 14 enzymes of terpene metabolism and three reference transcripts by qPCR, 18 primer pairs were synthesized as listed in Table 1. The efficiency of these primer pairs was calculated from the qPCR experiments and was found to be between 1.82 and 1.93 as summarized in Table 1.
Relative expression of genes of terpene metabolism in different tissues
The gene expression of the 14 enzymes of terpene metabolism was studied in different tissues of A. annua and the results are summarized in Table 2 and Figure 3. Transcripts of all the genes studied could be detected in all tissues at different amounts using qPCR. The only exception was that no ADS could be detected in roots. In Figure 3, expression levels in the different tissues are all compared separately to old leaves, which were arbitrary chosen as reference tissue. The results of the different experiments are presented below.
Relative expression of genes encoding HMGR and FDS in different tissues
The mevalonate pathway produces IDP and DMADP for production of different types of terpenoids as summarized in Figure 1. One enzyme of the mevalonate pathway (3-hydroxy-3-methyl-glutaryl-CoA reductase; HMGR) has been studied by qPCR in this paper. Our results show that HMGR was expressed at higher levels in biosynthetically active tissues such as flower buds, young leaves, roots and hairy roots. The CT-values for this transcript in these tissues were relatively low (23.8-24.7) indicating that HMGR was expressed at relatively high levels (Table 2). HMGR is a rate-limiting enzyme of the mevalonate pathway [24, 25], and recently it has been shown that HMGR expression limits artemisinin formation in A. annua . The regulation of the mevalonate pathway by HMGR was reflected in the result that around 20- and 8-fold more expression was seen in highly biosynthetically active flower buds and young leaves, respectively, than in old leaves (Figure 3). Stermer et al. showed that the HMGR activity is much reduced in mature tissues, correlating well with the results presented here . It is interesting to note that HMGR also was enhanced 20-fold in roots. This might be due to the relatively rapid growth of roots and the requirement for sterols and protective terpenoids .
Farnesyl diphosphate synthase (FDS) plays a central role in the conversion of IDP and DMADP, produced by the mevalonate pathway, to various terpenoids such as sterols and sesquiterpenes. A small gene family encodes isoenzymes of FDS in plants. FDS1 is a housekeeping gene as it is expressed in all tissues at all developmental stages . It is involved in the synthesis of isoprenoids serving basic plant cell functions such as production of sterols for various membranes. Overexpression of FDS1 in A. annua results in an increase in artemisinin production [29, 30] indicating that it is a rate-limiting enzyme in this pathway. FDS2 is often inducible and found in special tissues at particular stages of development  and is also involved in phytoalexin biosynthesis. FDS1 [31, 32] and FDS2 (GenBank accession number AF136602) have been cloned from A. annua. The two FDS genes from A. annua show 96.6% nucleotide sequence identity and the segment that was amplified with the primers used (Table 1) is identical in the two genes. Consequently, the expression of both genes was monitored in our qPCR experiments. From Figure 3 it is evident that the total expression of FDS1/FDS2 is essentially the same in all tissues, and the relatively small variation observed may be explained by the fact that the housekeeping gene FDS1 constitutes a major part of total FDS1/FDS2 and that the eventual different level of FDS2 expression in different tissues thereby becomes masked.
Relative expression of genes of the MEP pathway in different tissues
The MEP pathway produces the precursors for many different terpenoids as shown in Figure 1. Three of the eight enzymes of the MEP pathway, 1-deoxy-D-xylulose 5-phosphate synthase (DXS), 1-deoxy-D-xylulose 5-phosphate reductase (DXR) and hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase (HDR), have been studied by qPCR in this paper. Our results show that the expression of DXS was 7-8-fold higher in flower buds and young leaves compared to old leaves (Figure 3). It is well established that DXS, the first enzyme of the MEP pathway, plays a major role in the overall regulation of the pathway . Higher expression of DXS leads to an enhanced production of IDP and DMADP, which may influence the synthesis of terpenoids in the plastid. However, an enhanced DXS activity may also influence the production of artemisinin since it has recently been shown that one C5-unit of FDP, the precursor of artemisinin, is produced by the MEP pathway . In the other tissues, the expression of DXS was equal to that in old leaves. DXR and HDR may also have rate-limiting roles for the production of IDP and DMADP . The expression of DXR appears to be somewhat enhanced in flower buds and hairy roots but lower in stems and roots. The expression of HDR was 10-30-fold higher in old leaves than in the other tissues.
Relative expression of genes of artemisinin biosynthesis in different tissues
During the last decade, a number of enzymes of artemisinin biosynthesis have been cloned and a putative pathway has been constructed as outlined in Figure 1. We have monitored the expression of four genes (ADS, CYP71AV1, DBR2 and ALDH1) involved in the conversion of FDP to dihydroartemisinic acid, which is a late precursor of artemisinin . We have also included RED1 in these studies as it has been suggested that this enzyme may have an influence on the yield of artemisinin by withdrawing dihydroartemisinic aldehyde from further conversion to artemisinin .
The first committed enzyme of artemisinin biosynthesis is ADS , which has been cloned by several investigators [6, 35–37]. The next enzyme, CYP71AV1, was cloned by two groups independently [12, 38]. The following two enzymes, DBR2 and ALDH1, were cloned and characterized by Zhang et al.  and Teoh et al. , respectively. Finally, RED1 was cloned and characterized by Rydén et al. .
In Table 3 we summarize previous studies of gene expression of enzymes involved in artemisinin biosynthesis in A. annua. The specific expression of CYP71AV1, DBR2 and ALDH1 in glandular trichomes of A. annua has been shown by PCR, and ADS has been shown to be expressed in glandular trichomes by immuno gold staining in combination with silver enhancement . The relative expression of these different genes cannot be compared since the data in Table 3 is from several different studies and on different varieties of A. annua. It may, however, be concluded from Table 3 that the biosynthetic enzymes are expressed in flower buds, leaves and stems that carry glandular trichomes but not in roots.
In A. annua, DBR2 belongs to a small family of enolate reductases and a second enzyme of this family is 12-oxophytodienoate reductase (OPR3), which in Arabidopsis is involved in jasmonic acid biosynthesis . The OPR3-like gene of A. annua (Genbank accession number EU848577) shows a very high similarity to DBR2, and in our initial attempt to determine the relative expression of DBR2, a primer pair (see Table 1) with identical nucleotide sequence in the DBR2 and OPR3-like sequences was used. Consequently, the two transcripts were amplified simultaneously. The combined expression of these two genes is essentially the same in flower buds, young leaves and old leaves (Figure 3; panels A and B). However, a 21 bp deletion in the DBR2 sequence was thereafter used to design, by spanning this deletion, a reversed primer that is specific for DBR2. Using this primer, the expression of DBR2 could be determined as discussed below.
The qPCR data presented here clearly demonstrate that the expression levels of ADS, CYP71AV1, DBR2 and ALDH1 were very high in flower buds and young leaves compared to old leaves (Figure 3; panel A and B), i.e. ~500 and ~150 times more ADS, ~150 and ~150 times more CYP71AV1, ~40 and ~70 times more DBR2 and ~250 and ~130 times more ALDH1 in flower buds and young leaves, respectively. These high expression levels indicate a high capacity to produce artemisinin precursors in these tissues carrying biosynthetically active trichomes. These findings are supported by a study on the induction of enzymes involved in artemisinin biosynthesis by jasmonic acid . It was shown by qPCR that the levels of ADS, CYP71AV1, DBR2 and ALDH1 were around 50, 13, 160 and 55 times higher in isolated trichomes (from flower buds) than in leaves, respectively , which once again shows that the genes encoding enzymes of artemisinin biosynthesis are specifically expressed in trichomes.
Even though a 70-150-fold difference in expression of the four genes was observed between young and old leaves, the trichome density on these tissues only differs around three times. Obviously, the relative amount of the biosynthetic enzymes is much higher in young trichomes, and the regulation of expression of these enzymes appears to be linked to the developmental stage of the trichomes.
The qPCR data may be used to estimate the relative amount of transcripts using the 2-ΔΔCTmethod  and the β-actin as reference gene. As an example, we have made such an estimation for flower buds. For this crude estimation, we assume that the number of active sites is proportional to the level of transcription and that the enzymes are working at substrate saturation with an optimal NADPH/NADP+ ratio. Under such conditions the kcat-value is a good indicator of the conversion of substrate to product. Consequently, we may calculate the relative turnover potential for the different enzymes. The values were normalized to RED1 (=1) due to its low abundance in flower buds. Using published kcat-values for the different enzymes, we may calculate their relative turnover potential in planta as summarized in Table 4. For CYP71AV1 no kcat-value is available and therefore this enzyme is not included in this estimation. It may be concluded that ADS is a rate-limiting enzyme due to its low kcat-value and that the other trichome-specific enzymes DBR2 and ALDH1 (Figure 3) are present in excess assuming that they are acting on the same pool of intermediates. Due to differences in the kcat-value, ALDH1 exhibits a 5-fold higher potential conversion of dihydroartemisinic aldehyde as compared to artemisinic aldehyde (Table 4).
From Figure 3, we may conclude that only part of the FDS1/FDS2 activity is present in trichomes and involved in artemisinin biosynthesis. In fact, it may be speculated that FDS1 is a critical enzyme for the yield of artemisinin. If we assume that only one percent of the FDS1 activity is found in trichome cells, the relative turnover potential is in the same range as ADS. This assumption is supported by the fact that upregulation of FDS by genetic transformation of A. annua using constitutive promoter resulted in an increased artemisinin production  indicating that FDS1 also may be a rate-limiting enzyme in artemisinin biosynthesis. However, due to the relative high potential conversion capacity of DBR2 and ALDH1 it may be assumed that as amorpha-4,11-diene is formed by ADS, it is efficiently converted all the way to dihydroartemisinic acid, which accumulates in this chemotype of A. annua. This is supported by the fact that essentially no amorpha-4,11-diene is detected by GC-MS of an hexane extract from this variety of A. annua (unpublished). Finally, since dihydroartemisinic acid accumulates in this chemotype, we may conclude that the activity of DBR2 is significantly higher than that of CYP71AV1. The fact that DBR2 is a soluble cytosolic protein while CYP71AV1 is embedded in the ER may influence the metabolic flow and the rate of intracellular transport of intermediates has to be considered. At present, we do not have any information on the transport of intermediates of artemisinin biosynthesis within cells. It would be interesting to study the ratio of DBR2 and CYP71AV1 in the two different chemotypes of A. annua since these enzymes are the key enzymes for the formation of dihydroartemisinic acid and artemisinic acid, respectively.
Next it may be concluded that RED1 does not appear to have any significant influence on the biosynthesis of artemisinin, which was suggested by Rydén et al.  due to a high Km for dihydroartemisinic aldehyde, low relative turnover potential and only partial localization to trichomes. Expression of RED1 in flower buds, young leaves and stems was relatively low and just a fraction of that observed in old leaves (Figure 3A-3C). Furthermore, it is interesting to note that a high expression of RED1 was observed in hairy roots of A. annua, i.e. around 50 times more than in old leaves (Figure 3E). The function of RED1 in hairy roots is not known but in this case the high RED1 activity together with an extremely low ALDH1 expression may influence the yield of artemisinin in hairy roots, which is relatively low [43, 44].
Finally, no ADS could be detected in roots of A. annua while small amounts were found in hairy roots. Hairy roots produce as mentioned above small amounts of artemisinin. However, ALDH1 was almost not detectable in hairy roots, which may lead to the formation of artemisinic acid and arteannuin B instead of artemisinin. These compounds have been isolated from hairy root cultures of A. annua .
Relative expression of genes of sesquiterpene biosynthesis in different tissues
In addition to ADS, four other sesquiterpene synthases have been cloned from A. annua. These enzymes with very similar kinetic properties compete for the FDP available and may therefore influence the production of artemisinin precursors. Three of these sesquiterpene synthases, i.e. CPS, ECS and GAS have been monitored by qPCR in this study.
The expression of ECS was much lower in flower buds and young leaves than in old leaves (Figure 3A and 3B), while expression of CPS and GAS was considerably higher in these tissues as compared to old leaves. In old leaves, ECS appears to be the only sesquiterpene synthase highly expressed. The function of ECS in A. annua has not yet been established.
The relative amounts of the sesquiterpene synthases has been calculated using the 2-ΔΔCTmethod  for some tissues as summarized in Table 5 using β-actin as reference. The values were normalized to GAS expression in stems (=1) due to its low abundance there. As can be seen in Table 5, ADS is the dominating sesquiterpene synthase in flower buds and young leaves, i.e. the tissues where the biosynthesis of artemisinin precursors take place. It may be concluded that ADS is capable to channel a major part of the available substrate into artemisinin biosynthesis. Ma et al  showed that upregulation of ADS by genetic transformation of A. annua resulted in an increased formation of artemisinin showing that ADS is a rate-limiting enzyme .
Relative expression of genes of sterol and triterpene biosynthesis in different tissues
A key enzyme in the biosynthesis of sterols and triterpenes is squalene synthase (SQS), which condensates two molecules of FDP to squalene. The expression of SQS varied only modest between different tissues in our qPCR study. SQS is a microsomal enzyme and difficult to purify and characterize. Consequently, there is no reported kcat-value for any plant enzyme, but a kcat = 3.3 sec-1 has been reported for a soluble recombinant yeast enzyme . The A. annua enzyme has been cloned and is similar to the enzymes from A. thaliana and Nicotiana tabacum . As seen in Table 5, a significant amount of SQS was expressed in the tissues with high artemisinin biosynthesis. Assuming a kcat-value for the A. annua SQS in the order of 1 sec-1 (= 2 FDPs used per second), results in about 45- and 25-fold higher turnover potential, of FDP to products, for SQS, as compared to the turnover potential of ADS in flower buds and young leaves, respectively. The fraction of SQS expressed in trichome cells may channel part of available FDP into sterols and thereby lower the yield of artemisinin. Zhang et al. showed that down-regulation of SQS by hairpin-RNA-mediated gene silencing in A. annua resulted in a 3-fold increased artemisinin production .
The aim of our studies was to increase our understanding of terpene metabolism in the plant A. annua. The qPCR data presented in this paper demonstrate that four genes of the artemisinin biosynthetic pathway (ADS, CYP71AV1, DBR2 and ALDH1) showed remarkably higher expression (between 40- to 500-fold) in flower buds and young leaves compared to other tissues (old leaves, stems, roots and hairy root cultures) (Figure 3). These high expression levels indicate a much higher capacity to produce artemisinin precursors in flower buds and young leaves, which is partly due to the considerably higher density of trichomes on these tissues (Figure 2).
Our aim was also to evaluate the competition for precursors, which may influence the yield of artemisinin in the plant. The expression of other sesquiterpene synthases was much lower than ADS in tissues producing artemisinin precursors (i.e. flower buds and young leaves). Consequently, their influence on artemisinin yield appear to be relatively limited and downregulation of other sesquiterpene synthase(s) will most likely not affect artemisinin production in A. annua.
However, squalene synthesis may influence the yield of artemisinin, since the potential utilization of FDP by SQS is most likely higher than that of ADS. The fact that ADS is a cytosolic enzyme and SQS is localized to the ER may be favorable for the sesquiterpene synthase. Further, the ratio of CYP71AV1 and DBR2 may be critical for the type of end product formed; high DBR2 activity will result in an efficient formation of dihydroartemisinic aldehyde, dihydroartemisinic acid and artemisinin while high CYP71AV1 activity may result in artemisinic acid and arteannuin B. A direct comparison of the capacity of these two enzymes is not possible due to the fact that no kinetic data is available for CYP71AV1. However in this chemotype of A. annua, the turnover capacity of DBR2 is apparently sufficient to convert most of the artemisinic aldehyde to dihydroartemisinic aldehyde for further conversion to artemisinin. It is notable that the RED1 activity was very low in artemisinin-producing tissues and consequently this enzyme appears not to influence the yield of the artemisinin precursor dihydroartemisinic acid to any significant extent.
Seeds of Artemisia annua L. cv Artemis were obtained from Anamed, Germany (http://www.anamed.net). This variety is a high dihydroartemisinic acid chemotype (the ratio artemisinic acid to dihydroartemisinic acid is 0.36 ). Plants were grown under 16 h days and 8 h nights at 22°C to a height of approximately 1 m followed by flower bud induction at 8 h days and 16 h nights at 22°C.
Artemisia annua L. hairy root cultures, kindly provided by Dr. Kanyaratt Supaibulwatana, Mahidol University, Bangkok, were grown in MSMO medium (Sigma-Adrich, Stockholm, Sweden) pH 5.8, supplemented with 3% (w/v) sucrose at room temperature on a rotary shaker at 75 rpm.
Flower buds, young leaves, old leaves, stems, roots and hairy roots were collected separately from 5-6 months old plants after induction of flower buds, frozen in liquid nitrogen, ground to a fine powder in a mortar and used for RNA extraction.
RNA extraction was performed using Purelink™ Plant RNA Reagent kit (Invitrogen, Carlsbad, California, USA) for small scale RNA isolation according to the manufacturer's instructions. Frozen plant tissue powder (100 mg) was used for each RNA extraction. The RNA was DNase treated using DNase I (Fermentas, St Leo-Roth, Germany) to remove remaining genomic DNA.
First strand cDNA synthesis
RNA (1 μg) was reverse transcribed using RevertAid™ H Minus-MuLV reverse transcriptase (Fermentas) primed with 0.5 μg oligo(dT)18 primer. The RNA was removed from the first strand cDNA by RNase treatment using RNase H (Fermentas) according to the manufacturer's instructions.
Quantitative PCR (qPCR)
The qPCR was performed using specific primers (Table 1) on a 7500 qPCR equipment (Applied Biosystems, Foster City, USA). First strand cDNA was used as template in 20 μl reactions including 10 μl Power SYBR®Green PCR Master Mix (Applied Biosystems) and 2 pmol of each primer. The qPCR cycling was performed at 50 °C (2 min), 95 °C (10 min), 40 cycles at 95 °C (15 s), 60 °C (1 min) and finally a dissociation stage at 95 °C (15 s), 60 °C (1 min), 95 °C (15 s). The dissociation stage was performed to determine the PCR product size and to detect possible primer dimers. Triplets of all samples were run, and a negative control of the Master Mix in addition of primers was performed in all qPCR runs.
Any outliers in the triplets were detected by Grubbs test and removed. The mean efficiency of the amplicons was calculated by the program Linreg v. 12.1 based on the log linear phase of the amplification curve . Baseline corrected data was imported set between cycles 3-15, from the 7500 qPCR software (Applied Biosystems) to the Linreg software. Linreg calculated window of linearity for each amplicon and the efficiency from the curve fit showing highest correlation value between 4-6 points within the window of linearity. The cycle threshold (CT) values and efficiency values given by the Linreg software were used for further analysis. The BestKeeper software  was used to search for stable reference genes among all genes tested. Based on BestKeeper, the genes CPR, β-actin and PAL were selected as reference genes. Relative expression levels were calculated using the REST 2009 software V. 2.0.13 (Qiagen, Hilden, Germany) .
Flower buds, young leaves, old leaves, stems and roots were selected from the same A. annua plant and microscopy studies were performed to compare the amount of trichomes on different plant tissue. All micrographs were captured using a Nikon e400 C-SHG1 fluorescence microscope equipped with digital camera, using light microscopy and filter settings for FITC (λex 480 nm; λem = 535 nm) and Texas Red (λex 570 nm; λem = 625 nm, images not shown). Autofluorescence was separated in the red and green channels and brightness adjusted using NIS-elements imaging software v. 2.20 (Nikon, Badhoevedorp, The Netherlands). All tissue images have been captured using the same magnification (4x objective; 10x on ocular).
aldehyde dehydrogenase 1
cytochrome P450 reductase
artemisinic aldehyde Δ11(13) reductase
1-deoxy-D-xylulose 5-phosphate reductase
1-deoxy-D-xylulose 5-phosphate synthase
farnesyl diphosphate synthase
germacrene A synthase
hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase
Murashige and Skoog medium with minimal organics
phenylalanine ammonia lyase
dihydroartemisinic aldehyde reductase
Bosman A, Mendis KN: A major transition in malaria treatment: The adoption and deployment of artemisinin-based combination therapies. Am J Trop Med Hyg. 2007, 77: 193-197.
Covello PS: Making artemisinin. Phytochemistry. 2008, 69: 2881-2885. 10.1016/j.phytochem.2008.10.001.
Duke SO, Paul RN: Development and fine structure of the glandular trichomes of Artemisia annua L. Int J Plant Sci. 1993, 154: 107-118. 10.1086/297096.
Simon JE, Charles D, Cebert E, Grant L, Janick J, Whipkey A: Artemisia annua L.: a promising aromatic and medicinal. Advances in new crops. Edited by: Janick J and Simon JE. Portland: Timber Press; 1990, 522-526.
Mercke P, Crock J, Croteau R, Brodelius PE: Cloning, expression and characterization of epi-cedrol synthase, a sesquiterpene cyclase from Artemisia annua L. Arch Biochem Biophys. 1999, 369: 213-222. 10.1006/abbi.1999.1358.
Mercke P, Bengtsson M, Bouwmeester HJ, Posthumus MA, Brodelius PE: Molecular cloning, expression, and characterization of amorpha-4,11-diene synthase, a key enzyme of artemisinin biosynthesis in Artemisia annua L. Arch Biochem Biophys. 2000, 381: 173-180. 10.1006/abbi.2000.1962.
Cai Y, Jia J-W, Crock J, Lin Z-X, Chen X-Y, Croteau R: A cDNA clone for β-caryophyllene synthase from Artemisia annua. Phytochemistry. 2002, 61: 523-529. 10.1016/S0031-9422(02)00265-0.
Bertea CM, Voster A, Verstappen FW, Maffei M, Beekwilder J, Bouwmeester HJ: Isoprenoid biosynthesis in Artemisia annua: cloning and heterologous expression of a germacrene A synthase from a glandular trichome cDNA library. Arch Biochem Biophys. 2006, 448: 3-12. 10.1016/j.abb.2006.02.026.
Picaud S, Brodelius M, Brodelius PE: Expression, purification and characterization of β-farnesene synthase from Artemisia annua L. Phytochemistry. 2005, 66: 961-967. 10.1016/j.phytochem.2005.03.027.
Brown GD: The biosynthesis of artemisinin and the phytochemistry of Artemisia annua L. Molecules. 2010, 15: 7603-7698. 10.3390/molecules15117603.
Bouwmeester HJ, Wallaart TE, Janssen MH, van Loo B, Jansen BJ, Posthumus MA, Schmidt CO, de Kraker JW, Konig WA, Franssen MC: Amorpha-4, 11-diene synthase catalyses the first probable step in artemisinin biosynthesis. Phytochemistry. 1999, 52: 843-854. 10.1016/S0031-9422(99)00206-X.
Teoh KH, Polichuk DR, Reed DW, Nowak G, Covello PS: Artemisia annua L. (Asteraceae) trichome-specific cDNAs reveal CYP71AV1, a cytochrome P450 with a key role in the biosynthesis of the antimalarial sesquiterpene lactone artemisinin. FEBS Lett. 2006, 580: 1411-1416. 10.1016/j.febslet.2006.01.065.
Brown GD, Sy LK: In vivo transformations of artemisinic acid in Artemisia annua plants. Tetrahedron. 2006, 63: 9548-9566. 10.1016/j.tet.2007.06.062.
Brown GD, Sy LK: In vivo transformations of dihydroartemisinic acid in Artemisia annua plants. Tetrahedron. 2004, 60: 1139-1159. 10.1016/j.tet.2003.11.070.
Zhang Y, Teoh KH, Reed DW, Maes L, Goossens A, Olson DJ, Ross AR, Covello PS: The molecular cloning of artemisinic aldehyde Δ11(13) reductase and its role in glandular trichome-dependent biosynthesis of artemisinin in Artemisia annua. J Biol Chem. 2008, 283: 21501-21508. 10.1074/jbc.M803090200.
Teoh KH, Polichuk DR, Reed DW, Covello PS: Molecular cloning of an aldehyde dehydrogenase implicated in artemisinin biosynthesis in Artemsia annua. Botany. 2009, 87: 635-642. 10.1139/B09-032.
Rydén A-M: Identification, characterization and expression of early biosynthetic genes from Artemisia annua. PhD Thesis. University of Groningen, Department of Pharmaceutical Biology; 2010.
Rydén A-M, Ruyter-Spira C, Quax WJ, Osada H, Muranaka T, Kayser O, Bouwmeester H: The Molecular cloning of dihydroartemisinic aldehyde reductase and its implication in artemisinin biosynthesis in Artemisia annua. Planta Med. 2010, 76 (15): 1778-1783.
Wallaart TE, Pras N, Beekman AC, Quax WJ: Seasonal variation of artemisinin and its biosynthetic precursors in plants of Artemisia annua of different geographical origin; proof for the existence of chemotypes. Planta Med. 2000, 66: 57-62. 10.1055/s-2000-11115.
Dhingra V, Narasu ML: Purification and characterization of an enzyme involved in biochemical transformation of arteannuin B to artemisinin from Artemisia annua. Biochem Biophys Res Commun. 2001, 281: 558-561. 10.1006/bbrc.2000.4197.
Ferreira JFS, Janick J: Distribution of artemisinin in Artemisia annua. Progress in new crops. Edited by: Janick J. Arlington: ASHS Press; 1996,579-584.
Tellez MR, Canel C, Rimando AM, Duke SO: Differential accumulation of isoprenoids in glanded and glandless Artemisia annua L. Phytochemistry. 1999, 52: 1035-1040. 10.1016/S0031-9422(99)00308-8.
Biswas KK, Foster AJ, Aung T, Mahmoud SS: Essential oil production: relationship with abundance of glandular trichomes in aerial surface of plants. Acta Physiol Plant. 2009, 31: 13-19. 10.1007/s11738-008-0214-y.
Chappell J, Wolf F, Proulx J, Cuellar R, Saunders C: Is the reaction catalyzed by 3-hydroxy-3-methylglutaryl coenzyme A reductase a rate-limiting step for isoprenoid biosynthesis in plants?. Plant Physiol. 1995, 109: 1337-1343.
Schaller H, Grausem B, Benveniste P, Chye M-L, Tan CT, Song Y-H, Chua N-H: Expression of the Hevea brasiliensis (H.B.K.) müll. Arg. 3-hydroxy-3-methylglutaryl coenzyme A reductase 1 in tobacco results in sterol overproduction. Plant Physiol. 1995, 109: 761-770.
Ram M, Khan MA, Jha P, Khan S, Kiran U, Ahmad MM, Javed S, Abdin MZ: HMG-CoA reductase limits artemisinin biosynthesis and accumulation in Artemisia annua L. plants. Acta Physiol Plant. 2010, 32: 859-866. 10.1007/s11738-010-0470-5.
Stermer BA, Bianchini GM, Korth KL: Regulation of HMG-CoA reductase activity in plants. J Lipid Res. 1994, 35: 1133-1140.
Cunillera N, Boronat A, Ferrer A: Spatial and temporal patterns of GUS expression directed by 5' regions of the Arabidopsis thaliana farnesyl diphosphate synthase genes FPS1 and FPS2. Plant Mol Biol. 2000, 44: 747-758. 10.1023/A:1026588708849.
Han JL, Liu BY, Ye HC, Wang H, Li ZQ, Li GF: Effects of overexpression of the endogenouse farnesyl diphosphate synthase on the artemisinin content in Artemisia annua L. J Integr Plant Biol. 2006, 48: 482-487. 10.1111/j.1744-7909.2006.00208.x.
Banyai W, Kirdmanee C, Mii M, Supaibulwatana K: Overexpression of farnesyl pyrophosphate synthase (FPS) gene affected artemisinin content and growth of Artemisia annuaL. Plant Cell Tiss Organ Cult.
Matsushita Y, Kang WK, Charlwood BV: Cloning and analysis of a cDNA encoding farnesyl diphosphate synthase from Artemisia annua. Gene. 1996, 172: 207-209. 10.1016/0378-1119(96)00054-6.
Zhao Y, Ye H, Li G, Chen D, Liu Y: Cloning and enzymology analysis of farnesyl pyrophosphate synthase gene from a superior strain of Artemisia annua L. Chi Sci Bull. 2003, 48: 63-67.
Cordoba E, Salmi M, León P: Unravelling the regulatory mechanisms that modulate the MEP pathway in higher plants. J Exp Bot. 2009, 60: 2933-2943. 10.1093/jxb/erp190.
Schramek N, Wang H, Römisch-Margl W, Keil B, Radykewicz T, Winzenhörlein B, Beerhues L, Bacher A, Rohdich F, Gershenzon J, Liu B, Eisenreich W: Artemisinin biosynthesis in growing plants of Artemisia annua. A 13CO2 study. Phytochemistry. 2010, 71: 179-187. 10.1016/j.phytochem.2009.10.015.
Chang Y-J, Song S-H, Park S-H, Kim S-U: Amorpha-4,11-diene synthase of Artemisia annua: cDNA isolation and bacterial expression of a terpene synthase involved in artemisinin biosynthesis. Arch Biochem Biophys. 2000, 383: 178-184. 10.1006/abbi.2000.2061.
Wallaart TE, Bouwmeester HJ, Hille J, Poppinga L, Maijers NC: Amorpha-4,11-diene synthase: cloning and functional expression of a key enzyme in the biosynthetic pathway of the novel antimalarial drug artemisinin. Planta. 2001, 212: 460-465. 10.1007/s004250000428.
Li ZQ, Liu Y, Liu BY, Wang H, Ye HC, Li GF: Cloning, E. coli expression and molecular analysis of amorpha-4,11-diene synthase from a high-yield strain of Artemisia annua L. J Integ Plant Biol. 2006, 48: 1486-1492. 10.1111/j.1744-7909.2006.00381.x.
Ro DK, Paradise EM, Ouellet M, Fisher KJ, Newman KL, Ndungu JM, Ho KA, Eachus RA, Ham TS, Kirby J, Chang MC, Withers ST, Shiba Y, Sarpong R, Keasling JD: Production of the antimalarial drug precursor artemisinic acid in engineered yeast. Nature. 2006, 440: 940-943. 10.1038/nature04640.
Olsson ME, Olofsson LM, Lindahl A-L, Lundgren A, Brodelius M, Brodelius PE: Localization of enzymes of artemisinin biosynthesis to the apical cells of glandular secretory trichomes of Artemisia annua L. Phytochemistry. 2009, 68: 1864-1871.
Schaller F, Biesgen C, Müssig C, Altmann T, Weiler EW: 12-Oxophytodienoate reductase 3 (OPR3) is the isoenzyme involved in jasmonate biosynthesis. Planta. 2000, 210: 979-984. 10.1007/s004250050706.
Maes L, Van Nieuwerburgh FCW, Zhang Y, Reed DW, Pollier J, Vande Casteele SRF, Inzé D, Covello PS, Deforce DLD LD, Goossens A: Dissection of the phytohormonal regulation of trichome formation and biosynthesis of the antimalarial compound artemisinin in Artemisia annua plants. New Phytologist. 2011, 189: 176-189. 10.1111/j.1469-8137.2010.03466.x.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real- time quantitative PCR and the 2-ΔΔCTmethod. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.
Liu CZ, Wang YC, Ouyang F, Ye HC, Li GF: Production of artemisinin by hairy root cultures of Artemisia annua L. Biotechnol Lett. 1997, 19: 927-929. 10.1023/A:1018362309677.
Weathers PJ, Bunk G, McCoy MC: The effect of phytohormones on growth and artemisinin production in Artemisia annua hairy root. In Vitro Cell Dev Biol - Plant. 2005, 41: 47-53. 10.1079/IVP2004604.
Banerjee S, Zehra M, Gupta MM, Kumar S: Agrobacterium rhizogenes-mediated transformation of Artemisia annua: production of transgenic plants. Planta Med. 1997, 63: 467-469. 10.1055/s-2006-957737.
Ma C, Wang H, Lu X, Wang H, Xu G, Liu B: Terpenoid metabolic profiling analysis of transgenic Artemisia annua L. by comprehensive two-dimensional gas chromatography time-of-flight mass spectrometry. Metabolomics. 2009, 5: 497-506. 10.1007/s11306-009-0170-6.
Zhang D, Jennings SM, Robinson GW, Poulter CD: Yeast squalene synthase: Expression, purification and characterization of soluble recombinant enzyme. Arch Biochem Biophys. 1993, 304: 133-143. 10.1006/abbi.1993.1331.
Liu CZ, Ye HC, Wang H, Li GF: Molecular cloning, Escherichia coli expression and genomic organization of squalene synthase gene from Artemisia annua. Acta Bot Sin. 2003, 45: 608-613.
Zhang L, Jing F, Li F, Li M, Wang Y, Wang G, Sun X, Tang K: Development of transgenic Artemisia annua (Chinese wormwood) plants with an enhanced content of artemisinin, an effective anti-malarial drug, by hairpin-RNA-mediated gene silencing. Biotechnol Appl Biochem. 2009, 52: 199-207. 10.1042/BA20080068.
Ferreira JFS, Gonzalez JM: Analysis of underivatized artemisinin and related sesquiterpene lactones by high-performance liquid chromatography with ultraviolet detection. Phytochem Anal. 2009, 20: 91-97. 10.1002/pca.1101.
Ruijter JM, Ramakers C, Hoogaars WMH, Karlen Y, Bakker O, van den Hoff MJB, Moorman AFM: Amplification efficiancy, linking baseline and bias in the analysis of quantitative PCR data. Nucl Acids Res. 2009, 37: e45-10.1093/nar/gkp045.
Pfaffl MW, Tichopád A, Prgomet C, Neuvians TP: Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper - Excel-based tool using pair-wise correlations. Biotechnol Lett. 2004, 26: 509-515. 10.1023/B:BILE.0000019559.84305.47.
Pfaffl MW, Horgan GW, Dempfle L: Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real-time PCR. Nucl Acids Res. 2002, 30: e36-10.1093/nar/30.9.e36.
Picaud S, Olofsson L, Brodelius M, Brodelius PE: Expression, purification and characterization of amorpha-4,11-diene synthase from Artemisia annua L. Arch Biochem Biophys. 2005, 436 (2): 215-226.
We are grateful to the reviewers of this paper who made us aware of the high similarity between DBR2 and OPR3. We are grateful to Dr. Kanyaratt Supaibulwatana, Mahidol University, Bangkok who supplied the hairy root culture of Artemisia annua L. The financial support by the Faculty of Science and Engineering is acknowledged.
LO planned the experimental setup and PB was involved in the study design. AL prepared the plants. LO and AE carried out the RNA extraction, cDNA synthesis and qPCR analysis and analyzed the data. LO and AL performed the fluorescence microscopy. PB, LO, AL drafted and wrote the manuscript. All authors have read and approved the final manuscript.
About this article
Cite this article
Olofsson, L., Engström, A., Lundgren, A. et al. Relative expression of genes of terpene metabolism in different tissues of Artemisia annuaL. BMC Plant Biol 11, 45 (2011). https://doi.org/10.1186/1471-2229-11-45
- Hairy Root
- Hairy Root Culture
- Glandular Trichome
- Artemisinic Acid