Something Old, Something New: Conserved Enzymes and the Evolution of Novelty in Plant Specialized Metabolism

Genome-wide association identifies a BAHD acyltransferase activity that assembles an ester of glucuronosylglycerol and phenylacetic acid

Genome-wide association studies (GWAS) are an effective approach to identify new specialized metabolites and the genes involved in their biosynthesis and regulation. In this study, GWAS of Arabidopsis thaliana soluble leaf and stem metabolites identified alleles of an uncharacterized BAHD-family acyltransferase (AT5G57840) associated with natural variation in three structurally related metabolites. These metabolites were esters of glucuronosylglycerol, with one metabolite containing phenylacetic acid as the acyl component of the ester. Knockout and overexpression of AT5G57840 in Arabidopsis and heterologous overexpression in Nicotiana benthamiana and Escherichia coli demonstrated that it is capable of utilizing phenylacetyl-CoA as an acyl donor and glucuronosylglycerol as an acyl acceptor. We, thus, named the protein Glucuronosylglycerol Ester Synthase (GGES). Additionally, phenylacetyl glucuronosylglycerol increased in Arabidopsis CYP79A2 mutants that overproduce phenylacetic acid and was lost in knockout mutants of UDP-sulfoquinovosyl: diacylglycerol sulfoquinovosyl transferase, an enzyme required for glucuronosylglycerol biosynthesis and associated with glycerolipid metabolism under phosphate-starvation stress. GGES is a member of a well-supported clade of BAHD family acyltransferases that arose by duplication and neofunctionalized during the evolution of the Brassicales within a larger clade that includes HCT as well as enzymes that synthesize other plant-specialized metabolites. Together, this work extends our understanding of the catalytic diversity of BAHD acyltransferases and uncovers a pathway that involves contributions from both phenylalanine and lipid metabolism.

Overexpression of the Ginkgo biloba dihydroflavonol 4-reductase gene GbDFR6 results in the self-incompatibility-like phenotypes in transgenic tobacco

Although flavonoids play multiple roles in plant growth and development, the involvement in plant self-incompatibility (SI) have not been reported. In this research, the fertility of transgenic tobacco plants overexpressing the Ginkgo biloba dihydroflavonol 4-reductase gene, GbDFR6, were investigated. To explore the possible physiological defects leading to the failure of embryo development in transgenic tobacco plants, functions of pistils and pollen grains were examined. Transgenic pistils pollinated with pollen grains from another tobacco plants (either transgenic or wild-type), developed full of well-developed seeds. In contrast, in self-pollinated transgenic tobacco plants, pollen-tube growth was arrested in the upper part of the style, and small abnormal seeds developed without fertilization. Although the mechanism remains unclear, our research may provide a valuable method to create SI tobacco plants for breeding.

Cardiac glycosides protect wormseed wallflower (Erysimum cheiranthoides) against some, but not all, glucosinolate-adapted herbivores

The chemical arms race between plants and insects is foundational to the generation and maintenance of biological diversity. We asked how the evolution of a novel defensive compound in an already well-defended plant lineage impacts interactions with diverse herbivores. Erysimum cheiranthoides (Brassicaceae), which produces both ancestral glucosinolates and novel cardiac glycosides, served as a model.We analyzed gene expression to identify cardiac glycoside biosynthetic enzymes in E. cheiranthoides and characterized these enzymes via heterologous expression and CRISPR/Cas9 knockout. Using E. cheiranthoides cardiac glycoside-deficient lines, we conducted insect experiments in both the laboratory and field.EcCYP87A126 initiates cardiac glycoside biosynthesis via sterol side chain cleavage, and EcCYP716A418 has a role in cardiac glycoside hydroxylation. In EcCYP87A126 knockout lines, cardiac glycoside production was eliminated. Laboratory experiments with these lines revealed that cardiac glycosides were highly effective defenses against two species of glucosinolate-tolerant specialist herbivores but did not protect against all crucifer-feeding specialist herbivores in the field. Cardiac glycosides had lesser to no effect on two broad generalist herbivores.These results begin elucidation of the E. cheiranthoides cardiac glycoside biosynthetic pathway and demonstrate in vivo that cardiac glycoside production allows Erysimum to escape from some, but not all, specialist herbivores.

QT-GWAS: A novel method for unveiling biosynthetic loci affecting qualitative metabolic traits

Although the plant kingdom provides an enormous diversity of metabolites with potentially beneficial applications for humankind, a large fraction of these metabolites and their biosynthetic pathways remain unknown. Resolving metabolite structures and their biosynthetic pathways is key to gaining biological understanding and to allow metabolic engineering. In order to retrieve novel biosynthetic genes involved in specialized metabolism, we developed a novel untargeted method designated as qualitative trait GWAS (QT-GWAS) that subjects qualitative metabolic traits to a genome-wide association study, while the conventional metabolite GWAS (mGWAS) mainly considers the quantitative variation of metabolites. As a proof of the validity of QT-GWAS, 23 and 15 of the retrieved associations identified in Arabidopsis thaliana by QT-GWAS and mGWAS, respectively, were supported by previous research. Furthermore, seven gene-metabolite associations retrieved by QT-GWAS were confirmed in this study through reverse genetics combined with metabolomics and/or in vitro enzyme assays. As such, we established that CYTOCHROME P450 706A5 (CYP706A5) is involved in the biosynthesis of chroman derivatives, UDP-GLYCOSYLTRANSFERASE 76C3 (UGT76C3) is able to hexosylate guanine in vitro and in planta, and SULFOTRANSFERASE 202B1 (SULT202B1) catalyzes the sulfation of neolignans in vitro. Collectively, our study demonstrates that the untargeted QT-GWAS method can retrieve valid gene-metabolite associations at the level of enzyme-encoding genes, even new associations that cannot be found by the conventional mGWAS, providing a new approach for dissecting qualitative metabolic traits.

Specialized metabolism by trichome-enriched Rubisco and fatty acid synthase components

Acylsugars, specialized metabolites with defense activities, are secreted by trichomes of many solanaceous plants. Several acylsugar metabolic genes (AMGs) remain unknown. We previously reported multiple candidate AMGs. Here, using multiple approaches, we characterized additional AMGs. First, we identified differentially expressed genes between high- and low-acylsugar-producing F2 plants derived from a cross between cultivated tomato (Solanum lycopersicum) and a wild relative (Solanum pennellii), which produce acylsugars that are ∼1% and ∼20% of leaf dry weight, respectively. Expression levels of many known and candidate AMGs positively correlated with acylsugar amounts in F2 individuals. Next, we identified lycopersicum-pennellii putative orthologs with higher nonsynonymous to synonymous substitutions. These analyses identified four candidate genes, three of which showed enriched expression in stem trichomes compared to underlying tissues (shaved stems). Virus-induced gene silencing confirmed two candidates, Sopen05g009610 [beta-ketoacyl-(acyl-carrier-protein) reductase; fatty acid synthase component] and Sopen07g006810 (Rubisco small subunit), as AMGs. Phylogenetic analysis indicated that Sopen05g009610 is distinct from specialized metabolic cytosolic reductases but closely related to two capsaicinoid biosynthetic reductases, suggesting evolutionary relationship between acylsugar and capsaicinoid biosynthesis. Analysis of publicly available datasets revealed enriched expression of Sopen05g009610 orthologs in trichomes of several acylsugar-producing species. Similarly, orthologs of Sopen07g006810 were identified as solanaceous trichome-enriched members, which form a phylogenetic clade distinct from those of mesophyll-expressed "regular" Rubisco small subunits. Furthermore, δ13C analyses indicated recycling of metabolic CO2 into acylsugars by Sopen07g006810 and showed how trichomes support high levels of specialized metabolite production. These findings have implications for genetic manipulation of trichome-specialized metabolism in solanaceous crops.

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.

Protocorm-like-body extract of Phalaenopsis aphrodite combats watermelon fruit blotch disease

Bacterial fruit blotch, caused by the seedborne gram-negative bacterium Acidovorax citrulli, is one of the most destructive bacterial diseases of cucurbits (gourds) worldwide. Despite its prevalence, effective and reliable means to control bacterial fruit blotch remain limited. Transcriptomic analyses of tissue culture-based regeneration processes have revealed that organogenesis-associated cellular reprogramming is often associated with upregulation of stress- and defense-responsive genes. Yet, there is limited evidence supporting the notion that the reprogrammed cellular metabolism of the regenerated tissued confers bona fide antimicrobial activity. Here, we explored the anti-bacterial activity of protocorm-like-bodies (PLBs) of Phalaenopsis aphrodite. Encouragingly, we found that the PLB extract was potent in slowing growth of A. citrulli, reducing the number of bacteria attached to watermelon seeds, and alleviating disease symptoms of watermelon seedlings caused by A. citrulli. Because the anti-bacterial activity can be fractionated chemically, we predict that reprogrammed cellular activity during the PLB regeneration process produces metabolites with antibacterial activity. In conclusion, our data demonstrated the antibacterial activity in developing PLBs and revealed the potential of using orchid PLBs to discover chemicals to control bacterial fruit blotch disease.

Polyketide reductases in defense-related parasorboside biosynthesis in Gerbera hybrida share processing strategies with microbial polyketide synthase systems

Plant polyketides are well-known for their crucial functions in plants and their importance in the context of human health. They are synthesized by type III polyketide synthases (PKSs) and their final functional diversity is determined by post-PKS tailoring enzymes. Gerbera hybrida is rich in two defense-related polyketides: gerberin and parasorboside. Their synthesis is known to be initiated by GERBERA 2-PYRONE SYNTHASE 1 (G2PS1), but the polyketide reductases (PKRs) that determine their final structure have not yet been identified. We identified two PKR candidates in the pathway, GERBERA REDUCTASE 1 (GRED1) and GRED2. Gene expression and metabolite analysis of different gerbera tissues, cultivars, and transgenic gerbera plants, and in vitro enzyme assays, were performed for functional characterization of the enzymes. GRED1 and GRED2 catalyze the second reduction step in parasorboside biosynthesis. They reduce the proximal keto domain of the linear CoA bound intermediate before lactonization. We identified a crucial tailoring step in an important gerbera PKS pathway and show that plant polyketide biosynthesis shares processing strategies with fungi and bacteria. The two tailoring enzymes are recruited from the ancient sporopollenin biosynthetic pathway to a defense-related PKS pathway in gerbera. Our data provide an example of how plants recruit conserved genes to new functions in secondary metabolism that are important for environmental adaptation.

Consecutive action of two BAHD acyltransferases promotes tetracoumaroyl spermine accumulation in chicory

Fully substituted phenolamide accumulation in the pollen coat of Eudicotyledons is a conserved evolutionary chemical trait. Interestingly, spermidine derivatives are replaced by spermine derivatives as the main phenolamide accumulated in the Asteraceae family. Here, we show that the full substitution of spermine in chicory (Cichorium intybus) requires the successive action of two enzymes, that is spermidine hydroxycinnamoyl transferase-like proteins 1 and 2 (CiSHT1 and CiSHT2), two members of the BAHD enzyme family. Deletion of these genes in chicory using CRISPR/Cas9 gene editing technology evidenced that CiSHT2 catalyzes the first N-acylation steps, whereas CiSHT1 fulfills the substitution to give rise to tetracoumaroyl spermine. Additional experiments using Nicotiana benthamiana confirmed these findings. Expression of CiSHT2 alone promoted partially substituted spermine accumulation, and coexpression of CiSHT2 and CiSHT1 promoted synthesis and accumulation of the fully substituted spermine. Structural characterization of the main product of CiSHT2 using nuclear magnetic resonance revealed that CiSHT2 preferentially catalyzed N-acylation of secondary amines to form N5,N10-dicoumaroyl spermine, whereas CiSHT1 used this substrate to synthesize tetracoumaroyl spermine. We showed that spermine availability may be a key determinant toward preferential accumulation of spermine derivatives over spermidine derivatives in chicory. Our results reveal a subfunctionalization among the spermidine hydroxycinnamoyl transferase that was accompanied by a modification of free polyamine metabolism that has resulted in the accumulation of this new phenolamide in chicory and most probably in all Asteraceae. Finally, genetically engineered yeast (Saccharomyces cerevisiae) was shown to be a promising host platform to produce these compounds.

Using interdisciplinary, phylogeny-guided approaches to understand the evolution of plant metabolism

To cope with relentless environmental pressures, plants produce an arsenal of structurally diverse chemicals, often called specialized metabolites. These lineage-specific compounds are derived from the simple building blocks made by ubiquitous core metabolic pathways. Although the structures of many specialized metabolites are known, the underlying metabolic pathways and the evolutionary events that have shaped the plant chemical diversity landscape are only beginning to be understood. However, with the advent of multi-omics data sets and the relative ease of studying pathways in previously intractable non-model species, plant specialized metabolic pathways are now being systematically identified. These large datasets also provide a foundation for comparative, phylogeny-guided studies of plant metabolism. Comparisons of metabolic traits and features like chemical abundances, enzyme activities, or gene sequences from phylogenetically diverse plants provide insights into how metabolic pathways evolved. This review highlights the power of studying evolution through the lens of comparative biochemistry, particularly how placing metabolism into a phylogenetic context can help a researcher identify the metabolic innovations enabling the evolution of structurally diverse plant metabolites.

Predictive metabolomics of multiple Atacama plant species unveils a core set of generic metabolites for extreme climate resilience

Current crop yield of the best ideotypes is stagnating and threatened by climate change. In this scenario, understanding wild plant adaptations in extreme ecosystems offers an opportunity to learn about new mechanisms for resilience. Previous studies have shown species specificity for metabolites involved in plant adaptation to harsh environments. Here, we combined multispecies ecological metabolomics and machine learning-based generalized linear model predictions to link the metabolome to the plant environment in a set of 24 species belonging to 14 families growing along an altitudinal gradient in the Atacama Desert. Thirty-nine common compounds predicted the plant environment with 79% accuracy, thus establishing the plant metabolome as an excellent integrative predictor of environmental fluctuations. These metabolites were independent of the species and validated both statistically and biologically using an independent dataset from a different sampling year. Thereafter, using multiblock predictive regressions, metabolites were linked to climatic and edaphic stressors such as freezing temperature, water deficit and high solar irradiance. These findings indicate that plants from different evolutionary trajectories use a generic metabolic toolkit to face extreme environments. These core metabolites, also present in agronomic species, provide a unique metabolic goldmine for improving crop performances under abiotic pressure.

Generalizable approaches for genomic prediction of metabolites in plants

Plant metabolites are important traits for plant breeders seeking to improve nutrition and agronomic performance yet integrating selection for metabolomic traits can be limited by phenotyping expense and degree of genetic characterization, especially of uncommon metabolites. As such, developing generalizable genomic selection methods based on biochemical pathway biology for metabolites that are transferable across plant populations would benefit plant breeding programs. We tested genomic prediction accuracy for >600 metabolites measured by gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) in oat (Avena sativa L.) seed. Using a discovery germplasm panel, we conducted metabolite genome-wide association study (mGWAS) and selected loci to use in multikernel models that encompassed metabolome-wide mGWAS results or mGWAS from specific metabolite structures or biosynthetic pathways. Metabolite kernels developed from LC-MS metabolites in the discovery panel improved prediction accuracy of LC-MS metabolite traits in the validation panel consisting of more advanced breeding lines. No approach, however, improved prediction accuracy for GC-MS metabolites. We ranked model performance by metabolite and found that metabolites with similar polarity had consistent rankings of models. Overall, testing biological rationales for developing kernels for genomic prediction across populations contributes to developing frameworks for plant breeding for metabolite traits.

Origin, Expansion, and Divergence of ETHYLENE-INSENSITIVE 3 (EIN3)/EIN3-LIKE Transcription Factors During Streptophytes Evolution

The transition of plants to land required several regulatory adaptive mechanisms. Little is known about these mechanisms, but they no doubt involved the evolution of transcription factor (TF) families. ETHYLENE-INSENSITIVE 3 (EIN3)/EIN3-LIKE (EIL) transcription factors (TFs) are core components of the ethylene signaling pathway that play important roles in almost every aspect of plant development and environmental responses by regulating the transcription of numerous genes. However, the evolutionary history of EIN3/EIL TFs, which are present in a wide range of streptophytes, is still not clear. Here, to explore the evolution and functions of EIN3/EIL TFs, we performed phylogenetic analysis of these TFs and investigated their gene and protein structures as well as sequence features. Our results suggest that the EIN3/EIL TF family was already was already present in the ancestor of streptophyte algae. Phylogenetic analysis divided the EIN3/EIL TFs into three groups (Group A-C). Analysis of gene and protein structure revealed that most of the structural features of these TFs had already formed in ancient lineages. Further investigation suggested that all groups have undergone several duplication events related to whole-genome duplications in plants, generating multiple, functionally diverse gene copies. Therefore, as plants colonized terrestrial habitats and evolved key traits, the EIN3/EIL TF family expanded broadly via multiple duplication events, which could have promoted their fundamental neo- and sub-functionalization to help plants adapt to terrestrial life. Our findings shed light on the functional evolution of the EIN3/EIL TF family in the streptophytes.

Interplay between gene expression and gene architecture as a consequence of gene and genome duplications: evidence from metabolic genes of Arabidopsis thaliana

Gene and genome duplications have been widespread during the evolution of flowering plant which resulted in the increment of biological complexity as well as creation of plasticity of a genome helping the species to adapt to changing environments. Duplicated genes with higher evolutionary rates can act as a mechanism of generating novel functions in secondary metabolism. In this study, we explored duplication as a potential factor governing the expression heterogeneity and gene architecture of Primary Metabolic Genes (PMGs) and Secondary Metabolic Genes (SMGs) of Arabidopsis thaliana. It is remarkable that different types of duplication processes controlled gene expression and tissue specificity differently in PMGs and SMGs. A complex relationship exists between gene architecture and expression patterns of primary and secondary metabolic genes. Our study reflects, expression heterogeneity and gene structure variation of primary and secondary metabolism in Arabidopsis thaliana are partly results of duplication events of different origins. Our study suggests that duplication has differential effect on PMGs and SMGs regarding expression pattern by controlling gene structure, epigenetic modifications, multifunctionality and subcellular compartmentalization. This study provides an insight into the evolution of metabolism in plants in the light of gene and genome scale duplication.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01188-2.

Computational prediction of plant metabolic pathways

Uncovering genes encoding enzymes responsible for the biosynthesis of diverse plant metabolites is essential for metabolic engineering and production of plant metabolite-derived medicine. With the availability of multi-omics data for an ever-increasing number of plant species and the development of computational approaches, the metabolic pathways of many important plant compounds can be predicted, complementing a more traditional genetic and/or biochemical approach. Here, we summarize recent progress in predicting plant metabolic pathways using genome, transcriptome, proteome, interactome, and/or metabolome data, and the utility of integrating these data with machine learning to further improve metabolic pathway predictions.

The ease and complexity of identifying and using specialized metabolites for crop engineering

Plants produce a broad variety of specialized metabolites with distinct biological activities and potential applications. Despite this potential, most biosynthetic pathways governing specialized metabolite production remain largely unresolved across the plant kingdom. The rapid advancement of genetics and biochemical tools has enhanced our ability to identify plant specialized metabolic pathways. Further advancements in transgenic technology and synthetic biology approaches have extended this to a desire to design new pathways or move existing pathways into new systems to address long-running difficulties in crop systems. This includes improving abiotic and biotic stress resistance, boosting nutritional content, etc. In this review, we assess the potential and limitations for (1) identifying specialized metabolic pathways in plants with multi-omics tools and (2) using these enzymes in synthetic biology or crop engineering. The goal of these topics is to highlight areas of research that may need further investment to enhance the successful application of synthetic biology for exploiting the myriad of specialized metabolic pathways.

Reciprocal mutations of two multifunctional β-amyrin synthases from Barbarea vulgaris shift α/β-amyrin ratios

In the wild cruciferous wintercress (Barbarea vulgaris), β-amyrin-derived saponins are involved in resistance against insect herbivores like the major agricultural pest diamondback moth (Plutella xylostella). Enzymes belonging to the 2,3-oxidosqualene cyclase family have been identified and characterized in B. vulgaris G-type and P-type plants that differ in their natural habitat, insect resistance and saponin content. Both G-type and P-type plants possess highly similar 2,3-oxidosqualene cyclase enzymes that mainly produce β-amyrin (Barbarea vulgaris Lupeol synthase 5 G-Type; BvLUP5-G) or α-amyrin (Barbarea vulgaris Lupeol synthase 5 P-Type; BvLUP5-P), respectively. Despite the difference in product formation, the two BvLUP5 enzymes are 98% identical at the amino acid level. This provides a unique opportunity to investigate determinants of product formation, using the B. vulgaris 2,3-oxidosqualene cyclase enzymes as a model for studying amino acid residues that determine differences in product formation. In this study, we identified two amino acid residues at position 121 and 735 that are responsible for the dominant changes in generated product ratios of β-amyrin and α-amyrin in both BvLUP5 enzymes. These amino acid residues have not previously been highlighted as directly involved in 2,3-oxidosqualene cyclase product specificity. Our results highlight the functional diversity and promiscuity of 2,3-oxidosqualene cyclase enzymes. These enzymes serve as important mediators of metabolic plasticity throughout plant evolution.

Knockout of Arabidopsis thaliana VEP1, Encoding a PRISE (Progesterone 5β-Reductase/Iridoid Synthase-Like Enzyme), Leads to Metabolic Changes in Response to Exogenous Methyl Vinyl Ketone (MVK)

Small or specialized natural products (SNAPs) produced by plants vary greatly in structure and function, leading to selective advantages during evolution. With a limited number of genes available, a high promiscuity of the enzymes involved allows the generation of a broad range of SNAPs in complex metabolic networks. Comparative metabolic studies may help to understand why—or why not—certain SNAPs are produced in plants. Here, we used the wound-induced, vein patterning regulating VEP1 (AtStR1, At4g24220) and its paralogue gene on locus At5g58750 (AtStR2) from Arabidopsis to study this issue. The enzymes encoded by VEP1-like genes were clustered under the term PRISEs (progesterone 5β-reductase/iridoid synthase-like enzymes) as it was previously demonstrated that they are involved in cardenolide and/or iridoid biosynthesis in other plants. In order to further understand the general role of PRISEs and to detect additional more “accidental” roles we herein characterized A. thaliana steroid reductase 1 (AtStR1) and compared it to A. thaliana steroid reductase 2 (AtStR2). We used A. thaliana Col-0 wildtype plants as well as VEP1 knockout mutants and VEP1 knockout mutants overexpressing either AtStR1 or AtStR2 to investigate the effects on vein patterning and on the stress response after treatment with methyl vinyl ketone (MVK). Our results added evidence to the assumption that AtStR1 and AtStR2, as well as PRISEs in general, play specific roles in stress and defense situations and may be responsible for sudden metabolic shifts.

It happened again: Convergent evolution of acylglucose specialized metabolism in black nightshade and wild tomato

Plants synthesize myriad phylogenetically restricted specialized (aka “secondary”) metabolites with diverse structures. Metabolism of acylated sugar esters in epidermal glandular secreting trichomes across the Solanaceae (nightshade) family is ideal for investigating the mechanisms of evolutionary metabolic diversification. We developed methods to structurally analyze acylhexose mixtures by 2D NMR, which led to the insight that the Old World species black nightshade (Solanum nigrum) accumulates acylglucoses and acylinositols in the same tissue. Detailed in vitro biochemistry, cross-validated by in vivo virus-induced gene silencing, revealed two unique features of the four-step acylglucose biosynthetic pathway: A trichome-expressed, neofunctionalized invertase-like enzyme, SnASFF1, converts BAHD-produced acylsucroses to acylglucoses, which, in turn, are substrates for the acylglucose acyltransferase, SnAGAT1. This biosynthetic pathway evolved independently from that recently described in the wild tomato Solanum pennellii, reinforcing that acylsugar biosynthesis is evolutionarily dynamic with independent examples of primary metabolic enzyme cooption and additional variation in BAHD acyltransferases.

Comparisons within the Rice GA 2-Oxidase Gene Family Revealed Three Dominant Paralogs and a Functional Attenuated Gene that Led to the Identification of Four Amino Acid Variants Associated with GA Deactivation Capability

BACKGROUND

GA 2-oxidases (GA2oxs) are involved in regulating GA homeostasis in plants by inactivating bioactive GAs through 2β-hydroxylation. Rice GA2oxs are encoded by a family of 10 genes; some of them have been characterized, but no comprehensive comparisons for all these genes have been conducted.

RESULTS

Rice plants with nine functional GA2oxs were demonstrated in the present study, and these genes not only were differentially expressed but also revealed various capabilities for GA deactivation based on their height-reducing effects in transgenic plants. Compared to that of wild-type plants, the relative plant height (RPH) of transgenic plants was scored to estimate their reducing effects, and 8.3% to 59.5% RPH was observed. Phylogenetic analysis of class I GA2ox genes revealed two functionally distinct clades in the Poaceae. The OsGA2ox3, 4, and 8 genes belonging to clade A showed the most severe effect (8.3% to 8.7% RPH) on plant height reduction, whereas the OsGA2ox7 gene belonging to clade B showed the least severe effect (59.5% RPH). The clade A OsGA2ox3 gene contained two conserved C186/C194 amino acids that were crucial for enzymatic activity. In the present study, these amino acids were replaced with OsGA2ox7-conserved arginine (C186R) and proline (C194P), respectively, or simultaneously (C186R/C194P) to demonstrate their importance in planta. Another two amino acids, Q220 and Y274, conserved in OsGA2ox3 were substituted with glutamic acid (E) and phenylalanine (F), respectively, or simultaneously to show their significance in planta. In addition, through sequence divergence, RNA expression profile and GA deactivation capability analyses, we proposed that OsGA2ox1, OsGA2ox3 and OsGA2ox6 function as the predominant paralogs in each of their respective classes.

CONCLUSIONS

This study demonstrates rice has nine functional GA2oxs and the class I GA2ox genes are divided into two functionally distinct clades. Among them, the OsGA2ox7 of clade B is a functional attenuated gene and the OsGA2ox1, OsGA2ox3 and OsGA2ox6 are the three predominant paralogs in the family.

Evolutionary Implications of a Peroxidase with High Affinity for Cinnamyl Alcohols from Physcomitrium patens, a Non-Vascular Plant

Physcomitrium (Physcomitrella) patens is a bryophyte highly tolerant to different stresses, allowing survival when water supply is a limiting factor. This moss lacks a true vascular system, but it has evolved a primitive water-conducting system that contains lignin-like polyphenols. By means of a three-step protocol, including ammonium sulfate precipitation, adsorption chromatography on phenyl Sepharose and cationic exchange chromatography on SP Sepharose, we were able to purify and further characterize a novel class III peroxidase, PpaPrx19, upregulated upon salt and H₂O₂ treatments. This peroxidase, of a strongly basic nature, shows surprising homology to angiosperm peroxidases related to lignification, despite the lack of true lignins in P. patens cell walls. Moreover, PpaPrx19 shows catalytic and kinetic properties typical of angiosperm peroxidases involved in oxidation of monolignols, being able to efficiently use hydroxycinnamyl alcohols as substrates. Our results pinpoint the presence in P. patens of peroxidases that fulfill the requirements to be involved in the last step of lignin biosynthesis, predating the appearance of true lignin.

Metabolic networks of the Nicotiana genus in the spotlight: content, progress and outlook

Manually curated metabolic databases residing at the Sol Genomics Network comprise two taxon-specific databases for the Solanaceae family, i.e. SolanaCyc and the genus Nicotiana, i.e. NicotianaCyc as well as six species-specific databases for Nicotiana tabacum TN90, N. tabacum K326, Nicotiana benthamiana, N. sylvestris, N. tomentosiformis and N. attenuata. New pathways were created through the extraction, examination and verification of related data from the literature and the aid of external database guided by an expert-led curation process. Here we describe the curation progress that has been achieved in these databases since the first release version 1.0 in 2016, the curation flow and the curation process using the example metabolic pathway for cholesterol in plants. The current content of our databases comprises 266 pathways and 36 superpathways in SolanaCyc and 143 pathways plus 21 superpathways in NicotianaCyc, manually curated and validated specifically for the Solanaceae family and Nicotiana genus, respectively. The curated data have been propagated to the respective Nicotiana-specific databases, which resulted in the enrichment and more accurate presentation of their metabolic networks. The quality and coverage in those databases have been compared with related external databases and discussed in terms of literature support and metabolic content.

Genetic variation, environment and demography intersect to shape Arabidopsis defense metabolite variation across Europe

Plants produce diverse metabolites to cope with the challenges presented by complex and ever-changing environments. These challenges drive the diversification of specialized metabolites within and between plant species. However, we are just beginning to understand how frequently new alleles arise controlling specialized metabolite diversity and how the geographic distribution of these alleles may be structured by ecological and demographic pressures. Here, we measure the variation in specialized metabolites across a population of 797 natural Arabidopsis thaliana accessions. We show that a combination of geography, environmental parameters, demography and different genetic processes all combine to influence the specific chemotypes and their distribution. This showed that causal loci in specialized metabolism contain frequent independently generated alleles with patterns suggesting potential within-species convergence. This provides a new perspective about the complexity of the selective forces and mechanisms that shape the generation and distribution of allelic variation that may influence local adaptation.

Since plants cannot move, they have evolved chemical defenses to help them respond to changes in their surroundings. For example, where animals run from predators, plants may produce toxins to put predators off. This approach is why plants are such a rich source of drugs, poisons, dyes and other useful substances. The chemicals plants produce are known as specialized metabolites, and they can change a lot between, and even within, plant species. The variety of specialized metabolites is a result of genetic changes and evolution over millions of years.

Evolution is a slow process, yet plants are able to rapidly develop new specialized metabolites to protect them from new threats. Even different populations of the same species produce many distinct metabolites that help them survive in their surroundings. However, the factors that lead plants to produce new metabolites are not well understood, and it is not known how this affects genetic variation.

To gain a better understanding of this process, Katz et al. studied 797 European variants of a common weed species called Arabidopsis thaliana, which is widely studied. The investigation found that many factors affect the range of specialized metabolites in each variant. These included local geography and environment, as well as genetics and population history (demography). Katz et al. revealed a pattern of relationships between the variants that could mirror their evolutionary history as the species spread and adapted to new locations.

These results highlight the complex network of factors that affect plant evolution. Rapid diversification is key to plant survival in new and changing environments and has resulted in a wide range of specialized metabolites. As such they are of interest both for studying plant evolution and for understanding their ecology. Expanding similar work to more populations and other species will broaden the scope of our ability to understand how plants adapt to their surroundings.

Kingdom-wide analysis of the evolution of the plant type III polyketide synthase superfamily

The emergence of type III polyketide synthases (PKSs) was a prerequisite for the conquest of land by the green lineage. Within the PKS superfamily, chalcone synthases (CHSs) provide the entry point reaction to the flavonoid pathway, while LESS ADHESIVE POLLEN 5 and 6 (LAP5/6) provide constituents of the outer exine pollen wall. To study the deep evolutionary history of this key family, we conducted phylogenomic synteny network and phylogenetic analyses of whole-genome data from 126 species spanning the green lineage including Arabidopsis thaliana, tomato (Solanum lycopersicum), and maize (Zea mays). This study thereby combined study of genomic location and context with changes in gene sequences. We found that the two major clades, CHS and LAP5/6 homologs, evolved early by a segmental duplication event prior to the divergence of Bryophytes and Tracheophytes. We propose that the macroevolution of the type III PKS superfamily is governed by whole-genome duplications and triplications. The combined phylogenetic and synteny analyses in this study provide insights into changes in the genomic location and context that are retained for a longer time scale with more recent functional divergence captured by gene sequence alterations.

Strategies to Modulate Specialized Metabolism in Mediterranean Crops: From Molecular Aspects to Field

Plant specialized metabolites (SMs) play an important role in the interaction with the environment and are part of the plant defense response. These natural products are volatile, semi-volatile and non-volatile compounds produced from common building blocks deriving from primary metabolic pathways and rapidly evolved to allow a better adaptation of plants to environmental cues. Specialized metabolites include terpenes, flavonoids, alkaloids, glucosinolates, tannins, resins, etc. that can be used as phytochemicals, food additives, flavoring agents and pharmaceutical compounds. This review will be focused on Mediterranean crop plants as a source of SMs, with a special attention on the strategies that can be used to modulate their production, including abiotic stresses, interaction with beneficial soil microorganisms and novel genetic approaches.

Integration of mass spectral fingerprinting analysis with precursor ion (MS1) quantification for the characterisation of botanical extracts: application to extracts of Centella asiatica (L.) Urban

INTRODUCTION

The phytochemical composition of plant material governs the bioactivity and potential health benefits as well as the outcomes and reproducibility of laboratory studies and clinical trials.

OBJECTIVE

The objective of this work was to develop an efficient method for the in-depth characterisation of plant extracts and quantification of marker compounds that can be potentially used for subsequent product integrity studies. Centella asiatica (L.) Urb., an Ayurvedic herb with potential applications in enhancing mental health and cognitive function, was used as a case study.

METHODS

A quadrupole time-of-flight analyser in conjunction with an optimised high-performance liquid chromatography (HPLC) separation was used for in-depth untargeted fingerprinting and post-acquisition precursor ion quantification to determine levels of distinct phytochemicals in various C. asiatica extracts.

RESULTS

We demonstrate the utility of this workflow for the characterisation of extracts of C. asiatica. This integrated workflow allowed the identification or tentative identification of 117 compounds, chemically interconnected based on Tanimoto chemical similarity, and the accurate quantification of 24 phytochemicals commonly found in C. asiatica extracts.

CONCLUSION

We report a phytochemical analysis method combining liquid chromatography, high resolution mass spectral data acquisition, and post-acquisition interrogation that allows chemical fingerprints of botanicals to be obtained in conjunction with accurate quantification of distinct phytochemicals. The variability in the composition of specialised metabolites across different C. asiatica accessions was substantial, demonstrating that detailed characterisation of plant extracts is a prerequisite for reproducible use in laboratory studies, clinical trials and safe consumption. The methodological approach is generally applicable to other botanical products.

Critical enzymes for biosynthesis of cucurbitacin derivatives in watermelon and their biological significance

Various cucurbitacins have been isolated, and their structures have been elucidated. Owing to their economic potential and importance as active pharmacological compounds, their cytotoxicity in various cancer cells has been assessed. Here, we mined several candidate genes with potential involvement in cucurbitacin biosynthesis in watermelon (Citrullus lanatus) and performed in vitro enzymatic assays and instrumental analyses using various substrates to identify cucurbitacin functions and products. Enzymatic activities of two acetyltransferases (ACTs) and one UDP-glucosyltransferase (UGT) against cucurbitacins were confirmed, resulting in the synthesis of novel cucurbitacins in vivo and/or in vitro to our knowledge. As ACTs and UGT are involved in the dynamic conversion of cucurbitacins by catalyzing acetylation and glucosylation at moieties in the cucurbitacins skeleton, these findings improve our knowledge on how these genes contribute to the diversity of cucurbitacins.

Kim et al. use RNAseq of two watermelons to select candidate genes coding for enzymes that catalyze modifications of cucurbitacins. They characterise four of the 16 candidate enzymes (3 different acetyltransferases and one UDP-glucosyltransferase) by HPLC, LC-MS, NMR, and in vitro enzymatic assay. They further show with in vivo assay in Drosophila, that acetylation of cucurbitacin increases neuronal activity in insects.

Drivers of metabolic diversification: how dynamic genomic neighbourhoods generate new biosynthetic pathways in the Brassicaceae

Plants produce an array of specialized metabolites with important ecological functions. The mechanisms underpinning the evolution of new biosynthetic pathways are not well-understood. Here, we exploit available genome sequence resources to investigate triterpene biosynthesis across the Brassicaceae. Oxidosqualene cyclases (OSCs) catalyze the first committed step in triterpene biosynthesis. Systematic analysis of 13 sequenced Brassicaceae genomes was performed to identify all OSC genes. The genome neighbourhoods (GNs) around a total of 163 OSC genes were investigated to identify Pfam domains significantly enriched in these regions. All-vs-all comparisons of OSC neighbourhoods and phylogenomic analysis were used to investigate the sequence similarity and evolutionary relationships of the numerous candidate triterpene biosynthetic gene clusters (BGCs) observed. Functional analysis of three representative BGCs was carried out and their triterpene pathway products were elucidated. Our results indicate that plant genomes are remarkably plastic, and that dynamic GNs generate new biosynthetic pathways in different Brassicaceae lineages by shuffling the genes encoding a core palette of triterpene-diversifying enzymes, presumably in response to strong environmental selection pressure. These results illuminate a genomic basis for diversification of plant-specialized metabolism through natural combinatorics of enzyme families, which can be mimicked using synthetic biology to engineer diverse bioactive molecules.

Evolution of a plant gene cluster in Solanaceae and emergence of metabolic diversity

Plants produce phylogenetically and spatially restricted, as well as structurally diverse specialized metabolites via multistep metabolic pathways. Hallmarks of specialized metabolic evolution include enzymatic promiscuity and recruitment of primary metabolic enzymes and examples of genomic clustering of pathway genes. Solanaceae glandular trichomes produce defensive acylsugars, with sidechains that vary in length across the family. We describe a tomato gene cluster on chromosome 7 involved in medium chain acylsugar accumulation due to trichome specific acyl-CoA synthetase and enoyl-CoA hydratase genes. This cluster co-localizes with a tomato steroidal alkaloid gene cluster and is syntenic to a chromosome 12 region containing another acylsugar pathway gene. We reconstructed the evolutionary events leading to this gene cluster and found that its phylogenetic distribution correlates with medium chain acylsugar accumulation across the Solanaceae. This work reveals insights into the dynamics behind gene cluster evolution and cell-type specific metabolite diversity.

Structural Studies of Glutamate Dehydrogenase (Isoform 1) From Arabidopsis thaliana, an Important Enzyme at the Branch-Point Between Carbon and Nitrogen Metabolism

Glutamate dehydrogenase (GDH) releases ammonia in a reversible NAD(P)⁺-dependent oxidative deamination of glutamate that yields 2-oxoglutarate (2OG). In current perception, GDH contributes to Glu homeostasis and plays a significant role at the junction of carbon and nitrogen assimilation pathways. GDHs are members of a superfamily of ELFV (Glu/Leu/Phe/Val) amino acid dehydrogenases and are subdivided into three subclasses, based on coenzyme specificity: NAD⁺-specific, NAD⁺/NADP⁺ dual-specific, and NADP⁺-specific. We determined in this work that the mitochondrial AtGDH1 isozyme from A. thaliana is NAD⁺-specific. Altogether, A. thaliana expresses three GDH isozymes (AtGDH1-3) targeted to mitochondria, of which AtGDH2 has an extra EF-hand motif and is stimulated by calcium. Our enzymatic assays of AtGDH1 established that its sensitivity to calcium is negligible. In vivo the AtGDH1-3 enzymes form homo- and heterohexamers of varied composition. We solved the crystal structure of recombinant AtGDH1 in the apo-form and in complex with NAD⁺ at 2.59 and 2.03 Å resolution, respectively. We demonstrate also that both in the apo form and in 1:1 complex with NAD⁺, it forms D ₃-symmetric homohexamers. A subunit of AtGDH1 consists of domain I, which is involved in hexamer formation and substrate binding, and of domain II which binds coenzyme. Most of the subunits in our crystal structures, including those in NAD⁺ complex, are in open conformation, with domain II forming a large (albeit variable) angle with domain I. One of the subunits of the AtGDH1-NAD⁺ hexamer contains a serendipitous 2OG molecule in the active site, causing a dramatic (∼25°) closure of the domains. We provide convincing evidence that the N-terminal peptide preceding domain I is a mitochondrial targeting signal, with a predicted cleavage site for mitochondrial processing peptidase (MPP) at Leu17-Leu18 that is followed by an unexpected potassium coordination site (Ser27, Ile30). We also identified several MPD [(+/-)-2-methyl-2,4-pentanediol] binding sites with conserved sequence. Although AtGDH1 is insensitive to MPD in our assays, the observation of druggable sites opens a potential for non-competitive herbicide design.

Evolutionary significance of amino acid permease transporters in 17 plants from Chlorophyta to Angiospermae

Background

Nitrogen is an indispensable nutrient for plant growth. It is used and transported in the form of amino acids in living organisms. Transporting amino acids to various parts of plants requires relevant transport proteins, such as amino acid permeases (AAPs), which were our focus in this study.

Results

We found that 5 AAP genes were present in Chlorophyte species and more AAP genes were predicted in Bryophyta and Lycophytes. Two main groups were defined and group I comprised 5 clades. Our phylogenetic analysis indicated that the origin of clades 2, 3, and 4 is Gymnospermae and that these clades are closely related. The members of clade 1 included Chlorophyta to Gymnospermae. Group II, as a new branch consisting of non-seed plants, is first proposed in our research. Our results also indicated that the AAP family was already present in Chlorophyta and then expanded accompanying the development of vasculature. Concurrently, the AAP family experienced multiple duplication events that promoted the generation of new functions and differentiation of sub-functions.

Conclusions

Our findings suggest that the AAP gene originated in Chlorophyta, and some non-seed AAP genes clustered in one group. A second group, which contained plants of all evolutionary stages, indicated the evolution of AAPs. These new findings can be used to guide future research.

The evolutionary origins of the cat attractant nepetalactone in catnip

Catnip or catmint (Nepeta spp.) is a flowering plant in the mint family (Lamiaceae) famed for its ability to attract cats. This phenomenon is caused by the compound nepetalactone, a volatile iridoid that also repels insects. Iridoids are present in many Lamiaceae species but were lost in the ancestor of the Nepetoideae, the subfamily containing Nepeta. Using comparative genomics, ancestral sequence reconstructions, and phylogenetic analyses, we probed the re-emergence of iridoid biosynthesis in Nepeta. The results of these investigations revealed mechanisms for the loss and subsequent re-evolution of iridoid biosynthesis in the Nepeta lineage. We present evidence for a chronology of events that led to the formation of nepetalactone biosynthesis and its metabolic gene cluster. This study provides insights into the interplay between enzyme and genome evolution in the origins, loss, and re-emergence of plant chemical diversity.

Enhancing flavonoid production by promiscuous activity of prenyltransferase, BrPT2 from Boesenbergia rotunda

Flavonoids and prenylated flavonoids are active components in medicinal plant extracts which exhibit beneficial effects on human health. Prenylated flavonoids consist of a flavonoid core with a prenyl group attached to it. This prenylation process is catalyzed by prenyltranferases (PTs). At present, only a few flavonoid-related PT genes have been identified. In this study, we aimed to investigate the roles of PT in flavonoid production. We isolated a putative PT gene (designated as BrPT2) from a medicinal ginger, Boesenbergia rotunda. The deduced protein sequence shared highest gene sequence homology (81%) with the predicted homogentisate phytyltransferase 2 chloroplastic isoform X1 from Musa acuminata subsp. Malaccensis. We then cloned the BrPT2 into pRI vector and expressed in B. rotunda cell suspension cultures via Agrobacterium-mediated transformation. The BrPT2-expressing cells were fed with substrate, pinostrobin chalcone, and their products were analyzed by liquid chromatography mass spectrometry. We found that the amount of flavonoids, namely alpinetin, pinostrobin, naringenin and pinocembrin, in BrPT2-expressing cells was higher than those obtained from the wild type cells. However, we were unable to detect any targeted prenylated flavonoids. Further in-vitro assay revealed that the reaction containing the BrPT2 protein produced the highest accumulation of pinostrobin from the substrate pinostrobin chalcone compared to the reaction without BrPT2 protein, suggesting that BrPT2 was able to accelerate the enzymatic reaction. The finding of this study implied that the isolated BrPT2 may not be involved in the prenylation of pinostrobin chalcone but resulted in high yield and production of other flavonoids, which is likely related to enzyme promiscuous activities.

Role of cytosolic, tyrosine-insensitive prephenate dehydrogenase in Medicago truncatula

l-Tyrosine (Tyr) is an aromatic amino acid synthesized de novo in plants and microbes downstream of the shikimate pathway. In plants, Tyr and a Tyr pathway intermediate, 4-hydroxyphenylpyruvate (HPP), are precursors to numerous specialized metabolites, which are crucial for plant and human health. Tyr is synthesized in the plastids by a TyrA family enzyme, arogenate dehydrogenase (ADH/TyrAa), which is feedback inhibited by Tyr. Additionally, many legumes possess prephenate dehydrogenases (PDH/TyrAp), which are insensitive to Tyr and localized to the cytosol. Yet the role of PDH enzymes in legumes is currently unknown. This study isolated and characterized Tnt1-transposon mutants of MtPDH1 (pdh1) in Medicago truncatula to investigate PDH function. The pdh1 mutants lacked PDH transcript and PDH activity, and displayed little aberrant morphological phenotypes under standard growth conditions, providing genetic evidence that MtPDH1 is responsible for the PDH activity detected in M. truncatula. Though plant PDH enzymes and activity have been specifically found in legumes, nodule number and nitrogenase activity of pdh1 mutants were not significantly reduced compared with wild-type (Wt) during symbiosis with nitrogen-fixing bacteria. Although Tyr levels were not significantly different between Wt and mutants under standard conditions, when carbon flux was increased by shikimate precursor feeding, mutants accumulated significantly less Tyr than Wt. These data suggest that MtPDH1 is involved in Tyr biosynthesis when the shikimate pathway is stimulated and possibly linked to unidentified legume-specific specialized metabolism.

The evolution of metabolism: How to test evolutionary hypotheses at the genomic level

The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.

A conserved strategy of chalcone isomerase-like protein to rectify promiscuous chalcone synthase specificity

Land plants produce diverse flavonoids for growth, survival, and reproduction. Chalcone synthase is the first committed enzyme of the flavonoid biosynthetic pathway and catalyzes the production of 2',4,4',6'-tetrahydroxychalcone (THC). However, it also produces other polyketides, including p-coumaroyltriacetic acid lactone (CTAL), because of the derailment of the chalcone-producing pathway. This promiscuity of CHS catalysis adversely affects the efficiency of flavonoid biosynthesis, although it is also believed to have led to the evolution of stilbene synthase and p-coumaroyltriacetic acid synthase. In this study, we establish that chalcone isomerase-like proteins (CHILs), which are encoded by genes that are ubiquitous in land plant genomes, bind to CHS to enhance THC production and decrease CTAL formation, thereby rectifying the promiscuous CHS catalysis. This CHIL function has been confirmed in diverse land plant species, and represents a conserved strategy facilitating the efficient influx of substrates from the phenylpropanoid pathway to the flavonoid pathway.

S-adenosylmethionine synthases in plants: Structural characterization of type I and II isoenzymes from Arabidopsis thaliana and Medicago truncatula

S-adenosylmethionine synthases (MATs) are responsible for production of S-adenosylmethionine, the cofactor essential for various methylation reactions, production of polyamines and phytohormone ethylene, etc. Plants have two distinct MAT types (I and II). This work presents the structural analysis of MATs from Arabidopsis thaliana (AtMAT1 and AtMAT2, both type I) and Medicago truncatula (MtMAT3a, type II), which, unlike most MATs from other domains of life, are dimers where three-domain subunits are sandwiched flat with one another. Although MAT types are very similar, their subunits are differently oriented within the dimer. Structural snapshots along the enzymatic reaction reveal the exact conformation of precatalytic methionine in the active site and show a binding niche, characteristic only for plant MATs, that may serve as a lock of the gate loop. Nevertheless, plants, in contrary to mammals, lack the MAT regulatory subunit, and the regulation of plant MAT activity is still puzzling. Our structures open a possibility of an allosteric activity regulation of type I plant MATs by linear compounds, like polyamines, which would tighten the relationship between S-adenosylmethionine and polyamine biosynthesis.

Quaternary Structure of the Tryptophan Synthase α-Subunit Homolog BX1 from Zea mays

BX1 from Zea mays (zmBX1) is an enzyme of plant secondary metabolism that generates indole for the synthesis of plant defensins. It is a homologue of the tryptophan synthase α-subunit, TrpA. Whereas TrpA itself is a monomer in solution, zmBX1 is dimeric, confirmed in our work by native MS. Using cross-linking and mutagenesis, we identified the physiological dimerization interface of zmBX1. We found that homodimerization has only minor effects on catalysis and stability. A comparison of the zmBX1-zmBX1 homodimer and zmTrpA-zmTrpB heterodimer interfaces suggest that homodimerization in zmBX1 might, at an early point in evolution, have served as a mechanism to exclude the interaction with the tryptophan synthase β-subunit (zmTrpB), marking its transition from primary to secondary metabolism.

Genomic Origin and Diversification of the Glucosinolate MAM Locus

Glucosinolates are a diverse group of plant metabolites that characterize the order Brassicales. The MAM locus is one of the most significant QTLs for glucosinolate diversity. However, most of what we understand about evolution at the locus is focused on only a few species and not within a phylogenetic context. In this study, we utilize a micro-synteny network and phylogenetic inference to investigate the origin and diversification of the MAM/IPMS gene family. We uncover unique MAM-like genes found at the orthologous locus in the Cleomaceae that shed light on the transition from IPMS to MAM. In the Brassicaceae, we identify six distinct MAM clades across Lineages I, II, and III. We characterize the evolutionary impact and consequences of local duplications, transpositions, whole genome duplications, and gene fusion events, generating several new hypothesizes on the function and diversity of the MAM locus.

The Phytotoxicity of Meta-Tyrosine Is Associated With Altered Phenylalanine Metabolism and Misincorporation of This Non-Proteinogenic Phe-Analog to the Plant's Proteome

Plants produce a myriad of specialized (secondary) metabolites that are highly diverse chemically, and exhibit distinct biological functions. Here, we focus on meta-tyrosine (m-tyrosine), a non-proteinogenic byproduct that is often formed by a direct oxidation of phenylalanine (Phe). Some plant species (e.g., Euphorbia myrsinites and Festuca rubra) produce and accumulate high levels of m-tyrosine in their root-tips via enzymatic pathways. Upon its release to soil, the Phe-analog, m-tyrosine, affects early post-germination development (i.e., altered root development, cotyledon or leaf chlorosis, and retarded growth) of nearby plant life. However, the molecular basis of m-tyrosine-mediated (phyto)toxicity remains, to date, insufficiently understood and are still awaiting their functional characterization. It is anticipated that upon its uptake, m-tyrosine impairs key metabolic processes, or affects essential cellular activities in the plant. Here, we provide evidences that the phytotoxic effects of m-tyrosine involve two distinct molecular pathways. These include reduced steady state levels of several amino acids, and in particularly altered biosynthesis of the phenylalanine (Phe), an essential α-amino acid, which is also required for the folding and activities of proteins. In addition, proteomic studies indicate that m-tyrosine is misincorporated in place of Phe, mainly into the plant organellar proteomes. These data are supported by analyses of adt mutants, which are affected in Phe-metabolism, as well as of var2 mutants, which lack FtsH2, a major component of the chloroplast FtsH proteolytic machinery, which show higher sensitivity to m-tyrosine. Plants treated with m-tyrosine show organellar biogenesis defects, reduced respiration and photosynthetic activities and growth and developmental defect phenotypes.

Metabolic diversification of nitrogen-containing metabolites by the expression of a heterologous lysine decarboxylase gene in Arabidopsis

Lysine decarboxylase converts l-lysine to cadaverine as a branching point for the biosynthesis of plant Lys-derived alkaloids. Although cadaverine contributes towards the biosynthesis of Lys-derived alkaloids, its catabolism, including metabolic intermediates and the enzymes involved, is not known. Here, we generated transgenic Arabidopsis lines by expressing an exogenous lysine/ornithine decarboxylase gene from Lupinus angustifolius (La-L/ODC) and identified cadaverine-derived metabolites as the products of the emerged biosynthetic pathway. Through untargeted metabolic profiling, we observed the upregulation of polyamine metabolism, phenylpropanoid biosynthesis and the biosynthesis of several Lys-derived alkaloids in the transgenic lines. Moreover, we found several cadaverine-derived metabolites specifically detected in the transgenic lines compared with the non-transformed control. Among these, three specific metabolites were identified and confirmed as 5-aminopentanal, 5-aminopentanoate and δ-valerolactam. Cadaverine catabolism in a representative transgenic line (DC29) was traced by feeding stable isotope-labeled [α-¹⁵ N]- or [ε-¹⁵ N]-l-lysine. Our results show similar ¹⁵ N incorporation ratios from both isotopomers for the specific metabolite features identified, indicating that these metabolites were synthesized via the symmetric structure of cadaverine. We propose biosynthetic pathways for the metabolites on the basis of metabolite chemistry and enzymes known or identified through catalyzing specific biochemical reactions in this study. Our study shows that this pool of enzymes with promiscuous activities is the driving force for metabolite diversification in plants. Thus, this study not only provides valuable information for understanding the catabolic mechanism of cadaverine but also demonstrates that cadaverine accumulation is one of the factors to expand plant chemodiversity, which may lead to the emergence of Lys-derived alkaloid biosynthesis.

Multiple genes recruited from hormone pathways partition maize diterpenoid defences

Duplication and divergence of primary pathway genes underlie the evolution of plant specialized metabolism; however, mechanisms partitioning parallel hormone and defence pathways are often speculative. For example, the primary pathway intermediate ent-kaurene is essential for gibberellin biosynthesis and is also a proposed precursor for maize antibiotics. By integrating transcriptional coregulation patterns, genome-wide association studies, combinatorial enzyme assays, proteomics and targeted mutant analyses, we show that maize kauralexin biosynthesis proceeds via the positional isomer ent-isokaurene formed by a diterpene synthase pair recruited from gibberellin metabolism. The oxygenation and subsequent desaturation of ent-isokaurene by three promiscuous cytochrome P450s and a new steroid 5α reductase indirectly yields predominant ent-kaurene-associated antibiotics required for Fusarium stalk rot resistance. The divergence and differential expression of pathway branches derived from multiple duplicated hormone-metabolic genes minimizes dysregulation of primary metabolism via the circuitous biosynthesis of ent-kaurene-related antibiotics without the production of growth hormone precursors during defence.

Leaf metabolic signatures induced by real and simulated herbivory in black mustard (Brassica nigra)

INTRODUCTION

The oxylipin methyl jasmonate (MeJA) is a plant hormone active in response signalling and defence against herbivores. Although MeJA is applied experimentally to mimic herbivory and induce plant defences, its downstream effects on the plant metabolome are largely uncharacterized, especially in the context of primary growth and tissue-specificity of the response.

OBJECTIVES

We investigated the effects of MeJA-simulated and real caterpillar herbivory on the foliar metabolome of the wild plant Brassica nigra and monitored the herbivore-induced responses in relation to leaf ontogeny.

METHODS

As single or multiple herbivory treatments, MeJA- and mock-sprayed plants were consecutively exposed to caterpillars or left untreated. Gas chromatography (GC) and liquid chromatography (LC) time-of-flight mass-spectrometry (TOF-MS) were combined to analyse foliar compounds, including central primary and specialized defensive plant metabolites.

RESULTS

Plant responses were stronger in young leaves, which simultaneously induced higher chlorophyll levels. Both MeJA and caterpillar herbivory induced similar, but not identical, accumulation of tricarboxylic acids (TCAs), glucosinolates (GSLs) and phenylpropanoids (PPs), but only caterpillar feeding led to depletion of amino acids. MeJA followed by caterpillars caused higher induction of defence compounds, including a three-fold increase in the major defence compound allyl-GSL (sinigrin). When feeding on MeJA-treated plants, caterpillars gained less weight indicative of the reduced host-plant quality and enhanced resistance.

CONCLUSIONS

The metabolomics approach showed that plant responses induced by herbivory extend beyond the regulation of defence metabolism and are tightly modulated throughout leaf development. This leads to a new understanding of the plant metabolic potential that can be exploited for future plant protection strategies.

Harnessing evolutionary diversification of primary metabolism for plant synthetic biology

Plants produce numerous natural products that are essential to both plant and human physiology. Recent identification of genes and enzymes involved in their biosynthesis now provides exciting opportunities to reconstruct plant natural product pathways in heterologous systems through synthetic biology. The use of plant chassis, although still in infancy, can take advantage of plant cells' inherent capacity to synthesize and store various phytochemicals. Also, large-scale plant biomass production systems, driven by photosynthetic energy production and carbon fixation, could be harnessed for industrial-scale production of natural products. However, little is known about which plants could serve as ideal hosts and how to optimize plant primary metabolism to efficiently provide precursors for the synthesis of desirable downstream natural products or specialized (secondary) metabolites. Although primary metabolism is generally assumed to be conserved, unlike the highly-diversified specialized metabolism, primary metabolic pathways and enzymes can differ between microbes and plants and also among different plants, especially at the interface between primary and specialized metabolisms. This review highlights examples of the diversity in plant primary metabolism and discusses how we can utilize these variations in plant synthetic biology. I propose that understanding the evolutionary, biochemical, genetic, and molecular bases of primary metabolic diversity could provide rational strategies for identifying suitable plant hosts and for further optimizing primary metabolism for sizable production of natural and bio-based products in plants.

Dynamic metabolic solutions to the sessile life style of plants

Plants are sessile organisms. To compensate for not being able to escape when challenged by unfavorable growth conditions, pests or herbivores, plants have perfected their metabolic plasticity by having developed the capacity for on demand dynamic biosynthesis and storage of a plethora of phytochemicals.

Covering: up to 2018

Plants are sessile organisms. To compensate for not being able to escape when challenged by unfavorable growth conditions, pests or herbivores, plants have perfected their metabolic plasticity by having developed the capacity for on demand synthesis of a plethora of phytochemicals to specifically respond to the challenges arising during plant ontogeny. Key steps in the biosynthesis of phytochemicals are catalyzed by membrane-bound cytochrome P450 enzymes which in plants constitute a superfamily. In planta, the P450s may be organized in dynamic enzyme clusters (metabolons) and the genes encoding the P450s and other enzymes in a specific pathway may be clustered. Metabolon formation facilitates transfer of substrates between sequential enzymes and therefore enables the plant to channel the flux of general metabolites towards biosynthesis of specific phytochemicals. In the plant cell, compartmentalization of the operation of specific biosynthetic pathways in specialized plastids serves to avoid undesired metabolic cross-talk and offers distinct storage sites for molar concentrations of specific phytochemicals. Liquid–liquid phase separation may lead to formation of dense biomolecular condensates within the cytoplasm or vacuole allowing swift activation of the stored phytochemicals as required upon pest or herbivore attack. The molecular grid behind plant plasticity offers an endless reservoir of functional modules, which may be utilized as a synthetic biology tool-box for engineering of novel biological systems based on rational design principles. In this review, we highlight some of the concepts used by plants to coordinate biosynthesis and storage of phytochemicals.

SuCComBase: a manually curated repository of plant sulfur-containing compounds

Plants produce a wide range of secondary metabolites that play important roles in plant defense and immunity, their interaction with the environment and symbiotic associations. Sulfur-containing compounds (SCCs) are a group of important secondary metabolites produced in members of the Brassicales order. SCCs constitute various groups of phytochemicals, but not much is known about them. Findings from previous studies on SCCs were scattered in published literatures, hence SuCComBase was developed to store all molecular information related to the biosynthesis of SCCs. Information that includes genes, proteins and compounds that are involved in the SCC biosynthetic pathway was manually identified from databases and published scientific literatures. Sets of co-expression data was analyzed to search for other possible (previously unknown) genes that might be involved in the biosynthesis of SCC. These genes were named as potential SCC-related encoding genes. A total of 147 known and 92 putative Arabidopsis thaliana SCC-related genes from literatures were used to identify other potential SCC-related encoding genes. We identified 778 potential SCC-related encoding genes, 4026 homologs to the SCC-related encoding genes and 116 SCCs as shown on SuCComBase homepage. Data entries are searchable from the Main page, Search, Browse and Datasets tabs. Users can easily download all data stored in SuCComBase. All publications related to SCCs are also indexed in SuCComBase, which is currently the first and only database dedicated to plant SCCs. SuCComBase aims to become a manually curated and au fait knowledge-based repository for plant SCCs.

Flavor-cyber-agriculture: Optimization of plant metabolites in an open-source control environment through surrogate modeling

Food production in conventional agriculture faces numerous challenges such as reducing waste, meeting demand, maintaining flavor, and providing nutrition. Contained environments under artificial climate control, or cyber-agriculture, could in principle be used to meet many of these challenges. Through such environments, phenotypic expression of the plant—mass, edible yield, flavor, and nutrients—can be actuated through a “climate recipe,” where light, water, nutrients, temperature, and other climate and ecological variables are optimized to achieve a desired result. This paper describes a method for doing this optimization for the desired result of flavor by combining cyber-agriculture, metabolomic phenotype (chemotype) measurements, and machine learning. In a pilot experiment, (1) environmental conditions, i.e. photoperiod and ultraviolet (UV) light (known to affect production of flavor-active molecules in edible plants) were applied under different regimes to basil plants (Ocimum basilicum) growing inside a hydroponic farm with an open-source design; (2) flavor-active volatile molecules were measured in each plant using gas chromatography-mass spectrometry (GC-MS); and (3) symbolic regression was used to construct a surrogate model of this chemistry from the input environmental variables, and this model was used to discover new combinations of photoperiod and UV light to increase this chemistry. These new combinations, or climate recipes, were then implemented in the hydroponic farm, and several of them resulted in a marked increase in volatiles over control. The process also led to two important insights: it demonstrated a “dilution effect”, i.e. a negative correlation between weight and desirable chemical species, and it discovered the surprising effect that a 24-hour photoperiod of photosynthetic-active radiation, the equivalent of all-day light, induces the most flavor molecule production in basil. In this manner, surrogate optimization through machine learning can be used to discover effective recipes for cyber-agriculture that would be difficult and time-consuming to find using hand-designed experiments.

Evolution of metabolic novelty: A trichome-expressed invertase creates specialized metabolic diversity in wild tomato

Plants produce a myriad of taxonomically restricted specialized metabolites. This diversity-and our ability to correlate genotype with phenotype-makes the evolution of these ecologically and medicinally important compounds interesting and experimentally tractable. Trichomes of tomato and other nightshade family plants produce structurally diverse protective compounds termed acylsugars. While cultivated tomato (Solanum lycopersicum) strictly accumulates acylsucroses, the South American wild relative Solanum pennellii produces copious amounts of acylglucoses. Genetic, transgenic, and biochemical dissection of the S. pennellii acylglucose biosynthetic pathway identified a trichome gland cell-expressed invertase-like enzyme that hydrolyzes acylsucroses (Sopen03g040490). This enzyme acts on the pyranose ring-acylated acylsucroses found in the wild tomato but not on the furanose ring-decorated acylsucroses of cultivated tomato. These results show that modification of the core acylsucrose biosynthetic pathway leading to loss of furanose ring acylation set the stage for co-option of a general metabolic enzyme to produce a new class of protective compounds.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

Robust predictions of specialized metabolism genes through machine learning

Specialized metabolites are critical for plant–environment interactions, e.g., attracting pollinators or defending against herbivores, and are important sources of plant-based pharmaceuticals. However, it is unclear what proportion of enzyme-encoding genes play a role in specialized metabolism (SM) as opposed to general metabolism (GM) in any plant species. This is because of the diversity of specialized metabolites and the considerable number of incompletely characterized pathways responsible for their production. In addition, SM gene ancestors frequently played roles in GM. We evaluate features distinguishing SM and GM genes and build a computational model that accurately predicts SM genes. Our predictions provide candidates for experimental studies, and our modeling approach can be applied to other species that produce medicinally or industrially useful compounds.

Plant specialized metabolism (SM) enzymes produce lineage-specific metabolites with important ecological, evolutionary, and biotechnological implications. Using Arabidopsis thaliana as a model, we identified distinguishing characteristics of SM and GM (general metabolism, traditionally referred to as primary metabolism) genes through a detailed study of features including duplication pattern, sequence conservation, transcription, protein domain content, and gene network properties. Analysis of multiple sets of benchmark genes revealed that SM genes tend to be tandemly duplicated, coexpressed with their paralogs, narrowly expressed at lower levels, less conserved, and less well connected in gene networks relative to GM genes. Although the values of each of these features significantly differed between SM and GM genes, any single feature was ineffective at predicting SM from GM genes. Using machine learning methods to integrate all features, a prediction model was established with a true positive rate of 87% and a true negative rate of 71%. In addition, 86% of known SM genes not used to create the machine learning model were predicted. We also demonstrated that the model could be further improved when we distinguished between SM, GM, and junction genes responsible for reactions shared by SM and GM pathways, indicating that topological considerations may further improve the SM prediction model. Application of the prediction model led to the identification of 1,220 A. thaliana genes with previously unknown functions, each assigned a confidence measure called an SM score, providing a global estimate of SM gene content in a plant genome.