Insights from a century of data reveal global trends in ex situ living plant collections

Ex situ living plant collections play a crucial role in providing nature-based solutions to twenty-first century global challenges. However, the complex dynamics of these artificial ecosystems are poorly quantified and understood, affecting biodiversity storage, conservation and utilization. To evaluate the management of ex situ plant diversity, we analysed a century of data comprising 2.2 million records, from a meta-collection currently holding ~500,000 accessions and 41% of global ex situ species diversity. Our study provides critical insights into the historical evolution, current state and future trajectory of global living collections. We reveal sigmoidal growth of a meta-collection that has reached capacity in both total accessions and total diversity, and identify intrinsic constraints on biodiversity management, including a median survival probability of 15 years. We explore the impact of external constraints and quantify the influence of the Convention on Biological Diversity, which we link to reduced acquisition of wild-origin and internationally sourced material by 44% and 38%, respectively. We further define the impact of these constraints on ex situ conservation but highlight targeted initiatives that successfully mitigate these challenges. Ultimately, our study underscores the urgent need for strategic prioritization and the re-evaluation of ex situ biodiversity management to achieve both scientific and conservation goals.

Characterization and functional analysis of type III polyketide synthases in Selaginella moellendorffii

MAIN CONCLUSION

The evolutionary conservation of type III polyketide synthases (PKS) in Selaginella has been elucidated, and the critical amino acid residues of the anther-specific chalcone synthase-like enzyme (SmASCL) have been identified. Selaginella species are the oldest known vascular plants and a valuable resource for the study of metabolic evolution in land plants. Polyketides, especially flavonoids and sporopollenin precursors, are essential prerequisites for plant land colonization. Although type III polyketide synthases (PKS) are widely studied in seed plants, the related enzymes in Selaginella remain poorly characterized. Here, eight type III PKSs were identified in the Selaginella moellendorffii genome and classified into three clusters. Two PKSs were selected for further research based on their phylogenetic relationships and protein sequence similarity. Functional studies revealed that they were chalcone synthase (SmCHS) and anther-specific CHS-like enzyme (SmASCL). These enzymes are involved in the biosynthesis of flavonoids and sporopollenin, respectively. Their sequence information and enzymatic activity are similar to the orthologs in other plants. Phylogenetic analysis revealed that the ASCL and CHS enzymes were separated into two clades from the Bryophyta. These results suggest that CHS and ASCL emerged in the first land plants and then remained conserved during plant evolution. To study the structural basis of the enzymatic function of SmASCL, a series of mutants were constructed. The number of condensation reactions catalyzed by the P210L/Y211D and I200V/G201T double mutants exceeds that of the wild-type enzyme. Our study provides insight into the characteristics and functions of type III PKSs in S. moellendorffii. It also offers clues for a deeper understanding of the relationship between active sites and the enzymatic function of ASCLs.

Carbene-catalyzed chirality-controlled site-selective acylation of saccharides

Acylation stands as a fundamental process in both biological pathways and synthetic chemical reactions, with acylated saccharides and their derivatives holding diverse applications ranging from bioactive agents to synthetic building blocks. A longstanding objective in organic synthesis has been the site-selective acylation of saccharides without extensive pre-protection of alcohol units. In this study, we demonstrate that by simply altering the chirality of N-heterocyclic carbene (NHC) organic catalysts, the site-selectivity of saccharide acylation reactions can be effectively modulated. Our investigation reveals that this intriguing selectivity shift stems from a combination of factors, including chirality match/mismatch and inter- / intramolecular hydrogen bonding between the NHC catalyst and saccharide substrates. These findings provide valuable insights into catalyst design and reaction engineering, highlighting potential applications in glycoside analysis, such as fluorescent labelling, α/β identification, orthogonal reactions, and selective late-stage modifications.

Multi-omics analysis reveals the role of UGT72 family genes in arbutin biosynthesis in Pyrus and evolution driven by whole genome duplication

The UGT72 gene family encodes proteins that glycosylate phenylpropanoids, and thus contribute to the synthesis of various phenolic substances. However, their functional role and evolutionary history in Pyrus spp. remains poorly understood. Here we explored the evolution, amplification, coding region structural variation, and functional divergence of the UGT72 gene family and its subfamilies. Further, we identified functional genes involved in arbutin synthesis and functionally validated the key genes. 15 UGT72 genes were identified in the complete genome sequence and classified into two subfamilies of Pyrus betulifolia. Significant expansion of the UGT72 gene family occurred after genome duplication in P. betulifolia. 53.33 % of all UGT72 family genes were found to have undergone expansion via WGD/segmental duplication. A noteworthy discovery was that the amplification of functional genes such as PbUGT72B1714 during polyploidization, combined with the loss of vital motifs and variations at important sites within these genes, significantly impacted the diversification of arbutin metabolism. These findings offer novel insights into how gene gains and losses caused by WGDs have contributed to metabolic diversification and evolutionary adaptation in Pyrus, as well as a groundwork for more detailed investigations into the mechanisms of arbutin metabolism.

Analysis of substrate specificity of cytochrome P450 monooxygenases involved in trichothecene toxin biosynthesis

Trichothecenes are a structurally diverse family of toxic secondary metabolites produced by certain species of multiple fungal genera. All trichothecene analogs share a core 12,13-epoxytrichothec-9-ene (EPT) structure but differ in presence, absence and types of substituents attached to various positions of EPT. Formation of some of the structural diversity begins early in the biosynthetic pathway such that some producing species have few trichothecene biosynthetic intermediates in common. Cytochrome P450 monooxygenases (P450s) play critical roles in formation of trichothecene structural diversity. Within some species, relaxed substrate specificities of P450s allow individual orthologs of the enzymes to modify multiple trichothecene biosynthetic intermediates. It is not clear, however, whether the relaxed specificity extends to biosynthetic intermediates that are not produced by the species in which the orthologs originate. To address this knowledge gap, we used a mutant complementation-heterologous expression analysis to assess whether orthologs of three trichothecene biosynthetic P450s (TRI11, TRI13 and TRI22) from Fusarium sporotrichioides, Trichoderma arundinaceum, and Paramyrothecium roridum can modify trichothecene biosynthetic intermediates that they do not encounter in the organism in which they originated. The results indicate that TRI13 and TRI22 could not modify the intermediates that they do not normally encounter, whereas TRI11 could modify an intermediate that it does not normally encounter. These findings indicate that substrate promiscuity varies among trichothecene biosynthetic P450s. One structural feature that likely impacts the ability of the P450s to use biosynthetic intermediates as substrates is the presence and absence of an oxygen atom attached to carbon atom 3 of EPT.

Understanding metabolic diversification in plants: branchpoints in the evolution of specialized metabolism

Plants are chemical engineers par excellence. Collectively they make a vast array of structurally diverse specialized metabolites. The raw materials for building new pathways (genes encoding biosynthetic enzymes) are commonly recruited directly or indirectly from primary metabolism. Little is known about how new metabolic pathways and networks evolve in plants, or what key nodes contribute to branches that lead to the biosynthesis of diverse chemicals. Here we review the molecular mechanisms underlying the generation of biosynthetic branchpoints. We also consider examples in which new metabolites are formed through the joining of precursor molecules arising from different biosynthetic routes, a scenario that greatly increases both the diversity and complexity of specialized metabolism. Given the emerging importance of metabolic gene clustering in helping to identify new enzymes and pathways, we further cover the significance of biosynthetic gene clusters in relation to metabolic networks and dedicated biosynthetic pathways. In conclusion, an improved understanding of the branchpoints between metabolic pathways will be key in order to be able to predict and illustrate the complex structure of metabolic networks and to better understand the plasticity of plant metabolism. This article is part of the theme issue 'The evolution of plant metabolism'.

Applications of ancestral sequence reconstruction for understanding the evolution of plant specialized metabolism

Studies of enzymes in modern-day plants have documented the diversity of metabolic activities retained by species today but only provide limited insight into how those properties evolved. Ancestral sequence reconstruction (ASR) is an approach that provides statistical estimates of ancient plant enzyme sequences which can then be resurrected to test hypotheses about the evolution of catalytic activities and pathway assembly. Here, I review the insights that have been obtained using ASR to study plant metabolism and highlight important methodological aspects. Overall, studies of resurrected plant enzymes show that (i) exaptation is widespread such that even low or undetectable levels of ancestral activity with a substrate can later become the apparent primary activity of descendant enzymes, (ii) intramolecular epistasis may or may not limit evolutionary paths towards catalytic or substrate preference switches, and (iii) ancient pathway flux often differs from modern-day metabolic networks. These and other insights gained from ASR would not have been possible using only modern-day sequences. Future ASR studies characterizing entire ancestral metabolic networks as well as those that link ancient structures with enzymatic properties should continue to provide novel insights into how the chemical diversity of plants evolved. This article is part of the theme issue 'The evolution of plant metabolism'.

Multiplexed perturbation of yew reveals cryptic proteins that enable a total biosynthesis of baccatin III and Taxol precursors

Plants make complex and potent therapeutic molecules, but difficulties in sourcing from natural producers or chemical synthesis can challenge their use in the clinic. A prominent example is the anti-cancer therapeutic paclitaxel (Taxol ® ). Identification of the full paclitaxel biosynthetic pathway would enable heterologous drug production, but it has eluded discovery despite a half century of intensive research. Within the search space of Taxus' large, enzyme-rich genome, we suspected the complex paclitaxel pathway would be difficult to resolve using conventional gene co-expression analysis and small sample sets. To improve the resolution of gene set identification, we developed a multiplexed perturbation strategy to transcriptionally profile cell states spanning tissues, cell types, developmental stages, and elicitation conditions. This approach revealed a set of paclitaxel biosynthetic genes that segregate into expression modules that suggest consecutive biosynthetic sub-pathways. These modules resolved seven new genes that, when combined with previously known enzymes, are sufficient for the de novo biosynthesis and isolation of baccatin III, an industrial precursor for Taxol, in Nicotiana benthamiana leaves at levels comparable to the natural abundance in Taxus needles. Included are taxane 1β-hydroxylase (T1βH), taxane 9α-hydroxylase (T9αH), taxane 7β- O -acyltransferase (T7ΑΤ), taxane 7β- O -deacetylase (T7dA), taxane 9α- O -deacetylase (T9dA), and taxane 9-oxidase (T9ox). Importantly, the T9αH we discovered is distinct and independently evolved from those recently reported, which failed to yield baccatin III with downstream enzymes. Unexpectedly, we also found a nuclear transport factor 2 (NTF2)-like protein (FoTO1) crucial for high yields of taxanes; this gene promotes the formation of the desired product during the first taxane oxidation step, resolving a longstanding bottleneck in paclitaxel pathway reconstitution. Together with a new β-phenylalanine-CoA-ligase, the eight genes discovered in this study enables the complete reconstitution of 3'- N -debenzoyl-2'-deoxy-paclitaxel with a 20-enzyme pathway in Nicotiana plants. More broadly, we establish a generalizable approach for pathway discovery that scales the power of co-expression studies to match the complexity of specialized metabolism, enabling discovery of gene sets responsible for high-value biological functions.

The HOMO-LUMO Gap as Discriminator of Biotic from Abiotic Chemistries

Low-molecular-mass organic chemicals are widely discussed as potential indicators of life in extraterrestrial habitats. However, demarcation lines between biotic chemicals and abiotic chemicals have been difficult to define. Here, we have analyzed the potential utility of the quantum chemical property, HOMO-LUMO gap (HLG), as a novel proxy variable of life, since a significant trend towards incrementally smaller HLGs has been described in the genetically encoded amino acids. The HLG is a zeroth-order predictor of chemical reactivity. Comparing a set of 134 abiotic organic molecules recovered from meteorites, with 570 microbial and plant secondary metabolites thought to be exclusively biotic, we found that the average HLG of biotic molecules was significantly narrower (-10.4 ± 0.9 eV versus -12.4 ± 1.6 eV), with an effect size of g = 1.87. Limitation to hydrophilic molecules (XlogP < 2) improved the separation of biotic from abiotic compounds (g = 2.52). The "hydrophilic reactivity" quadrant defined by |HLG| < 11.25 eV and XlogP < 2 was populated exclusively by 183 biotic compounds and 6 abiotic compounds, 5 of which were nucleobases. We conclude that hydrophilic molecules with small HLGs represent valuable indicators of biotic activity, and we discuss the evolutionary plausibility of this inference.

Decoding resilience: ecology, regulation, and evolution of biosynthetic gene clusters

Secondary metabolism is crucial for plant survival and can generate chemistry with nutritional, therapeutic, and industrial value. Biosynthetic genes of selected secondary metabolites cluster within localised chromosomal regions. The arrangement of these biosynthetic gene clusters (BGCs) challenges the long-held model of random gene order in eukaryotes, raising questions about their regulation, ecological significance, and evolution. In this review, we address these questions by exploring the contribution of BGCs to ecologically relevant plant-biotic interactions, while also evaluating the molecular-(epi)genetic mechanisms controlling their coordinated stress- and tissue-specific expression. Based on evidence that BGCs have distinct chromatin signatures and are enriched with transposable elements (TEs), we integrate emerging hypotheses into an updated evolutionary model emphasising how stress-induced epigenetic processes have shaped BGC formation.

A developmental gradient reveals biosynthetic pathways to eukaryotic toxins in monocot geophytes

Numerous eukaryotic toxins that accumulate in geophytic plants are valuable in the clinic, yet their biosynthetic pathways have remained elusive. A notable example is the >150 Amaryllidaceae alkaloids (AmAs), including galantamine, an FDA-approved treatment for Alzheimer's disease. We show that while AmAs accumulate to high levels in many daffodil tissues, biosynthesis is localized to nascent, growing tissue at the leaf base. A similar trend is found in the production of steroidal alkaloids (e.g., cyclopamine) in corn lily. This model of active biosynthesis enabled the elucidation of a complete set of biosynthetic genes that can be used to produce AmAs. Taken together, our work sheds light on the developmental and enzymatic logic of diverse alkaloid biosynthesis in daffodils. More broadly, it suggests a paradigm for biosynthesis regulation in monocot geophytes, where plants are protected from herbivory through active charging of newly formed cells with eukaryotic toxins that persist as above-ground tissue develops.

Tissue ontogeny and chemical composition influence bacterial biodiversity in the wood and shoot tip of Populus nigra

Plant-microbe interactions significantly influence plant growth dynamics and adaptability. This study explores the impact of metabolites on microbial biodiversity in shoot tips and wood of Populus nigra under greenhouse conditions, using high-throughput sequencing and metabolite profiling. Branches from P. nigra were harvested, rooted, and transplanted into pots for growth. After 3 months, tissue samples from shoot tips and wood were collected, and metabolites extracted and analysed using GC-MS and LC-MS. Genomic DNA was extracted and subjected to high-throughput sequencing for bacterial biodiversity profiling. Both datasets were analysed using bioinformatic and statistical pipelines. Metabolite profiling indicated that shoot tips had a higher relative abundance of primary and secondary metabolites, including sugars, fatty acids, organic acids, phenolic acid derivatives and salicinoids, while wood was enriched in flavonoids. Bacterial biodiversity also differed significantly between these tissues, with Clostridiales, Bacteroidales and Bacillales dominating in shoot tips, associated with rapid growth and anaerobic fermentation, while wood tissues were characterized by diazotrophs from Rhizobiales, Sphingomonadales and Frankiales. PCoA clustering confirmed tissue-specific microbial differences. Functional analysis revealed an enrichment of fundamental cellular processes in shoot tips, while wood exhibited pathways related to degradation and mortality. Metabolite profiling revealed significant variations in primary and secondary metabolites, highlighting their influence on microbial biodiversity across plant tissues. The dominance of specific bacterial orders and distinct functional pathways in each tissue suggests a tailored microbial response to the unique environments of shoot tips and wood.

The evolutionary advantage of artemisinin production by Artemisia annua

Artemisinin, a potent antimalarial compound, is predominantly derived from Artemisia annua. The uniqueness of artemisinin production in A. annua lies in its complex biochemical pathways and genetic composition, distinguishing it from other plant species, even within the Asteraceae family. In this review, we investigate the potential of A. annua for artemisinin production, drawing evidence from natural populations and mutants. Leveraging high-quality whole-genome sequence analyses, we offer insights into the evolution of artemisinin biosynthesis. We also highlight current understanding of the protective functions of artemisinin in A. annua in response to both biotic and abiotic stresses. In addition, we explore the mechanisms used by A. annua to mitigate the phytotoxicity generated by artemisinin catabolism.

Mapping of specialized metabolite terms onto a plant phylogeny using text mining and large language models

Plants produce a staggering array of chemicals that are the basis for organismal function and important human nutrients and medicines. However, it is poorly defined how these compounds evolved and are distributed across the plant kingdom, hindering a systematic view and understanding of plant chemical diversity. Recent advances in plant genome/transcriptome sequencing have provided a well-defined molecular phylogeny of plants, on which the presence of diverse natural products can be mapped to systematically determine their phylogenetic distribution. Here, we built a proof-of-concept workflow where previously reported diverse tyrosine-derived plant natural products were mapped onto the plant tree of life. Plant chemical-species associations were mined from literature, filtered, evaluated through manual inspection of over 2500 scientific articles, and mapped onto the plant phylogeny. The resulting "phylochemical" map confirmed several highly lineage-specific compound class distributions, such as betalain pigments and Amaryllidaceae alkaloids. The map also highlighted several lineages enriched in dopamine-derived compounds, including the orders Caryophyllales, Liliales, and Fabales. Additionally, the application of large language models, using our manually curated data as a ground truth set, showed that post-mining processing can largely be automated with a low false-positive rate, critical for generating a reliable phylochemical map. Although a high false-negative rate remains a challenge, our study demonstrates that combining text mining with language model-based processing can generate broader phylochemical maps, which will serve as a valuable community resource to uncover key evolutionary events that underlie plant chemical diversity and enable system-level views of nature's millions of years of chemical experimentation.

Analysis of the Mechanism of Wood Vinegar and Butyrolactone Promoting Rapeseed Growth and Improving Low-Temperature Stress Resistance Based on Transcriptome and Metabolomics

Rapeseed is an important oil crop in the world. Wood vinegar could increase the yield and abiotic resistance of rapeseed. However, little is known about the underlying mechanisms of wood vinegar or its valid chemical components on rapeseed. In the present study, wood vinegar and butyrolactone (γ-Butyrolactone, one of the main components of wood vinegar) were applied to rapeseed at the seedling stage, and the molecular mechanisms of wood vinegar that affect rapeseed were studied by combining transcriptome and metabolomic analyses. The results show that applying wood vinegar and butyrolactone increases the biomass of rapeseed by increasing the leaf area and the number of pods per plant, and enhances the tolerance of rapeseed under low temperature by reducing membrane lipid oxidation and improving the content of chlorophyll, proline, soluble sugar, and antioxidant enzymes. Compared to the control, 681 and 700 differentially expressed genes were in the transcriptional group treated with wood vinegar and butyrolactone, respectively, and 76 and 90 differentially expressed metabolites were in the metabolic group. The combination of transcriptome and metabolomic analyses revealed the key gene-metabolic networks related to various pathways. Our research shows that after wood vinegar and butyrolactone treatment, the amino acid biosynthesis pathway of rapeseed may be involved in mediating the increase in rapeseed biomass, the proline metabolism pathway of wood vinegar treatment may be involved in mediating rapeseed's resistance to low-temperature stress, and the sphingolipid metabolism pathway of butyrolactone treatment may be involved in mediating rapeseed's resistance to low-temperature stress. It is suggested that the use of wood vinegar or butyrolactone are new approaches to increasing rapeseed yield and low-temperature resistance.

Metabolites-induced co-evolutionary warfare between plants, viruses, and their associated vectors: So close yet so far away

Agriculture and global food security encounter significant challenges due to viral threats. In the following decades, several molecular studies have focused on discovering biosynthetic pathways of numerous defensive and signaling compounds, as key regulators of plant interactions, either with viruses or their associated vectors. Nevertheless, the complexities of specialized metabolites mediated plant-virus-vector tripartite viewpoint and the identification of their co-evolutionary crossroads toward antiviral defense system, remain elusive. The current study reviews the various roles of plant-specialized metabolites (PSMs) and how plants use these metabolites to defend against viruses. It discusses recent examples of specialized metabolites that have broad-spectrum antiviral properties. Additionally, the study presents the co-evolutionary basis of metabolite-mediated plant-virus-insect interactions as a potential bioinspired approach to combat viral threats. The prospects also show promising metabolic engineering strategies aimed at discovering a wide range of PSMs that are effective in fending off viruses and their related vectors. These advances in understanding the potential role of PSMs in plant-virus interactions not only serve as a cornerstone for developing plant antiviral systems, but also highlight essential principles of biological control.

The immediate metabolomic effects of whole-genome duplication in the greater duckweed, Spirodela polyrhiza

PREMISE

In plants, whole-genome duplication (WGD) is a common mutation with profound evolutionary potential. Given the costs associated with a superfluous genome copy, polyploid establishment is enigmatic. However, in the right environment, immediate phenotypic changes following WGD can facilitate establishment. Metabolite abundances are the direct output of the cell's regulatory network and determine much of the impact of environmental and genetic change on the phenotype. While it is well known that an increase in the bulk amount of genetic material can increase cell size, the impact of gene dosage multiplication on the metabolome remains largely unknown.

METHODS

We used untargeted metabolomics on four genetically distinct diploid-neoautotetraploid pairs of the greater duckweed, Spirodela polyrhiza, to investigate how WGD affects metabolite abundances per cell and per biomass.

RESULTS

Autopolyploidy increased metabolite levels per cell, but the response of individual metabolites varied considerably. However, the impact on metabolite level per biomass was restricted because the increased cell size reduced the metabolite concentration per cell. Nevertheless, we detected both quantitative and qualitative effects of WGD on the metabolome. Many effects were strain-specific, but some were shared by all four strains.

CONCLUSIONS

The nature and impact of metabolic changes after WGD depended strongly on the genotype. Dosage effects have the potential to alter the plant metabolome qualitatively and quantitatively, but were largely balanced out by the reduction in metabolite concentration due to an increase in cell size in this species.

Target control of linear directed networks based on the path cover problem

Securing complete control of complex systems comprised of tens of thousands of interconnected nodes holds immense significance across various fields, spanning from cell biology and brain science to human-engineered systems. However, depending on specific functional requirements, it can be more practical and efficient to focus on a pre-defined subset of nodes for control, a concept known as target control. While some methods have been proposed to find the smallest driver node set for target control, they either rely on heuristic approaches based on k-walk theory, lacking a guarantee of optimal solutions, or they are overly complex and challenging to implement in real-world networks. To address this challenge, we introduce a simple and elegant algorithm, inspired by the path cover problem, which efficiently identifies the nodes required to control a target node set within polynomial time. To practically apply the algorithm in real-world systems, we have selected several networks in which a specific set of nodes with functional significance can be designated as a target control set. The analysed systems include the complete connectome of the nematode worm C. elegans, the recently disclosed connectome of the Drosophila larval brain, as well as dozens of genome-wide metabolic networks spanning major plant lineages. The target control analysis shed light on distinctions between neural systems in nematode worms and larval brain insects, particularly concerning the number of nodes necessary to regulate specific functional systems. Furthermore, our analysis uncovers evolutionary trends within plant lineages, notably when examining the proportion of nodes required to control functional pathways.

SCPL acyltransferases catalyze the metabolism of chlorogenic acid during purple coneflower seed germination

The metabolism of massively accumulated chlorogenic acid is crucial for the successful germination of purple coneflower (Echinacea purpurea (L.) Menoch). A serine carboxypeptidase-like (SCPL) acyltransferase (chicoric acid synthase, CAS) utilizes chlorogenic acid to produce chicoric acid during germination. However, it seems that the generation of chicoric acid lags behind the decrease in chlorogenic acid, suggesting an earlier route of chlorogenic acid metabolism. We discovered another chlorogenic acid metabolic product, 3,5-dicaffeoylquinic acid, which is produced before chicoric acid, filling the lag phase. Then, we identified two additional typical clade IA SCPL acyltransferases, named chlorogenic acid condensing enzymes (CCEs), that catalyze the biosynthesis of 3,5-dicaffeoylquinic acid from chlorogenic acid with different kinetic characteristics. Chlorogenic acid inhibits radicle elongation in a dose-dependent manner, explaining the potential biological role of SCPL acyltransferases-mediated continuous chlorogenic acid metabolism during germination. Both CCE1 and CCE2 are highly conserved among Echinacea species, supporting the observed metabolism of chlorogenic acid to 3,5-dicaffeoylquinic acid in two Echinacea species without chicoric acid accumulation. The discovery of SCPL acyltransferase involved in the biosynthesis of 3,5-dicaffeoylquinic acid suggests convergent evolution. Our research clarifies the metabolism strategy of chlorogenic acid in Echinacea species and provides more insight into plant metabolism.

Using Synthetic Biology to Understand the Function of Plant Specialized Metabolites

Plant specialized metabolites (PSMs) are variably distributed across taxa, tissues, and ecological contexts; this variability has inspired many theories about PSM function, which, to date, remain poorly tested because predictions have outpaced the available data. Advances in mass spectrometry-based metabolomics have enabled unbiased PSM profiling, and molecular biology techniques have produced PSM-free plants; the combination of these methods has accelerated our understanding of the complex ecological roles that PSMs play in plants. Synthetic biology techniques and workflows are producing high-value, structurally complex PSMs in quantities and purities sufficient for both medicinal and functional studies. These workflows enable the reengineering of PSM transport, externalization, structural diversity, and production in novel taxa, facilitating rigorous tests of long-standing theoretical predictions about why plants produce so many different PSMs in particular tissues and ecological contexts. Plants use their chemical prowess to solve ecological challenges, and synthetic biology workflows are accelerating our understanding of these evolved functions.

Multisubstrate specificity shaped the complex evolution of the aminotransferase family across the tree of life

Aminotransferases (ATs) are an ancient enzyme family that play central roles in core nitrogen metabolism, essential to all organisms. However, many of the AT enzyme functions remain poorly defined, limiting our fundamental understanding of the nitrogen metabolic networks that exist in different organisms. Here, we traced the deep evolutionary history of the AT family by analyzing AT enzymes from 90 species spanning the tree of life (ToL). We found that each organism has maintained a relatively small and constant number of ATs. Mapping the distribution of ATs across the ToL uncovered that many essential AT reactions are carried out by taxon-specific AT enzymes due to wide-spread nonorthologous gene displacements. This complex evolutionary history explains the difficulty of homology-based AT functional prediction. Biochemical characterization of diverse aromatic ATs further revealed their broad substrate specificity, unlike other core metabolic enzymes that evolved to catalyze specific reactions today. Interestingly, however, we found that these AT enzymes that diverged over billion years share common signatures of multisubstrate specificity by employing different nonconserved active site residues. These findings illustrate that AT family enzymes had leveraged their inherent substrate promiscuity to maintain a small yet distinct set of multifunctional AT enzymes in different taxa. This evolutionary history of versatile ATs likely contributed to the establishment of robust and diverse nitrogen metabolic networks that exist throughout the ToL. The study provides a critical foundation to systematically determine diverse AT functions and underlying nitrogen metabolic networks across the ToL.

Molecular-level carbon traits underlie the multidimensional fine root economics space

Carbon influences the evolution and functioning of plants and their roots. Previous work examining a small number of commonly measured root traits has revealed a global multidimensionality of the resource economics traits in fine roots considering carbon as primary currency but without considering the diversity of carbon-related traits. To address this knowledge gap, we use data from 66 tree species from a tropical forest to illustrate that root economics space co-varies with a novel molecular-level traits space based on nuclear magnetic resonance. Thinner fine roots exhibit higher proportions of carbohydrates and lower diversity of molecular carbon than thicker roots. Mass-denser fine roots have more lignin and aromatic carbon compounds but less bioactive carbon compounds than lighter roots. Thus, the transition from thin to thick fine roots implies a shift in the root carbon economy from 'do-it-yourself' soil exploration to collaboration with mycorrhizal fungi, while the shift from light to dense fine roots emphasizes a shift from acquisitive to conservative root strategy. We reveal a previously undocumented role of molecular-level carbon traits that potentially undergird the multidimensional root economics space. This finding offers new molecular insight into the diversity of root form and function, which is fundamental to our understanding of plant evolution, species coexistence and adaptations to heterogeneous environments.

Arogenate dehydratase isoforms strategically deregulate phenylalanine biosynthesis in Akebia trifoliata

Arogenate dehydratase (ADT) is key for phenylalanine (Phe) biosynthesis in plants. To examine ADT components and function in Akebia trifoliata, a representative of Ranunculaceae, we first identified eight ADTs (AktADT1-8, encoding sequences varying from 1032 to 1962 bp) in the A. trifoliata reference genome and five proteins (AktADT1, AktADT4, AktADT7, AktADT8 and AktADT8s) with moonlighting prephenate dehydratase (PDT) activity and Km values varying from 0.43 to 2.17 mM. Structurally, two basic residue combinations (Val314/Ala317 and Ala314/Val317) in the PAC domain are essential for the moonlighting PDT activity of ADTs. Functionally, AktADT4 and AktADT8 successfully restored the wild-type phenotype of pha2, a knockout mutant of Saccharomyces cerevisiae. In addition, AktADTs are ubiquitously expressed, but their expression levels are tissue specific, and the half maximal inhibitory concentration (IC50) of Phe for AktADTs ranged from 49.81 to 331.17 μM. Both AktADT4 and AktADT8 and AktADT8s localized to chloroplast stromules and the cytosol, respectively, while the remaining AktADTs localized to the chloroplast stroma. These findings suggest that various strategies exist for regulating Phe biosynthesis in A. trifoliata. This provides a reasonable explanation for the high Phe content and insights for further genetic improvement of the edible fruits of A. trifoliata.

Mechanistic docking in terpene synthases using EnzyDock

Terpene Synthases (TPS) catalyze the formation of multicyclic, complex terpenes and terpenoids from linear substrates. Molecular docking is an important research tool that can further our understanding of TPS multistep mechanisms and guide enzyme design. Standard docking programs are not well suited to tackle the unique challenges of TPS, like the many chemical steps which form multiple stereo-centers, the weak dispersion interactions between the isoprenoid chain and the hydrophobic region of the active site, description of carbocation intermediates, and finding mechanistically meaningful sets of docked poses. To address these and other unique challenges, we developed the multistate, multiscale docking program EnzyDock and used it to study many TPS and other enzymes. In this review we discuss the unique challenges of TPS, the special features of EnzyDock developed to address these challenges and demonstrate its successful use in ongoing research on the bacterial TPS CotB2.

Molecular evolution and characterization of type III polyketide synthase gene family in Aquilaria sinensis

2-(2-Phenylethyl) chromone (PEC) and its derivatives are markers of agarwood formation and are also related to agarwood quality. However, the biosynthetic and regulatory mechanisms of PECs still remain mysterious. Several studies suggested that type III polyketide synthases (PKSs) contribute to PEC biosynthesis in Aquilaria sinensis. Furthermore, systematic studies on the evolution of PKSs in A. sinensis have rarely been reported. Herein, we comprehensively analyzed PKS genes from 12 plant genomes and characterized the AsPKSs in detail. A unique branch contained only AsPKS members was identified through evolutionary analysis, including AsPKS01 that was previously indicated to participate in PEC biosynthesis. AsPKS07 and AsPKS08, two tandem-duplicated genes of AsPKS01 and lacking orthologous genes in evolutionary models, were selected for their transient expression in the leaves of Nicotiana benthamiana. Subsequently, PECs were detected in the extracts of N. benthamiana leaves, suggesting that AsPKS07 and AsPKS08 promote PEC biosynthesis. The interaction between the promoters of AsPKS07, AsPKS08 and five basic leucine zippers (bZIPs) from the S subfamily indicated that their transcripts could be regulated by these transcription factors (TFs) and might further contribute to PECs biosynthesis in A. sinensis. Our findings provide valuable insights into the molecular evolution of the PKS gene family in A. sinensis and serve as a foundation for advancing PEC production through the bioengineering of gene clusters. Ultimately, this contribution is expected to shed light on the mechanism underlying agarwood formation.

Three types of enzymes complete the furanocoumarins core skeleton biosynthesis in Angelica sinensis

Furanocoumarins (FCs) are widely distributed secondary metabolites found in higher plants, including Apiaceae, Rutaceae, Moraceae, and Fabaceae. They play a crucial role in the physiological functions of plants and are well-known for their diverse pharmacological activities. As a representative plant of the Apiaceae family, Angelica sinensis is highly valued for its medicinal properties and FCs are one of the main ingredients of A. sinensis. However, the biosynthetic mechanism of FCs in A. sinensis remains poorly understood. In this study, we successfully cloned and verified three types of enzymes using genome analysis and in vitro functional verification, which complete the biosynthesis of the FCs core skeleton in A. sinensis. It includes a p-coumaroyl CoA 2'-hydroxylase (AsC2'H) responsible for umbelliferone formation, two UbiA prenyltransferases (AsPT1 and AsPT2) that convert umbelliferone to demethylsuberosin (DMS) and osthenol, respectively, and two CYP736 subfamily cyclases (AsDC and AsOD) that catalyze the formation of FCs core skeleton. Interestingly, AsOD was demonstrated to be a bifunctional cyclase and could catalyze both DMS and osthenol, but had a higher affinity to osthenol. The characterization of these enzymes elucidates the molecular mechanism of FCs biosynthesis, providing new insights and technologies for understanding the diverse origins of FCs biosynthesis.

An intramolecular macrocyclase in plant ribosomal peptide biosynthesis

The biosynthetic dogma of ribosomally synthesized and posttranslationally modified peptides (RiPP) involves enzymatic intermolecular modification of core peptide motifs in precursor peptides. The plant-specific BURP-domain protein family, named after their four founding members, includes autocatalytic peptide cyclases involved in the biosynthesis of side-chain-macrocyclic plant RiPPs. Here we show that AhyBURP, a representative of the founding Unknown Seed Protein-type BURP-domain subfamily, catalyzes intramolecular macrocyclizations of its core peptide during the sequential biosynthesis of monocyclic lyciumin I via glycine-tryptophan crosslinking and bicyclic legumenin via glutamine-tyrosine crosslinking. X-ray crystallography of AhyBURP reveals the BURP-domain fold with two type II copper centers derived from a conserved stapled-disulfide and His motif. We show the macrocyclization of lyciumin-C(sp3)-N-bond formation followed by legumenin-C(sp3)-O-bond formation requires dioxygen and radical involvement based on enzyme assays in anoxic conditions and isotopic labeling. Our study expands enzymatic intramolecular modifications beyond catalytic moiety and chromophore biogenesis to RiPP biosynthesis.

Functional analysis of CYP71AV1 reveals the evolutionary landscape of artemisinin biosynthesis

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

The evolutionary origin of naturally occurring intermolecular Diels-Alderases from Morus alba

Biosynthetic enzymes evolutionarily gain novel functions, thereby expanding the structural diversity of natural products to the benefit of host organisms. Diels-Alderases (DAs), functionally unique enzymes catalysing [4 + 2] cycloaddition reactions, have received considerable research interest. However, their evolutionary mechanisms remain obscure. Here, we investigate the evolutionary origins of the intermolecular DAs in the biosynthesis of Moraceae plant-derived Diels-Alder-type secondary metabolites. Our findings suggest that these DAs have evolved from an ancestor functioning as a flavin adenine dinucleotide (FAD)-dependent oxidocyclase (OC), which catalyses the oxidative cyclisation reactions of isoprenoid-substituted phenolic compounds. Through crystal structure determination, computational calculations, and site-directed mutagenesis experiments, we identified several critical substitutions, including S348L, A357L, D389E and H418R that alter the substrate-binding mode and enable the OCs to gain intermolecular DA activity during evolution. This work provides mechanistic insights into the evolutionary rationale of DAs and paves the way for mining and engineering new DAs from other protein families.

A trait-based framework linking the soil metabolome to plant-soil feedbacks

By modifying the biotic and abiotic properties of the soil, plants create soil legacies that can affect vegetation dynamics through plant-soil feedbacks (PSF). PSF are generally attributed to reciprocal effects of plants and soil biota, but these interactions can also drive changes in the identity, diversity and abundance of soil metabolites, leading to more or less persistent soil chemical legacies whose role in mediating PSF has rarely been considered. These chemical legacies may interact with microbial or nutrient legacies to affect species coexistence. Given the ecological importance of chemical interactions between plants and other organisms, a better understanding of soil chemical legacies is needed in community ecology. In this Viewpoint, we aim to: highlight the importance of belowground chemical interactions for PSF; define and integrate soil chemical legacies into PSF research by clarifying how the soil metabolome can contribute to PSF; discuss how functional traits can help predict these plant-soil interactions; propose an experimental approach to quantify plant responses to the soil solution metabolome; and describe a testable framework relying on root economics and seed dispersal traits to predict how plant species affect the soil metabolome and how they could respond to soil chemical legacies.

Analysing a Group of Homologous BAHD Enzymes Provides Insights into the Evolutionary Transition of Rosmarinic Acid Synthases from Hydroxycinnamoyl-CoA:Shikimate/Quinate Hydroxycinnamoyl Transferases

The interplay of various enzymes and compounds gives rise to the intricate secondary metabolic networks observed today. However, the current understanding of their formation and expansion remains limited. BAHD acyltransferases play important roles in the biosynthesis of numerous significant secondary metabolites. In plants, they are widely distributed and exhibit a diverse range of activities. Among them, rosmarinic acid synthase (RAS) and hydroxycinnamoyl-CoA:shikimate/quinate hydroxycinnamoyl transferase (HCT) have gained significant recognition and have been extensively investigated as prominent members of the BAHD acyltransferase family. Here, we conducted a comprehensive study on a unique group of RAS homologous enzymes in Mentha longifolia that display both catalytic activities and molecular features similar to HCT and Lamiaceae RAS. Subsequent phylogenetic and comparative genome analyses revealed their derivation from expansion events within the HCT gene family, indicating their potential as collateral branches along the evolutionary trajectory, leading to Lamiaceae RAS while still retaining certain ancestral vestiges. This discovery provides more detailed insights into the evolution from HCT to RAS. Our collective findings indicate that gene duplication is the driving force behind the observed evolutionary pattern in plant-specialized enzymes, which probably originated from ancestral enzyme promiscuity and were subsequently shaped by principles of biological adaptation.

The evolution of chemodiversity in plants-From verbal to quantitative models

Plants harbour a great chemodiversity, that is diversity of specialised metabolites (SMs), at different scales. For instance, individuals can produce a large number of SMs, and populations can differ in their metabolite composition. Given the ecological and economic importance of plant chemodiversity, it is important to understand how it arises and is maintained over evolutionary time. For other dimensions of biodiversity, that is species diversity and genetic diversity, quantitative models play an important role in addressing such questions. Here, we provide a synthesis of existing hypotheses and quantitative models, that is mathematical models and computer simulations, for the evolution of plant chemodiversity. We describe each model's ingredients, that is the biological processes that shape chemodiversity, the scales it considers and whether it has been formalized as a quantitative model. Although we identify several quantitative models, not all are dynamic and many influential models have remained verbal. To fill these gaps, we outline our vision for the future of chemodiversity modelling. We identify quantitative models used for genetic variation that may be adapted for chemodiversity, and we present a flexible framework for the creation of individual-based models that address different scales of chemodiversity and combine different ingredients that bring this chemodiversity about.

Insights into the differences related to the resistance mechanisms to the highly toxic fruit Hippomane mancinella (Malpighiales: Euphorbiaceae) between the larvae of the sister species Anastrepha acris and Anastrepha ludens (Diptera: Tephritidae) through comparative transcriptomics

The Manchineel, Hippomane mancinella ("Death Apple Tree") is one of the most toxic fruits worldwide and nevertheless is the host plant of the monophagous fruit fly species Anastrepha acris (Diptera: Tephritidae). Here we aimed at elucidating the detoxification mechanisms in larvae of A. acris reared on a diet enriched with the toxic fruit (6% lyophilizate) through comparative transcriptomics. We compared the performance of A. acris larvae with that of the sister species A. ludens, a highly polyphagous pest species that is unable to infest H. mancinella in nature. The transcriptional alterations in A. ludens were significantly greater than in A. acris. We mainly found two resistance mechanisms in both species: structural, activating cuticle protein biosynthesis (chitin-binding proteins likely reducing permeability to toxic compounds in the intestine), and metabolic, triggering biosynthesis of serine proteases and xenobiotic metabolism activation by glutathione-S-transferases and cytochrome P450 oxidoreductase. Some cuticle proteins and serine proteases were not orthologous between both species, suggesting that in A. acris, a structural resistance mechanism has been selected allowing specialization to the highly toxic host plant. Our results represent a nice example of how two phylogenetically close species diverged over recent evolutionary time related to resistance mechanisms to plant secondary metabolites.

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.

Genome-wide identification and expression of BAHD acyltransferase gene family shed novel insights into the regulation of linalyl acetate and lavandulyl acetate in lavender

The BAHD acyltransferase superfamily has a variety of biological functions, especially in catalyzing the synthesis of ester compounds and improving plant stress resistance. Linalyl acetate and lavandulyl acetate, the most important volatile esters in lavender, are generated by LaBAHDs. However, the systematic identification, expression characteristics of LaBAHD genes and their correlations with ester formation remain elusive. Here, 166 LaBAHD genes were identified from the lavender genome. Based on detailed phylogenetic analysis, the LaBAHD family genes were divided into five groups, among which the LaBAHDs involved in volatile ester biosynthesis belong to the IIIa and Va clades. Whole-genome duplications (WGDs) and tandem duplications (TDs) jointly drive the expansion of LaBAHD superfamily. The promoter regions of LaBAHDs contained a variety of stress- and hormone-related motifs, as well as binding sites with five types of transcription factors (TFs). Then, linalyl acetate- and lavandulyl acetate-regulated coexpression modules were established and some candidate TFs that may function in inducing ester formation were identified. Based on the correlation analysis between the ester contents and expression profiles of BAHD genes in different tissues, five candidate genes were screened for further examination. Drought, salt and MeJA treatments increased the accumulation of linalyl acetate and lavandulyl acetate, and induced the expression of LaBAHDs. Our results indicated that LaBAHD57, LaBAHD63, LaBAHD104, LaBAHD105 and LaBAHD119 are crucial candidate genes involved in linalyl acetate and lavandulyl acetate biosynthesis. Our findings offer a theoretical foundation for further studying the specific biological functions of LaBAHD family and improving the quality of lavender essential oil.

Studying Plant Specialized Metabolites Using Computational Metabolomics Strategies

Plant specialized metabolites have diversified vastly over the course of plant evolution, and they are considered key players in complex interactions between plants and their environment. The chemical diversity of these metabolites has been widely explored and utilized in agriculture and crop enhancement, the food industry, and drug development, among other areas. However, the immensity of the plant metabolome can make its exploration challenging. Here we describe a protocol for exploring plant specialized metabolites that combines high-resolution mass spectrometry and computational metabolomics strategies, including molecular networking, identification of structural motifs, as well as prediction of chemical structures and metabolite classes.

The R2R3-MYB transcription factor EVER controls the emission of petunia floral volatiles by regulating epicuticular wax biosynthesis in the petal epidermis

The epidermal cells of petunia (Petunia × hybrida) flowers are the main site of volatile emission. However, the mechanisms underlying the release of volatiles into the environment are still being explored. Here, using cell-layer-specific transcriptomic analysis, reverse genetics by virus-induced gene silencing and clustered regularly interspaced short palindromic repeat (CRISPR), and metabolomics, we identified EPIDERMIS VOLATILE EMISSION REGULATOR (EVER)-a petal adaxial epidermis-specific MYB activator that affects the emission of volatiles. To generate ever knockout lines, we developed a viral-based CRISPR/Cas9 system for efficient gene editing in plants. These knockout lines, together with transient-suppression assays, revealed EVER's involvement in the repression of low-vapor-pressure volatiles. Internal pools and annotated scent-related genes involved in volatile production and emission were not affected by EVER. RNA-Seq analyses of petals of ever knockout lines and EVER-overexpressing flowers revealed enrichment in wax-related biosynthesis genes. Liquid chromatography/gas chromatography-MS analyses of petal epicuticular waxes revealed substantial reductions in wax loads in ever petals, particularly of monomers of fatty acids and wax esters. These results implicate EVER in the emission of volatiles by fine-tuning the composition of petal epicuticular waxes. We reveal a petunia MYB regulator that interlinks epicuticular wax composition and volatile emission, thus unraveling a regulatory layer in the scent-emission machinery in petunia flowers.

Systematic thermal analysis of the Arabidopsis proteome: Thermal tolerance, organization, and evolution

Understanding the thermal stability of the plant proteome in the context of the native cellular environment would aid the design of crops with high thermal tolerance, but only limited such data are available. Here, we applied quantitative mass spectrometry to profile the thermal stability of the Arabidopsis proteome and identify thermo-sensitive and thermo-resilient protein networks in Arabidopsis, providing a basis for understanding heat-induced damage. We also show that the similarities of the protein-melting curves can be used as a proxy to evaluate system-wide protein-protein interactions in non-engineered plants and enable the identification of transient interactions exhibited by metabolons in the context of the cellular environment. Finally, we report a systematic comparison of the thermal stability of paralogs in Arabidopsis to aid the investigation and understanding of gene duplication and protein evolution. Taken together, our results could have broad implications for the fields of plant thermal tolerance, plant protein assemblies, and evolution.

Jasmonates Coordinate Secondary with Primary Metabolism

Jasmonates (JAs), including jasmonic acid (JA), its precursor 12-oxo-phytodienoic acid (OPDA) and its derivatives jasmonoyl-isoleucine (JA-Ile), methyl jasmonate (MeJA), cis-jasmone (CJ) and other oxylipins, are important in the regulation of a range of ecological interactions of plants with their abiotic and particularly their biotic environments. Plant secondary/specialized metabolites play critical roles in implementing these ecological functions of JAs. Pathway and transcriptional regulation analyses have established a central role of JA-Ile-mediated core signaling in promoting the biosynthesis of a great diversity of secondary metabolites. Here, we summarized the advances in JAs-induced secondary metabolites, particularly in secondary metabolites induced by OPDA and volatile organic compounds (VOCs) induced by CJ through signaling independent of JA-Ile. The roles of JAs in integrating and coordinating the primary and secondary metabolism, thereby orchestrating plant growth-defense tradeoffs, were highlighted and discussed. Finally, we provided perspectives on the improvement of the adaptability and resilience of plants to changing environments and the production of valuable phytochemicals by exploiting JAs-regulated secondary metabolites.

Structural and functional diversity in plant specialized metabolism signals and products: The case of oxylipins and triterpenes

Metabolic enzymes tend to evolve towards catalytic efficacy, precision and speed. This seems particularly true for ancient and conserved enzymes involved in fundamental cellular processes that are present virtually in every cell and organism and converting and producing relatively limited metabolite numbers. Nevertheless, sessile organisms like plants have an astonishing repertoire of specific (specialized) metabolites that, by numbers and chemical complexity, by far exceed primary metabolites. Most theories agree that early gene duplication, subsequent positive selection and diversifying evolution have allowed relaxed selection of duplicated metabolic genes, thus facilitating the accumulation of mutations that could broaden substrate/product specificity and lower activation barriers and kinetics. Here, we use oxylipins, oxygenated fatty acids of plastidial origin to which the phytohormone jasmonate belongs, and triterpenes, a large group of specialized metabolites whose biosynthesis is often elicited by jasmonates, to showcase the structural and functional diversity of chemical signals and products in plant metabolism.

Transcriptomic and metabolomic perspectives for the growth of alfalfa (Medicago sativa L.) seedlings with the effect of vanadium exposure

Hitherto, the effect of vanadium on higher plant growth remains an open topic. Therefore, nontargeted metabolomic and RNA-Seq profiling were implemented to unravel the possible alteration in alfalfa seedlings subjected to 0.1 mg L^-1 (B group) and 0.5 mg L^-1 (C group) pentavalent vanadium [(V(V)] versus control (A group) in this study. Results revealed that vanadium exposure significantly altered some pivotal transcripts and metabolites. The number of differentially expressed genes (DEGs) markedly up- and down-regulated was 21 and 23 in B_vs_A, 27 and 33 in C_vs_A, and 24 and 43 in C_vs_B, respectively. The number for significantly up- and down-regulated differential metabolites was 17 and 15 in B_vs_A, 43 and 20 in C_vs_A, and 24 and 16 in C_vs_B, respectively. Metabolomics and transcriptomics co-analysis characterized three significantly enriched metabolic pathways in C_vs_A comparing group, viz., α-linolenic acid metabolism, flavonoid biosynthesis, and phenylpropanoid biosynthesis, from which some differentially expressed genes and differential metabolites participated. The metabolite of traumatic acid in α-linolenic acid metabolism and apigenin in flavonoid biosynthesis were markedly upregulated, while phenylalanine in phenylpropanoid biosynthesis was remarkably downregulated. The genes of allene oxide cyclase (AOC) and acetyl-CoA acyltransferase (fadA) in α-linolenic acid metabolism, and chalcone synthase (CHS), flavonoid 3'-monooxygenase (CYP75B1), and flavonol synthase (FLS) in flavonoid biosynthesis, and caffeoyl-CoA O-methyltransferase (CCoAOMT) in phenylpropanoid biosynthesis were significantly downregulated. While shikimate O-hydroxycinnamoyltransferase (HCT) in flavanoid and phenylpropanoid biosynthesis were conspicuously upregulated. Briefly, vanadium exposure induces a readjustment yielding in metabolite and the correlative synthetic precursors (transcripts/unigenes) in some branched metabolic pathways. This study provides a practical and in-depth perspective from transcriptomics and metabolomics in investigating the effects conferred by vanadium on plant growth and development.

Plant terpene specialized metabolism: complex networks or simple linear pathways?

From the perspectives of pathway evolution, discovery and engineering of plant specialized metabolism, the nature of the biosynthetic routes represents a critical aspect. Classical models depict biosynthesis typically from an end-point angle and as linear, for example, connecting central and specialized metabolism. As the number of functionally elucidated routes increased, the enzymatic foundation of complex plant chemistries became increasingly well understood. The perception of linear pathway models has been severely challenged. With a focus on plant terpenoid specialized metabolism, we review here illustrative examples supporting that plants have evolved complex networks driving chemical diversification. The completion of several diterpene, sesquiterpene and monoterpene routes shows complex formation of scaffolds and their subsequent functionalization. These networks show that branch points, including multiple sub-routes, mean that metabolic grids are the rule rather than the exception. This concept presents significant implications for biotechnological production.

Elucidating the Mesocarp Drupe Transcriptome of Açai (Euterpe oleracea Mart.): An Amazonian Tree Palm Producer of Bioactive Compounds

Euterpe oleracea palm, endemic to the Amazon region, is well known for açai, a fruit violet beverage with nutritional and medicinal properties. During E. oleracea fruit ripening, anthocyanin accumulation is not related to sugar production, contrarily to grape and blueberry. Ripened fruits have a high content of anthocyanins, isoprenoids, fibers, and proteins, and are poor in sugars. E. oleracea is proposed as a new genetic model for metabolism partitioning in the fruit. Approximately 255 million single-end-oriented reads were generated on an Ion Proton NGS platform combining fruit cDNA libraries at four ripening stages. The de novo transcriptome assembly was tested using six assemblers and 46 different combinations of parameters, a pre-processing and a post-processing step. The multiple k-mer approach with TransABySS as an assembler and Evidential Gene as a post-processer have shown the best results, with an N50 of 959 bp, a read coverage mean of 70x, a BUSCO complete sequence recovery of 36% and an RBMT of 61%. The fruit transcriptome dataset included 22,486 transcripts representing 18 Mbp, of which a proportion of 87% had significant homology with other plant sequences. Approximately 904 new EST-SSRs were described, and were common and transferable to Phoenix dactylifera and Elaeis guineensis, two other palm trees. The global GO classification of transcripts showed similar categories to that in P. dactylifera and E. guineensis fruit transcriptomes. For an accurate annotation and functional description of metabolism genes, a bioinformatic pipeline was developed to precisely identify orthologs, such as one-to-one orthologs between species, and to infer multigenic family evolution. The phylogenetic inference confirmed an occurrence of duplication events in the Arecaceae lineage and the presence of orphan genes in E. oleracea. Anthocyanin and tocopherol pathways were annotated entirely. Interestingly, the anthocyanin pathway showed a high number of paralogs, similar to in grape, whereas the tocopherol pathway exhibited a low and conserved gene number and the prediction of several splicing forms. The release of this exhaustively annotated molecular dataset of E. oleracea constitutes a valuable tool for further studies in metabolism partitioning and opens new great perspectives to study fruit physiology with açai as a model.

A developmental gradient reveals biosynthetic pathways to eukaryotic toxins in monocot geophytes

Numerous eukaryotic toxins that accumulate in geophytic plants are valuable in the clinic, yet their biosynthetic pathways have remained elusive. A lead example is the >150 Amaryllidaceae alkaloids (AmAs) including galantamine, an FDA-approved treatment for Alzheimer's disease. We show that while AmAs accumulate to high levels in many tissues in daffodils, biosynthesis is localized to nascent, growing tissue at the base of leaves. A similar trend is found for the production of steroidal alkaloids (e.g. cyclopamine) in corn lily. This model of active biosynthesis enabled elucidation of a complete set of biosynthetic genes for the production of AmAs. Taken together, our work sheds light on the developmental and enzymatic logic of diverse alkaloid biosynthesis in daffodil. More broadly, it suggests a paradigm for biosynthesis regulation in monocot geophytes where plants are protected from herbivory through active charging of newly formed cells with eukaryotic toxins that persist as aboveground tissue develops.

Biochemistry and genetics of floral scent: a historical perspective

Floral scent plays a crucial role in the reproductive process of many plants. Humans have been fascinated by floral scents throughout history, and have transported and traded floral scent products for which they have found multiple uses, such as in food additives, hygiene and perfume products, and medicines. Yet the scientific study of how plants synthesize floral scent compounds began later than studies on most other major plant metabolites, and the first report of the characterization of an enzyme responsible for the synthesis of a floral scent compound, namely linalool in Clarkia breweri, a California annual, appeared in 1994. In the almost 30 years since, enzymes and genes involved in the synthesis of hundreds of scent compounds from multiple plant species have been described. This review recapitulates this history and describes the major findings relating to the various aspects of floral scent biosynthesis and emission, including genes and enzymes and their evolution, storage and emission of scent volatiles, and the regulation of the biochemical processes.

Perplexing Polyphenolics: The Isolations, Syntheses, Reappraisals, and Bioactivities of Flavonoids, Isoflavonoids, and Neoflavonoids from 2016 to 2022

Flavonoids, isoflavonoids, neoflavonoids, and their various subcategories are polyphenolics-an extensive class of natural products. These compounds are bioactive and display multiple activities, including anticancer, antibacterial, antiviral, antioxidant, and neuroprotective activities. Thus, these compounds can serve as leads for therapeutic agents or targets for complex synthesis; they are coveted and routinely isolated, characterized, biologically evaluated, and synthesized. However, data regarding the compounds' sources, isolation procedures, structural novelties, bioactivities, and synthetic schemes are often dispersed and complex, a dilemma this review aims to address. To serve as an easily accessible guide for researchers wanting to apprise themselves of the latest advancements in this subfield, this review summarizes seventy-six (76) articles published between 2016 and 2022 that detail the isolation and characterization of two hundred and forty-nine (249) novel compounds, the total and semisyntheses of thirteen (13) compounds, and reappraisals of the structures of twenty (20) previously reported compounds and their bioactivities. This article also discusses new synthetic methods and enzymes capable of producing or modifying flavonoids, isoflavonoids, or neoflavonoids.

Loss of MYB34 Transcription Factor Supports the Backward Evolution of Indole Glucosinolate Biosynthesis in a Subclade of the Camelineae Tribe and Releases the Feedback Loop in This Pathway in Arabidopsis

Glucosinolates are specialized defensive metabolites characteristic of the Brassicales order. Among them, aliphatic and indolic glucosinolates (IGs) are usually highly abundant in species from the Brassicaceae family. The exceptions this trend are species representing a subclade of the Camelineae tribe, including Capsella and Camelina genera, which have reduced capacity to produce and metabolize IGs. Our study addresses the contribution of specific glucosinolate-related myeloblastosis (MYB) transcription factors to this unprecedented backward evolution of IG biosynthesis. To this end, we performed phylogenomic and functional studies of respective MYB proteins. The obtained results revealed weakened conservation of glucosinolate-related MYB transcription factors, including loss of functional MYB34 protein, in the investigated species. We showed that the introduction of functional MYB34 from Arabidopsis thaliana partially restores IG biosynthesis in Capsella rubella, indicating that the loss of this transcription factor contributes to the backward evolution of this metabolic pathway. Finally, we performed an analysis of the impact of particular myb mutations on the feedback loop in IG biosynthesis, which drives auxin overproduction, metabolic dysregulation and strong growth retardation caused by mutations in IG biosynthetic genes. This uncovered the unique function of MYB34 among IG-related MYBs in this feedback regulation and consequently in IG conservation in Brassicaceae plants.

Two O-Methyltransferases from Phylogenetically Unrelated Cow Parsley (Anthriscus sylvestris) and Hinoki-Asunaro (Thujopsis dolabrata var. hondae) as a Signature of Lineage-Specific Evolution in Lignan Biosynthesis

O-Methyltransferases (OMTs) play important roles in antitumor lignan biosynthesis. To date, six OMTs catalyzing the methylation of dibenzylbutyrolactone lignans as biosynthetic precursors of antitumor lignans have been identified. However, there is still no systematic understanding of the diversity and regularity of the biosynthetic mechanisms among various plant lineages. Herein, we report the characterization of two OMTs from Anthriscus sylvestris and Thujopsis dolabrata var. hondae [designated as AsSecoNorYatein (SNY) OMT and TdSNYOMT] together with the six known OMTs to evaluate their diversity and regularity. Although A. sylvestris 5-O-methylthujaplicatin (SecoNorYatein) and 4-O-demethylyatein (NorYatein) OMT (AsSNYOMT) and TdSNYOMT accept 5-O-methylthujaplicatin and 4-O-demethylyatein as substrates, phylogenetic analysis indicated that these two OMTs shared low amino acid sequence identity, 33.8%, indicating a signature of parallel evolution. The OMTs and the six previously identified OMTs were found to be diverse in terms of their substrate specificity, regioselectivity and amino acid sequence identity, indicating independent evolution in each plant species. Meanwhile, two-entropy analysis detected four amino acid residues as being specifically acquired by dibenzylbutyrolactone lignan OMTs. Site-directed mutation of AsSNYOMT indicated that two of them contributed specifically to 5-O-methylthujaplicatin methylation. The results provide a new example of parallel evolution and the diversity and regularity of OMTs in plant secondary (specialized) metabolism.

Phylogenomic analyses across land plants reveals motifs and coexpression patterns useful for functional prediction in the BAHD acyltransferase family

The BAHD acyltransferase family is one of the largest enzyme families in flowering plants, containing dozens to hundreds of genes in individual genomes. Highly prevalent in angiosperm genomes, members of this family contribute to several pathways in primary and specialized metabolism. In this study, we performed a phylogenomic analysis of the family using 52 genomes across the plant kingdom to gain deeper insights into its functional evolution and enable function prediction. We found that BAHD expansion in land plants was associated with significant changes in various gene features. Using pre-defined BAHD clades, we identified clade expansions in different plant groups. In some groups, these expansions coincided with the prominence of metabolite classes such as anthocyanins (flowering plants) and hydroxycinnamic acid amides (monocots). Clade-wise motif-enrichment analysis revealed that some clades have novel motifs fixed on either the acceptor or the donor side, potentially reflecting historical routes of functional evolution. Co-expression analysis in rice and Arabidopsis further identified BAHDs with similar expression patterns, however, most co-expressed BAHDs belonged to different clades. Comparing BAHD paralogs, we found that gene expression diverges rapidly after duplication, suggesting that sub/neo-functionalization of duplicate genes occurs quickly via expression diversification. Analyzing co-expression patterns in Arabidopsis in conjunction with orthology-based substrate class predictions and metabolic pathway models led to the recovery of metabolic processes of most of the already-characterized BAHDs as well as definition of novel functional predictions for some uncharacterized BAHDs. Overall, this study provides new insights into the evolution of BAHD acyltransferases and sets up a foundation for their functional characterization.