Search

Acylsugars are defensive, trichome-synthesized sugar esters produced in plants across the Solanaceae (nightshade) family. Although assembled from simple metabolites and synthesized by a relatively short core biosynthetic pathway, tremendous within- and across-species acylsugar structural variation is documented across the family. To advance our understanding of the diversity and the synthesis of acylsugars within the Nicotiana genus, trichome extracts were profiled across the genus coupled with transcriptomics-guided enzyme discovery and in vivo and in vitro analysis. Differences in the types of sugar cores, numbers of acylations, and acyl chain structures contributed to over 300 unique annotated acylsugars throughout Nicotiana. Placement of acyl chain length into a phylogenetic context revealed that an unsaturated acyl chain type was detected in a few closely related species. A comparative transcriptomics approach identified trichome-enriched Nicotiana acuminata acylsugar biosynthetic candidate enzymes. More than 25 acylsugar variants could be produced in a single enzyme assay with four N. acuminata acylsugar acyltransferases (NacASAT1–4) together with structurally diverse acyl-CoAs and sucrose. Liquid chromatography coupled with mass spectrometry screening of in vitro products revealed the ability of these enzymes to make acylsugars not present in Nicotiana plant extracts. In vitro acylsugar production also provided insights into acyltransferase acyl donor promiscuity and acyl acceptor specificity as well as regiospecificity of some ASATs. This study suggests that promiscuous Nicotiana acyltransferases can be used as synthetic biology tools to produce novel and potentially useful metabolites.

Plant derived bioactive small molecules have attracted attention of scientists across fundamental and applied scientific disciplines. We seek to understand the influence of these phytochemicals on functional phytobiomes. Increased knowledge of specialized metabolite bioactivities could inform strategies for sustainable crop production. Our recent investigations of metabolomes of the upland and lowland ecotypes of the bioenergy crop switchgrass (Panicum virgatum) revealed large differences in types and abundances of terpenoid specialized metabolites. We hypothesized that – consistent with accumulating evidence that switchgrass genotype impacts microbiome assembly – differential terpenoid accumulation contributes to switchgrass ecotype-specific microbiome composition. An initialin vitroplate-based disc diffusion screen of 18 switchgrass root derived fungal isolates revealed differential responses to upland- and lowland-isolated metabolites. To identify specific fungal growth-modulating metabolites, we tested fractions from root extracts on three ecologically important fungal isolates –Linnemania elongata, Trichodermasp. andFusariumsp. Saponins and diterpenoids were identified as the most prominent antifungal metabolites. Finally, analysis of liquid chromatography-purified terpenoids revealed fungal inhibition structure – activity relationships (SAR). Saponin antifungal activity was primarily determined by the number of sugar moieties – saponins glycosylated at a single core position were inhibitory whereas saponins glycosylated at two core positions were inactive. Saponin core hydroxylation and acetylation were also associated with reduced activity. Diterpenoid activity required the presence of an intact furan ring for strong fungal growth inhibition. This study demonstrates that switchgrass genotype-specific metabolites differentially inhibit fungal isolates from the root and rhizosphere, supporting the hypothesis that these small molecules contribute to microbiome assembly and function.

SUMMARYSwitchgrass (Panicum virgatum) is a bioenergy model crop valued for its energy efficiency and drought tolerance resilience. The related monocot species rice (Oryza sativa) and maize (Zea mays) deploy species-specific, specialized metabolites as core stress defenses. By contrast, specialized chemical defenses in switchgrass are largely unknown.To investigate specialized metabolic drought responses in switchgrass, we integrated tissue-specific transcriptome and metabolite analyses of the genotypes Alamo and Cave-in-Rock that feature different drought tolerance.The more drought-susceptible Cave-in-Rock featured an earlier onset of transcriptomic changes and significantly more differentially expressed genes in response to drought compared to Alamo. Specialized pathways showed moderate differential expression compared to pronounced transcriptomic alterations in carbohydrate and amino acid metabolism. However, diterpenoid-biosynthetic genes showed drought-inducible expression in Alamo roots, contrasting largely unaltered triterpenoid and phenylpropanoid pathways. Metabolomic analyses identified common and genotype-specific flavonoids and terpenoids. Consistent with transcriptomic alterations, several root diterpenoids showed significant drought-induced accumulation, whereas triterpenoid abundance remained predominantly unchanged. Structural analysis of drought-responsive root diterpenoids verified these metabolites as oxygenated furanoditerpenoids.Drought-dependent transcriptome and metabolite profiles provide the foundation to understand the molecular mechanisms underlying switchgrass environmental resilience. Accumulation of specialized root diterpenoids and corresponding pathway transcripts supports a role in drought stress tolerance for these compounds.Significance statementWith an increasing demand for renewable energy opposed by rising climate-driven crop losses, understanding, and leveraging plant natural defenses can enable the development of sustainable crop production systems. Here, we integrated comparative transcriptomics and metabolomics analyses to gain a detailed understanding of the diversity and physiological relevance of specialized metabolites in upland and lowland switchgrass ecotypes and provide resources for future investigations of drought response mechanisms in switchgrass.

The hormone salicylic acid (SA) plays crucial roles in plant defense, stress responses, and in the regulation of plant growth and development. Whereas the biosynthetic pathways and biological functions of SA have been extensively studied, SA catabolism is less well understood. In this study, we report the identification and functional characterization of an FAD/NADH-dependent SA 1-hydroxylase from tomato (Solanum lycopersicum; SlSA1H), which catalyzes the oxidative decarboxylation of SA to catechol. Transcript levels of SlSA1H were highest in stems and its expression was correlated with the formation of the methylated catechol derivatives guaiacol and veratrole. Consistent with a role in SA catabolism, SlSA1H RNAi plants accumulated lower amounts of guaiacol and failed to produce any veratrole. Two O-methyltransferases involved in the conversion of catechol to guaiacol and guaiacol to veratrole were also functionally characterized. Subcellular localization analyses revealed the cytosolic localization of this degradation pathway. Phylogenetic analysis and functional characterization of SA1H homologs from other species indicated that this type of FAD/NADH-dependent SA 1-hydroxylases evolved recently within the Solanaceae family.

Given the potential health benefits (and adverse effects), of polyphenolic and steroidal glycoalkaloids in the diet there is a growing interest in fully elucidating the genetic control of their levels in foodstuffs. Here we carried out profiling of the specialized metabolites in the seeds of the Solanum pennellii introgression lines identifying 338 putative metabolite quantitative trait loci (mQTL) for flavonoids, steroidal glycoalkaloids and further specialized metabolites. Two putative mQTL for flavonols and one for steroidal glycoalkaloids were cross-validated by evaluation of the metabolite content of recombinants harboring smaller introgression in the corresponding QTL interval or by analysis of lines from an independently derived backcross inbred line population. The steroidal glycoalkaloid mQTL was localized to a chromosomal region spanning 14 genes, including a previously defined steroidal glycoalkaloid gene cluster. The flavonoid mQTL was further vali dated via the use of transient and stable overexpression of the Solyc12g098600 and Solyc12g096870 genes, which encode seed-specific uridine 50 -diphosphate-glycosyltransferases. The results are discussed in the context of our understanding of the accumulation of polyphenols and steroidal glycoalkaloids, and how this knowledge may be incorporated into breeding strategies aimed at improving nutritional aspects of plants as well as in fortifying them against abiotic stress.

Plants make a variety of specialized metabolites that can mediate interactions with animals, microbes and competitor plants. Understanding how plants synthesize these compounds enables studies of their biological roles by manipulating their synthesisin vivoas well as producing themin vitro. Acylsugars are a group of protective metabolites that accumulate in the trichomes of many Solanaceae family plants. Acylinositol biosynthesis is of interest because it appears to be restricted to a subgroup of species within the Solanum genus. Previous work characterized a triacylinositol acetyltransferase involved in acylinositol biosynthesis in the Andean fruit plantSolanum quitoense(lulo or naranjillo). We characterized three additionalS. quitoensetrichome expressed enzymes, and found that virus induced gene silencing of each caused changes in acylinositol accumulation. Surprisingly, thein vitrotriacylinositol products of these enzymes are distinct from those that accumulatein planta. These enzymes, nonetheless, provide an opportunities to test the biological impact and properties of these triacylinositolsin vitro.

Acylsugars are defensive, trichome-synthesized sugar esters produced in plants across the Solanaceae (nightshade) family. Although assembled from simple metabolites and synthesized by a relatively short core biosynthetic pathway, tremendous within- and across-species acylsugar structural variation is documented across the family. To advance our understanding of the diversity and the synthesis of acylsugars within the Nicotiana genus, trichome extracts were profiled across the genus coupled with transcriptomics-guided enzyme discovery and in vivo and in vitro analysis. Differences in the types of sugar cores, numbers of acylations, and acyl chain structures contributed to over 300 unique annotated acylsugars throughout Nicotiana. Placement of acyl chain length into a phylogenetic context revealed that an unsaturated acyl chain type was detected in a few closely-related species. A comparative transcriptomics approach identified trichome-enriched Nicotiana acuminata acylsugar biosynthetic candidate enzymes. > 25 acylsugar variants could be produced in a single enzyme assay with four acylsugar acyltransferases (NacASAT1-4) together with structurally diverse acyl-CoAs and sucrose. Liquid chromatography coupled with mass spectrometry screening of in vitro products revealed the ability of these enzymes to make acylsugars not present in Nicotiana plant extracts. In vitro acylsugar production also provided insights into acyltransferase acyl donor promiscuity and acyl acceptor specificity as well as regiospecificity of some ASATs. This study suggests that promiscuous Nicotiana acyltransferases can be used as synthetic biology tools to produce novel and potentially useful metabolites.ONE SENTENCE SUMMARYAnalysis of Nicotiana glandular trichome metabolites and BAHD acyltransferases revealed diverse sucrose and glucose based acylesters.

Description, Analysis of plant protective surface hair chemistry revealed evolutionary mechanisms leading to metabolic innovation., 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.

Correction for ‘Fruity, sticky, stinky, spicy, bitter, addictive, and deadly: evolutionary signatures of metabolic complexity in the Solanaceae’ by Paul D. Fiesel et al., Nat. Prod. Rep., 2022, 39, 1438–1464, https://doi.org/10.1039/D2NP00003B.

Switchgrass (Panicum virgatum L.) is a bioenergy crop that grows productively on lands not suitable for food production, and is an excellent target for low-pesticide input biomass production. We hypothesize that resistance to insect pests and microbial pathogens is influenced by low molecular weight compounds known as specialized metabolites. We employed untargeted liquid chromatography-mass spectrometry (LC-MS), quantitative gas chromatography-mass spectrometry (GC-MS) and nuclear magnetic resonance (NMR) spectroscopy to identify differences in switchgrass ecotype metabolomes. This analysis revealed striking differences between upland and lowland switchgrass metabolomes as well as distinct developmental profiles. Terpenoid and polyphenol derived specialized metabolites were identified, including steroidal saponins, di- and sesqui-terpenoids and flavonoids. The saponins are especially abundant in switchgrass extracts and have diverse aglycone cores and sugar moieties. We report seven structurally distinct steroidal saponin classes with unique steroidal cores and glycosylated at one or two positions. Quantitative GC-MS revealed differences in total saponin concentrations in leaf blade, leaf sheath, stem, rhizome and root. The quantitative data also demonstrated that saponin concentrations is higher in roots of lowland than upland ecotype plants, suggesting ecotypic specific biosynthesis and/or biological functions. These results enable future testing of these specialized metabolites on biotic and abiotic stress tolerance and can inform development of low-input bioenergy crops.One sentence summaryIntegrated mass-spectrometry and nuclear magnetic resonance spectroscopy based metabolomics reveal that switchgrass accumulates structurally diverse terpenoids and phenolics, which vary in abundance and structure in a tissue- and ecotype-specific manner.

Natural pyrethrin insecticides produced by Dalmatian pyrethrum (Tanacetum cinerariifolium) have low mammalian toxicity and short environmental persistence, providing an alternative to widely used synthetic agricultural insecticides that pose a threat to human health and the environment. A recent surge of interest in the use of pyrethrins as agricultural insecticides coincides with the discovery of several new genes in the pyrethrin biosynthetic pathway. Elucidation of this pathway facilitates efforts to breed improved pyrethrum varieties and to engineer plants with improved endogenous defenses or hosts for heterologous pyrethrin production. We describe the current state of knowledge related to global pyrethrum production, the pyrethrin biosynthetic pathway and its regulation, and recent efforts to engineer the pyrethrin pathway in diverse plant hosts.

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., eLife digest Plants produce a vast variety of different molecules known as secondary or specialized metabolites to attract pollinating insects, such as bees, or protect themselves against herbivores and pests. The secondary metabolites are made from simple building blocks that are readily available in plants, including amino acids, fatty acids and sugars. Different species of plant, and even different parts of the same plant, produce their own sets of secondary metabolites. For example, the hairs on the surface of tomatoes and other members of the nightshade family of plants make metabolites known as acylsugars. These chemicals deter herbivores and pests from damaging the plants. To make acylsugars, the plants attach long chains known as fatty acyl groups to molecules of sugar, such as sucrose. Some members of the nightshade family produce acylsugars with longer chains than others. In particular, acylsugars with long chains are only found in tomatoes and other closely-related species. It remained unclear how the nightshade family evolved to produce acylsugars with chains of different lengths. To address this question, Fan et al. used genetic and biochemical approaches to study tomato plants and other members of the nightshade family. The experiments identified two genes known as AACS and AECH in tomatoes that produce acylsugars with long chains. These two genes originated from the genes of older enzymes that metabolize fatty acids – the building blocks of fats – in plant cells. Unlike the older genes, AACS and AECH were only active at the tips of the hairs on the plant’s surface. Fan et al. then investigated the evolutionary relationship between 11 members of the nightshade family and two other plant species. This revealed that AACS and AECH emerged in the nightshade family around the same time that longer chains of acylsugars started appearing. These findings provide insights into how plants evolved to be able to produce a variety of secondary metabolites that may protect them from a broader range of pests. The gene cluster identified in this work could be used to engineer other species of crop plants to start producing acylsugars as natural pesticides.

Plants make hundreds of thousands of biologically active specialized metabolites varying widely in structure, biosynthesis and the processes that they influence. An increasing number of these compounds are documented to protect plants from harmful insects, pathogens, or herbivores, or mediate interactions with beneficial organisms including pollinators and nitrogen fixing microbes. Acylsugars – one class of protective compounds – are made in glandular trichomes of plants across the Solanaceae family. While most described acylsugars are acylsucroses, published examples also include acylsugars with hexose cores. The South American fruit crop Solanum quitoense (Naranjilla) produces acylsugars that contain a myo-inositol core. We identified an enzyme that acetylates triacylinositols, a function homologous to the last step in the Solanum lycopersicum acylsucrose biosynthetic pathway. Our analysis reveals parallels between S. lycopersicum acylsucrose and S. quitoense acylinositol biosynthesis, suggesting a common evolutionary origin.Material availabilityThe author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantphysiol.org) is: Robert L. Last (lastr@msu.edu).One sentence summaryEvidence that the final step in Solanum quitoense acylinositol biosynthesis evolved from an acylsucrose acetyltransferase enzyme.

Plants make many biologically active, specialized metabolites, which vary in structure, biosynthesis, and the processes they influence. An increasing number of these compounds are documented to protect plants from insects, pathogens, or herbivores or to mediate interactions with beneficial organisms, including pollinators and nitrogen-fixing microbes. Acylsugars, one class of protective compounds, are made in glandular trichomes of plants across the Solanaceae family. While most described acylsugars are acylsucroses, published examples also include acylsugars with hexose cores. The South American fruit crop naranjilla (lulo; Solanum quitoense) produces acylsugars containing a myoinositol core. We identified an enzyme that acetylates triacylinositols, a function homologous to the last step in the acylsucrose biosynthetic pathway of tomato (Solanum lycopersicum). Our analysis reveals parallels between S. lycopersicum acylsucrose and S. quitoense acylinositol biosynthesis, suggesting a common evolutionary origin.

Collectively, plants produce hundreds of thousands of specialized metabolites from simple building blocks such as amino acids, fatty acids, and isoprenoids. As additional specialized metabolic enzymes are described, there is increasing recognition of the importance of cooption of general metabolic enzymes to specialized metabolism by gene duplication, narrowing of expression, and alteration of enzymatic activities. Here, we examine how several classes of enzymes were each coopted multiple times. We demonstrate the simplicity of achieving the synthesis of analogous chemicals by coopting existing enzymes and summarize emerging insights that could inform rational metabolic engineering of both general and specialized metabolic enzymes.

The pyrethrum plant, Tanacetum cineramfolium (Asteraceae) synthesizes a class of compounds called pyrethrins that have strong insecticidal properties but are safe to humans. Class I pyrethrins are esters of the monoterpenoid trans-chrysanthemic acid with one of three jasmonic-acid derived alcohols. We reconstructed the trans-chry-santhemic acid biosynthetic pathway in tomato fruits, which naturally produce high levels of the tetraterpene pigment lycopene, an isoprenoid which shares a common precursor, dimethylallyl diphosphate (DMAPP), with trans-chrysanthemic acid. trans-Chrysanthemic acid biosynthesis in tomato fruit was achieved by expressing the chrysanthemyl diphosphate synthase gene from T. cinerariifolium, encoding the enzyme that uses DMAPP to make trans-chrysanthemol, under the control of the fruit specific promoter PG, as well as an alcohol dehy-drogenease (ADH) gene and aldehyde dehydrogenase (ALDH) gene from a wild tomato species, also under the control of the PG promoter. Tomato fruits expressing all three genes had a concentration of trans-chrysanthemic acid that was about 1.7-fold higher (by weight) than the levels of lycopene present in non-transgenic fruit, while the level of lycopene in the transgenic plants was reduced by 68%. Ninety seven percent of the diverted DMAPP was converted to trans-chrysanthemic acid, but 62% of this acid was further glycosylated. We conclude that the tomato fruit is an alternative platform for the biosynthesis of trans-chrysanthemic acid by metabolic engineering.

Specialized metabolic enzymes and metabolite diversity evolve through a variety of mechanisms including promiscuity, changes in substrate specificity, modifications of gene expression and gene duplication. For example, gene duplication and substrate binding site changes led to the evolution of the glucosinolate biosynthetic enzyme, AtIPMDH1, from a Leu biosynthetic enzyme. BAHD acyltransferases illustrate how enzymatic promiscuity leads to metabolite diversity. The examples 4-coumarate:CoA ligase and aromatic acid methyltransferases illustrate how promiscuity can potentiate the evolution of these specialized metabolic enzymes.

Flowers of Tanacetum cinerariifolium produce a set of compounds known collectively as pyrethrins, which are commercially important pesticides that are strongly toxic to flying insects but not to most vertebrates. A pyrethrin molecule is an ester consisting of either trans-chrysanthemic acid or its modified form, pyrethric acid, and one of three alcohols, jasmolone, pyrethrolone, and cinerolone, that appear to be derived from jasmonic acid. Chrysanthemyl diphosphate synthase (CDS), the first enzyme involved in the synthesis of trans-chrysanthemic acid, was characterized previously and its gene isolated. TcCDS produces free trans-chrysanthemol in addition to trans-chrysanthemyl diphosphate, but the enzymes responsible for the conversion of trans-chrysanthemol to the corresponding aldehyde and then to the acid have not been reported. We used an RNA sequencing-based approach and coexpression correlation analysis to identify several candidate genes encoding putative trans-chrysanthemol and trans-chrysanthemal dehydrogenases. We functionally characterized the proteins encoded by these genes using a combination of in vitro biochemical assays and heterologous expression in planta to demonstrate that TcADH2 encodes an enzyme that oxidizes trans-chrysanthemol to trans-chrysanthemal, while TcALDH1 encodes an enzyme that oxidizes trans-chrysanthemal into trans-chrysanthemic acid. Transient coexpression of TcADH2 and TcALDH1 together with TcCDS in Nicotiana benthamiana leaves results in the production of trans-chrysanthemic acid as well as several other side products. The majority (58%) of trans-chrysanthemic acid was glycosylated or otherwise modified. Overall, these data identify key steps in the biosynthesis of pyrethrins and demonstrate the feasibility of metabolic engineering to produce components of these defense compounds in a heterologous host.

The branched-chain amino acids (BCAAs) Ile, Val and Leu are essential nutrients that humans and other animals obtain from plants. However, total and relative amounts of plant BCAAs rarely match animal nutritional needs, and improvement requires a better understanding of the mechanistic basis for BCAA homeostasis. We present an in vivo regulatory model of BCAA homeostasis derived from analysis of feedback-resistant Arabidopsis thaliana mutants for the three allosteric committed enzymes in the biosynthetic network: threonine deaminase (also named L-O-methylthreonine resistant 1, OMR1), acetohydroxyacid synthase small subunit 2 (AHASS2) and isopropylmalate synthase 1 (IPMS1). In this model, OMR1 exerts primary control on Ile accumulation and functions independently of AHAS and IPMS. AHAS and IPMS regulate Val and Leu homeostasis, where AHAS affects total Val+Leu and IPMS controls partitioning between these amino acids. In addition, analysis of feedback-resistant and loss-of-function single and double mutants revealed that each AHAS and IPMS isoenzyme contributes to homeostasis rather than being functionally redundant. The characterized feedback-resistant mutations caused increased free BCAA levels in both seedlings and seeds. These results add to our understanding of the basis of in vivo BCAA homeostasis and inform approaches to improve the amount and balance of these essential nutrients in crops.

The plant pyrethrum (Tanacetum cinerariifolium) synthesizes highly effective natural pesticides known as pyrethrins. Pyrethrins are esters consisting of an irregular monoterpenoid acid and an alcohol derived from jasmonic acid (JA). These alcohols, referred to as rethrolones, can be jasmolone, pyrethrolone, or cinerolone. We recently showed that jasmolone is synthesized from jasmone, a degradation product of JA, in a single hydroxylation step catalyzed by jasmone hydroxylase (TcJMH). TcJMH belongs to the CYP71 clade of the cytochrome P450 oxidoreductase family. Here, we used coexpression analysis, heterologous gene expression, and in vitro biochemical assays to identify the enzyme responsible for conversion of jasmolone to pyrethrolone. A further T. cinerariifolium cytochrome P450 family member, CYP82Q3 (designated Pyrethrolone Synthase; TcPYS), appeared to catalyze the direct desaturation of the C1–C2 bond in the pentyl side chain of jasmolone to produce pyrethrolone. TcPYS is highly expressed in the trichomes of the ovaries in pyrethrum flowers, similar to TcJMH and other T. cinerariifolium genes involved in JA biosynthesis. Thus, as previously shown for biosynthesis of the monoterpenoid acid moiety of pyrethrins, rethrolones are synthesized in the trichomes. However, the final assembly of pyrethrins occurs in the developing achenes. Our data provide further insight into pyrethrin biosynthesis, which could ultimately be harnessed to produce this natural pesticide in a heterologous system.

Evolution of a cell specific invertase enzyme creates phenotypic diversity in wild tomato insecticidal metabolite production., 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.

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

Glandular trichomes of cultivated tomato (Solanum lycopersicum) and many other species throughout the Solanaceae produce and secrete mixtures of sugar esters (acylsugars) on the plant aerial surfaces. In wild and cultivated tomato, these metabolites consist of a sugar backbone, typically glucose or sucrose, and two to five acyl chains esterified to various positions on the sugar core. The aliphatic acyl chains vary in length and branching and are transferred to the sugar by a series of reactions catalyzed by acylsugar acyltransferases. A phenotypic screen of a set of S. lycopersicum M82 × Solanum pennellii LA0716 introgression lines identified a dominant genetic locus on chromosome 5 from the wild relative that affected total acylsugar levels. Genetic mapping revealed that the reduction in acylsugar levels was consistent with the presence and increased expression of two S. pennellii genes (Sopen05g030120 and Sopen05g030130) encoding putative carboxylesterase enzymes of the α/β-hydrolase superfamily. These two enzymes, named ACYLSUGAR ACYLHYDROLASE1 (ASH1) and ASH2, were shown to remove acyl chains from specific positions of certain types of acylsugars in vitro. A survey of related genes in M82 and LA0716 identified another trichome-expressed ASH gene on chromosome 9 (M82, Solyc09g075710; LA0716, Sopen09g030520) encoding a protein with similar activity. Characterization of the in vitro activities of the SpASH enzymes showed reduced activities with acylsugars produced by LA0716, presumably contributing to the high-level production of acylsugars in the presence of highly expressed SpASH genes.

The diversity of life on Earth is a result of continual innovations in molecular networks influencing morphology and physiology. Plant specialized metabolism produces hundreds of thousands of compounds, offering striking examples of these innovations. To understand how this novelty is generated, we investigated the evolution of the Solanaceae family-specific, trichome-localized acylsugar biosynthetic pathway using a combination of mass spectrometry, RNA-seq, enzyme assays, RNAi and phylogenomics in different non-model species. Our results reveal hundreds of acylsugars produced across the Solanaceae family and even within a single plant, built on simple sugar cores. The relatively short biosynthetic pathway experienced repeated cycles of innovation over the last 100 million years that include gene duplication and divergence, gene loss, evolution of substrate preference and promiscuity. This study provides mechanistic insights into the emergence of plant chemical novelty, and offers a template for investigating the ~300,000 non-model plant species that remain underexplored., eLife digest There are about 300,000 species of plant on Earth, which together produce over a million different small molecules called metabolites. Plants use many of these molecules to grow, to communicate with each other or to defend themselves against pests and disease. Humans have co-opted many of the same molecules as well; for example, some are important nutrients while others are active ingredients in medicines. Many plant metabolites are found in almost all plants, but hundreds of thousands of them are more specialized and only found in small groups of related plant species. These specialized metabolites have a wide variety of structures, and are made by different enzymes working together to carry out a series of biochemical reactions. Acylsugars are an example of a group of specialized metabolites with particularly diverse structures. These small molecules are restricted to plants in the Solanaceae family, which includes tomato and tobacco plants. Moghe et al. have now focused on acylsugars to better understand how plants produce the large diversity of chemical structures found in specialized metabolites, and how these processes have evolved over time. An analysis of over 35 plant species from across the Solanaceae family revealed hundreds of acylsugars, with some plants accumulating 300 or more different types of these specialized metabolites. Moghe et al. then looked at the enzymes that make acylsugars from a poorly studied flowering plant called Salpiglossis sinuata, partly because it produces a large diversity of these small molecules and partly because it sits in a unique position in the Solanaceae family tree. The activities of the enzymes were confirmed both in test tubes and in plants. This suggested that many of the enzymes were “promiscuous”, meaning that they could likely use a variety of molecules as starting points for their chemical reactions. This finding could help to explain how this plant species can make such a wide variety of acylsugars. Moghe et al. also discovered that many of the enzymes that make acylsugars are encoded by genes that were originally copies of other genes and that have subsequently evolved new activities. Plant scientists and plant breeders value tomato plants that produce acylsugars because these natural chemicals protect against pests like whiteflies and spider mites. A clearer understanding of the diversity of acylsugars in the Solanaceae family, as well as the enzymes that make these specialized metabolites, could help efforts to breed crops that are more resistant to pests. Some of the enzymes related to those involved in acylsugar production could also help to make chemicals with pharmaceutical value. These new findings might also eventually lead to innovative ways to produce these chemicals on a large scale.

The diversity of life on Earth is a result of continual innovations in molecular networks influencing morphology and physiology. Plant specialized metabolism produces hundreds of thousands of compounds, offering striking examples of these innovations. To understand how this novelty is generated, we investigated the evolution of the Solanaceae family-specific, trichome-localized acylsugar biosynthetic pathway using a combination of mass spectrometry, RNA-seq, enzyme assays, RNAi and phylogenetics in non-model species. Our results reveal that hundreds of acylsugars are produced across the Solanaceae family and even within a single plant, revealing this phenotype to be hyper-diverse. The relatively short biosynthetic pathway experienced repeated cycles of innovation over the last 100 million years that include gene duplication and divergence, gene loss, evolution of substrate preference and promiscuity. This study provides mechanistic insights into the emergence of plant chemical novelty, and offers a template for investigating the ∼300,000 non-model plant species that remain underexplored.

Flavonoids are ubiquitous plant aromatic specialized metabolites found in a variety of cell types and organs. Methylated flavonoids are detected in secreting glandular trichomes of various Solanum species, including the cultivated tomato (Solanum lycopersicum). Inspection of the sequenced S. lycopersicum Heinz 1706 reference genome revealed a close homolog of Solanum habrochaites MOMT1 3'/5' myricetin O-methyltransferase gene, but this gene (Solyc06g083450) is missing the first exon, raising the question of whether cultivated tomato has a distinct 3' or 3'/5' O-methyltransferase. A combination of mining genome and cDNA sequences from wild tomato species and S. lycopersicum cultivar M82 led to the identification of Sl-MOMT4 as a 3' O-methyltransferase. In parallel, three independent ethyl methanesulfonate mutants in the S. lycopersicum cultivar M82 background were identified as having reduced amounts of di- and trimethylated myricetins and increased monomethylated myricetin. Consistent with the hypothesis that Sl-MOMT4 is a 3' O-methyltransferase gene, all three myricetin methylation defective mutants were found to have defects in MOMT4 sequence, transcript accumulation, or 3'-O-methyltransferase enzyme activity. Surprisingly, no MOMT4 sequence is found in the Heinz 1706 reference genome sequence, and this cultivar accumulates 3-methyl myricetin and is deficient in 3'-methyl myricetins, demonstrating variation in this gene among cultivated tomato varieties.

Acylsugars are important insect defense compounds produced at high levels by glandular trichomes of the wild tomato, Solanum pennellii. The ability to produce acylsugars at elevated levels was bred into the tomato line CU071026. This study utilized a marker-assisted backcross approach to individually introgress into CU071026 and to fine map the three quantitative trait loci (QTL) fatty acid 5 (FA5QTL), fatty acid 7 (FA7QTL), and fatty acid 8 (FA8QTL), which were previously associated with changes in acylsugar chemistry. Additional breeding with and fine mapping the previously introgressed QTL, fatty acid 2 (FA2QTL), was also conducted. The effect of these four QTL on acylsugar quality and quantity in the presence of the five introgressions of CU071026 was evaluated. Incorporation of the QTL altered acylsugar chemotype by modulating the length, orientation, and/or relative proportion of fatty acid acyl groups. The resulting quantities of acylsugar produced in most of the new lines were similar to that of CU071026; however, introgression of FA5QTL reduced acylsugar levels. The acylsugar lines containing each QTL were characterized for acylsugar level, trichome abundance, and acylsugar chemistry through gas chromatography/mass spectrometry and liquid chromatography/mass spectrometry. The novel acylsugar chemotype lines created can contribute to elucidation of the mechanism of insect resistance mediated by acylsugars and help with identification of yet-unknown genes contributing to acylsugar synthesis and diversity.

A systematic screen of volatile terpene production in the glandular trichomes of 79 accessions of Solanum habrochaites was conducted and revealed the presence of 21 mono- and sesquiterpenes that exhibit a range of qualitative and quantitative variation. Hierarchical clustering identified distinct terpene phenotypic modules with shared patterns of terpene accumulation across accessions. Several terpene modules could be assigned to previously identified terpene synthase (TPS) activities that included members of the TPS-e/f subfamily that utilize the unusual cis-prenyl diphosphate substrates neryl diphosphate and 2z,6z-farnesyl diphosphate. DNA sequencing and in vitro enzyme activity analysis of TPS-e/f members from S. habrochaites identified three previously unassigned enzyme activities that utilize these cisoid substrates. These produce either the monoterpenes α-pinene and limonene, or the sesquiterpene 7-epizingiberene, with the in vitro analyses that recapitulated the trichome chemistry found in planta. Comparison of the distribution of S. habrochaites accessions with terpene content revealed a strong preference for the presence of particular TPS20 alleles at distinct geographic locations. This study reveals that the unusually high intra-specific variation of volatile terpene synthesis in glandular trichomes of S. habrochaites is due at least in part to evolution at the TPS20 locus.

This work identifies LOW QUANTUM YIELD OF PHOTOSYSTEM II1 (LQY1), a Zn finger protein that shows disulfide isomerase activity, interacts with the photosystem II (PSII) core complex, and may act in repair of photodamaged PSII complexes. Two mutants of an unannotated small Zn finger containing a thylakoid membrane protein of Arabidopsis thaliana (At1g75690; LQY1) were found to have a lower quantum yield of PSII photochemistry and reduced PSII electron transport rate following high-light treatment. The mutants dissipate more excess excitation energy via nonphotochemical pathways than wild type, and they also display elevated accumulation of reactive oxygen species under high light. After high-light treatment, the mutants have less PSII–light-harvesting complex II supercomplex than wild-type plants. Analysis of thylakoid membrane protein complexes showed that wild-type LQY1 protein comigrates with the PSII core monomer and the CP43-less PSII monomer (a marker for ongoing PSII repair and reassembly). PSII repair and reassembly involve the breakage and formation of disulfide bonds among PSII proteins. Interestingly, the recombinant LQY1 protein demonstrates a protein disulfide isomerase activity. LQY1 is more abundant in stroma-exposed thylakoids, where key steps of PSII repair and reassembly take place. The absence of the LQY1 protein accelerates turnover and synthesis of PSII reaction center protein D1. These results suggest that the LQY1 protein may be involved in maintaining PSII activity under high light by regulating repair and reassembly of PSII complexes.

Large-scale phenotypic screening presents challenges and opportunities not encountered in typical forward or reverse genetics projects. We describe a modular database and laboratory information management system that was implemented in support of the Chloroplast 2010 Project, an Arabidopsis (Arabidopsis thaliana) reverse genetics phenotypic screen of more than 5,000 mutants (http://bioinfo.bch.msu.edu/2010_LIMS; www.plastid.msu.edu). The software and laboratory work environment were designed to minimize operator error and detect systematic process errors. The database uses Ruby on Rails and Flash technologies to present complex quantitative and qualitative data and pedigree information in a flexible user interface. Examples are presented where the database was used to find opportunities for process changes that improved data quality. We also describe the use of the data-analysis tools to discover mutants defective in enzymes of leucine catabolism (heteromeric mitochondrial 3-methylcrotonyl-coenzyme A carboxylase [At1g03090 and At4g34030] and putative hydroxymethylglutaryl-coenzyme A lyase [At2g26800]) based upon a syndrome of pleiotropic seed amino acid phenotypes that resembles previously described isovaleryl coenzyme A dehydrogenase (At3g45300) mutants. In vitro assay results support the computational annotation of At2g26800 as hydroxymethylglutaryl-coenzyme A lyase.

Shotgun proteomics analysis allows hundreds of proteins to be identified and quantified from a single sample at relatively low cost. Extensive DNA sequence information is a prerequisite for shotgun proteomics, and it is ideal to have sequence for the organism being studied rather than from related species or accessions. While this requirement has limited the set of organisms that are candidates for this approach, next generation sequencing technologies make it feasible to obtain deep DNA sequence coverage from any organism. As part of our studies of specialized (secondary) metabolism in tomato (Solanum lycopersicum) trichomes, 454 sequencing of cDNA was combined with shotgun proteomics analyses to obtain in-depth profiles of genes and proteins expressed in leaf and stem glandular trichomes of 3-week-old plants. The expressed sequence tag and proteomics data sets combined with metabolite analysis led to the discovery and characterization of a sesquiterpene synthase that produces β-caryophyllene and α-humulene from E,E-farnesyl diphosphate in trichomes of leaf but not of stem. This analysis demonstrates the utility of combining high-throughput cDNA sequencing with proteomics experiments in a target tissue. These data can be used for dissection of other biochemical processes in these specialized epidermal cells.

In 1998 Cereon Genomics LLC, a subsidiary of Monsanto Co., performed a shotgun sequencing of the Arabidopsis thaliana Landsberg erecta genome to a depth of twofold coverage using 'classic' Sanger sequencing. This sequence was assembled and aligned to the Columbia ecotype sequence produced by the Arabidopsis Genome Initiative. The analysis provided tens of thousands of high-confidence predictions of polymorphisms between these two varieties of A. thaliana, and the predicted polymorphisms and Landsberg erecta sequence were subsequently made available to the not-for-profit research community by Monsanto. These data have been used for a wide variety of published studies, including map-based gene identification from forward genetic screens, studies of recombination and organelle genetics, and gene expression studies. The combination of resequencing approaches with next-generation sequencing technology has led to an increasing number of similar studies of genome-wide genetic diversity in A. thaliana, including the 1001 genomes project (http://1001genomes.org). Similar approaches are becoming possible in any number of crop species as DNA sequencing costs plummet and throughput rapidly increases, promising to lay the groundwork for revolutionizing our understanding of the relationship between genotype and phenotype in plants.

Summary Glandular secreting trichomes of cultivated tomato (Solanum lycopersicum) and close relatives produce a variety of structurally diverse volatile and non-volatile specialized (‘secondary’) metabolites, including terpenes, flavonoids and acyl sugars. A genetic screen is described here to profile leaf trichome and surface metabolite extracts of nearly isogenic chromosomal substitution lines covering the tomato genome. These lines contain specific regions of the Solanum pennellii LA0716 genome in an otherwise ‘wild-type’ M82 tomato genetic background. Regions that have an impact on the total amount of extractable mono- and sesquiterpenes (IL2-2) or only sesquiterpenes (IL10-3) or specifically influence accumulation of the monoterpene α-thujene (IL1-3 and IL1-4) were identified using GC-MS. A rapid LC-TOF-MS method was developed and used to identify changes in non-volatile metabolites through non-targeted analysis. Metabolite profiles generated using this approach led to the discovery of introgression lines producing different acyl chain substitutions on acyl sugar metabolites (IL1-3/1-4 and IL8-1/8-1-1), as well as two regions that influence the quantity of acyl sugars (IL5-3 and IL11-3). Chromosomal region 1-1/1-1-3 was found to influence the types of glycoalkaloids that are detected in leaf surface extracts. These results show that direct chemical screening is a powerful way to characterize genetic diversity in trichome specialized metabolism.

Polar membrane glycerolipids occur in a mixture of molecular species defined by a polar head group and characteristic acyl groups esterified to a glycerol backbone. A molecular species of phosphatidylglycerol specific to chloroplasts of plants carries a Delta(3-trans) hexadecenoic acid in the sn-2 position of its core glyceryl moiety. The fad4-1 mutant of Arabidopsis thaliana missing this particular phosphatidylglycerol molecular species lacks the necessary fatty acid desaturase, or a component thereof. The overwhelming majority of acyl groups associated with membrane lipids in plants contains double bonds with a cis configuration. However, FAD4 is unusual because it is involved in the formation of a trans double bond introduced close to the carboxyl group of palmitic acid, which is specifically esterified to the sn-2 glyceryl carbon of phosphatidylglycerol. As a first step towards the analysis of this unusual desaturase reaction, the FAD4 gene was identified by mapping of the FAD4 locus and coexpression analysis with known lipid genes. FAD4 encodes a predicted integral membrane protein that appears to be unrelated to classic membrane bound fatty acid desaturases based on overall sequence conservation. However, the FAD4 protein contains two histidine motifs resembling those of metalloproteins such as fatty acid desaturases. FAD4 is targeted to the plastid. Overexpression of the cDNA in transgenic Arabidopsis led to increased accumulation of the Delta(3-trans) hexadecanoyl group in phosphatidylglycerol relative to wild type. Taken together these results are consistent with the hypothesis that FAD4 is the founding member of a novel class of fatty acid desaturases.

Traditionally, phenotype-driven forward genetic plant mutant studies have been among the most successful approaches to revealing the roles of genes and their products and elucidating biochemical, developmental, and signaling pathways. A limitation is that it is time consuming, and sometimes technically challenging, to discover the gene responsible for a phenotype by map-based cloning or discovery of the insertion element. Reverse genetics is also an excellent way to associate genes with phenotypes, although an absence of detectable phenotypes often results when screening a small number of mutants with a limited range of phenotypic assays. The Arabidopsis Chloroplast 2010 Project (www.plastid.msu.edu) seeks synergy between forward and reverse genetics by screening thousands of sequence-indexed Arabidopsis (Arabidopsis thaliana) T-DNA insertion mutants for a diverse set of phenotypes. Results from this project are discussed that highlight the strengths and limitations of the approach. We describe the discovery of altered fatty acid desaturation phenotypes associated with mutants of At1g10310, previously described as a pterin aldehyde reductase in folate metabolism. Data are presented to show that growth, fatty acid, and chlorophyll fluorescence defects previously associated with antisense inhibition of synthesis of the family of acyl carrier proteins can be attributed to a single gene insertion in Acyl Carrier Protein4 (At4g25050). A variety of cautionary examples associated with the use of sequence-indexed T-DNA mutants are described, including the need to genotype all lines chosen for analysis (even when they number in the thousands) and the presence of tagged and untagged secondary mutations that can lead to the observed phenotypes.

In traditional mutant screening approaches, genetic variants are tested for one or a small number of phenotypes. Once bona fide variants are identified, they are typically subjected to a limited number of secondary phenotypic screens. Although this approach is excellent at finding genes involved in specific biological processes, the lack of wide and systematic interrogation of phenotype limits the ability to detect broader syndromes and connections between genes and phenotypes. It could also prevent detection of the primary phenotype of a mutant. As part of a systems biology approach to understand plastid function, large numbers of Arabidopsis thaliana homozygous T-DNA lines are being screened with parallel morphological, physiological, and chemical phenotypic assays (www.plastid.msu.edu). To refine our approaches and validate the use of this high-throughput screening approach for understanding gene function and functional networks, approximately 100 wild-type plants and 13 known mutants representing a variety of phenotypes were analyzed by a broad range of assays including metabolite profiling, morphological analysis, and chlorophyll fluorescence kinetics. Data analysis using a variety of statistical approaches showed that such industrial approaches can reliably identify plant mutant phenotypes. More significantly, the study uncovered previously unreported phenotypes for these well-characterized mutants and unexpected associations between different physiological processes, demonstrating that this approach has strong advantages over traditional mutant screening approaches. Analysis of wild-type plants revealed hundreds of statistically robust phenotypic correlations, including metabolites that are not known to share direct biosynthetic origins, raising the possibility that these metabolic pathways have closer relationships than is commonly suspected.

A modified LC-MS/MS method for large-scale screening of metabolic phenotypes was developed and validated. Twenty amino acids and 5 metabolically related compounds are measured within 4 min using multiple reaction monitoring (MRM) transitions selective for each compound. Separation with a short C18 column and rapid gradient using the ion-pairing reagent perfluoroheptanoic acid allows chromatographic resolution of the isomers Ile and Leu and improved chromatographic peak shapes for Lys, Arg, and His. MRM transitions were established with capability to distinguish isomers Leu from Ile and Thr from homoserine even when chromatographic resolution is incomplete. The reproducibility of the assay was tested by adding eight stable isotope-labeled amino acid standards (AA*) to extracts of Arabidopsis thaliana seeds. In intra- and interday assay comparisons, mean coefficients of variation of the peak area ratios of AA/AA* were less than 5% for all but Gly/Gly-D2-15N1. Recoveries of these eight amino acids ranged from 62 to 94% and suppression of signal by the matrix were 31-65%. Dilution of seed extracts reduced ion suppression for early-eluting amino acids but had minimal effects for those eluting later. The intra- and interday accuracies for these eight amino acids were 77-131 and 88-133% of nominal concentrations, respectively.

Glucosinolates are secondary metabolites involved in pathogen and insect defense of cruciferous plants. Although seeds and vegetative tissue often have very different glucosinolate profiles, few genetic factors that determine seed glucosinolate accumulation have been identified. An HPLC-based screen of 5500 mutagenized Arabidopsis thaliana lines produced 33 glucosinolate mutants, of which 21 have seed-specific changes. Five of these mutant lines, representing three genetic loci, are compromised in the biosynthesis of benzoyloxyglucosinolates, which are only found in seeds and young seedlings of A. thaliana. Genetic mapping and analysis of T-DNA insertions in candidate genes identified BZO1 (At1g65880), which encodes an enzyme with benzoyl-CoA ligase activity, as being required for the accumulation of benzoyloxyglucosinolates. Long-chain aliphatic glucosinolates are elevated in bzo1 mutants, suggesting substrate competition for the common short-chain aliphatic glucosinolate precursors. Whereas bzo1 mutations have seed-specific effects on benzoyloxyglucosinolate accumulation, the relative abundance of 3-benzoyloxypropyl- and 4-benzoyloxybutylglucosinolates depends on the maternal genotype.