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.

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.

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.

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.

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

Illuminating a plant's tissue-specific metabolic diversity using computational metabolomics and information theory.

Secondary metabolite diversity is considered an important fitness determinant for plants' biotic and abiotic interactions in nature. This diversity can be examined in two dimensions. The first one considers metabolite diversity across plant species. A second way of looking at this diversity is by considering the tissue-specific localization of pathways underlying secondary metabolism within a plant. Although these cross-tissue metabolite variations are increasingly regarded as important readouts of tissue-level gene function and regulatory processes, they have rarely been comprehensively explored by nontargeted metabolomics. As such, important questions have remained superficially addressed. For instance, which tissues exhibit prevalent signatures of metabolic specialization? Reciprocally, which metabolites contribute most to this tissue specialization in contrast to those metabolites exhibiting housekeeping characteristics? Here, we explore tissue-level metabolic specialization in Nicotiana attenuata, an ecological model with rich secondary metabolism, by combining tissue-wide nontargeted mass spectral data acquisition, information theory analysis, and tandem MS (MS/MS) molecular networks. This analysis was conducted for two different methanolic extracts of 14 tissues and deconvoluted 895 nonredundant MS/MS spectra. Using information theory analysis, anthers were found to harbor the most specialized metabolome, and most unique metabolites of anthers and other tissues were annotated through MS/MS molecular networks. Tissue-metabolite association maps were used to predict tissue-specific gene functions. Predictions for the function of two UDP-glycosyltransferases in flavonoid metabolism were confirmed by virus-induced gene silencing. The present workflow allows biologists to amortize the vast amount of data produced by modern MS instrumentation in their quest to understand gene function.

Beyond the Canon: Within-Plant and Population-Level Heterogeneity in Jasmonate Signaling Engaged by Plant-Insect Interactions.

Plants have evolved sophisticated communication and defense systems with which they interact with insects. Jasmonates are synthesized from the oxylipin pathway and act as pivotal cellular orchestrators of many of the metabolic and physiological processes that mediate these interactions. Many of these jasmonate-dependent responses are tissue-specific and translate from modulations of the canonical jasmonate signaling pathway. Here we provide a short overview of within-plant heterogeneities in jasmonate signaling and dependent responses in the context of plant-insect interactions as illuminated by examples from recent work with the ecological model, Nicotiana attenuata. We then discuss means of manipulating jasmonate signaling by creating tissue-specific jasmonate sinks, and the micrografting of different transgenic plants. The metabolic phenotyping of these manipulations provides an integrative understanding of the functional significance of deviations from the canonical model of this hormonal pathway. Additionally, natural variation in jasmonate biosynthesis and signaling both among and within species can explain polymorphisms in resistance to insects in nature. In this respect, insect-guided explorations of population-level variations in jasmonate metabolism have revealed more complexity than previously realized and we discuss how different "omic" techniques can be used to exploit the natural variation that occurs in this important signaling pathway.

Using the knowns to discover the unknowns: MS-based dereplication uncovers structural diversity in 17-hydroxygeranyllinalool diterpene glycoside production in the Solanaceae.

Exploring the diversity of plant secondary metabolism requires efficient methods to obtain sufficient structural insights to discriminate previously known from unknown metabolites. De novo structure elucidation and confirmation of known metabolites (dereplication) remain a major bottleneck for mass spectrometry-based metabolomic workflows, and few systematic dereplication strategies have been developed for the analysis of entire compound classes across plant families, partly due to the complexity of plant metabolic profiles that complicates cross-species comparisons. 17-hydroxygeranyllinalool diterpene glycosides (HGL-DTGs) are abundant defensive secondary metabolites whose malonyl and glycosyl decorations are induced by jasmonate signaling in the ecological model plant Nicotiana attenuata. The multiple labile glycosidic bonds of HGL-DTGs result in extensive in-source fragmentation (IS-CID) during ionization. To reconstruct these IS-CID clusters from profiling data and identify precursor ions, we applied a deconvolution algorithm and created an MS/MS library from positive-ion spectra of purified HGL-DTGs. From this library, 251 non-redundant fragments were annotated, and a workflow to characterize leaf, flower and fruit extracts of 35 solanaceous species was established. These analyses predicted 105 novel HGL-DTGs that were restricted to Nicotiana, Capsicum and Lycium species. Interestingly, malonylation is a highly conserved step in HGL-DTG metabolism, but is differentially affected by jasmonate signaling among Nicotiana species. This MS-based workflow is readily applicable for cross-species re-identification/annotation of other compound classes with sufficient fragmentation knowledge, and therefore has the potential to support hypotheses regarding secondary metabolism diversification.

Navigating natural variation in herbivory-induced secondary metabolism in coyote tobacco populations using MS/MS structural analysis.

Natural variation can be extremely useful in unraveling the determinants of phenotypic trait evolution but has rarely been analyzed with unbiased metabolic profiling to understand how its effects are organized at the level of biochemical pathways. Native populations of Nicotiana attenuata, a wild tobacco species, have been shown to be highly genetically diverse for traits important for their interactions with insects. To resolve the chemodiversity existing in these populations, we developed a metabolomics and computational pipeline to annotate leaf metabolic responses to Manduca sexta herbivory. We selected seeds from 43 accessions of different populations from the southwestern United States--including the well-characterized Utah 30th generation inbred accession--and grew 183 plants in the glasshouse for standardized herbivory elicitation. Metabolic profiles were generated from elicited leaves of each plant using a high-throughput ultra HPLC (UHPLC)-quadrupole TOFMS (qTOFMS) method, processed to systematically infer covariation patterns among biochemically related metabolites, as well as unknown ones, and finally assembled to map natural variation. Navigating this map revealed metabolic branch-specific variations that surprisingly only partly overlapped with jasmonate accumulation polymorphisms and deviated from canonical jasmonate signaling. Fragmentation analysis via indiscriminant tandem mass spectrometry (idMS/MS) was conducted with 10 accessions that spanned a large proportion of the variance found in the complete accession dataset, and compound spectra were computationally assembled into spectral similarity networks. The biological information captured by this networking approach facilitates the mining of the mass spectral data of unknowns with high natural variation, as demonstrated by the annotation of a strongly herbivory-inducible phenolic derivative, and can guide pathway analysis.

Jasmonoyl-L-isoleucine coordinates metabolic networks required for anthesis and floral attractant emission in wild tobacco (Nicotiana attenuata).

Jasmonic acid and its derivatives (jasmonates [JAs]) play central roles in floral development and maturation. The binding of jasmonoyl-L-isoleucine (JA-Ile) to the F-box of CORONATINE INSENSITIVE1 (COI1) is required for many JA-dependent physiological responses, but its role in anthesis and pollinator attraction traits remains largely unexplored. Here, we used the wild tobacco Nicotiana attenuata, which develops sympetalous flowers with complex pollination biology, to examine the coordinating function of JA homeostasis in the distinct metabolic processes that underlie flower maturation, opening, and advertisement to pollinators. From combined transcriptomic, targeted metabolic, and allometric analyses of transgenic N. attenuata plants for which signaling deficiencies were complemented with methyl jasmonate, JA-Ile, and its functional homolog, coronatine (COR), we demonstrate that (1) JA-Ile/COR-based signaling regulates corolla limb opening and a JA-negative feedback loop; (2) production of floral volatiles (night emissions of benzylacetone) and nectar requires JA-Ile/COR perception through COI1; and (3) limb expansion involves JA-Ile-induced changes in limb fresh mass and carbohydrate metabolism. These findings demonstrate a master regulatory function of the JA-Ile/COI1 duet for the main function of a sympetalous corolla, that of advertising for and rewarding pollinator services. Flower opening, by contrast, requires JA-Ile signaling-dependent changes in primary metabolism, which are not compromised in the COI1-silenced RNA interference line used in this study.

Revealing insect herbivory-induced phenolamide metabolism: from single genes to metabolic network plasticity analysis.

The phenylpropanoid metabolic space comprises a network of interconnected metabolic branches that contribute to the biosynthesis of a large array of compounds with functions in plant development and stress adaptation. During biotic challenges, such as insect attack, a major rewiring of gene networks associated with phenylpropanoid metabolism is observed. This rapid reconfiguration of gene expression allows prioritized production of metabolites that help the plant solve ecological problems. Phenolamides are a group of phenolic derivatives that originate from diversion of hydroxycinnamoyl acids from the main phenylpropanoid pathway after N-acyltransferase-dependent conjugation to polyamines or aryl monoamines. These structurally diverse metabolites are abundant in the reproductive organs of many plants, and have recently been shown to play roles as induced defenses in vegetative tissues. In the wild tobacco, Nicotiana attenuata, in which herbivory-induced regulation of these metabolites has been studied, rapid elevations of the levels of phenolamides that function as induced defenses result from a multi-hormonal signaling network that re-shapes connected metabolic pathways. In this review, we summarize recent findings in the regulation of phenolamides obtained by mass spectrometry-based metabolomics profiling, and outline a conceptual framework for gene discovery in this pathway. We also introduce a multifactorial approach that is useful in deciphering metabolic pathway reorganizations among tissues in response to stress.

The roots of plant defenses: integrative multivariate analyses uncover dynamic behaviors of gene and metabolic networks of roots elicited by leaf herbivory.

High-throughput analyses have frequently been used to characterize herbivory-induced reconfigurations in plant primary and secondary metabolism in above- and below-ground tissues, but the conclusions drawn from these analyses are often limited by the univariate methods used to analyze the data. Here we use our previously described multivariate time-series data analysis to evaluate leaf herbivory-elicited transcriptional and metabolic dynamics in the roots of Nicotiana attenuata. We observed large, but transient, systemic responses in the roots that contrasted with the pattern of co-linearity observed in the up- and downregulation of genes and metabolites across the entire time series in treated and systemic leaves. Using this newly developed approach for the analysis of whole-plant molecular responses in a time-course multivariate data set, we simultaneously analyzed stress responses in leaves and roots in response to the elicitation of a leaf. We found that transient systemic responses in roots resolved into two principal trends characterized by: (i) an inversion of root-specific semi-diurnal (12 h) transcript oscillations and (ii) transcriptional changes with major amplitude effects that translated into a distinct suite of root-specific secondary metabolites (e.g. alkaloids synthesized in the roots of N. attenuata). These findings underscore the importance of understanding tissue-specific stress responses in the correct day-night phase context and provide a holistic framework for the important role played by roots in above-ground stress responses.

An integrative statistical method to explore herbivory-specific responses in plants.

Spatial-temporal coordination between multiple processes/pathways is a key determinant of whole-organism transcriptome and metabolome reconfigurations in plants' response to biotic stresses. To explore tissue-based interdependencies in Nicotiana attenuata's resistance to insect attack, we performed time course analyses of the plant's transcriptome and metabolome in herbivory-elicited source leaves and unelicited sink leaves and roots. To dissect the multidimensionality of these responses, we have recently designed a novel approach of constructing interactive motifs by combining an extended self-organizing maps (SOM) based dimensionality reduction method with bootstrap-based non-parametric AN OVA models. In this previous study, we used this method to study nonlinearities in gene-metabolite associations involved in the acyclic diterpene glucoside pathway. Here, we extend the application of this method to the extraction of genes showing herbivory-elicitation specifically in systemic (distal from the treatment sites) tissues using motif analysis for different combinations of treatment applied to Nicotiana attenuata.