Plant metabolic clusters - from genetics to genomics

Hunting for the Intermolecular Diels-Alderase

ConspectusThe Diels-Alder reaction is well known as a concerted [4 + 2] cycloaddition governed by the Woodward-Hoffmann rules. Since Prof. Otto Diels and his student Kurt Alder initially reported the intermolecular [4 + 2] cycloaddition between cyclopentadiene and quinone in 1928, it has been recognized as one of the most powerful chemical transformations to build C-C bonds and construct cyclic structures. This named reaction has been widely used in synthesizing natural products and drug molecules. Driven by the synthetic importance of the Diels-Alder reaction, identifying the enzyme that stereoselectively catalyzes the Diels-Alder reaction has become an intriguing research area in natural product biosynthesis and biocatalysis. With significant progress in sequencing and bioinformatics, dozens of Diels-Alderases have been characterized in microbial natural product biosynthesis. However, few are evolutionally dedicated to catalyzing an intermolecular Diels-Alder reaction with a concerted mechanism.This Account summarizes our endeavors to hunt for the naturally occurring intermolecular Diels-Alderase from plants. Our research journey started from the biomimetic syntheses of D-A-type terpenoids and flavonoids, showing that plants use both nonenzymatic and enzymatic intermolecular [4 + 2] cycloadditions to create complex molecules. Inspired by the biomimetic syntheses, we identify an intermolecular Diels-Alderase hidden in the biosynthetic pathway of mulberry Diels-Alder-type cycloadducts using a biosynthetic intermediate probe-based target identification strategy. This enzyme, MaDA, is an endo-selective Diels-Alderase and is then functionally characterized as a standalone intermolecular Diels-Alderase with a concerted but asynchronous mechanism. We also discover the exo-selective intermolecular Diels-Alderases in Morus plants. Both the endo- and exo-selective Diels-Alderases feature a broad substrate scope, but their mechanisms for controlling the endo/exo pathway are different. These unique intermolecular Diels-Alderases phylogenetically form a subgroup of FAD-dependent enzymes that can be found only in moraceous plants, explaining why this type of [4 + 2] cycloadduct is unique to moraceous plants. Further studies of the evolutionary mechanism reveal that an FAD-dependent oxidocyclase could acquire the Diels-Alderase activity via four critical amino acid mutations and then gradually lose its original oxidative activity to become a standalone Diels-Alderase during the natural evolution. Based on these insights, we designed new Diels-Alderases and achieved the diversity-oriented chemoenzymatic synthesis of D-A products using either naturally occurring or engineered Diels-Alderases.Overall, this Account describes our decade-long efforts to discover the intermolecular Diels-Alderases in Morus plants, particularly highlighting the importance of biomimetic synthesis and chemical proteomics in discovering new intermolecular Diels-Alderases from plants. Meanwhile, this Account also covers the evolutionary and catalytic mechanism study of intermolecular Diels-Alderases that may provide new insights into how to discover and design new Diels-Alderases as powerful biocatalysts for organic synthesis.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Plant antibacterials: The challenges and opportunities

Nature possesses an inexhaustible reservoir of agents that could serve as alternatives to combat the growing threat of antimicrobial resistance (AMR). While some of the most effective drugs for treating bacterial infections originate from natural sources, they have predominantly been derived from fungal and bacterial species. However, a substantial body of literature is available on the promising antibacterial properties of plant-derived compounds. In this comprehensive review, we address the major challenges associated with the discovery and development of plant-derived antimicrobial compounds, which have acted as obstacles preventing their clinical use. These challenges encompass limited sourcing, the risk of agent rediscovery, suboptimal drug metabolism, and pharmacokinetics (DMPK) properties, as well as a lack of knowledge regarding molecular targets and mechanisms of action, among other pertinent issues. Our review underscores the significance of these challenges and their implications in the quest for the discovery and development of effective plant-derived antimicrobial agents. Through a critical examination of the current state of research, we give valuable insights that will advance our understanding of these classes of compounds, offering potential solutions to the global crisis of AMR. © 2017 Elsevier Inc. All rights reserved.

Uptake, translocation, and metabolism of organophosphate esters (OPEs) in plants and health perspective for human: A review

Plant uptake, accumulation, and transformation of organophosphate esters (OPEs) play vital roles in their geochemical cycles and exposure risks. Here we reviewed the recent research advances in OPEs in plants. The mean OPE concentrations based on dry/wet/lipid weight varied in 4.80-3,620/0.287-26.8/12,000-315,000 ng g-1 in field plants, and generally showed positive correlations with those in plant habitats. OPEs with short-chain substituents and high hydrophilicity, particularly the commonly used chlorinated OPEs, showed dominance in most plant samples, whereas some tree barks, fruits, seeds, and roots demonstrated dominance of hydrophobic OPEs. Both hydrophilic and hydrophobic OPEs can enter plants via root and foliar uptake, and the former pathway is mainly passively mediated by various membrane proteins. After entry, different OPEs undergo diverse subcellular distributions and acropetal/basipetal/intergenerational translocations, depending on their physicochemical properties. Hydrophilic OPEs mainly exist in cell sap and show strong transferability, hydrophobic OPEs demonstrate dominant distributions in cell wall and limited migrations owing to the interception of Casparian strips and cell wall. Additionally, plant species, transpiration capacity, growth stages, commensal microorganisms, and habitats also affect OPE uptake and transfer in plants. OPE metabolites derived from various Phase I transformations and Phase II conjugations are increasingly identified in plants, and hydrolysis and hydroxylation are the most common metabolic processes. The metabolisms and products of OPEs are closely associated with their structures and degradation resistance and plant species. In contrast, plant-derived food consumption contributes considerably to the total dietary intakes of OPEs by human, particularly the cereals, and merits specifical attention. Based on the current research limitations, we proposed the research perspectives regarding OPEs in plants, with the emphases on their behavior and fate in field plants, interactions with plant-related microorganisms, multiple uptake pathways and mechanisms, and comprehensive screening analysis and risk evaluation.

Machine learning assists prediction of genes responsible for plant specialized metabolite biosynthesis by integrating multi-omics data

BACKGROUND

Plant specialized (or secondary) metabolites (PSM), also known as phytochemicals, natural products, or plant constituents, play essential roles in interactions between plants and environment. Although many research efforts have focused on discovering novel metabolites and their biosynthetic genes, the resolution of metabolic pathways and identified biosynthetic genes was limited by rudimentary analysis approaches and enormous number of candidate genes.

RESULTS

Here we integrated state-of-the-art automated machine learning (ML) frame AutoGluon-Tabular and multi-omics data from Arabidopsis to predict genes encoding enzymes involved in biosynthesis of plant specialized metabolite (PSM), focusing on the three main PSM categories: terpenoids, alkaloids, and phenolics. We found that the related features of genomics and proteomics were the top two crucial categories of features contributing to the model performance. Using only these key features, we built a new model in Arabidopsis, which performed better than models built with more features including those related with transcriptomics and epigenomics. Finally, the built models were validated in maize and tomato, and models tested for maize and trained with data from two other species exhibited either equivalent or superior performance to intraspecies predictions.

CONCLUSIONS

Our external validation results in grape and poppy on the one hand implied the applicability of our model to the other species, and on the other hand showed enormous potential to improve the prediction of enzymes synthesizing PSM with the inclusion of valid data from a wider range of species.

Tomato root specialized metabolites evolved through gene duplication and regulatory divergence within a biosynthetic gene cluster

Tremendous plant metabolic diversity arises from phylogenetically restricted specialized metabolic pathways. Specialized metabolites are synthesized in dedicated cells or tissues, with pathway genes sometimes colocalizing in biosynthetic gene clusters (BGCs). However, the mechanisms by which spatial expression patterns arise and the role of BGCs in pathway evolution remain underappreciated. In this study, we investigated the mechanisms driving acylsugar evolution in the Solanaceae. Previously thought to be restricted to glandular trichomes, acylsugars were recently found in cultivated tomato roots. We demonstrated that acylsugars in cultivated tomato roots and trichomes have different sugar cores, identified root-enriched paralogs of trichome acylsugar pathway genes, and characterized a key paralog required for root acylsugar biosynthesis, SlASAT1-LIKE (SlASAT1-L), which is nested within a previously reported trichome acylsugar BGC. Last, we provided evidence that ASAT1-L arose through duplication of its paralog, ASAT1, and was trichome-expressed before acquiring root-specific expression in the Solanum genus. Our results illuminate the genomic context and molecular mechanisms underpinning metabolic diversity in plants.

Genetic control of abiotic stress-related specialized metabolites in sunflower

BACKGROUND

Abiotic stresses in plants include all the environmental conditions that significantly reduce yields, like drought and heat. One of the most significant effects they exert at the cellular level is the accumulation of reactive oxygen species, which cause extensive damage. Plants possess two mechanisms to counter these molecules, i.e. detoxifying enzymes and non-enzymatic antioxidants, which include many classes of specialized metabolites. Sunflower, the fourth global oilseed, is considered moderately drought resistant. Abiotic stress tolerance in this crop has been studied using many approaches, but the control of specialized metabolites in this context remains poorly understood. Here, we performed the first genome-wide association study using abiotic stress-related specialized metabolites as molecular phenotypes in sunflower. After analyzing leaf specialized metabolites of 450 hybrids using liquid chromatography-mass spectrometry, we selected a subset of these compounds based on their association with previously known abiotic stress-related quantitative trait loci. Eventually, we characterized these molecules and their associated genes.

RESULTS

We putatively annotated 30 compounds which co-localized with abiotic stress-related quantitative trait loci and which were associated to seven most likely candidate genes. A large proportion of these compounds were potential antioxidants, which was in agreement with the role of specialized metabolites in abiotic stresses. The seven associated most likely candidate genes, instead, mainly belonged to cytochromes P450 and glycosyltransferases, two large superfamilies which catalyze greatly diverse reactions and create a wide variety of chemical modifications. This was consistent with the high plasticity of specialized metabolism in plants.

CONCLUSIONS

This is the first characterization of the genetic control of abiotic stress-related specialized metabolites in sunflower. By providing hints concerning the importance of antioxidant molecules in this biological context, and by highlighting some of the potential molecular mechanisms underlying their biosynthesis, it could pave the way for novel applications in breeding. Although further analyses will be required to better understand this topic, studying how antioxidants contribute to the tolerance to abiotic stresses in sunflower appears as a promising area of research.

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.

Differential gene expression and potential regulatory network of fatty acid biosynthesis during fruit and leaf development in yellowhorn (Xanthoceras sorbifolium), an oil-producing tree with significant deployment values

Xanthoceras sorbifolium (yellowhorn) is a woody oil plant with super stress resistance and excellent oil characteristics. The yellowhorn oil can be used as biofuel and edible oil with high nutritional and medicinal value. However, genetic studies on yellowhorn are just in the beginning, and fundamental biological questions regarding its very long-chain fatty acid (VLCFA) biosynthesis pathway remain largely unknown. In this study, we reconstructed the VLCFA biosynthesis pathway and annotated 137 genes encoding relevant enzymes. We identified four oleosin genes that package triacylglycerols (TAGs) and are specifically expressed in fruits, likely playing key roles in yellowhorn oil production. Especially, by examining time-ordered gene co-expression network (TO-GCN) constructed from fruit and leaf developments, we identified key enzymatic genes and potential regulatory transcription factors involved in VLCFA synthesis. In fruits, we further inferred a hierarchical regulatory network with MYB-related (XS03G0296800) and B3 (XS02G0057600) transcription factors as top-tier regulators, providing clues into factors controlling carbon flux into fatty acids. Our results offer new insights into key genes and transcriptional regulators governing fatty acid production in yellowhorn, laying the foundation for efforts to optimize oil content and fatty acid composition. Moreover, the gene expression patterns and putative regulatory relationships identified here will inform metabolic engineering and molecular breeding approaches tailored to meet biofuel and bioproduct demands.

The evolutionary trajectories of specialized metabolites towards antiviral defense system in plants

Viral infections in plants pose major challenges to agriculture and global food security in the twenty-first century. Plants have evolved a diverse range of specialized metabolites (PSMs) for defenses against pathogens. Although, PSMs-mediated plant-microorganism interactions have been widely discovered, these are mainly confined to plant-bacteria or plant-fungal interactions. PSM-mediated plant-virus interaction, however, is more complicated often due to the additional involvement of virus spreading vectors. Here, we review the major classes of PSMs and their emerging roles involved in antiviral resistances. In addition, evolutionary scenarios for PSM-mediated interactions between plant, virus and virus-transmitting vectors are presented. These advancements in comprehending the biochemical language of PSMs during plant-virus interactions not only lay the foundation for understanding potential co-evolution across life kingdoms, but also open a gateway to the fundamental principles of biological control strategies and beyond.

Selection and Verification of Standardized Reference Genes of Angelica dahurica under Various Abiotic Stresses by Real-Time Quantitative PCR

In traditional Chinese medicine, Angelica dahurica is a valuable herb with numerous therapeutic applications for a range of ailments. There have not yet been any articles on the methodical assessment and choice of the best reference genes for A. dahurica gene expression studies. Real-time quantitative PCR (RT-qPCR) is widely employed as the predominant method for investigating gene expression. In order to ensure the precise determination of target gene expression outcomes in RT-qPCR analysis, it is imperative to employ stable reference genes. In this study, a total of 11 candidate reference genes including SAND family protein (SAND), polypyrimidine tract-binding protein (PTBP), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), actin (ACT), TIP41-like protein (TIP41), cyclophilin 2 (CYP2), elongation factor 1 α (EF1α), ubiquitin-protein ligase 9 (UBC9), tubulin β-6 (TUB6), thioredoxin-like protein YLS8 (YLS8), and tubulin-α (TUBA) were selected from the transcriptome of A. dahurica. Subsequently, three statistical algorithms (geNorm, NormFinder, and BestKeeper) were employed to assess the stability of their expression patterns across seven distinct stimulus treatments. The outcomes obtained from these analyses were subsequently amalgamated into a comprehensive ranking using RefFinder. Additionally, one target gene, phenylalanine ammonia-lyase (PAL), was used to confirm the effectiveness of the selected reference genes. According to the findings of this study, the two most stable reference genes for normalizing the expression of genes in A. dahurica are TIP41 and UBC9. Overall, our research has determined the appropriate reference genes for RT-qPCR in A. dahurica and provides a crucial foundation for gene screening and identifying genes associated with the biosynthesis of active ingredients in A. dahurica.

Single-cell RNA sequencing reveals a hierarchical transcriptional regulatory network of terpenoid biosynthesis in cotton secretory glandular cells

Plants can synthesize a wide range of terpenoids in response to various environmental cues. However, the specific regulatory mechanisms governing terpenoid biosynthesis at the cellular level remain largely elusive. In this study, we employed single-cell RNA sequencing to comprehensively characterize the transcriptome profile of cotton leaves and established a hierarchical transcriptional network regulating cell-specific terpenoid production. We observed substantial expression levels of genes associated with the biosynthesis of both volatile terpenes (such as β-caryophyllene and β-myrcene) and non-volatile gossypol-type terpenoids in secretory glandular cells. Moreover, two novel transcription factors, namely GoHSFA4a and GoNAC42, are identified to function downstream of the Gossypium PIGMENT GLAND FORMATION genes. Both transcription factors could directly regulate the expression of terpenoid biosynthetic genes in secretory glandular cells in response to developmental and environmental stimuli. For convenient retrieval of the single-cell RNA sequencing data generated in this study, we developed a user-friendly web server . Our findings not only offer valuable insights into the precise regulation of terpenoid biosynthesis genes in cotton leaves but also provide potential targets for cotton breeding endeavors.

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.

Chromosome-level genome assembly of Niphotrichum japonicum provides new insights into heat stress responses in mosses

With a diversity of approximately 22,000 species, bryophytes (hornworts, liverworts, and mosses) represent a major and diverse lineage of land plants. Bryophytes can thrive in many extreme environments as they can endure the stresses of drought, heat, and cold. The moss Niphotrichum japonicum (Grimmiaceae, Grimmiales) can subsist for extended periods under heat and drought conditions, providing a good candidate for studying the genetic basis underlying such high resilience. Here, we de novo assembled the genome of N. japonicum using Nanopore long reads combined with Hi-C scaffolding technology to anchor the 191.61 Mb assembly into 14 pseudochromosomes. The genome structure of N. japonicum's autosomes is mostly conserved and highly syntenic, in contrast to the sparse and disordered genes present in its sex chromosome. Comparative genomic analysis revealed the presence of 10,019 genes exclusively in N. japonicum. These genes may contribute to the species-specific resilience, as demonstrated by the gene ontology (GO) enrichment. Transcriptome analysis showed that 37.44% (including 3,107 unique genes) of the total annotated genes (26,898) exhibited differential expression as a result of heat-induced stress, and the mechanisms that respond to heat stress are generally conserved across plants. These include the upregulation of HSPs, LEAs, and reactive oxygen species (ROS) scavenging genes, and the downregulation of PPR genes. N. japonicum also appears to have distinctive thermal mechanisms, including species-specific expansion and upregulation of the Self-incomp_S1 gene family, functional divergence of duplicated genes, structural clusters of upregulated genes, and expression piggybacking of hub genes. Overall, our study highlights both shared and species-specific heat tolerance strategies in N. japonicum, providing valuable insights into the heat tolerance mechanism and the evolution of resilient plants.

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.

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.

Tetraose steroidal glycoalkaloids from potato provide resistance against Alternaria solani and Colorado potato beetle

Plants with innate disease and pest resistance can contribute to more sustainable agriculture. Natural defence compounds produced by plants have the potential to provide a general protective effect against pathogens and pests, but they are not a primary target in resistance breeding. Here, we identified a wild relative of potato, Solanum commersonii, that provides us with unique insight in the role of glycoalkaloids in plant immunity. We cloned two atypical resistance genes that provide resistance to Alternaria solani and Colorado potato beetle through the production of tetraose steroidal glycoalkaloids (SGA). Moreover, we provide in vitro evidence to show that these compounds have potential against a range of different (potato pathogenic) fungi. This research links structural variation in SGAs to resistance against potato diseases and pests. Further research on the biosynthesis of plant defence compounds in different tissues, their toxicity, and the mechanisms for detoxification, can aid the effective use of such compounds to improve sustainability of our food production.

Genome assembly of Polygala tenuifolia provides insights into its karyotype evolution and triterpenoid saponin biosynthesis

Polygala tenuifolia is a perennial medicinal plant that has been widely used in traditional Chinese medicine for treating mental diseases. However, the lack of genomic resources limits the insight into its evolutionary and biological characterization. In the present work, we reported the P. tenuifolia genome, the first genome assembly of the Polygalaceae family. We sequenced and assembled this genome by a combination of Illumnina, PacBio HiFi, and Hi-C mapping. The assembly includes 19 pseudochromosomes covering ~92.68% of the assembled genome (~769.62 Mb). There are 36 463 protein-coding genes annotated in this genome. Detailed comparative genome analysis revealed that P. tenuifolia experienced two rounds of whole genome duplication that occurred ~39-44 and ~18-20 million years ago, respectively. Accordingly, we systematically reconstructed ancestral chromosomes of P. tenuifolia and inferred its chromosome evolution trajectories from the common ancestor of core eudicots to the present species. Based on the transcriptomics data, enzyme genes and transcription factors involved in the synthesis of triterpenoid saponin in P. tenuifolia were identified. Further analysis demonstrated that whole-genome duplications and tandem duplications play critical roles in the expansion of P450 and UGT gene families, which contributed to the synthesis of triterpenoid saponins. The genome and transcriptome data will not only provide valuable resources for comparative and functional genomic researches on Polygalaceae, but also shed light on the synthesis of triterpenoid saponin.

A redundant transcription factor network steers spatiotemporal Arabidopsis triterpene synthesis

Plant specialized metabolites modulate developmental and ecological functions and comprise many therapeutic and other high-value compounds. However, the mechanisms determining their cell-specific expression remain unknown. Here we describe the transcriptional regulatory network that underlies cell-specific biosynthesis of triterpenes in Arabidopsis thaliana root tips. Expression of thalianol and marneral biosynthesis pathway genes depends on the phytohormone jasmonate and is limited to outer tissues. We show that this is promoted by the activity of redundant bHLH-type transcription factors from two distinct clades and coactivated by homeodomain factors. Conversely, the DOF-type transcription factor DAG1 and other regulators prevent expression of the triterpene pathway genes in inner tissues. We thus show how precise expression of triterpene biosynthesis genes is determined by a robust network of transactivators, coactivators and counteracting repressors.

Lavandula Species, Their Bioactive Phytochemicals, and Their Biosynthetic Regulation

Lavandula species are one of the most useful aromatic and medicinal plants and have great economic potential. The phytopharmaceutical contribution of the secondary metabolites of the species is unquestionable. Most recent studies have been focusing on the elucidation of the genetic background of secondary metabolite production in lavender species. Therefore, knowledge of not only genetic but especially epigenetic mechanisms for the regulation of secondary metabolites is necessary for the modification of those biosynthesis processes and the understanding of genotypic differences in the content and compositional variability of these products. The review discusses the genetic diversity of Lavandula species in relation to the geographic area, occurrence, and morphogenetic factors. The role of microRNAs in secondary-metabolites biosynthesis is described.

Single-cell multi-omics in the medicinal plant Catharanthus roseus

Advances in omics technologies now permit the generation of highly contiguous genome assemblies, detection of transcripts and metabolites at the level of single cells and high-resolution determination of gene regulatory features. Here, using a complementary, multi-omics approach, we interrogated the monoterpene indole alkaloid (MIA) biosynthetic pathway in Catharanthus roseus, a source of leading anticancer drugs. We identified clusters of genes involved in MIA biosynthesis on the eight C. roseus chromosomes and extensive gene duplication of MIA pathway genes. Clustering was not limited to the linear genome, and through chromatin interaction data, MIA pathway genes were present within the same topologically associated domain, permitting the identification of a secologanin transporter. Single-cell RNA-sequencing revealed sequential cell-type-specific partitioning of the leaf MIA biosynthetic pathway that, when coupled with a single-cell metabolomics approach, permitted the identification of a reductase that yields the bis-indole alkaloid anhydrovinblastine. We also revealed cell-type-specific expression in the root MIA pathway.

Vacuolar MATE/DTX protein-mediated cucurbitacin C transport is co-regulated with bitterness biosynthesis in cucumber

Membrane-localized transporters constitute important components for specialized metabolism in plants. However, due to the vast array of specialized metabolites produced by plants, and the large families of transporter genes, knowledge about the intracellular and intercellular transport of plant metabolites is still in its infancy. Cucurbitacins are bitter and defensive triterpenoids produced mainly in the cucurbits. Using a comparative genomics and multi-omics approach, a MATE gene (CsMATE1), physically clustered with cucurbitacin C (CuC) biosynthetic genes, was identified and functionally shown to sequester CuC in cucumber leaf mesophyll cells. Notably, the CuC transport process is strictly co-regulated with CuC biosynthesis. CsMATE1 clustering with bitterness biosynthesis genes may provide benefits and a basis for this feedback regulation on CuC sequestration and biosynthesis. Identification of transport systems for plant-specialized metabolites can accelerate the metabolic engineering of high-value-added compounds by simplifying their purification process.

WAX INDUCER 1 Regulates β-Diketone Biosynthesis by Mediating Expression of the Cer-cqu Gene Cluster in Barley

Plant surface properties are crucial determinants of resilience to abiotic and biotic stresses. The outer layer of the plant cuticle consists of chemically diverse epicuticular waxes. The WAX INDUCER1/SHINE subfamily of APETALA2/ETHYLENE RESPONSIVE FACTORS regulates cuticle properties in plants. In this study, four barley genes homologous to the Arabidopsis thaliana AtWIN1 gene were mutated using RNA-guided Cas9 endonuclease. Mutations in one of them, the HvWIN1 gene, caused a recessive glossy sheath phenotype associated with β-diketone deficiency. A complementation test for win1 knockout (KO) and cer-x mutants showed that Cer-X and WIN1 are allelic variants of the same genomic locus. A comparison of the transcriptome from leaf sheaths of win1 KO and wild-type plants revealed a specific and strong downregulation of a large gene cluster residing at the previously known Cer-cqu locus. Our findings allowed us to postulate that the WIN1 transcription factor in barley is a master mediator of the β-diketone biosynthesis pathway acting through developmental stage- and organ-specific transactivation of the Cer-cqu gene cluster.

A reference-grade genome assembly for Astragalus mongholicus and insights into the biosynthesis and high accumulation of triterpenoids and flavonoids in its roots

Astragalus membranaceus var. mongholicus (AMM), a member of the Leguminosae, is one of the most important medicinal plants worldwide. The dried roots of AMM have a wide range of pharmacological effects and are a traditional Chinese medicine. Here, we report the first chromosome-level reference genome of AMM, comprising nine pseudochromosomes with a total size of 1.47 Gb and 27 868 protein-encoding genes. Comparative genomic analysis reveals that AMM has not experienced an independent whole-genome duplication (WGD) event after the WGD event shared by the Papilionoideae species. Analysis of long terminal repeat retrotransposons suggests a recent burst of these elements at approximately 0.13 million years ago, which may explain the large size of the AMM genome. Multiple gene families involved in the biosynthesis of triterpenoids and flavonoids were expanded, and our data indicate that tandem duplication has been the main driver for expansion of these families. Among the expanded families, the phenylalanine ammonia-lyase gene family was primarily expressed in the roots of AMM, suggesting their roles in the biosynthesis of phenylpropanoid compounds. The functional versatility of 2,3-oxidosqualene cyclase genes in cluster III may play a critical role in the diversification of triterpenoids in AMM. Our findings provide novel insights into triterpenoid and flavonoid biosynthesis and can facilitate future research on the genetics and medical applications of AMM.

A chromosome-level genome assembly reveals that a bipartite gene cluster formed via an inverted duplication controls monoterpenoid biosynthesis in Schizonepeta tenuifolia

Biosynthetic gene clusters (BGCs) are regions of a genome where genes involved in a biosynthetic pathway are in proximity. The origin and evolution of plant BGCs as well as their role in specialized metabolism remain largely unclear. In this study, we have assembled a chromosome-scale genome of Japanese catnip (Schizonepeta tenuifolia) and discovered a BGC that contains multiple copies of genes involved in four adjacent steps in the biosynthesis of p-menthane monoterpenoids. This BGC has an unprecedented bipartite structure, with mirrored biosynthetic regions separated by 260 kilobases. This bipartite BGC includes identical copies of a gene encoding an old yellow enzyme, a type of flavin-dependent reductase. In vitro assays and virus-induced gene silencing revealed that this gene encodes the missing isopiperitenone reductase. This enzyme evolved from a completely different enzyme family to isopiperitenone reductase from closely related Mentha spp., indicating convergent evolution of this pathway step. Phylogenomic analysis revealed that this bipartite BGC has emerged uniquely in the S. tenuifolia lineage and through insertion of pathway genes into a region rich in monoterpene synthases. The cluster gained its bipartite structure via an inverted duplication. The discovered bipartite BGC for p-menthane biosynthesis in S. tenuifolia has similarities to the recently described duplicated p-menthane biosynthesis gene pairs in the Mentha longifolia genome, providing an example of the convergent evolution of gene order. This work expands our understanding of plant BGCs with respect to both form and evolution, and highlights the power of BGCs for gene discovery in plant biosynthetic pathways.

Chromosome-level Genome Assembly of Euphorbia peplus, a Model System for Plant Latex, Reveals that Relative Lack of Ty3 Transposons Contributed to Its Small Genome Size

Euphorbia peplus (petty spurge) is a small, fast-growing plant that is native to Eurasia and has become a naturalized weed in North America and Australia. Euphorbia peplus is not only medicinally valuable, serving as a source for the skin cancer drug ingenol mebutate, but also has great potential as a model for latex production owing to its small size, ease of manipulation in the laboratory, and rapid reproductive cycle. To help establish E. peplus as a new model, we generated a 267.2-Mb Hi-C-anchored PacBio HiFi nuclear genome assembly with a BUSCO score of 98.5%, a genome annotation based on RNA-seq data from six organs, and publicly accessible tools including a genome browser and an interactive organ-specific expression atlas. Chromosome number is highly variable across Euphorbia species. Using a comparative analysis of our newly sequenced E. peplus genome with other Euphorbiaceae genomes, we show that variation in Euphorbia chromosome number between E. peplus and Euphorbia lathyris is likely due to fragmentation and rearrangement rather than chromosomal duplication followed by diploidization of the duplicated sequence. Moreover, we found that the E. peplus genome is relatively compact compared with related members of the genus in part due to restricted expansion of the Ty3 transposon family. Finally, we identify a large gene cluster that contains many previously identified enzymes in the putative ingenol mebutate biosynthesis pathway, along with additional gene candidates for this biosynthetic pathway. The genomic resources we have created for E. peplus will help advance research on latex production and ingenol mebutate biosynthesis in the commercially important Euphorbiaceae family.

Oryzalexin S biosynthesis: a cross-stitched disappearing pathway

Rice produces many diterpenoid phytoalexins and, reflecting the importance of these natural products in this important cereal crop plant, its genome contains three biosynthetic gene clusters (BGCs) for such metabolism. The chromosome 4 BGC (c4BGC) is largely associated with momilactone production, in part due to the presence of the initiating syn-copalyl diphosphate (CPP) synthase gene (OsCPS4). Oryzalexin S is also derived from syn-CPP. However, the relevant subsequently acting syn-stemarene synthase gene (OsKSL8) is not located in the c4BGC. Production of oryzalexin S further requires hydroxylation at carbons 2 and 19 (C2 and C19), presumably catalyzed by cytochrome P450 (CYP) monooxygenases. Here it is reported the closely related CYP99A2 and CYP99A3, whose genes are also found in the c4BGC catalyze the necessary C19-hydroxylation, while the closely related CYP71Z21 and CYP71Z22, whose genes are found in the recently reported chromosome 7 BGC (c7BGC), catalyze subsequent hydroxylation at C2α. Thus, oryzalexin S biosynthesis utilizes two distinct BGCs, in a pathway cross-stitched together by OsKSL8. Notably, in contrast to the widely conserved c4BGC, the c7BGC is subspecies (ssp.) specific, being prevalent in ssp. japonica and only rarely found in the other major ssp. indica. Moreover, while the closely related syn-stemodene synthase OsKSL11 was originally considered to be distinct from OsKSL8, it has now been reported to be a ssp. indica derived allele at the same genetic loci. Intriguingly, more detailed analysis indicates that OsKSL8(j) is being replaced by OsKSL11 (OsKSL8i), suggesting introgression from ssp. indica to (sub)tropical japonica, with concurrent disappearance of oryzalexin S production.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42994-022-00092-3.

Genetic and epigenetic control of the plant metabolome

Plant metabolites are mainly produced through chemical reactions catalysed by enzymes encoded in the genome. Mutations in enzyme-encoding or transcription factor-encoding genes can alter the metabolome by changing the enzyme's catalytic activity or abundance, respectively. Insertion of transposable elements into non-coding regions has also been reported to affect transcription and ultimately metabolite content. In addition to genetic mutations, transgenerational epigenetic variations have also been found to affect metabolic content by controlling the transcription of metabolism-related genes. However, the majority of cases reported so far, in which epigenetic mechanisms are associated with metabolism, are non-transgenerational, and are triggered by developmental signals or environmental stress. Although, accumulating research has provided evidence of strong genetic control of the metabolome, epigenetic control has been largely untouched. Here, we provide a review of the genetic and epigenetic control of metabolism with a focus on epigenetics. We discuss both transgenerational and non-transgenerational epigenetic marks regulating metabolism as well as prospects of the field of metabolic control where intricate interactions between genetics and epigenetics are involved.

Making small molecules in plants: A chassis for synthetic biology-based production of plant natural products

Plant natural products have been extensively exploited in food, medicine, flavor, cosmetic, renewable fuel, and other industrial sectors. Synthetic biology has recently emerged as a promising means for the cost-effective and sustainable production of natural products. Compared with engineering microbes for the production of plant natural products, the potential of plants as chassis for producing these compounds is underestimated, largely due to challenges encountered in engineering plants. Knowledge in plant engineering is instrumental for enabling the effective and efficient production of valuable phytochemicals in plants, and also paves the way for a more sustainable future agriculture. In this manuscript, we briefly recap the biosynthesis of plant natural products, focusing primarily on industrially important terpenoids, alkaloids, and phenylpropanoids. We further summarize the plant hosts and strategies that have been used to engineer the production of natural products. The challenges and opportunities of using plant synthetic biology to achieve rapid and scalable production of high-value plant natural products are also discussed.

Uncovering a miltiradiene biosynthetic gene cluster in the Lamiaceae reveals a dynamic evolutionary trajectory

The spatial organization of genes within plant genomes can drive evolution of specialized metabolic pathways. Terpenoids are important specialized metabolites in plants with diverse adaptive functions that enable environmental interactions. Here, we report the genome assemblies of Prunella vulgaris, Plectranthus barbatus, and Leonotis leonurus. We investigate the origin and subsequent evolution of a diterpenoid biosynthetic gene cluster (BGC) together with other seven species within the Lamiaceae (mint) family. Based on core genes found in the BGCs of all species examined across the Lamiaceae, we predict a simplified version of this cluster evolved in an early Lamiaceae ancestor. The current composition of the extant BGCs highlights the dynamic nature of its evolution. We elucidate the terpene backbones generated by the Callicarpa americana BGC enzymes, including miltiradiene and the terpene (+)-kaurene, and show oxidization activities of BGC cytochrome P450s. Our work reveals the fluid nature of BGC assembly and the importance of genome structure in contributing to the origin of metabolites.

Genome assembly of Musa beccarii shows extensive chromosomal rearrangements and genome expansion during evolution of Musaceae genomes

BACKGROUND

Musa beccarii (Musaceae) is a banana species native to Borneo, sometimes grown as an ornamental plant. The basic chromosome number of Musa species is x = 7, 10, or 11; however, M. beccarii has a basic chromosome number of x = 9 (2n = 2x = 18), which is the same basic chromosome number of species in the sister genera Ensete and Musella. Musa beccarii is in the section Callimusa, which is sister to the section Musa. We generated a high-quality chromosome-scale genome assembly of M. beccarii to better understand the evolution and diversity of genomes within the family Musaceae.

FINDINGS

The M. beccarii genome was assembled by long-read and Hi-C sequencing, and genes were annotated using both long Iso-seq and short RNA-seq reads. The size of M. beccarii was the largest among all known Musaceae assemblies (∼570 Mbp) due to the expansion of transposable elements and increased 45S ribosomal DNA sites. By synteny analysis, we detected extensive genome-wide chromosome fusions and fissions between M. beccarii and the other Musa and Ensete species, far beyond those expected from differences in chromosome number. Within Musaceae, M. beccarii showed a reduced number of terpenoid synthase genes, which are related to chemical defense, and enrichment in lipid metabolism genes linked to the physical defense of the cell wall. Furthermore, type III polyketide synthase was the most abundant biosynthetic gene cluster (BGC) in M. beccarii. BGCs were not conserved in Musaceae genomes.

CONCLUSIONS

The genome assembly of M. beccarii is the first chromosome-scale genome assembly in the Callimusa section in Musa, which provides an important genetic resource that aids our understanding of the evolution of Musaceae genomes and enhances our knowledge of the pangenome.

Chromosome-scale genome assembly of Glycyrrhiza uralensis revealed metabolic gene cluster centred specialized metabolites biosynthesis

A high-quality genome assembly is imperative to explore the evolutionary basis of characteristic attributes that define chemotype and provide essential resources for a molecular breeding strategy for enhanced production of medicinal metabolites. Here, using single-molecule high-fidelity (HiFi) sequencing reads, we report chromosome-scale genome assembly for Chinese licorice (Glycyrrhiza uralensis), a widely used herbal and natural medicine. The entire genome assembly was achieved in eight chromosomes, with contig and scaffold N50 as 36.02 and 60.2 Mb, respectively. With only 17 assembly gaps and half of the chromosomes having no or one assembly gap, the presented genome assembly is among the best plant genomes to date. Our results showed an advantage of using highly accurate long-read HiFi sequencing data for assembling a highly heterozygous genome including its complexed repeat content. Additionally, our analysis revealed that G. uralensis experienced a recent whole-genome duplication at approximately 59.02 million years ago post a gamma (γ) whole-genome triplication event, which contributed to its present chemotype features. The metabolic gene cluster analysis identified 355 gene clusters, which included the entire biosynthesis pathway of glycyrrhizin. The genome assembly and its annotations provide an essential resource for licorice improvement through molecular breeding and the discovery of valuable genes for engineering bioactive components and understanding the evolution of specialized metabolites biosynthesis.

The Vinca minor genome highlights conserved evolutionary traits in monoterpene indole alkaloid synthesis

Vinca minor, also known as the lesser periwinkle, is a well-known species from the Apocynaceae, native to central and southern Europe. This plant synthesizes monoterpene indole alkaloids, which are a class of specialized metabolites displaying a wide range of bioactive- and pharmacologically important properties. Within the almost 50 monoterpene indole alkaloids it produces, V. minor mainly accumulates vincamine, which is commercially used as a nootropic. Using a combination of Oxford Nanopore Technologies long read- and Illumina short-read sequencing, a 679,098 Mb V. minor genome was assembled into 296 scaffolds with an N50 scaffold length of 6 Mb, and encoding 29,624 genes. These genes were functionally annotated and used in a comparative genomic analysis to establish gene families and to investigate gene family expansion and contraction across the phylogenetic tree. Furthermore, homology-based monoterpene indole alkaloid gene predictions together with a metabolic analysis across 4 different V. minor tissue types guided the identification of candidate monoterpene indole alkaloid genes. These candidates were finally used to identify monoterpene indole alkaloid gene clusters, which combined with synteny analysis allowed for the discovery of a functionally validated vincadifformine-16-hydroxylase, reinforcing the potential of this dataset for monoterpene indole alkaloids gene discovery. It is expected that access to these resources will facilitate the elucidation of unknown monoterpene indole alkaloid biosynthetic routes with the potential of transferring these pathways to heterologous expression systems for large-scale monoterpene indole alkaloid production.

Comparative transcriptome analysis reveals key genes for polyphyllin difference in five Paris species

Paris species accumulate a large amount of steroidal saponins, which have numerous pharmacological activities and have become an essential component in many patented drugs. However, only two among all Paris species. Paris are identified as official sources due to high level of bioactive compounds. To clarify the composition of steroidal saponins and the molecular basis behind the differences between species, we investigated transcriptome and metabolic profiles of leaves and rhizomes in Paris polyphylla var. chinensis (PPC), Paris polyphylla var. yunnanensis (PPY), Paris polyphylla var. stenophylla (PPS), Paris fargesii (PF), and Paris mairei (PM). Phytochemical results displayed that the accumulation of steroidal saponins was tissue- and species-specific. PF and PPS contained more steroidal saponins in leaves than rhizomes, while PPY accumulated more steroidal saponins in rhizomes than leaves. PPC and PM contained similar amounts of steroidal saponins in leaves and rhizomes. Transcriptome analysis illustrated that most differentially expressed genes related to the biosynthesis of steroidal saponins were abundantly expressed in rhizomes than leaves. Meanwhile, more biosynthetic genes had significant correlations with steroidal saponins in rhizomes than in leaves. The result of CCA indicated that ACAT, DXS, DWF1, and CYP90 constrained 97.35% of the variance in bioactive compounds in leaves, whereas CYP72, UGT73, ACAT, and GPPS constrained 98.61% of the variance in phytochemicals in rhizomes. This study provided critical information for enhancing the production of steroidal saponins by biotechnological approaches and methodologies.

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.

Plant metabolic gene clusters in the multi-omics era

Secondary metabolism in plants gives rise to a vast array of small-molecule natural products. The discovery of operon-like gene clusters in plants has provided a new perspective on the evolution of specialized metabolism and the opportunity to rapidly advance the metabolic engineering of natural product production. Here, we review historical aspects of the study of plant metabolic gene clusters as well as general strategies for identifying plant metabolic gene clusters in the multi-omics era. We also emphasize the exploration of their natural variation and evolution, as well as new strategies for the prospecting of plant metabolic gene clusters and a deeper understanding of how their structure influences their function.

Chromosome-scale assembly of the Moringa oleifera Lam. genome uncovers polyploid history and evolution of secondary metabolism pathways through tandem duplication

The African Orphan Crops Consortium (AOCC) selected the highly nutritious, fast growing and drought tolerant tree crop moringa (Moringa oleifera Lam.) as one of the first of 101 plant species to have its genome sequenced and a first draft assembly was published in 2019. Given the extensive uses and culture of moringa, often referred to as the multipurpose tree, we generated a significantly improved new version of the genome based on long-read sequencing into 14 pseudochromosomes equivalent to n = 14 haploid chromosomes. We leveraged this nearly complete version of the moringa genome to investigate main drivers of gene family and genome evolution that may be at the origin of relevant biological innovations including agronomical favorable traits. Our results reveal that moringa has not undergone any additional whole-genome duplication (WGD) or polyploidy event beyond the gamma WGD shared by all core eudicots. Moringa duplicates retained following that ancient gamma events are also enriched for functions commonly considered as dosage balance sensitive. Furthermore, tandem duplications seem to have played a prominent role in the evolution of specific secondary metabolism pathways including those involved in the biosynthesis of bioactive glucosinolate, flavonoid, and alkaloid compounds as well as of defense response pathways and might, at least partially, explain the outstanding phenotypic plasticity attributed to this species. This study provides a genetic roadmap to guide future breeding programs in moringa, especially those aimed at improving secondary metabolism related traits.

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.

Biological Dark Matter Exploration using Data Mining for the Discovery of Antimicrobial Natural Products

The discovery of novel antimicrobials has significantly slowed down over the last three decades. At the same time, humans rely increasingly on antimicrobials because of the progressive antimicrobial resistance in medical practices, human communities, and the environment. Data mining is currently considered a promising option in the discovery of new antibiotics. Some of the advantages of data mining are the ability to predict chemical structures from sequence data, anticipation of the presence of novel metabolites, the understanding of gene evolution, and the corroboration of data from multiple omics technologies. This review analyzes the state-of-the-art for data mining in the fields of bacteria, fungi, and plant genomic data, as well as metabologenomics. It also summarizes some of the most recent research accomplishments in the field, all pinpointing to innovation through uncovering and implementing the next generation of antimicrobials.

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.

Chromosome-scale genome assembly of Rhododendron molle provides insights into its evolution and terpenoid biosynthesis

BACKGROUND

Rhododendron molle (Ericaceae) is a traditional Chinese medicine, which has been used to treat rheumatism and relieve pain since ancient times. The characteristic grayanoids of this plant have been demonstrated to be the chemical basis for the analgesic activity. Moreover, unlike morphine, these diterpenoids are non-addictive. Grayanoids mainly distribute in the leaves, flowers, roots, and fruits of R. molle, with low content. Currently the research on the biosynthesis of grayanoids is hindered, partially due to lack of the genomic information.

RESULTS

In the present study, a total of 744 Mb sequences were generated and assembled into 13 chromosomes. An ancient whole-genome duplication event (Ad-β) was discovered that occurred around 70 million years ago. Tandem and segmental gene duplications led to specific gene expansions in the terpene synthase and cytochrome P450 (CYP450) gene families. Two diterpene synthases were demonstrated to be responsible for the biosynthesis of 16α-hydroxy-ent-kaurane, the key precursor for grayanoids. Phylogenetic analysis revealed a species-specific bloom of the CYP71AU subfamily, which may involve the candidate CYP450s responsible for the biosynthesis of grayanoids. Additionally, three putative terpene biosynthetic gene clusters were found.

CONCLUSIONS

We reported the first genome assembly of R. molle and investigated the molecular basis underpinning terpenoids biosynthesis. Our work provides a foundation for elucidating the complete biosynthetic pathway of grayanoids and studying the terpenoids diversity in R. molle.

Convergent gene clusters underpin hyperforin biosynthesis in St John's wort

The meroterpenoid hyperforin is responsible for the antidepressant activity of St John's wort extracts, but the genes controlling its biosynthesis are unknown. Using genome mining and biochemical work, we characterize two biosynthetic gene clusters (BGCs) that encode the first three steps in the biosynthesis of hyperforin precursors. The findings of syntenic and phylogenetic analyses reveal the parallel assembly of the two BGCs. The syntenous BGC in Mesua ferrea indicates that the first cluster was assembled before the divergence of the Hypericaceae and Calophyllaceae families. The assembly of the second cluster is the result of a coalescence of genomic fragments after a major duplication event. The differences between the two BGCs - in terms of gene expression, response to methyl jasmonate, substrate specificity and subcellular localization of key enzymes - suggest that the presence of the two clusters could serve to generate separate pools of precursors. The parallel assembly of two BGCs with similar compositions in a single plant species is uncommon, and our work provides insights into how and when these gene clusters form. Our discovery helps to advance our understanding of the evolution of plant specialized metabolism and its genomic organization. Additionally, our results offer a foundation from which hyperforin biosynthesis can be more fully understood, and which can be used in future metabolic engineering applications.

Selection and validation of reference genes for quantitative real-time PCR normalization in Psoralea corylifolia (Babchi) under various abiotic stress

Psoralea corylifolia L. is a popular herb and has long been used in traditional Ayurvedic and Chinese medicine owing to its extensive pharmacological activities, especially in the treatment of various shin diseases. To date, the systematic evaluation and selection of the optimum reference genes for gene expression analysis of P. corylifolia were not reported. Reverse transcription-quantitative polymerase chain reaction (RT-qPCR) is a method for gene expression quantification. Selecting appropriate reference genes is the prerequisite for accurate normalization of RT-qPCR results. This study assessed the expression stability of 10 candidate reference genes under different abiotic stresses. First, amplification primers for RT-qPCR were designed and received testing and optimization. Then, expression data from each candidate gene were evaluated based on three statistical algorithms, and their results were further integrated into a comprehensive ranking based on the geometric mean. Additionally, one target gene, i.e., 1-aminocyclopropane-1-carboxylate oxidase (ACO), was used to validate the effectiveness of the selected reference. Our analysis suggested that thioredoxin-like protein YLS8 (YLS8), TIP41-like family protein (TIP41), and cyclophilin 2 (CYP2) genes provided superior expression normalization under different abiotic stresses. Overall, this work constitutes the first effort to select optimal endogenous controls for RT-qPCR studies of gene expression in P. corylifolia. It also provides a reasonable normalization standard and basis for further analysis of the gene expression of bioactive components in P. corylifolia.

Mapping‑by‑Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus

Rapeseed (Brassica napus L.) is an important oil crop and has the potential to serve as a highly productive source of protein. This protein exhibits an excellent amino acid composition and has high nutritional value for humans. Seed protein content (SPC) and seed oil content (SOC) are two complex quantitative and polygenic traits which are negatively correlated and assumed to be controlled by additive and epistatic effects. A reduction in seed glucosinolate (GSL) content is desired as GSLs cause a stringent and bitter taste. The goal here was the identification of genomic intervals relevant for seed GSL content and SPC/SOC. Mapping by sequencing (MBS) revealed 30 and 15 new and known genomic intervals associated with seed GSL content and SPC/SOC, respectively. Within these intervals, we identified known but also so far unknown putatively causal genes and sequence variants. A 4 bp insertion in the MYB28 homolog on C09 shows a significant association with a reduction in seed GSL content. This study provides insights into the genetic architecture and potential mechanisms underlying seed quality traits, which will enhance future breeding approaches in B. napus.

Omics-guided metabolic pathway discovery in plants: Resources, approaches, and opportunities

Plants produce a vast array of metabolites, the biosynthetic routes of which remain largely undetermined. Genome-scale enzyme and pathway annotations and omics technologies have revolutionized research to decrypt plant metabolism and produced a growing list of functionally characterized metabolic genes and pathways. However, what is known is still a tiny fraction of the metabolic capacity harbored by plants. Here, we review plant enzyme and pathway annotation resources and cutting-edge omics approaches to guide discovery and characterization of plant metabolic pathways. We also discuss strategies for improving enzyme function prediction by integrating protein 3D structure information and single cell omics. This review aims to serve as a primer for plant biologists to leverage omics datasets to facilitate understanding and engineering plant metabolism.

Gene discovery and virus-induced gene silencing reveal branched pathways to major classes of bioactive diterpenoids in Euphorbia peplus

SignificanceEuphorbia peplus, a member of the Euphorbia genus, is rich in jatrophane and ingenane diterpenoids. Using a metabolomics-guided transcriptomic approach to gene candidate identification, we have discovered a short-chain dehydrogenase gene involved in the production of the lathyrane jolkinol E. We have developed a virus-induced gene-silencing method in E. peplus that has allowed us to demonstrate the direct relationship between casbene and polycyclic diterpenoids and that jolkinol C acts as a key branch point intermediate in the production of ingenanes and jatrophanes. This work contributes both knowledge and tools for engineering production of bioactive diterpenoids in heterologous host systems, thus enabling their further evaluation and development.

Ancient plant-like terpene biosynthesis in corals

Octocorals are major contributors of terpenoid chemical diversity in the ocean. Natural products from other sessile marine animals are primarily biosynthesized by symbiotic microbes rather than by the host. Here, we challenge this long-standing paradigm by describing a monophyletic lineage of animal-encoded terpene cyclases (TCs) ubiquitous in octocorals. We characterized 15 TC enzymes from nine genera, several of which produce precursors of iconic coral-specific terpenoids, such as pseudopterosin, lophotoxin and eleutherobin. X-ray crystallography revealed that coral TCs share conserved active site residues and structural features with bacterial TCs. The identification of coral TCs enabled the targeted identification of the enzyme that constructs the coral-exclusive capnellane scaffold. Several TC genes are colocalized with genes that encode enzymes known to modify terpenes. This work presents an example of biosynthetic capacity in the kingdom Animalia that rivals the chemical complexity generated by plants, unlocking the biotechnological potential of octocorals for biomedical applications.

Horizontal transfer and evolution of the biosynthetic gene cluster for benzoxazinoids in plants

Benzoxazinoids are a class of protective and allelopathic plant secondary metabolites that have been identified in multiple grass species and are encoded by the Bx biosynthetic gene cluster (BGC) in maize. Data mining of 41 high-quality grass genomes identified complete Bx clusters (containing genes Bx1-Bx5 and Bx8) in three genera (Zea, Echinochloa, and Dichanthelium) of Panicoideae and partial clusters in Triticeae. The Bx cluster probably originated from gene duplication and chromosomal translocation of native homologs of Bx genes. An ancient Bx cluster that included additional Bx genes (e.g., Bx6) is presumed to have been present in ancestral Panicoideae. The ancient Bx cluster was putatively gained by the Triticeae ancestor via horizontal transfer (HT) from the ancestral Panicoideae and later separated into multiple segments on different chromosomes. Bx6 appears to have been under less constrained selection compared with the Bx cluster during the evolution of Panicoideae, as evidenced by the fact that it was translocated away from the Bx cluster in Zea mays, moved to other chromosomes in Echinochloa, and even lost in Dichanthelium. Further investigations indicate that purifying selection and polyploidization have shaped the evolutionary trajectory of Bx clusters in the grass family. This study provides the first candidate case of HT of a BGC between plants and sheds new light on the evolution of BGCs.

The lncRNA MARS modulates the epigenetic reprogramming of the marneral cluster in response to ABA

Clustered organization of biosynthetic non-homologous genes is emerging as a characteristic feature of plant genomes. The co-regulation of clustered genes seems to largely depend on epigenetic reprogramming and three-dimensional chromatin conformation. In this study, we identified the long non-coding RNA (lncRNA) MARneral Silencing (MARS), localized inside the Arabidopsis marneral cluster, which controls the local epigenetic activation of its surrounding region in response to abscisic acid (ABA). MARS modulates the POLYCOMB REPRESSIVE COMPLEX 1 (PRC1) component LIKE HETEROCHROMATIN PROTEIN 1 (LHP1) binding throughout the cluster in a dose-dependent manner, determining H3K27me3 deposition and chromatin condensation. In response to ABA, MARS decoys LHP1 away from the cluster and promotes the formation of a chromatin loop bringing together the MARNERAL SYNTHASE 1 (MRN1) locus and a distal ABA-responsive enhancer. The enrichment of co-regulated lncRNAs in clustered metabolic genes in Arabidopsis suggests that the acquisition of novel non-coding transcriptional units may constitute an additional regulatory layer driving the evolution of biosynthetic pathways.