[Simultaneous determination of seven artemisinin-related compounds in Artemisia annua by UPLC-QQQ-MS/MS].

A variety of compounds in Artemisia annua were simultaneously determined to evaluate the quality of A. annua from multiple perspectives. A method based on ultra-high performance liquid chromatography-triple quadrupole tandem mass spectrometry(UPLC-QQQ-MS/MS) was established for the simultaneous determination of seven compounds: amorpha-4,11-diene, artemisinic aldehyde, dihydroartemisinic acid, artemisinic acid, artemisinin B, artemisitene, and artemisinin, in A. annua. The content of the seven compounds in different tissues(roots, stems, leaves, and lateral branches) of A. annua were compared. The roots, stems, leaves, and lateral branches of four-month-old A. annua were collected and the content of seven artemisinin-related compounds in different tissues was determined. A multi-reaction monitoring(MRM) acquisition mode of UPLC-QQQ-MS/MS was used, with a positive ion mode of atmospheric pressure chemical ion source(APCI). Chromatographic separation was achieved on an Eclipse Plus RRHD C_(18) column(2.1 mm×50 mm, 1.8 µm). The gradient elution was performed with the mobile phase consisted of formic acid(0.1%)-ammonium formate(5 mmol·L~(-1))(A) and the methanol(B) gradient program of 0-8 min, 55%-100% B, 8-11 min, 100% B, and equilibrium for 3 min, the flow rate of 0.6 mL·min~(-1), the column temperature of 40 ℃, the injection volume of 5 µL, and the detection time of 8 min. Through methodological investigation, a method based on UPLC-QQQ-MS/MS was established for the simultaneous quantitative determination of seven representative compounds involved in the biosynthesis of artemisinin. The content of artemisinin in A. annua was higher than that of artemisinin B, and the content of artemisinin and dihydroartemisinic acid were high in all the tissues of A. annua. The content of the seven compounds varied considerably in different tissues, with the highest levels in the leaves and neither artemisinene nor artemisinic aldehyde was detected in the roots. In this study, a quantitative method based on UPLC-QQQ-MS/MS for the simultaneous determination of seven representative compounds involved in the biosynthesis of artemisinin was established, which was accurate, sensitive, and highly efficient, and can be used for determining the content of artemisinin-related compounds in A. annua, breeding new varieties, and controlling the quality of Chinese medicinal materials.

Multi-omics analysis reveals promiscuous O-glycosyltransferases involved in the diversity of flavonoid glycosides in Periploca forrestii (Apocynaceae).

Flavonoid glycosides are widespread in plants, and are of great interest owing to their diverse biological activities and effectiveness in preventing chronic diseases. Periploca forrestii, a renowned medicinal plant of the Apocynaceae family, contains diverse flavonoid glycosides and is clinically used to treat rheumatoid arthritis and traumatic injuries. However, the mechanisms underlying the biosynthesis of these flavonoid glycosides have not yet been elucidated. In this study, we used widely targeted metabolomics and full-length transcriptome sequencing to identify flavonoid diversity and biosynthetic genes in P. forrestii. A total of 120 flavonoid glycosides, including 21 C-, 96 O-, and 3 C/O-glycosides, were identified and annotated. Based on 24,123 full-length coding sequences, 99 uridine diphosphate sugar-utilizing glycosyltransferases (UGTs) were identified and classified into 14 groups. Biochemical assays revealed that four UGTs exhibited O-glycosyltransferase activity toward apigenin and luteolin. Among them, PfUGT74B4 and PfUGT92A8 were highly promiscuous and exhibited multisite O-glycosylation or consecutive glycosylation activities toward various flavonoid aglycones. These four glycosyltransferases may significantly contribute to the diversity of flavonoid glycosides in P. forrestii. Our findings provide a valuable genetic resource for further studies on P. forrestii and insights into the metabolic engineering of bioactive flavonoid glycosides.

[Genome-wide identification of bZIP family genes and screening of candidate AarbZIPs involved in terpenoid biosynthesis in Artemisia argyi].

Artemisia argyi is an important medicinal and economic plant in China, with the effects of warming channels, dispersing cold, and relieving pain, inflammation, and allergy. The essential oil of this plant is rich in volatile terpenoids and widely used in moxi-bustion and healthcare products, with huge market potential. The bZIP transcription factors compose a large family in plants and are involved in the regulation of plant growth and development, stress response, and biosynthesis of secondary metabolites such as terpenoids. However, little is known about the bZIPs and their roles in A. argyi. In this study, the bZIP transcription factors in the genome of A. argyi were systematically identified, and their physicochemical properties, phylogenetic relationship, conserved motifs, and promoter-binding elements were analyzed. Candidate AarbZIP genes involved in terpenoid biosynthesis were screened out. The results showed that a total of 156 AarbZIP transcription factors were identified at the genomic level, with the lengths of 99-618 aa, the molecular weights of 11.7-67.8 kDa, and the theoretical isoelectric points of 4.56-10.16. According to the classification of bZIPs in Arabidopsis thaliana, the 156 AarbZIPs were classified into 12 subfamilies, and the members in the same subfamily had similar conserved motifs. The cis-acting elements of promoters showed that AarbZIP genes were possibly involved in light and hormonal pathways. Five AarbZIP genes that may be involved in the regulation of terpenoid biosynthesis were screened out by homologous alignment and phylogenetic analysis. The qRT-PCR results showed that the expression levels of the five AarbZIP genes varied significantly in different tissues of A. argyi. Specifically, AarbZIP29 and AarbZIP55 were highly expressed in the leaves and AarbZIP81, AarbZIP130, and AarbZIP150 in the flower buds. This study lays a foundation for the functional study of bZIP genes and their regulatory roles in the terpenoid biosynthesis in A. argyi.

Characterization and molecular evolution analysis of Periploca forrestii inferred from its complete chloroplast genome sequence.

Periploca forrestii, a medicinal plant of the family Apocynaceae, is known as an effective and widely used clinical prescription for the treatment of rheumatoid diseases. In this study, we de novo sequenced and assembled the completement chloroplast (cp) genome of P. forrestii based on combined Oxford Nanopore PromethION and Illumina data. The cp genome was 153 724 bp in length and had four subregions. Moreover, an 84 433 bp large single-copy and a 17 731 bp small single-copy were separated by 25 780 bp inverted repeats (IRs). The cp genome included 132 genes with 18 duplicates in the IRs. A total of 45 repeat structures and 183 simple sequence repeats were detected. Codon usage showed a bias toward A/T-ending codons. A comparative study of Apocynaceae revealed that an IR expansion occurred on P. forrestii. The Ka/Ks values of eight species of Apocynaceae suggested that positive selection was exerted on the psaI and ycf2 genes, which might reflect specific adaptions to the P. forrestii particular growth environment. Phylogenetic analysis indicated that Periplocoideae was a sister to Asclepiadoideae, forming a monophyletic group in the family Apocynaceae. This study provided an important P. forrestii genomic resource for future evolutionary studies and the phylogenetic reconstruction of the family Apocynaceae.

Comparative and phylogenetic analysis of complete chloroplast genomes from five Artemisia species.

Artemisia Linn. is a large genus within the family Asteraceae that includes several important medicinal plants. Because of their similar morphology and chemical composition, traditional identification methods often fail to distinguish them. Therefore, developing an effective identification method for Artemisia species is an urgent requirement. In this study, we analyzed 15 chloroplast (cp) genomes, including 12 newly sequenced genomes, from 5 Artemisia species. The cp genomes from the five Artemisia species had a typical quadripartite structure and were highly conserved across species. They had varying lengths of 151,132-151,178 bp, and their gene content and codon preferences were similar. Mutation hotspot analysis identified four highly variable regions, which can potentially be used as molecular markers to identify Artemisia species. Phylogenetic analysis showed that the five Artemisia species investigated in this study were sister branches to each other, and individuals of each species formed a monophyletic clade. This study shows that the cp genome can provide distinguishing features to help identify closely related Artemisia species and has the potential to serve as a universal super barcode for plant identification.

[Chloroplast genome structure characteristics and phylogenetic analysis of Artemisia indica].

Artemisia indica is an important medicinal plant in the Asteraceae family, but its molecular genetic information has been rarely reported. In this study, the chloroplast genome of A. indica was sequenced, assembled, and annotated by the high-throughput sequencing technology, and its sequence characteristics, repeat sequences, codon usage bias, and phylogeny were analyzed. The results showed that the length of the chloroplast genome for A. indica was 151 161 bp, which was a typical circular four-segment structure, including two inverted repeat regions(IRs), a large single-copy(LSC) region, and a small single-copy(SSC) region, with a GC content of 37.47%. A total of 132 genes were annotated, and 114 were obtained after de-duplication, including 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Fifty long repeat sequences and 191 SSRs were detected in the chloroplast genome of A. indica, and SSRs were mainly single nucleotides. Codon usage bias analysis showed that leucine was the most frequently used amino acid(10.77%) in the chloroplast genome, and there were 30 codons with relative synonymous codon usage(RSCU)>1 and all ended with A/U. The phylogenetic tree constructed based on the chloroplast genomes of the 19 species from the Asteraceae family showed that A. indica and A. argyi were closest in the genetic relationship, and Artemisia species clustered into separate evolutionary branches. The results of this study are expected to provide a theoretical basis for the genetic diversity and resource conservation of Artemisia medicinal plants.

Identification of Specific Glycosyltransferases Involved in Flavonol Glucoside Biosynthesis in Ginseng Using Integrative Metabolite Profiles, DIA Proteomics, and Phylogenetic Analysis.

Ginseng contains a variety of flavonol glycosides that possess diverse biological activities; however, scant information of flavonoid glycosylation was reported in ginseng. We found that panasenoside and kaempferol 3-O-glucoside were commonly accumulated along with cultivation years in leaves. In order to explore the procedure of flavonol glycosylation in ginseng, 50 UDP-glycosyltransferases (UGTs) were screened out using differentiated data-independent acquisition (DIA) proteomics and phylogenetic analysis. UGT92A10 and UGT94Q4 were found contributing to the formation of kaempferol 3-O-glucoside. UGT73A18, UGT74T4, and UGT75W1 could catalyze galactosylation of kaempferol 3-O-glucoside. Ser278, Trp335, Gln338, and Val339 were found forming hydrogen bonds with UDP-galactose in UGT75W1 by docking. MeJA induced transcripts of UGT73A18 and UGT74T4 by over fourfold, consistent with the decrease of kaempferol 3-O-glucoside, which indicated that these genes may be related to resisting adversity stress in ginseng. These results highlight the significance of integrative metabolite profiles, proteomics, and phylogenetic analysis for exploring flavonol glycosylation in ginseng.

Genome-wide analyses reveals a glucosyltransferase involved in rutin and emodin glucoside biosynthesis in tartary buckwheat.

With people's increasing needs for health concern, rutin and emodin in tartary buckwheat have attracted much attention for their antioxidant, anti-diabetic and reducing weight function. However, the biosynthesis of rutin and emodin in tartary buckwheat is still unclear; especially their later glycosylation contributing to make them more stable and soluble is uncovered. Based on tartary buckwheat' genome, the gene structures of 106 UGTs were analyzed; 21 candidate FtUGTs were selected to enzymatic test by comparing their transcript patterns. Among them, FtUGT73BE5 and other 4 FtUGTs were identified to glucosylate flavonol or emodin in vitro; especially rFtUGT73BE5 could catalyze the glucosylation of all tested flavonoids and emodin. Furthermore, the identical in vivo functions of FtUGT73BE5 were demonstrated in tartary buckwheat hairy roots. The transcript profile of FtUGT73BE5 was consistent with the accumulation trend of rutin in plant; this gene may relate to anti-adversity for its transcripts were up-regulated by MeJA, and repressed by ABA.

[Molecular genetics of medicinal plants, past achievement, and new era].

Unraveling the genetic basis of medicinal plant metabolism and developmental traits is a long-standing goal for pharmacologists and plant biologists. This paper discusses the definition of molecular genetics of medicinal plants, which is an integrative discipline with medicinal plants as the research object. This discipline focuses on the heredity and variation of medicinal plants, and elucidates the relationship between the key traits of medicinal plants(active compounds, yield, resistance, etc.) and genotype, studies the structure and function, heredity and variation of medicinal plant genes mainly at molecular level, so as to reveal the molecular mechanisms of transmission, expression and regulation of genetic information of medicinal plants. Specifically, we emphasize on three major aspects of this discipline.(1)Individual and population genetics of medicinal plants, this part mainly highlights the genetic mechanism of the domestication, the individual genomics at the species level, and the formation of genetic diversity of medicinal plants.(2)Elucidation of biosynthetic pathways of active compounds and their evolutionary significance. This part summarizes the biosynthesis, diversity and molecular evolution of active compounds in medicinal plants.(3) Molecular mechanisms that shaping the key agronomic traits by internal and external factors. This part focuses on the accumulation and distribution of active compounds within plants and the regulation of metabolic network by environmental factors. Finally, we prospect the future direction of molecular genetics of medicinal plants based on the rapid development of multi-omics technology, as well as the application of molecular genetics in the future strategies to achieve conservation and breeding of medicinal plants and efficient biosynthesis of active compounds.

[Synthesis of triterpenoid saponins from Aesculus chinensis based on transcriptome data].

Aesculus chinensis belongs to Hippocastanaceae family,bears medicinal and ornamental values. The oleanane type triterpenoid saponin aescin is regarded as active ingredient and accumulated in seed. In order to understand its molecular basis of the triterpenoid biosynthesis,we used high-throughput sequencing under Illumina Hi Seq 2000 platform to obtain the transcriptome data of seed and flower from A. chinensis to further mine the genes involved in its metabolic pathway. Unigene's de novo splicing was performed using Trinity software; the transcriptome results were annotated with KEGG database to predict the specific pathways of the aescin triterpenoid metabolism. Terpenoid and triterpenoid pathways were found from transcriptome data,and forty seven and twenty seven corresponding genes were uncovered respectively. It was found that there are eight kinds of enzymes related to the terpenoid metabolism pathway precursors and three kinds of enzymes related to the triterpenoid metabolism pathway. In this study,five genes corresponding to triterpene cyclase were analyzed in A. chinensis for the first time,which may participate in the synthesis of triterpenoid. It' s revealed that there were thirty three differential genes associated with the ko00900 and ko00909 pathways by analysis on the difference in transcriptome expression between seeds and flowers; seventeen unigenes were up-regulated and sixteen unigenes were down-regulated in the seeds relative to flowers. In this study, qRT-PCR experiments were used to verify the expression of three key enzyme genes of SQE( Unigene25806),HMGS( Unigene36710),and ß-AS( Unigene33291). The results of qRT-PCR were consistent with the transcriptome data. The candidate genes related to triterpenoid saponin aescin synthesis in A. chinensis found in this study can provide theoretical basis for the metabolism synthesis and regulation of aescin.

The genome of the medicinal plant Andrographis paniculata provides insight into the biosynthesis of the bioactive diterpenoid neoandrographolide.

Andrographis paniculata is a herbaceous dicot plant widely used for its anti-inflammatory and anti-viral properties across its distribution in China, India and other Southeast Asian countries. A. paniculata was used as a crucial therapeutic treatment during the influenza epidemic of 1919 in India, and is still used for the treatment of infectious disease in China. A. paniculata produces large quantities of the anti-inflammatory diterpenoid lactones andrographolide and neoandrographolide, and their analogs, which are touted to be the next generation of natural anti-inflammatory medicines for lung diseases, hepatitis, neurodegenerative disorders, autoimmune disorders and inflammatory skin diseases. Here, we report a chromosome-scale A. paniculata genome sequence of 269 Mb that was assembled by Illumina short reads, PacBio long reads and high-confidence (Hi-C) data. Gene annotation predicted 25 428 protein-coding genes. In order to decipher the genetic underpinning of diterpenoid biosynthesis, transcriptome data from seedlings elicited with methyl jasmonate were also obtained, which enabled the identification of genes encoding diterpenoid synthases, cytochrome P450 monooxygenases, 2-oxoglutarate-dependent dioxygenases and UDP-dependent glycosyltransferases potentially involved in diterpenoid lactone biosynthesis. We further carried out functional characterization of pairs of class-I and -II diterpene synthases, revealing the ability to produce diversified labdane-related diterpene scaffolds. In addition, a glycosyltransferase able to catalyze O-linked glucosylation of andrograpanin, yielding the major active product neoandrographolide, was also identified. Thus, our results demonstrate the utility of the combined genomic and transcriptomic data set generated here for the investigation of the production of the bioactive diterpenoid lactone constituents of the important medicinal herb A. paniculata.

Panax ginseng genome examination for ginsenoside biosynthesis.

Ginseng, which contains ginsenosides as bioactive compounds, has been regarded as an important traditional medicine for several millennia. However, the genetic background of ginseng remains poorly understood, partly because of the plant's large and complex genome composition. We report the entire genome sequence of Panax ginseng using next-generation sequencing. The 3.5-Gb nucleotide sequence contains more than 60% repeats and encodes 42 006 predicted genes. Twenty-two transcriptome datasets and mass spectrometry images of ginseng roots were adopted to precisely quantify the functional genes. Thirty-one genes were identified to be involved in the mevalonic acid pathway. Eight of these genes were annotated as 3-hydroxy-3-methylglutaryl-CoA reductases, which displayed diverse structures and expression characteristics. A total of 225 UDP-glycosyltransferases (UGTs) were identified, and these UGTs accounted for one of the largest gene families of ginseng. Tandem repeats contributed to the duplication and divergence of UGTs. Molecular modeling of UGTs in the 71st, 74th, and 94th families revealed a regiospecific conserved motif located at the N-terminus. Molecular docking predicted that this motif captures ginsenoside precursors. The ginseng genome represents a valuable resource for understanding and improving the breeding, cultivation, and synthesis biology of this key herb.

Identification of UDP-glycosyltransferases involved in the biosynthesis of astringent taste compounds in tea (Camellia sinensis).

Galloylated catechins and flavonol 3-O-glycosides are characteristic astringent taste compounds in tea (Camellia sinensis). The mechanism involved in the formation of these metabolites remains unknown in tea plants. In this paper, 178 UGT genes (CsUGTs) were identified inC. sinensis based on an analysis of tea transcriptome data. Phylogenetic analysis revealed that 132 of these genes were clustered into 15 previously established phylogenetic groups (A to M, O and P) and a newly identified group R. Three of the 11 recombinant UGT proteins tested were found to be involved in the in vitro biosynthesis of ß-glucogallin and glycosylated flavonols. CsUGT84A22 exhibited catalytic activity toward phenolic acids, in particular gallic acid, to produce ß-glucogallin, which is the immediate precursor of galloylated catechin biosynthesis in tea plants. CsUGT78A14 and CsUGT78A15 were found to be responsible for the biosynthesis of flavonol 3-O-glucosides and flavonol 3-O-galactosides, respectively. Site-directed mutagenesis of the Q373H substitution for CsUGT78A14 indicated that the Q (Gln) residue played a catalytically crucial role for flavonoid 3-O-glucosyltransferase activity. The expression profiles of the CsUGT84A22, CsUGT78A14, and CsUGT78A15 genes were correlated with the accumulation patterns of ß-glucogallin and the glycosylated flavonols which indicated that these three CsUGT genes were involved in the biosynthesis of astringent compounds inC. sinensis.