A Genome-Wide Scenario of Terpene Pathways in Self-pollinated Artemisia annua

Scenarios of genes to metabolites in Artemisia annua remain uninvestigated. Here, we report the use of an integrated approach combining metabolomics, transcriptomics, and gene function analyses to characterize gene-to-terpene and terpene pathway scenarios in a self-pollinating variety of this species. Eighty-eight metabolites including 22 sesquiterpenes (e.g., artemisinin), 26 monoterpenes, two triterpenes, one diterpene and 38 other non-polar metabolites were identified from 14 tissues. These metabolites were differentially produced by leaves and flowers at lower to higher positions. Sequences from cDNA libraries of six tissues were assembled into 18 871 contigs and genome-wide gene expression profiles in tissues were strongly associated with developmental stages and spatial specificities. Sequence mining identified 47 genes that mapped to the artemisinin, non-amorphadiene sesquiterpene, monoterpene, triterpene, 2-C-methyl-D-erythritol 4-phosphate and mevalonate pathways. Pearson correlation analysis resulted in network integration that characterized significant correlations of gene-to-gene expression patterns and gene expression-to-metabolite levels in six tissues simultaneously. More importantly, manipulations of amorpha-4,11-diene synthase gene expression not only affected the activity of this pathway toward artemisinin, artemisinic acid, and arteannuin b but also altered non-amorphadiene sesquiterpene and genome-wide volatile profiles. Such gene-to-terpene landscapes associated with different tissues are fundamental to the metabolic engineering of artemisinin.