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

Functional analysis of CYP71AV1 reveals the evolutionary landscape of artemisinin biosynthesis

Artemisinin biosynthesis, unique to Artemisia annua, is suggested to have evolved from the ancestral costunolide biosynthetic pathway commonly found in the Asteraceae family. However, the evolutionary landscape of this process is not fully understood. The first oxidase in artemisinin biosynthesis, CYP71AV1, also known as amorpha-4,11-diene oxidase (AMO), has specialized from ancestral germacrene A oxidases (GAOs). Unlike GAO, which exhibits catalytic promiscuity toward amorpha-4,11-diene, the natural substrate of AMO, AMO has lost its ancestral activity on germacrene A. Previous studies have suggested that the loss of the GAO copy in A. annua is responsible for the abolishment of the costunolide pathway. In the genome of A. annua, there are two copies of AMO, each of which has been reported to be responsible for the different product profiles of high- and low-artemisinin production chemotypes. Through analysis of their tissue-specific expression and comparison of their sequences with those of other GAOs, it was discovered that one copy of AMO (AMOHAP) exhibits a different transcript compared to the reported artemisinin biosynthetic genes and shows more sequence similarity to other GAOs in the catalytic regions. Furthermore, in a subsequent in vitro enzymatic assay, the recombinant protein of AMOHAP unequivocally demonstrated GAO activity. This result clearly indicates that AMOHAP is a GAO rather than an AMO and that its promiscuous activity on amorpha-4,11-diene has led to its misidentification as an AMO in previous studies. In addition, the divergent expression pattern of AMOHAP compared to that of the upstream germacrene A synthase may have contributed to the abolishment of costunolide biosynthesis in A. annua. Our findings reveal a complex evolutionary landscape in which the emergence of a new metabolic pathway replaces an ancestral one.

Graphene enhances artemisinin production in the traditional medicinal plant Artemisia annua via dynamic physiological processes and miRNA regulation

We investigated the effects of graphene on the model herb Artemisia annua, which is renowned for producing artemisinin, a widely used pharmacological compound. Seedling growth and biomass were promoted when A. annua was cultivated with low concentrations of graphene, an effect which was attributed to a 1.4-fold increase in nitrogen uptake, a 15%-22% increase in chlorophyll fluorescence, and greater abundance of carbon cycling-related bacteria. Exposure to 10 or 20 mg/L graphene resulted in a ∼60% increase in H2O2, and graphene could act as a catalyst accelerator, leading to a 9-fold increase in catalase (CAT) activity in vitro and thereby maintaining reactive oxygen species (ROS) homeostasis. Importantly, graphene exposure led to an 80% increase in the density of glandular secreting trichomes (GSTs), in which artemisinin is biosynthesized and stored. This contributed to a 5% increase in artemisinin content in mature leaves. Interestingly, expression of miR828 was reduced by both graphene and H2O2 treatments, resulting in induction of its target gene AaMYB17, a positive regulator of GST initiation. Subsequent molecular and genetic assays showed that graphene-induced H2O2 inhibits micro-RNA (miRNA) biogenesis through Dicers and regulates the miR828-AaMYB17 module, thus affecting GST density. Our results suggest that graphene may contribute to yield improvement in A. annua via dynamic physiological processes together with miRNA regulation, and it may thus represent a new cultivation strategy for increasing yield capacity through nanobiotechnology.

Molecular regulation of the key specialized metabolism pathways in medicinal plants

The basis of modern pharmacology is the human ability to exploit the production of specialized metabolites from medical plants, for example, terpenoids, alkaloids, and phenolic acids. However, in most cases, the availability of these valuable compounds is limited by cellular or organelle barriers or spatio-temporal accumulation patterns within different plant tissues. Transcription factors (TFs) regulate biosynthesis of these specialized metabolites by tightly controlling the expression of biosynthetic genes. Cutting-edge technologies and/or combining multiple strategies and approaches have been applied to elucidate the role of TFs. In this review, we focus on recent progress in the transcription regulation mechanism of representative high-value products and describe the transcriptional regulatory network, and future perspectives are discussed, which will help develop high-yield plant resources.

Advanced metabolic engineering strategies for increasing artemisinin yield in Artemisia annua L

Artemisinin, also known as 'Qinghaosu', is a chemically sesquiterpene lactone containing an endoperoxide bridge. Due to the high activity to kill Plasmodium parasites, artemisinin and its derivatives have continuously served as the foundation for antimalarial therapies. Natural artemisinin is unique to the traditional Chinese medicinal plant Artemisia annua L., and its content in this plant is low. This has motivated the synthesis of this bioactive compound using yeast, tobacco, and Physcomitrium patens systems. However, the artemisinin production in these heterologous hosts is low and cannot fulfil its increasing clinical demand. Therefore, A. annua plants remain the major source of this bioactive component. Recently, the transcriptional regulatory networks related to artemisinin biosynthesis and glandular trichome formation have been extensively studied in A. annua. Various strategies including (i) enhancing the metabolic flux in artemisinin biosynthetic pathway; (ii) blocking competition branch pathways; (iii) using transcription factors (TFs); (iv) increasing peltate glandular secretory trichome (GST) density; (v) applying exogenous factors; and (vi) phytohormones have been used to improve artemisinin yields. Here we summarize recent scientific advances and achievements in artemisinin metabolic engineering, and discuss prospects in the development of high-artemisinin yielding A. annua varieties. This review provides new insights into revealing the transcriptional regulatory networks of other high-value plant-derived natural compounds (e.g., taxol, vinblastine, and camptothecin), as well as glandular trichome formation. It is also helpful for the researchers who intend to promote natural compounds production in other plants species.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

AaABCG20 transporter involved in cutin and wax secretion affects the initiation and development of glandular trichomes in Artemisia annua

Glandular trichomes are specialized structures found on the surface of plants to produce specific compounds, including terpenes, alkaloids, and other organic substances. Artemisia annua, commonly known as sweet wormwood, synthesizes and stores the antimalarial drug artemisinin in glandular trichomes. Previous research indicated that increasing the glandular trichome density could enhance artemisinin production, and the cuticle synthesis affected the initiation and development of glandular trichomes in A. annua. In this study, AaABCG12 and AaABCG20 were isolated from A. annua that exhibited similar expression patterns to artemisinin biosynthetic genes. Of the two, AaABCG20 acted as a specific transporter in glandular trichomes. Downregulating the expression of AaABCG20 resulted in a notable reduction in the density of glandular trichome, while overexpressing AaABCG20 resulted in an increase in glandular trichome density. GC-MS analysis demonstrated that AaABCG20 was responsible for the transport of cutin and wax in A. annua. These findings indicated that AaABCG20 influenced the initiation and development of glandular trichomes through transporting cutin and wax in A. annua. This glandular trichome specific half-size ABCG-type transporter is crucial in facilitating the transportation of cutin and wax components, ultimately contributing to the successful initiation and development of glandular trichomes.

IMP: bridging the gap for medicinal plant genomics

Medicinal plants have garnered significant attention in ethnomedicine and traditional medicine due to their potential antitumor, anti-inflammatory and antioxidant properties. Recent advancements in genome sequencing and synthetic biology have revitalized interest in natural products. Despite the availability of sequenced genomes and transcriptomes of these plants, the absence of publicly accessible gene annotations and tabular formatted gene expression data has hindered their effective utilization. To address this pressing issue, we have developed IMP (Integrated Medicinal Plantomics), a freely accessible platform at https://www.bic.ac.cn/IMP. IMP curated a total of 8 565 672 genes for 84 high-quality genome assemblies, and 2156 transcriptome sequencing samples encompassing various organs, tissues, developmental stages and stimulations. With the integrated 10 analysis modules, users could simply examine gene annotations, sequences, functions, distributions and expressions in IMP in a one-stop mode. We firmly believe that IMP will play a vital role in enhancing the understanding of molecular metabolic pathways in medicinal plants or plants with medicinal benefits, thereby driving advancements in synthetic biology, and facilitating the exploration of natural sources for valuable chemical constituents like drug discovery and drug production.

The genomes of Dahlia pinnata, Cosmos bipinnatus, and Bidens alba in tribe Coreopsideae provide insights into polyploid evolution and inulin biosynthesis

BACKGROUND

The Coreopsideae tribe, a subset of the Asteraceae family, encompasses economically vital genera like Dahlia, Cosmos, and Bidens, which are widely employed in medicine, horticulture, ecology, and food applications. Nevertheless, the lack of reference genomes hinders evolutionary and biological investigations in this tribe.

RESULTS

Here, we present 3 haplotype-resolved chromosome-level reference genomes of the tribe Coreopsideae, including 2 popular flowering plants (Dahlia pinnata and Cosmos bipinnatus) and 1 invasive weed plant (Bidens alba), with assembled genome sizes 3.93 G, 1.02 G, and 1.87 G, respectively. We found that Gypsy transposable elements contribute mostly to the larger genome size of D. pinnata, and multiple chromosome rearrangements have occurred in tribe Coreopsideae. Besides the shared whole-genome duplication (WGD-2) in the Heliantheae alliance, our analyses showed that D. pinnata and B. alba each underwent an independent recent WGD-3 event: in D. pinnata, it is more likely to be a self-WGD, while in B. alba, it is from the hybridization of 2 ancestor species. Further, we identified key genes in the inulin metabolic pathway and found that the pseudogenization of 1-FEH1 and 1-FEH2 genes in D. pinnata and the deletion of 3 key residues of 1-FFT proteins in C. bipinnatus and B. alba may probably explain why D. pinnata produces much more inulin than the other 2 plants.

CONCLUSIONS

Collectively, the genomic resources for the Coreopsideae tribe will promote phylogenomics in Asteraceae plants, facilitate ornamental molecular breeding improvements and inulin production, and help prevent invasive weeds.

Strategies on biosynthesis and production of bioactive compounds in medicinal plants

Medicinal plants are a valuable source of essential medicines and herbal products for healthcare and disease therapy. Compared with chemical synthesis and extraction, the biosynthesis of natural products is a very promising alternative for the successful conservation of medicinal plants, and its rapid development will greatly facilitate the conservation and sustainable utilization of medicinal plants. Here, we summarize the advances in strategies and methods concerning the biosynthesis and production of natural products of medicinal plants. The strategies and methods mainly include genetic engineering, plant cell culture engineering, metabolic engineering, and synthetic biology based on multiple "OMICS" technologies, with paradigms for the biosynthesis of terpenoids and alkaloids. We also highlight the biosynthetic approaches and discuss progress in the production of some valuable natural products, exemplifying compounds such as vindoline (alkaloid), artemisinin and paclitaxel (terpenoids), to illustrate the power of biotechnology in medicinal plants.

AaABI5 transcription factor mediates light and abscisic acid signaling to promote anti-malarial drug artemisinin biosynthesis in Artemisia annua

Artemisia annua, a member of the Asteraceae family, remains the primary source of artemisinin. However, the artemisinin content in the existing varieties of this plant is very low. In this study, we found that the environmental factors light and phytohormone abscisic acid (ABA) could synergistically promote the expression of artemisinin biosynthetic genes. Notably, the increased expression levels of those genes regulated by ABA depended on light. Gene expression analysis found that AaABI5, a transcription factor belonging to the basic leucine zipper (bZIP) family, was inducible by the light and ABA treatment. Analysis of AaABI5-overexpressing and -suppressing lines suggested that AaABI5 could enhance artemisinin biosynthesis and activate the expression of four core biosynthetic genes. In addition, the key regulator of light-induced artemisinin biosynthesis, AaHY5, could bind to the promoter of AaABI5 and activate its expression. In conclusion, our results demonstrated that AaABI5 acts as an important molecular junction for the synergistic promotion of artemisinin biosynthesis by light and ABA signals, which provides a candidate gene for developing new germplasms of high-quality A. annua.

Applications of some advanced sequencing, analytical, and computational approaches in medicinal plant research: a review

The potential active chemicals found in medicinal plants, which have long been employed as natural medicines, are abundant. Exploring the genes responsible for producing these compounds has given new insights into medicinal plant research. Previously, the authentication of medicinal plants was done via DNA marker sequencing. With the advancement of sequencing technology, several new techniques like next-generation sequencing, single molecule sequencing, and fourth-generation sequencing have emerged. These techniques enshrined the role of molecular approaches for medicinal plants because all the genes involved in the biosynthesis of medicinal compound(s) could be identified through RNA-seq analysis. In several research insights, transcriptome data have also been used for the identification of biosynthesis pathways. miRNAs in several medicinal plants and their role in the biosynthesis pathway as well as regulation of the disease-causing genes were also identified. In several research articles, an in silico study was also found to be effective in identifying the inhibitory effect of medicinal plant-based compounds against virus' gene(s). The use of advanced analytical methods like spectroscopy and chromatography in metabolite proofing of secondary metabolites has also been reported in several recent research findings. Furthermore, advancement in molecular and analytic methods will give new insight into studying the traditionally important medicinal plants that are still unexplored.

LcMYB106 suppresses monoterpene biosynthesis by negatively regulating LcTPS32 expression in Litsea cubeba

Litsea cubeba, the core species of the Lauraceae family, is valuable for the production of essential oils due to its high concentration of monoterpenes (90%). The key monoterpene synthase and metabolic regulatory network of monoterpene biosynthesis have provided new insights for improving essential oil content. However, there are few studies on the regulation mechanism of monoterpenes in L. cubeba. In this study, we investigated LcTPS32, a member of the TPS-b subfamily, and identified its function as an enzyme for the synthesis of monoterpenes, including geraniol, α-pinene, β-pinene, β-myrcene, linalool and eucalyptol. The quantitative real-time PCR analysis showed that LcTPS32 was highly expressed in the fruits of L. cubeba and contributed to the characteristic flavor of its essential oil. Overexpression of LcTPS32 resulted in a significant increase in the production of monoterpenes in L. cubeba by activating both the MVA and MEP pathways. Additionally, the study revealed that LcMYB106 played a negative regulatory role in monoterpenes biosynthesis by directly binding to the promoter of LcTPS32. Our study indicates that LcMYB106 could serve as a crucial target for metabolic engineering endeavors, aiming at enhancing the monoterpene biosynthesis in L. cubeba.

Progress and prospect: Biosynthesis of plant natural products based on plant chassis

Plant-derived natural products are a specific class of active substances with numerous applications in the medical, energy, and industrial fields. Many of these substances are in high demand and have become the fundamental materials for various purposes. Recently, the use of synthetic biology to produce plant-derived natural products has become a significant trend. Plant chassis, in particular, offer unique advantages over microbial chassis in terms of cell structure, product affinity, safety, and storage. The development of the plant hairy root tissue culture system has accelerated the commercialization and industrialization of synthetic biology in the production of plant-derived natural products. This paper will present recent progress in the synthesis of various plant natural products using plant chassis, organized by the types of different structures. Additionally, we will summarize the four primary types of plant chassis used for synthesizing natural products from plant sources and review the enabling technologies that have contributed to the development of synthetic biology in recent years. Finally, we will present the role of isolated and combined use of different optimization strategies in breaking the upper limit of natural product production in plant chassis. This review aims to provide practical references for synthetic biologists and highlight the great commercial potential of plant chassis biosynthesis, such as hairy roots.

Functions of Representative Terpenoids and Their Biosynthesis Mechanisms in Medicinal Plants

Terpenoids are the broadest and richest group of chemicals obtained from plants. These plant-derived terpenoids have been extensively utilized in various industries, including food and pharmaceuticals. Several specific terpenoids have been identified and isolated from medicinal plants, emphasizing the diversity of biosynthesis and specific functionality of terpenoids. With advances in the technology of sequencing, the genomes of certain important medicinal plants have been assembled. This has improved our knowledge of the biosynthesis and regulatory molecular functions of terpenoids with medicinal functions. In this review, we introduce several notable medicinal plants that produce distinct terpenoids (e.g., Cannabis sativa, Artemisia annua, Salvia miltiorrhiza, Ginkgo biloba, and Taxus media). We summarize the specialized roles of these terpenoids in plant-environment interactions as well as their significance in the pharmaceutical and food industries. Additionally, we highlight recent findings in the fields of molecular regulation mechanisms involved in these distinct terpenoids biosynthesis, and propose future opportunities in terpenoid research, including biology seeding, and genetic engineering in medicinal plants.

Drought risk assessment on the eastern part of Indian peninsula-a study on Purulia district, West Bengal

This study focuses on measuring the spatial nature of drought risk which is conceived as the product of drought severity, drought vulnerability, and drought exposure in the Purulia district, located in the eastern part of the Indian peninsula. Drought severity is measured using the Standard Precipitation Index and drought vulnerability is calculated as the average condition of meteorological, hydrological, agricultural, and socio-economic drought. The drought types and drought exposure conditions are the outcome of multi-criteria analysis where the Fuzzy Analytical Hierarchy Process is used for assigning weights to the respective parameters and the Analytical Hierarchy Process is used for determining the class ranks. 31.46% of the total district area has registered moderate to high and high vulnerability to drought situations, while 16.57% of the entire district area has been found moderate to high and highly exposed to drought situations. Similarly, 39.39% of the district's total area is under a significant drought risk. Blocks like Barabazar (75.49%), Jhalda-I (71.85%), and Purulia-II (52.66%) have the majority of their area under extreme drought risk conditions. The modeled outcome of drought vulnerability was found significant while being tested with phenomena highly correlated to drought events, land surface temperature, and aridity index. The computed spatial profile of the districts' drought risk condition is of substantial help for the policymakers in preparing effective drought mitigation measures to restrict drought impacts reasonably.

Molecular understanding of anthocyanin biosynthesis activated by PAP1 and regulated by 2, 4-dichlorophenoxyacetic acid in engineered red Artemisia annua cells

MAIN CONCLUSION

Eight promoters were cloned, from which AC and G-box cis-elements were identified. PAP1 enhanced the promoter activity. 2,4-D reduced the anthocyanin biosynthesis via downregulating the expression of the PAP1 transgene. Artemisia annua is an effective antimalarial medicinal crop. We have established anthocyanin-producing red cell cultures from this plant with the overexpression of Production of Anthocyanin Pigment 1 (PAP1) encoding a R2R3MYB transcription factor. To understand the molecular mechanism by which PAP1 activated the entire anthocyanin pathway, we mined the genomic sequences of A. annua and obtained eight promoters of the anthocyanin pathway genes. Sequence analysis identified four types of AC cis-elements from six promoters, the MYB response elements (MRE) bound by PAP1. In addition, six promoters were determined to have at least one G-box cis-element. Eight promoters were cloned for activity analysis. Dual luciferase assays showed that PAP1 significantly enhanced the promoting activity of seven promoters, indicating that PAP1 turned on the biosynthesis of anthocyanins via the activation of these pathway gene expression. To understand how 2,4-dichlorophenoxyacetic acid (2,4-D), an auxin, regulates the PAP1-activated anthocyanin biosynthesis, five different concentrations (0, 0.05, 0.5, 2.5, and 5 µM) were tested to characterize anthocyanin production and profiles. The resulting data showed that the concentrations tested decreased the fresh weight of callus growth, anthocyanin levels, and the production of anthocyanins per Petri dish. HPLC-qTOF-MS/MS-based profiling showed that these concentrations did not alter anthocyanin profiles. Real-time RT-PCR was completed to characterize the expression PAP1 and four representative pathway genes. The results showed that the five concentrations reduced the expression levels of the constitutive PAP1 transgene and three pathway genes significantly and eliminated the expression of the chalcone synthase gene either significantly or slightly. These data indicate that the constitutive PAP1 expression depends on gradients added in the medium. Based on these findings, the regulation of 2,4-D is discussed for anthocyanin engineering in red cells of A. annua.

Promoter variations in DBR2-like affect artemisinin production in different chemotypes of Artemisia annua

Artemisia annua is the only known plant source of the potent antimalarial artemisinin, which occurs as the low- and high-artemisinin producing (LAP and HAP) chemotypes. Nevertheless, the different mechanisms of artemisinin producing between these two chemotypes were still not fully understood. Here, we performed a comprehensive analysis of genome resequencing, metabolome, and transcriptome data to systematically compare the difference in the LAP chemotype JL and HAP chemotype HAN. Metabolites analysis revealed that 72.18% of sesquiterpenes was highly accumulated in HAN compared to JL. Integrated omics analysis found a DBR2-Like (DBR2L) gene may be involved in artemisinin biosynthesis. DBR2L was highly homologous with DBR2, belonged to ORR3 family, and had the DBR2 activity of catalyzing artemisinic aldehyde to dihydroartemisinic aldehyde. Genome resequencing and promoter cloning revealed that complicated variations existed in DBR2L promoters among different varieties of A. annua and were clustered into three variation types. The promoter activity of diverse variant types showed obvious differences. Furthermore, the core region (-625 to 0) of the DBR2L promoter was identified and candidate transcription factors involved in DBR2L regulation were screened. Thus, the result indicates that DBR2L is another key enzyme involved in artemisinin biosynthesis. The promoter variation in DBR2L affects its expression level, and thereby may result in the different yield of artemisinin in varieties of A. annua. It provides a novel insight into the mechanism of artemisinin-producing difference in LAP and HAP chemotypes of A. annua, and will assist in a high yield of artemisinin in A. annua.

Chromosome-level reference genome of stinkwort, Dittrichia graveolens (L.) Greuter: A resource for studies on invasion, range expansion, and evolutionary adaptation under global change

Dittrichia graveolens (L.) Greuter, or stinkwort, is a weedy annual plant within the family Asteraceae. The species is recognized for the rapid expansion of both its native and introduced ranges: in Europe, it has expanded its native distribution northward from the Mediterranean basin by nearly 7 °C latitude since the mid-20th century, while in California and Australia the plant is an invasive weed of concern. Here, we present the first de novo D. graveolens genome assembly (1N = 9 chromosomes), including complete chloroplast (151,013 bp) and partial mitochondrial genomes (22,084 bp), created using Pacific Biosciences HiFi reads and Dovetail Omni-C data. The final primary assembly is 835 Mbp in length, of which 98.1% are represented by 9 scaffolds ranging from 66 to 119 Mbp. The contig N50 is 74.9 Mbp and the scaffold N50 is 96.9 Mbp, which, together with a 98.8% completeness based on the BUSCO embryophyta10 database containing 1,614 orthologs, underscores the high quality of this assembly. This pseudo-molecule-scale genome assembly is a valuable resource for our fundamental understanding of the genomic consequences of range expansion under global change, as well as comparative genomic studies in the Asteraceae.

The Light- and Jasmonic Acid-Induced AaMYB108-like Positive Regulates the Initiation of Glandular Secretory Trichome in Artemisia annua L

The plant Artemisia annua L. is famous for producing "artemisinin", which is an essential component in the treatment of malaria. The glandular secretory trichomes (GSTs) on the leaves of A. annua secrete and store artemisinin. Previous research has demonstrated that raising GST density can effectively raise artemisinin content. However, the molecular mechanism of GST initiation is not fully understood yet. In this study, we identified an MYB transcription factor, the AaMYB108-like, which is co-induced by light and jasmonic acid, and positively regulates glandular secretory trichome initiation in A. annua. Overexpression of the AaMYB108-like gene in A. annua increased GST density and enhanced the artemisinin content, whereas anti-sense of the AaMYB108-like gene resulted in the reduction in GST density and artemisinin content. Further experiments demonstrated that the AaMYB108-like gene could form a complex with AaHD8 to promote the expression of downstream AaHD1, resulting in the initiation of GST. Taken together, the AaMYB108-like gene is a positive regulator induced by light and jasmonic acid for GST initiation in A. annua.

Medicinal plant genomics

Recent developments in plant genomics have enabled a comprehensive analysis of the medicinal potential of plants based on their gene repertoire. Genes of biosynthesis pathways can be discovered through comparative genomics and through integration of transcriptomic data. Data-driven discovery of specialized metabolites could accelerate research.

Multi-omics reveal key enzymes involved in the formation of phenylpropanoid glucosides in Artemisia annua

Although mainly known for producing artemisinin, Artemisia annua is enriched in phenylpropanoid glucosides (PGs) with significant bioactivities. However, the biosynthesis of A. annua PGs is insufficiently investigated. Different A. annua ecotypes from distinct growing environments accumulate varying amounts of metabolites, including artemisinin and PGs such as scopolin. UDP-glucose:phenylpropanoid glucosyltransferases (UGTs) transfers glucose from UDP-glucose in PG biosynthesis. Here, we found that the low-artemisinin ecotype GS produces a higher amount of scopolin, compared to the high-artemisinin ecotype HN. By combining transcriptome and proteome analyses, we selected 28 candidate AaUGTs from 177 annotated AaUGTs. Using AlphaFold structural prediction and molecular docking, we determined the binding affinities of 16 AaUGTs. Seven of the AaUGTs enzymatically glycosylated phenylpropanoids. AaUGT25 converted scopoletin to scopolin and esculetin to esculin. The lack of accumulation of esculin in the leaf and the high catalytic efficiency of AaUGT25 on esculetin suggest that esculetin is methylated to scopoletin, the precursor of scopolin. We also discovered that AaOMT1, a previously uncharacterized O-methyltransferase, converts esculetin to scopoletin, suggesting an alternative route for producing scopoletin, which contributes to the high-level accumulation of scopolin in A. annua leaves. AaUGT1 and AaUGT25 responded to induction of stress-related phytohormones, implying the involvement of PGs in stress responses.

Genome-wide identification of bZIP transcription factor family in Artemisia annua, its transcriptional profiling and regulatory role in phenylpropanoid metabolism under different light conditions

The basic leucine zipper (bZIP) protein transcription factors are known to modulate development, plant growth, metabolic response, and resistance to several biotic and abiotic stressors and have been widely studied in the model plant Arabidopsis thaliana. However, no comprehensive information about the bZIP transcription factor family in Artemisia annua has been explored to date. In this genome-wide study, we identified 61 bZIP TFs after removing false positives and incomplete sequences from Artemisia annua. Seven highly expressed homolog AabZIP TF genes under UV-B and differential light conditions in different tissues were identified from the publicly available microarray dataset as having their cis-regulatory elements involved in, flavonoids biosynthesis, seed-specific gene regulation, stress responses, and metabolic regulation. In-silico analysis and electrophoretic mobility shift assay (EMSA) confirmed the interaction of AabZIP19 TF over the AaPAL1 promoter in order to regulate the phenolics and flavonoid biosynthesis via the phenylpropanoid pathway. Further, RT-PCR analysis has been carried out to validate the transcript levels of selected AabZIP genes under white light, red light, blue light (45 min), and UV-B exposure (12 and 24 h). These genes have their highest expression levels under UV-B and blue light exposure, in contrast with white light. Therefore, the detection of ROS through staining confirms the accumulation of superoxide radicals and H2O2, and in addition to reducing ROS accumulation under UV-B and blue light irradiation, total phenols and flavonoids are significantly enhanced. This study laid the groundwork for deciphering the possible role of AabZIP TFs under different light stress-responsive conditions and in the regulation of secondary metabolism.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-023-01338-0.

Chemical profile and analysis of biosynthetic pathways and genes of volatile terpenes in Pityopsis ruthii, a rare and endangered flowering plant

It is critical to gather biological information about rare and endangered plants to incorporate into conservation efforts. The secondary metabolism of Pityopsis ruthii, an endangered flowering plant that only occurs along limited sections of two rivers (Ocoee and Hiwassee) in Tennessee, USA was studied. Our long-term goal is to understand the mechanisms behind P. ruthii's adaptation to restricted areas in Tennessee. Here, we profiled the secondary metabolites, specifically in flowers, with a focus on terpenes, aiming to uncover the genomic and molecular basis of terpene biosynthesis in P. ruthii flowers using transcriptomic and biochemical approaches. By comparative profiling of the nonpolar portion of metabolites from various tissues, P. ruthii flowers were rich in terpenes, which included 4 monoterpenes and 10 sesquiterpenes. These terpenes were emitted from flowers as volatiles with monoterpenes and sesquiterpenes accounting for almost 68% and 32% of total emission of terpenes, respectively. These findings suggested that floral terpenes play important roles for the biology and adaptation of P. ruthii to its limited range. To investigate the biosynthesis of floral terpenes, transcriptome data for flowers were produced and analyzed. Genes involved in the terpene biosynthetic pathway were identified and their relative expressions determined. Using this approach, 67 putative terpene synthase (TPS) contigs were detected. TPSs in general are critical for terpene biosynthesis. Seven full-length TPS genes encoding putative monoterpene and sesquiterpene synthases were cloned and functionally characterized. Three catalyzed the biosynthesis of sesquiterpenes and four catalyzed the biosynthesis of monoterpenes. In conclusion, P. ruthii plants employ multiple TPS genes for the biosynthesis of a mixture of floral monoterpenes and sesquiterpenes, which probably play roles in chemical defense and attracting insect pollinators alike.

Future antimalarials from Artemisia? A rationale for natural product mining against drug-refractory Plasmodium stages

Covering: up to 2023Infusions of the plants Artemisia annua and A. afra are gaining broad popularity to prevent or treat malaria. There is an urgent need to address this controversial public health question by providing solid scientific evidence in relation to these uses. Infusions of either species were shown to inhibit the asexual blood stages, the liver stages including the hypnozoites, but also the sexual stages, the gametocytes, of Plasmodium parasites. Elimination of hypnozoites and sterilization of mature gametocytes remain pivotal elements of the radical cure of P. vivax, and the blockage of P. vivax and P. falciparum transmission, respectively. Drugs active against these stages are restricted to the 8-aminoquinolines primaquine and tafenoquine, a paucity worsened by their double dependence on the host genetic to elicit clinical activity without severe toxicity. Besides artemisinin, these Artemisia spp. contain many natural products effective against Plasmodium asexual blood stages, but their activity against hypnozoites and gametocytes was never investigated. In the context of important therapeutic issues, we provide a review addressing (i) the role of artemisinin in the bioactivity of these Artemisia infusions against specific parasite stages, i.e., alone or in association with other phytochemicals; (ii) the mechanisms of action and biological targets in Plasmodium of ca. 60 infusion-specific Artemisia phytochemicals, with an emphasis on drug-refractory parasite stages (i.e., hypnozoites and gametocytes). Our objective is to guide the strategic prospecting of antiplasmodial natural products from these Artemisia spp., paving the way toward novel antimalarial "hit" compounds either naturally occurring or Artemisia-inspired.

Two-step model of paleohexaploidy, ancestral genome reshuffling and plasticity of heat shock response in Asteraceae

An ancient hexaploidization event in the most but not all Asteraceae plants, may have been responsible for shaping the genomes of many horticultural, ornamental, and medicinal plants that promoting the prosperity of the largest angiosperm family on the earth. However, the duplication process of this hexaploidy, as well as the genomic and phenotypic diversity of extant Asteraceae plants caused by paleogenome reorganization, are still poorly understood. We analyzed 11 genomes from 10 genera in Asteraceae, and redated the Asteraceae common hexaploidization (ACH) event ~70.7-78.6 million years ago (Mya) and the Asteroideae specific tetraploidization (AST) event ~41.6-46.2 Mya. Moreover, we identified the genomic homologies generated from the ACH, AST and speciation events, and constructed a multiple genome alignment framework for Asteraceae. Subsequently, we revealed biased fractionations between the paleopolyploidization produced subgenomes, suggesting the ACH and AST both are allopolyplodization events. Interestingly, the paleochromosome reshuffling traces provided clear evidence for the two-step duplications of ACH event in Asteraceae. Furthermore, we reconstructed ancestral Asteraceae karyotype (AAK) that has 9 paleochromosomes, and revealed a highly flexible reshuffling of Asteraceae paleogenome. Of specific significance, we explored the genetic diversity of Heat Shock Transcription Factors (Hsfs) associated with recursive whole-genome polyploidizations, gene duplications, and paleogenome reshuffling, and revealed that the expansion of Hsfs gene families enable heat shock plasticity during the genome evolution of Asteraceae. Our study provides insights on polyploidy and paleogenome remodeling for the successful establishment of Asteraceae, and is helpful for further communication and exploration of the diversification of plant families and phenotypes.

Phylogenomics and morphological evolution of the mega-diverse genus Artemisia (Asteraceae: Anthemideae): implications for its circumscription and infrageneric taxonomy

BACKGROUND AND AIMS

Artemisia is a mega-diverse genus consisting of ~400 species. Despite its medicinal importance and ecological significance, a well-resolved phylogeny for global Artemisia, a natural generic delimitation and infrageneric taxonomy remain missing, owing to the obstructions from limited taxon sampling and insufficient information on DNA markers. Its morphological characters, such as capitulum, life form and leaf, show marked variations and are widely used in its infrageneric taxonomy. However, their evolution within Artemisia is poorly understood. Here, we aimed to reconstruct a well-resolved phylogeny for global Artemisia via a phylogenomic approach, to infer the evolutionary patterns of its key morphological characters and to update its circumscription and infrageneric taxonomy.

METHODS

We sampled 228 species (258 samples) of Artemisia and its allies from both fresh and herbarium collections, covering all the subgenera and its main geographical areas, and conducted a phylogenomic analysis based on nuclear single nucleotide polymorphisms (SNPs) obtained from genome skimming data. Based on the phylogenetic framework, we inferred the possible evolutionary patterns of six key morphological characters widely used in its previous taxonomy.

KEY RESULTS

The genus Kaschgaria was revealed to be nested in Artemisia with strong support. A well-resolved phylogeny of Artemisia consisting of eight highly supported clades was recovered, two of which were identified for the first time. Most of the previously recognized subgenera were not supported as monophyletic. Evolutionary inferences based on the six morphological characters showed that different states of these characters originated independently more than once.

CONCLUSIONS

The circumscription of Artemisia is enlarged to include the genus Kaschgaria. The morphological characters traditionally used for the infrageneric taxonomy of Artemisia do not match the new phylogenetic tree. They experienced a more complex evolutionary history than previously thought. We propose a revised infrageneric taxonomy of the newly circumscribed Artemisia, with eight recognized subgenera to accommodate the new results.

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.

Novel players in organogenesis and flavonoid biosynthesis in cucumber glandular trichomes

Glandular trichomes (GTs) are outgrowths of plant epidermal cells that secrete and store specialized secondary metabolites that protect plants against biotic and abiotic stresses and have economic importance for human use. While extensive work has been done to understand the molecular mechanisms of trichome organogenesis in Arabidopsis (Arabidopsis thaliana), which forms unicellular, non-glandular trichomes (NGTs), little is known about the mechanisms of GT development or regulation of secondary metabolites in plants with multicellular GTs. Here, we identified and functionally characterized genes associated with GT organogenesis and secondary metabolism in GTs of cucumber (Cucumis sativus). We developed a method for effective separation and isolation of cucumber GTs and NGTs. Transcriptomic and metabolomic analyses showed that flavonoid accumulation in cucumber GTs is positively associated with increased expression of related biosynthesis genes. We identified 67 GT development-related genes, the functions of 7 of which were validated by virus-induced gene silencing. We further validated the role of cucumber ECERIFERUM1 (CsCER1) in GT organogenesis by overexpression and RNA interference transgenic approaches. We further show that the transcription factor TINY BRANCHED HAIR (CsTBH) serves as a central regulator of flavonoid biosynthesis in cucumber glandular trichomes. Work from this study provides insight into the development of secondary metabolite biosynthesis in multi-cellular glandular trichomes.

Analyses of a chromosome-scale genome assembly reveal the origin and evolution of cultivated chrysanthemum

Chrysanthemum (Chrysanthemum morifolium Ramat.) is a globally important ornamental plant with great economic, cultural, and symbolic value. However, research on chrysanthemum is challenging due to its complex genetic background. Here, we report a near-complete assembly and annotation for C. morifolium comprising 27 pseudochromosomes (8.15 Gb; scaffold N50 of 303.69 Mb). Comparative and evolutionary analyses reveal a whole-genome triplication (WGT) event shared by Chrysanthemum species approximately 6 million years ago (Mya) and the possible lineage-specific polyploidization of C. morifolium approximately 3 Mya. Multilevel evidence suggests that C. morifolium is likely a segmental allopolyploid. Furthermore, a combination of genomics and transcriptomics approaches demonstrate the C. morifolium genome can be used to identify genes underlying key ornamental traits. Phylogenetic analysis of CmCCD4a traces the flower colour breeding history of cultivated chrysanthemum. Genomic resources generated from this study could help to accelerate chrysanthemum genetic improvement.

Transcriptional regulatory network of high-value active ingredients in medicinal plants

High-value active ingredients in medicinal plants have attracted research attention because of their benefits for human health, such as the antimalarial artemisinin, anticardiovascular disease tanshinones, and anticancer Taxol and vinblastine. Here, we review how hormones and environmental factors promote the accumulation of active ingredients, thereby providing a strategy to produce high-value drugs at a low cost. Focusing on major hormone signaling events and environmental factors, we review the transcriptional regulatory network mediating biosynthesis of representative active ingredients. In this network, many transcription factors (TFs) simultaneously control multiple synthase genes; thus, understanding the molecular mechanisms affecting transcriptional regulation of active ingredients will be crucial to developing new breeding possibilities.

AaWIN1, an AP2/ERF protein, positively regulates glandular secretory trichome initiation in Artemisia annua

Exploring the genetic network of glandular trichomes and manipulating genes relevant to secondary metabolite biosynthesis are of great importance and value. Artemisinin, a key antimalarial drug ingredient, is synthesized and stored in glandular secretory trichomes (GSTs) in Artemisia annua. WIN/SHN proteins, a clade of AP2/ERF family, are known as regulators for cuticle biosynthesis. However, their function in glandular trichome development is less unknown. In this study, we identified a WIN/SHN gene from A. annua and named it as AaWIN1. AaWIN1 was predominantly expressed in buds, flowers and trichomes, and encoded a nuclear-localized protein. Overexpressing AaWIN1 in A. annua significantly increased the density of GST as well as the artemisinin content. Furthermore, AaGSW2 was reported to play an important role in promoting GST initiation, and the expression of AaGSW2 was induced in AaWIN1-overexpression lines. AaMIXTA1, a MYB protein positively regulating trichome initiation and cuticle biosynthesis, was confirmed to interact with AaWIN1. In addition, the ectopic expression of AaWIN1 resulted in slender and curled leaves, fewer trichomes, and rising expressions of cuticle biosynthesis genes in Arabidopsis thaliana. Taken together, based on phenotype observations, content measurements and gene expression detections, AaWIN1 was considered as a positive regulator for GST initiation in A. annua.

Discovering the Repeatome of Five Species Belonging to the Asteraceae Family: A Computational Study

Genome divergence by repeat proliferation and/or loss is a process that plays a crucial role in species evolution. Nevertheless, knowledge of the variability related to repeat proliferation among species of the same family is still limited. Considering the importance of the Asteraceae family, here we present a first contribution towards the metarepeatome of five Asteraceae species. A comprehensive picture of the repetitive components of all genomes was obtained by genome skimming with Illumina sequence reads and by analyzing a pool of full-length long terminal repeat retrotransposons (LTR-REs). Genome skimming allowed us to estimate the abundance and variability of repetitive components. The structure of the metagenome of the selected species was composed of 67% repetitive sequences, of which LTR-REs represented the bulk of annotated clusters. The species essentially shared ribosomal DNA sequences, whereas the other classes of repetitive DNA were highly variable among species. The pool of full-length LTR-REs was retrieved from all the species and their age of insertion was established, showing several lineage-specific proliferation peaks over the last 15-million years. Overall, a large variability of repeat abundance at superfamily, lineage, and sublineage levels was observed, indicating that repeats within individual genomes followed different evolutionary and temporal dynamics, and that different events of amplification or loss of these sequences may have occurred after species differentiation.

Application of third-generation sequencing to herbal genomics

There is a long history of traditional medicine use. However, little genetic information is available for the plants used in traditional medicine, which limits the exploitation of these natural resources. Third-generation sequencing (TGS) techniques have made it possible to gather invaluable genetic information and develop herbal genomics. In this review, we introduce two main TGS techniques, PacBio SMRT technology and Oxford Nanopore technology, and compare the two techniques against Illumina, the predominant next-generation sequencing technique. In addition, we summarize the nuclear and organelle genome assemblies of commonly used medicinal plants, choose several examples from genomics, transcriptomics, and molecular identification studies to dissect the specific processes and summarize the advantages and disadvantages of the two TGS techniques when applied to medicinal organisms. Finally, we describe how we expect that TGS techniques will be widely utilized to assemble telomere-to-telomere (T2T) genomes and in epigenomics research involving medicinal plants.

AaMYB108 is the core factor integrating light and jasmonic acid signaling to regulate artemisinin biosynthesis in Artemisia annua

Artemisinin, a sesquiterpene compound synthesized and stored in the glandular trichome of Artemisia annua leaves, has been used to treat malaria. Previous studies have shown that both light and jasmonic acid (JA) can promote the biosynthesis of artemisinin, and the promotion of artemisinin by JA is dependent on light. However, the specific molecular mechanism remains unclear. Here, we report a MYB transcription factor, AaMYB108, identified from transcriptome analysis of light and JA treatment, as a positive regulator of artemisinin biosynthesis in A. annua. AaMYB108 promotes artemisinin biosynthesis by interacting with a previously characterized positive regulator of artemisinin, AaGSW1. Then, we found that AaMYB108 interacted with AaCOP1 and AaJAZ8, respectively. The function of AaMYB108 was influenced by AaCOP1 and AaJAZ8. Through the treatment of AaMYB108 transgenic plants with light and JA, it was found that the promotion of artemisinin by light and JA depends on the presence of AaMYB108. Taken together, our results reveal the molecular mechanism of JA regulating artemisinin biosynthesis depending on light in A. annua. This study provides new insights into the integration of light and phytohormone signaling to regulate terpene biosynthesis in plants.

Profiling of phytohormone-specific microRNAs and characterization of the miR160-ARF1 module involved in glandular trichome development and artemisinin biosynthesis in Artemisia annua

MicroRNAs (miRNAs) play crucial roles in plant development and secondary metabolism through different modes of sequence-specific interaction with their targets. Artemisinin biosynthesis is extensively regulated by phytohormones. However, the function of phytohormone-responsive miRNAs in artemisinin biosynthesis remains enigmatic. Thus, we combined the analysis of transcriptomics, small RNAs, and the degradome to generate a comprehensive resource for identifying key miRNA-target circuits involved in the phytohormone-induced process of artemisinin biosynthesis in Artemisia annua. In total, 151 conserved and 52 novel miRNAs and their 4132 targets were determined. Based on the differential expression analysis, miR160 was selected as a potential miRNA involved in artemisinin synthesis. Overexpressing MIR160 significantly impaired glandular trichome formation and suppressed artemisinin biosynthesis in A. annua, while repressing its expression resulted in the opposite effect, indicating that miR160 negatively regulates glandular trichome development and artemisinin biosynthesis. RNA ligase-mediated 5' RACE and transient transformation assays showed that miR160 mediates the RNA cleavage of Auxin Response Factor 1 (ARF1) in A. annua. Furthermore, ARF1 was shown to increase artemisinin synthesis by activating AaDBR2 expression. Taken together, our results reveal the intrinsic link between the miR160-ARF1 module and artemisinin biosynthesis, and may expedite the innovation of metabolic engineering approaches for high and stable production of artemisinin in the future.

Elevation of artemisinin content by co-transformation of artemisinin biosynthetic pathway genes and trichome-specific transcription factors in Artemisia annua

Artemisinin, derived from Artemisia annua, is currently used as the first-line treatment for malaria. However, wild-type plants have a low artemisinin biosynthesis rate. Although yeast engineering and plant synthetic biology have shown promising results, plant genetic engineering is considered the most feasible strategy, but it is also constrained by the stability of progeny development. Here we constructed three independent unique overexpressing vectors harboring three mainstream artemisinin biosynthesis enzymes HMGR, FPS, and DBR2, as well as two trichomes-specific transcription factors AaHD1 and AaORA. The simultaneous co-transformation of these vectors by Agrobacterium resulted in the successful increase of the artemisinin content in T0 transgenic lines by up to 3.2-fold (2.72%) leaf dry weight compared to the control plants. We also investigated the stability of transformation in progeny T1 lines. The results indicated that the transgenic genes were successfully integrated, maintained, and overexpressed in some of the T1 progeny plants' genomes, potentially increasing the artemisinin content by up to 2.2-fold (2.51%) leaf dry weight. These results indicated that the co-overexpression of multiple enzymatic genes and transcription factors via the constructed vectors provided promising results, which could be used to achieve the ultimate goal of a steady supply of artemisinin at affordable prices around the world.

Identification of drought-responsive phenolic compounds and their biosynthetic regulation under drought stress in Ligularia fischeri

Ligularia fischeri, a leafy edible plant found in damp shady regions, has been used as an herbal medicine and is also consumed as a horticultural crop. In this study, we investigated the physiological and transcriptomic changes, especially those involved in phenylpropanoid biosynthesis, induced by severe drought stress in L. fischeri plants. A distinguishing characteristic of L. fischeri is a color change from green to purple due to anthocyanin biosynthesis. We chromatographically isolated and identified two anthocyanins and two flavones upregulated by drought stress using liquid chromatography-mass spectrometry and nuclear magnetic resonance analyses in this plant for the first time. In contrast, all types of caffeoylquinic acids (CQAs) and flavonol contents were decreased under drought stress. Further, we performed RNA sequencing to examine the molecular changes in these phenolic compounds at the transcriptome level. In an overview of drought-inducible responses, we identified 2,105 hits for 516 distinct transcripts as drought-responsive genes. Moreover, differentially expressed genes (DEGs) associated with phenylpropanoid biosynthesis accounted for the greatest number of both up- and downregulated DEGs by Kyoto Encyclopedia of Genes and Genomes enrichment analysis. We identified 24 meaningful DEGs based on the regulation of phenylpropanoid biosynthetic genes. Potential drought-responsive genes included upregulated flavone synthase (LfFNS, TRINITY DN31661 c0 g1 i1) and anthocyanin 5-O-glucosyltransferase (LfA5GT1, TRINITY DN782 c0 g1 i1), which could contribute to the high levels of flavones and anthocyanins under drought stress in L. fischeri. In addition, the downregulated shikimate O-hydroxycinnamolytransferase (LfHCT, TRINITY DN31661 c0 g1 i1) and hydroxycinnamoyl-CoA quinate/shikimate transferase (LfHQT4, TRINITY DN15180 c0 g1 i1) genes led to a reduction in CQAs. Only one or two BLASTP hits for LfHCT were obtained for six different Asteraceae species. It is possible that the HCT gene plays a crucial role in CQAs biosynthesis in these species. These findings expand our knowledge of the response mechanisms to drought stress, particularly regarding the regulation of key phenylpropanoid biosynthetic genes in L. fischeri.

The burdock database: a multi-omic database for Arctium lappa, a food and medicinal plant

BACKGROUND

Burdock is a biennial herb of Asteraceae found in Northern Europe, Eurasia, Siberia, and China. Its mature dry fruits, called Niu Bang Zi, are recorded in various traditional Chinese medicine books. With the development of sequencing technology, the mitochondrial, chloroplast, and nuclear genomes, transcriptome, and sequence-related amplified polymorphism (SRAP) fingerprints of burdock have all been reported. To make better use of this data for further research and analysis, a burdock database was constructed.

RESULTS

This burdock multi-omics database contains two burdock genome datasets, two transcriptome datasets, eight burdock chloroplast genomes, one burdock mitochondrial genome, one A. tomentosum chloroplast genome, one A. tomentosum mitochondrial genome, 26 phenotypes of burdock varieties, burdock rhizosphere-associated microorganisms, and chemical constituents of burdock fruit, pericarp, and kernel at different growth stages (using UPLC-Q-TOF-MS). The wild and cultivation distribution of burdock in China was summarized, and the main active components and pharmacological effects of burdock currently reported were concluded. The database contains ten central functional modules: Home, Genome, Transcriptome, Jbrowse, Search, Tools, SRAP fingerprints, Associated microorganisms, Chemical, and Publications. Among these, the "Tools" module can be used to perform sequence homology alignment (Blast), multiple sequence alignment analysis (Muscle), homologous protein prediction (Genewise), primer design (Primer), large-scale genome analysis (Lastz), and GO and KEGG enrichment analyses (GO Enrichment and KEGG Enrichment).

CONCLUSIONS

The database URL is http://210.22.121.250:41352/ . This burdock database integrates molecular and chemical data to provide a comprehensive information and analysis platform for interested researchers and can be of immense help to the cultivation, breeding, and molecular pharmacognosy research of burdock.

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.

Genome-wide identification and characterisation of bHLH transcription factors in Artemisia annua

BACKGROUND

A. annua (also named Artemisia annua, sweet wormwood) is the main source of the anti-malarial drug artemisinin, which is synthesised and stored in its trichomes. Members of the basic Helix-Loop-Helix (bHLH) family of transcription factors (TFs) have been implicated in artemisinin biosynthesis in A. annua and in trichome development in other plant species.

RESULTS

Here, we have systematically identified and characterised 226 putative bHLH TFs in A. annua. All of the proteins contain a HLH domain, 213 of which also contain the basic motif that mediates DNA binding of HLH dimers. Of these, 22 also contained a Myc domain that permits dimerisation with other families of TFs; only two proteins lacking the basic motif contained a Myc domain. Highly conserved GO annotations reflected the transcriptional regulatory role of the identified TFs, and suggested conserved roles in biological processes such as iron homeostasis, and guard cell and endosperm development. Expression analysis revealed that three genes (AabHLH80, AabHLH96, and AaMyc-bHLH3) exhibited spatiotemporal expression patterns similar to genes encoding key enzymes in artemisinin synthesis.

CONCLUSIONS

This comprehensive analysis of bHLH TFs provides a new resource to direct further analysis into key molecular mechanisms underlying and regulating artemisinin biosynthesis and trichome development, as well as other biological processes, in the key medicinal plant A. annua.

Intraspecific variation in genome size in Artemisia argyi determined using flow cytometry and a genome survey

Different collections and accessions of Artemisia argyi (Chinese mugwort) harbour considerable diversity in morphology and bioactive compounds, but no mechanisms have been reported that explain these variations. We studied genome size in A. argyi accessions from different regions of China by flow cytometry. Genome size was significantly distinct among origins of these 42 Chinese mugwort accessions, ranging from 8.428 to 11.717 pg. There were no significant intraspecific differences among the 42 accessions from the five regions of China. The clustering analysis showed that these 42 A. argyi accessions could be divided into three groups, which had no significant relationship with geographical location. In a genome survey, the total genome size of A. argyi (A15) was estimated to be 7.852 Gb (or 8.029 pg) by K-mer analysis. This indicated that the results from the two independent methods are consistent, and that the genome survey can be used as an adjunct to flow cytometry to compensate for its deficiencies. In addition, genome survey can provide the information about heterozygosity, repeat sequences, GC content and ploidy of A. argyi genome. The nuclear DNA contents determined here provide a new reference for intraspecific variation in genome size in A. argyi, and may also be a potential resource for the study of genetic diversity and for breeding new cultivar.

Towards greenhouse cultivation of Artemisia annua: The application of LEDs in regulating plant growth and secondary metabolism

Artemisinin is a sesquiterpene lactone produced in glandular trichomes of Artemisia annua, and is extensively used in the treatment of malaria. Growth and secondary metabolism of A. annua are strongly regulated by environmental conditions, causing unstable supply and quality of raw materials from field grown plants. This study aimed to bring A. annua into greenhouse cultivation and to increase artemisinin production by manipulating greenhouse light environment using LEDs. A. annua plants were grown in a greenhouse compartment for five weeks in vegetative stage with either supplemental photosynthetically active radiation (PAR) (blue, green, red or white) or supplemental radiation outside PAR wavelength (far-red, UV-B or both). The colour of supplemental PAR hardly affected plant morphology and biomass, except that supplemental green decreased plant biomass by 15% (both fresh and dry mass) compared to supplemental white. Supplemental far-red increased final plant height by 23% whereas it decreased leaf area, plant fresh and dry weight by 30%, 17% and 7%, respectively, compared to the treatment without supplemental radiation. Supplemental UV-B decreased plant leaf area and dry weight (both by 7%). Interestingly, supplemental green and UV-B increased leaf glandular trichome density by 11% and 9%, respectively. However, concentrations of artemisinin, arteannuin B, dihydroartemisinic acid and artemisinic acid only exhibited marginal differences between the light treatments. There were no interactive effects of far-red and UV-B on plant biomass, morphology, trichome density and secondary metabolite concentrations. Our results illustrate the potential of applying light treatments in greenhouse production of A. annua to increase trichome density in vegetative stage. However, the trade-off between light effects on plant growth and trichome initiation needs to be considered. Moreover, the underlying mechanisms of light spectrum regulation on artemisinin biosynthesis need further clarification to enhance artemisinin yield in greenhouse production of A. annua.

Molecular insights into the genetic diversity and population structure of Artemisia annua L. as revealed by insertional polymorphisms

The knowledge about the level of genetic diversity and population structure in natural populations of Artemisia annua L. is a primary step in breeding programs for development of new cultivars with higher artemisinin level and better quality of secondary metabolites composition. We used PCR-based "retrotransposon-microsatellite amplified polymorphisms" (REMAPs) to study insertional polymorphism in A. annua genome to assess genetic variability and population structure in a collection of 118 accessions collected from north and northwest of Iran. Twenty-five primer combinations of 10 retrotransposon and seven ISSR primers amplified a total of 693 clear and unambiguous fragments in the studied accessions. The average number of bands, polymorphic bands, polymorphism, effective number of alleles, Shannon's information index and expected heterozygosity were 27.72, 24.76, 88.14%, 1.47, 0.42 and 0.28, respectively. The analysis of molecular variance revealed high genetic variation present within sampled geographical regions. Distance-based cluster analysis assigned the studied accessions into four clusters according to their geographical origin, which were also confirmed by principal coordinate analysis. In model-based Bayesian clustering, the maximum value of Δk was obtained when the collection of 118 assayed A. annua accessions assigned into two subgroups (K = 2). The results showed the high genetic variation in the collection of Iranian sweet wormwood which revealed by REMAP markers indicating the reliability and efficiency of this marker system for analysis of genetic diversity and population structure of A. annua.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s40415-022-00860-x.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

Co-option of a carotenoid cleavage dioxygenase gene (CCD4a) into the floral symmetry gene regulatory network contributes to the polymorphic floral shape-color combinations in Chrysanthemum sensu lato

Morphological novelties, including formation of trait combinations, may result from de novo gene origination and/or co-option of existing genes into other developmental contexts. A variety of shape-color combinations of capitular florets occur in Chrysanthemum and its allies. We hypothesized that co-option of a carotenoid cleavage dioxygenase gene into the floral symmetry gene network would generate a white zygomorphic ray floret. We tested this hypothesis in an evolutionary context using species in Chrysanthemum sensu lato, a monophyletic group with diverse floral shape-color combinations, based on morphological investigation, interspecific crossing, molecular interaction and transgenic experiments. Our results showed that white color was significantly associated with floret zygomorphy. Specific expression of the carotenoid cleavage dioxygenase gene CCD4a in marginal florets resulted in white color. Crossing experiments between Chrysanthemum lavandulifolium and Ajania pacifica indicated that expression of CCD4a is trans-regulated. The floral symmetry regulator CYC2g can activate expression of CCD4a with a dependence on TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING (TCP) binding element 8 on the CCD4a promoter. Based on all experimental findings, we propose that gene co-option of carotenoid degradation into floral symmetry regulation, and the subsequent dysfunction or loss of either CYC2g or CCD4a, may have led to evolution of capitular shape-color patterning in Chrysanthemum sensu lato.

Emerging trends in genomic and epigenomic regulation of plant specialised metabolism

Regulation of specialised metabolism genes is multilayered and complex, influenced by an array of genomic, epigenetic and epigenomic mechanisms. Here, we review the most recent knowledge in this field, drawing from discoveries in several plant species. Our aim is to improve understanding of how plant genome structure and function influence specialised metabolism. We also highlight key areas for future exploration. Gene regulatory mechanisms influencing specialised metabolism include gene duplication and neo-functionalization, conservation of operon-like clusters of specialised metabolism genes, local chromatin modifications, and the organisation of higher order chromatin structures within the nucleus. Genomic and epigenomic research to-date in the discipline have focused on a relatively small number of plant species, primarily at whole organ or tissue level. This is largely due to the technical demands of the experimental methods needed. However, a high degree of cell-type specificity of function exists in specialised metabolism, driven by similarly specific gene regulation. In this review we focus on the genomic characteristics of genes that are found in different types of clusters within the genome. We propose that acquisition of cell-resolution epigenomic datasets in emerging models, such as the glandular trichomes of Cannabis sativa, will yield important advances. Data such as chromatin accessibility and histone modification profiles can pinpoint which regulatory sequences are active in individual cell types and at specific times in development. These could provide fundamental biological insight as well as novel targets for genetic engineering and crop improvement.

The genomes of chicory, endive, great burdock and yacon provide insights into Asteraceae palaeo-polyploidization history and plant inulin production

Inulin is an important reserve polysaccharide in Asteraceae plants, and is also widely used as a sweetener, a source of dietary fibre and prebiotic. Nevertheless, a lack of genomic resources for inulin-producing plants has hindered extensive studies on inulin metabolism and regulation. Here, we present chromosome-level reference genomes for four inulin-producing plants: chicory (Cichorium intybus), endive (Cichorium endivia), great burdock (Arctium lappa) and yacon (Smallanthus sonchifolius), with assembled genome sizes of 1.28, 0.89, 1.73 and 2.72 Gb, respectively. We found that the chicory, endive and great burdock genomes were shaped by whole genome triplication (WGT-1), and the yacon genome was shaped by WGT-1 and two subsequent whole genome duplications (WGD-2 and WGD-3). A yacon unique whole genome duplication (WGD-3) occurred 5.6-5.8 million years ago. Our results also showed the genome size difference between chicory and endive is largely due to LTR retrotransposons, and rejected a previous hypothesis that chicory is an ancestor of endive. Furthermore, we identified fructan-active-enzyme and transcription-factor genes, and found there is one copy in chicory, endive and great burdock but two copies in yacon for most of these genes, except for the 1-FEH II gene which is significantly expanded in chicory. Interestingly, inulin synthesis genes 1-SST and 1-FFT are located close to each other, as are the degradation genes 1-FEH I and 1-FEH II. Finally, we predicted protein structures for 1-FFT genes to explore the mechanism determining inulin chain length.

From Plant to Yeast-Advances in Biosynthesis of Artemisinin

Malaria is a life-threatening disease. Artemisinin-based combination therapy (ACT) is the preferred choice for malaria treatment recommended by the World Health Organization. At present, the main source of artemisinin is extracted from Artemisia annua; however, the artemisinin content in A. annua is only 0.1-1%, which cannot meet global demand. Meanwhile, the chemical synthesis of artemisinin has disadvantages such as complicated steps, high cost and low yield. Therefore, the application of the synthetic biology approach to produce artemisinin in vivo has magnificent prospects. In this review, the biosynthesis pathway of artemisinin was summarized. Then we discussed the advances in the heterologous biosynthesis of artemisinin using microorganisms (Escherichia coli and Saccharomyces cerevisiae) as chassis cells. With yeast as the cell factory, the production of artemisinin was transferred from plant to yeast. Through the optimization of the fermentation process, the yield of artemisinic acid reached 25 g/L, thereby producing the semi-synthesis of artemisinin. Moreover, we reviewed the genetic engineering in A. annua to improve the artemisinin content, which included overexpressing artemisinin biosynthesis pathway genes, blocking key genes in competitive pathways, and regulating the expression of transcription factors related to artemisinin biosynthesis. Finally, the research progress of artemisinin production in other plants (Nicotiana, Physcomitrella, etc.) was discussed. The current advances in artemisinin biosynthesis may help lay the foundation for the remarkable up-regulation of artemisinin production in A. annua through gene editing or molecular design breeding in the future.

Genome-wide transcriptome analysis reveals the regulatory network governing terpene trilactones biosynthesis in Ginkgo biloba

Ginkgo biloba L. is currently the only remaining gymnosperm of the Ginkgoaceae Ginkgo genus, and its history can be traced back to the Carboniferous 200 million years ago. Terpene trilactones (TTLs) are one of the main active ingredients in G. biloba, including ginkgolides and bilobalide. They have a good curative effect on cardiovascular and cerebrovascular diseases because of their special antagonistic effect on platelet-activating factors. Therefore, it is necessary to deeply mine genes related to TTLs and to analyze their transcriptional regulation mechanism, which will hold vitally important scientific and practical significance for quality improvement and regulation of G. biloba. In this study, we performed RNA-Seq on the root, stem, immature leaf, mature leaf, microstrobilus, ovulate strobilus, immature fruit and mature fruit of G. biloba. The TTL regulatory network of G. biloba in different organs was revealed by different transcriptomic analysis strategies. Weighted gene co-expression network analysis (WGCNA) revealed that the five modules were closely correlated with organs. The 12 transcription factors, 5 structural genes and 24 Cytochrome P450 (CYP450) were identified as candidate regulators for TTL accumulation by WGCNA and cytoscape visualization. Finally, 6 APETALA2/ethylene response factors, 2 CYP450s and bHLH were inferred to regulate the metabolism of TTLs by correlation analysis. This study is the comprehensive in authenticating transcription factors, structural genes and CYP450 involved in TTL biosynthesis, thereby providing molecular evidence for revealing the comprehensive regulatory network involved in TTL metabolism in G. biloba.