Evolutionary metabolomics of specialized metabolism diversification in the genus Nicotiana highlights N-acylnornicotine innovations.

Specialized metabolite (SM) diversification is a core process to plants' adaptation to diverse ecological niches. Here, we implemented a computational mass spectrometry-based metabolomics approach to exploring SM diversification in tissues of 20 species covering Nicotiana phylogenetics sections. To markedly increase metabolite annotation, we created a large in silico fragmentation database, comprising >1 million structures, and scripts for connecting class prediction to consensus substructures. Together, the approach provides an unprecedented cartography of SM diversity and section-specific innovations in this genus. As a case study and in combination with nuclear magnetic resonance and mass spectrometry imaging, we explored the distribution of N-acylnornicotines, alkaloids predicted to be specific to Repandae allopolyploids, and revealed their prevalence in the genus, albeit at much lower magnitude, as well as a greater structural diversity than previously thought. Together, the data integration approaches provided here should act as a resource for future research in plant SM evolution.

TOR acts as a metabolic gatekeeper for auxin-dependent lateral root initiation in Arabidopsis thaliana.

Plant organogenesis requires matching the available metabolic resources to developmental programs. In Arabidopsis, the root system is determined by primary root-derived lateral roots (LRs), and adventitious roots (ARs) formed from non-root organs. Lateral root formation entails the auxin-dependent activation of transcription factors ARF7, ARF19, and LBD16. Adventitious root formation relies on LBD16 activation by auxin and WOX11. The allocation of shoot-derived sugar to the roots influences branching, but how its availability is sensed for LRs formation remains unknown. We combine metabolic profiling with cell-specific interference to show that LRs switch to glycolysis and consume carbohydrates. The target-of-rapamycin (TOR) kinase is activated in the lateral root domain. Interfering with TOR kinase blocks LR initiation while promoting AR formation. The target-of-rapamycin inhibition marginally affects the auxin-induced transcriptional response of the pericycle but attenuates the translation of ARF19, ARF7, and LBD16. TOR inhibition induces WOX11 transcription in these cells, yet no root branching occurs as TOR controls LBD16 translation. TOR is a central gatekeeper for root branching that integrates local auxin-dependent pathways with systemic metabolic signals, modulating the translation of auxin-induced genes.

BAHD Company: The Ever-Expanding Roles of the BAHD Acyltransferase Gene Family in Plants.

Plants' ability to chemically modify core structures of specialized metabolites is the main reason why the plant kingdom contains such a wide and rich array of diverse compounds. One of the most important types of chemical modifications of small molecules is the addition of an acyl moiety to produce esters and amides. Large-scale phylogenomics analyses have shown that the enzymes that perform acyl transfer reactions on the myriad small molecules synthesized by plants belong to only a few gene families. This review is focused on describing the biochemistry, evolutionary origins, and chemical ecology implications of one of these families-the BAHD acyltransferases. The growth of advanced metabolomic studies coupled with next-generation sequencing of diverse plant species has confirmed that the BAHD family plays critical roles in modifying nearly all known classes of specialized metabolites. The current and future outlook for research on BAHDs includes expanding their roles in synthetic biology and metabolic engineering.

Synthesis of Non-natural Aza-Iridoids via Ynamides and Molecular Networking-Based Tracing of Their In Planta Bioconversion.

A synthesis of new-to-nature aza-iridoids via ynamides is presented. ZrCl4 proved to be the best acid to perform this transformation. Various ynamides were accommodated, and seco-iridoids could be obtained as well. Aza-iridoids were infiltrated into leaves of Scrophularia Nodosa, an iridoid-producing plant species. High-resolution mass spectrometry coupled to computational metabolomic approaches was employed for the detection of aza-iridoid bioconversion products.

Natural history-guided omics reveals plant defensive chemistry against leafhopper pests.

Although much is known about plant traits that function in nonhost resistance against pathogens, little is known about nonhost resistance against herbivores, despite its agricultural importance. Empoasca leafhoppers, serious agricultural pests, identify host plants by eavesdropping on unknown outputs of jasmonate (JA)-mediated signaling. Forward- and reverse-genetics lines of a native tobacco plant were screened in native habitats with native herbivores using high-throughput genomic, transcriptomic, and metabolomic tools to reveal an Empoasca-elicited JA-JAZi module. This module induces an uncharacterized caffeoylputrescine-green leaf volatile compound, catalyzed by a polyphenol oxidase in a Michael addition reaction, which we reconstitute in vitro; engineer in crop plants, where it requires a berberine bridge enzyme-like 2 (BBL2) for its synthesis; and show that it confers resistance to leafhoppers. Natural history-guided forward genetics reveals a conserved nonhost resistance mechanism useful for crop protection.

Evolutionary footprints of a cold relic in a rapidly warming world.

With accelerating global warming, understanding the evolutionary dynamics of plant adaptation to environmental change is increasingly urgent. Here, we reveal the enigmatic history of the genus Cochlearia (Brassicaceae), a Pleistocene relic that originated from a drought-adapted Mediterranean sister genus during the Miocene. Cochlearia rapidly diversified and adapted to circum-Arctic regions and other cold-characterized habitat types during the Pleistocene. This sudden change in ecological preferences was accompanied by a highly complex, reticulate polyploid evolution, which was apparently triggered by the impact of repeated Pleistocene glaciation cycles. Our results illustrate that two early diversified Arctic-alpine diploid gene pools contributed differently to the evolution of this young polyploid genus now captured in a cold-adapted niche. Metabolomics revealed central carbon metabolism responses to cold in diverse species and ecotypes, likely due to continuous connections to cold habitats that may have facilitated widespread adaptation to alpine and subalpine habitats, and which we speculate were coopted from existing drought adaptations. Given the growing scientific interest in the adaptive evolution of temperature-related traits, our results provide much-needed taxonomic and phylogenomic resolution of a model system as well as first insights into the origins of its adaptation to cold.

Toward a Systematic Nomenclature for (Neo)Lignanamides.

Phenylalkenoic acid amides, often referred to as phenol amides or hydroxycinnamic acid amides, are bioactive phytochemicals, whose bioactivity can be enhanced by coupling to form dimers or oligomers. Phenylalkenoic acid amides consist of a (hydroxy)cinnamic acid derivative (i.e., the phenylalkenoic acid subunit) linked to an amine-containing compound (i.e., the amine subunit) via an amide bond. The phenylalkenoic acid moiety can undergo oxidative coupling, either catalyzed by oxidative enzymes or due to autoxidation, which leads to the formation of (neo)lignanamides. Dimers described in the literature are often named after the species in which the compound was first discovered; however, the naming of these compounds lacks a systematic approach. We propose a new nomenclature, inspired by the existing system used for hydroxycinnamic acid dimers and lignin. In the proposed systematic nomenclature for (neo)lignanamides, compound names will be composed of three-letter codes and prefixes denoting the subunits, and numbers that indicate the carbon atoms involved in the linkage between the monomeric precursors. The proposed nomenclature is consistent, future-proof, and systematic.

Next-Generation Mass Spectrometry Metabolomics Revives the Functional Analysis of Plant Metabolic Diversity.

The remarkable diversity of specialized metabolites produced by plants has inspired several decades of research and nucleated a long list of theories to guide empirical ecological studies. However, analytical constraints and the lack of untargeted processing workflows have long precluded comprehensive metabolite profiling and, consequently, the collection of the critical currencies to test theory predictions for the ecological functions of plant metabolic diversity. Developments in mass spectrometry (MS) metabolomics have revolutionized the large-scale inventory and annotation of chemicals from biospecimens. Hence, the next generation of MS metabolomics propelled by new bioinformatics developments provides a long-awaited framework to revisit metabolism-centered ecological questions, much like the advances in next-generation sequencing of the last two decades impacted all research horizons in genomics. Here, we review advances in plant (computational) metabolomics to foster hypothesis formulation from complex metabolome data. Additionally, we reflect on how next-generation metabolomics could reinvigorate the testing of long-standing theories on plant metabolic diversity.

Specific decorations of 17-hydroxygeranyllinalool diterpene glycosides solve the autotoxicity problem of chemical defense in Nicotiana attenuata.

The native diploid tobacco Nicotiana attenuata produces abundant, potent anti-herbivore defense metabolites known as 17-hydroxygeranyllinalool diterpene glycosides (HGL-DTGs) whose glycosylation and malonylation biosynthetic steps are regulated by jasmonate signaling. To characterize the biosynthetic pathway of HGL-DTGs, we conducted a genome-wide analysis of uridine diphosphate glycosyltransferases (UGTs) and identified 107 family-1 UGT members. The transcript levels of three UGTs were highly correlated with the transcript levels two key HGL-DTG biosynthetic genes: geranylgeranyl diphosphate synthase (NaGGPPS) and geranyllinalool synthase (NaGLS). NaGLS's role in HGL-DTG biosynthesis was confirmed by virus-induced gene silencing. Silencing the Uridine diphosphate (UDP)-rhamnosyltransferase gene UGT91T1 demonstrated its role in the rhamnosylation of HGL-DTGs. In vitro enzyme assays revealed that UGT74P3 and UGT74P4 use UDP-glucose for the glucosylation of 17-hydroxygeranyllinalool (17-HGL) to lyciumoside I. Plants with stable silencing of UGT74P3 and UGT74P5 were severely developmentally deformed, pointing to a phytotoxic effect of the aglycone. The application of synthetic 17-HGL and silencing of the UGTs in HGL-DTG-free plants confirmed this phytotoxic effect. Feeding assays with tobacco hornworm (Manduca sexta) larvae revealed the defensive functions of the glucosylation and rhamnosylation steps in HGL-DTG biosynthesis. Glucosylation of 17-HGL is therefore a critical step that contributes to the resulting metabolites' defensive function and solves the autotoxicity problem of this potent chemical defense.

Ethylene is a local modulator of jasmonate-dependent phenolamide accumulation during Manduca sexta herbivory in Nicotiana attenuata.

Rapid reconfigurations of interconnected phytohormone signalling networks allow plants to tune their physiology to constantly varying ecological conditions. During insect herbivory, most of the induced changes in defence-related leaf metabolites are controlled by jasmonate (JA) signalling, which, in the wild tobacco Nicotiana attenuata, recruits MYB8, a transcription factor controlling the accumulation of phenolic-polyamine conjugates (phenolamides). In this and other plant species, herbivory also locally triggers ethylene (ET) production but the outcome of the JA-ET cross-talk at the level of secondary metabolism regulation has remained only superficially investigated. Here, we analysed local and systemic herbivory-induced changes by mass spectrometry-based metabolomics in leaves of transgenic plants impaired in JA, ET and MYB8 signalling. Parsing deregulations in this factorial data-set identified a network of JA/MYB8-dependent phenolamides for which impairment of ET signalling attenuated their accumulation only in locally damaged leaves. Further experiments revealed that ET, albeit biochemically interrelated to polyamine metabolism via the intermediate S-adenosylmethionine, does not alter the free polyamine levels, but instead significantly modulates phenolamide levels with marginal modulations of transcript levels. The work identifies ET as a local modulator of phenolamide accumulations and provides a metabolomics data-platform with which to mine associations among herbivory-induced signalling and specialized metabolites in N. attenuata.

TOC1 in Nicotiana attenuata regulates efficient allocation of nitrogen to defense metabolites under herbivory stress.

The circadian clock contextualizes plant responses to environmental signals. Plants use temporal information to respond to herbivory, but many of the functional roles of circadian clock components in these responses, and their contribution to fitness, remain unknown. We investigate the role of the central clock regulator TIMING OF CAB EXPRESSION 1 (TOC1) in Nicotiana attenuata's defense responses to the specialist herbivore Manduca sexta under both field and glasshouse conditions. We utilize 15 N pulse-labeling to quantify nitrogen incorporation into pools of three defense compounds: caffeoylputrescine (CP), dicaffeoyl spermidine (DCS) and nicotine. Nitrogen incorporation was decreased in CP and DCS and increased in nicotine pools in irTOC1 plants compared to empty vector (EV) under control conditions, but these differences were abolished after simulated herbivory. Differences between EV and irTOC1 plants in nicotine, but not phenolamide production, were abolished by treatment with the ethylene agonist 1-methylcyclopropene. Using micrografting, TOC1's effect on nicotine was isolated to the root and did not affect the fitness of heterografts under field conditions. These results suggest that the circadian clock contributes to plant fitness by balancing production of metabolically expensive nitrogen-rich defense compounds and mediating the allocation of resources between vegetative biomass and reproduction.

Information theory tests critical predictions of plant defense theory for specialized metabolism.

Different plant defense theories have provided important theoretical guidance in explaining patterns in plant specialized metabolism, but their critical predictions remain to be tested. Here, we systematically explored the metabolomes of Nicotiana attenuata, from single plants to populations, as well as of closely related species, using unbiased tandem mass spectrometry (MS/MS) analyses and processed the abundances of compound spectrum-based MS features within an information theory framework to test critical predictions of optimal defense (OD) and moving target (MT) theories. Information components of plant metabolomes were consistent with the OD theory but contradicted the main prediction of the MT theory for herbivory-induced dynamics of metabolome compositions. From micro- to macroevolutionary scales, jasmonate signaling was confirmed as the master determinant of OD, while ethylene signaling provided fine-tuning for herbivore-specific responses annotated via MS/MS molecular networks.

Annotation of Specialized Metabolites from High-Throughput and High-Resolution Mass Spectrometry Metabolomics.

High-throughput mass spectrometry (MS) metabolomics profiling of highly complex samples allows the comprehensive detection of hundreds to thousands of metabolites under a given condition and point in time and produces information-rich data sets on known and unknown metabolites. One of the main challenges is the identification and annotation of metabolites from these complex data sets since the number of authentic standards available for specialized metabolites is far lower than an account for the number of mass spectral features. Previously, we reported two novel tools, MetNet and MetCirc, for putative annotation and structural prediction on unknown metabolites using known metabolites as baits. MetNet employs differences between m/z values of MS1 features, which correspond to metabolic transformations, and statistical associations, while MetCirc uses MS/MS features as input and calculates similarity scores of aligned spectra between features to guide the annotation of metabolites. Here, we showcase the use of MetNet and MetCirc to putatively annotate metabolites and provide detailed instructions as to how those can be used. While our case studies are from plants, the tools find equal utility in studies on bacterial, fungal, or mammalian xenobiotic samples.

The Active Jasmonate JA-Ile Regulates a Specific Subset of Plant Jasmonate-Mediated Resistance to Herbivores in Nature.

The jasmonate hormones are essential regulators of plant defense against herbivores and include several dozen derivatives of the oxylipin jasmonic acid (JA). Among these, the conjugate jasmonoyl isoleucine (JA-Ile) has been shown to interact directly with the jasmonate co-receptor complex to regulate responses to jasmonate signaling. However, functional studies indicate that some aspects of jasmonate-mediated defense are not regulated by JA-Ile. Thus, it is not clear whether JA-Ile is best characterized as the master jasmonate regulator of defense, or if it regulates more specific aspects. We investigated possible functions of JA-Ile in anti-herbivore resistance of the wild tobacco Nicotiana attenuata, a model system for plant-herbivore interactions. We first analyzed the soluble and volatile secondary metabolomes of irJAR4xirJAR6, asLOX3, and WT plants, as well as an RNAi line targeting the jasmonate co-receptor CORONATINE INSENSITIVE 1 (irCOI1), following a standardized herbivory treatment. irJAR4xirJAR6 were the most similar to WT plants, having a ca. 60% overlap in differentially regulated metabolites with either asLOX3 or irCOI1. In contrast, while at least 25 volatiles differed between irCOI1 or asLOX3 and WT plants, there were few or no differences in herbivore-induced volatile emission between irJAR4xirJAR6 and WT plants, in glasshouse- or field-collected samples. We then measured the susceptibility of jasmonate-deficient vs. JA-Ile-deficient plants in nature, in comparison to wild-type (WT) controls, and found that JA-Ile-deficient plants (irJAR4xirJAR6) are much better defended even than a mildly jasmonate-deficient line (asLOX3). The differences among lines could be attributed to differences in damage from specific herbivores, which appeared to prefer either one or the other jasmonate-deficient phenotype. We further investigated the elicitation of one herbivore-induced volatile known to be jasmonate-regulated and to mediate resistance to herbivores: (E)-α-bergamotene. We found that JA was a more potent elicitor of (E)-α-bergamotene emission than was JA-Ile, and when treated with JA, irJAR4xirJAR6 plants emitted 20- to 40-fold as much (E)-α-bergamotene than WT. We conclude that JA-Ile regulates specific aspects of herbivore resistance in N. attenuata. This specificity may allow plants flexibility in their responses to herbivores and in managing trade-offs between resistance, vs. growth and reproduction, over the course of ontogeny.

Current Challenges in Plant Eco-Metabolomics.

The relatively new research discipline of Eco-Metabolomics is the application of metabolomics techniques to ecology with the aim to characterise biochemical interactions of organisms across different spatial and temporal scales. Metabolomics is an untargeted biochemical approach to measure many thousands of metabolites in different species, including plants and animals. Changes in metabolite concentrations can provide mechanistic evidence for biochemical processes that are relevant at ecological scales. These include physiological, phenotypic and morphological responses of plants and communities to environmental changes and also interactions with other organisms. Traditionally, research in biochemistry and ecology comes from two different directions and is performed at distinct spatiotemporal scales. Biochemical studies most often focus on intrinsic processes in individuals at physiological and cellular scales. Generally, they take a bottom-up approach scaling up cellular processes from spatiotemporally fine to coarser scales. Ecological studies usually focus on extrinsic processes acting upon organisms at population and community scales and typically study top-down and bottom-up processes in combination. Eco-Metabolomics is a transdisciplinary research discipline that links biochemistry and ecology and connects the distinct spatiotemporal scales. In this review, we focus on approaches to study chemical and biochemical interactions of plants at various ecological levels, mainly plant⁻organismal interactions, and discuss related examples from other domains. We present recent developments and highlight advancements in Eco-Metabolomics over the last decade from various angles. We further address the five key challenges: (1) complex experimental designs and large variation of metabolite profiles; (2) feature extraction; (3) metabolite identification; (4) statistical analyses; and (5) bioinformatics software tools and workflows. The presented solutions to these challenges will advance connecting the distinct spatiotemporal scales and bridging biochemistry and ecology.

Wild tobacco genomes reveal the evolution of nicotine biosynthesis.

Nicotine, the signature alkaloid of Nicotiana species responsible for the addictive properties of human tobacco smoking, functions as a defensive neurotoxin against attacking herbivores. However, the evolution of the genetic features that contributed to the assembly of the nicotine biosynthetic pathway remains unknown. We sequenced and assembled genomes of two wild tobaccos, Nicotiana attenuata (2.5 Gb) and Nicotiana obtusifolia (1.5 Gb), two ecological models for investigating adaptive traits in nature. We show that after the Solanaceae whole-genome triplication event, a repertoire of rapidly expanding transposable elements (TEs) bloated these Nicotiana genomes, promoted expression divergences among duplicated genes, and contributed to the evolution of herbivory-induced signaling and defenses, including nicotine biosynthesis. The biosynthetic machinery that allows for nicotine synthesis in the roots evolved from the stepwise duplications of two ancient primary metabolic pathways: the polyamine and nicotinamide adenine dinucleotide (NAD) pathways. In contrast to the duplication of the polyamine pathway that is shared among several solanaceous genera producing polyamine-derived tropane alkaloids, we found that lineage-specific duplications within the NAD pathway and the evolution of root-specific expression of the duplicated Solanaceae-specific ethylene response factor that activates the expression of all nicotine biosynthetic genes resulted in the innovative and efficient production of nicotine in the genus Nicotiana Transcription factor binding motifs derived from TEs may have contributed to the coexpression of nicotine biosynthetic pathway genes and coordinated the metabolic flux. Together, these results provide evidence that TEs and gene duplications facilitated the emergence of a key metabolic innovation relevant to plant fitness.

MetCirc: navigating mass spectral similarity in high-resolution MS/MS metabolomics data.

SUMMARY: Among the main challenges in metabolomics are the rapid dereplication of previously characterized metabolites across a range of biological samples and the structural prediction of unknowns from MS/MS data. Here, we developed MetCirc to comprehensively align and calculate pairwise similarity scores among MS/MS spectral data and visualize these across a range of biological samples. MetCirc comprises functionalities to interactively organize these data according to compound familial groupings and to accelerate the discovery of shared metabolites and hypothesis formulation for unknowns. As such, MetCirc provides a significant advance to address biological questions in areas where chemodiversity plays a role. AVAILABILITY AND IMPLEMENTATION: MetCirc , implemented in the open-source R language, together with its vignette are available in the Bioconductor project and at https://github.com/PlantDefenseMetabolism/MetCirc . CONTACT: thomasnaake@googlemail.com or emmanuel.gaquerel@cos.uni-heidelberg.de. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

O-Acyl Sugars Protect a Wild Tobacco from Both Native Fungal Pathogens and a Specialist Herbivore.

O-Acyl sugars (O-AS) are abundant trichome-specific metabolites that function as indirect defenses against herbivores of the wild tobacco Nicotiana attenuata; whether they also function as generalized direct defenses against herbivores and pathogens remains unknown. We characterized natural variation in O-AS among 26 accessions and examined their influence on two native fungal pathogens, Fusarium brachygibbosum U4 and Alternaria sp. U10, and the specialist herbivore Manduca sexta At least 15 different O-AS structures belonging to three classes were found in N. attenuata leaves. A 3-fold quantitative variation in total leaf O-AS was found among the natural accessions. Experiments with natural accessions and crosses between high- and low-O-AS accessions revealed that total O-AS levels were associated with resistance against herbivores and pathogens. Removing O-AS from the leaf surface increased M. sexta growth rate and plant fungal susceptibility. O-AS supplementation in artificial diets and germination medium reduced M. sexta growth and fungal spore germination, respectively. Finally, silencing the expression of a putative branched-chain α-ketoacid dehydrogenase E1 ß-subunit-encoding gene (NaBCKDE1B) in the trichomes reduced total leaf O-AS by 20% to 30% and increased susceptibility to Fusarium pathogens. We conclude that O-AS function as direct defenses to protect plants from attack by both native pathogenic fungi and a specialist herbivore and infer that their diversification is likely shaped by the functional interactions among these biotic stresses.

Nicotiana attenuata Data Hub (NaDH): an integrative platform for exploring genomic, transcriptomic and metabolomic data in wild tobacco.

BACKGROUND: Nicotiana attenuata (coyote tobacco) is an ecological model for studying plant-environment interactions and plant gene function under real-world conditions. During the last decade, large amounts of genomic, transcriptomic and metabolomic data have been generated with this plant which has provided new insights into how native plants interact with herbivores, pollinators and microbes. However, an integrative and open access platform that allows for the efficient mining of these -omics data remained unavailable until now. DESCRIPTION: We present the Nicotiana attenuata Data Hub (NaDH) as a centralized platform for integrating and visualizing genomic, phylogenomic, transcriptomic and metabolomic data in N. attenuata. The NaDH currently hosts collections of predicted protein coding sequences of 11 plant species, including two recently sequenced Nicotiana species, and their functional annotations, 222 microarray datasets from 10 different experiments, a transcriptomic atlas based on 20 RNA-seq expression profiles and a metabolomic atlas based on 895 metabolite spectra analyzed by mass spectrometry. We implemented several visualization tools, including a modified version of the Electronic Fluorescent Pictograph (eFP) browser, co-expression networks and the Interactive Tree Of Life (iTOL) for studying gene expression divergence among duplicated homologous. In addition, the NaDH allows researchers to query phylogenetic trees of 16,305 gene families and provides tools for analyzing their evolutionary history. Furthermore, we also implemented tools to identify co-expressed genes and metabolites, which can be used for predicting the functions of genes. Using the transcription factor NaMYB8 as an example, we illustrate that the tools and data in NaDH can facilitate identification of candidate genes involved in the biosynthesis of specialized metabolites. CONCLUSION: The NaDH provides interactive visualization and data analysis tools that integrate the expression and evolutionary history of genes in Nicotiana, which can facilitate rapid gene discovery and comparative genomic analysis. Because N. attenuata shares many genome-wide features with other Nicotiana species including cultivated tobacco, and hence NaDH can be a resource for exploring the function and evolution of genes in Nicotiana species in general. The NaDH can be accessed at: http://nadh.ice.mpg.de/ .