Next Article in Journal
Detection and Identification of Estrogen Based on Surface-Enhanced Resonance Raman Scattering (SERRS)
Previous Article in Journal
Replacing Standard Reporters from Molecular Cloning Plasmids with Chromoproteins for Positive Clone Selection
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 23
Issue 6
10.3390/molecules23061329
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Red and Blue Light Promote the Accumulation of Artemisinin in Artemisia annua L.

by
Dong Zhang
1,2,
Wei Sun
2,
Yuhua Shi
1,2,
Lan Wu
1,2,
Tianyuan Zhang
2 and
Li Xiang
1,2,*
1
Artemisinin Reserch Center, China Academy of Chinese Medical Sciences, Beijing 100700, China
2
Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, Beijing 100700, China
*
Author to whom correspondence should be addressed.
Molecules 2018, 23(6), 1329; https://doi.org/10.3390/molecules23061329
Submission received: 27 April 2018 / Revised: 30 May 2018 / Accepted: 30 May 2018 / Published: 31 May 2018
(This article belongs to the Section Metabolites)
Browse Figures
Versions Notes

Abstract

:
Artemisinin, which has been isolated from Artemisia annua L., is the most effective antimalarial drug and has saved millions of lives. In addition, artemisinin and its derivatives have anti-tumor, anti-parasitic, anti-fibrosis, and anti-arrhythmic properties, which enhances the demand for these compounds. Improving the content of artemisinin in A. annua is therefore becoming an increasing research interest, as the chemical synthesis of this metabolite is not viable. Ultraviolet B and C irradiation have been reported to improve the artemisinin content in A. annua, but they are harmful to plant growth and development. Therefore, we screened other light sources to examine if they could promote artemisinin content without affecting plant growth and development. We found that red and blue light could enhance artemisinin accumulation by promoting the expression of the genes that were involved in artemisinin biosynthesis, such as amorpha-4,11-diene synthase (ADS) and cytochrome P450 monooxygenase (CYP71AV1) genes. Thus, in addition to being the main light sources for photosynthesis, red and blue light play a key role in plant secondary metabolism, and optimizing the combination of these light might allow for the productionof artemisinin-rich A. annua.

1. Introduction

Malaria, a serious disease that is caused by Plasmodium sp. infection, threatens human health and safety in 108 countries and regions. About 3300 million people are exposed to malaria, more than 100 million are infected, and nearly 800,000 die every year. The World Health Organization ranks it alongside acquired immune deficiency syndrome (AIDS) and tuberculosis as one of the world's top three public health problems [1]. Artemisinin-based therapies are the main and best treatments for malaria [2], and the discovery of artemisinin as well as its contribution to malaria treatment granted Professor Youyou Tu the 2015 Nobel Prize in Physiology or Medicine. Artemisinin and its derivatives have also manifested anti-tumor, anti-parasitic, anti-fibrosis, and anti-arrhythmic properties [3]. Therefore, an increasing demand for artemisinin and its derivatives is expected [4]. Presently, the main source of artemisinin is through direct extraction from Artemisia annua L. Although artemisinin has been chemically synthesized since 1983, this method has complex steps and thus is far from industrialization, as is the synthesis of artemisinin by cell culture. In addition, artemisinin can be obtained by heterologous synthesis in Escherichia coli or yeast [5]. Because the artemisinin content is very low in most varieties of A. annua, improving its content in this plant is an urgent requirement.
Plant secondary metabolites are greatly affected by various environmental factors, such as light, temperature, humidity, and pathogenic bacteria [6,7,8]. Generally, plants under biotic or abiotic stress produce a large number of secondary metabolites as a defense against such stress. Light is an important environmental factor that regulates a plant’s life cycle, including plant growth and development [9,10,11]. Several studies have reported that ultraviolet (UV)-B irradiation can promote the accumulation of flavonoids in buckwheat [12,13], by enhancing the expression of key metabolic genes, such as phe ammonia lyase (PAL), cinnamate-4-hydroxylase (C4H), chalcone synthase (CHS), chalcone isomerase (CHI), and flavonol synthase (FLS). In addition, resveratrol accumulation in mature grape berries is induced by UV-C irradiation [14,15,16,17]. It was also found that UV-B and -C irradiation could improve the content of artemisinin in A. annua. In the plants that were treated with UV-B and -C irradiation for a short time (14 days) before flowering, the artemisinin content increased by 1.63 and 1.75 times, respectively, compared with the control group [18].
Artemisinin is a sesquiterpene lactone with a peroxy bridge structure (i.e., 1,2,4-three dioxane ring) and its biosynthetic pathway has been clearly identified (Figure 1). In plants, this five carbon-skeleton compound has two independent biosynthetic pathways, which are the mevalonate (MVA) pathway and are located in the cytoplasm, and the methylerythritol phosphate (MEP) pathway, which occur in plastids [19]. The metabolic steps before the farnesyl diphosphate (FPP) synthesis are considered as an upstream of artemisinin biosynthesis, and FPP is formed by isopentenyl diphosphate (IPP), which is provided by the MVA and MEP pathways [19]. Under the catalysis of amorpha-4,11-diene synthase (ADS), FPP closed-loop produces amorpha-4,11-diene [20], which is oxidized to sequentially generate artenimol, artemisinic aldehyde, and artemisinic acid, which is catalyzed by the cytochrome P450 monooxygenase CYP71AV1 [21,22,23]. Under the catalysis of artemisinic aldehyde delta11(13)-reductase (DBR2) and aldehyde dehydrogenase 1 (ALDH1), artemisinic aldehyde is converted to dihydroartemisinic acid [24], which is converted to artemisinin.
Over the past decades, the biosynthetic pathway of artemisinin has been clearly understood. However, the transcriptional regulation of artemisinin is relatively slow and the synthesis and accumulation of secondary metabolites have specific temporal and spatial patterns [25]. Transcription factors (TFs), which are the proteins that bind to specific sequences of downstream target genes, and that activate or inhibit gene expression, play an important role in this type of regulation [26] and, in recent years, the discovery of some TFs has promoted the study of artemisinin transcriptional regulation [27,28,29,30,31]. Therefore, studying their function is crucial for understanding artemisinin biosynthesis and its regulation. It has been found that WRKY (named after its WRKY domain, which binds to the DNA sequence), APETALA2/ethylene-responsive (AP2/ERF), basic helix–loop–helix (bHLH), and basic leucine zipper (bZIP) are involved in the regulation of artemisinin biosynthesis. Ma et al. [27] cloned AaWRKY1, which can promote the expression of ADS in A. annua. Lu et al. [28] cloned the AP2/ERF TF AaORA, which was specifically expressed in A. annua trichomes, and by overexpressing the AaORA gene, this also promoted the overexpression of the artemisinin biosynthetic pathway genes ADS, CYP71AV1, and DBR2, resulting in an increase of artemisinin contents by 40–53% in overexpressing transgenic plants. Ji et al. [29] cloned AabHLH1, which also can promote the expression of artemisinin biosynthesis genes ADS and CYP71AV1. Recently, the bHLH TF AaMYC2, which was cloned by Tang et al. [30], was found to promote the expression of CYP71AV1 and DBR2 genes. Zhang et al. [31] found that the TF AabZIP1, which mediates abscisic acid (ABA) signaling, was able to promote the synthesis of artemisinin by promoting the expression of ADS and CYP71AV1 genes. Compared with the wild-type, the artemisinin content was increased from 0.7- to 1.5-fold in AabZIP1-OX lines.
The present study comprehensively evaluated the impact of different light waves (far-red, white, blue, and red) on the synthesis of artemisinin. Using transcriptome analysis, several candidate genes that transcriptionally regulate the synthesis of artemisinin were found, providing important information for further research regarding the regulation of artemisinin synthesis. The use of other light sources to promote the content of artemisinin not affecting the growth and development of A. annua have been discussed regarding the establishment of A. annua plantations for artemisinin intensive production.

2. Results and Discussion

2.1. Analysis of Artemisinin Content under Various Treatments

The internal transcribed spacer 2 (ITS2) ribosomal DNA sequence was characterized by high nucleotide variation and thus was widely used in phylogenetic and taxonomic studies. The ITS2 sequences of the samples that were collected in Hainan, China, for the purposes of the present study, completely matched the ITS2 sequence of A. annua (data not shown), which confirmed the identity of the collected plants. After sowing the A. annua seeds into the soil for four weeks under white light, the seedlings were transferred to darkness, white light, far-red light, red light, and blue light conditions for two days. After the light treatments, the chlorophyll content was detected, and the seedlings under red and blue light accumulated increased levels of chlorophyll compared with the white light treatment (Figure S1). In addition, the artemisinin and artemisinic acid content was analyzed by liquid chromatography–mass spectrometry (LC–MS). The results showed that light had an important effect on the content of artemisinin and artemisinic acid, as blue and red light notably increased the content of these compounds, compared with white light. Under far-red light and darkness conditions, the content of artemisinin was equal to and obviously lower than under the white light, respectively (Figure 2). Overall, these results indicated that light, especially blue and red light, could enhance the content of artemisinin and artemisinic acid.
Light not only provided energy for plant photosynthesis, but it also acted as an important environmental signal so as to regulate several types of plant growth and development processes, such as the transition from heterotrophic to autotrophic growth, which involves chlorophyll synthesis, chloroplast development, and photomorphogenesis [9,10,11]. However, light could be a stress factor regulating plant secondary metabolism. Although UV-B and -C light could improve the content of artemisinin in A. annua [18], UV light was harmful to organisms, causing DNA mutations. Thus, examining the effects of other light on the synthesis of artemisinin was crucial for cultivating A. annua plants that were able to produce high amounts of this metabolic compound. In plants, there were several important light receptors, including phytochromes, which perceived far-red and red light signals [32]; cryptochromes, which perceived blue light signals [33,34,35]; and UVR8, which perceived UV-B signals [36]. Far-red, red, and blue light were used to treat the A. annua seedlings. Our results showed that blue and red light notably increased the content of artemisinin compared with the white light (Figure 1). This was of great significance, as red and blue light did not harm plant growth [37].

2.2. Transcriptome Sequencing and Assembly

The samples were subjected to Illumina Xten paired-end sequencing. After removing the adaptor sequences and filtering the low quality and ambiguous reads, 713 million clean reads that contained 107 Gb of valid data, were acquired. A summary of the sequencing statistics is shown in Table 1. These reads were assembled into 651,000 transcripts and 271,932 unigenes using Trinity. The average length of the transcripts was 922 bp and the shortest sequence length, at 50% of the genome (N50), was 1501 bp. The average length of the unigenes was 711 bp and that of N50 was 1176 bp. The maximum length of both transcripts and unigenes was 18,284 bp, and their minimum length was 201 bp. However, the transcripts had 600,292,436 nucleotides in total, whereas the unigenes had only 193,418,027 nucleotides. The sequencing data have been deposited in the National Center for Biotechnology Information (NCBI) database (Accession number: SRP133983)

2.3. Functional Annotation and Classification

For functional annotation, all of the assembled unigenes were searched against the NCBI non-redundant protein database (NR), NCBI nucleotide sequences (NT), Swiss-Prot, Protein family (PFAM), euKaryotic Ortholog Groups (KOG), gene ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases (Table 2).
For the global functional analysis, 87,236 unigenes were annotated in the GO database and were classified into 56 functional groups, including 25 groups in the biological process, 21 in the cellular component, and 10 in the molecular function (Figure 3). Within the biological process, ‘cellular process’ (GO: 0009987; 47,389 unigenes) and ‘metabolic process’ (GO: 0008152; 43,972 unigenes) were predominant. In the cellular component category, the two main groups were the ‘cell’ (GO: 0005623; 26,159 unigenes) and ‘cell part’ (GO: 0044464; 26,129 unigenes). The terms ‘binding’ (GO: 0005488; 48,030 unigenes) and ‘catalytic activity’ (GO: 0003824; 36,550 unigenes) were the most common in the molecular function category. Comparing GO annotations of up-regulated genes under red and blue light treatments with those that were obtained under darkness, revealed that these genes were closely related to photosynthesis, which was consistent with previous reports stating that red and blue light were the main light sources for photosynthesis [37].
The KEGG pathway analysis, which was performed to identify metabolic pathways (Figure S2), revealed that 37,459 unigenes were grouped into 19 KEGG pathways, mainly including translation (5046), metabolism (4818), carbohydrate metabolism (3414), as well as folding, sorting, and degradation (3176).
The KOG pathway analysis (Figure S3) revealed that 36,681 unigenes were grouped into 26 pathways, mainly including posttranslational modification, protein turnover, chaperones (5141), translation, ribosomal structure and biogenesis (4750), general function prediction only (4728), and signal transduction mechanisms (3611).

2.5. Comparative Analysis of Transcriptional Profiles

In order to identify differentially expressed genes (DEGs) in different treatment groups, the fragments per kilobase of exon per million mapped fragments (FPKM) values of the assembled unigenes were calculated. Under red light and blue light treatments, 9855 and 7295 unigenes, accounting for 3.6% and 2.7% of the total unigenes (271,932), respectively, were significantly overexpressed by at least 2-fold (Figure 4A). In addition, under red light and blue light, 22,236 and 7791 unigenes, accounting for 8.2% and 2.9% of the total unigenes, respectively, were significantly underexpressed by at least 2-fold (Figure 4B). The Venn diagram that was constructed for the number of up- and down-regulated unigenes among the red and blue light treatments displays overlapping and unique unigenes to each condition.

2.6. Expression Analysis of the Genes Involved in Artemisinin Biosynthesis

The transcriptome profiles revealed that many of the differentially expressed unigenes were annotated as artemisinin biosynthesis genes (Figure 5). In the mevalonate (MVA) pathway, two steps that were catalyzed by 3-hydroxy-3-methyl-glutaryl CoA synthase (HMGS) and 3-hydroxy-3-methyl-glutaryl CoA reductase (HMGR) converted acetoacetyl-CoA to mevalonic acid. This was reported when the overexpressed HMGR in A. annua, artemisinin content in the transgenic lines was 38.9% higher than that in non-transgenic plants [38]. In our results, the expression of HMGR was almost identical under dark or white light (control) conditions, which suggested that light could not promote the HMGR expression. However, the expression of HMGS was 3.1-fold higher under white light than in darkness, and red and blue light had almost the same effect as white light. 1-deoxy-d-xylulose-5-phosphate synthase (DXS) and 1-deoxy-d-xylulose-5-phosphate reductoisomerase (DXR), which were the first and second enzymes of the methylerythritol 4-phosphate (MEP) pathway and played a major role in the overall regulation of this route [39], catalyzed the conversion of pyruvate and d-glyceraldehyde-3-phosphate to 2-C-methyl-d-erythritol-4-phosphate, and hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase (HDR) catalyzed this compound to isopentenyl-diphosphate and dimethylallyl-diphosphate in darkness. Overexpression of HDR1 increased the contents of artemisinin, and the suppression of HDR1 led to opposite results [40]. The expression of DXS was slightly higher under white light than in darkness. However, the DXS and HDR expressions under white light were 3.3- and 2-fold higher than in darkness, respectively. Isopentenyl diphosphate and dimethylallyl diphosphate were utilized as substrates by FPS, to form FPP. Banyai et al. overexpressed FPS in A. annua and gained a 2.5-fold increase of artemisinin content at best [41]. Under white light, the expression of FPS was 2.5-fold higher than in darkness. Overall, the above results indicated that, by upstream FPP synthesis, light could have promoted the accumulation of artemisinin precursors and that white light was better than monochromatic light for enhancing such accumulation.
Downstream FPP synthesis, ADS, and CYP71AV1 were the most important enzymes. The content of artemisinin was increased by about 82% in the ADS overexpressing transgenic plant lines [27]. Jing et al. [42] overexpressed CYP71AV1 and CPR, and the artemisinin content was increased by as much as 2.4-fold compared with the control plant. Compared to darkness, the red light, white light, and blue light led to ADS overexpression by 2.5-, 5.9-, and 14.7-fold, respectively. Red light, white light, and blue light also led to the overexpression of CYP71AV1 by 24.3-, 16.7-, and 46.3-fold, respectively. Interestingly, the far-red light inhibited both of the enzymes, indicating that light promoted artemisinin synthesis in a complex process, as different light might have had different effects on the artemisinin biosynthesis gene expression. Blue light, in particular, promoted the expression of ALDH1 and DBR2. The contents of artemisinin were remarkably increased in the transgenic plants of A. annua with DBR2 overexpression [43]. Blue light inhibited the expression of ESC and SQS, which catalyzed FPP to 8-epi-cedrol and squalene, respectively. A study showed that suppressing the expression of SQS significantly increased the artemisinin content in transgenic A. annua [44]. Hence, light, especially blue light, not only promoted the expression of artemisinin synthetase, but it could have also inhibited the expression of synthetases in the branching pathway, which allowed more metabolites to flow to artemisinin.

2.7. Validation of Differential Expression via Quantitative Real-Time PCR (qRT-PCR)

The levels of DEGs were validated using qRT-PCR and eight selected genes. The qRT-PCR expression profile showed a positive correlation with the transcriptome data (Figure 6 and Table S1).

2.8. Identification of Co-Expressed TFs and Artemisinin Biosynthesis Genes

For more information about the regulation of the genes that were involved in the artemisinin biosynthesis, transcription factors were identified using the iTAK software. We identified a total of 5106 TFs. Then, the co-expression analysis was carried out to obtain the specific TFs that might have regulated the artemisinin biosynthesis gene expression (Table S2), which revealed that 365 TFs were co-expressed with TRINITY_DN75917_c2_g1 (ADS) and TRINITY_DN83022_c1_g1 (CYP71AV1) (|r| > 0.95). Within these 365 TFs, 29 were MYB TFs, 18 were bHLH TFs, 15 were bZIP TFs, and 12 were WRKY TFs (Figure 7). These TF families had been reported to play important roles in regulating artemisinin biosynthesis gene expression.
In order to validate the TF gene expression results, eight TFs that were co-expressed with ADS and CYP71AV1 were selected for qRT-PCR analysis (Figure 8), and the qRT-PCR expression profile showed a positive correlation with the transcriptome data. Three MYB TFs, three bZIP TFs, one IAA TF, and one ERF TF showed the highest expression under blue light, followed by red light, except MYB4, which showed the highest expression under red light.
In the plants, the synthesis and accumulation of secondary metabolites were under strict control from the TFs. Transcriptional regulation was a complex network, as a single TF might have regulated the expression of many genes that were involved in a specific metabolic pathway, and many TFs regulated the expression of a specific gene [25]. Several TF families, including WRKY, AP2/ERF, bHLH, and bZIP, were involved in the regulation of artemisinin biosynthesis. Most of these TFs were involved in hormone signals, but the TF that was regulated by light was still unknown. To explore which TFs might have been regulated by light, a co-expression analysis was carried out. The results revealed that MYB TFs were the most frequent TFs in the process, and that their expression was induced by light. In addition, 18 bHLH TFs and 15 bZIP TFs were found (Figure 7). In the light signal transduction pathway, the well-known bZIP TF HY5 was the central point of the pathway [45,46]. The phytochrome interacting factors (PIFs) belonging to the bHLH family were also very important light signal TFs, which regulated many of the development processes [47]. Twelve WRKY TFs were also found, and the WRKY TFs were usually related with light stress response and were induced by light when plants transitioned from heterotrophic to autotrophic growth [48]. Many TFs were found by co-expression analysis and were induced by light (Figure 8), and thus might have been involved in the light signal pathway, interacting with HY5 or PIFs. Overall, the results that were presented here provided candidate TFs for the transcriptional regulation of artemisinin by light, which was of great significance for further in-depth studies regarding artemisinin synthesis regulation.

3. Materials and Methods

3.1. Plant Material and Growth Conditions

The Artemisia annua that was used in this study was a wild type, which we named ‘Haiqing 1’, was obtained from Haikou, Hainan, China. The Artemisia annua seeds were sown to the soil for four weeks under white light at 25 °C. The photoperiod in the growth chamber was 16/8 h (light/dark) (16 h light and 8 h darkness per day). Then, the seedlings were transferred to four continuous light treatments for two days before harvesting. After different light treatments, we collected the aboveground parts of 10 seedlings randomly, as a sample approximately 0.8 g. Then, the sample was grinded into powder and put into liquid nitrogen. The content of artemisinin was detected using 0.1 g of the powder, and the other 0.1 g was used to extract RNA for sequencing. Light-emitting diode (LED) red light (670 nm), LED blue light (470 nm), and LED far-red light (735 nm). The light intensity was 50 ± 5 µmol/ m2 s in all of the treatments. Each treatment comprised three biological replicates (10 seedlings per repeat).

3.2. Chlorophyll Measurement

The chlorophyll measurement was conducted as previously described [49]. Chlorophyll content was expressed as 7.18 × OD663 + 17.32 × OD646 per 20 seedlings. All of the experiments were performed in triplicate.

3.3. Transcriptome Sequencing

The process was described previously [50], the plant total RNA was isolated using an RNA Extraction Kit (Tiangen, Beijing, China) according to the manufacturer’s instructions. To build the sequencing library, two micrograms of total RNA, with 28S/18S RNA ratio ≥1.8 were used by the NEBNext® Ultra RNA Library Prep Kit (New England Biolabs Inc., Ipswich, MA, USA). The sequencing library quality was determined with Agilent 2100 Bioanalyzer, which had a minimum RNA integrity number (RIN) value of 7 and 250–300 bp insertion elements. The RNA sequencing was performed using the Illumina Xten sequencing system (Illumina Inc., San Diego, CA, USA).

3.4. De Novo Assembly

In order to obtain high quality clean reads, the adapter sequences, low-quality sequences with N (undetermined bases) percentage over 10%, and the sequences containing more than 50% bases with q-value ≤5 and shorter than 35 bases, were removed. Clean reads were assembled into contigs de novo by the Trinity software with the min_kmer_cov set to 2 and all of the other parameters set to default values [51]. The program TransDecoder was used to predict the complete open reading frames (ORFs) of the unigenes (https://transdecoder.github.io/).

3.5. Functional Annotation

For the functional annotation, the resulting contigs were blasted against the National Center for Biotechnology Information (NCBI, http://www.ncbi.nlm.nih.gov) non-redundant protein database (NR) and nucleotide sequences (NT) (e-value = 1e−5), the Pfam database hmmscan (e-value = 0.01) (http://pfam.sanger.ac.uk/), the Swiss-Prot (e-value = 1e−5) (http://www.ebi.ac.uk/uniprot/), euKaryotic Ortholog Groups (KOG) (e-value = 1e−5) (http://www.ncbi.nlm.nih.gov/COG/), and the Kyoto Encyclopedia of Genes and Genomes (KEGG), using the KEGG Automatic Annotation Server (e-value = 1e−10) (http://www.genome.jp/kegg/). Subsequently, the gene ontology (GO) annotation was performed using the Blast2GO v2.5 (e-value =1e−6) (http://www.geneontology.org/) [52].

3.6. Transcription Factor Identification and Co-Expression Analysis

The transcription factor families were identified using the iTAKsoftware. In order to select the possible TFs regulating artemisinin biosynthesis, two highly expressed genes—TRINITY_DN75917_c2_g1 (ADS) and TRINITY_DN83022_c1_g1 (CYP71AV1)—were used in co-expression analysis, with a default value 0.6 [53]. The paired genes showing a Pearson correlation coefficient (r) greater than 0.95 were considered as significantly co-expressed and were selected to build a co-expression network, using the Perl script [50]. The data correlation and visualization were performed using Cytoscape v3.4.10 (National Institute of General Medical Sciences, Bethesda, MD, U.S.) [54].

3.7. Quantitative Real-Time PCR (qRT-PCR) Analysis of Gene Expression

The plant total RNA was isolated using an RNA Extraction Kit (Tiangen) and the first-strand cDNA was synthesized from 2 mg of RNA by reverse transcriptase (Invitrogen, Carlsbad, CA, USA), and then diluted (1:100) for use in qRT-PCR, with SYBR Premix ExTaq Mix (Takara, Dalian, Liaoning, China) in a total volume of 15 mL. The reactions were performed in a LightCycler 480 thermal cycler (Roche, Basel, Switzerland), following the manufacturer’s instructions. Three biological replicates were performed for each sample, and the expression level was normalized to that of actin, which was used as the control gene.

3.8. Liquid Chromatography–Mass Spectrometry (LC–MS) Analysis of Secondary Metabolites

Of the sample powder, 0.1 g was extracted overnight at 4 °C with 1.0 mL pure methanol (or 70% aqueous methanol) containing 0.1 mg·L−1 lidocaine. The lipid-solubility extracts were absorbed and 0.4 mL of each extract was mixed and filtrated (SCAA-104, 0.22 µm pore size) before LC–MS analysis. Artemisinin and artemisinic acid reference substances were accurately weighed (0.001 g) and used to obtain a stock solution (1 mg/mL), which was dissolved by adding the corresponding volume of methanol. The stock solutions of artemisinin and artemisinic acid standards were diluted to 100 µg/mL, 1600 ng/mL, 800 ng/mL, 400 ng/mL, 200 ng/mL, 100 ng/mL, 50 ng/mL, 25 ng/mL, and 10 ng/mL. The standard curves were determined as Y = 600.39X + 6012.6, and the linear range was 10–800 ng/mL.

3.9. High-Performance Liquid Chromatography (HPLC) Conditions

The analytical conditions were as follows: column was Agilent Elipse Plus C18 column (Santa Clara, CA, USA., pore size 1.8 µm, length 2.1 × 50 mm); solvent system was water (0.1% formic acid), acetonitrile (0.1% formic acid); gradient program was 100:0 v/v at 0 min, 5:95 v/v, at 20.0 min, 5:95 v/v at 22.0 min, 95:5 v/v at 22.1 min, and 95:5 v/v at 28 min; flow rate was 0.2 mL·min−1; temperature was 40 °C; injection volume was 1 µL. The effluent was alternatively connected to an electrospray ionization (ESI)-triple quadrupole-linear ion trap (Q-TRAP)-MS or ESI-QqTOF-MS.

4. Conclusions

Using LC–MS analysis, the present study revealed that red and blue light promoted artemisinin and artemisinic acid content in A. annua seedlings. Additionally, GO annotations revealed that the up-regulated genes under red and blue light treatments were closely related to photosynthesis, and the heatmap and qRT-PCR analyses revealed that red and blue light led to the up-regulation of many artemisinin biosynthesis genes and inhibited the expression of synthetases in the branching pathway. The co-expression was analysis carried out to reveal the TFs that might have regulated the artemisinin biosynthesis gene expression, and the qRT-PCR analysis confirmed that the expression of MYB39 and bZIP9, among others, were induced by red and blue light. Overall, the results that have been presented here revealed that red and blue light promoted the accumulation of artemisinin and artemisinic acid in A. annua.

Supplementary Materials

The chlorophyll content under different light treatments, Figure S2: Functional classification and pathway assignment of assembled unigenes by KEGG, Figure S3: The KOG annotation for all assembled unigenes. Table S1: List of primers used in qRT-PCR, Table S2: List of co-expressed transcript factors.

Author Contributions

D.Z. and T.Z. analyzed the data; D.Z. wrote the manuscript draft; W.S. and L.W. prepared the material for sequencing; Y.S. and W.S. revised the manuscript; and L.X. designed the experiments and revised the manuscript.

Funding

This work was funded by funding from the University of Macau (MYRG2016-00055-ICMS-QRCM), Hong Kong, Macao and the Taiwan Science and Technology Cooperation Project (2015DFM30030), the Fundamental Research Funds for the Central public welfare research institutes (ZZ10-024) and Major New Drug Creation ‘Study on the Innovative Drugs of Artemisinin and its Derivatives’ (2017ZX09101002-003-001).

Acknowledgments

The authors thank Qi Shen and Mingkun Huang, as well as the Qinggang Yin Institute of Chinese Materia Medica and the China Academy of Chinese Medical Sciences for the helpful discussion and language polishing.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. World Health Organization. World Malaria Report 2015. 2015. Available online: http://www.who.int/mediacentre/factsheets/fs094/en/ (accessed on 9 December 2015).
  2. Dondorp, A.M.; Nosten, F.; Yi, P.; Das, D.; Phyo, A.P.; Tarning, J.; Lwin, K.M.; Ariey, F.; Hanpithakpong, W.; Lee, S.J.; et al. Artemisinin Resistance in Plasmodium falciparum Malaria. N. Engl. J. Med. 2009, 361, 455–467. [Google Scholar] [CrossRef] [PubMed]
  3. Bin, L.; Hong, Z. Advances on the study of the pharmacological effects of artemisinin and its derivatives. Chin. J. Clin. Pharmacol. Ther. 2010, 15, 572–576. [Google Scholar]
  4. Peplow, M. Malaria drug made in yeast causes market ferment. Nature 2013, 494, 160–161. [Google Scholar] [CrossRef] [PubMed]
  5. Paddon, C.J.; Westfall, P.J.; Pitera, D.J.; Benjamin, K.; Fisher, K.; McPhee, D.; Leavell, M.D.; Tai, A.; Main, A.; Eng, D.; et al. High-level semi-synthetic production of the potent antimalarial artemisinin. Nature 2013, 496, 528–532. [Google Scholar] [CrossRef] [PubMed]
  6. Aziz, A.; Poinssot, B.; Daire, X.; Adrian, M.; Bézier, A.; Lambert, B.; Joubert, J.M.; Pugin, A. Laminarin elicits defense responses in grapevine and induces protection against Botrytis cinerea and Plasmopara viticola. Mol. Plant-Microbe. Interact. 2003, 16, 1118–1128. [Google Scholar] [CrossRef] [PubMed]
  7. Korn, M.; Peterek, S.; Mock, H.P.; Heyer, A.G.; Hincha, D.K. Heterosis in the freezing tolerance, and sugar and flavonoid contents of crosses between Arabidopsis thaliana accessions of widely varying freezing tolerance. Plant Cell Environ. 2008, 31, 813–827. [Google Scholar] [CrossRef] [PubMed]
  8. Tsurunaga, Y.; Takahashi, T.; Katsube, T.; Kudo, A.; Kuramitsu, O.; Ishiwata, M.; Matsumoto, S. Effects of UV-B irradiation on the levels of anthocyanin, rutin and radical scavenging activity of buckwheat sprouts. Food Chem. 2013, 141, 552–556. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, H.Y.; Deng, X.W. Dissecting the phytochrome A-dependent signaling network in higher plants. Trends Plant Sci. 2003, 8, 172–178. [Google Scholar] [CrossRef]
  10. Jiao, Y.; Lau, O.S.; Deng, X.W. Light-regulated transcriptional networks in higher plants. Genetics 2007, 8, 217–230. [Google Scholar] [CrossRef] [PubMed]
  11. Lau, O.S.; Deng, X.W. Plant hormone signaling lightens up: Integrators of light and hormones. Curr. Opin. Plant Biol. 2010, 13, 571–577. [Google Scholar] [CrossRef] [PubMed]
  12. Wilson, M.I.; Greenberg, B.M. Protection of the D1 photosystem II reaction center protein from degradation in ultraviolet radiation following adaptation of Brassica napus L. to growth in ultraviolet radiation. Photochem. Photobiol. 1993, 57, 556–563. [Google Scholar] [CrossRef]
  13. Landrey, L.G.; Chapple, C.C.S.; Last, R.L. Arabidopsis mutants lacking phenolic sunscreens exhibit enhanced ultraviolet-B injury and oxidative damage. Plant Physiol. 1995, 109, 1159–1166. [Google Scholar] [CrossRef]
  14. Adrian, M.; Jeandet, P.; Douillet-Breuil, A.C.; Tesson, L.; Bessis, R. Stilbene content of 3 mature Vitis vinifera berries in response to UV-C elicitation. J. Agric. Food Chem. 2000, 48, 6103–6105. [Google Scholar] [CrossRef] [PubMed]
  15. Versari, A.; Parpinello, G.P.; Tornielli, G.B.; Ferrarini, R.; Giulivo, C. Stilbene compounds and stilbene synthase expression during ripening, wilting, and UV treatment in grape cv. Corvina. J. Agric. Food Chem. 2001, 49, 5531–5536. [Google Scholar] [CrossRef] [PubMed]
  16. Takayanagi, T.; Okuda, T.; Mine, Y.; Yokotsuka, K. Induction of resveratrol biosynthesis in skins of three grape cultivars by ultraviolet irradiation. J. Jpn. Soc. Hortic. Sci. 2004, 73, 193–199. [Google Scholar] [CrossRef]
  17. Nishikawa, Y.; Tominori, S.; Wada, K.; Kondo, H. Effect of cultivation practices on resveratrol content in grape berry skins. Hortic. Res. 2011, 10, 249–253. [Google Scholar] [CrossRef]
  18. Rai, R.; Meena, R.P.; Smita, S.S.; Shukla, A.; Rai, K.; Pandey-Rai, S. UV-B and UV-C pre-treatments induce physiological changes and artemisinin biosynthesis in Artemisia annua L.—An antimalarial plant. J. Photochem. Photobiol. B Biol. 2011, 105, 216–225. [Google Scholar] [CrossRef] [PubMed]
  19. Schramek, N.; Wang, H.; Römisch-Margl, W.; Keil, B.; Radykewicz, T.; Winzenhörlein, B.; Beerhues, L.; Bacher, A.; Rohdich, F.; Gershenzon, J.; et al. Artemisinin biosynthesis in growing plants of Artemisia annua. A 13CO2 study. Phytochemistry 2010, 71, 179–187. [Google Scholar] [CrossRef] [PubMed]
  20. Nguyen, K.T.; Arsenault, P.R.; Weathers, P.J. Trichomes + roots + ROS = artemisinin: Regulating artemisinin biosynthesis in Artemisia annua L. In Vitro Cell. Dev. Biol.-Plant 2011, 47, 329–338. [Google Scholar] [CrossRef] [PubMed]
  21. Wang, H.; Han, J.; Kanagarajan, S.; Lundgren, A.; Brodelius, P.E. Trichome-specific expression of the amorpha-4,11-diene 12-hydroxylase (cyp71av1) gene, encoding a key enzyme of artemisinin biosynthesis in Artemisia annua, as reported by a promoter-GUS fusion. Plant Mol. Biol. 2013, 81, 119–138. [Google Scholar] [CrossRef] [PubMed]
  22. Wang, Y.; Yang, K.; Jing, F.; Li, M.; Deng, T.; Huang, R.; Wang, B.; Wang, G.; Sun, X.; Tang, K.X. Cloning and Characterization of Trichome Specific Promoter of cpr71av1 Gene Involved in Artemisinin Biosynthesis in Artemisia annua L. Mol. Biol. 2011, 45, 751–758. [Google Scholar] [CrossRef]
  23. Teoh, K.H.; Polichuk, D.R.; Reed, D.W.; Nowak, G.; Covello, P.S. 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. [Google Scholar] [CrossRef] [PubMed]
  24. Fang, X.; Lu, S.; Yu, Z.X.; Chen, X.Y. Research in biosynthesis of artemisinin (in Chinese). Sci. Technol. Rev. 2015, 33, 31–35. [Google Scholar]
  25. Yang, C.Q.; Fang, X.; Wu, X.M.; Mao, Y.; Wang, L.; Chen, X. Transcriptional regulation of plant secondary metabolism. J. Integr. Plant Biol. 2012, 54, 703–712. [Google Scholar] [CrossRef] [PubMed]
  26. Patra, B.; Schluttenhofer, C.; Wu, Y.; Pattanaik, S.; Yuan, L. Transcriptional regulation of secondary metabolite biosynthesis in plants. BBA-Gene Regul. Mech. 2013, 1829, 1236–1247. [Google Scholar] [CrossRef] [PubMed]
  27. Ma, D.; Pu, G.; Lei, C.; Ma, L.; Wang, H.; Guo, Y.; Chen, J.; Du, Z.; Wang, H.; Li, G.; et al. Isolation and characterization of AaWRKY1, an Artemisia annua transcription factor that regulates the amorpha-4,11-diene synthase gene, a key gene of artemisinin biosynthesis. Plant Cell Physiol. 2009, 50, 2146–2161. [Google Scholar] [CrossRef] [PubMed]
  28. Lu, X.; Zhang, L.; Zhang, F.; Jiang, W.; Shen, Q.; Zhang, L.; Lv, Z.; Wang, G.; Tang, K. AaORA, a trichome-specific AP2/ERF transcription factor of Artemisia annua, is a positive regulator in the artemisinin biosynthetic pathway and in disease resistance to Botrytis cinerea. New Phytol. 2013, 198, 1191–1202. [Google Scholar] [CrossRef] [PubMed]
  29. Ji, Y.; Xiao, J.; Shen, Y.; Ma, D.; Li, Z.; Pu, G.; Li, X.; Huang, L.; Liu, B.; Ye, H.; et al. Cloning and Characterization of AabHLH1, a bHLH Transcription Factor that Positively Regulates Artemisinin Biosynthesis in Artemisia annua. Plant Cell Physiol. 2014, 55, 1–13. [Google Scholar] [CrossRef] [PubMed]
  30. Tang, K.X.; Shen, Q.; Lu, X.; Fu, X.Q.; Yan, T.X.; Sun, X.F.; Wang, G.F. Sweet Wormwood MYC2 Transcription Factor Protein Coding Sequence and Applications Thereof. China Patent ZL201310413155.3, 11 September 2013. [Google Scholar]
  31. Zhang, F.; Fu, X.; Lv, Z.; Lu, X.; Shen, Q.; Zhang, L.; Zhu, M.; Wang, G.; Sun, X.; Liao, Z.; et al. A Basic Leucine Zipper Transcription Factor, AabZIP1, Connects Abscisic Acid Signaling with Artemisinin Biosynthesis in Artemisia annua L. Mol. Plant 2015, 8, 163–175. [Google Scholar] [CrossRef] [PubMed]
  32. Rockwell, N.C.; Su, Y.S.; Lagarias, J.C. Phytochrome structure and signaling mechanisms. Annu. Rev. Plant Biol. 2006, 57, 837–858. [Google Scholar] [CrossRef] [PubMed]
  33. Lin, C.; Robertson, D.E.; Ahmad, M.; Raibekas, A.A.; Jorns, M.S.; Dutton, P.L.; Cashmore, A.R. Association of flavin adenine dinucleotide with the Arabidopsis blue light receptor CRY1. Science 1995, 269, 968–970. [Google Scholar] [CrossRef] [PubMed]
  34. Malhotra, K.; Kim, S.T.; Batschauer, A.; Dawut, L.; Sancar, A. Putative blue-light photoreceptors from Arabidopsis thaliana and Sinapis alba with a high degree of sequence homology to DNA photolyase contain the two photolyase cofactors but lack DNA repair activity. Biochemistry 1995, 34, 6892–6899. [Google Scholar] [CrossRef] [PubMed]
  35. Hsu, D.S.; Zhao, X.; Zhao, S.; Kazantsev, A.; Wang, R.P.; Todo, T.; Wei, Y.F.; Sancar, A. Putative human blue-light photoreceptors hCRY1 and hCRY2 are flavoproteins. Biochemistry 1996, 35, 13871–13877. [Google Scholar] [CrossRef] [PubMed]
  36. Brown, B.A.; Cloix, C.; Jiang, G.H.; Kaiserli, E.; Herzyk, P.; Kliebenstein, D.J.; Jenkins, G.I. A UV-B-specific signaling component orchestrates plant UV protection. Proc. Natl. Acad. Sci. USA 2005, 102, 18225–18230. [Google Scholar] [CrossRef] [PubMed]
  37. Moore, R.; Clark, W.D.; Kingsley, R.S.; Vodopich, D. Plant physiology and growth. In Botany; Wm. C. Brown Pub: Dubuque, IA, USA, 1995; Chapter 7; ISBN 9780697037756. [Google Scholar]
  38. Nafis, T.; Akmal, M.; Ram, M.; Alam, P.; Ahlawat, S.; Mohd, A.; Abdin, M.Z. Enhancement of artemisinin content by constitutive expression of the HMG-CoA reductase gene in highyielding strain of Artemisia annua L. Plant Biotechnol. Rep. 2011, 5, 53–60. [Google Scholar] [CrossRef]
  39. 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. [Google Scholar] [CrossRef] [PubMed]
  40. Ma, D.; Li, G.; Zhu, Y.; Xie, D.Y. Overexpression and Suppression of Artemisia annua 4-Hydroxy- 3-Methylbut-2-enyl Diphosphate Reductase 1 Gene (AaHDR1) Differentially Regulate Artemisinin and Terpenoid Biosynthesis. Front. Plant Sci. 2017, 8, 77. [Google Scholar] [CrossRef] [PubMed]
  41. Banyai, W.; Kirdmanee, C.; Mii, M.; Supaibulwatana, K. Overexpression of farnesyl pyrophosphate synthase (FPS) gene affected artemisinin content and growth of Artemisia annua L. Plant Cell Tissue Organ. Cult. 2010, 103, 255–265. [Google Scholar] [CrossRef]
  42. Jing, F.; Zhang, L.; Li, M.; Tang, K. Over-expressing CYP71AV1 and CPR genes enhances artemisinin content in Artemisia annua L. J. Agric. Sci. Technol. 2008, 10, 64–70. [Google Scholar]
  43. Yuan, Y.; Liu, W.; Zhang, Q.; Xiang, L.; Liu, X.; Chen, M.; Lin, Z.; Wang, Q.; Liao, Z. Overexpression of artemisinic aldehyde Δ11(13) reductase gene-enhanced artemisinin and its relative metabolite biosynthesis in transgenic Artemisia annua L. Biotechnol. Appl. Biochem. 2015, 62, 17–23. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, L.; Jing, F.Y.; Li, F.P.; Li, M.Y.; Wang, Y.L.; Wang, G.F.; Sun, X.F.; Tang, K.X. 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. [Google Scholar] [CrossRef] [PubMed]
  45. Osterlund, M.T.; Hardtke, C.S.; Wei, N.; Deng, X.W. Targeted destabilization of HY5 during light-regulated development of Arabidopsis. Nature 2000, 405, 462–466. [Google Scholar] [CrossRef] [PubMed]
  46. Saijo, Y.; Sullivan, J.A.; Wang, H.; Yang, J.; Shen, Y.; Rubio, V.; Ma, L.; Hoecker, U.; Deng, X.W. The COP1-SPA1 interaction defines a critical step in phytochrome A-mediated regulation of HY5 activity. Genes Dev. 2003, 17, 2642–2647. [Google Scholar] [CrossRef] [PubMed]
  47. Heim, M.A.; Jakoby, M.; Werber, M.; Martin, C.; Weisshaar, B.; Bailey, P.C. The basic helix-loop-helix transcription factor family in plants: A genome-wide study of protein structure and functional diversity. Mol. Biol. Evol. 2003, 20, 735–747. [Google Scholar] [CrossRef] [PubMed]
  48. Rushton, D.L.; Tripathi, P.; Rabara, R.; Liu, J.; Ringler, P.; Boken, A.K.; Langum, T.J.; Smidt, L.; Boomsma, D.D.; Emme, N.J.; et al. WRKY transcription factors: Key components in abscisic acid signalling. Plant Biotechnol. J. 2012, 10, 2–11. [Google Scholar] [CrossRef] [PubMed]
  49. Lin, R.; Ding, L.; Casola, C.; Ripoll, D.R.; Feschotte, C.; Wang, H. Transposase-derived transcription factors regulate light signaling in Arabidopsis. Science 2007, 318, 1302–1305. [Google Scholar] [CrossRef] [PubMed]
  50. Zhang, T.; Song, C.; Song, L.; Shang, Z.; Yang, S.; Zhang, D.; Sun, W.; Shen, Q.; Zhao, D. RNA Sequencing and Coexpression Analysis Reveal Key Genes Involved in α-Linolenic Acid Biosynthesis in Perilla frutescens Seed. Int. J. Mol. Sci. 2017, 18, 2433. [Google Scholar] [CrossRef] [PubMed]
  51. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q. Trinity: Reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat. Biotechnol. 2013, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  52. Götz, S.; García-Gómez, J.M.; Terol, J.; Williams, T.D.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talón, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef] [PubMed]
  53. Ariani, A.; Gepts, P. Genome-wide identification and characterization of aquaporin gene family in common bean (Phaseolus vulgaris L.). Mol. Genet. Genom. 2015, 290, 1771. [Google Scholar] [CrossRef] [PubMed]
  54. Smoot, M.; Ono, K.; Ideker, T.; Maere, S. PiNGO: A Cytoscape plugin to find candidate genes in biological networks. Bioinformatics 2011, 27, 1030–1031. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Samples of the compounds are available from the authors.
Molecules 23 01329 g001 550
Figure 1. Schematic representation of artemisinin biosynthesis in Artemisia annua L. HMGS—3-hydroxy-3-methyl-glutaryl coenzyme A synthase; HMGR—3-hydroxy-3-methyl-glutaryl coenzyme A reductase; DXS—1-deoxy-d-xylulose-5-phosphate synthase; DXR—1-deoxy-Dxylulose-5-phosphate reductoisomerase; HDR—hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase; FDS—farnesyl diphosphate synthase; SQS—squalene synthase; ECS—epi-cedrol synthase; ADS—amorpha-4,11-diene synthase; CYP71AV1—amorphadiene-12-hydroxylase; DBR2—artemisinic aldehyde Δ11(13) reductase; ALDH1—aldehyde dehydrogenase 1; MEP—methylerythritol phosphate.
Figure 1. Schematic representation of artemisinin biosynthesis in Artemisia annua L. HMGS—3-hydroxy-3-methyl-glutaryl coenzyme A synthase; HMGR—3-hydroxy-3-methyl-glutaryl coenzyme A reductase; DXS—1-deoxy-d-xylulose-5-phosphate synthase; DXR—1-deoxy-Dxylulose-5-phosphate reductoisomerase; HDR—hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase; FDS—farnesyl diphosphate synthase; SQS—squalene synthase; ECS—epi-cedrol synthase; ADS—amorpha-4,11-diene synthase; CYP71AV1—amorphadiene-12-hydroxylase; DBR2—artemisinic aldehyde Δ11(13) reductase; ALDH1—aldehyde dehydrogenase 1; MEP—methylerythritol phosphate.
Molecules 23 01329 g001
Molecules 23 01329 g002 550
Figure 2. Artemisinin (A) and artemisinic acid (B) content under different light treatments. In (A) and (B), data represent means ± standard deviation of biological triplicates (10 seedlings per repeat). D—darkness; FR—far red light; WL—white light; R—red light; B—blue light; FW—fresh weight. Double asterisks indicate significant differences at p < 0.01 using Student’s t-test.
Figure 2. Artemisinin (A) and artemisinic acid (B) content under different light treatments. In (A) and (B), data represent means ± standard deviation of biological triplicates (10 seedlings per repeat). D—darkness; FR—far red light; WL—white light; R—red light; B—blue light; FW—fresh weight. Double asterisks indicate significant differences at p < 0.01 using Student’s t-test.
Molecules 23 01329 g002
Molecules 23 01329 g003 550
Figure 3. Functional annotation of all of the assembled unigenes into biological process, cellular component, and molecular function categories within the gene ontology (GO) database.
Figure 3. Functional annotation of all of the assembled unigenes into biological process, cellular component, and molecular function categories within the gene ontology (GO) database.
Molecules 23 01329 g003
Molecules 23 01329 g004 550
Figure 4. The Venn diagram shows the number of up-regulated (A) and down-regulated (B) unigenes in the two treatment groups. The numbers in the circles indicate the number of unigenes regulated by one or two treatment.
Figure 4. The Venn diagram shows the number of up-regulated (A) and down-regulated (B) unigenes in the two treatment groups. The numbers in the circles indicate the number of unigenes regulated by one or two treatment.
Molecules 23 01329 g004
Molecules 23 01329 g005 550
Figure 5. Heatmap showing variations in the expression of the genes involved in artemisinin biosynthesis. D—darkness; FR—far red light; WL—white light; R—red light; B—blue light. HMGS—3-hydroxy-3-methyl-glutaryl coenzyme A synthase; HMGR—3-hydroxy-3-methyl-glutaryl coenzyme A reductase; DXS—1-deoxy-d-xylulose-5-phosphate synthase; DXR—1-deoxy-Dxylulose-5-phosphate reductoisomerase; HDR—hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase; FPS—farnesyl diphosphate synthase; SQS—squalene synthase; ECS—epi-cedrol synthase; ADS—amorpha-4,11-diene synthase; CYP71AV1—amorphadiene-12-hydroxylase; DBR2—artemisinic aldehyde Δ11(13) reductase; ALDH1—aldehyde dehydrogenase 1. Color indicates the strength of the expression signal, with gene expression becoming stronger from green to red.
Figure 5. Heatmap showing variations in the expression of the genes involved in artemisinin biosynthesis. D—darkness; FR—far red light; WL—white light; R—red light; B—blue light. HMGS—3-hydroxy-3-methyl-glutaryl coenzyme A synthase; HMGR—3-hydroxy-3-methyl-glutaryl coenzyme A reductase; DXS—1-deoxy-d-xylulose-5-phosphate synthase; DXR—1-deoxy-Dxylulose-5-phosphate reductoisomerase; HDR—hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase; FPS—farnesyl diphosphate synthase; SQS—squalene synthase; ECS—epi-cedrol synthase; ADS—amorpha-4,11-diene synthase; CYP71AV1—amorphadiene-12-hydroxylase; DBR2—artemisinic aldehyde Δ11(13) reductase; ALDH1—aldehyde dehydrogenase 1. Color indicates the strength of the expression signal, with gene expression becoming stronger from green to red.
Molecules 23 01329 g005
Molecules 23 01329 g006 550
Figure 6. Validation of RNA sequencing results by real-time quantitative PCR. D—darkness; FR—far red light; WL—white light; R—red light; B—blue light. ADS—amorpha-4,11-diene synthase; CYP71AV1—amorphadiene-12-hydroxylase; DBR2—artemisinic aldehyde Δ11(13) reductase; ALDH1—aldehyde dehydrogenase 1; HDR—hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase; HMGS—3-hydroxy-3-methyl-glutaryl coenzyme A synthase; SQS—squalene synthase; ECS—epi-cedrol synthase. The relative expression levels were normalized to that of the actin control. Data represent means ± standard deviation of biological triplicates.
Figure 6. Validation of RNA sequencing results by real-time quantitative PCR. D—darkness; FR—far red light; WL—white light; R—red light; B—blue light. ADS—amorpha-4,11-diene synthase; CYP71AV1—amorphadiene-12-hydroxylase; DBR2—artemisinic aldehyde Δ11(13) reductase; ALDH1—aldehyde dehydrogenase 1; HDR—hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase; HMGS—3-hydroxy-3-methyl-glutaryl coenzyme A synthase; SQS—squalene synthase; ECS—epi-cedrol synthase. The relative expression levels were normalized to that of the actin control. Data represent means ± standard deviation of biological triplicates.
Molecules 23 01329 g006
Molecules 23 01329 g007 550
Figure 7. Distribution of transcription factor (TF) families in the co-expression results. There were 365 TFs that were co-expressed with TRINITY_DN75917_c2_g1 (ADS) and TRINITY_DN83022_c1_g1 (CYP71AV1) (r > 0.95).
Figure 7. Distribution of transcription factor (TF) families in the co-expression results. There were 365 TFs that were co-expressed with TRINITY_DN75917_c2_g1 (ADS) and TRINITY_DN83022_c1_g1 (CYP71AV1) (r > 0.95).
Molecules 23 01329 g007
Molecules 23 01329 g008 550
Figure 8. Quantitative real time PCR results obtained for selected transcription factors. D—darkness; FR—far red light; WL—white light; R—red light; B—blue light. The relative expression levels were normalized to that of the actin control. Data represent means ± standard deviation of biological triplicates.
Figure 8. Quantitative real time PCR results obtained for selected transcription factors. D—darkness; FR—far red light; WL—white light; R—red light; B—blue light. The relative expression levels were normalized to that of the actin control. Data represent means ± standard deviation of biological triplicates.
Molecules 23 01329 g008
Table
Table 1. Summary of the Artemisia annua transcriptome data obtained by Illumina sequencing.
Table 1. Summary of the Artemisia annua transcriptome data obtained by Illumina sequencing.
SampleClean Reads (n)Clean Bases (n)Q20 (%)Q30 (%)GC Content (%)
D141,315,65439,620,94497.6393.6642.76
D244,4941,9842,460,92897.3292.9642.92
D345,997,16244,196,85697.6293.6342.85
FR144,583,75842,993,87097.6093.5842.92
FR246,616,43644,924,62897.6593.8042.92
FR353,563,70851,433,96697.1692.7942.93
WL148,061,22246,239,32697.5693.6242.67
WL241,232,12239,771,77297.6393.7742.75
WL351,498,24249,412,74497.6893.8742.82
R146,278,94044,459,75896.9092.1242.62
R250,660,67848,157,44697.3293.0142.55
R350,543,01247,291,92097.2192.9042.61
B163,937,45861,502,15497.9594.6442.50
B260,062,10257,833,94097.9194.5442.52
B355,204,49852,916,96497.5993.6042.51
D—darkness; FR—far-red light; WL—white light; R—red light; B—blue light. Numbers correspond to three biological duplicates. Q20 and Q30: the base value of Phred is greater than 20 and 30, respectively, accounting for the total base percentage.
Table
Table 2. Unigenes annotated per database.
Table 2. Unigenes annotated per database.
DatabaseNumber of UnigenesPercentage (%)
NR90,04933.11
NT67,31124.75
KEGG37,45913.77
Swiss-Prot83,50630.70
PFAM85,70231.51
GO87,23632.08
KOG36,68113.48
Annotated in all of the databases13,3264.90
Annotated in at least one database136,36550.14
Total Unigenes271,932100.00
NCBI—National Center for Biotechnology Information; NCBI NR—NCBI non-redundant database; NCBI NT—NCBI nucleotide sequences; KEGG—Kyoto Encyclopedia of Genes and Genomes; SwissProt—a manually annotated and reviewed protein sequence database; PFAM—Protein family; GO—Gene Ontology; KOG—euKaryotic Ortholog Groups.

Share and Cite

MDPI and ACS Style

Zhang, D.; Sun, W.; Shi, Y.; Wu, L.; Zhang, T.; Xiang, L. Red and Blue Light Promote the Accumulation of Artemisinin in Artemisia annua L. Molecules 2018, 23, 1329. https://doi.org/10.3390/molecules23061329

AMA Style

Zhang D, Sun W, Shi Y, Wu L, Zhang T, Xiang L. Red and Blue Light Promote the Accumulation of Artemisinin in Artemisia annua L. Molecules. 2018; 23(6):1329. https://doi.org/10.3390/molecules23061329

Chicago/Turabian Style

Zhang, Dong, Wei Sun, Yuhua Shi, Lan Wu, Tianyuan Zhang, and Li Xiang. 2018. "Red and Blue Light Promote the Accumulation of Artemisinin in Artemisia annua L." Molecules 23, no. 6: 1329. https://doi.org/10.3390/molecules23061329

APA Style

Zhang, D., Sun, W., Shi, Y., Wu, L., Zhang, T., & Xiang, L. (2018). Red and Blue Light Promote the Accumulation of Artemisinin in Artemisia annua L. Molecules, 23(6), 1329. https://doi.org/10.3390/molecules23061329

Article Metrics

Back to TopTop