Introduction
Secondary metabolites in plants are the products of long-term evolution in a particular environment; these compounds play important roles in many vital activities [
1,
2]. The development of an ecotype within a habitat could significantly affect the accumulation of secondary metabolites. Heavy metal pollution in soil and environmental stresses such as drought, pest, and disease damage, are vital factors that affect secondary metabolite accumulation. Pollution of cadmium (Cd), one of the most poisonous heavy metals, causes physiological and biochemical damages to plants; these damages include inhibited seed germination and plant growth [
3,
4], decreased photosynthetic efficiency [
5,
6], and varied uptake rates of soil ions [
7]. In addition, Cd can alter the biosynthesis of secondary metabolites and their related gene expression. Vitti
et al. showed that the amount of IAA in
Arabidopsis thaliana roots treated with Cd was significantly higher than that in the control group; the expression of
AtIPT in the shoot and
AtCKX in the root was significantly up-regulated after Cd treatment [
8]. Daud
et al. showed that the malondialdehyde (MDA) concentrations in the leaves, stems, and roots of two transgenic cotton cultivars treated with 1000 μmol/L Cd were higher than those in their respective controls, except in the roots of BR001; the highest MDA concentrations were found in the roots of GK30 [
9]. Thus far, the effects of Cd on medicinal plants, particularly on the biosynthesis of active substances, have been rarely investigated [
10,
11].
The above-ground parts of
Artemisia annua, a therophyte from the Asteraceae family, is a traditional Chinese medicine, known as
qinghao, used to clear vacuity heat, cool the blood, and prevent malaria. Modern pharmacological research shows that artemisinin (AN) is the main ingredient for malaria prevention. In 2001, the World Health Organization recommended artemisinin combination therapies (ACT) as the first-line therapy for malaria worldwide. Since its discovery, commercial preparations of artemisinin have been mainly isolated from plants. Studies show that wounding [
12], water [
13], chitosan [
14], sugars [
15], salicylic acid [
16], and MeJA [
17] enhanced the biosynthesis of artemisinin metabolites by upregulating the expression of key enzyme genes in the biosynthesis pathways. However, gene expression and increased rate differ when different elicitors are used. For example, MeJA upregulated
CYP71AV1, and miconazole enhanced
CPR and
DBR2 expression [
17].
AN is mainly produced and stored in the glandular secretory trichomes on leaves, stems, and inflorescences [
18] of
A. annua. Fig.1 shows the possible pathways of AN biosynthesis. In the plant, AN is biosynthesized from farnesyl pyrophosphate (FPP) produced from isopentenyl pyrophosphate (IPP), which is derived from the mevalonate (MVA) and the 2C-methyl-D-erythritol4-phosphate (MEP) pathway [
16,
19,
20]. FPP is converted into amorpha-4,11-diene cyclized by amorpha-4,11-diene synthase (ADS) [
21]. Amorpha-4,11-diene is then converted into artemisinic acid (AA). In this step, cytochrome P450 monooxygenase (CYP71AV1) hydroxylate amorpha-4,11-dieneto produces artemisinic alcohol, which oxidizes to generate AA via artemisinic aldehyde [
22,
23]. The final step is the generation of dihydroartemisinic acid (DHAA) formed from artemisinic aldehyde through dihydroartemisinic aldehyde. Scholars reported that the final step involves two genes, namely, a double bond reductase (DBR2) [
24] and an aldehyde dehydrogenase (ALDH1) [
25]. Furthermore, DHAA scavenges reactive oxygen species (ROS) and is converted into AN through a non-enzymatic spontaneous reaction [
26,
27]. In addition, arteannuin B (AB), which is generated from artemisinic acid (AA), may be converted into AN [
14].
Our previous work showed that Cd enhanced the accumulation of artemisinin and its precursors in
A. annua; however, the underlying mechanism remains unclear [
10,
28]. Our team aims to understand the effects of Cd on the accumulation of secondary metabolites in
A. annua and provide a basis for further studies on other pharmacologically active medicinal plants. In the present study, the effects and the mechanisms of Cd on the accumulation of other secondary metabolites such as AA, AB, and DHAA, were investigated using different Cd concentrations.
A. annua was grown in a controlled hydroponic environment and treated with different concentrations of Cd. Leaves were collected at different time periods after the treatment. Four secondary metabolites and eight key enzyme genes involved in AN synthesis were analyzed.
Materials and methods
Plant material and experimental design
A. annua seeds (guihao third) were obtained from the National Engineering Laboratory of South-west Endangered Medicinal Resources Development. The seeds were sown with 2/3 sterilized vermiculite and watered with distilled water and Hoagland solution every day. Thereafter, the plants were watered every 3 days. When seedlings reached a height of 5 cm, they were transferred to hydroponic pots with 1 L of Hoagland solution (pH 5.8). Three seedlings of the same height were transplanted into each pot. The Hoagland solution was changed every 2 days.
Plants were grown in a greenhouse with a 12 h/12 h (light/dark) photo-period with sodium lighting at 25 °C. The seedlings were grown for 40 days. Four different Cd (Cd(NO3)2·4H2O) concentrations was added to the Hoagland solution: 0 (control group), 20, 60, and 120 μmol/L. Each Cd level was replicated nine times, with each group consisting of three pots containing nine seedlings.
Solution preparation
Briefly, 20.0 mg of the AN standard was weighed and dissolved in 10 ml of methyl alcohol to obtain a standard stock solution (2000 μg/ml). The stock solution was diluted to 500.0, 200.0, 100.0, 50.0, 20.0, 10.0, and 5.0 μg/ml by adding methyl alcohol and stored at 4 °C. Moreover, 2 mg of AA, AB, and DHAA standards were weighed and dissolved in 10 ml of methyl alcohol to obtain stock solutions of 200 μg/ml. The solutions were diluted to 50.0, 20.0, 10.0, 5.0, 2.0, 2.0, and 0.5 μg/ml with methyl alcohol and stored at 4 °C.
One plant from each hydroponic pot was harvested at 0, 1, and 3 days (d) after Cd treatment and dried to a constant weight in a ventilated oven at 45 °C.
A. annua was extracted according to previously reported method with slight modifications [
29]. Dry crushed plant material (0.20 g) was weighed in an Eppendorf tube. The samples were sonicated in 2.5 ml of methanol for 20 min at 30 °C and centrifuged at 12 000 r/min for 5 min. The supernatant was transferred to a 10-ml volumetric flask. The procedure was repeated twice, and respective supernatants were combined. The final volume was adjusted to volume with methanol and mixed thoroughly. After filtration with 0.22 μm nylon membrane filter (Agela Technologies, Tianjin, China), the extract was used to quantify AA, AB, DHAA, and AN in
A. annua.
HPLC-ESI-MS/MS analysis
Chromatographic analyses were performed on an Agilent 1220 Infinity System (Agilent Corporation, CA, USA). HPLC separation was carried out on ZORBAX RRHD Eclipse Plus C18 column (2.1 mm × 50 mm, 1.8 μm, Agilent Technologies, MD, USA) with the thermostat set at 30 °C. The injection volume for all samples was 5 μl, and the flow rate was 0.200 ml/min. The binary gradient consisted of solvent system A (formic acid/acetonitrile 0.1:99.9 v/v) and solvent system B (formic acid/water 0.1:99.9 v/v). The conditions were as follows: 0 min, 25% A; 7 min, 98% A; 9.5 min, 98% A; 10 min, 25% A; 15 min, 25% A; and isocratic at 25% A from 15 to 20 min to create equilibrium in the column for the next injection. For ESI-MS/MS analysis, an Agilent QQQ 6420 Mass Spectrometer (Santa Clara, CA, USA) was connected on the same Agilent 1220 HPLC instrument via an electrospray ionization (ESI) interface. MS analysis was carried out in positive ionization mode (ESI+) by monitoring the protonated molecular ions under the following operating conditions: nebulizer gas pressure of 30.00 psi; dry gas flow rate of 10.00 L/min; ion source electrospray voltage of 4000 V; and capillary temperature of 350 °C. The MS parameters were manually optimized (Table 1). Quantification was performed using multiple reaction monitoring (MRM) modes for the above transitions. Data were acquired with Agilent MassHunter Workstation Acquisition software and analyzed using Agilent MassHunter Workstation Qual software.
Extraction of total RNA and real-time PCR analysis
Leaves were collected from A. annua at 0, 4, 8, 12, 24, 48, and 72 h after the treatment, immediately frozen in liquid nitrogen, and stored at -80 °C until further use. Total RNA was extracted from the samples using RNeasy Plant Mini Kit. DNase treatment was performed on an RNeasy Mini spin column with RNase-free DNase set (QIAGEN, Germany) according to the manufacturer’s protocol. Approximately 0.5 μg of total RNA was reverse transcribed using Primer ScriptTMII 1st Strand cDNA Synthesis Kit (Takara, Japan) according to the manufacturer’s instructions.
Real-time PCR (qPCR) was performed using specific primers designed by the online PrimerQuest Tool (Table 2) on 7500 qPCR equipment (Applied Biosystems, USA). First-strand cDNA was used as template in 10-μl reactions including 5 μl of Power SYBR Green PCR Master Mix (Applied Biosystems, USA) and 0.025 μg of each primer. 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), and 95 °C (15 s). Dissociation was performed to detect possible primer dimers. Each sample was performed in triplicate, and a negative control was included in all qPCR runs. Relative gene expression was calculated based on 2
-ΔΔCT comparative method [
30].
Statistical analysis
All experiments were conducted at least three times, and all the values are shown as mean±SE. Significant differences between the control and treatment groups were analyzed using one-way ANOVA procedure with SPSS 20.0 software.
Results
Effects of Cd on biosynthesis of AA, AB, DHAA, and AN in A. annua
Fig.2 shows changes in AA, AB, DHAA, and AN concentrations after treatment with Cd at different concentrations. 1 d after treatment, the concentrations of the four metabolites showed no significant difference in all the treated groups compared with the control group. However, after 3 d of treatment, the 60 and 120 μmol/L concentration groups showed a significant increase in AA, AB, and AN concentrations compared with the control. After 3 d of treatment with 120 μmol/L of Cd, AA, AB, and AN concentrations were 2.26 (P<0.05), 102.08 (P<0.01) and 33.63 (P<0.01) times that of the control. After 3 d of treatment with 60 μmol/L Cd, AB and AN concentrations were 61.10 (P<0.01) and 26.40 (P<0.01) times that of the control; the concentration of AA was slightly higher than that of the control (P>0.05). After 3 d of treatment with 20 μmol/L Cd, AA, AB, and AN concentrations were slightly higher than that of the control (P>0.05). After 3 d, the concentrations of DHAA in all treatment groups were slightly lower than that of the control (P>0.05).
Effects of Cd on AN biosynthesis genes
Studies have shown that the biosynthesis and accumulation of secondary metabolites in plants are closely associated with the expression of key enzyme genes in their biosynthesis pathways [
31–
34]. The results above show that Cd can significantly increase the concentrations of AA, AB, and AN in
A. annua. qPCR was used to further investigate the effects of different concentrations of Cd on the expression of key enzyme genes in the AN biosynthesis pathways in
A. annua, and the relationship between the key enzyme gene expression levels and the concentrations of the four metabolites. This method was used to determine the relative expression levels of eight genes (
HMGR,
FPS,
ADS,
CYP71AV1,
DBR2,
ALDH1,
DXS, and
DXR).
Expression trends of HMGR and FPS, key enzyme genes in the upstream of AN biosynthesis pathway, were similar in the control and in all the treatment groups. Relative expression levels increased after 4 h of treatment, reaching its peak at 8 h before decreasing, and remaining consistent at 24 h, 48 h, and 72 h (Fig. 3A and 3C). Comparing the control with the treatment groups at any time point, the relative expression levels of HMGR and FPS were highest in the 120 μmol/L group; the relative expression level of HMGR was 4.33 times that of the control at 8 h, and that of FPS was 6.95 times higher at 48 h (Fig. 3B and 3D).
ADS is widely acknowledged as the first key enzyme gene in the AN biosynthesis pathway. Relative expression level of ADS was highest at 4 h in both the control and treatment groups, then declined until reaching another peak at 48 h (Fig. 4A). The relative expression level of ADS was highest in the 120 μmol/L group compared with the control and the other treatment groups at any time point, with a relative expression level of ADS 112.68 times that of the control at 8 h. However, the relative expression of ADS was not detected at 24 h (Fig. 4B).
The expression trends of CYP71AV1, a multiple functional enzyme which catalyzes amorpha-4,11-diene successively into AA, were similar to that of the ADS in the control and in all treatment groups, with the highest level detected at 4 h, but then decreased before reaching another peak at 48 h (Fig. 4C). The relative expression of CYP71AV1 was highest in the 120 μmol/L group compared with that of the control and the other treatment groups at any time point, with its relative expression level 997.98 times that of the control at 48 h (Fig. 4D).
The expression trends of DBR2, a key enzyme specifically catalyzing artemisinic aldehyde into dihydroartemisinic aldehyde, were similar to that of HMGR and FPS, which increased at 4 h, reaching its peak at 8 h before declining, and remained constant at 24 h, 48 h, and 72 h (Fig. 5A). The relative expression level of DBR2 was highest in the 120 μmol/L group compared with the control and other treatment groups at any time point, with its relative expression level 62.73 times that of the control at 72 h (Fig. 5B).
ALDH1 simultaneously catalyzes artemisinic aldehyde into AA and dihydroartemisinic aldehyde into DHAA. The relative expression level of ALDH1 was highest in the 20 μmol/L group compared with the control and the other treatment groups at 4 h, 8 h, and 12 h, with its relative expression level 533.76 times that of the control at 8 h. The relative expression level of ALDH1 was highest in the 120 μmol/L group compared with the control and other treatment groups at 24 h, 48 h, and 72 h, with its relative expression level 188.70 times that of control at 72 h (Fig. 5D).
DXS and DXR, genes of key enzymes in the upstream MEP pathway, are involved in AN biosynthesis by catalyzing IPP, a precursor of AN. The expression trends of DXS and DXR were similar in the control and treatment groups, with relative expression level increasing at 4 h, reached a peak at 8 h, then decreased; remained stable at 24, 48, and 72 h (Fig. 6).
Discussion
Changes in the concentrations of AA, AB, and AN in
A. annua between the control and Cd-treated groups were almost the same at 1 d and 3 d after treatment. At 1 d, concentrations of all three secondary metabolites decreased, increased, and finally decreased with increasing Cd concentration. However, at day 3, the concentrations of these metabolites increased with increasing Cd concentration. The results indicated that AA and AB showed the same trends with AN in
A. annua after Cd treatment, but the accumulation trend of DHAA slightly differed from that of AN. Similar results were reported by Wang
et al. [
35]. To ensure that Cd influences the transcription of genes involved in artemisinin biosynthesis, we used real-time PCR for analyzing the relative expression levels of the following genes:
HMGR,
FPS,
ADS,
CYP71AV1,
DBR2,
ALDH1,
DXR, and DXS. The relative expression levels of these eight genes reached their peak at 4 h and 8 h and then declined. In addition, the abundance of these transcripts in response to different concentrations of Cd are different during this time period. Cd might be the specific signal that causes changes in the transcript levels of
HMGR,
FPS,
ADS,
CYP71AV1,
DBR2,
ALDH1, and
DXR in 120 μmol/L treatment group. However, in the 20 μmol/L and 60 μmol/L treatment groups,
ADS,
CYP71AV1,
ALDH1, and
DXS are the genes that showed changes in transcript levels.
Cd can increase the accumulation of secondary metabolites such as AN by up-regulating key enzyme genes in the AN biosynthesis pathways. The relative expression levels of HMGR, FPS, ADS, CYP71AV1, DBR2, ALDH1, and DXR were highest after treatment with 120 μmol/L of Cd. These observations indicate that the concentrations of AN, AA, and AB were the highest in the group treated with 120 μmol/L Cd for 3 d; whereas the relative expression levels of ADS, CYP71AV1, ALDH1, and DXS were up-regulated in the 20 and 60 μmol/L -treated groups. These findings suggest that different concentrations of Cd showed varied effects on the expression levels of key enzyme genes in AN biosynthesis pathways, which may be the reason why concentrations of metabolites like AN differed in A. annua after treatment with different concentrations of Cd.
Many strains of
A. annua, such as SP18 [
14], 001 strain [
16], and YU strain [
15], which produce a high yield of artemisinin, and can be divided into two chemotypes, AA/AB and AN/DHAA, based on the concentrations of AA, AB, DHAA, and AN. YU and 001 strains belong to AA/AB chemotype, whereas SP18 belongs to AN/DHAA [
36]. In our study, concentrations of AA, AB, DHAA, and AN in
A. annua (guihao third) were determined by HPLC-ESI-MS/MS, and results showed that the concentrations of AN and DHAA were higher than that of AB and AA. This observation suggests that, like SP18,
A. annua (guihao third) also belongs to the AN/DHAA chemotype. Research has showed that AN biosynthesis pathways differed between these two chemotypes [
36].
Interestingly, concentrations of DHAA in the control and treatment groups showed no significant difference. However, some studies have reported [
14,
31] that concentrations of DHAA significantly increased in
A. annua (SP18) after chitosan and MeJA treatment.
A. annua (SP18) and
A. annua (guihao third) are both the AN/DHAA chemotype, and this difference in DHAA concentrations between these two cultivars may lie in different elicitors. DHAA, the accepted direct precursor of AN, could turn into AN and eliminate ROS simultaneously by quenching singlet oxygen in a non-enzymatic fashion [
27]. This finding suggested that the presence of ROS is conducive to the transformation of DHAA to AN. Previous studies in our laboratory have shown that Cd could induce the generation of ROS in
A. annua [
10]. This result suggested that DHAA could turn into AN in
A. annua under Cd treatment. DHAA concentration did not significantly differ between the control and treatment groups. The concentrations of AN in groups treated with 60 and 120 μmol/L Cd were significantly higher than those in the control, suggesting that DHAA was rapidly and efficiently converted into the product AN by ROS produced after Cd treatment.
Higher Education Press and Springer-Verlag Berlin Heidelberg