Differential transcriptome analysis of glandular and filamentous trichomes in Artemisia annua

Morpho-anatomical characterization and DNA barcoding of Artemesia vulgaris L

Artemisia vulgaris L. belongs to Asteraceae, is a herbal plant that has various benefits in the medical field, so that its use in the medical field can be explored optimally, the plant must be thoroughly identified. This study aims to identify A. vulgaris both in terms of descriptive morpho-anatomy and DNA barcoding using BLAST and phylogenetic tree reconstruction. The morpho-anatomical character was observed on root, stem, and leaf. DNA barcoding analysis was carried out through amplification and alignment of the rbcL and matK genes. All studies were conducted on three samples from Taman Husada (Medicinal Plant Garden) Graha Famili Surabaya, Indonesia. The anatomical slide was prepared by the paraffin method. Morphological studies revealed that the leaves of A. vulgaris both on the lower-middle part and on the upper part of the stem have differences, especially in the character of the stipules, petioles, and incisions they have. Meanwhile, from the study of anatomy, A. vulgaris has an anomocytic type of stomata and its distribution is mostly on the ventral part of the leaves. Through the BLAST process and phylogenetic tree reconstruction, the plant sequences being studied are closely related to several species of the genus Artemisia as indicated by a percentage identity above 98% and branch proximity between taxa in the reconstructed phylogenetic tree.

A Cotyledon-based Virus-Induced Gene Silencing (Cotyledon-VIGS) approach to study specialized metabolism in medicinal plants

BACKGROUND

Virus-induced gene silencing (VIGS) is widely used in plant functional genomics. However, the efficiency of VIGS in young plantlets varies across plant species. Additionally, VIGS is not optimized for many plant species, especially medicinal plants that produce valuable specialized metabolites.

RESULTS

We evaluated the efficacy of five-day-old, etiolated seedlings of Catharanthus roseus (periwinkle) for VIGS. The seedlings were vacuum-infiltrated with Agrobacterium tumefaciens GV3101 cells carrying the tobacco rattle virus (TRV) vectors. The protoporphyrin IX magnesium chelatase subunit H (ChlH) gene, a key gene in chlorophyll biosynthesis, was used as the target for VIGS, and we observed yellow cotyledons 6 days after infiltration. As expected, the expression of CrChlH and the chlorophyll contents of the cotyledons were significantly decreased after VIGS. To validate the cotyledon based-VIGS method, we silenced the genes encoding several transcriptional regulators of the terpenoid indole alkaloid (TIA) biosynthesis in C. roseus, including two activators (CrGATA1 and CrMYC2) and two repressors (CrGBF1 and CrGBF2). Silencing CrGATA1 led to downregulation of the vindoline pathway genes (T3O, T3R, and DAT) and decreased vindoline contents in cotyledons. Silencing CrMYC2, followed by elicitation with methyl jasmonate (MeJA), resulted in the downregulation of ORCA2 and ORCA3. We also co-infiltrated C. roseus seedlings with TRV vectors that silence both CrGBF1 and CrGBF2 and overexpress CrMYC2, aiming to simultaneous silencing two repressors while overexpressing an activator. The simultaneous manipulation of repressors and activator resulted in significant upregulation of the TIA pathway genes. To demonstrate the broad application of the cotyledon-based VIGS method, we optimized the method for two other valuable medicinal plants, Glycyrrhiza inflata (licorice) and Artemisia annua (sweet wormwood). When TRV vectors carrying the fragments of the ChlH genes were infiltrated into the seedlings of these plants, we observed yellow cotyledons with decreased chlorophyll contents.

CONCLUSIONS

The widely applicable cotyledon-based VIGS method is faster, more efficient, and easily accessible to additional treatments than the traditional VIGS method. It can be combined with transient gene overexpression to achieve simultaneous up- and down-regulation of desired genes in non-model plants. This method provides a powerful tool for functional genomics of medicinal plants, facilitating the discovery and production of valuable therapeutic compounds.

Exploring the co-operativity of secretory structures for defense and pollination in flowering plants

MAIN CONCLUSION

In flowers multiple secretory systems cooperate to deliver specialized metabolites to support specific roles in defence and pollination. The collective roles of cell types, enzymes, and transporters are discussed. The interplay between reproductive strategies and defense mechanisms in flowering plants has long been recognized, with trade-offs between investment in defense and reproduction predicted. Glandular trichomes and secretory cavities or ducts, which are epidermal and internal structures, play a pivotal role in the secretion, accumulation, and transport of specialized secondary metabolites, and contribute significantly to defense and pollination. Recent investigations have revealed an intricate connection between these two structures, whereby specialized volatile and non-volatile metabolites are exchanged, collectively shaping their respective ecological functions. However, a comprehensive understanding of this profound integration remains largely elusive. In this review, we explore the secretory systems and associated secondary metabolism primarily in Asteraceous species to propose potential shared mechanisms facilitating the directional translocation of these metabolites to diverse destinations. We summarize recent advances in our understanding of the cooperativity between epidermal and internal secretory structures in the biosynthesis, secretion, accumulation, and emission of terpenes, providing specific well-documented examples from pyrethrum (Tanacetum cinerariifolium). Pyrethrum is renowned for its natural pyrethrin insecticides, which accumulate in the flower head, and more recently, for emitting an aphid alarm pheromone. These examples highlight the diverse specializations of secondary metabolism in pyrethrum and raise intriguing questions regarding the regulation of production and translocation of these compounds within and between its various epidermal and internal secretory systems, spanning multiple tissues, to serve distinct ecological purposes. By discussing the cooperative nature of secretory structures in flowering plants, this review sheds light on the intricate mechanisms underlying the ecological roles of terpenes in defense and pollination.

Trichome-Specific Analysis and Weighted Gene Co-Expression Correlation Network Analysis (WGCNA) Reveal Potential Regulation Mechanism of Artemisinin Biosynthesis in Artemisia annua

Trichomes are attractive cells for terpenoid biosynthesis and accumulation in Artemisia annua. However, the molecular process underlying the trichome of A. annua is not yet fully elucidated. In this study, an analysis of multi-tissue transcriptome data was performed to examine trichome-specific expression patterns. A total of 6646 genes were screened and highly expressed in trichomes, including artemisinin biosynthetic genes such as amorpha-4,11-diene synthase (ADS) and cytochrome P450 monooxygenase (CYP71AV1). Mapman and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that trichome-specific genes were mainly enriched in lipid metabolism and terpenoid metabolism. These trichome-specific genes were analyzed by a weighted gene co-expression network analysis (WGCNA), and the blue module linked to terpenoid backbone biosynthesis was determined. Hub genes correlated with the artemisinin biosynthetic genes were selected based on TOM value. ORA, Benzoate carboxyl methyltransferase (BAMT), Lysine histidine transporter-like 8 (AATL1), Ubiquitin-like protease 1 (Ulp1) and TUBBY were revealed as key hub genes induced by methyl jasmonate (MeJA) for regulating artemisinin biosynthesis. In summary, the identified trichome-specific genes, modules, pathways and hub genes provide clues and shed light on the potential regulatory mechanisms of artemisinin biosynthesis in trichomes in A. annua.

Synergistic interaction of two bHLH transcription factors positively regulates artemisinin biosynthetic pathway in Artemisia annua L

The wonder drug artemisinin, a sesquiterpene lactone endoperoxide from Artemisia annua is the million-dollar molecule required to curb the deadliest disease, Malaria. One of the major challenges even today is to increase the concentration of artemisinin within plants. The transcription factors are important regulators of plant secondary metabolites and have the potential to regulate key steps or the whole biosynthetic pathway. In this study, we have identified and characterised two bHLH transcription factors (Aa6119 and Aa7162) from A. annua. Both the transcription factors turned out to be transcriptionally active and nuclear-localised typical bHLH proteins. In our study, we found that Aa6119 specifically binds to the E-box element present on the promoter of artemisinin biosynthetic gene, AMORPHA-4,11-DIENE SYNTHASE (ADS). The protein-DNA interaction confirmed by Yeast one-hybrid assay was specific as Aa6119 was unable to bind to the mutated E-boxes of ADS. Further, Aa6119 interacted physically with Aa7162, which was confirmed in vitro by Yeast two-hybrid assay and in vivo by Bimolecular Fluorescent complementation assay. Our quantitative expression studies have confirmed that Aa6119 and Aa7162 act synergistically in the regulation of artemisinin biosynthetic and trichome developmental genes. The higher accumulation of artemisinin content in the transient co-transformed transgenic plants than in the individual over-expression transgenic plants has further validated that Aa6119 and Aa7162 act positively and synergistically to regulate artemisinin accumulation.

Xanthanolides in Xanthium L.: Structures, Synthesis and Bioactivity

Xanthanolides were particularly characteristic of the genus Xanthium, which exhibited broad biological effects and have drawn much attention in pharmacological application. The review surveyed the structures and bioactivities of the xanthanolides in the genus Xanthium, and summarized the synthesis tactics of xanthanolides. The results indicated that over 30 naturally occurring xanthanolides have been isolated from the genus Xanthium in monomeric, dimeric and trimeric forms. The bioassay-guided fractionation studies suggested that the effective fractions on antitumor activities were mostly from weak polar solvents, and xanthatin (1) was the most effective and well-studied xanthanolide. The varieties of structures and structure-activity relationships of the xanthanolides had provided the promising skeleton for the further study. The review aimed at providing guidance for the efficient preparation and the potential prospects of the xanthanolides in the medicinal industry.

Allele-aware chromosome-level genome assembly of Artemisia annua reveals the correlation between ADS expansion and artemisinin yield

Artemisia annua is the major natural source of artemisinin, an anti-malarial medicine commonly used worldwide. Here, we present chromosome-level haploid maps for two A. annua strains with different artemisinin contents to explore the relationships between genomic organization and artemisinin production. High-fidelity sequencing, optical mapping, and chromatin conformation capture sequencing were used to assemble the heterogeneous and repetitive genome and resolve the haplotypes of A. annua. Approximately 50,000 genes were annotated for each haplotype genome, and a triplication event that occurred approximately 58.12 million years ago was examined for the first time in this species. A total of 3,903,467-5,193,414 variants (SNPs, indels, and structural variants) were identified in the 1.5-Gb genome during pairwise comparison between haplotypes, consistent with the high heterozygosity of this species. Genomic analyses revealed a correlation between artemisinin concents and the copy number of amorpha-4,11-diene synthase genes. This correlation was further confirmed by resequencing of 36 A. annua samples with varied artemisinin contents. Circular consensus sequencing of transcripts facilitated the detection of paralog expression. Collectively, our study provides chromosome-level allele-aware genome assemblies for two A. annua strains and new insights into the biosynthesis of artemisinin and its regulation, which will contribute to conquering malaria worldwide.

The transcription factors TLR1 and TLR2 negatively regulate trichome density and artemisinin levels in Artemisia annua

The important antimalarial drug artemisinin is biosynthesized and stored in Artemisia annua glandular trichomes and the artemisinin content correlates with trichome density; however, the factors affecting trichome development are largely unknown. Here, we demonstrate that the A. annua R2R3 MYB transcription factor TrichomeLess Regulator 1 (TLR1) negatively regulates trichome development. In A. annua, TLR1 overexpression lines had 44.7%-64.0% lower trichome density and 11.5%-49.4% lower artemisinin contents and TLR1-RNAi lines had 33%-93.3% higher trichome density and 32.2%-84.0% higher artemisinin contents compared with non-transgenic controls. TLR1 also negatively regulates the expression of anthocyanin biosynthetic pathway genes in A. annua. When heterologously expressed in Arabidopsis thaliana, TLR1 interacts with GLABROUS3a, positive regulator of trichome development, and represses trichome development. Yeast two-hybrid and pull-down assays indicated that TLR1 interacts with the WUSCHEL homeobox (WOX) protein AaWOX1, which interacts with the LEAFY-like transcription factor TLR2. TLR2 overexpression in Arabidopsis and A. annua showed that TLR2 reduces trichome development by reducing gibberellin levels. Furthermore, artemisinin contents were 19%-43% lower in TLR2-overexpressing A. annua plants compared to controls. These data indicate that TLR1 and TLR2 negatively regulate trichome density by lowering gibberellin levels and may enable approaches to enhance artemisinin yields.

Transcriptomic and functional analysis provides molecular insights into multicellular trichome development

Trichomes, the hair-like structures located on aerial parts of most vascular plants, are associated with a wide array of biological processes and affect the economic value of certain species. The processes involved in unicellular trichome formation have been well-studied in Arabidopsis (Arabidopsis thaliana). However, our understanding of the morphological changes and the underlying molecular processes involved in multicellular trichome development is limited. Here, we studied the dynamic developmental processes involved in glandular and nonglandular multicellular trichome formation in cucumber (Cucumis sativus L.) and divided these processes into five sequential stages. To gain insights into the underlying mechanisms of multicellular trichome formation, we performed a time-course transcriptome analysis using RNA-sequencing analysis. A total of 711 multicellular trichome-related genes were screened and a model for multicellular trichome formation was developed. The transcriptome and co-expression datasets were validated by reverse transcription-quantitative PCR and in situ hybridization. In addition, virus-induced gene silencing analysis revealed that CsHOMEOBOX3 (CsHOX3) and CsbHLH1 are involved in nonglandular trichome elongation and glandular trichome formation, respectively, which corresponds with the transcriptome data. This study presents a transcriptome atlas that provides insights into the molecular processes involved in multicellular trichome formation in cucumber and can be an important resource for future functional studies.

Coenzymes and the primary and specialized metabolism interface

In plants, primary and specialized metabolism have classically been distinguished as either essential for growth or required for survival in a particular environment. Coenzymes (organic cofactors) are essential for growth but their importance to specialized metabolism is often not considered. In line with the recent proposal of viewing primary and specialized metabolism as an integrated whole rather than segregated lots with a defined interface, we highlight here the importance of collating information on the regulation of coenzyme supply with metabolic demands using examples of vitamin B derived coenzymes. We emphasize that coenzymes can have enormous influence on the outcome of metabolic as well as engineered pathways and should be taken into account in the era of synthetic biology.

Methyl jasmonate treatment, aphid resistance assay, and transcriptomic analysis revealed different herbivore defensive roles between tobacco glandular and non-glandular trichomes

KEY MESSAGE

Methyl jasmonate treatment and aphid resistance assays reveal different roles in herbivore defensive responses between tobacco glandular and non-glandular trichomes. These roles correlate with trichome gene expression patterns. In plants, trichomes greatly contribute to biotic stress resistance. To better understand the different defensive functions between glandular and non-glandular trichomes, we used Nicotiana tabacum as a model. This species bears three types of trichomes: long and short stalk glandular trichomes (LGT and SGT, respectively), and non-glandular trichomes (NGT). Tobacco accession T.I.1068 (lacking NGT) and T.I.1112 (lacking LGT) were used for the experiment. After methyl jasmonate (MeJA) treatment, LGT formation was promoted not only in T.I.1068, but also in T.I.1112, whereas NGT remained absent in T.I.1068, and was slightly reduced in T.I.1112. Diterpenoids, which play important roles in herbivore resistance, accumulated abundantly in T.I.1068 and were elevated by MeJA; however, they were not found in T.I.1112 but became detectable after MeJA treatment. The aphid resistance of T.I.1068 was higher than that of T.I.1112, and both were enhanced by MeJA, which was closely correlated with LGT density. Trichomes detached from T.I.1068 and T.I.1112 were used for RNA-Seq analysis, the results showed that pentose phosphate, photosynthesis, and diterpenoid biosynthesis genes were much more expressed in T.I.1068 than in T.I.1112, which was consistent with the vigorous diterpenoid biosynthesis in T.I.1068. In T.I.1112, citrate cycle, propanoate, and glyoxylate metabolism processes were enriched, and some defensive protein genes were expressed at higher levels than those in T.I.1068.These results suggested that LGT plays a predominant role in aphid resistance, whereas NGT could strengthen herbivore resistance by accumulating defensive proteins, and the roles of LGT and NGT are associated with their gene expression patterns.

Plant Secondary Metabolite Transporters: Diversity, Functionality, and Their Modulation

Secondary metabolites (SMs) play crucial roles in the vital functioning of plants such as growth, development, defense, and survival via their transportation and accumulation at the required site. However, unlike primary metabolites, the transport mechanisms of SMs are not yet well explored. There exists a huge gap between the abundant presence of SM transporters, their identification, and functional characterization. A better understanding of plant SM transporters will surely be a step forward to fulfill the steeply increasing demand for bioactive compounds for the formulation of herbal medicines. Thus, the engineering of transporters by modulating their expression is emerging as the most viable option to achieve the long-term goal of systemic metabolic engineering for enhanced metabolite production at minimum cost. In this review article, we are updating the understanding of recent advancements in the field of plant SM transporters, particularly those discovered in the past two decades. Herein, we provide notable insights about various types of fully or partially characterized transporters from the ABC, MATE, PUP, and NPF families including their diverse functionalities, structural information, potential approaches for their identification and characterization, several regulatory parameters, and their modulation. A novel perspective to the concept of "Transporter Engineering" has also been unveiled by highlighting its potential applications particularly in plant stress (biotic and abiotic) tolerance, SM accumulation, and removal of anti-nutritional compounds, which will be of great value for the crop improvement program. The present study creates a roadmap for easy identification and a better understanding of various transporters, which can be utilized as suitable targets for transporter engineering in future research.

Transcriptome analysis reveals important candidate gene families related to oligosaccharides biosynthesis in Morindaofficinalis

Morinda officinalis How (MO) is one of the best-known traditional herbs and is widely cultivated in subtropical and tropical areas for many years, especially in southern China. Oligosaccharides are the major constituents in the roots of MO, which is well known for its therapeutic effects with anti-depression, anti-osteoporosis, memory-enhancing, ect. To date, the main gene families that regulate the biosynthetic pathway of MO oligosaccharides metabolism yet have been published. In our study, six cDNA libraries generated from six plants of MO were sequenced utilizing an Illumina HiSeq 4000 platform. Corresponding totals of more than 132.60 million clean reads were obtained from the six libraries and assembled into 25,812 unigenes with an average length of 1288 bp. Moreover, 6036 unigenes were found to be allocated to 26 pathways maps using several public databases, and 2538 differential expression genes (DEGs) were screened. Among them, 25 genes from three families were selected as the mainly candidate genes related to MO oligosaccharides biosynthesis. Then, the expression patterns of six DEGs closely related to MO oligosaccharides biosynthesis were verified by quantitative real-time PCR (qRT-PCR). Besides, the MO was clustered more closely to Coffea arabica of Rubiaceae. In summary, the transcriptomic analysis was used to investigate the differences in expression genes of oligosaccharides biosynthesis, with the notable outcome that several key gene families were closely linked to oligosaccharides biosynthesis.

Transcriptome and Small RNA Combined Sequencing Analysis of Cold Tolerance in Non-heading Chinese Cabbage

Background

Non-heading Chinese cabbage (Brassica rapa ssp. chinensis) is an important leaf vegetable grown worldwide. However, there has currently been not enough transcriptome and small RNA combined sequencing analysis of cold tolerance, which hinders further functional genomics research.

Results

In this study, 63.43 Gb of clean data was obtained from the transcriptome analysis. The clean data of each sample reached 6.99 Gb, and the basic percentage of Q30 was 93.68% and above. The clean reads of each sample were sequence aligned with the designated reference genome (Brassica rapa, IVFCAASv1), and the efficiency of the alignment varied from 81.54 to 87.24%. According to the comparison results, 1,860 new genes were discovered in Pak-choi, of which 1,613 were functionally annotated. Among them, 13 common differentially expressed genes were detected in all materials, including seven upregulated and six downregulated. At the same time, we used quantitative real-time PCR to confirm the changes of these gene expression levels. In addition, we sequenced miRNA of the same material. Our findings revealed a total of 34,182,333 small RNA reads, 88,604,604 kinds of small RNAs, among which the most common size was 24 nt. In all materials, the number of common differential miRNAs is eight. According to the corresponding relationship between miRNA and its target genes, we carried out Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis on the set of target genes on each group of differentially expressed miRNAs. Through the analysis, it is found that the distributions of candidate target genes in different materials are different. We not only used transcriptome sequencing and small RNA sequencing but also used experiments to prove the expression levels of differentially expressed genes that were obtained by sequencing. Sequencing combined with experiments proved the mechanism of some differential gene expression levels after low-temperature treatment.

Conclusion

In all, this study provides a resource for genetic and genomic research under abiotic stress in Pak-choi.

Single Wavelengths of LED Light Supplement Promote the Biosynthesis of Major Cyclic Monoterpenes in Japanese Mint

Environmental light conditions influence the biosynthesis of monoterpenes in the mint plant. Cyclic terpenes, such as menthol, menthone, pulegone, and menthofuran, are major odor components synthesized in mint leaves. However, it is unclear how light for cultivation affects the contents of these compounds. Artificial lighting using light-emitting diodes (LEDs) for plant cultivation has the advantage of preferential wavelength control. Here, we monitored monoterpene contents in hydroponically cultivated Japanese mint leaves under blue, red, or far-red wavelengths of LED light supplements. Volatile cyclic monoterpenes, pulegone, menthone, menthol, and menthofuran were quantified using the head-space solid phase microextraction method. As a result, all light wavelengths promoted the biosynthesis of the compounds. Remarkably, two weeks of blue-light supplement increased all compounds: pulegone (362% increase compared to the control), menthofuran (285%), menthone (223%), and menthol (389%). Red light slightly promoted pulegone (256%), menthofuran (178%), and menthol (197%). Interestingly, the accumulation of menthone (229%) or menthofuran (339%) was observed with far-red light treatment. The quantification of glandular trichomes density revealed that no increase under light supplement was confirmed. Blue light treatment even suppressed the glandular trichome formation. No promotion of photosynthesis was observed by pulse-amplitude-modulation (PAM) fluorometry. The present result indicates that light supplements directly promoted the biosynthetic pathways of cyclic monoterpenes.

The WRKY transcription factor AaGSW2 promotes glandular trichome initiation in Artemisia annua

Glandular secreting trichomes (GSTs) synthesize and secrete large quantities of secondary metabolites, some of which have well-established commercial value. An example is the anti-malarial compound artemisinin, which is synthesized in the GSTs of Artemisia annua. Accordingly, there is considerable interest in understanding the processes that regulate GST density as a strategy to increase artemisinin production. In this study, we identified a GST-specific WRKY transcription factor from A. annua, AaGSW2, which is positively regulated by the direct binding of the homeodomain proteins AaHD1 and AaHD8 to the L1-box of the AaGSW2 promoter. Overexpression of AaGSW2 in A. annua significantly increased GST density, while AaGSW2 knockdown lines showed impaired GST initiation. Ectopic expression of AaGSW2 homologs from two mint cultivars, Mentha spicata and Mentha haplocalyx, in A. annua also induced GST formation. These results reveal a molecular mechanism involving homeodomain and WRKY proteins that controls glandular trichome initiation, at least part of which is shared by A. annua and mint.

An R2R3-MYB Transcription Factor Positively Regulates the Glandular Secretory Trichome Initiation in Artemisia annua L

Artemisia annua L. is known for its specific product "artemisinin" which is an active ingredient for curing malaria. Artemisinin is secreted and accumulated in the glandular secretory trichomes (GSTs) on A. annua leaves. Earlier studies have shown that increasing GST density is effective in increasing artemisinin content. However, the mechanism of GST initiation is not fully understood. To this end, we isolated and characterized an R2R3-MYB gene, AaMYB17, which is expressed specifically in the GSTs of shoot tips. Overexpression of AaMYB17 in A. annua increased GST density and enhanced the artemisinin content, whereas RNA interference of AaMYB17 resulted in the reduction of GST density and artemisinin content. Additionally, neither overexpression lines nor RNAi lines showed an abnormal phenotype in plant growth and the morphology of GSTs. Our study demonstrates that AaMYB17 is a positive regulator of GSTs' initiation, without influencing the trichome morphology.

Transcriptome analysis for molecular landscaping of genes controlling diterpene andrographolide biosynthesis in Andrographis paniculata (Burm . f.) Nees

Kalmegh [Andrographis paniculata (Burm. f.) Nees.] is one of the essential medicinal plants due to an important terpenoid, i.e. andrographolide possesses immense therapeutic and pharmacological uses. The experiment was performed to elucidate the expression of candidate genes associated with andrographolide biosynthesis. Based on results obtained in chromatography for andrographolide content analysis of six genotypes, two contrast genotypes, i.e. IC-520361 (maximum andrographolide content-2.33%) and Anand Local (lowest andrographolide content-1.01%) were selected for the transcriptome analysis. A total of 1.04 Gb of raw data were produced using MiSeq Illumina platform, in which IC 520361 generated 645 million base pairs sequence along with 4,524,251 raw reads and Anand Local produced 419 million base pairs sequence along with 3,021,316 raw reads. The combined assembly of high quality reads generated for both the samples had 33,247,454 bp of total assembled bases and 38,292 of transcripts. The GC percent of assembled transcripts was 44.79%, an average read length was 800 bp and N50 value was 1186 bp. Species-specific distribution using BLAST X (Nr), showed the highest Blast hits with Sesamum indicum. Out of 23,346 transcripts, 87% of transcripts annotated in UniProt KB (Universal Protein Resource KnowledgeBase) database and only 0.21% of transcripts were annotated in TAIR (The Arabidopsis Information Resources). Biological processes gene ontology classified based on Blast2GO showed, out of 6853 transcripts, 1370 of transcripts were represented by terpenoid biosynthetic pathway, which involved in secondary metabolite andrographolide biosynthesis. The heat map showed 1016 transcripts were differentially expressed between two kalmegh genotypes, in which nine important differentially expressed transcripts related to MEP (2C methyl-d-erythritol 4-phosphate) and MVA (Mevalonic acid) andrographolide biosynthesis pathways such as, geranyl diphosphate synthase small subunit, Isopentenyl-diphosphate delta-isomerase i-like, 4, 13-hydroxy-3-methylglutaryl-coenzyme a reductase etc. were upregulated in IC 520361 as compared to Anand Local, which were validated through RT-qPCR. The highest expression of gene 13-hydroxy-3-methylglutaryl-coenzyme a reductase (HMGR) was reported, which is responsible for accumulation of andrographolide in leaf. This comparative transcriptome analysis confirmed the expression level of genes were higher in accession IC 520361 as compare to Anand Local related to andrographolide biosynthesis pathways i.e. MEP and MVA. These up-regulated genes could be over-expressed to enhance the andrographolide content using genetic engineering of these metabolic pathways. It will also give an idea to the breeder for development of molecular markers for direct screening of the genotypes.

Comparative transcriptome analysis of root, stem, and leaf tissues of Entada phaseoloides reveals potential genes involved in triterpenoid saponin biosynthesis

Background

Entada phaseoloides (L.) Merr. is an important traditional medicinal plant. The stem of Entada phaseoloides is popularly used as traditional medicine because of its significance in dispelling wind and dampness and remarkable anti-inflammatory activities. Triterpenoid saponins are the major bioactive compounds of Entada phaseoloides. However, genomic or transcriptomic technologies have not been used to study the triterpenoid saponin biosynthetic pathway in this plant.

Results

We performed comparative transcriptome analysis of the root, stem, and leaf tissues of Entada phaseoloides with three independent biological replicates and obtained a total of 53.26 Gb clean data and 116,910 unigenes, with an average N50 length of 1218 bp. Putative functions could be annotated to 42,191 unigenes (36.1%) based on BLASTx searches against the Non-redundant, Uniprot, KEGG, Pfam, GO, KEGG and COG databases. Most of the unigenes related to triterpenoid saponin backbone biosynthesis were specifically upregulated in the stem. A total of 26 cytochrome P450 and 17 uridine diphosphate glycosyltransferase candidate genes related to triterpenoid saponin biosynthesis were identified. The differential expressions of selected genes were further verified by qPT-PCR.

Conclusions

The dataset reported here will facilitate the research about the functional genomics of triterpenoid saponin biosynthesis and genetic engineering of Entada phaseoloides.

Overexpression of blue light receptor AaCRY1 improves artemisinin content in Artemisia annua L.

Artemisinin, an effective antimalarial compound, is isolated from the medicinal plant Artemisia annua L. However, because of the low content of artemisinin in A. annua, the demand of artemisinin exceeds supply. Previous studies show that the artemisinin biosynthesis is promoted by light in A. annua. Cryptochrome1 (CRY1) is involved in many processes in the light response. In this study, AaCRY1 was cloned from A. annua. Overexpressing AaCRY1 in Arabidopsis thaliana cry1 mutant resulted in blue-light-dependent short hypocotyl phenotype and short coleoptile under blue light. Yeast two-hybrid and subcellular colocalization showed that AaCRY1 interacted with AtCOP1 (ubiquitin E3 ligase CONSTITUTIVE PHOTOMORPHOGENIC1). Overexpression of AaCRY1 in transgenic A. annua increased the artemisinin content. When AaCRY1 was overexpressed in A. annua driven by the CYP71AV1 (cytochrome P450 dependent amorpha-4,11-diene 12-hydroxylase) promoter, the artemisinin content was 1.6 times higher than that of the control. Furthermore, we expressed the C terminal of AaCRY1(CCT) involved a GUS-CCT fusion protein in A. annua. The results showed that the artemisinin content was increased to 1.7- to 2.4-fold in GUS-CCT transgenic A. annua plants. These results demonstrate that overexpression of GUS-CCT is an effective strategy to increase artemisinin production in A. annua.

Cross-Species Annotation of Expressed Genes and Detection of Different Functional Gene Modules Between 10 Cold- and 10 Hot-Propertied Chinese Herbal Medicines

According to the traditional Chinese medicine (TCM) system, Chinese herbal medicines (HMs) can be divided into four categories: hot, warm, cold, and cool. A cool nature usually is categorized as a cold nature, and a warm nature is classified as a hot nature. However, the detectable characteristics of the gene expression profile associated with the cold and hot properties have not been studied. To address this question, a strategy for the cross-species annotation of conserved genes was established in the present study by using transcriptome data of 20 HMs with cold and hot properties. Functional enrichment analysis was performed on group-specific expressed genes inferred from the functional genome of the reference species (i.e., Arabidopsis). Results showed that metabolic pathways relevant to chrysoeriol, luteolin, paniculatin, and wogonin were enriched for cold-specific genes, and pathways of inositol, heptadecane, lauric acid, octanoic acid, hexadecanoic acid, and pentadecanoic acid were enriched for hot-specific genes. Six functional modules were identified in the HMs with the cold property: nucleotide biosynthetic process, peptidy-L-cysteine S-palmitoylation, lipid modification, base-excision repair, dipeptide transport, and response to endoplasmic reticulum stress. For the hot HMs, another six functional modules were identified: embryonic meristem development, embryonic pattern specification, axis specification, regulation of RNA polymerase II transcriptional preinitiation complex assembly, mitochondrial RNA modification, and cell redox homeostasis. The research provided a new insight into HMs’ cold and hot properties from the perspective of the gene expression profile of plants.

Genetic Control of Glandular Trichome Development

Plant glandular trichomes are epidermal secretory structures producing various specialized metabolites. These metabolites are involved in plant adaptation to its environment and many of them have remarkable properties exploited by fragrance, flavor, and pharmaceutical industries. The identification of genes controlling glandular trichome development is of high interest to understand how plants produce specialized metabolites. Our knowledge about this developmental process is still limited, but genes controlling glandular trichome initiation and morphogenesis have recently been identified. In particular, R2R3-MYB and HD-ZIP IV transcription factors appear to play essential roles in glandular trichome initiation in Artemisia annua and tomato. In this review, we focus on the results obtained in these two species and we propose genetic regulation models integrating these data.

An oxidatively stressful situation: a case of Artemisia annua L

Artemisinin (ART) is an antimalarial compound that possesses a variety of novel biological activities. Due to the low abundance of ART in natural sources, agricultural supply has been erratic, and prices are highly volatile. While heterologous biosynthesis and semi-synthesis are advantageous in certain aspects, these approaches remained disadvantageous in terms of productivity and cost-effectiveness. Therefore, further improvement in ART production calls for approaches that should supplement the agricultural production gap, while reducing production costs and stabilising supply. The present review offers a discussion on the elicitation of plants and/or in vitro cultures as an economically feasible yield enhancement strategy to address the global problem of access to affordable ART. Deemed critical for the manipulation of biosynthetic potential, the mechanism of ART biosynthesis is reviewed. It includes a discussion on the current biotechnological solutions to ART production, focusing on semi-synthesis and elicitation. A brief commentary on the possible aspects that influence elicitation efficiency and how oxidative stress modulates ART synthesis is also presented. Based on the critical analysis of current literature, a hypothesis is put forward to explain the possible involvement of enzymes in assisting the final non-enzymatic transformation step leading to ART formation. This review highlights the critical factors limiting the success of elicitor-induced modulation of ART metabolism, that will help inform strategies for future improvement of ART production. Additionally, new avenues for future research based on the proposed hypothesis will lead to exciting perspectives in this research area and continue to enhance our understanding of this intricate metabolic process.

Comprehensive transcriptome analysis provides insights into metabolic and gene regulatory networks in trichomes of Nicotiana tabacum

KEY MESSAGE

Comprehensive transcriptome analysis suggested that the primary metabolism is modulated to augment the supply of substrates towards secondary metabolism operating in the glandular trichomes of Nicotiana tabacum. The comparative gene expression and co-expression network analysis revealed that certain members of transcription factor genes belonging to the MYB, HD-ZIP, ERF, TCP, SRS, WRKY and DOF families may be involved in the regulation of metabolism and/other aspects in the glandular trichomes of N. tabacum The glandular trichomes of Nicotiana tabacum are highly productive in terms of secondary metabolites and therefore have been projected to be used as a prognostic platform for metabolic engineering of valuable natural products. For obvious reasons, detailed studies pertaining to the metabolic and gene regulatory networks operating in the glandular trichomes of N. tabacum are of pivotal significance to be undertaken. We have carried out next-generation sequencing of glandular trichomes of N. tabcaum and investigated differential gene expression among different tissues, including trichome-free leaves. We identified a total of 37,269 and 37,371 genes, expressing in trichome free leaf and glandular trichomes, respectively, at a cutoff of FPKM ≥ 1. The analysis revealed that different pathways involved with the primary metabolism are modulated in glandular trichomes of N. tabacum, providing a plausible explanation for the enhanced biosynthesis of secondary metabolism in the glandular trichomes. Further, comparative gene expression analysis revealed several genes, which display preferential expression in the glandular trichomes and thereby seem to be potential candidate genes for future studies in connection to the discovery of novel trichome specific promoters. The present study also led to the comprehensive identification of 1750 transcription factor genes expressing at a cutoff of FPKM ≥ 1 in the glandular trichomes of N. tabacum. The clustering and co-expression analysis suggested that transcription factor genes belonging to HD-ZIP, ERF, WRKY, MYB, TCP, SRS and DOF families may be the major players in the regulation of gene expression in the glandular trichomes of N. tabacum. To the best of our knowledge, the present work is the first effort towards detailed identification of genes, especially regulatory genes expressing in the glandular trichomes of N. tabacum. The data resource and the empirical findings from present work in all probability must, therefore, provide a reference and background context for future work aiming at deciphering molecular mechanism of regulation of secondary metabolism and gene expression in the glandular trichomes of N. tabacum.

Integrative Transcriptomic and Metabolic Analyses Provide Insights into the Role of Trichomes in Tea Plant (Camellia Sinensis)

Trichomes, which develop from epidermal cells, are regarded as one of the key features that are involved in the evaluation of tea quality and tea germplasm resources. The metabolites from trichomes have been well characterized in tea products. However, little is known regarding the metabolites in fresh tea trichomes and the molecular differences in trichomes and tea leaves per se. In this study, we developed a method to collect trichomes from tea plant tender shoots, and their main secondary metabolites, including catechins, caffeine, amino acids, and aroma compounds, were determined. We found that the majority of these compounds were significantly less abundant in trichomes than in tea leaves. RNA-Seq was used to investigate the differences in the molecular regulatory mechanism between trichomes and leaves to gain further insight into the differences in trichomes and tea leaves. In total, 52.96 Gb of clean data were generated, and 6560 differentially expressed genes (DEGs), including 4471 upregulated and 2089 downregulated genes, were identified in the trichomes vs. leaves comparison. Notably, the structural genes of the major metabolite biosynthesis pathways, transcription factors, and other key DEGs were identified and comparatively analyzed between trichomes and leaves, while trichome-specific genes were also identified. Our results provide new insights into the differences between tea trichomes and leaves at the metabolic and transcriptomic levels, and open up new doors to further recognize and re-evaluate the role of trichomes in tea quality formation and tea plant growth and development.

Synthesis of [3,3-2H2]-Dihydroartemisinic Acid to Measure the Rate of Nonenzymatic Conversion of Dihydroartemisinic Acid to Artemisinin

Dihydroartemisinic acid is the biosynthetic precursor to artemisinin, the endoperoxide-containing natural product used to treat malaria. The conversion of dihydroartemisinic acid to artemisinin is a cascade reaction that involves C-C bond cleavage, hydroperoxide incorporation, and polycyclization to form the endoperoxide. Whether or not this reaction is enzymatically controlled has been controversial. A method was developed to quantify the nonenzymatic conversion of dihydroartemisinic acid to artemisinin using LC-MS. A seven-step synthesis of 3,3-dideuterodihydroartemisinic acid (23) was accomplished beginning with dihydroartemisinic acid (1). The nonenzymatic rates of formation of 3,3-dideuteroartemisinin (24) from 3,3-dideuterodihydroartemisinic acid (23) were 1400 ng/day with light and 32 ng/day without light. Moreover, an unexpected formation of nondeuterated artemisinin (3) from 3,3-dideuterodihydroartemisinic acid (23) was detected in both the presence and absence of light. This formation of nondeuterated artemisinin (3) from its dideuterated precursor (23) suggests an alternative mechanistic pathway that operates independent of light to form artemisinin, involving the loss of the two C-3 deuterium atoms.

Flavonoid Versus Artemisinin Anti-malarial Activity in Artemisia annua Whole-Leaf Extracts

Artemisinin, a sesquiterpene lactone produced by Artemisia annua glandular secretory trichomes, is the active ingredient in the most effective treatment for uncomplicated malaria caused by Plasmodium falciparum parasites. Other metabolites in A. annua or related species, particularly flavonoids, have been proposed to either act as antimalarials on their own or act synergistically with artemisinin to enhance antimalarial activity. We identified a mutation that disrupts the CHALCONE ISOMERASE 1 (CHI1) enzyme that is responsible for the second committed step of flavonoid biosynthesis. Detailed metabolite profiling revealed that chi1-1 lacks all major flavonoids but produces wild-type artemisinin levels, making this mutant a useful tool to test the antiplasmodial effects of flavonoids. We used whole-leaf extracts from chi1-1 and mutant lines impaired in artemisinin production in bioactivity in vitro assays against intraerythrocytic P. falciparum Dd2. We found that chi1-1 extracts did not differ from wild-type extracts in antiplasmodial efficacy nor initial rate of cytocidal action. Furthermore, extracts from the A. annua cyp71av1-1 mutant and RNAi lines impaired in amorpha-4,11-diene synthase gene expression, which are both severely compromised in artemisinin biosynthesis but unaffected in flavonoid metabolism, showed very low or no antiplasmodial activity. These results demonstrate that in vitro bioactivity against P. falciparum of flavonoids is negligible when compared to that of artemisinin.

Transcriptome analysis and identification of genes related to terpenoid biosynthesis in Cinnamomum camphora

BACKGROUND

Cinnamomum camphora has been cultivated as an economically important tree for its medicinal and aromatic properties. Selective breeding has produced Cinnamomum plants for special uses, including spice strains with characteristic flavors and aromas and high-potency medicinal cultivars. The molecular biology underlying terpenoid biosynthesis is still unexplored.

RESULTS

Gas chromatography-mass spectrometry was used to analyze the differences in contents and compositions of essential oil terpenoids in linalool- and borneol-type chemotypes of C. camphora. The data revealed that the essential oils consist primarily of monoterpenes with only very minor quantities of sesquiterpenes and diterpenes and that the essential oil differs in different chemotypes of C. camphora, with higher yields of (-)-borneol from the borneol-type than from the linalool-type. To study the terpenoid biosynthesis of signature compounds of the major monoterpenes, we performed RNA sequencing to profile the leaf transcriptomes of the two chemotypes of C. camphora. A total of 23.76 Gb clean data was generated from two chemotypes and assembled into 156,184 unigenes. The total length, average length, N50 and GC content of unigenes were 155,645,929 bp, 997 bp, 1430 bp, and 46.5%, respectively. Among them, 76,421 unigenes were annotated by publicly available databases, of which 67 candidate unigenes were identified to be involved in terpenoid biosynthesis in C. camphora. A total of 2863 unigenes were identified to be differentially expression between borneol-type and linalool-type, including 1714 up-regulated and 1149 down-regulated unigenes. Most genes encoding proteins involved in terpenoid precursor MVA and MEP pathways were expressed in similar levels in both chemotypes of C. camphora. In addition, 10 and 17 DEGs were significantly enriched in the terpene synthase activity and oxidoreductase activity terms of their directed acyclic graphs (DAG), respectively. Three monoterpene synthase genes, TPS14-like1, TPS14-like2 and TPS14-like3 were up-regulated in the borneol-type compared to the linalool-type, and their expression levels were further verified using quantitative real-time PCR.

CONCLUSIONS

This study provides a global overview of gene expression patterns related to terpenoid biosynthesis in C. camphora, and could contribute to a better understanding of the differential accumulation of terpenoids in different C. camphora chemotypes.

The Genome of Artemisia annua Provides Insight into the Evolution of Asteraceae Family and Artemisinin Biosynthesis

Artemisia annua, commonly known as sweet wormwood or Qinghao, is a shrub native to China and has long been used for medicinal purposes. A. annua is now cultivated globally as the only natural source of a potent anti-malarial compound, artemisinin. Here, we report a high-quality draft assembly of the 1.74-gigabase genome of A. annua, which is highly heterozygous, rich in repetitive sequences, and contains 63 226 protein-coding genes, one of the largest numbers among the sequenced plant species. We found that, as one of a few sequenced genomes in the Asteraceae, the A. annua genome contains a large number of genes specific to this large angiosperm clade. Notably, the expansion and functional diversification of genes encoding enzymes involved in terpene biosynthesis are consistent with the evolution of the artemisinin biosynthetic pathway. We further revealed by transcriptome profiling that A. annua has evolved the sophisticated transcriptional regulatory networks underlying artemisinin biosynthesis. Based on comprehensive genomic and transcriptomic analyses we generated transgenic A. annua lines producing high levels of artemisinin, which are now ready for large-scale production and thereby will help meet the challenge of increasing global demand of artemisinin.

De novo assembly and analysis of the Artemisia argyi transcriptome and identification of genes involved in terpenoid biosynthesis

Artemisia argyi Lev. et Vant. (A. argyi) is widely utilized for moxibustion in Chinese medicine, and the mechanism underlying terpenoid biosynthesis in its leaves is suggested to play an important role in its medicinal use. However, the A. argyi transcriptome has not been sequenced. Herein, we performed RNA sequencing for A. argyi leaf, root and stem tissues to identify as many as possible of the transcribed genes. In total, 99,807 unigenes were assembled by analysing the expression profiles generated from the three tissue types, and 67,446 of those unigenes were annotated in public databases. We further performed differential gene expression analysis to compare leaf tissue with the other two tissue types and identified numerous genes that were specifically expressed or up-regulated in leaf tissue. Specifically, we identified multiple genes encoding significant enzymes or transcription factors related to terpenoid synthesis. This study serves as a valuable resource for transcriptome information, as many transcribed genes related to terpenoid biosynthesis were identified in the A. argyi transcriptome, providing a functional genomic basis for additional studies on molecular mechanisms underlying the medicinal use of A. argyi.

Detailed Phytochemical Analysis of High- and Low Artemisinin-Producing Chemotypes of Artemisia annua

Chemical derivatives of artemisinin, a sesquiterpene lactone produced by Artemisia annua, are the active ingredient in the most effective treatment for malaria. Comprehensive phytochemical analysis of two contrasting chemotypes of A. annua resulted in the characterization of over 80 natural products by NMR, more than 20 of which are novel and described here for the first time. Analysis of high- and low-artemisinin producing (HAP and LAP) chemotypes of A. annua confirmed the latter to have a low level of DBR2 (artemisinic aldehyde Δ¹¹⁽¹³⁾ reductase) gene expression. Here we show that the LAP chemotype accumulates high levels of artemisinic acid, arteannuin B, epi-deoxyarteannuin B and other amorpha-4,11-diene derived sesquiterpenes which are unsaturated at the 11,13-position. By contrast, the HAP chemotype is rich in sesquiterpenes saturated at the 11,13-position (dihydroartemisinic acid, artemisinin and dihydro-epi-deoxyarteannunin B), which is consistent with higher expression levels of DBR2, and also with the presence of a HAP-chemotype version of CYP71AV1 (amorpha-4,11-diene C-12 oxidase). Our results indicate that the conversion steps from artemisinic acid to arteannuin B, epi-deoxyarteannuin B and artemisitene in the LAP chemotype are non-enzymatic and parallel the non-enzymatic conversion of DHAA to artemisinin and dihyro-epi-deoxyarteannuin B in the HAP chemotype. Interestingly, artemisinic acid in the LAP chemotype preferentially converts to arteannuin B rather than the endoperoxide bridge containing artemisitene. In contrast, in the HAP chemotype, DHAA preferentially converts to artemisinin. Broader metabolomic and transcriptomic profiling revealed significantly different terpenoid profiles and related terpenoid gene expression in these two morphologically distinct chemotypes.

Silencing amorpha-4,11-diene synthase Genes in Artemisia annua Leads to FPP Accumulation

Artemisia annua is established as an efficient crop for the production of the anti-malarial compound artemisinin, a sesquiterpene lactone synthesized and stored in Glandular Secretory Trichomes (GSTs) located on the leaves and inflorescences. Amorpha-4,11-diene synthase (AMS) catalyzes the conversion of farnesyl pyrophosphate (FPP) to amorpha-4,11-diene and diphosphate, which is the first committed step in the synthesis of artemisinin. FPP is the precursor for sesquiterpene and sterol biosynthesis in the plant. This work aimed to investigate the effect of blocking the synthesis of artemisinin in the GSTs of a high artemisinin yielding line, Artemis, by down regulating AMS. We determined that there are up to 12 AMS gene copies in Artemis, all expressed in GSTs. We used sequence homology to design an RNAi construct under the control of a GST specific promoter that was predicted to be effective against all 12 of these genes. Stable transformation of Artemis with this construct resulted in over 95% reduction in the content of artemisinin and related products, and a significant increase in the FPP pool. The Artemis AMS silenced lines showed no morphological alterations, and metabolomic and gene expression analysis did not detect any changes in the levels of other major sesquiterpene compounds or sesquiterpene synthase genes in leaf material. FPP also acts as a precursor for squalene and sterol biosynthesis but levels of these compounds were also not altered in the AMS silenced lines. Four unknown oxygenated sesquiterpenes were produced in these lines, but at extremely low levels compared to Artemis non-transformed controls (NTC). This study finds that engineering A. annua GSTs in an Artemis background results in endogenous terpenes related to artemisinin being depleted with the precursor FPP actually accumulating rather than being utilized by other endogenous enzymes. The challenge now is to establish if this precursor pool can act as substrate for production of alternative sesquiterpenes in A. annua.

Plant Glandular Trichomes: Natural Cell Factories of High Biotechnological Interest

Multicellular glandular trichomes are epidermal outgrowths characterized by the presence of a head made of cells that have the ability to secrete or store large quantities of specialized metabolites. Our understanding of the transcriptional control of glandular trichome initiation and development is still in its infancy. This review points to some central questions that need to be addressed to better understand how such specialized cell structures arise from the plant protodermis. A key and unique feature of glandular trichomes is their ability to synthesize and secrete large amounts, relative to their size, of a limited number of metabolites. As such, they qualify as true cell factories, making them interesting targets for metabolic engineering. In this review, recent advances regarding terpene metabolic engineering are highlighted, with a special focus on tobacco (Nicotiana tabacum). In particular, the choice of transcriptional promoters to drive transgene expression and the best ways to sink existing pools of terpene precursors are discussed. The bioavailability of existing pools of natural precursor molecules is a key parameter and is controlled by so-called cross talk between different biosynthetic pathways. As highlighted in this review, the exact nature and extent of such cross talk are only partially understood at present. In the future, awareness of, and detailed knowledge on, the biology of plant glandular trichome development and metabolism will generate new leads to tap the largely unexploited potential of glandular trichomes in plant resistance to pests and lead to the improved production of specialized metabolites with high industrial or pharmacological value.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of the Medicinal Plant Artemisia annua

The complete chloroplast genome of Artemisia annua (Asteraceae), the primary source of artemisinin, was sequenced and analyzed. The A. annua cp genome is 150,995 bp, and harbors a pair of inverted repeat regions (IRa and IRb), of 24,850 bp each that separate large (LSC, 82,988 bp) and small (SSC, 18,267 bp) single-copy regions. Our annotation revealed that the A. annua cp genome contains 113 genes and 18 duplicated genes. The gene order in the SSC region of A. annua is inverted; this fact is consistent with the sequences of chloroplast genomes from three other Artemisia species. Fifteen (15) forward and seventeen (17) inverted repeats were detected in the genome. The existence of rich SSR loci in the genome suggests opportunities for future population genetics work on this anti-malarial medicinal plant. In A. annua cpDNA, the rps19 gene was found in the LSC region rather than the IR region, and the rps19 pseudogene was absent in the IR region. Sequence divergence analysis of five Asteraceae species indicated that the most highly divergent regions were found in the intergenic spacers, and that the differences between A. annua and A. fukudo were very slight. A phylogenetic analysis revealed a sister relationship between A. annua and A. fukudo. This study identified the unique characteristics of the A. annua cp genome. These results offer valuable information for future research on Artemisia species identification and for the selective breeding of A. annua with high pharmaceutical efficacy.

Comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri identifies potential genes related to triterpenoid saponin biosynthesis

Background

Bacopa monnieri commonly known as Brahmi is utilized in Ayurveda to improve memory and many other human health benefits. Bacosides enriched standardized extract of Bacopa monnieri is being marketed as a memory enhancing agent. In spite of its well known pharmacological properties it is not much studied in terms of transcripts involved in biosynthetic pathway and its regulation that controls the secondary metabolic pathway in this plant. The aim of this study was to identify the potential transcripts and provide a framework of identified transcripts involved in bacosides production through transcriptome assembly.

Results

We performed comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri in two independent biological replicate and obtained 22.48 million and 22.0 million high quality processed reads in shoot and root respectively. After de novo assembly and quantitative assessment total 26,412 genes got annotated in root and 18,500 genes annotated in shoot sample. Quality of raw reads was determined by using SeqQC-V2.2. Assembled sequences were annotated using BLASTX against public database such as NR or UniProt. Searching against the KEGG pathway database indicated that 37,918 unigenes from root and 35,130 unigenes from shoot were mapped to 133 KEGG pathways. Based on the DGE data we found that most of the transcript related to CYP450s and UDP-glucosyltransferases were specifically upregulated in shoot tissue as compared to root tissue. Finally, we have selected 43 transcripts related to secondary metabolism including transcription factor families which are differentially expressed in shoot and root tissues were validated by qRT-PCR and their expression level were monitored after MeJA treatment and wounding for 1, 3 and 5 h.

Conclusions

This study not only represents the first de novo transcriptome analysis of Bacopa monnieri but also provides information about the identification, expression and differential tissues specific distribution of transcripts related to triterpenoid sapogenin which is one of the most important pharmacologically active secondary metabolite present in Bacopa monnieri. The identified transcripts in this study will establish a foundation for future studies related to carrying out the metabolic engineering for increasing the bacosides biosynthesis and its regulation for human health benefits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3865-5) contains supplementary material, which is available to authorized users.

Multi-Omics of Tomato Glandular Trichomes Reveals Distinct Features of Central Carbon Metabolism Supporting High Productivity of Specialized Metabolites

Glandular trichomes are metabolic cell factories with the capacity to produce large quantities of secondary metabolites. Little is known about the connection between central carbon metabolism and metabolic productivity for secondary metabolites in glandular trichomes. To address this gap in our knowledge, we performed comparative metabolomics, transcriptomics, proteomics, and ¹³C-labeling of type VI glandular trichomes and leaves from a cultivated (Solanum lycopersicum LA4024) and a wild (Solanum habrochaites LA1777) tomato accession. Specific features of glandular trichomes that drive the formation of secondary metabolites could be identified. Tomato type VI trichomes are photosynthetic but acquire their carbon essentially from leaf sucrose. The energy and reducing power from photosynthesis are used to support the biosynthesis of secondary metabolites, while the comparatively reduced Calvin-Benson-Bassham cycle activity may be involved in recycling metabolic CO₂ Glandular trichomes cope with oxidative stress by producing high levels of polyunsaturated fatty acids, oxylipins, and glutathione. Finally, distinct mechanisms are present in glandular trichomes to increase the supply of precursors for the isoprenoid pathways. Particularly, the citrate-malate shuttle supplies cytosolic acetyl-CoA and plastidic glycolysis and malic enzyme support the formation of plastidic pyruvate. A model is proposed on how glandular trichomes achieve high metabolic productivity.

The 2-oxoglutarate-dependent dioxygenase superfamily participates in tanshinone production in Salvia miltiorrhiza

Highly oxidized tanshinones are pharmacological ingredients extracted from the medicinal model plant Salvia miltiorrhiza and are mainly used to treat cardiovascular diseases. Previous studies have confirmed that cytochrome P450 mono-oxygenases (CYP450s) have a key function in the biosynthesis of tanshinones; however, no solid evidence links oxidation to the 2-oxoglutarate-dependent dioxygenase (2OGD) superfamily. Here, we identified 132 members of the DOXB and DOXC subfamilies of 2OGD by scanning the 2OG-FeII Oxy domain using a genome-wide strategy in S. miltiorrhiza. The DOXC class was phylogenetically divided into twelve clades. Combining phylogenetic relationships, differential expression and co-expression from various organs and tissues revealed that two 2OGDs were directly related to flavonoid metabolism, and that 13 2OGDs from different clades were predicted to be involved in tanshinone biosynthesis. Based on this insight into tanshinone production, we experimentally detected significant decreases in miltirone, cryptotanshinone, and tanshinone IIA (0.16-, 0.56-, and 0.56-fold, respectively) in 2OGD5 RNAi transgenic lines relative to the control lines using a metabonomics analysis. 2OGD5 was found to play a crucial role in the downstream biosynthesis of tanshinones following the hydroxylation of CYPs. Our results highlight the evolution and diversification of 2OGD superfamily members and suggest that they contribute to the complexity of tanshinone metabolites.

Making Plants Break a Sweat: the Structure, Function, and Evolution of Plant Salt Glands

Salt stress is a complex trait that poses a grand challenge in developing new crops better adapted to saline environments. Some plants, called recretohalophytes, that have naturally evolved to secrete excess salts through salt glands, offer an underexplored genetic resource for examining how plant development, anatomy, and physiology integrate to prevent excess salt from building up to toxic levels in plant tissue. In this review we examine the structure and evolution of salt glands, salt gland-specific gene expression, and the possibility that all salt glands have originated via evolutionary modifications of trichomes. Salt secretion via salt glands is found in more than 50 species in 14 angiosperm families distributed in caryophyllales, asterids, rosids, and grasses. The salt glands of these distantly related clades can be grouped into four structural classes. Although salt glands appear to have originated independently at least 12 times, they share convergently evolved features that facilitate salt compartmentalization and excretion. We review the structural diversity and evolution of salt glands, major transporters and proteins associated with salt transport and secretion in halophytes, salt gland relevant gene expression regulation, and the prospect for using new genomic and transcriptomic tools in combination with information from model organisms to better understand how salt glands contribute to salt tolerance. Finally, we consider the prospects for using this knowledge to engineer salt glands to increase salt tolerance in model species, and ultimately in crops.

AaPDR3, a PDR Transporter 3, Is Involved in Sesquiterpene β-Caryophyllene Transport in Artemisia annua

Artemisinin, a sesquiterpenoid endoperoxide, isolated from the plant Artemisia annua L., is widely used in the treatment of malaria. Another sesquiterpenoid, β-caryophyllene having antibiotic, antioxidant, anticarcinogenic and local anesthetic activities, is also presented in A. annua. The role played by sesquiterpene transporters in trichomes and accumulation of these metabolites is poorly understood in A. annua and in trichomes of other plant species. We identified AaPDR3, encoding a pleiotropic drug resistance (PDR) transporter located to the plasma membrane from A. annua. Expression of AaPDR3 is tissue-specifically and developmentally regulated in A. annua. GUS activity is primarily restricted to T-shaped trichomes of old leaves and roots of transgenic A. annua plants expressing proAaPDR3: GUS. The level of β-caryophyllene was decreased in transgenic A. annua plants expressing AaPDR3-RNAi while transgenic A. annua plants expressing increased levels of AaPDR3 accumulated higher levels of β-caryophyllene. When AaPDR3 was expressed in transformed yeast, yeasts expressing AaPDR3 accumulated more β-caryophyllene, rather than germacrene D and β-farnesene, compared to the non-expressing control.

Transcriptome Analysis of Genes Associated with the Artemisinin Biosynthesis by Jasmonic Acid Treatment under the Light in Artemisia annua

Artemisinin is a sesquiterpene lactone endoperoxide extracted from a traditional Chinese medicinal plant Artemisia annua. Artemisinin-based combination therapies (ACTs) are recommended as the best treatment of malaria by the World Health Organization (WHO). Both the phytohormone jasmonic acid (JA) and light promote artemisinin biosynthesis in A. annua. Interestingly, we found that the increase of artemisinin biosynthesis by JA was dependent on light. However, the relationship between the two signal pathways mediated by JA and light remains unclear. Here, we collected the A. annua seedlings of 24 h continuous light (Light), 24 h dark treatment (Dark), 4 h MeJA treatment under the continuous light conditions (Light-MeJA-4h) and 4 h MeJA treatment under the dark conditions (Dark-MeJA-4h) and performed the transcriptome sequencing using Illumina HiSeq 4000 System. A total of 266.7 million clean data were produced and assembled into 185,653 unigenes, with an average length of 537 bp. Among them, 59,490 unigenes were annotated and classified based on the public information. Differential expression analyses were performed between Light and Dark, Light and Light-MeJA-4h, Dark and Dark-MeJA-4h, Light-MeJA-4h, and Dark-MeJA-4h, respectively. Furthermore, transcription factor (TF) analysis revealed that 1588 TFs were identified and divided into 55 TF families, with 284 TFs down-regulated in the Dark relative to Light and 96 TFs up-regulated in the Light-MeJA-4h relative to Light. 8 TFs were selected as candidates for regulating the artemisinin biosynthesis and one of them was validated to be involved in artemisinin transcriptional regulation by Dual-Luciferase (Dual-LUC) assay. The transcriptome data shown in our study offered a comprehensive transcriptional expression pattern influenced by the MeJA and light in A. annua seedling, which will serve as a valuable resource for further studies on transcriptional regulation mechanisms underlying artemisinin biosynthesis.

Artemisia annua mutant impaired in artemisinin synthesis demonstrates importance of nonenzymatic conversion in terpenoid metabolism

Artemisinin, a sesquiterpene lactone produced by Artemisia annua glandular secretory trichomes, is the active ingredient in the most effective treatment for malaria currently available. We identified a mutation that disrupts the amorpha-4,11-diene C-12 oxidase (CYP71AV1) enzyme, responsible for a series of oxidation reactions in the artemisinin biosynthetic pathway. Detailed metabolic studies of cyp71av1-1 revealed that the consequence of blocking the artemisinin biosynthetic pathway is the redirection of sesquiterpene metabolism to a sesquiterpene epoxide, which we designate arteannuin X. This sesquiterpene approaches half the concentration observed for artemisinin in wild-type plants, demonstrating high-flux plasticity in A. annua glandular trichomes and their potential as factories for the production of novel alternate sesquiterpenes at commercially viable levels. Detailed metabolite profiling of leaf maturation time-series and precursor-feeding experiments revealed that nonenzymatic conversion steps are central to both artemisinin and arteannuin X biosynthesis. In particular, feeding studies using ¹³C-labeled dihydroartemisinic acid (DHAA) provided strong evidence that the final steps in the synthesis of artemisinin are nonenzymatic in vivo. Our findings also suggest that the specialized subapical cavity of glandular secretory trichomes functions as a location for both the chemical conversion and the storage of phytotoxic compounds, including artemisinin. We conclude that metabolic engineering to produce high yields of novel secondary compounds such as sesquiterpenes is feasible in complex glandular trichomes. Such systems offer advantages over single-cell microbial hosts for production of toxic natural products.

Root regulation of artemisinin production in Artemisia annua: trichome and metabolite evidence

Roots of plants with high artemisinin-producing leaves increased leaf production of artemisinin in low-producing plants and vice versa indicating roots are involved in controlling artemisinin biosynthesis in shoots. The anti-malarial sesquiterpene, artemisinin, is produced and stored in glandular trichomes (GLTs) of Artemisia annua. Evidence suggested roots, which produce no significant artemisinin nor precursor compounds, regulate production of artemisinin biosynthesis in the leaves. Using grafting, we studied the relationship between rootstock and scion by measuring GLTs and five artemisinic metabolites (artemisinin, deoxyartemisinin, dihydroartemisinic acid, artemisinic acid, arteannuin B) in scions of ungrafted, self-grafted, and cross-grafted plants among three cultivars: S and 15 both having GLTs with artemisinin at 1.49 and 0.57 %, respectively, and G producing neither GLTs nor detectable artemisinin. All artemisinin-producing self-grafts, e.g., S/S (scion/rootstock) and 15/15, produced more artemisinin than ungrafted plants, likely from grafting stress. S/S grafts also produced more GLTs. The 15/S grafts produced more artemisinin than S/15, suggesting rootstocks from high producing S plants stimulated artemisinin production in 15 scions. S/15 grafts yielded less artemisinin than S/S, but more than either 15/15 or ungrafted n15 and nS; S/15 grafts also had a lower density of GLTs than S/S, suggesting rootstock inhibition of the scion. The S rootstock induced trace artemisinin production in G scions, but did not induce GLT formation in G/S grafts. Different grafts exhibited different trichome morphologies and effects on artemisinic pathway flux. This study provides new information regarding the role of roots in GLT development and artemisinin production in this important medicinal plant.

Transient production of artemisinin in Nicotiana benthamiana is boosted by a specific lipid transfer protein from A. annua

Our lack of full understanding of transport and sequestration of the heterologous products currently limit metabolic engineering in plants for the production of high value terpenes. For instance, although all genes of the artemisinin/arteannuin B (AN/AB) biosynthesis pathway (AN-PW) from Artemisia annua have been identified, ectopic expression of these genes in Nicotiana benthamiana yielded mostly glycosylated pathway intermediates and only very little free (dihydro)artemisinic acid [(DH)AA]. Here we demonstrate that Lipid Transfer Protein 3 (AaLTP3) and the transporter Pleiotropic Drug Resistance 2 (AaPDR2) from A. annua enhance accumulation of (DH)AA in the apoplast of N. benthamiana leaves. Analysis of apoplast and cell content and apoplast exclusion assays show that AaLTP3 and AaPDR2 prevent reflux of (DH)AA from the apoplast back into the cells and enhances overall flux through the pathway. Moreover, AaLTP3 is stabilized in the presence of AN-PW activity and co-expression of AN-PW+AaLTP3+AaPDR2 genes yielded AN and AB in necrotic N. benthamiana leaves at 13 days post-agroinfiltration. This newly discovered function of LTPs opens up new possibilities for the engineering of biosynthesis pathways of high value terpenes in heterologous expression systems.

Jasmonates: signal transduction components and their roles in environmental stress responses

Jasmonates, oxylipin-type plant hormones, are implicated in diverse aspects of plant growth development and interaction with the environment. Following diverse developmental and environmental cues, jasmonate is produced, conjugated to the amino acid isoleucine and perceived by a co-receptor complex composed of the Jasmonate ZIM-domain (JAZ) repressor proteins and an E3 ubiquitin ligase complex containing the F-box CORONATINE INSENSITIVE 1 (COI1). This event triggers the degradation of the JAZ proteins and the release of numerous transcription factors, including MYC2 and its homologues, which are otherwise bound and inhibited by the JAZ repressors. Here, we will review the role of the COI1, JAZ and MYC2 proteins in the interaction of the plant with its environment, illustrating the significance of jasmonate signalling, and of the proteins involved, for responses to both biotic stresses caused by insects and numerous microbial pathogens and abiotic stresses caused by adverse climatic conditions. It has also become evident that crosstalk with other hormone signals, as well as light and clock signals, plays an important role in the control and fine-tuning of these stress responses. Finally, we will discuss how several pathogens exploit the jasmonate perception and early signalling machinery to decoy the plants defence systems.

Genome-Scale Transcriptome Analysis of the Desert Shrub Artemisia sphaerocephala

BACKGROUND

Artemisia sphaerocephala, a semi-shrub belonging to the Artemisia genus of the Compositae family, is an important pioneer plant that inhabits moving and semi-stable sand dunes in the deserts and steppes of northwest and north-central China. It is very resilient in extreme environments. Additionally, its seeds have excellent nutritional value, and the abundant lipids and polysaccharides in the seeds make this plant a potential valuable source of bio-energy. However, partly due to the scarcity of genetic information, the genetic mechanisms controlling the traits and environmental adaptation capacity of A. sphaerocephala are unknown.

RESULTS

Here, we present the first in-depth transcriptomic analysis of A. sphaerocephala. To maximize the representation of conditional transcripts, mRNA was obtained from 17 samples, including living tissues of desert-growing A. sphaerocephala, seeds germinated in the laboratory, and calli subjected to no stress (control) and high and low temperature, high and low osmotic, and salt stresses. De novo transcriptome assembly performed using an Illumina HiSeq 2500 platform resulted in the generation of 68,373 unigenes. We analyzed the key genes involved in the unsaturated fatty acid synthesis pathway and identified 26 A. sphaerocephala fad2 genes, which is the largest fad2 gene family reported to date. Furthermore, a set of genes responsible for resistance to extreme temperatures, salt, drought and a combination of stresses was identified.

CONCLUSION

The present work provides abundant genomic information for functional dissection of the important traits of A. sphaerocephala and contributes to the current understanding of molecular adaptive mechanisms of A. sphaerocephala in the desert environment. Identification of the key genes in the unsaturated fatty acid synthesis pathway could increase understanding of the biological regulatory mechanisms of fatty acid composition traits in plants and facilitate genetic manipulation of the fatty acid composition of oil crops.

Overexpression of AaWRKY1 Leads to an Enhanced Content of Artemisinin in Artemisia annua

Artemisinin is an effective component of drugs against malaria. The regulation of artemisinin biosynthesis is at the forefront of artemisinin research. Previous studies showed that AaWRKY1 can regulate the expression of ADS, which is the first key enzyme in artemisinin biosynthetic pathway. In this study, AaWRKY1 was cloned, and it activated ADSpro and CYPpro in tobacco using dual-LUC assay. To further study the function of AaWRKY1, pCAMBIA2300-AaWRKY1 construct under 35S promoter was generated. Transgenic plants containing AaWRKY1 were obtained, and four independent lines with high expression of AaWRKY1 were analyzed. The expression of ADS and CYP, the key enzymes in artemisinin biosynthetic pathway, was dramatically increased in AaWRKY1-overexpressing A. annua plants. Furthermore, the artemisinin yield increased significantly in AaWRKY1-overexpressing A. annua plants. These results showed that AaWRKY1 increased the content of artemisinin by regulating the expression of both ADS and CYP. It provides a new insight into the mechanism of regulation on artemisinin biosynthesis via transcription factors in the future.

Expression of β-glucosidase increases trichome density and artemisinin content in transgenic Artemisia annua plants

Artemisinin is highly effective against multidrug-resistant strains of Plasmodium falciparum, the aetiological agent of the most severe form of malaria. However, a low level of accumulation of artemisinin in Artemisia annua is a major limitation for its production and delivery to malaria endemic areas of the world. While several strategies to enhance artemisinin have been extensively explored, enhancing storage capacity in trichome has not yet been considered. Therefore, trichome density was increased with the expression of β-glucosidase (bgl1) gene in A. annua through Agrobacterium-mediated transformation. Transgene (bgl1) integration and transcript were confirmed by molecular analysis. Trichome density increased up to 20% in leaves and 66% in flowers of BGL1 transgenic plants than Artemisia control plants. High-performance liquid chromatography, time of flight mass spectrometer data showed that artemisinin content increased up to 1.4% in leaf and 2.56% in flowers (per g DW), similar to the highest yields achieved so far through metabolic engineering. Artemisinin was enhanced up to five-fold in BGL1 transgenic flowers. This study opens the possibility of increasing artemisinin content by manipulating trichomes' density, which is a major reservoir of artemisinin. Combining biosynthetic pathway engineering with enhancing trichome density may further increase artemisinin yield in A. annua. Because oral feeding of Artemisia plant cells reduced parasitemia more efficiently than the purified drug, reduced drug resistance and cost of prohibitively expensive purification process, enhanced expression should play a key role in making this valuable drug affordable to treat malaria in a large global population that disproportionally impacts low-socioeconomic areas and underprivileged children.

Identification and Functional Characterization of Sesquiterpene Synthases from Xanthium strumarium

Xanthium strumarium synthesizes various pharmacologically active sesquiterpenes. The molecular characterization of sesquiterpene biosynthesis in X. strumarium has not been reported so far. In this study, the cDNAs coding for three sesquiterpene synthases (designated as XsTPS1, XsTPS2 and XsTPS3) were isolated using the X. strumarium transcriptome that we recently constructed. XsTPS1, XsTPS2 and XsTPS3 were revealed to have primary activities forming germacrene D, guaia-4,6-diene and germacrene A, respectively, by either ectopic expression in yeast cells or purified recombinant protein-based in vitro assays. Quantitative real-time PCRs and metabolite analysis for the different plant parts showed that the transcript abundance of XsTPS1-XsTPS3 is consistent with the accumulation pattern of their enzymatic products, supporting their biochemical functions in vivo. In particular, we discovered that none of the XsTPS2 product, guaia-4,6-diene, can be detected in one of the X. strumarium cultivars used in this study (it was named the Hubei-cultivar), in which a natural deletion of two A bases in the XsTPS2 cDNA disrupts its activity, which further confirmed the proposed biochemical role of XsTPS2 in X. strumarium in vivo.

Proteomic analysis of the medicinal plant Artemisia annua: Data from leaf and trichome extracts

This article contains raw and processed data related to research published by Bryant et al.[1]. Data was obtained by MS-based proteomics, analysing trichome-enriched, trichome-depleted and whole leaf samples taken from the medicinal plant Artemisia annua and searching the acquired MS/MS data against a recently published contig database [2] and other genomic and proteomic sequence databases for comparison. The processed data shows that an order-of-magnitude more proteins have been identified from trichome-enriched Artemisia annua samples in comparison to previously published data. Proteins known to have a role in the biosynthesis of artemisinin and other highly abundant proteins were found which imply additional enzymatically driven processes occurring within the trichomes that are significant for the biosynthesis of artemisinin.