Genome of Tripterygium wilfordii and identification of cytochrome P450 involved in triptolide biosynthesis

Demethylzeylasteral exerts potent efficacy against non-small-cell lung cancer via the P53 signaling pathway

Lung cancer has one of the highest mortality rates worldwide, with non-small-cell lung cancer (NSCLC) constituting approximately 85% of all cases. Demethylzeylasteral (DEM), extracted from Tripterygium wilfordii Hook F, exhibits notable anti-tumor properties. In this study, we revealed that DEM could effectively induce NSCLC cell apoptosis. Specifically, DEM can dose-dependently suppress the viability and migration of human NSCLC cells. RNA-seq analysis revealed that DEM regulates the P53-signaling pathway, which was further validated by assessing crucial proteins involved in this pathway. Biacore analysis indicated that DEM has high affinity with the P53 protein. The CDX model demonstrated DEM's anti-tumor actions. This work provided evidence that DEM-P53 interaction stabilizes P53 protein and triggers downstream anti-tumor activities. These findings indicate that DEM treatment holds promise as a potential therapeutic approach for NSCLC, which warrants further clinical assessment in patients with NSCLC.

Biosynthetic pathway of prescription cucurbitacin IIa and high-level production of key triterpenoid intermediates in engineered yeast and tobacco

Cucurbitacin IIa is a triterpenoid isolated exclusively from Hemsleya plants and a non-steroidal anti-inflammatory drug that functions as the main ingredient of prescription Hemslecin capsules and tablets in China. Synthetic biology provides new strategies for production of such valuable cucurbitacins at a large scale; however, the biosynthetic pathway of cucurbitacin IIa has been unknown, and the heterologous production of cucurbitacins in galactose medium has been expensive and low yielding. In this study, we characterized the functions of genes encoding two squalene epoxidases (HcSE1-2), six oxidosqualene cyclases (HcOSC1-6), two CYP450s (HcCYP87D20 and HcCYP81Q59), and an acyltransferase (HcAT1) in cucurbitacin IIa biosynthesis by heterologous expression in Saccharomyces cerevisiae and Nicotiana benthamiana. We achieved high-level production of the key cucurbitacin precursor 11-carbonyl-20β-hydroxy-Cuol from glucose in yeast via modular engineering of the mevalonate pathway and optimization of P450 expression levels. The resulting yields of 46.41 mg/l 11-carbonyl-20β-hydroxy-Cuol and 126.47 mg/l total cucurbitacin triterpenoids in shake flasks are the highest yields yet reported from engineered microbes. Subsequently, production of 11-carbonyl-20β-hydroxy-Cuol by transient gene expression in tobacco resulted in yields of 1.28 mg/g dry weight in leaves. This work reveals the key genes involved in biosynthesis of prescription cucurbitacin IIa and demonstrates that engineered yeast cultivated with glucose can produce high yields of key triterpenoid intermediates. We describe a low-cost and highly efficient platform for rapid screening of candidate genes and high-yield production of pharmacological triterpenoids.

The chromosome-level genome and functional database accelerate research about biosynthesis of secondary metabolites in Rosa roxburghii

Rosa roxburghii Tratt, a valuable plant in China with long history, is famous for its fruit. It possesses various secondary metabolites, such as L-ascorbic acid (vitamin C), alkaloids and poly saccharides, which make it a high nutritional and medicinal value. Here we characterized the chromosome-level genome sequence of R. roxburghii, comprising seven pseudo-chromosomes with a total size of 531 Mb and a heterozygosity of 0.25%. We also annotated 45,226 coding gene loci after masking repeat elements. Orthologs for 90.1% of the Complete Single-Copy BUSCOs were found in the R. roxburghii annotation. By aligning with protein sequences from public platform, we annotated 85.89% genes from R. roxburghii. Comparative genomic analysis revealed that R. roxburghii diverged from Rosa chinensis approximately 5.58 to 13.17 million years ago, and no whole-genome duplication event occurred after the divergence from eudicots. To fully utilize this genomic resource, we constructed a genomic database RroFGD with various analysis tools. Otherwise, 69 enzyme genes involved in L-ascorbate biosynthesis were identified and a key enzyme in the biosynthesis of vitamin C, GDH (L-Gal-1-dehydrogenase), is used as an example to introduce the functions of the database. This genome and database will facilitate the future investigations into gene function and molecular breeding in R. roxburghii.

Haplotype-resolved genome of Mimosa bimucronata revealed insights into leaf movement and nitrogen fixation

BACKGROUND

Mimosa bimucronata originates from tropical America and exhibits distinctive leaf movement characterized by a relative slow speed. Additionally, this species possesses the ability to fix nitrogen. Despite these intriguing traits, comprehensive studies have been hindered by the lack of genomic resources for M. bimucronata.

RESULTS

To unravel the intricacies of leaf movement and nitrogen fixation, we successfully assembled a high-quality, haplotype-resolved, reference genome at the chromosome level, spanning 648 Mb and anchored in 13 pseudochromosomes. A total of 32,146 protein-coding genes were annotated. In particular, haplotype A was annotated with 31,035 protein-coding genes, and haplotype B with 31,440 protein-coding genes. Structural variations (SVs) and allele specific expression (ASE) analyses uncovered the potential role of structural variants in leaf movement and nitrogen fixation in M. bimucronata. Two whole-genome duplication (WGD) events were detected, that occurred ~ 2.9 and ~ 73.5 million years ago. Transcriptome and co-expression network analyses revealed the involvement of aquaporins (AQPs) and Ca2+-related ion channel genes in leaf movement. Moreover, we also identified nodulation-related genes and analyzed the structure and evolution of the key gene NIN in the process of symbiotic nitrogen fixation (SNF).

CONCLUSION

The detailed comparative genomic and transcriptomic analyses provided insights into the mechanisms governing leaf movement and nitrogen fixation in M. bimucronata. This research yielded genomic resources and provided an important reference for functional genomic studies of M. bimucronata and other legume species.

Genetic diversities in wild and cultivated populations of the two closely-related medical plants species, Tripterygium Wilfordii and T. Hypoglaucum (Celastraceae)

BACKGROUND

The sustainable supply of medicinal plants is important, and cultivating and domesticating them has been suggested as an optimal strategy. However, this can lead to a loss of genetic diversity. Tripterygium wilfordii Hook. f. is a medicinal plant commonly used in traditional Chinese medicine, but its wild populations are dwindling due to excessive harvesting. To protect the species and meet the increasing demand, it is urgent to cultivate it on a large scale. However, distinguishing between T. wilfordii and T. hypoglaucum, two similar species with different medicinal properties, is challenging. Therefore, it is crucial to understand the genetic diversity and population structure of these species for their sustainable utilization.

RESULTS

In this study, we investigated the genetic diversity and population structure of the two traditional medicinal semiwoody vines plant species, Tripterygium wilfordii and T. hypoglaucum, including wild and cultivated populations using chloroplast DNA (cpDNA) sequences and microsatellite loci. Our results indicated that the two species maintain a high level of genetic divergence, indicating possible genetic bases for the different contents of bioactive compounds of the two species. T. wilfordii showed lower genetic diversity and less subdivided population structures of both markers than T. hypoglaucum. The potential factors in shaping these interesting differences might be differentiated pollen-to-seed migration rates, interbreeding, and history of population divergence. Analyses of cpDNA and microsatellite loci supported that the two species are genetically distinct entities. In addition, a significant reduction of genetic diversity was observed for cultivated populations of the two species, which mainly resulted from the small initial population size and propagated vegetative practice during their cultivation.

CONCLUSION

Our findings indicate significant genetic divergence between T. wilfordii and T. hypoglaucum. The genetic diversity and population structure analyses provide important insights into the sustainable cultivation and utilization of these medicinal plants. Accurate identification and conservation efforts are necessary for both species to ensure the safety and effectiveness of crude drug use. Our study also highlighted the importance of combined analyses of different DNA markers in addressing population genetics of medicinal plants because of the contrasts of inheritance and rates of gene flow. Large-scale cultivation programs should consider preserving genetic diversity to enhance the long-term sustainability of T. wilfordii and T. hypoglaucum. Our study proposed that some populations showed higher genetic diversity and distinctness, which can be considered with priority for conservation and as the sources for future breeding and genetic improvement.

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.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Integrating Transcriptome and Chemical Analyses to Provide Insights into Biosynthesis of Terpenoids and Flavonoids in the Medicinal Industrial Crop Andrographis paniculate and Its Antiviral Medicinal Parts

Andrographis paniculata is a medicinal plant traditionally used to produce diterpene lactones and flavonoids, which possess various biological activities. Widely distributed in China, India, and other Southeast Asia countries, A. paniculata has become an important economic crop, significantly treating SARS-CoV-2, and is being cultivated on a large scale in southern China. The biosynthesis of active ingredients in A. paniculata are regulated and controlled by genes, but their specific roles are still not fully understood. To further explore the growth regulation factors and utilization of its medicinal parts of this industrial crop, chemical and transcriptome analyses were conducted on the roots, stems, and leaves of A. paniculata to identify the biosynthesis pathways and related candidate genes of the active ingredients. The chemical analysis revealed that the main components of A. paniculata were diterpene lactones and flavonoids, which displayed potential ability to treat SARS-CoV-2 through molecular docking. Moreover, the transcriptome sequencing annotated a total of 40,850 unigenes, including 7962 differentially expressed genes. Among these, 120 genes were involved in diterpene lactone biosynthesis and 60 genes were involved in flavonoid biosynthesis. The expression of diterpene lactone-related genes was the highest in leaves and the lowest in roots, consistent with our content determination results. It is speculated that these highly expressed genes in leaves may be involved in the biosynthesis pathway of diterpenes. Furthermore, two class Ⅰ terpene synthases in A. paniculata transcriptome were also annotated, providing reference for the downstream pathway of the diterpene lactone biosynthesis. With their excellent market value, our experiments will promote the study of the biosynthetic genes for active ingredients in A. paniculata and provide insights for subsequent in vitro biosynthesis.

From genomics to metabolomics: Deciphering sanguinarine biosynthesis in Dicranostigma leptopodum

Dicranostigma leptopodum (Maxim) Fedde (DLF) is a renowned medicinal plant in China, known to be rich in alkaloids. However, the unavailability of a reference genome has impeded investigation into its plant metabolism and genetic breeding potential. Here we present a high-quality chromosomal-level genome assembly for DLF, derived using a combination of Nanopore long-read sequencing, Illumina short-read sequencing and Hi-C technologies. Our assembly genome spans a size of 621.81 Mb with an impressive contig N50 of 93.04 Mb. We show that the species-specific whole-genome duplication (WGD) of DLF and Papaver somniferum corresponded to two rounds of WGDs of Papaver setigerum. Furthermore, we integrated comprehensive homology searching, gene family analyses and construction of a gene-to-metabolite network. These efforts led to the discovery of co-expressed transcription factors, including NAC and bZIP, alongside sanguinarine (SAN) pathway genes CYP719 (CFS and SPS). Notably, we identified P6H as a promising gene for enhancing SAN production. By providing the first reference genome for Dicranostigma, our study confirms the genomic underpinning of SAN biosynthesis and establishes a foundation for advancing functional genomic research on Papaveraceae species. Our findings underscore the pivotal role of high-quality genome assemblies in elucidating genetic variations underlying the evolutionary origin of secondary metabolites.

Gene identification and semisynthesis of the anti-inflammatory oleanane-type triterpenoid wilforlide A

Wilforlide A is one of the main active constituents produced in trace amounts in Tripterygium wilfordii Hook F, which has excellent anti-inflammatory and immune suppressive effects. Despite the seeming structural simplicity of the compound, the biosynthetic pathway of wilforlide A remains unknown. Gene-specific expression analysis and genome mining were used to identify the gene candidates, and their functions were studied in vitro and in vivo. A modularized two-step (M2S) technique and CRISPR-Cas9 methods were used to construct engineering yeast. Here, we identified a cytochrome P450, TwCYP82AS1, that catalyses C-22 hydroxylation during wilforlide A biosynthesis. We also found that TwCYP712K1 to K3 can further oxidize the C-29 carboxylation of oleanane-type triterpenes in addition to friedelane-type triterpenes. Reconstitution of the biosynthetic pathway in engineered yeast increased the precursor supply, and combining TwCYP82AS1 and TwCYP712Ks produced abrusgenic acid, which was briefly acidified to achieve the semisynthesis of wilforlide A. Our work presents an alternative metabolic engineering approach for obtaining wilforlide A without relying on extraction from plants.

The haplotype-resolved genome assembly of autotetraploid rhubarb Rheum officinale provides insights into its genome evolution and massive accumulation of anthraquinones

Rheum officinale, a member of the Polygonaceae family, is an important medicinal plant that is widely used in traditional Chinese medicine. Here, we report a 7.68-Gb chromosome-scale assembly of R. officinale with a contig N50 of 3.47 Mb, which was clustered into 44 chromosomes across four homologous groups. Comparative genomics analysis revealed that transposable elements have made a significant contribution to its genome evolution, gene copy number variation, and gene regulation and expression, particularly of genes involved in metabolite biosynthesis, stress resistance, and root development. We placed the recent autotetraploidization of R. officinale at ∼0.58 mya and analyzed the genomic features of its homologous chromosomes. Although no dominant monoploid genomes were observed at the overall expression level, numerous allele-differentially-expressed genes were identified, mainly with different transposable element insertions in their regulatory regions, suggesting that they functionally diverged after polyploidization. Combining genomics, transcriptomics, and metabolomics, we explored the contributions of gene family amplification and tetraploidization to the abundant anthraquinone production of R. officinale, as well as gene expression patterns and differences in anthraquinone content among tissues. Our report offers unprecedented genomic resources for fundamental research on the autopolyploid herb R. officinale and guidance for polyploid breeding of herbs.

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.

Multiomics analyses of two Leonurus species illuminate leonurine biosynthesis and its evolution

The Lamiaceae family is renowned for its terpenoid-based medicinal components, but Leonurus, which has traditional medicinal uses, stands out for its alkaloid-rich composition. Leonurine, the principal active compound found in Leonurus, has demonstrated promising effects in reducing blood lipids and treating strokes. However, the biosynthetic pathway of leonurine remains largely unexplored. Here, we present the chromosome-level genome sequence assemblies of Leonurus japonicus, known for its high leonurine production, and Leonurus sibiricus, characterized by very limited leonurine production. By integrating genomics, RNA sequencing, metabolomics, and enzyme activity assay data, we constructed the leonurine biosynthesis pathway and identified the arginine decarboxylase (ADC), uridine diphosphate glucosyltransferase (UGT), and serine carboxypeptidase-like (SCPL) acyltransferase enzymes that catalyze key reactions in this pathway. Further analyses revealed that the UGT-SCPL gene cluster evolved by gene duplication in the ancestor of Leonurus and neofunctionalization of SCPL in L. japonicus, which contributed to the accumulation of leonurine specifically in L. japonicus. Collectively, our comprehensive study illuminates leonurine biosynthesis and its evolution in Leonurus.

Elucidation of the melitidin biosynthesis pathway in pummelo

Specialized plant metabolism is a rich resource of compounds for drug discovery. The acylated flavonoid glycoside melitidin is being developed as an anti-cholesterol statin drug candidate, but its biosynthetic route in plants has not yet been fully characterized. Here, we describe the gene discovery and functional characterization of a new flavonoid gene cluster (UDP-glucuronosyltransferases (CgUGTs), 1,2 rhamnosyltransferase (Cg1,2RhaT), acyltransferases (CgATs)) that is responsible for melitidin biosynthesis in pummelo (Citrus grandis (L.) Osbeck). Population variation analysis indicated that the tailoring of acyltransferases, specific for bitter substrates, mainly determine the natural abundance of melitidin. Moreover, 3-hydroxy-3-methylglutaryl-CoA reductase enzyme inhibition assays showed that the product from this metabolic gene cluster, melitidin, may be an effective anti-cholesterol statin drug candidate. Co-expression of these clustered genes in Nicotiana benthamiana resulted in the formation of melitidin, demonstrating the potential for metabolic engineering of melitidin in a heterologous plant system. This study establishes a biosynthetic pathway for melitidin, which provides genetic resources for the breeding and genetic improvement of pummelo aimed at fortifying the content of biologically active metabolites.

CYP72D19 from Tripterygium wilfordii catalyzes C-2 hydroxylation of abietane-type diterpenoids

KEY MESSAGE

CYP72D19, the first functional gene of the CYP72D subfamily, catalyzes the C-2 hydroxylation of abietane-type diterpenoids. The abietane-type diterpenoids, e.g., triptolide, tripdiolide, and 2-epitripdiolide, are the main natural products for the anti-tumor, anti-inflammatory, and immunosuppressive activities of Tripterygium wilfordii, while their biosynthetic pathways are not resolved. Here, we cloned and characterized the CYP72D19-catalyzed C-2 hydroxylation of dehydroabietic acid, a compound that has been proven to be a biosynthetic intermediate in triptolide biosynthesis. Through molecular docking and site-directed mutagenesis, L386, L387, and I493 near the active pocket were found to have an important effect on the enzyme activity, which also indicates that steric hindrance of residues plays an important role in function. In addition, CYP72D19 also catalyzed a variety of abietane-type diterpenoids with benzene ring, presumably because the benzene ring of the substrate molecule stabilized the C-ring, allowing the protein and the substrate to form a relatively stable spatial structure. This is the first demonstration of CYP72D subfamily gene function. Our research provides important genetic elements for the structural modification of active ingredients and the heterologous production of other 2-hydroxyl abietane-type natural products.

Comparative genomics of the medicinal plants Lonicera macranthoides and L. japonica provides insight into genus genome evolution and hederagenin-based saponin biosynthesis

Lonicera macranthoides (LM) and L. japonica (LJ) are medicinal plants widely used in treating viral diseases, such as COVID-19. Although the two species are morphologically similar, their secondary metabolite profiles are significantly different. Here, metabolomics analysis showed that LM contained ~86.01 mg/g hederagenin-based saponins, 2000-fold higher than LJ. To gain molecular insights into its secondary metabolite production, a chromosome-level genome of LM was constructed, comprising 9 pseudo-chromosomes with 40 097 protein-encoding genes. Genome evolution analysis showed that LM and LJ were diverged 1.30-2.27 million years ago (MYA). The two plant species experienced a common whole-genome duplication event that occurred ∼53.9-55.2 MYA before speciation. Genes involved in hederagenin-based saponin biosynthesis were arranged in clusters on the chromosomes of LM and they were more highly expressed in LM than in LJ. Among them, oleanolic acid synthase (OAS) and UDP-glycosyltransferase 73 (UGT73) families were much more highly expressed in LM than in LJ. Specifically, LmOAS1 was identified to effectively catalyse the C-28 oxidation of β-Amyrin to form oleanolic acid, the precursor of hederagenin-based saponin. LmUGT73P1 was identified to catalyse cauloside A to produce α-hederin. We further identified the key amino acid residues of LmOAS1 and LmUGT73P1 for their enzymatic activities. Additionally, comparing with collinear genes in LJ, LmOAS1 and LmUGT73P1 had an interesting phenomenon of 'neighbourhood replication' in LM genome. Collectively, the genomic resource and candidate genes reported here set the foundation to fully reveal the genome evolution of the Lonicera genus and hederagenin-based saponin biosynthetic pathway.

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.

Multi-omics analysis reveals the evolutionary origin of diterpenoid alkaloid biosynthesis pathways in Aconitum

Diterpenoid alkaloids (DAs) have been often utilized in clinical practice due to their analgesic and anti-inflammatory properties. Natural DAs are prevalent in the family Ranunculaceae, notably in the Aconitum genus. Nevertheless, the evolutionary origin of the biosynthesis pathway responsible for DA production remains unknown. In this study, we successfully assembled a high-quality, pseudochromosome-level genome of the DA-rich species Aconitum vilmorinianum (A. vilmorinianum) (5.76 Gb). An A. vilmorinianum-specific whole-genome duplication event was discovered using comparative genomic analysis, which may aid in the evolution of the DA biosynthesis pathway. We identified several genes involved in DA biosynthesis via integrated genomic, transcriptomic, and metabolomic analyses. These genes included enzymes encoding target ent-kaurene oxidases and aminotransferases, which facilitated the activation of diterpenes and insertion of nitrogen atoms into diterpene skeletons, thereby mediating the transformation of diterpenes into DAs. The divergence periods of these genes in A. vilmorinianum were further assessed, and it was shown that two major types of genes were involved in the establishment of the DA biosynthesis pathway. Our integrated analysis offers fresh insights into the evolutionary origin of DAs in A. vilmorinianum as well as suggestions for engineering the biosynthetic pathways to obtain desired DAs.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

Biosynthesis and biotechnological production of the anti-obesity agent celastrol

Obesity is a major health risk still lacking effective pharmacological treatment. A potent anti-obesity agent, celastrol, has been identified in the roots of Tripterygium wilfordii. However, an efficient synthetic method is required to better explore its biological utility. Here we elucidate the 11 missing steps for the celastrol biosynthetic route to enable its de novo biosynthesis in yeast. First, we reveal the cytochrome P450 enzymes that catalyse the four oxidation steps that produce the key intermediate celastrogenic acid. Subsequently, we show that non-enzymatic decarboxylation-triggered activation of celastrogenic acid leads to a cascade of tandem catechol oxidation-driven double-bond extension events that generate the characteristic quinone methide moiety of celastrol. Using this acquired knowledge, we have developed a method for producing celastrol starting from table sugar. This work highlights the effectiveness of combining plant biochemistry with metabolic engineering and chemistry for the scalable synthesis of complex specialized metabolites.

Pharmacological inhibition of demethylzeylasteral on JAK-STAT signaling ameliorates vitiligo

BACKGROUND

The activation of CD8⁺ T cells and their trafficking to the skin through JAK-STAT signaling play a central role in the development of vitiligo. Thus, targeting this key disease pathway with innovative drugs is an effective strategy for treating vitiligo. Natural products isolated from medicinal herbs are a useful source of novel therapeutics. Demethylzeylasteral (T-96), extracted from Tripterygium wilfordii Hook F, possesses immunosuppressive and anti-inflammatory properties.

METHODS

The efficacy of T-96 was tested in our mouse model of vitiligo, and the numbers of CD8⁺ T cells infiltration and melanocytes remaining in the epidermis were quantified using whole-mount tail staining. Immune regulation of T-96 in CD8⁺ T cells was evaluated using flow cytometry. Pull-down assay, mass spectrum analysis, molecular docking, knockdown and overexpression approaches were utilized to identify the target proteins of T-96 in CD8⁺ T cells and keratinocytes.

RESULTS

Here, we found that T-96 reduced CD8⁺ T cell infiltration in the epidermis using whole-mount tail staining and alleviated the extent of depigmentation to a comparable degree of tofacitinib (Tofa) in our vitiligo mouse model. In vitro, T-96 decreased the proliferation, CD69 membrane expression, and IFN-γ, granzyme B, (GzmB), and perforin (PRF) levels in CD8⁺ T cells isolated from patients with vitiligo. Pull-down assays combined with mass spectrum analysis and molecular docking showed that T-96 interacted with JAK3 in CD8⁺ T cell lysates. Furthermore, T-96 reduced JAK3 and STAT5 phosphorylation following IL-2 treatment. T-96 could not further reduce IFN-γ, GzmB and PRF expression following JAK3 knockdown or inhibit increased immune effectors expression upon JAK3 overexpression. Additionally, T-96 interacted with JAK2 in IFN-γ-stimulated keratinocytes, inhibiting the activation of JAK2, decreasing the total and phosphorylated protein levels of STAT1, and reducing the production and secretion of CXCL9 and CXCL10. T-96 did not significantly inhibit STAT1 and CXCL9/10 expression following JAK2 knockdown, nor did it suppress upregulated STAT1-CXCL9/10 signaling upon JAK2 overexpression. Finally, T-96 reduced the membrane expression of CXCR3, and the culture supernatants pretreated with T-96 under IFN-γ stressed keratinocytes markedly blocked the migration of CXCR3⁺CD8⁺ T cells, similarly to Tofa in vitro.

CONCLUSION

Our findings demonstrated that T-96 might have positive therapeutic responses to vitiligo by pharmacologically inhibiting the effector functions and skin trafficking of CD8⁺ T cells through JAK-STAT signaling.

Functional Study and Efficient Catalytic Element Mining of CYP76AHs in Salvia Plants

Salvia is a large genus with hundreds of species used in traditional Chinese medicine. Tanshinones are a highly representative class of exclusive compounds found in the Salvia genus that exhibit significant biological activity. Tanshinone components have been identified in 16 Salvia species. The CYP76AH subfamily (P450) is crucial for the synthesis of tanshinone due to its catalytic generation of polyhydroxy structures. In this study, a total of 420 CYP76AH genes were obtained, and phylogenetic analysis showed their clear clustering relationships. Fifteen CYP76AH genes from 10 Salvia species were cloned and studied from the perspectives of evolution and catalytic efficiency. Three CYP76AHs with significantly improved catalytic efficiency compared to SmCYP76AH3 were identified, providing efficient catalytic elements for the synthetic biological production of tanshinones. A structure-function relationship study revealed several conserved residues that might be related to the function of CYP76AHs and provided a new mutation direction for the study of the directed evolution of plant P450.

Haplotype-resolved genome assembly of Coriaria nepalensis a non-legume nitrogen-fixing shrub

Coriaria nepalensis Wall. (Coriariaceae) is a nitrogen-fixing shrub which forms root nodules with the actinomycete Frankia. Oils and extracts of C. nepalensis have been reported to be bacteriostatic and insecticidal, and C. nepalensis bark provides a valuable tannin resource. Here, by combining PacBio HiFi sequencing and Hi-C scaffolding techniques, we generated a haplotype-resolved chromosome-scale genome assembly for C. nepalensis. This genome assembly is approximately 620 Mb in size with a contig N50 of 11 Mb, with 99.9% of the total assembled sequences anchored to 40 pseudochromosomes. We predicted 60,862 protein-coding genes of which 99.5% were annotated from databases. We further identified 939 tRNAs, 7,297 rRNAs, and 982 ncRNAs. The chromosome-scale genome of C. nepalensis is expected to be a significant resource for understanding the genetic basis of root nodulation with Frankia, toxicity, and tannin biosynthesis.

Identification and functional characterization of three cytochrome P450 genes for the abietane diterpenoid biosynthesis in Isodon lophanthoides

We identify two ferruginol synthases and a 11-hydroxyferruginol synthase from a traditional Chinese medicinal herb Isodon lophanthoides and propose their involvement in two independent abietane diterpenoids biosynthetic pathways. Isodon lophanthoides is a traditional Chinese medicinal herb rich in highly oxidized abietane-type diterpenoids. These compounds exhibit a wide range of pharmaceutical activities, yet the biosynthesis is barely known. Here, we describe the screening and functional characterization of P450s that oxidize the abietane skeleton abietatriene. We mainly focused on CYP76 family and identified 12 CYP76AHs by mining the RNA-seq data of I. lophanthoides. Among the 12 CYP76AHs, 6 exhibited similar transcriptional expression features as upstream diterpene synthases, including root or leaf-preferential expression pattern and highly MeJA inducibility. These six P450s were considered as first-tier candidates and functionally characterized in yeast and plant cells. In yeast assays showed that both CYP76AH42 and CYP76AH43 were ferruginol synthases hydroxylating the C12 position of abietatriene, whereas CYP76AH46 was characterized as a 11-hydroxyferruginol synthase which catalyzes two successive oxidations at C12 and C11 of abietatriene. Heterologous expression of three CYP76AHs in Nicotiana benthamiana resulted in the formation of ferruginol. qPCR analysis showed CYP76AH42 and CYP76AH43 were mainly expressed in the root, which was consistent with the distribution of ferruginol in the root periderms. CYP76AH46 was primarily expressed in the leaves where barely ferruginol or 11-hydroxyferruginol was detected. In addition to distinct organ-specific expression pattern, three CYP76AHs exhibited different genomic structures (w or w/o introns), low protein sequence identities (51-63%) and were placed in separate subclades in the phylogenetic tree. These results suggest that the identified CYP76AHs may be involved in at least two independent abietane biosynthetic pathways in the aerial and underground parts of I. lophanthoides.

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.

Efficient heterologous expression of cytochrome P450 enzymes in microorganisms for the biosynthesis of natural products

Natural products, a chemically and structurally diverse class of molecules, possess a wide spectrum of biological activities, have been used therapeutically for millennia, and have provided many lead compounds for the development of synthetic drugs. Cytochrome P450 enzymes (P450s, CYP) are widespread in nature and are involved in the biosynthesis of many natural products. P450s are heme-containing enzymes that use molecular oxygen and the hydride donor NAD(P)H (coupled via enzymic redox partners) to catalyze the insertion of oxygen into C-H bonds in a regio- and stereo-selective manner, effecting hydroxylation and several other reactions. With the rapid development of systems biology, numerous novel P450s have been identified for the biosynthesis of natural products, but there are still several challenges to the efficient heterologous expression of active P450s. This review covers recent developments in P450 research and development, including the properties and functions of P450s, discovery and mining of novel P450s, modification and screening of P450 mutants, improved heterologous expression of P450s in microbial hosts, efficient whole-cell transformation with P450s, and current applications of P450s for the biosynthesis of natural products. This resource provides a solid foundation for the application of highly active and stable P450s in microbial cell factories to biosynthesize natural products.

Genome-wide analyzation and functional characterization on the TPS family provide insight into the biosynthesis of mono-terpenes in the camphor tree

Terpene synthase (TPS) plays an important role in terpenoids biosynthesis. Cinnamomum camphora (camphor tree) contains dozens of terpenoids with medicinal value, especially borneol, which has been widely used since ancient times. However, limited information is available regarding the genome-wide identification and characterization of the TPS family in the C. camphora. In this study, 82 CcTPS genes were identified from the camphor tree genome (CTG). Gene cluster and sequence syntenic analysis suggested that tandem duplication occurred within the TPS family of the CTG, especially for the TPS-b subfamily. The chemotype-specific gene expression analysis showed significantly differential expression patterns among six chemotypes. It is worth noting that three genes (CcTPS26, CcTPS49 and CcTPS72) exhibited relatively high expression in the borneol-type camphor tree, compared to the other five chemotypes. Further functional characterization of them indicated that they were all bornyl diphosphate synthases (BPPSs), which function in catalyzing GPP into BPP and then undergoes dephosphorylation to yield borneol. This is the first report that multiple BPPSs exist within a single species. Intriguingly, CcTPS49 and CcTPS72 lead to the generation of dextral-borneol, while CcTPS26 contributes to the biosynthesis of levo-borneol. In addition, the functional characterization of another six CcTPSs suggested that they are responsible for the biosynthesis of linalool, eucalyptol and several other monoterpenes in camphor tree. In conclusion, these novel results provide a foundation for further exploration of the role of the CcTPS gene family and shed light on a better understanding of the biosynthesis and accumulation of monoterpenes in camphor tree.

Tandemly duplicated CYP82Ds catalyze 14-hydroxylation in triptolide biosynthesis and precursor production in Saccharomyces cerevisiae

Triptolide is a valuable multipotent antitumor diterpenoid in Tripterygium wilfordii, and its C-14 hydroxyl group is often selected for modification to enhance both the bioavailability and antitumor efficacy. However, the mechanism for 14-hydroxylation formation remains unknown. Here, we discover 133 kb of tandem duplicated CYP82Ds encoding 11 genes on chromosome 12 and characterize CYP82D274 and CYP82D263 as 14-hydroxylases that catalyze the metabolic grid in triptolide biosynthesis. The two CYP82Ds catalyze the aromatization of miltiradiene, which has been repeatedly reported to be a spontaneous process. In vivo assays and evaluations of the kinetic parameters of CYP82Ds indicate the most significant affinity to dehydroabietic acid among multiple intermediates. The precursor 14-hydroxy-dehydroabietic acid is successfully produced by engineered Saccharomyces cerevisiae. Our study provides genetic elements for further elucidation of the downstream biosynthetic pathways and heterologous production of triptolide and of the currently intractable biosynthesis of other 14-hydroxyl labdane-type secondary metabolites.

Cytochrome P450s in plant terpenoid biosynthesis: discovery, characterization and metabolic engineering

As the largest family of natural products, terpenoids play valuable roles in medicine, agriculture, cosmetics and food. However, the traditional methods that rely on direct extraction from the original plants not only produce low yields, but also result in waste of resources, and are not applicable at all to endangered species. Modern heterologous biosynthesis is considered a promising, efficient, and sustainable production method, but it relies on the premise of a complete analysis of the biosynthetic pathway of terpenoids, especially the functionalization processes involving downstream cytochrome P450s. In this review, we systematically introduce the biotech approaches used to discover and characterize plant terpenoid-related P450s in recent years. In addition, we propose corresponding metabolic engineering approaches to increase the effective expression of P450 and improve the yield of terpenoids, and also elaborate on metabolic engineering strategies and examples of heterologous biosynthesis of terpenoids in Saccharomyces cerevisiae and plant hosts. Finally, we provide perspectives for the biotech approaches to be developed for future research on terpenoid-related P450.

Transcriptome analysis reveals that jasmonic acid biosynthesis and signaling is associated with the biosynthesis of asperosaponin VI in Dipsacus asperoides

Dipsacus asperoides is a perennial herb, the roots of which are abundant in asperosaponin VI, which has important medicinal value. However, the molecular mechanism underlying the biosynthesis of asperosaponin VI in D. asperoides remains unclear. In present study, a comprehensive investigation of asperosaponin VI biosynthesis was conducted at the levels of metabolite and transcript during root development. The content of asperosaponin VI was significantly accumulated in two-leaf stage roots, and the spatial distribution of asperosaponin VI was localized in the xylem. The concentration of asperosaponin VI gradually increased in the root with the development process. Transcriptome analysis revealed 3916 unique differentially expressed genes (DEGs) including 146 transcription factors (TFs) during root development in D. asperoides. In addition, α-linolenic acid metabolism, jasmonic acid (JA) biosynthesis, JA signal transduction, sesquiterpenoid and triterpenoid biosynthesis, and terpenoid backbone biosynthesis were prominently enriched. Furthermore, the concentration of JA gradually increased, and genes involved in α-linolenic acid metabolism, JA biosynthesis, and triterpenoid biosynthesis were up-regulated during root development. Moreover, the concentration of asperosaponin VI was increased following methyl jasmonate (MeJA) treatment by activating the expression of genes in the triterpenoid biosynthesis pathway, including acetyl-CoA acetyltransferase (DaAACT), 3-hydroxy-3-methylglutaryl coenzyme A synthase (DaHMGCS), 3-hydroxy-3-methylglutaryl coenzyme-A reductase (DaHMGCR). We speculate that JA biosynthesis and signaling regulates the expression of triterpenoid biosynthetic genes and facilitate the biosynthesis of asperosaponin VI. The results suggest a regulatory network wherein triterpenoids, JA, and TFs co-modulate the biosynthesis of asperosaponin VI in D. asperoides.

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.

Chromosome-level genome assembly and functional characterization of terpene synthases provide insights into the volatile terpenoid biosynthesis of Wurfbainia villosa

Wurfbainia villosa is a well-known medicinal and edible plant that is widely cultivated in the Lingnan region of China. Its dried fruits (called Fructus Amomi) are broadly used in traditional Chinese medicine for curing gastrointestinal diseases and are rich in volatile terpenoids. Here, we report a high-quality chromosome-level genome assembly of W. villosa with a total size of approximately 2.80 Gb, 42 588 protein-coding genes, and a very high percentage of repetitive sequences (87.23%). Genome analysis showed that W. villosa likely experienced a recent whole-genome duplication event prior to the W. villosa-Zingiber officinale divergence (approximately 11 million years ago), and a recent burst of long terminal repeat insertions afterward. The W. villosa genome enabled the identification of 17 genes involved in the terpenoid skeleton biosynthesis pathway and 66 terpene synthase (TPS) genes. We found that tandem duplication events have an important contribution to the expansion of WvTPSs, which likely drove the production of volatile terpenoids. In addition, functional characterization of 18 WvTPSs, focusing on the TPS-a and TPS-b subfamilies, showed that most of these WvTPSs are multi-product TPS and are predominantly expressed in seeds. The present study provides insights into the genome evolution and the molecular basis of the volatile terpenoids diversity in W. villosa. The genome sequence also represents valuable resources for the functional gene research and molecular breeding of W. villosa.

Increased genetic diversity loss and genetic differentiation in a model marine diatom adapted to ocean warming compared to high CO2

Although high CO2 and warming could act interactively on marine phytoplankton, little is known about the molecular basis for this interaction on an evolutionary scale. Here we explored the adaptation to high CO2 in combination with warming in a model marine diatom Phaeodactylum tricornutum. Whole-genome re-sequencing identifies, in comparison to populations grown under control conditions, a larger genetic diversity loss and a higher genetic differentiation in the populations adapted for 2 years to warming than in those adapted to high CO2. However, this diversity loss was less under high CO2 combined with warming, suggesting that the evolution driven by warming was constrained by high CO2. By integrating genomics, transcriptomics, and physiological data, we found that the underlying molecular basis for this constraint is associated with the expression of genes involved in some key metabolic pathways or biological processes, such as the glyoxylate pathway, amino acid and fatty acid metabolism, and diel variability. Our results shed new light on the evolutionary responses of marine phytoplankton to multiple environmental changes in the context of global change and provide new insights into the molecular basis underpinning interactions among those multiple drivers.

The Chromosome-level genome of Aesculus wilsonii provides new insights into terpenoid biosynthesis and Aesculus evolution

Aesculus L. (buckeye and horse chestnut) are woody plant species with important horticultural and medicinal values. Aesculus seeds are widely used as biomedicine and cosmetic ingredients due to their saponins. We report a chromosomal-scale genome of Aesculus wilsonii. Sequences amounting to a total of 579.01 Mb were assembled into 20 chromosomes. More than half of the genome (54.46%) were annotated as repetitive sequences, and 46,914 protein-coding genes were predicted. In addition to the widespread gamma event with core eudicots, a unique whole-genome duplication (WGD) event (17.69 Mya) occurred in Aesculus after buckeye differentiated from longan. Due to WGD events and tandem duplications, the related synthetic genes of triterpene saponins unique to Aesculus increased significantly. Combined with transcriptome characterization, the study preliminarily resolved the biosynthetic pathway of triterpenoid saponins like aescin in A. wilsonii genome. Analyses of the resequencing of 104 buckeye accessions revealed clear relationship between the geographic distribution and genetic differentiation of buckeye trees in China. We found that the buckeye species found in southern Shaanxi is A. wilsonii rather than A. chinensis. Population dynamics analysis further suggests that the population size and evolution of existing buckeye species have been influenced by climate fluctuations during the Pleistocene and recent domestication events. The genome of A. wilsonii and population genomics of Aesculus provide a resource for future research on Hippocastanaceae. These findings will contribute to the utilization and diversity protection of Aesculus.

Achievements and perspectives of synthetic biology in botanical insecticides

Botanical insecticides are the origin of all insecticidal compounds. They have been widely used to control pests in crops for a long time. Currently, the commercial production of botanical insecticides extracted from plants is limited because of insufficient raw material supply. Synthetic biology is a promising and effective approach for addressing the current problems of the production of botanical insecticides. It is an emerging biological research hotspot in the field of botanical insecticides. However, the biosynthetic pathways of many botanical insecticides are not completely elucidated. On the other hand, the cytotoxicity of botanical pesticides and low efficiency of these biosynthetic enzymes in new hosts make it still challenging for their heterologous production. In the present review, we summarized the recent developments in the heterologous production of botanical insecticides, analyzed the current challenges, and discussed the feasible production strategies, focusing on elucidating biosynthetic pathways, enzyme engineering, host engineering, and cytotoxicity engineering. Looking to the future, synthetic biology promises to further advance heterologous production of more botanical pesticides.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Haplotype-resolved genome assembly of Bletilla striata (Thunb.) Reichb.f. to elucidate medicinal value

Bletilla striata, commonly known as baiji, is a species used in traditional Chinese medicine; it is highly regarded for its medicinal applications and therefore has high economic value. Here, we report a high-quality haplotype-resolved genome of B. striata, haplotype A (2.37 Gb, with a scaffold N50 of 146.39 Mb and a contig N50 of 1.65 Mb) and haplotype B (2.43 Gb, with a scaffold N50 of 150.22 Mb and a contig N50 of 1.66 Mb), assembled from high-fidelity (HiFi) reads and chromosome conformation capture (Hi-C) reads. We find evidence that B. striata has undergone two whole-genome duplication (WGD) events: an ancient WGD event shared by most monocots and a recent WGD event unique to all orchids. We also reconstructed the ancestral orchid karyotype (AOK) of 18 ancient chromosomes and the evolutionary trajectories of 16 modern B. striata chromosomes. Comparative genomic analysis suggests that the expanded gene families of B. striata might play important roles in secondary metabolite biosynthesis and environmental adaptation. By combining genomic and transcriptomic data, we identified the 10 core members from nine gene families that were probably involved in B. striata polysaccharide (BSP) biosynthesis. Based on virus-induced gene silencing (VIGS) and yeast two-hybrid experiments, we present an MYB transcription factor (TF), BsMYB2, that can regulate BSP biosynthesis by directly interacting with eight key BSP-related genes: sacA1, HK1, scrK1, scrK2, GPI1, manA1, GMPP1 and UGP2_1. Our study will enhance the understanding of orchid evolution and accelerate the molecular-assisted breeding of B. striata for improving traits of medicinal value.

Ancestral sequence reconstruction of the CYP711 family reveals functional divergence in strigolactone biosynthetic enzymes associated with gene duplication events in monocot grasses

The strigolactone (SL) class of phytohormones shows broad chemical diversity, the functional importance of which remains to be fully elucidated, along with the enzymes responsible for the diversification of the SL structure. Here we explore the functional evolution of the highly conserved CYP711A P450 family, members of which catalyze several key monooxygenation reactions in the strigolactone pathway. Ancestral sequence reconstruction was utilized to infer ancestral CYP711A sequences based on a comprehensive set of extant CYP711 sequences. Eleven ancestral enzymes, corresponding to key points in the CYP711A phylogenetic tree, were resurrected and their activity was characterized towards the native substrate carlactone and the pure enantiomers of the synthetic strigolactone analogue, GR24. The ancestral and extant CYP711As tested accepted GR24 as a substrate and catalyzed several diversifying oxidation reactions on the structure. Evidence was obtained for functional divergence in the CYP711A family. The monocot group 3 ancestor, arising from gene duplication events within monocot grasses, showed both increased catalytic activity towards GR24 and high stereoselectivity towards the GR24 isomer resembling strigol-type SLs. These results are consistent with a role for CYP711As in strigolactone diversification in early land plants, which may have extended to the diversification of strigol-type SLs.

Tripterygium wilfordii cytochrome P450s catalyze the methyl shift and epoxidations in the biosynthesis of triptonide

The diterpenoid triepoxides triptolide and triptonide from Tripterygium wilfordii (thunder god wine) exhibit unique bioactivities with potential uses in disease treatment and as a non-hormonal male contraceptives. Here, we show that cytochrome P450s (CYPs) from the CYP71BE subfamily catalyze an unprecedented 18(4→3) methyl shift required for biosynthesis of the abeo-abietane core structure present in diterpenoid triepoxides and in several other plant diterpenoids. In combination with two CYPs of the CYP82D subfamily, four CYPs from T. wilfordii are shown to constitute the minimal set of biosynthetic genes that enables triptonide biosynthesis using Nicotiana benthamiana and Saccharomyces cerevisiae as heterologous hosts. In addition, co-expression of a specific T. wilfordii cytochrome b₅ (Twcytb₅-A) increases triptonide output more than 9-fold in S. cerevisiae and affords isolation and structure elucidation by NMR spectroscopic analyses of 18 diterpenoids, providing insights into the biosynthesis of diterpenoid triepoxides. Our findings pave the way for diterpenoid triepoxide production via fermentation.

How triptonide is made in the medicinal plant Tripterygium wilfordii is largely unknown. Here, the authors report the identification and characterization of a suite of cytochrome P450s and show their function in catalyzing the formation of triptonide from miltriadiene in tobacco and baker’s yeast.

Genomic insights into the evolution of plant chemical defense

Plant trait evolution can be impacted by common mechanisms of genome evolution, including whole-genome and small-scale duplication, rearrangement, and selective pressures. With the increasing accessibility of genome sequencing for non-model species, comparative studies of trait evolution among closely related or divergent lineages have supported investigations into plant chemical defense. Plant defensive compounds include major chemical classes, such as terpenoids, alkaloids, and phenolics, and are used in primary and secondary plant functions. These include the promotion of plant health, facilitation of pollination, defense against pathogens, and responses to a rapidly changing climate. We discuss mechanisms of genome evolution and use examples from recent studies to impress a stronger understanding of the link between genotype and phenotype as it relates to the evolution of plant chemical defense. We conclude with considerations for how to leverage genomics, transcriptomics, metabolomics, and functional assays for studying the emergence and evolution of chemical defense systems.

Biosynthesis, total synthesis, and pharmacological activities of aryltetralin-type lignan podophyllotoxin and its derivatives

Covering: up to 2022Podophyllotoxin (PTOX, 1), a kind of aryltetralin-type lignan, was first discovered in the plant Podophyllum peltatum and its structure was clarified by W. Borsche and J. Niemann in 1932. Due to its potent anti-cancer and anti-viral activities, it is considered one of the molecules most likely to be developed into modern drugs. With the increasing market demand and insufficient storage of natural resources, it is crucial to expand the sources of PTOXs. The original extraction method from plants has gradually failed to meet the requirements, and the biosynthesis and total synthesis have become the forward-looking alternatives. As key enzymes in the biosynthetic pathway of PTOXs and their catalytic mechanisms being constantly revealed, it is possible to realize the heterogeneous biosynthesis of PTOXs in the future. Chemical and chemoenzymatic synthesis also provide schemes for strictly controlling the asymmetric configuration of the tetracyclic core. Currently, the pharmacological activities of some PTOX derivatives have been extensively studied, laying the foundation for clinical candidate drugs. This review focuses primarily on the latest research progress in the biosynthesis, total synthesis, and pharmacological activities of PTOX and its derivatives, providing a more comprehensive understanding of these widely used compounds and supporting the future search for clinical applications.

Plant biosynthetic gene clusters in the context of metabolic evolution

Covering: up to 2022Plants produce a wide range of structurally and biosynthetically diverse natural products to interact with their environment. These specialised metabolites typically evolve in limited taxonomic groups presumably in response to specific selective pressures. With the increasing availability of sequencing data, it has become apparent that in many cases the genes encoding biosynthetic enzymes for specialised metabolic pathways are not randomly distributed on the genome. Instead they are physically linked in structures such as arrays, pairs and clusters. The exact function of these clusters is debated. In this review we take a broad view of gene arrangement in plant specialised metabolism, examining types of structures and variation. We discuss the evolution of biosynthetic gene clusters in the wider context of metabolism, populations and epigenetics. Finally, we synthesise our observations to propose a new hypothesis for biosynthetic gene cluster formation in plants.

A Systematic Review and Meta-Analysis of the Inhibitory Effects of Plant-Derived Sterilants on Rodent Population Abundance

Owing to their low minimal environmental risk and other ethical considerations, plant-derived sterilants are used to control rodent populations. However, the effects of plant-derived sterilants are not immediate, and their efficacy on rodent control is controversial, which negatively affects sterilant research and application. Here, a meta-analysis of the available literature was conducted to evaluate the effects of two plant-derived sterilants, triptolide and curcumol, on rodent populations. Using a random-effects and a fixed-effects model, we calculated the weighted mean difference (WMD) and relative risk (RR) and their corresponding 95% confidence intervals (95% CIs). After the application of plant-derived sterilants, the rodent population density tended to decrease. Three outcome-related measures in rodents, i.e., capture rate (RR = 0.31, 95% CI [0.20, 0.47]), pregnancy rate (RR = 0.49, 95% CI [0.40, 0.61]), and sperm survival rate (WMD = −17.53, 95% CI [−28.96, −6.06]), significantly decreased, as shown by a significant reduction of ovarian, uterine, and testicular organ coefficients. However, the number of effective rodent holes did not change significantly after the interventions, indicating that the studied sterilants did not directly eradicate the rodent populations. This study provides a theoretical basis for elucidating the inhibitory mechanisms of plant-derived sterilants on rodent populations and for the rational use of these sterilants.

Genome-wide analysis of MYB family genes in Tripterygium wilfordii and their potential roles in terpenoid biosynthesis

Terpenoids are a class of significant bioactive components in the woody vine of Tripterygium wilfordii. Previous studies have shown that MYB transcription factors play important roles in plant secondary metabolism, growth, and developmental processes. However, the MYB involved in terpenoid biosynthesis in Tripterygium wilfordii are unknown. To identify Tripterygium wilfordii MYB (TwMYB) genes that are involved in terpenoid biosynthesis, we conducted the genome-wide analysis of the TwMYB gene family. A total of 207 TwMYBs were identified including 84 1R-TwMYB, 117 R2R3-TwMYB, four 3R-TwMYB, and two 4R-TwMYB genes. The most abundant R2R3-TwMYBs together with their Arabidopsis homologs were categorized into 26 subgroups. Intraspecific collinearity analysis found that the 74.9% of the TwMYBs may be generated by segmental duplication events, and 36.7% of duplicated gene pairs were derived from the specific whole genome duplication (WGD) event in Tripterygium wilfordii. In addition, interspecies collinearity analysis found that 16 TwMYB genes formed homologous gene pairs with MYB genes in seven representative species, which indicated they may have a key role in evolution. Notably, we found that the TwMYB genes were differentially expressed in various tissues by expression pattern analysis. In order to further select the candidate genes related to terpenoid biosynthesis, the assay of Methyl jasmonate (MeJA) induction and analysis of phylogenetic tree was conducted. It was speculated that six candidate TwMYB genes (TwMYB33, TwMYB34, TwMYB45, TwMYB67, TwMYB102, and TwMYB103) are involved in regulating terpenoid biosynthesis. This study is the first systematic analysis of the TwMYB gene family and will lay a foundation for the functional characterization of TwMYB genes in the regulation of terpenoid biosynthesis.

Genome sequence and population genomics provide insights into chromosomal evolution and phytochemical innovation of Hippophae rhamnoides

Plants of the Elaeagnaceae family are widely used to treat various health disorders owing to their natural phytochemicals. Seabuckthorn (Hippophae rhamnoides L.) is an economically and ecologically important species within the family with richness of biologically and pharmacologically active substances. Here, we present a chromosome-level genome assembly of seabuckthorn (http://hipp.shengxin.ren/), the first genome sequence of Elaeagnaceae, which has a total length of 849.04 Mb with scaffold N50 of 69.52 Mb and 30 864 annotated genes. Two sequential tetraploidizations with one occurring ~36-41 million years ago (Mya) and the last ~24-27 Mya were inferred, resulting in expansion of genes related to ascorbate and aldarate metabolism, lipid biosynthesis, and fatty acid elongation. Comparative genomic analysis reconstructed the evolutionary trajectories of the seabuckthorn genome with the predicted ancestral genome of 14 proto-chromosomes. Comparative transcriptomic and metabonomic analyses identified some key genes contributing to high content of polyunsaturated fatty acids and ascorbic acid (AsA). Additionally, we generated and analysed 55 whole-genome sequences of diverse accessions, and identified 9.80 million genetic variants in the seabuckthorn germplasms. Intriguingly, genes in selective sweep regions identified through population genomic analysis appeared to contribute to the richness of AsA and fatty acid in seabuckthorn fruits, among which GalLDH, GMPase and ACC, TER were the potentially major-effect causative genes controlling AsA and fatty acid content of the fruit, respectively. Our research offers novel insights into the molecular basis underlying phytochemical innovation of seabuckthorn, and provides valuable resources for exploring the evolution of the Elaeagnaceae family and molecular breeding.

The Manchurian Walnut Genome: Insights into Juglone and Lipid Biosynthesis

BACKGROUND

Manchurian walnut (Juglans mandshurica Maxim.) is a tree with multiple industrial uses and medicinal properties in the Juglandaceae family (walnuts and hickories). J. mandshurica produces juglone, which is a toxic allelopathic agent and has potential utilization value. Furthermore, the seed of J. mandshurica is rich in various unsaturated fatty acids and has high nutritive value.

FINDINGS

Here, we present a high-quality chromosome-scale reference genome assembly and annotation for J. mandshurica (n = 16) with a contig N50 of 21.4 Mb by combining PacBio high-fidelity reads with high-throughput chromosome conformation capture data. The assembled genome has an estimated sequence size of 548.7 Mb and consists of 657 contigs, 623 scaffolds, and 40,453 protein-coding genes. In total, 60.99% of the assembled genome consists of repetitive sequences. Sixteen super-scaffolds corresponding to the 16 chromosomes were assembled, with a scaffold N50 length of 33.7 Mb and a BUSCO complete gene percentage of 98.3%. J. mandshurica displays a close sequence relationship with Juglans cathayensis, with a divergence time of 13.8 million years ago. Combining the high-quality genome, transcriptome, and metabolomics data, we constructed a gene-to-metabolite network and identified 566 core and conserved differentially expressed genes, which may be involved in juglone biosynthesis. Five CYP450 genes were found that may contribute to juglone accumulation. NAC, bZip, NF-YA, and NF-YC are positively correlated with the juglone content. Some candidate regulators (e.g., FUS3, ABI3, LEC2, and WRI1 transcription factors) involved in the regulation of lipid biosynthesis were also identified.

CONCLUSIONS

Our genomic data provide new insights into the evolution of the walnut genome and create a new platform for accelerating molecular breeding and improving the comprehensive utilization of these economically important tree species.

Integrated Transcriptomic and Metabolomic Analyses Reveal the Mechanisms Underlying Anthocyanin Coloration and Aroma Formation in Purple Fennel

The color and aroma are the significant traits of vegetables and fruits, but the metabolic and molecular mechanisms underlying anthocyanin accumulation and aroma formation remain almost unknown in fennel (Anethum foeniculum L.), which is a crucial vegetable crop and grown widely for aromatic leaves and bulbs. Here, ten major anthocyanins identified and quantified by ultra-high performance liquid chromatography coupled with quadrupole Orbitrap high-resolution mass spectrometry (UHPLC-Q-Orbitrap HRMS) were mainly responsible for the coloration of purple fennel leaf. With the application of GC-MS, it was found that the reduced volatile phenylpropanoids including isoeugenol, trans-isoeugenol, and apiol chiefly account for the characteristic aroma changes of the purple fennel. Moreover, the characteristic anthocyanin coloration and aroma formation in purple fennel were systematically studied with the integrated transcriptomics and metabolomics. The critical genes associated with the biosynthesis and regulation of anthocyanins and volatile phenylpropanoids were isolated and studied carefully in transiently transfected tobacco cells and transgenic tomato plants. Together with the results of UHPLC-Q-Orbitrap HRMS, RT-qPCR, and yeast two hybrid (Y2H), it is proved that the metabolic flux redirection of phenylpropanoid pathway primarily regulated by a functional MYB-bHLH-WD40 complex consisting of AfTT8, AfMYB7, and AfTTG1 accounts for the characteristic anthocyanin coloration and aroma formation in purple fennel leaf. The systematic understanding of the anthocyanin accumulation and aroma formation will assist in the improvement of fennel resource utilization and breeding.

Research Progress on the Synthetic Biology of Botanical Biopesticides

The production and large-scale application of traditional chemical pesticides will bring environmental pollution and food safety problems. With the advantages of high safety and environmental friendliness, botanical biopesticides are in line with the development trend of modern agriculture and have gradually become the mainstream of modern pesticide development. However, the traditional production of botanical biopesticides has long been faced with prominent problems, such as limited source and supply, complicated production processes, and excessive consumption of resources. In recent years, the rapid development of synthetic biology will break through these bottlenecks, and many botanical biopesticides are produced using synthetic biology, such as emodin, celangulin, etc. This paper reviews the latest progress and application prospect of synthetic biology in the development of botanical pesticides so as to provide new ideas for the analysis of synthetic pathways and heterologous and efficient production of botanical biopesticides and accelerate the research process of synthetic biology of natural products.