Overexpression of SQUAMOSA promoter binding protein-like 4a (NtSPL4a) alleviates Cd toxicity in Nicotiana tabacum.

Squamosa Promoter Binding Protein-Like (SPL) plays a crucial role in regulating plant development and combating stress, yet its mechanism in regulating resistance to Cd toxicity remains unclear. In this study, we cloned a nuclear-localized transcription factor, NtSPL4a, from the tobacco cultivar TN90. Transient co-expression results showed that miR156 significantly reduced the expression of NtSPL4a by binding to the 3'-UTR of its transcript. We obtained transgenic tobacco overexpressing NtSPL4a (including the 3'-UTR) and NtSPL4aΔ (lacking the 3'-UTR) through Agrobacterium-mediated genetic transformation. Compared to the wild type (WT), overexpression of NtSPL4a/NtSPL4aΔ shortened the flowering time and exhibited a more developed root system. The transgenic tobacco showed significantly reduced Cd content, being 85.1% (OE-NtSPL4a) and 46.7% (OE-NtSPL4aΔ) of WT, respectively. Moreover, the upregulation of NtSPL4a affected the mineral nutrient homeostasis in transgenic tobacco. Additionally, overexpression of NtSPL4a/NtSPL4aΔ effectively alleviated leaf chlorosis and oxidative stress induced by Cd toxicity. One possible reason is that the overexpression of NtSPL4a/NtSPL4aΔ can effectively promote the accumulation of non-enzymatic antioxidants. A comparative transcriptomic analysis was performed between transgenic tobacco and WT to further unravel the global impacts brought by NtSPL4a. The tobacco overexpressing NtSPL4a had 183 differentially expressed genes (77 upregulated, 106 downregulated), while the tobacco overexpressing NtSPL4aΔ had 594 differentially expressed genes (244 upregulated, 350 downregulated) compared to WT. These differentially expressed genes mainly included transcription factors, metal transport proteins, flavonoid biosynthesis pathway genes, and plant stress-related genes. Our study provides new insights into the role of the transcript factor SPL in regulating Cd tolerance.

Biosynthesis of Scopoletin in Sweet Potato Confers Resistance against Fusarium oxysporum.

Fusarium wilt is a severe fungal disease caused by Fusarium oxysporum in sweet potato. We conducted transcriptome analysis to explore the resistance mechanism of sweet potato against F. oxysporum. Our findings highlighted the role of scopoletin, a hydroxycoumarin, in enhancing resistance. In vitro experiments confirmed that scopoletin and umbelliferone had inhibitory effects on the F. oxysporum growth. We identified hydroxycoumarin synthase genes IbF6'H2 and IbCOSY that are responsible for scopoletin production in sweet potatoes. The co-overexpression of IbF6'H2 and IbCOSY in tobacco plants produced the highest scopoletin levels and disease resistance. This study provides insights into the molecular basis of sweet potato defense against Fusarium wilt and identifies valuable genes for breeding wilt-resistant cultivars.

Overexpression of NtGPX8a Improved Cadmium Accumulation and Tolerance in Tobacco (Nicotiana tabacum L.).

Cadmium (Cd)-induced oxidative stress detrimentally affects hyperaccumulator growth, thereby diminishing the efficacy of phytoremediation technology aimed at Cd pollution abatement. In the domain of plant antioxidant mechanisms, the role of glutathione peroxidase (GPX) in conferring Cd tolerance to tobacco (Nicotiana tabacum) remained unclear. Our investigation employed genome-wide analysis to identify 14 NtGPX genes in tobacco, revealing their organization into seven subgroups characterized by analogous conserved domain patterns. Notably, qPCR analysis highlighted NtGPX8a as markedly responsive to Cd2+ stress. Subsequent exploration through yeast two-hybridization unveiled NtGPX8a's utilization of thioredoxins AtTrxZ and AtTrxm2 as electron donors, and without interaction with AtTrx5. Introduction of NtGPX8a into Escherichia coli significantly ameliorated Cd-induced adverse effects on bacterial growth. Transgenic tobacco overexpressing NtGPX8a demonstrated significantly augmented activities of GPX, SOD, POD, and CAT under Cd2+ stress compared to the wild type (WT). Conversely, these transgenic plants exhibited markedly reduced levels of MDA, H2O2, and proline. Intriguingly, the expression of NtGPX8a in both E. coli and transgenic tobacco led to increased Cd accumulation, confirming its dual role in enhancing Cd tolerance and accumulation. Consequently, NtGPX8a emerges as a promising candidate gene for engineering transgenic hyperaccumulators endowed with robust tolerance for Cd-contaminated phytoremediation.

Heterologous expression of the tobacco metallothionein gene NtMT2F confers enhanced tolerance to Cd stress in Escherichia coli and Arabidopsis thaliana.

Heavy metal pollution in the soil is a serious threat to crop growth and human health. Metallothionein (MT) is a low molecular weight protein that is rich in cysteine, which can effectively alleviate the toxicity of heavy metals in plants. In this study, a novel metallothionein encoding gene, NtMT2F, was cloned from the Cd-hyperaccumulator tobacco and heterologously expressed in E. coli and A. thaliana to verify its biological function. Recombinant E. coli incubated with NtMT2F effectively resisted heavy metal stress, particularly Cd. The recombinant strain grew significantly faster and had a higher content of Cd than the control. Mutations in the C-terminal Cys residues of NtMT2F significantly reduced its ability to chelate heavy metals. The overexpression of NtMT2F significantly enhanced resistance to Cd toxicity in transgenic A. thaliana. The germination rate, root length, and fresh weight of transgenic plants under Cd stress were higher than those of the wild type (WT). The contents of hydrogen peroxide (H2O2) and malondialdehyde (MDA) were lower than those of the WT. In addition, the activities of anti-peroxidase enzymes including glutathione reductase (GR), catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD), were significantly increased in the transgenic plants. The results of this study indicate that NtMT2F significantly improved the tolerance of microorganisms and plants to Cd and could be an important candidate protein for phytoremediation.

The NtSPL Gene Family in Nicotiana tabacum: Genome-Wide Investigation and Expression Analysis in Response to Cadmium Stress.

The SQUAMOSA promoter binding protein-like (SPL)SPL family genes play an important role in regulating plant growth and development, synthesis of secondary metabolites, and resistance to stress. Understanding of the role of the SPL family in tobacco is still limited. In this study, 42 NtSPL genes were identified from the genome of the tobacco variety TN90. According to the results of the conserved motif and phylogenetic tree, the NtSPL genes were divided into eight subgroups, and the genes in the same subgroup showed similar gene structures and conserved domains. The cis-acting element analysis of the NtSPL promoters showed that the NtSPL genes were regulated by plant hormones and stresses. Twenty-eight of the 42 NtSPL genes can be targeted by miR156. Transcriptome data and qPCR results indicated that the expression pattern of miR156-targeted NtSPL genes was usually tissue specific. The expression level of miR156 in tobacco was induced by Cd stress, and the expression pattern of NtSPL4a showed a significant negative correlation with that of miR156. These results suggest that miR156-NtSPL4a may mediate the tobacco response to Cd stress. This study lays a foundation for further research on the function of the NtSPL gene and provides new insights into the involvement of NtSPL genes in the plant response to heavy metal stress.

AabHLH113 integrates jasmonic acid and abscisic acid signaling to positively regulate artemisinin biosynthesis in Artemisia annua.

Artemisinin, a sesquiterpene lactone isolated from Artemisia annua, is in huge market demand due to its efficient antimalarial action, especially after the COVID-19 pandemic. Many researchers have elucidated that phytohormones jasmonic acid (JA) and abscisic acid (ABA) positively regulate artemisinin biosynthesis via types of transcription factors (TFs). However, the crosstalk between JA and ABA in regulating artemisinin biosynthesis remains unclear. Here, we identified a novel ABA- and JA-induced bHLH TF, AabHLH113, which positively regulated artemisinin biosynthesis by directly binding to the promoters of artemisinin biosynthetic genes, DBR2 and ALDH1. The contents of artemisinin and dihydroartemisinic acid increased by 1.71- to 2.06-fold and 1.47- to 2.23-fold, respectively, in AabHLH1113 overexpressed A. annua, whereas they decreased by 14-36% and 26-53%, respectively, in RNAi-AabHLH113 plants. Furthermore, we demonstrated that AabZIP1 and AabHLH112, which, respectively, participate in ABA and JA signaling pathway to regulate artemisinin biosynthesis, directly bind to and activate the promoter of AabHLH113. Collectively, we revealed a complex network in which AabHLH113 plays a key interrelational role to integrate ABA- and JA-mediated regulation of artemisinin biosynthesis.

AabHLH112, a bHLH transcription factor, positively regulates sesquiterpenes biosynthesis in Artemisia annua.

The bHLH transcription factors play important roles in the regulation of plant growth, development, and secondary metabolism. ß-Caryophyllene, epi-cedrol, and ß-farnesene, three kinds of sesquiterpenes mainly found in plants, are widely used as spice in the food industry and biological pesticides in agricultural production. Furthermore, they also have a significant value in the pharmaceutical industry. However, there is currently a lack of knowledge on the function of bHLH family TFs in ß-caryophyllene, epi-cedrol, and ß-farnesene biosynthesis. Here, we found that AabHLH112 transcription factor had a novel function to positively regulate ß-carophyllene, epi-cedrol, and ß-farnesene biosynthesis in Artemisia annua. Exogenous MeJA enhanced the expression of AabHLH112 and genes of ß-caryophyllene synthase (CPS), epi-cedrol synthase (ECS), and ß-farnesene synthase (BFS), as well as sesquiterpenes content. Dual-LUC assay showed the activation of AaCPS, AaECS, and AaBFS promoters were enhanced by AabHLH112. Yeast one-hybrid assay showed AabHLH112 could bind to the G-box (CANNTG) cis-element in promoters of both AaCPS and AaECS. In addition, overexpression of AabHLH112 in A. annua significantly elevated the expression levels of AaCPS, AaECS, and AaBFS as well as the contents of ß-caryophyllene, epi-cedrol, and ß-farnesene, while suppressing AabHLH112 expression by RNAi reduced the expression of the three genes and the contents of the three sesquiterpenes. These results suggested that AabHLH112 is a positive regulator of ß-caryophyllene, epi-cedrol, and ß-farnesene biosynthesis in A. annua.

Molecular insights into AabZIP1-mediated regulation on artemisinin biosynthesis and drought tolerance in Artemisia annua.

Artemisia annua is the main natural source of artemisinin production. In A. annua, extended drought stress severely reduces its biomass and artemisinin production while short-term water-withholding or abscisic acid (ABA) treatment can increase artemisinin biosynthesis. ABA-responsive transcription factor AabZIP1 and JA signaling AaMYC2 have been shown in separate studies to promote artemisinin production by targeting several artemisinin biosynthesis genes. Here, we found AabZIP1 promote the expression of multiple artemisinin biosynthesis genes including AaDBR2 and AaALDH1, which AabZIP1 does not directly activate. Subsequently, it was found that AabZIP1 up-regulates AaMYC2 expression through direct binding to its promoter, and that AaMYC2 binds to the promoter of AaALDH1 to activate its transcription. In addition, AabZIP1 directly transactivates wax biosynthesis genes AaCER1 and AaCYP86A1. The biosynthesis of artemisinin and cuticular wax and the tolerance of drought stress were significantly increased by AabZIP1 overexpression, whereas they were significantly decreased in RNAi-AabZIP1 plants. Collectively, we have uncovered the AabZIP1-AaMYC2 transcriptional module as a point of cross-talk between ABA and JA signaling in artemisinin biosynthesis, which may have general implications. We have also identified AabZIP1 as a promising candidate gene for the development of A. annua plants with high artemisinin content and drought tolerance in metabolic engineering breeding.

The cold-induced transcription factor bHLH112 promotes artemisinin biosynthesis indirectly via ERF1 in Artemisia annua.

Basic helix-loop-helix (bHLH) proteins are the second largest family of transcription factors (TFs) involved in developmental and physiological processes in plants. In this study, 205 putative bHLH TF genes were identified in the genome of Artemisia annua and expression of 122 of these was determined from transcriptomes used to construct the genetic map of A. annua. Analysis of gene expression association allowed division of the 122 bHLH TFs into five groups. Group V, containing 15 members, was tightly associated with artemisinin biosynthesis genes. Phylogenetic analysis indicated that two bHLH TFs, AabHLH106 and AabHLH112, were clustered with Arabidopsis ICE proteins. AabHLH112 was induced by low temperature, while AabHLH106 was not. We therefore chose AabHLH112 for further examination. AabHLH112 was highly expressed in glandular secretory trichomes, flower buds, and leaves. Dual-luciferase assays demonstrated that AabHLH112 enhanced the promoter activity of artemisinin biosynthesis genes and AaERF1, an AP2/ERF TF that directly and positively regulates artemisinin biosynthesis genes. Yeast one-hybrid assays indicated that AabHLH112 could bind to the AaERF1 promoter, but not to the promoters of artemisinin biosynthesis genes. Overexpression of AabHLH112 significantly up-regulated the expression levels of AaERF1 and artemisinin biosynthesis genes and consequently promoted artemisinin production.

Molecular Characterization of the 1-Deoxy-D-Xylulose 5-Phosphate Synthase Gene Family in Artemisia annua.

Artemisia annua produces artemisinin, an effective antimalarial drug. In recent decades, the later steps of artemisinin biosynthesis have been thoroughly investigated; however, little is known about the early steps of artemisinin biosynthesis. Comparative transcriptomics of glandular and filamentous trichomes and 13CO2 radioisotope study have shown that the 2-C-methyl-D-erythritol-4-phosphate (MEP) pathway, rather than the mevalonate pathway, plays an important role in artemisinin biosynthesis. In this study, we have cloned three 1-deoxy-D-xylulose 5-phosphate synthase (DXS) genes from A. annua (AaDXS1, AaDXS2, and AaDXS3); the DXS enzyme catalyzes the first and rate-limiting enzyme of the MEP pathway. We analyzed the expression of these three genes in different tissues in response to multiple treatments. Phylogenetic analysis revealed that each of the three DXS genes belonged to a distinct clade. Subcellular localization analysis indicated that all three AaDXS proteins are targeted to chloroplasts, which is consistent with the presence of plastid transit peptides in their N-terminal regions. Expression analyses revealed that the expression pattern of AaDXS2 in specific tissues and in response to different treatments, including methyl jasmonate, light, and low temperature, was similar to that of artemisinin biosynthesis genes. To further investigate the tissue-specific expression pattern of AaDXS2, the promoter of AaDXS2 was cloned upstream of the ß-glucuronidase gene and was introduced in arabidopsis. Histochemical staining assays demonstrated that AaDXS2 was mainly expressed in the trichomes of Arabidopsis leaves. Together, these results suggest that AaDXS2 might be the only member of the DXS family in A. annua that is involved in artemisinin biosynthesis.

[Molecular cloning and characterization of CMK from Artemisia annua].

Artemisinin is a preferred medicine in the treatment of malaria. In this study, AaCMK, a key gene involved in the upstream pathway of artemisinin biosynthesis, was cloned and characterized from Artemisia annua for the first time. The full-length cDNA of AaCMK was 1 462 bp and contained an ORF of 1 197 bp that encoded a 399-anomo-acid polypeptide. Tissue expression pattern analysis showed that AaCMK was expressed in leaves, flowers, roots and stems, but with higher expression level in glandular secretory trichomes. In addition, the expression of AaCMK was markedly increased after MeJA treatment. Subcellular localization showed that the protein encoded by AaCMK was localized in chloroplast. Overexpression of AaCMK in Arabidopsis increased the contents of chlorophyll a, chlorophyll b and carotenoids. These results suggest that AaCMK plays an important role in the biosynthesis of terpenoids in A. annua and this research provids a candidate gene that could be used for engineering the artemisinin biosynthesis.

ARTEMISININ BIOSYNTHESIS PROMOTING KINASE 1 positively regulates artemisinin biosynthesis through phosphorylating AabZIP1.

The plant Artemisia annua produces the anti-malarial compound artemisinin. Although the transcriptional regulation of artemisinin biosynthesis has been extensively studied, its post-translational regulatory mechanisms, especially that of protein phosphorylation, remain unknown. Here, we report that an ABA-responsive kinase (AaAPK1), a member of the SnRK2 family, is involved in regulating artemisinin biosynthesis. The physical interaction of AaAPK1 with AabZIP1 was confirmed by multiple assays, including yeast two-hybrid, bimolecular fluorescence complementation, and pull-down. AaAPK1, mainly expressed in flower buds and leaves, could be induced by ABA, drought, and NaCl treatments. Phos-tag mobility shift assays indicated that AaAPK1 phosphorylated both itself and AabZIP1. As a result, the phosphorylated AaAPK1 significantly enhanced the transactivational activity of AabZIP1 on the artemisinin biosynthesis genes. Substituting the Ser37 with Ala37 of AabZIP1 significantly suppressed its phosphorylation, which inhibited the transactivational activity of AabZIP1. Consistent overexpression of AaAPK1 significantly increased the production of artemisinin, as well as the expression levels of the artemisinin biosynthesis genes. Our study opens a window into the regulatory network underlying artemisinin biosynthesis at the post-translational level. Importantly, and for the first time, we provide evidence for why the kinase gene AaAPK1 is a key candidate for the metabolic engineering of artemisinin biosynthesis.

[Molecular cloning and characterization of the 2-C-methyl-D-erythritol 4-phosphate cytidylyltransferase gene from Artemisia annua L.].

The plastidial methylerythritol phosphate(MEP) pathway provides 5-carbon precursors to the biosynthesis of isoprenoid (including artemisinin). 2-C-Methyl-D-erythritol-4-phosphate cytidylyltransferase (MCT) is the third enzyme of the MEP pathway, which catalyzes 2-C-methyl-D-erythritol-4-phosphate to form 4-(cytidine 5'-diphospho)-2-C-methyl-D-erythritol. The full-length MCT cDNA sequence (AaMCT) was cloned and characterized for the first time from Artemisia annua L. Analysis of tissue expression pattern revealed that AaMCT was highly expressed in glandular secretory trichome and poorly expressed in leaf, flower, root and stem. AaMCT was found to be a methyl jasmonate (Me JA)-induced genes, the expression of AaMCT was significantly increased after MeJA treatment. Subcellular localization indicated that the GFP protein fused with AaMCT was targeted specifically in chloroplasts. The transgenic plants of Arabidopsis thaliana with AaMCT overexpression exhibited a significantly increase in the content of chlorophyll a, chlorophyll b and carotenoids, demonstrating that AaMCT kinase plays an influential role in isoprenoid biosynthesis.

[Molecular cloning and functional characterization of the gene encoding hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase gene from Artemisia annua L.].

Artemisinin is the first choice for malaria treatment. The plastidial MEP pathway provides 5-carbon precursors (IPP and its isomer DMAPP) for the biosynthesis of isoprenoid (including artemisinin). Hydroxy-2-methyl-2-(E)-butenyl 4-diphosphate reductase (HDR) is the last enzyme involved in the MEP pathway, which catalyzes HMBPP to form IPP and DMAPP. In this study, we isolated the full-length cDNA of HDR from Artemisia annua L. (AaHDR2) and performed functional analysis. According to gene expression analysis of AaHDR2 (GenBank: KX058541) and AaHDR1 reported ever (GenBank: ADC84348.1) by qPCR, we found that AaHDR1 and AaHDR2 had much higher expression level in trichomes than that in roots, stems, leaves and flowers. AaHDR2 had much higher expression level in flowers than that in leaves. Further, the plant hormones such as Me JA and ABA respectively up-regulated the expression level of AaHDR1 and AaHDR2 significantly, but GA3 up-regulated the expression level of AaHDR2 only. The gene expression analysis of AaHDR1 and AaHDR2 showed that AaHDR2 had a greater contribution than AaHDR1 to isoprenoid biosynthesis(including artemisinin). We used AaHDR2 for the following experiments. Bioinformatic analysis indicated that AaHDR2 belonged to the HDR family and the functional complementation assay showed that AaHDR2 did have the enzymatic function of HDR, using E. coli mutant MG1655(ara)<>HDR as host cell. The subcellular localization assay showed that AaHDR2 fused with GFP at its N-terminal specifically targeted in chloroplasts. Finally, AaHDR2 was overexpressed in Arabidopsis thaliana. The AaHDR2-overexpressing plants produced the isoprenoids including chlorophyll a, chlorophyll b and carotenoids at significantly higher levels than the wild-type Arabidopsis plants. In summary, AaHDR2 might be a candidate gene for genetic improvement of the isoprenoid biosynthesis.

Overexpression of artemisinic aldehyde Δ11 (13) reductase gene-enhanced artemisinin and its relative metabolite biosynthesis in transgenic Artemisia annua L.

Artemisinic aldehyde Δ11 (13) reductase (DBR2) is the checkpoint enzyme catalyzing artemisinic aldehyde to form dihydroartemisinic aldehyde directly involved in artemisinin biosynthetic pathway. In the present study, DBR2 was employed to engineer the biosynthetic pathway of artemisinin in transgenic plants of Artemisia annua L. Seven independent transgenic plants of A. annua with DBR2 overexpression driven by the cauliflower mosaic virus 35S promoter were obtained by Agrobacterium-mediated genetic transformation and confirmed by genomic PCR. The results of real-time qPCR analysis showed that the expression levels of DBR2 gene in all the seven transgenic lines were significantly higher than in nontransgenic control. The high-performance liquid chromatography analysis of artemisinin and its relative metabolites demonstrated that the contents of artemisinin and its direct precursor dihydroartemisinic acid were remarkably increased in the transgenic plants of A. annua with DBR2 overexpression. Interestingly, it was also found that the contents of arteannuin B and its direct precursor artemisinic acid in the branch pathway competing against artemisinin biosynthesis were also improved in DBR2-overexpressed A. annua plants. The transgenic results in the present study indicated that DBR2 is a useful structural gene in engineering the artemisinin biosynthetic pathway to develop genetically modified A. annua with the higher yield of artemisinin.

[Enhancement of artemisinin biosynthesis in transgenic Artemisia annua L. by overexpressed HDR and ADS genes].

Artemisnin is a novel sesquiterpene lactone with an internal peroxide bridge structure, which is extracted from traditional Chinese herb Artemisia annua L. (Qinghao). Recommended by World Health Organization, artemisinin is the first-line drug in the treatment of encephalic and chloroquine-resistant malaria. In the present study, transgenic A. annua plants were developed by overexpressing the key enzymes involved in the biosynthetic pathway of artemisinin. Based on Agrobacterium-mediated transformation methods, transgenic plants of A. annua with overexpression of both HDR and ADS were obtained through hygromycin screening. The genomic PCR analysis confirmed six transgenic lines in which both HDR and ADS were integrated into genome. The gene expression analysis given by real-time quantitative PCR showed that all the transgenic lines had higher expression levels of HDR and ADS than the non-transgenic control (except ah3 in which the expression level of ADS showed no significant difference compared with control); and the HPLC analysis of artemisinin demonstrated that transgenic A. annua plants produced artemisinin at significantly higher level than non-transgenic plants. Especially, the highest content of artemisinin was found in transgenic line ah70, in which the artemisinin content was 3.48 times compared with that in non-transgenic lines. In summary, overexpression of HDR and ADS facilitated artemisinin biosynthesis and this method could be applied to develop transgenic plants of A. annua with higher yield of artemisinin.

[Relative expression of genes involved in artemisinin biosynthesis and artemisinin accumulation in different tissues of Artemisia annua].

OBJECTIVE: To study the relative expression of the genes involved in artemisinin biosynthesis in different tissues including roots, stems, leaves and flowers of Artemisia annua, and establish the relationship between gene expression and artemisinin accumulation, eventually leading to discover the mainly effective genes involved in artemisinin biosynthesis. METHOD: The 7 functional genes involved in artemisinin biosynthesis were detected at the level of expression by using qRT-PCR, and simultaneously the content of artemisinin in the 4 investigated tissues was detected in parallel. RESULT: The 3 genes including HMGR, DXR and FPS which were involved in the upstream pathway of artemisinin biosynthesis showed the highest expression levels in flowers, and the 4 functional genes including ADS, CYP71AV1, CPR and AAR which were involved in the artemisinin-specific biosynthetic pathway were found to be expressed in all the 4 detected tissues. The highest expression level of ADS was found in leaves, then followed by flowers, and the lowest expression level of ADS was found in roots and stems. CYP71AV1 had highest expression level in flowers and lowest in leaves. CPR showed highest expression level in flowers, and AAR had lower expression levels in the other 3 artemisinin-specific pathway genes in all the tissues. The highest content of artemisinin was found in leaves (0.343 mg x g(-1)), then followed by flowers (0.152 mg x g(-1)), roots (0.062 mg x g(-1)) and stems (0.060 mg x g(-1)). CONCLUSION: In the biosynthesis of artemisinin, the upstream genes including HMGR from the MVA pathway, DXR from the MEP pathway and the checkpoint gene FPS were much more active in flowers, and this suggested that flowers might be the tissues of artemisinin precursor biosynthesis, and further DXR contributed more to artemisinin biosynthesis. The positive correlation of ADS expression and artemisinin content in tissues demonstrated that ADS played a very important role in artemisinin biosynthesis, which was the ideal target for engineering the artemisinin biosynthetic pathway. In summary, the functional genes involved in artemisinin biosynthesis do not express at the same level but synergistically.