[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.

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.

[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.

Integrated multi-omics analysis and microbial recombinant protein system reveal hydroxylation and glycosylation involving nevadensin biosynthesis in Lysionotus pauciflorus.

BACKGROUND: Karst-adapted plant, Lysionotus pauciflours accumulates special secondary metabolites with a wide range of pharmacological effects for surviving in drought and high salty areas, while researchers focused more on their environmental adaptations and evolutions. Nevadensin (5,7-dihydroxy-6,8,4'-trimethoxyflavone), the main active component in L. pauciflours, has unique bioactivity of such as anti-inflammatory, anti-tubercular, and anti-tumor or cancer. Complex decoration of nevadensin, such as hydroxylation and glycosylation of the flavone skeleton determines its diversity and biological activities. The lack of omics data limits the exploration of accumulation mode and biosynthetic pathway. Herein, we integrated transcriptomics, metabolomics, and microbial recombinant protein system to reveal hydroxylation and glycosylation involving nevadensin biosynthesis in L. pauciflours. RESULTS: Up to 275 flavonoids were found to exist in L. pauciflorus by UPLC-MS/MS based on widely targeted metabolome analysis. The special flavone nevadensin (5,7-dihydroxy-6,8,4'-trimethoxyflavone) is enriched in different tissues, as are its related glycosides. The flavonoid biosynthesis pathway was drawn based on differential transcripts analysis, including 9 PAL, 5 C4H, 8 4CL, 6 CHS, 3 CHI, 1 FNSII, and over 20 OMTs. Total 310 LpCYP450s were classified into 9 clans, 36 families, and 35 subfamilies, with 56% being A-type CYP450s by phylogenetic evolutionary analysis. According to the phylogenetic tree with AtUGTs, 187 LpUGTs clustered into 14 evolutionary groups (A-N), with 74% being E, A, D, G, and K groups. Two LpCYP82D members and LpUGT95 were functionally identified in Saccharomyces cerevisiae and Escherichia coli, respectively. CYP82D-8 and CYP82D-1 specially hydroxylate the 6- or 8-position of A ring in vivo and in vitro, dislike the function of F6H or F8H discovered in basil which functioned depending on A-ring substituted methoxy. These results refreshed the starting mode that apigenin can be firstly hydroxylated on A ring in nevadensin biosynthesis. Furthermore, LpUGT95 clustered into the 7-OGT family was verified to catalyze 7-O glucosylation of nevadensin accompanied with weak nevadensin 5-O glucosylation function, firstly revealed glycosylation modification of flavones with completely substituted A-ring. CONCLUSIONS: Metabolomic and full-length transcriptomic association analysis unveiled the accumulation mode and biosynthetic pathway of the secondary metabolites in the karst-adapted plant L. pauciflorus. Moreover, functional identification of two LpCYP82D members and one LpUGT in microbe reconstructed the pathway of nevadensin biosynthesis.

[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.

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.

[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.

Genome-Wide Identification and Functional Characterization of UDP-Glucosyltransferase Genes Involved in Flavonoid Biosynthesis in Glycine max.

Flavonoids, natural products abundant in the model legume Glycine max, confer benefits to plants and to animal health. Flavonoids are present in soybean mainly as glycoconjugates. However, the mechanisms of biosynthesis of flavonoid glycosides are largely unknown in G. max. In the present study, 212 putative UDP-glycosyltransferase (UGT) genes were identified in G. max by genome-wide searching. The GmUGT genes were distributed differentially among the 20 chromosomes, and they were expressed in various tissues with distinct expression profiles. We further analyzed the enzymatic activities of 11 GmUGTs that are potentially involved in flavonoid glycosylation, and found that six of them (UGT72X4, UGT72Z3, UGT73C20, UGT88A13, UGT88E19 and UGT92G4) exhibited activity toward flavonol, isoflavone, flavone and flavanol aglycones with different kinetic properties. Among them, UGT72X4, UGT72Z3 and UGT92G4 are flavonol-specific UGTs, and UGT73C20 and UGT88E19 exhibited activity toward both flavonol and isoflavone aglycones. In particular, UGT88A13 exhibited activity toward epicatechin, but not for the flavonol aglycones kaempferol and quercetin. Overexpression of these six GmUGT genes significantly increased the contents of isoflavone and flavonol glucosides in soybean hairy roots. In addition, overexpression of these six GmUGT genes also affected flavonol glycoside contents differently in seedlings and seeds of transgenic Arabidopsis thaliana. We provide valuable information on the identification of all UGT genes in soybean, and candidate GmUGT genes for potential metabolic engineering of flavonoid compounds in both Escherichia coli and plants.

Involvement of three putative glucosyltransferases from the UGT72 family in flavonol glucoside/rhamnoside biosynthesis in Lotus japonicus seeds.

Flavonols are one of the largest groups of flavonoids that confer benefits for the health of plants and animals. Flavonol glycosides are the predominant flavonoids present in the model legume Lotus japonicus. The molecular mechanisms underlying the biosynthesis of flavonol glycosides as yet remain unknown in L. japonicus. In the present study, we identified a total of 188 UDP-glycosyltransferases (UGTs) in L. japonicus by genome-wide searching. Notably, 12 UGTs from the UGT72 family were distributed widely among L. japonicus chromosomes, expressed in all tissues, and showed different docking scores in an in silico bioinformatics docking analysis. Further enzymatic assays showed that five recombinant UGTs (UGT72AD1, UGT72AF1, UGT72AH1, UGT72V3, and UGT72Z2) exhibit activity toward flavonol, flavone, and isoflavone aglycones. In particular, UGT72AD1, UGT72AH1, and UGT72Z2 are flavonol-specific UGTs with different kinetic properties. In addition, the overexpression of UGT72AD1 and UGT72Z2 led to increased accumulation of flavonol rhamnosides in L. japonicus and Arabidopsis thaliana. Moreover, the increase of kaempferol 3-O-rhamnoside-7-O-rhamnoside in transgenic A. thaliana inhibited root growth as compared with the wild-type control. These results highlight the significance of the UGT72 family in flavonol glycosylation and the role of flavonol rhamnosides in plant growth.

GmMYB58 and GmMYB205 are seed-specific activators for isoflavonoid biosynthesis in Glycine max.

KEY MESSAGE: GmMYB58 and GmMYB205 are key positive regulators that are involved in isoflavonoid biosynthesis in seeds of Glycine max, and they activate the expression of several structural genes in the isoflavonoid pathway. MYB transcription factors (TFs) are major regulators involved in flavonoid/isoflavonoid biosynthesis in many plant species. However, functions of most MYB TFs remain unknown in flavonoid/isoflavonoid pathway in Glycine max. In this study, we identified 321 MYB TFs by genome-wide searching, and further isolated and functionally characterized two MYB TFs, GmMYB58 and GmMYB205. The deduced GmMYB58 and GmMYB205 proteins contain highly conserved R2R3 repeat domain at the N-terminal region that is the signature motif of R2R3-type MYB TFs. GmMYB58 and GmMYB205 were highly expressed in early seed development stages than in the other tested organs. GmMYB58 and GmMYB205 GFP fusion proteins were found to be localized in the nucleus when they were transiently expressed in Arabidopsis thaliana mesophyll protoplast. Both GmMYB58 and GmMYB205 can activate the promoter activities of GmCHS, GmIFS2, and GmHID in the transient trans-activation assays, and the activation of GmHID by both GmMYB58 and GmMYB205 was further confirmed by yeast one-hybrid assay. In addition, over-expression of GmMYB58 and GmMYB205 resulted in significant increases in expression levels of several pathway genes in soybean hairy roots, in particular, IFS2 by more than fivefolds in GmMYB205-over-expressing lines. Moreover, isoflavonoid contents were remarkably enhanced in the GmMYB58 and GmMYB205 over-expressing hairy roots than in the control. Our results suggest that GmMYB58 and GmMYB205 are seed-specific TFs, and they can enhance isoflavonoid biosynthesis mainly through the regulation of GmIFS2 and GmHID in G. max.

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.

Role of a chalcone isomerase-like protein in flavonoid biosynthesis in Arabidopsis thaliana.

Flavonoids are important natural products for plant defence and human health. Although almost all the flavonoid pathway genes have been well-documented by biochemical and/or genetic approaches, the role of the Arabidopsis chalcone isomerase-like (CHIL) gene remains unclear. Two chil mutants with a seed colour similar to that of wild-type Arabidopsis have been identified here, but in sharp contrast to the characteristic transparent testa seed phenotype associated with other known flavonoid pathway genes. CHIL loss-of-function mutations led to a strong reduction in the proanthocyanidin and flavonol levels in seeds, but not in the anthocyanin levels in leaves. CHIL over-expression could partially recover the mutant phenotype of the chil mutant and increased both proanthocyanidin and flavonol accumulation in wild-type Arabidopsis. However, the CHIL gene could not rescue the mutant phenotype of TT5 that encodes the intrinsic chalcone isomerase in Arabidopsis. Parallel phenotypical and metabolic analyses of the chil, tt5, chs, and f3h mutants revealed that, genetically, CHIL functions at the same step as TT5. Moreover, it is demonstrated that CHIL co-expresses, co-localizes, and interacts with TT5 in Arabidopsis for flavonoid production. Based on these genetic and metabolic studies, it is concluded that CHIL functions with TT5 to promote flavonoid production, which is a unique enhancer in the flavonoid pathway.

The APETALA2-MYBL2 module represses proanthocyanidin biosynthesis by affecting formation of the MBW complex in seeds of Arabidopsis thaliana.

Proanthocyanidins (PAs) are the second most abundant plant phenolic natural products. PA biosynthesis is regulated by the well-documented MYB/bHLH/WD40 (MBW) complex, but how this complex itself is regulated remains ill defined. Here, in situ hybridization and ß-glucuronidase staining show that APETALA2 (AP2), a well-defined regulator of flower and seed development, is strongly expressed in the seed coat endothelium, where PAs accumulate. AP2 negatively regulates PA content and expression levels of key PA pathway genes. AP2 activates MYBL2 transcription and interacts with MYBL2, a key suppressor of the PA pathway. AP2 exerts its function by directly binding to the AT-rich motifs near the promoter region of MYBL2. Molecular and biochemical analyses revealed that AP2 forms AP2-MYBL2-TT8/EGL3 complexes, disrupting the MBW complex and thereby repressing expression of ANR, TT12, TT19, and AHA10. Genetic analyses revealed that AP2 functions upstream of MYBL2, TT2, and TT8 in PA regulation. Our work reveals a new role of AP2 as a key regulator of PA biosynthesis in Arabidopsis. Overall, this study sheds new light on the comprehensive regulation network of PA biosynthesis as well as the dual regulatory roles of AP2 in seed development and accumulation of major secondary metabolites in Arabidopsis.

Integrated Full-Length Transcriptomics and Metabolomics Reveal Glycosyltransferase Involved in the Biosynthesis of Flavonol Glycosides in Laportea bulbifera.

Many species of the Urticaceae family are important cultivated fiber plants that are known for their economic and industrial values. However, their secondary metabolite profiles and associated biosynthetic mechanisms have not been well-studied. Using Laportea bulbifera as a model, we conducted widely targeted metabolomics, which revealed 523 secondary metabolites, including a unique accumulation of flavonol glycosides in bulblet. Through full-length transcriptomic and RNA-seq analyses, the related genes in the flavonoid biosynthesis pathway were identified. Finally, weighted gene correlation network analysis and functional characterization revealed four LbUGTs, including LbUGT78AE1, LbUGT72CT1, LbUGT71BX1, and LbUGT71BX2, can catalyze the glycosylation of flavonol aglycones (kaempferol, myricetin, gossypetin, and quercetagetin) using UDP-Gal and UDP-Glu as the sugar donors. LbUGT78AE1 and LbUGT72CT1 showed substrate promiscuity, whereas LbUGT71BX1 and LbUGT71BX2 exhibited different substrate and sugar donor selectivity. These results provide a genetic resource for studying Laportea in the Urticaceae family, as well as key enzymes responsible for the metabolism of valuable flavonoid glycosides.

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.

Identification of Laportea bulbifera using the complete chloroplast genome as a potentially effective super-barcode.

Laportea bulbifera, a Miao medicine grown in karst areas, has exerted a unique curative effect on skin itching in the elderly, with an annual sales of > 100 million Yuan. Owing to the shortage of resources and large morphological variations in L. bulbifera, it is difficult to identify the species correctly using only traditional methods, which seriously affects the safety of drug usage for patients. This study obtained the complete high-quality L. bulbifera chloroplast (cp) genome, using second- and third-generation high-throughput sequencing. The cp genome was 149,911 bp in length, with a typical quadripartite structure. A total of 127 genes were annotated, including 83 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. There was an inverted small single copy (SSC) structure in the L. bulbifera cp genome, one large-scale rearrangement of ~ 39 kb excised in the SSC and IR regions. The complete cp genome sequence is used as a potentially effective super-barcode and the highly variable regions (ycf1, matK, and ndhD) can be used as potentially specific barcodes to accurately distinguish L. bulbifera from counterfeits and closely related species. This study is important for the identification of L. bulbifera and lays a theoretical foundation for elucidating the phylogenetic relationship of the species.

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.

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.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis.

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.