# Unravelling Plant Responses to Stress-The Importance of Targeted and Untargeted Metabolomics

James William Allwood, Alex Williams, Henriette Uthe, Nicole M van Dam, Luis a J Mur, Murray R Grant, Pierre Pétriacq

► To cite this version:

James William Allwood, Alex Williams, Henriette Uthe, Nicole M van Dam, Luis a J Mur, et al. Unravelling Plant Responses to Stress-The Importance of Targeted and Untargeted Metabolomics. Metabolites, 2021, 11 (8), pp.558. 10.3390/metabo11080558. [https://doi.org/10.3390/metabo11080558](https://doi.org/10.3390/metabo11080558)

HAL Id: hal-03329826

https://hal.inrae.fr/hal-03329826

Submitted on 31 Aug 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Review

# Unravelling Plant Responses to Stress—The Importance of Targeted and Untargeted Metabolomics

James William Allwood <fnref n="1" />,* , Alex Williams <fnref n="2" />,<fnref n="3" /> , Henriette Uthe <fnref n="4" /> , Nicole M. van Dam <fnref n="4" /> ID , Luis A. J. Mur <fnref n="5" /> ID , Murray R. Grant <fnref n="6" /> and Pierre Pétriacq <fnref n="7" />,<fnref n="8" />,* ID

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Allwood, J.W.; Williams, A.; Uthe, H.; van Dam, N.M.; Mur, L.A.J.; Grant, M.R.; Pétriacq, P. Unravelling Plant Responses to Stress—The Importance of Targeted and Untargeted Metabolomics. *Metabolites* **2021**, *11*, 558. https://doi.org/10.3390/metabo11080558

Academic Editors:
Hirokazu Kawagishi and
Peter Meikle

Received: 29 June 2021
Accepted: 16 August 2021
Published: 22 August 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** Climate change and an increasing population, present a massive global challenge with respect to environmentally sustainable nutritious food production. Crop yield enhancements, through breeding, are decreasing, whilst agricultural intensification is constrained by emerging, re-emerging, and endemic pests and pathogens, accounting for ~30% of global crop losses, as well as mounting abiotic stress pressures, due to climate change. Metabolomics approaches have previously contributed to our knowledge within the fields of molecular plant pathology and plant-insect interactions. However, these remain incredibly challenging targets, due to the vast diversity in metabolite volatility and polarity, heterogeneous mixtures of pathogen and plant cells, as well as rapid rates of metabolite turn-over. Unravelling the systematic biochemical responses of plants to various individual and combined stresses, involves monitoring signaling compounds, secondary messengers, phytohormones, and defensive and protective chemicals. This demands both targeted and untargeted metabolomics approaches, as well as a range of enzymatic assays, protein assays, and proteomic and transcriptomic technologies. In this review, we focus upon the technical and biological challenges of measuring the metabolome associated with plant stress. We illustrate the challenges, with relevant examples from bacterial and fungal molecular pathologies, plant-insect interactions, and abiotic and combined stress in the environment. We also discuss future prospects from both the perspective of key innovative metabolomic technologies and their deployment in breeding for stress resistance.

**Keywords:** metabolomics; plant pathology; plant-insect interactions; biotic stress; abiotic stress; targeted; untargeted; systems biology; stress resistance; breeding

1. Introduction

Climate change and population growth pose a major global challenge to environmentally sustainable food production. Current agricultural challenges include the unsustainable carbon footprint of artificial fertilizers and pesticides, with growing public

<footnote n="1">1 Environmental and Biochemical Sciences, James Hutton Institute, Errol Road, Invergowrie,
Dundee DD2 5DA, UK</footnote>

<footnote n="2">2 School of Earth and Environmental Sciences, The University of Manchester, Oxford Road, Manchester M13 9PT, UK; alex.williams-4@sheffield.ac.uk</footnote>

<footnote n="3">3 Department of Animal and Plant Sciences, Biosciences, The University of Sheffield Western Bank, Sheffield S10 2TN, UK</footnote>

<footnote n="4">4 German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Molecular Interaction Ecology Group, Friedrich-Schiller University Jena, Puschstr. 4, 04103 Leipzig, Germany; henriette.uthe@idiv.de (H.U.); nicole.vandam@idiv.de (N.M.v.D.)</footnote>

<footnote n="5">5 Institute of Biological, Environmental and Rural Sciences (IBERS), Edward Llwyd Building,
Aberystwyth University, Aberystwyth SY23 3DA, UK; lum@aber.ac.uk</footnote>

<footnote n="6">6 Gibbet Hill Campus, School of Life Sciences, The University of Warwick, Coventry CV4 7AL, UK;
m.grant@warwick.ac.uk</footnote>

<footnote n="7">7 UMR 1332 Fruit Biology and Pathology, Centre INRAE de Nouvelle Aquitaine Bordeaux,
University of Bordeaux, 33140 Villenave d'Ornon, France</footnote>

<footnote n="8">8 Bordeaux Metabolome, MetaboHUB, PHENOME-EMPHASIS, Centre INRAE de Nouvelle Aquitaine-Bordeaux, 33140 Villenave d'Ornon, France</footnote>

<footnote>* Correspondence: will.allwood@hutton.ac.uk (J.W.A.); pierre.petriacq@inrae.fr (P.P.);
Tel.: +44-(0)1382-568-745 (J.W.A.); +33-(0)5-57-12-25-75 (P.P.)</footnote>
concern about their environmental and health impacts. Superimposed on this, yield enhancements, through crop breeding, are declining (below 1% year-on-year), while agricultural intensification continues to be constrained by emerging, re-emerging, and endemic pests and pathogens, which are responsible for 20–40% of crop losses worldwide (FAO: http://www.fao.org/news/story/en/item/1187738/icode/ (accessed 15 July 2021)). Climate change will further exacerbate these global issues as rising temperatures affect the geographical distribution of food production, change pest and pathogen distribution, and escalate the magnitude and frequency of abiotic stresses. Thus, feeding an increasing population, whilst limiting land use and environmental impact, represents one of the most significant challenges of the 21st century. It thus becomes increasingly important to deploy innovative genetic and metabolic GWAS (genome-wide association studies) approaches, to define agronomically important QTLs (quantitative trait loci), and adopt new breeding technologies to improve elite crops, focusing on robust disease, pest, and abiotic stress-resistant plants, with lower dependence on chemical fertilizers, while maintaining high yield and quality.

Plants face a multitude of abiotic (e.g., drought, salinity, heat, and water stress) and biotic (pathogenic bacteria, fungi, viruses, nematodes, and insect pests) stresses throughout their life cycles. Plants have evolved a remarkable range of mechanisms to respond rapidly to environmental changes. They can sense the presence of pests, parasites, and pathogens, through evolving effective molecular recognition mechanisms and the capacity to respond to physical changes [1]. The responses to pests and pathogens are complex, involving host-microbe/pest communication and distinct defensive mechanisms. Aside from the physical barriers, innate immune responses, initiated by tissue damage or recognition of microbe/pathogen/nematode-associated molecular patterns—DAMPs, MAMPs, PAMPs, and NAMPs (e.g., lipopolysaccharide—LPS, flagellin, and nematode specific pheromones, i.e., ascarosides)—by host membrane-localized pattern recognition receptors (PRRs), leads to microbe/pathogen-triggered immunity (M/PTI) [2,3]. Adapted pathogens have evolved virulence strategies, deploying proteinaceous or small-molecule “effectors” (virulence factors), which collectively act to suppress PTI and mobilize host nutrients for pathogen nutrition. In some cases, plants have evolved disease resistance (R) proteins that recognize effectors (avirulence factors), or their activity, to initiate a hypersensitive response (HR) that contains the invading pathogen, known as effector-triggered immunity (ETI) [3,4]. ETI also induces systemic-acquired resistance (SAR), a priming response in the whole plant that stimulates long-lasting immunity to a diverse range of pathogens [5–7].

The initial stress response of plants is through an integrated network of signaling molecules, inclusive of a wide range of secondary messengers (calcium, reactive oxygen species (ROS), nitric oxide (NO), and lipids) and phytohormones (e.g., salicylic acid (SA), jasmonic acid (JA), abscisic acid (ABA), and ethylene) [7–11]. The network is primarily mediated by mitogen and calcium-activated protein kinase (MAPK and CAPK) signal transduction cascades [12–14]. Understanding the timing, amplitude, and duration of signaling events during plant stress, as well as the interplay between phytohormones and transcriptional re-programming, leading to physical and chemical defenses, is a massive, but critically important, challenge within plant biology. Additionally, different subcellular organelles are known to play unique and pivotal roles in integrating biotic and abiotic stress signals, with, for example, chloroplasts orchestrating the production of a range of phytohormones and secondary metabolites [7,15]. This highlights the future importance of sub-cellular metabolomic analyses.

Over the past two decades, the emergence of metabolomics as a systems biology strategy has dramatically improved our knowledge of complex biological systems. Similarly to other “omics” sciences, metabolomics has made a prominent contribution to our understanding of plant immune responses [16]. A typical experimental design (Figure 1) includes targeted and untargeted metabolomics strategies, depending on the objective of the study, as follows: (i) the quantification of specific metabolites, due to their importance during the stress response, or (ii) the acquisition of global metabolic profiles from which a
search for defense response/resistance metabolic pathways, novel signaling molecules, or biomarkers can be conducted (Figure 2). In the face of microbial challenges, metabolomics has revealed specific and common metabolic signatures as hallmarks of infection processes. These include antimicrobial metabolite counterattacks or metabolic signals, to reorganize plant physiology and adapt to stress. Here, we will examine the major metabolomics techniques and approaches that have made, or have the potential to make, significant contributions to our understanding of plant stress responses, with a specific focus on bacterial, fungal, and pest interactions. We will further overview the metabolomics research into biotic stress responses within the natural environment, and discuss some challenges and perspectives of phytopathological metabolomics.

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Experimental design in phytopathological metabolomics. Plant tissues (here, tomato leaves) are inoculated without (i.e., mock) or with the pathogenic microbe (here, a fungus causing necrotic lesions), or exposed to other biotic challenges (e.g., elicitors). Subsequent sampling of control and infected material is usually performed following a time course expressed as day post-inoculation (DPI), sample metabolism is quenched by snap-freezing in liquid nitrogen, followed by storage at -80 °C until metabolomic analyses. Targeted metabolomics will allow measuring the levels of specific compounds (here, the defense hormone salicylate, SA). On the other hand, untargeted approaches can evaluate the global impacts of the microbial interactions on metabolic profiles using multivariate statistics (here, a principal component analysis score plot), where the best discriminators could be revealed (represented as arrows, e.g., SA).</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Typical LC-MS-based metabolomics workflow. Targeted metabolomics (bottom panel) aims to deliver quantitative concentrations of specific compounds using metabolite standards to optimize the extraction and LC-MS acquisition. Then, univariate statistical methods allow investigating the targeted data for statistical significance using adjustment of false discovery rate and descriptive analysis (e.g., adjusted *p*-value, single- and multiple-factor ANOVA), and comparative analysis between the control and infected samples at different time points using single or multiple groups (e.g., *t*-test, Kruskal–Wallis, ANOVA with post hoc tests). When using isotopically labelled standards, targeted metabolomics provide absolute quantification of the metabolite of interest. Alternatively, untargeted metabolomics (top panel) generates a comprehensive metabolic overview of several thousands of metabolic features from a global extraction, which are analyzed by untargeted LC-MS profiling, for instance. The processing of untargeted data relies on sophisticated bioinformatics workflow (including extraction, alignment, grouping and retention time correction) for all the detected LC-MS features. Multivariate statistical methods are subsequently applied to the dataset to visualize how metabolic profiles respond to the conditions (i.e., biotic stress) and unravel stress-responsive discriminant metabolic markers. These relevant markers are putatively annotated using metabolite databases, then confirmed by comparing MS/MS and retention time data with metabolite standards, thereby returning to a targeted approach for verification. Targeted and untargeted metabolomics are complementary in addressing how relevant metabolites and metabolism respond or comprehensively react to stress. Adapted from Williams et al. (2020).</figcaption></figure>
## 2. Plant Stress Responses Are Metabolically Diverse and Require a Suite of Technologies for Accurate Characterization

Chemical defense underpins survival strategies in the plant kingdom. The interactions between plants and their pathogens and herbivores have led to the evolution of an extensive diversity of defense compounds, representing many chemical classes, including terpenes, alkaloids, and glucosinolates. Defense responses in plants are dynamic, in terms of both the concentration and chemical modification of metabolites. In turn, pests and pathogens evolve new strategies to dampen host defense signaling and overcome the toxic effects of defense compounds, including evolving remarkable bioactive mimics, such as the phytotoxin coronatine, produced by certain strains of *Pseudomonas syringae* [17,18]. Signaling pathways, lipids and phytohormones, ROS generation, structural defenses, such as cell wall strengthening (callose deposition) and the deposition of leaf waxes, all add to this rich chemical tapestry and the technically challenging nature of measuring the plant stress response.

### 2.1. Initial Stress Responses Involve Complex Protein Kinase Signalling Networks and the Generation of Phytohormones

As the precursor to chemical defense, it is of crucial importance to monitor signaling processes in response to stress. MAPK and CAPK cascades are among the first steps downstream of PRR stress/pathogen recognition, and are subsequently deployed during the R protein recognition of intracellular effectors (avirulence factors). These processes are conserved and equally important in response to abiotic stress, such as drought and salinity [19–21]. Phosphorylation of the M/CAPK results in the activation of various effector proteins in the cytoplasm or nucleus, including other kinases, enzymes, cytoskeletal proteins, or transcription factors [19,20]. If the protein kinase target is known, its level of activity can be determined by applying Western blotting [22], whereas if the target is unknown, large-scale protein arrays can be screened for phosphorylation activity, in a high-throughput manner [21]. These targets help infer the signaling processes underpinning the metabolic changes that are induced during the plant immune or stress response.

Phytohormones are hugely significant, not just to plant stress responses, but with respect to all aspects of plant growth and development, inclusive of ABA, SA and derivatives, JA and derivatives, the auxins, brassinosteroids, cytokinins, gibberellins, strigolactones, ethylene, NO, pipecolic (PIP) and n-hydroxy pipecolic (NHP) acids, peptide hormones, and polyamines. Phytohormones serve as local and long-distance signaling molecules, and they can act both antagonistically and synergistically, altering defense-related gene expression and transcription levels, cellular division, and growth, both at the localized site of stress and systemically [23]. The phytohormones are generally low-molecular-weight (MW) compounds that provide activity at very low concentrations, $10^{-6}$ to $10^{-5}$ mol/L, their turnover rates can be extremely rapid, and once extracted, may be readily degraded, thus providing an analytical challenge. The chemical diversity of phytohormones further adds to this challenge, with phytohormones ranging from being highly polar to non-polar lipids, as well as varying in volatility. Certain hormones, such as SA and JA, are non-volatile, but the simple addition of a methyl group (i.e., MeJA and MeSA) confers volatility, allowing them to act as airborne signals, to communicate to neighboring leaves or plants that they are under attack by herbivores or pathogens. Targeted absolute-quantitative methods, applying ultra-high-performance liquid chromatography (UHPLC) and electro-spray ionization (ESI) triple-quadrupole mass spectrometry (QqQ-MS) in conjunction with surrogate-labelled internal standards, which cover a wide range of phytohormones, and from as little as 20 mg of fresh leaf tissue, have been successfully developed [24]. Volatile and semi-volatile phytohormones may also be analyzed via gas chromatography (GC)-MS analysis of methyl-chloroformate derivatives [25]. Ethylene and NO, being extremely low MW gases, require specialist approaches, such as targeted GC-flame ionization detection (FID) [26] or photo-acoustic laser spectrophotometry [27]. Pipecolic and hydroxy pipecolic acids have been determined via GC-MS [28] and liquid chromatography (LC)-MS [29]. Within samples collected at the localized site of infection or herbivore attack, the levels of
phytohormones, such as SA and JA, are induced at levels that are high enough to permit their detection through untargeted LC-ESI-MS or DI/FI (direct infusion/flow infusion)-ESI-MS metabolomics approaches [27,30–32]. However, with respect to the long-term priming of SAR, their concentrations are maintained at much lower levels, and therefore targeted LC-QqQ-MS methods are valuable [24].

## 2.2. Redox Carrier/Activation Signals Play a Key Role in the Homeostasis of Reactive Oxygen Species That Lead to Programmed Cell Death

Redox carrier/activation signals, including various pyridine nucleotides (i.e., nicotinamide adenine dinucleotide (NAD) and adenosine triphosphate (ATP)), are also of much importance in ROS and redox signaling against biotrophs and necrotrophs, and the activation of JA responses [33,34]. Likewise, NO has a pivotal role in regulating ROS [8,35]. ATP, NAD+, NADH, and NADPH have classically been assayed via enzymatic colorimetric assays [34,36], although LC-MS and nuclear magnetic resonance (NMR) methods are also available [37,38]. The methods of NO analysis have included the use of fluorescent probes, the Griess reaction, electron paramagnetic resonance (EPR), chemiluminescence, the oxyhemoglobin (Hb-O2) assay, membrane inlet-MS, amperometric methods, and laser photoacoustics [39]. Other important signaling compounds that modulate the expression of genes that are involved in plant signaling, transcriptional regulation, hormone biosynthesis, ROS generation, and polyamine metabolism, include the following: the 3'-phosphoadenosine 5'-phosphate (PAP) pathway, which interacts with ABA, especially during abiotic stress responses (e.g., stomatal regulation in drought and water stress), and which can be determined via HPLC-fluorescence or LC-MS approaches [40]; the non-proteinogenic amino acids, γ- and β-amino butyric acids (GABA and BABA), which are involved in defense priming responses to pathogens, insects, and abiotic stresses [41,42], and which can be determined via LC-MS, GC-MS, or NMR; and lipid-derived signals that are liberated from the phospho- and galacto-lipids, both enzymatically and non-enzymatically, including linolenic and linoleic acids, 18 carbon fatty acids that serve as precursors to 12-oxo-phytodienoic acid (OPDA), a C18 cyclopentenone that is converted to JA and its derivatives [30,43], as well as phosphatidic acids (PA's), which have been associated with the regulation of both SA's and JA's [44]. LC-MS lipidomics approaches present a powerful tool for the study of lipid signaling in plants. In conjunction with measuring stress-related and redox activation signals, it is also important to consider ROS generation via histochemical staining, or quantitative determination through the luminol chemiluminescence assay, and the levels of cell death via electrolyte leakage assays [34,45]. Genetically encoded ROS reporters are increasingly being deployed, to provide spatial-temporal information on subcellular ROS generation [46]. A range of genetically encoded reporters, with different subcellular locations, capable of measuring ATP [47], monitoring intracellular H₂O₂ dynamics [48], NAD redox status [49], or dynamic changes in NADPH [50], are being developed, and these will be invaluable in providing a spatial context to targeted metabolomics approaches.

## 2.3. Analysis of Central Metabolism, a Range of Defence Compound Chemistries, and Secondary Cell Wall Strengthening, Are Key to Successful Studies of Plant Pathogen or Herbivore-Induced Stress

Signaling events and phytohormonal regulation are, of course, just the tip of the 'plant stress' iceberg. Integrated studies, monitoring gene transcript expression, applying proteomics to monitor enzymatic regulation, as well as metabolomic analysis to study defense chemicals, physiological responses, such as cell wall strengthening, and the production of osmotic protectants, will be key to our systematic understanding of a given plant stress response [34,51,52]. Innate defenses are observed against both biotrophic and necrotrophic pathogens, as well as herbivores, and involve generalized physical defenses. The deposition of long-chain leaf cuticular waxes is a common response, both in the defense against herbivore feeding [53] and egg deposition, and with respect to desiccation resistance in abiotic stress [54]. Secondary cell wall strengthening, by the deposition of callose, serves to restrict the invasive spread of pathogens [34,52], and, in some cases, limits the effects of abiotic stresses, such as salinity [55]. Callose, a polysaccharide consisting of β-1,3-glucan
with some β-1,6-glucan branches, is commonly deposited in cell walls at the site of the pathogen ingress and is often imaged through fluorescence microscopy. Both the structural changes and spatial distribution of callose can be assessed through approaches such as Fourier transform infrared (FT-IR) and Raman microscopy, or MS imaging (MSi), such as matrix-assisted laser desorption/ionization (MALDI)-MS, desorption (D)ESI-MS, laser ablation (LA)ESI-MS, and time-of-flight secondary ionization (ToF-SI)MS [56,57]. Changes within the structure of callose, in response to stress, could also be assessed through multistage LC-MS/MS and MS<sup>n</sup>, as well as 2D NMR post purification. Long-chain cuticular waxes have been classically analyzed with GC-MS [58,59] and, recently, ToF-SIMS has also been applied to the spatial analysis of leaf cuticular wax distributions [60].

Dependent upon the host plant species under investigation, as well as the nature of the biotic stress (and in some examples, abiotic stress), a range of plant secondary metabolite classes, with known toxicity against microbes and/or insects, are present and should be considered for analysis. Classically, glucosinolates and isothiocyanates are important defensive compounds in Brassica responses to herbivory [31,61], as are alkaloids in a range of species, especially the Solanaceae [62,63] and benzoxazinoids, for example DIMBOA (2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one) in cereals such as maize [64]. With respect to both microbial pathogens and insect herbivores, phenolic compounds, such as the cyanogenic glycosides [65–67], flavonoids [68–70], and lignins [71,72], play important roles, as do carotenoids, in relation to stress response signaling and abiotic stress [73], and terpenes, which serve as volatile stress signals that prime neighboring leaves and plants [74,75]. All of these classes are amenable to untargeted metabolomics approaches, with C18 reverse-phase (RP) UHPLC high-resolution (HR)MS, with the exception of the terpenes and carotenoids. Terpenes require volatile amenable methods, such as thermal desorption (TD) or solid-phase micro extraction (SPME) interfaced with GC-MS [74], or, alternatively, proton-transfer reaction (PTR)-MS [75]. Carotenoids require targeted UHPLC-atmospheric pressure chemical ionization (APCI)-MS methods, applying C30 RP columns and high-strength non-polar gradients [76]. Finally, central metabolism, assayed via GC-MS of methoximated trimethylsilyl (MOX-TMS) derivatives, hydrophobic interaction LC (HILIC)-MS, or ¹H-NMR, is also of importance. Pathogens and insects target photosynthate and carbon storage forms, such as starch, and, in response, plant hosts mobilize primary metabolite resources away from the localized site of attack. Central metabolite hubs, such as phenylalanine, lysine, and glutamate, which divert resources to defensive compounds, play significant roles [77]. Moreover, the importance of osmoprotectants, such as a number of sugars, sugar alcohols, and amino acids, are becoming increasingly recognized, especially in abiotic stress [78,79].

This review has, so far, highlighted that studying plant stress metabolism is a technological challenge, requiring an array of extraction, purification and separation methods, and analytical technologies. To capture the full spectrum of key compound chemistries necessitates a range of conventional tissue extracts, volatile capture methods, and spatial imaging approaches. It is also challenging with respect to the complexity of sub-cellular (e.g., chloroplast) orchestrated defenses, and the differentiation of host and pathogen metabolomes [80]. Here, technologies such as cellular fractionation, cell sorting, and differential isotopic labelling, will be of significant value. In the following sections, we appraise the challenges and provide a review of metabolomics approaches in the field of plant stress, as well as discussing how metabolomics approaches can be deployed to improve our systematic understanding of plant stress and breed improved disease- and stress-resistant crop varieties.

### 3. Metabolomics Research in Molecular Plant Pathology—Challenges and Complexity

It is essential to recognize that, despite the wealth of technical approaches that have been described above, significant biological challenges to studying plant-pathogen interactions still persist. For example, defining the timing of in vivo sampling, to capture the critical early stages of pathogenesis, is one such difficulty. Pathogen effectors work
largely in collaboration and in a cell-autonomous manner, hence non-responding tissue
invariably consists of the majority of the sample, attenuating biological signal to noise.
Furthermore, the bioactive molecules underpinning early responses often do not exhibit
large changes in abundance and/or vary dynamically. For example, in any given infected
leaf, the establishment of fungal and oomycete haustoria are rare in the overall context
of total leaf tissue/cell counts. When sampling entire leaves, this would mean that only
the stronger metabolic changes can be elucidated. Though obvious symptoms develop,
and pathogen biomass increases later in the infection process, immunity has already been
suppressed. Besides, haustorial establishment is asynchronous, thus temporal dynamics
are not captured. By contrast, bacterial infections, particularly those that are apoplastically
localized, such as *P. syringae* pv. tomato (*Pst*), provide more tractable systems for both
targeted and untargeted defense metabolomics. Not only can one alter the inoculum
concentration, but synchronous infections across an entire leaf can be established with
high-resolution sampling [81]. Moreover, the availability of type III secretion-deficient
mutants, where virulence effector delivery into the cell is blocked, provides a readout
of the metabolic changes associated with PTI. Thus, bacterial pathosystems can provide
insight into the early metabolite changes captured between a non-pathogenic challenge
and a virulent challenge, which are modulated by effectors. Despite this, the actual number
of infection sites in the whole leaf remains low, even following a high inoculum challenge.
Later in the infection process, as the bacterial load increases, metabolomics approaches
are most likely capturing changes that are associated with primary metabolism repro-
gramming, to facilitate pathogen nutrition. Local and distal communication, imperative
in defining disease outcome, is predominately facilitated by small molecules. Untargeted
metabolomics provides the unparalleled opportunity to uncover new signaling molecules,
but will always be constrained by the great diversity and number of chemical structures
in plants.

Metabolomics approaches have revealed unexpected levels of complexity, with re-
spect to phytohormonal regulation and balance in biotic stress responses. Traditionally, SA
and JA are associated with positively regulating biotrophic and necrotrophic interactions,
respectively, and they are also assumed to be mutually antagonistic. The application of tran-
scriptome informed targeted phytohormone analyses, however, identified a critical early
role for ABA, revealing that biotrophic pathogens hijack ABA, to promote virulence [82,83].
It is becoming increasingly clear that regulatory changes affecting the bioactivity of phy-
tohormones also need to be analyzed. The conjugation of sugars and amino acids to
phytohormones generate a diversity of active and inactive compounds. For example,
SA is modified by SA glucosyltransferase enzymes, which are members of the uridine
diphosphate (UDP)-dependent glycosyl transferases (UGTs) family. UGT74F1 catalyzes
the addition of a glycoside to SA, to form SA-glycoside (SAG) [84]; UGT74F2 catalyzes the
formation of the SA glucose ester (SGE); [85] and UGT71C3 catalyzes the conversion of the
volatile MeSA to MeSA glycoside [86].

Defense metabolomics evolve. PIP, a lysine catabolite, was originally identified as a
key long-distance signaling molecule for SAR [87]. Subsequently, NHP was demonstrated
to be the bioactive inducer [88]. Further studies revealed that its signaling function was
more complex with glycosylated derivatives implicated [89]. Independent targeted and
untargeted metabolomics strategies [90,91] identified UGT76B1 as a glucosyltransferase
that modifies bioactive NHP, by catalyzing the formation of 1-O-glucosyl-pipecolic acid.
Interestingly, UGT76B1 can also multitask, using SA and dihydroxybenzoic acid derivatives
as substrates [92], implicating a pivotal role in modulating immune signaling.

Even with a genetically amenable pathosystem, analytical method development is
complex. This is probably best exemplified by the jasmonate family. It was only relatively
recently that (+)-7-iso-jasmonyl-L-isoleucine was recognized as the actual bioactive jas-
monate [17]. Indeed, jasmonates exist as a particularly large range of derivatives, which
contribute positively or negatively to plant defense. The recent development of a two-phase
extraction, coupled to LC-MS, for the absolute quantification of multiple jasmonates and
the semi-quantitative determination of associated jasmonate pathway metabolites [93],
represents major progress in this area.

Metabolomics approaches have also indicated a key role for nucleotides in multiple
defense responses. Plant extracellular adenosine-5-triphosphate (eATP) is a DAMP [94].
This recognition of damaged self is consistent with the reported release of cytoplasmic
ATP following wounding, leading to apoplastic eATP levels of ~40 uM [95]. Further,
eATP is recognized in Arabidopsis by the L-type lectin receptor kinases P2K1/DORN1
(does not respond to nucleotides 1) and P2K2. The activation of these PRRs induces
membrane depolarization, Ca²⁺ influx, and ROS formation, conferring enhanced immunity
against bacteria, necrotrophic fungi, and herbivores [96–99]. Paralleling a DAMP role
for ATP, the exogeneous application of NAD⁺ or NADP⁺ can also induce both local and
systemic pathogen resistance, underpinned by elevated SA and pathogenesis-related
gene (PR) expression. Similarly to eATP perception, eNAD⁺ or eNADP⁺ appear to be
recognized by lectin receptor kinases [100,101]. The establishment of pyridine nucleotides
in activated PTI responses leads to one of the most significant recent discoveries in the
field of plant immunity, which is the demonstration that the TIR (Toll-like, interleukin-1
receptor) N-terminal domains of plant resistance proteins (TNLs) undergo proximity-
induced dimerization, generating functional NADase enzymatic activity, which is essential
for the activation of ETI [102]. The plant TIRs utilize NAD/P⁺, generating a variant cyclic
ADP ribose (v-cADPR) with similar mass and mass spectra as cADPR, but a different
retention time. Notably, human TIR domains also dimerize to form functional NADases,
but their products are ADPr and cADPR. A significant immediate metabolomics challenge
is to better understand this “v-CADPR” and the near-identical derivates that are produced
by bacterial TIRs [102,103]. Thus, we are now entering a new area, where the absolute
quantitation of low-abundance (deoxy)ribonucleotides and (deoxy)ribonucleosides, and
their derivatives, such as v-cADPR, will be critical to dissecting the complex plant defense
metabolome. This will be greatly facilitated by a new method to quantify nucleotide and
nucleoside metabolism in plants [38].

It is quite remarkable that, despite our relatively comprehensive knowledge of the
molecular mechanisms underpinning PTI and ETI, and a developing understanding of how
pathogen effectors drive effector-triggered susceptibility (ETS), we still do not have a clear
understanding of the primary metabolic reprogramming that is associated with pathogen
nutrition, let alone know the key metabolites that are utilized by pathogens, as fundamental
carbon and nitrogen sources. Clearly, sugar and amino acid metabolism are modified, but
space considerations preclude detailed discussion. Instead, we provide two contrasting
examples of the effector modulation of primary metabolism for pathogen nutrition. The
activation of GABA metabolism has long been associated with plant-microbe interactions.
While it has the capacity to provide rich C and N sources, GABA appears to have contrasting
effects, depending on the pathosystem under study (reviewed: [104]). An elegant study
recently demonstrated that the *Ralstonia solanacearum* effector Rip interacts with host
glutamate decarboxylases (GADs), promoting their interaction with calmodulin, to enhance
GABA production [105]. In contrast, the indirect hijacking of host metabolism is elegantly
illustrated by the activation of host sugar transporters by pathogen effectors. For example,
*Xanthomonas oryzae* pv. oryzae, the causal agent of rice bacteria leaf streak, induces (*Oryzae
sativa*) clade III “sweet transporters” (OsSWEET11–15), to promote virulence by supporting
the translocation of sugars (OsSWEET11 transports sucrose) to the apoplast [106].

Integrating different technologies in a systems level approach will greatly improve
our understanding of plant defense responses. Untargeted metabolomics will play a critical
role in generating a large compendium of features that are significantly different in diseased
or resistant tissues, relative to their appropriate control. Integrating transcriptomics and
immunity mutants into such studies can help refine subsets of features to study further.
We have already briefly discussed the challenges of signal to noise, from the perspective
of disease and defense. While there is a potential scope to deploy cell/tissue-specific
markers, sampling remains a challenge, and is further confounded by different pathogen
virulence mechanisms and lifestyles. For example, while *Pst* forms discrete apoplastic colonies associated with mesophyll cells, many *Xanthomonas* spp. are vascular pathogens, reflecting the importance of developing methods for spatial sampling. The application of untargeted metabolomic approaches, to investigate stomatal immunity [107], is a recent example of such evolving spatial sampling methods.

4. Metabolomics Research in Fungal and Oomycete Plant Pathology

Metabolomic approaches are now well established in the assessment of plant pathogenic
fungal and oomycete interactions, and have been extensively considered elsewhere [16,108].
These studies have revealed similar changes in metabolism as previously described in this
review, although, unsurprisingly, each discrete plant–fungal/oomycete interaction shows
distinctive features. These features have been classified into changes in the following three
broad major biochemical groups: alkaloids, isoprenoids, and shikimates [16,109]. Exemplar
interactions include alkaloid changes in responses to *Fusarium gramininum*, isoprenoids
to *Magnaporthe oryzae*, and shikimates to *Puccinia coronta* [109]. However, metabolomic
analyses of fungal- and oomycete-plant interactions represent a considerable challenge,
due to the complex spatio-temporal changes that we have already referred to. For example,
fungal and oomycete pathogens will often be dispersed via asexually or sexually produced
spores, to germinate on the plant surface, and may germinate to form discrete infection
structures, which are initially deployed against only a few cells. The key examples are
outlined below.

*M. oryzae*—the causal agent of rice blast disease—can infect a range of cereal and grass hosts [110]. *M. oryzae* spore contact on a hydrophobic plant surface leads to germ-tube development and the formation of a melanized appressorium. Metabolite processing, of glycerol for example, with the appressorium, builds up enormous turgor pressures that contribute to the forcing of a penetration peg through the plant cuticle. After penetration, the primary invasive hyphae fill the cell, but still form an intact biotrophic interfacial complex (BIC) with the host, which is vital to the delivery of effectors [111]. Infection hyphae will move to neighboring host cells via plasmodesmata, and ultimately trigger a switch to a necrotrophic mode of pathogenesis [112]. Although phylogenetically distinct, oomycetes exhibit similar infection strategies. With the oomycete *Phytophthora infestans*, the dispersal of asexual zoospores leads to encystment on the plant surface. Cyst germination leads to the formation of an appressorium, penetration through the epidermis (occasionally via the stomata), and the formation of a haustorium, when effectors begin to be delivered. Further hyphal development leads to the formation of additional haustoria in cells in close proximity with the leaf [113,114]. Necrotrophic fungi tend to have less-sophisticated infection strategies, often focusing on the production of toxins and cell wall-degrading enzymes [115]. An important example here is *Fusarium graminearum*, which causes Fusarium head blight (FHB) on a range of cereal crops [116]. This species can secrete a series of mycotoxins, with the most well characterized being deoxynivalenol (DON). It should also be recognized that defense reactions exhibit spatial differences with the formations of cell wall appositions, as well as single or multicellular cell death being displayed at the infection site [117].

*M. oryzae* interactions have been well characterized using metabolomics approaches. Studying *M. oryzae* appressoria formation on a synthetic hydrophobic surface revealed ceramide processing as an important feature of appressorial maturation, linked to the protein kinase C-mediated cell wall integrity pathway [118]. Focusing on plant responses, three hosts of *M. oryzae* showed common changes in amino acids and sugars [110], which likely reflect nutritional reprogramming towards the host [109]. A more comprehensive approach, integrating GC-MS/MS, LC-MS/MS, and ¹H-NMR assessment of the different responses of rice to virulent and avirulent *M. oryzae* strains, revealed a significant correlation between alanine and disease development, potentially aiding infection by triggering host cell death [119]. By considering temporal changes, the metabolic "conflict" between pathogen nutrient mobilization and host defenses was revealed. During the early stages
of infection, there were changes in bioenergetic, polyamine and oxidative metabolism, which is indicative of a defense response, but by the stage of lesion formation, changes in mannitol and glycerol were evident, possibly supporting fungal growth [110]. Resistance to *M. oryzae* may be associated with wide-ranging phospholipid processing, leading to the production of JA, as occurs in the model grass *Brachypodium distachyon* [30]. Genetic evidence in rice has also suggested the importance of JAs, as well as SA and ethylene, in resistance [120].

The *Phytophthora infestans*–tomato interaction is a good example where the transcrip-
tome and metabolomes have been integrated, to derive a robust infection model [121].
Predictive genome-scale metabolic models (GEMs) were derived from *P. infestans* enzyme
sequences, which were related to KEGG metabolic networks [122]. This GEM was then
improved by incorporating further information, such as organellar location and, crucially,
transcriptomic and metabolomic data [121]. These models showed how *P. infestans* growth
was dependent on host thiamine biosynthesis. Moreover, separate consideration of the
host and oomycete transcriptomic data showed how the pathogen increasingly shuts down
some of its own biosynthetic pathways, becoming more dependent on host-sourced nutri-
ents. Although undoubtedly challenging, such multi-omic approaches offer the possibility
of major new insights.

Due to the multigenic nature of defense against necrotrophs, metabolomics is a power-
ful tool that can generate infection biomarkers or, through the use of metabolite (m)GWAS
approaches [123], accelerate the breeding of resistant varieties [124,125]. LC-MS has been
used to define metabolites associated with a QTL that is linked to resistance against *F.*
*graminearum* in barley [125]. The FHB1 QTL was linked to DON detoxification and, through
metabolomics, also to changes in phenylpropanoid and JA metabolism [126]. In wider
studies, metabolites linked to constitutive and induced resistance against head blight in-
fections, have been characterized. With DON, the possible biomarkers for this mycotoxin
included changes in JAs and the fatty alcohol dihydro-7-hydroxymyoporone [127]. Impor-
tantly, responses to the mycotoxin trichothecene were linked to indole acetic acid, picolinic
acid, and feruloyl glucosides [128]. *Botrytis cinerea* is an important necrotrophic disease
of, particularly, soft fruits, where the identification of biomarkers in presymptomatic ("la-
tent") infected fruit could help mitigate product spoilage. Recently, untargeted GC-MS
approaches identified hexadecanoic acid, octadecanoic acid, sucrose, β-lyxopyranose,
melibiose, and 1,1,4a-Trimethyl-5,6-dimethylenedecahydronaphthalene changes in the
strawberry latent phase [129].

To address the issue of spatial differences, our early research demonstrated the use of
FT-IR microscopy in Arabidopsis, following challenge with *Botrytis cinerea*. Using a series
of ethylene defense mutants, imaging at the infection site defined cell wall changes linked
to ferulate [27]. More powerfully, MSi approaches can be used where samples are scanned
across a predefined field and the intensity of analyte ions is generated at a single point. This
generates a pixel reflecting analyte intensity, which can be assembled to form a heat-map
image [130,131]. An elegant application of MSi technologies included the use of MALDI-
MS [132], to reveal changes in the distribution of glucosinolates in Arabidopsis, which
influenced the feeding patterns of cotton bollworm (*Helicoverpa armigera*) larvae [133]. Both
FT-IR and MALDI-MSi cannot be easily applied to living tissue, but this is not the case for
LAESI-MS. Studying the interaction of *Cladosporium fulvum* on tomato, LAESI-MS provided
spatial information on the fungal toxic glycoalkaloid α-tomatine, which suggested its
degradation at the site of the infection [134]. Metabolite imaging approaches have yet to be
widely applied to plant pathogen interactions, but their potential impact in understanding
metabolomic changes could be enormous.

**5. Metabolomics Research in Plant—Pest Resistance Breeding**

Despite all the phytosanitary measures taken by farmers, crops commonly will also
be confronted with insect herbivores. It is estimated that ~15% of all crop production is
lost to insect herbivore damage annually [135]. Synthetic pesticides have been used to
effectively reduce pest pressure and increase crop productivity since the 1950s. The current concerns on the non-target effects of synthetic pesticides, amongst others on pollinating insects, have led to the ban of several effective pesticides, such as neonicotinoids [136]. In addition, several herbivore species have developed pesticide tolerance. This has created the urge for more sustainable herbivore control measures. These strategies go hand in hand with a stronger ability of the plants to defend themselves. Plant breeders are thus left in the unenviable position of balancing the production of novel, pathogen, and pest-resistant varieties that are both tolerant to changing climate conditions and appropriate for the food market.

Metabolomics, combined with bioassays and genetic analyses, is a rapid method to identify specific leads for breeding efforts where natural resistance to pests is paramount. To identify (new) sources of chemical resistance, the large-scale panel screening of existing accessions and wild relatives, using bioassays, is a powerful tool [137]. With a focus on subsets with contrasting levels of resistance, metabolomic analyses can identify specific compounds whose concentrations differ between the susceptible and the resistant accessions. For example, an LC-MS-based metabolomics approach was used to pinpoint dimer acyclic diterpene glycosides (capsianosides) as resistance factors against thrips in Capsicum [53]. In addition, a GC-MS-based metabolomics approach was used to analyze which wax components were related to thrips resistance [53]. For plant breeders, linking these resistance traits to genetic markers is essential. The common approaches are QTL analysis or GWAS mapping [138]. With GWAS, natural accessions are genotyped by sequencing. These same accessions are subjected to metabolomics, which can be combined with QTL maps, yielding information on the genomic regions that certain metabolites, or mQTL, are associated with. Resistance factors have been identified by relating the gene sequence with the metabolome of the individuals [138]. Both mapping methods combined with metabolomics will identify accessions that are associated with resistance and particular metabolic traits. Only by further genetic fine mapping, the responsible regulatory or biosynthetic genes can be identified [138,139]. In addition, these approaches could be combined with targeted genome editing tools, such as CRISPR-Cas9, which precisely target specific DNA sequences [140]. Combined with bioassays, targeted mutations can provide proof-of-principle on the function of specific genes and metabolites. High-throughput metabolomic analyses can identify such desired metabolites in the breeding stock, and thus can be part of other speed-breeding approaches [141], to quickly yield natural insect-resistant accessions. One issue, however, is that plant defenses are often important determinants of the palatability and nutritional value of produce and fodder. Depending on the compound chemistry, plant defense may either be health-beneficial or toxic. Therefore, most crops have been selected to reduce toxicity and bitterness, to increase palatability. A famous example is the reduced concentrations of the glucosinolate progoitrin in oil seed rape. In mammals, progoitrin negatively effects thyroid metabolism, thus causing health issues [142]. The trade-off, however, is that low progoitrin accessions become more susceptible to herbivores [143]. Metabolomic analyses, combined with bioassays, may help to optimize the levels of chemical defenses in crops, without decreasing the nutritional value and marketability of the produce.

## 6. Unravelling the Response of Plants to Stress in the Natural and Agricultural Environment

As plants are adapted to surviving a multitude of climatic extremes over their growth season, which often occur simultaneously, within close proximity or at multiple points during development, their metabolic response needs to be plastic [144]. With increasingly frequent weather extremes in agriculturally valuable areas, understanding the metabolic underpinnings of the plant response to environmental stress is an important avenue for safeguarding natural and managed ecosystems [145,146]. Due to the vast suite of primary and secondary metabolites that plants employ to aide in combined (a)biotic stress responses, developing and employing metabolomics technologies is imperative to discern
the interactive impacts on system level metabolomics, and highlight potential agricultural exploitation for stress mitigation strategies, as discussed in the previous section.

Targeted metabolomics allows relative or absolute quantification of individual, or a small number, of metabolites. For example, phytohormones are often measured using targeted approaches [24,25]. Considering future CO₂ levels, GC-chemical ionization (CI)-MS revealed that reduced JA production in maize, at elevated CO₂ (eCO₂), coincided with reduced phytoalexin production and increased susceptibility to Fusarium wilt [147]. Similarly, HPLC-fluorescence was used to show that increases in SA, at eCO₂, coincided with increased scopoletin levels and resistance to *Pst* in Arabidopsis [148]. Interestingly, phytohormonal changes under eCO₂ link to changes in the structure of the secondary metabolome, which impacts the immune response. In eCO₂, untargeted GC-MS revealed that changes in the primary metabolite profile overlapped with increases in an oxidative stress background, and thus highlighted a role for redox metabolism of the plant during CO₂ stress [148]. Similarly, untargeted UHPLC-ESI-MS was used to show that plants that were grown in low CO₂ had a secondary metabolite profile that was consistent with greater ROS production, oxidative stress, and resistance to biotrophic pathogens [149]. Low CO₂ altered the photorespiration machinery, with knock-on impacts on the sugar metabolism and antioxidant status of the plant. This also suggests that photorespiration played a greater role in immunity in past CO₂ climates.

Metabolomics has found substantial application in drought research too. From directed studies, such as GC demonstrating that mannitol accumulates in the roots of droughted olives [78], to targeted LC-MS methods, illustrating major changes in many primary and secondary pathways, including amino acid metabolism and flavonoid biosynthesis in JA mutant backgrounds [150]. In drought-tolerant sesame, untargeted GC-MS identified accumulation of ABA, proline, arginine, lysine, GABA, saccharopine, 2-amino adipate, and allantoin [151], many of which are implicated and are of importance in plant-microbe responses. In this case, by coupling metabolome and transcriptome data, a strong argument was demonstrated that ABA itself is likely altering metabolite profiles, to generate a more osmotically stable environment. This is complementary with ABAs role in drought-induced stomatal closure, increased stomatal defense [152], and it serving as a direct signal in the rhizosphere, to recruit drought-protective microbes [153]. Combining different omics approaches is a powerful tool for better understanding and unravelling complex tripartite interactions. For instance, combined transcriptome and metabolome data, measured via LC-MS, in Arabidopsis experiencing combined drought, warming and viral infection, indicated that signaling networks shifted in a stress-specific manner, and resulted in a breakdown in plant defenses under multi-stress [154]. Further, the impact of heat stress can impact the pathogen response, by requiring a contrasting metabolomic strategy, rendering the plant more susceptible to disease [155]. In maize, the metabolomes under heat stress and disease were shown to be divergent, using both untargeted and targeted UHPLC-MS. In a combined stress scenario, hydroxycinnamic acid and p-coumaric acid accumulation were specifically identified to contribute to the heat-induced susceptibility of maize to the necrotrophic leaf pathogen *Cochliobolus heterostrophus* [155]. Previous examples illustrate that individual stresses are inadequate for determining the mechanisms of high-temperature-induced disease susceptibility [156], which is likely true across all (a)biotic stresses that require contrasting plant responses.

The benefits of the application of metabolomics technologies, to investigate complex tripartite stress interactions, are clear, although challenges remain. Future innovation will permit an integrated view of changes in dynamic plant metabolomes, across the tissue spectrum of the plant, in natural field systems with specific dissection of host and microbial metabolomic responses to a suite of environmental stresses.

## 7. Future Prospects: Which New Tools and Approaches Will Make an Impact?

Metabolomics is a valuable tool for discovering metabolic signals for which the compound identification procedure is underprovided and tedious [157]. Advances in
bioinformatics have dramatically improved metabolomic data processing [158,159], but, despite this, thousands of signals remain unidentified [160]. Thus, metabolomics could have a much greater impact, with the development of improved annotation workflows and more comprehensive plant metabolite libraries. Machine learning now plays a significant role in metabolomics (LC-MS being routinely utilized [161]). It predicts a vast array of molecular properties and processes, such as adduct species, MS fragmentation, MS/MS deconvolution, and chromatographic retention time [158,159,162]. However, it must still be noted that, due to the vast numbers of closely related isomeric compounds in plants, as well as their derivatives decorated with varying isomeric sugar moieties or fused with amino acids, which produce near-identical MS/MS and MS<sup>n</sup> fragmentation patterns, and which co-elute even with long (60 > min) UHPLC separations, there needs to be a focus upon developing higher resolution separation methodologies and technologies that can be combined with both MS and 2D NMR. The development of 2D LC separations, as well as time-series retention index standards that are applicable to LC, as n-alkanes and fatty acid methyl esters (FAMES) are applied in GC, would provide a massive step change towards developing retention time (RT)-MS/MS and MS<sup>n</sup> libraries that are transferable between different LC-MS systems and laboratories. Overcoming such hurdles and enhancing our metabolite annotation capabilities will be particularly valuable to the discovery of antimicrobial and viral molecules in plant-pathogen research. A recent study reports the metabolic setting of plant-virus interactions in a natural habitat and its predictive link to susceptibility, unravelling new insight into plant-virus interactions and providing attractive antimicrobial targets [163]. Likewise, machine learning-based predictive metabolomics is particularly promising, to reveal and validate metabolic predictors under biotic stress [164]. Future bioinformatics developments for metabolomics will reduce human intervention and increase data analysis throughput, thanks to robotized assays, for instance [157,164].

As previously mentioned, MSi is prominent among the new tools that have been developed for metabolome analysis [165], and offers promising perspectives for the spatial dissection of metabolic responses to infection, although without chromatographic separation it does lack the ability to differentiate structural isomers that display identical MS/MS and MS<sup>n</sup> fragmentation patterns. However, by applying ion mobility separations prior to HRMS detection, there is the potential to resolve isomeric compounds, without prior chromatographic separation. Recent salient MSi examples include the barley defense metabolite serotonin and other metabolites of the melanin pathway against *M. oryzae* [166], and the sulfur-containing glucosinolates in Arabidopsis, which intervene in fungal and pest responses [133]. Besides MSi, other approaches allowing the isolation of cells and organelles (e.g., flow cytometry, cell and organelle sorting, microfluidics) will make an outstanding contribution to the future impacts of metabolomics on plant stress studies.

Genome-scale metabolic modelling also holds great promise to study plant-pathogen interactions. This computational systems biology perspective implies structural and comparative genomics, transcriptomics, and protein-protein interactions [167]. The calculation of metabolite coefficients for genome-scale metabolic flux predictions requires high-resolution metabolome data (easily provided by today's Orbitrap and Q-TOF technologies) and high-resolution temporal sampling, as emphasized by a typical experimental design of phytopathological metabolomics (Figure 1). In this context, genome-scale metabolic models of potato and tomato leaves and *Phytophthora* have simulated the metabolic fluxes during infection [121,168]. Interesting molecular aspects, involving photosynthesis suppression by the pathogen, were demonstrated, as well as the uptake of nutrients by *Phytophthora* following the infection cycle phases. It should be noted that stage-specific profiles of the joint metabolism of the host and pathogen could be refined by integrating high-resolution data from the tomato infection metabolome [121].

To echo this last point, one of the main challenges of metabolomics in plant pathology is distinguishing plant and microbial metabolomes [80]. Metabolites, most particularly primary compounds, are often ubiquitous in the diverse plant, microbial and animal organisms, or exchanged between the host and pathogens, as signals or metabolic responses. In
pathological studies (Figure 1), plants and microbes are jointly incubated, but the pathogens
are rarely removed before metabolite extraction, thereby mixing metabolomes of different
origins. This probably seems to be of minor importance for other omics, such as transcrip-
tomics and proteomics, which exploit species-specific databases. However, metabolomics
suffers from cross-contamination of the host and non-host metabolites; therefore, it is often
difficult to draw precise conclusions as to the origin of the metabolic variations found.
Whilst one can search data for metabolites of non-plant origin, for example, arachidonic
acid with respect to fungal pathogens, one cannot assume that their absence suggests
that metabolites of fungal origin are not being detected. Implementing stable isotope
labelling of microbes could address this pressing concern. A recent metabolomics study
of the early stage of plant cell–microbe interactions employed an incubation method of
¹³C- and ¹⁵N-labelled (98%) *Pst* DC3000 with Arabidopsis cells, followed by a washing
step to remove the excess bacterial cells [169]. Subsequently, “light” plant metabolites were
analyzed by targeted LC-MS, highlighting an infection reprogramming of plant signaling
and primary metabolic pathways, including carbohydrates, nucleotides, and amino and
organic acids. Alternatively, uniformly ¹³C- and ¹⁵N-labelled plants can be prepared, which
will offer an accurate tool for metabolic analysis under biotic stress conditions [170]. The
growing interest in isotopically labelled biological material coincides with the development
of additional analytical tools, such as NMR and high-throughput GC- and LC-HRMS,
which allow the discrimination between native and labelled metabolites (Figure 3). It is
conceivable that elegant studies of this type will be carried out less sporadically, to help
clarify the origin of changes in the host and pathogen metabolome. In summary, future
prospects that will greatly contribute to the success of metabolomics in the study of plant
immunity will include (i) increasing the capacity to chromatopgraphically resolve and
 annotate plant and pathogen metabolomes, (ii) new tools and approaches to analyze the
metabolome at the subcellular level, and (iii) the development and implementation of
modelling, to predict and/or calculate the dynamics of metabolic variations.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** LC-HRMS discrimination between native and labelled metabolites. Unlabelled (light) plant tissues (e.g., leaves) are inoculated with fully labelled (heavy) pathogenic microbe (e.g., ¹³C isotopically labelled bacteria) then sampled at the appropriate time point for metabolome analysis (left). The resulting metabolite extracts thus contains a mixture of light and heavy biochemical (middle) that could be discriminated by LC-HRMS technique (right). Light and heavy metabolites display similar retention time (RT) on the LC-HRMS chromatogram, but natural isotope clusters of light and heavy forms are well separated on the mass spectra.</figcaption></figure>
## 8. Concluding Remarks

Metabolomics has revolutionized the understanding of plant-microbe interactions, by providing a valuable biochemical phenotypic screen of the chemical diversity of plants and their metabolic responses to biotic challenges. Targeted approaches have permitted a better insight into the role of key molecules that are involved in these phytopathological mechanisms, such as phytohormones, and other defense and signaling biochemicals. Untargeted metabolomics has also pivotally unveiled new and unsuspected metabolic players, such as pyridine cofactors and ATP, for instance. However, the use of metabolomics in the study of plant-microbe interactions still lags other omics approaches, notably transcriptomics. This problem is even more imperative for plant-virus interactions, for which there has been very little research. At a time when humanity is suffering from an unprecedented pandemic of COVID-19, the study of the metabolic pathways that are involved in viral resistance, as well as the discovery of plant antiviral molecules, has a bright future [162]. Beyond metabolite profiling, the forthcoming trend is likely towards more holistic and comprehensive large-scale multi-omics strategies, critically providing necessary insight at the systems level [171]. A closer association between metabolomics and chemical ecology will also facilitate discovering a plethora of metabolites with ecological and defensive roles [172].

**Author Contributions:** J.W.A. lead on conceptualisation with contributions from all authors, A.W., H.U., N.M.v.D., L.A.J.M., M.R.G. and P.P. contributed equally to writing—original draft preparation, writing—review and editing. P.P. and A.W. contributed the visualization. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding. J.W.A. acknowledges support from the Scottish Government Rural and Environmental Science and Analytical Services. P.P. is grateful for support from MetaboHUB (ANR-11-INBS-0010) and PHENOME (ANR-11-INBS-0012) projects. H.U. and N.M.v.D. gratefully acknowledge support of iDiv funded by the German Research Foundation (DFG-FZT 118, 202548816).

**Conflicts of Interest:** The authors declare no conflict of interest.

## Abbreviations

<table><tr><td><sup>1</sup>H-NMR</td><td>proton nuclear magnetic resonance spectroscopy</td></tr><tr><td>2D NMR</td><td>two-dimensional nuclear magnetic resonance spectroscopy</td></tr><tr><td>ABA</td><td>abscisic acid</td></tr><tr><td>ADP</td><td>adenosine diphosphate</td></tr><tr><td>APCI</td><td>atmospheric pressure chemical ionization</td></tr><tr><td>ATP</td><td>adenosine triphosphate</td></tr><tr><td>BABA</td><td>β-amino butyric acid</td></tr><tr><td>cADPR</td><td>variant cyclic ADP ribose</td></tr><tr><td>CAPK</td><td>calcium-activated protein kinase</td></tr><tr><td>CI</td><td>chemical ionization</td></tr><tr><td>CO<sub>2</sub></td><td>carbon dioxide</td></tr><tr><td>DAMP</td><td>damage-associated molecular patterns</td></tr><tr><td>DESI</td><td>desorption electrospray ionization</td></tr><tr><td>DI</td><td>direct infusion</td></tr><tr><td>DON</td><td>deoxynivalenol</td></tr><tr><td>eATP</td><td>extracellular adenosine-5-triphosphate</td></tr><tr><td>εCO<sub>2</sub></td><td>elevated carbon dioxide</td></tr><tr><td>EI</td><td>electron ionization</td></tr><tr><td>EPR</td><td>electron paramagnetic resonance</td></tr><tr><td>ESI</td><td>electrospray ionization</td></tr><tr><td>ETI</td><td>effector-triggered immunity</td></tr><tr><td>ETS</td><td>effector-triggered susceptibility</td></tr><tr><td>FAMES</td><td>fatty acid methyl esters</td></tr><tr><td>FHB</td><td>Fusarium head blight</td></tr></table>
<table>
  <tr>
    <td>FI</td>
    <td>flow infusion</td>
  </tr>
  <tr>
    <td>FID</td>
    <td>flame ionization detection</td>
  </tr>
  <tr>
    <td>FT-IR</td>
    <td>Fourier transform infrared</td>
  </tr>
  <tr>
    <td>GABA</td>
    <td>γ-amino butyric acid</td>
  </tr>
  <tr>
    <td>GAD</td>
    <td>glutamate decarboxylase</td>
  </tr>
  <tr>
    <td>GC</td>
    <td>gas chromatography</td>
  </tr>
  <tr>
    <td>GEM</td>
    <td>genome-scale metabolic models</td>
  </tr>
  <tr>
    <td>GWAS</td>
    <td>genome-wide association studies</td>
  </tr>
  <tr>
    <td>HILIC</td>
    <td>hydrophobic interaction liquid chromatography</td>
  </tr>
  <tr>
    <td>HPLC</td>
    <td>high-performance liquid chromatography</td>
  </tr>
  <tr>
    <td>HR</td>
    <td>hypersensitive response</td>
  </tr>
  <tr>
    <td>HRMS</td>
    <td>high-resolution mass spectrometry</td>
  </tr>
  <tr>
    <td>JA</td>
    <td>jasmonic acid</td>
  </tr>
  <tr>
    <td>LAESI</td>
    <td>laser ablation electrospray ionization</td>
  </tr>
  <tr>
    <td>LC</td>
    <td>liquid chromatography</td>
  </tr>
  <tr>
    <td>LPS</td>
    <td>lipopolysaccharide</td>
  </tr>
  <tr>
    <td>MALDI</td>
    <td>matrix-assisted laser desorption ionization</td>
  </tr>
  <tr>
    <td>MAMP</td>
    <td>microbe-associated molecular patterns</td>
  </tr>
  <tr>
    <td>MAPK</td>
    <td>mitogen-activated protein kinase</td>
  </tr>
  <tr>
    <td>MeJA</td>
    <td>methyl jasmonic acid</td>
  </tr>
  <tr>
    <td>MeSA</td>
    <td>methyl salicylic acid</td>
  </tr>
  <tr>
    <td>MOX-TMS</td>
    <td>methoximated trimethylsilyl</td>
  </tr>
  <tr>
    <td>mQTL</td>
    <td>metabolite quantitative trait loci</td>
  </tr>
  <tr>
    <td>MS</td>
    <td>mass spectrometry</td>
  </tr>
  <tr>
    <td>MS/MS</td>
    <td>tandem MS</td>
  </tr>
  <tr>
    <td>MSi</td>
    <td>mass spectrometry imaging</td>
  </tr>
  <tr>
    <td>MS<sup>n</sup></td>
    <td>multi-stage mass spectrometry</td>
  </tr>
  <tr>
    <td>MTI</td>
    <td>microbe-triggered immunity</td>
  </tr>
  <tr>
    <td>MW</td>
    <td>molecular weight</td>
  </tr>
  <tr>
    <td>NAD</td>
    <td>nicotinamide adenine dinucleotide</td>
  </tr>
  <tr>
    <td>NADH</td>
    <td>nicotinamide adenine dinucleotide + hydrogen</td>
  </tr>
  <tr>
    <td>NADPH</td>
    <td>nicotinamide adenine dinucleotide phosphate</td>
  </tr>
  <tr>
    <td>NAMP</td>
    <td>nematode-associated molecular patterns</td>
  </tr>
  <tr>
    <td>NHP</td>
    <td>n-hydroxy pipecolic acid</td>
  </tr>
  <tr>
    <td>NMR</td>
    <td>nuclear magnetic resonance spectroscopy</td>
  </tr>
  <tr>
    <td>NO</td>
    <td>nitric oxide</td>
  </tr>
  <tr>
    <td>OPDA</td>
    <td>12-oxo-phytodienoic acid</td>
  </tr>
  <tr>
    <td>PA</td>
    <td>phosphatidic acid</td>
  </tr>
  <tr>
    <td>PAMP</td>
    <td>pathogen-associated molecular patterns</td>
  </tr>
  <tr>
    <td>PAP</td>
    <td>3'-phosphoadenosine 5'-phosphate</td>
  </tr>
  <tr>
    <td>PIP</td>
    <td>pipecolic acid</td>
  </tr>
  <tr>
    <td>PRR</td>
    <td>pattern recognition receptors</td>
  </tr>
  <tr>
    <td>Pst</td>
    <td>Pseudomonas syringae pv. Tomato</td>
  </tr>
  <tr>
    <td>PTD</td>
    <td>pathogen-triggered defense</td>
  </tr>
  <tr>
    <td>PTI</td>
    <td>pathogen-triggered immunity</td>
  </tr>
  <tr>
    <td>PTR</td>
    <td>proton transfer reaction</td>
  </tr>
  <tr>
    <td>QqQ</td>
    <td>triple-quadrupole mass spectrometry</td>
  </tr>
  <tr>
    <td>QTL</td>
    <td>quantitative trait loci</td>
  </tr>
  <tr>
    <td>Q-ToF</td>
    <td>quadrupole time-of-flight</td>
  </tr>
  <tr>
    <td>R protein</td>
    <td>resistance protein</td>
  </tr>
  <tr>
    <td>ROS</td>
    <td>reactive oxygen species</td>
  </tr>
  <tr>
    <td>RT</td>
    <td>retention time</td>
  </tr>
  <tr>
    <td>SA</td>
    <td>salicylic acid</td>
  </tr>
  <tr>
    <td>SAG</td>
    <td>salicylic acid glycoside</td>
  </tr>
  <tr>
    <td>SAR</td>
    <td>systemic acquired resistance</td>
  </tr>
  <tr>
    <td>SIMS</td>
    <td>secondary ionization mass spectrometry</td>
  </tr>
</table>
<table><tr><td>SPME</td><td>solid-phase micro extraction</td></tr><tr><td>TD</td><td>thermal desorption</td></tr><tr><td>TIR</td><td>Toll-like, interleukin-1 receptor</td></tr><tr><td>TNL</td><td>N-terminal domains</td></tr><tr><td>ToF</td><td>time-of-flight</td></tr><tr><td>UHPLC</td><td>ultra-high-performance liquid chromatography</td></tr><tr><td>v-cADPR</td><td>variant cyclic ADP ribose</td></tr></table>

References

1. Piasecka, A.; Kachlicki, P.; Stobiecki, M. Analytical Methods for Detection of Plant Metabolomes Changes in Response to Biotic and Abiotic Stresses. *Int. J. Mol. Sci.* **2019**, *20*, 379. [CrossRef]

2. Choi, H.W.; Klessig, D.F. DAMPs, MAMPs, and NAMPs in plant innate immunity. *BMC Plant Biol.* **2016**, *16*, 232. [CrossRef] [PubMed]

3. Peng, Y.; van Wersch, R.; Zhang, Y. Convergent and Divergent Signaling in PAMP-Triggered Immunity and Effector-Triggered Immunity. *Mol. Plant-Microbe Interact.* **2018**, *31*, 403–409. [CrossRef]

4. Jones, J.D.G.; Dangl, J.L. The plant immune system. *Nature* **2006**, *444*, 323–329. [CrossRef]

5. Van der Burgh, A.M.; Joosten, M.H.A.J. Plant Immunity: Thinking Outside and Inside the Box. *Trends Plant Sci.* **2019**, *24*, 587–601. [CrossRef] [PubMed]

6. Vlot, A.C.; Sales, J.H.; Lenk, M.; Bauer, K.; Brambilla, A.; Sommer, A.; Chen, Y.; Wenig, M.; Nayem, S. Systemic propagation of immunity in plants. *New Phytol.* **2021**, *229*, 1234–1250. [CrossRef] [PubMed]

7. Kachroo, P.; Burch-Smith, T.; Grant, M. An Emerging Role for Chloroplasts in Disease and Defense. *Annu. Rev. Phytopathol.* **2021**, *59*, 432–445.

8. Mur, L.; Prats, E.; Pierre, S.; Hall, M.; Hebelstrup, K.H. Integrating nitric oxide into salicylic acid and jasmonic acid/ethylene plant defense pathways. *Front. Plant Sci.* **2013**, *4*, 215. [CrossRef] [PubMed]

9. Xing, J.; Zhang, L.; Duan, Z.; Lin, J. Coordination of Phospholipid-Based Signaling and Membrane Trafficking in Plant Immunity. *Trends Plant Sci.* **2021**, *26*, 407–420. [CrossRef]

10. Vlot, A.C.; Pabst, E.; Riedlmeier, M. Systemic Signalling in Plant Defence. In *eLS*; John Wiley & Sons: Hoboken, NJ, USA, 2017; pp. 1–9. [CrossRef]

11. Takahashi, F.; Shinozaki, K. Long-distance signaling in plant stress response. *Curr. Opin. Plant Biol.* **2019**, *47*, 106–111. [CrossRef]

12. Lin, L.; Wu, J.; Jiang, M.; Wang, Y. Plant Mitogen-Activated Protein Kinase Cascades in Environmental Stresses. *Int. J. Mol. Sci.* **2021**, *22*, 1543. [CrossRef] [PubMed]

13. Roychoudhury, A.; Banerjee, A. Abscisic Acid Signaling and Involvement of Mitogen Activated Protein Kinases and Calcium-Dependent Protein Kinases during Plant Abiotic Stress. In *Mechanism of Plant Hormone Signaling under Stress*; Pandey, G.K., Ed.; John and Wiley and Sons: Hoboken, NJ, USA, 2017; pp. 197–241.

14. Su, J.; Yang, L.; Zhu, Q.; Wu, H.; He, Y.; Liu, Y.; Xu, J.; Jiang, D.; Zhang, S. Active photosynthetic inhibition mediated by MPK3/MPK6 is critical to effector-triggered immunity. *PLoS Biol.* **2018**, *16*, e2004122. [CrossRef]

15. Littlejohn, G.R.; Breen, S.; Smirnoff, N.; Grant, M. Chloroplast immunity illuminated. *New Phytol.* **2021**, *229*, 3088–3107. [CrossRef] [PubMed]

16. Castro-Moretti, F.R.; Gentzel, I.N.; Mackey, D.; Alonso, A.P. Metabolomics as an Emerging Tool for the Study of Plant-Pathogen Interactions. *Metabolites* **2020**, *10*, 52. [CrossRef] [PubMed]

17. Fonseca, S.; Chini, A.; Hamberg, M.; Adie, B.; Porzel, A.; Kramell, R.; Miersch, O.; Wasternack, C.; Solano, R. (+)-7-iso-Jasmonoyl-L-isoleucine is the endogenous bioactive jasmonate. *Nat. Chem. Biol.* **2009**, *5*, 344–350. [CrossRef]

18. Sakata, N.; Ishiga, T.; Masuo, S.; Hashimoto, Y.; Ishiga, Y. Coronatine contributes to *Pseudomonas cannabina* pv. *alisalensis* virulence by overcoming both stomatal and apoplastic defenses in dicot and monocot plants. *Mol. Plant-Microbe Interact.* **2021**. [CrossRef]

19. Schulz, P.; Herde, M.; Romeis, T. Calcium-Dependent Protein Kinases: Hubs in Plant Stress Signaling and Development. *Plant Physiol.* **2013**, *163*, 523–530. [CrossRef]

20. Jagodzik, P.; Tajdel-Zielinska, M.; Ciesla, A.; Marczak, M.; Ludwikow, A. Mitogen-Activated Protein Kinase Cascades in Plant Hormone Signaling. *Front. Plant Sci.* **2018**, *9*, 1387. [CrossRef]

21. Wang, P.; Hsu, C.-C.; Du, Y.; Zhu, P.; Zhao, C.; Fu, X.; Zhang, C.; Paez, J.S.; Macho, A.P.; Tao, W.A.; et al. Mapping proteome-wide targets of protein kinases in plant stress responses. *Proc. Natl. Acad. Sci. USA* **2020**, *117*, 3270. [CrossRef]

22. Eto, M.; Katsuki, S.; Tanaka, Y.; Takeya, K. Kinase activity-tagged western blotting assay. *BioTechniques* **2020**, *68*, 211–213. [CrossRef]

23. Bürger, M.; Chory, J. Stressed Out About Hormones: How Plants Orchestrate Immunity. *Cell Host Microbe* **2019**, *26*, 163–172. [CrossRef]

24. Šimura, J.; Antoniadi, I.; Široká, J.; Tarkowská, D.E.; Strnad, M.; Ljung, K.; Novák, O. Plant Hormonomics: Multiple Phytohormone Profiling by Targeted Metabolomics. *Plant Physiol.* **2018**, *177*, 476–489. [CrossRef] [PubMed]

25. Rawlinson, C.; Kamphuis, L.G.; Gummer, J.P.A.; Singh, K.B.; Trengove, R.D. A rapid method for profiling of volatile and semi-volatile phytohormones using methyl chloroformate derivatisation and GC-MS. *Metabolomics* **2015**, *11*, 1922–1933. [CrossRef] [PubMed]
26. BASSI, P.K.; SPENCER, M.S. Comparative evaluation of photoionization and flame ionization detectors for ethylene analysis. *Plant Cell Environ.* **1985**, *8*, 161–165. [CrossRef]

27. Lloyd, A.J.; William Allwood, J.; Winder, C.L.; Dunn, W.B.; Heald, J.K.; Cristescu, S.M.; Sivakumaran, A.; Harren, F.J.M.; Mulema, J.; Denby, K.; et al. Metabolomic approaches reveal that cell wall modifications play a major role in ethylene-mediated resistance against *Botrytis cinerea*. *Plant J.* **2011**, *67*, 852–868. [CrossRef] [PubMed]

28. Hartmann, M.; Kim, D.; Bernsdorff, F.; Ajami-Rashidi, Z.; Scholten, N.; Schreiber, S.; Zeier, T.; Schuck, S.; Reichel-Deland, V.; Zeier, J. Biochemical Principles and Functional Aspects of Pipecolic Acid Biosynthesis in Plant Immunity. *Plant Physiol.* **2017**, *174*, 124–153. [CrossRef]

29. Ding, P.; Rekhter, D.; Ding, Y.; Feussner, K.; Busta, L.; Haroth, S.; Xu, S.; Li, X.; Jetter, R.; Feussner, I.; et al. Characterization of a Pipecolic Acid Biosynthesis Pathway Required for Systemic Acquired Resistance. *Plant Cell* **2016**, *28*, 2603–2615. [CrossRef]

30. William Allwood, J.; Ellis, D.I.; Heald, J.K.; Goodacre, R.; Mur, L.A.J. Metabolomic approaches reveal that phosphatidic and phosphatidyl glycerol phospholipids are major discriminatory non-polar metabolites in responses by *Brachypodium distachyon* to challenge by *Magnaporthe grisea*. *Plant J.* **2006**, *46*, 351–368. [CrossRef]

31. Jansen, J.J.; Allwood, J.W.; Marsden-Edwards, E.; van der Putten, W.H.; Goodacre, R.; van Dam, N.M. Metabolomic analysis of the interaction between plants and herbivores. *Metabolomics* **2008**, *5*, 150. [CrossRef]

32. Allwood, J.W.; Heald, J.; Lloyd, A.J.; Goodacre, R.; Mur, L.A.J. Separating the Inseparable: The Metabolomic Analysis of Plant–Pathogen Interactions. In *Plant Metabolomics: Methods and Protocols*; Hardy, N.W., Hall, R.D., Eds.; Humana Press: Totowa, NJ, USA, 2012; pp. 31–49.

33. Gakière, B.; Fernie, A.R.; Pétriacq, P. More to NAD+ than meets the eye: A regulator of metabolic pools and gene expression in *Arabidopsis*. *Free Radic. Biol. Med.* **2018**, *122*, 86–95. [CrossRef]

34. Pétriacq, P.; Ton, J.; Patrit, O.; Tcherkez, G.; Gakière, B. NAD Acts as an Integral Regulator of Multiple Defense Layers. *Plant Physiol.* **2016**, *172*, 1465–1479. [CrossRef] [PubMed]

35. Baxter, A.; Mittler, R.; Suzuki, N. ROS as key players in plant stress signalling. *J. Exp. Bot.* **2014**, *65*, 1229–1240. [CrossRef]

36. Zhang, Y.; Krahnert, I.; Bolze, A.; Gibon, Y.; Fernie, A.R. Adenine Nucleotide and Nicotinamide Adenine Dinucleotide Measurements in Plants. *Curr. Protoc. Plant Biol.* **2020**, *5*, e20115. [CrossRef]

37. Guérard, F.; Pétriacq, P.; Gakière, B.; Tcherkez, G. Liquid chromatography/time-of-flight mass spectrometry for the analysis of plant samples: A method for simultaneous screening of common cofactors or nucleotides and application to an engineered plant line. *Plant Physiol. Biochem.* **2011**, *49*, 1117–1125. [CrossRef]

38. Straube, H.; Niehaus, M.; Zwittian, S.; Witte, C.-P.; Herde, M. Enhanced nucleotide analysis enables the quantification of deoxynucleotides in plants and algae revealing connections between nucleoside and deoxynucleoside metabolism. *Plant Cell* **2021**, *33*, 270–289. [CrossRef]

39. Vishwakarma, A.; Wany, A.; Pandey, S.; Bulle, M.; Kumari, A.; Kishorekumar, R.; Igamberdiev, A.U.; Mur, L.A.J.; Gupta, K.J. Current approaches to measure nitric oxide in plants. *J. Exp. Bot.* **2019**, *70*, 4333–4343. [CrossRef] [PubMed]

40. Estavillo, G.M.; Crisp, P.A.; Pornsiriwong, W.; Wirtz, M.; Collinge, D.; Carrie, C.; Giraud, E.; Whelan, J.; David, P.; Javot, H.; et al. Evidence for a SAL1-PAP Chloroplast Retrograde Pathway That Functions in Drought and High Light Signaling in Arabidopsis. *Plant Cell* **2011**, *23*, 3992–4012. [CrossRef] [PubMed]

41. Thevenet, D.; Pastor, V.; Baccelli, I.; Balmer, A.; Vallat, A.; Neier, R.; Glauser, G.; Mauch-Mani, B. The priming molecule β-aminobutyric acid is naturally present in plants and is induced by stress. *New Phytol.* **2017**, *213*, 552–559. [CrossRef]

42. Podlešáková, K.; Ugena, L.; Spíchal, L.; Doležal, K.; De Diego, N. Phytohormones and polyamines regulate plant stress responses by altering GABA pathway. *New Biotechnol.* **2019**, *48*, 53–65. [CrossRef] [PubMed]

43. Howe, G.A. Metabolic End Run to Jasmonate. *Nat. Chem. Biol.* **2018**, *14*, 109–110. [CrossRef]

44. Li, J.; Wang, X. Phospholipase D and phosphatidic acid in plant immunity. *Plant Sci.* **2019**, *279*, 45–50. [CrossRef] [PubMed]

45. Pétriacq, P.; de Bont, L.; Hager, J.; Didierlaurent, L.; Mauve, C.; Guérard, F.; Noctor, G.; Pelletier, S.; Renou, J.-P.; Tcherkez, G.; et al. Inducible NAD overproduction in Arabidopsis alters metabolic pools and gene expression correlated with increased salicylate content and resistance to Pst-AvrRpm1. *Plant J.* **2012**, *70*, 650–665. [CrossRef] [PubMed]

46. Kostyuk, A.I.; Panova, A.S.; Kokova, A.D.; Kotova, D.A.; Maltsev, D.I.; Podgorny, O.V.; Belousov, V.V.; Bilan, D.S. In Vivo Imaging with Genetically Encoded Redox Biosensors. *Int. J. Mol. Sci.* **2020**, *21*, 8164. [CrossRef]

47. De Col, V.; Fuchs, P.; Nietzel, T.; Elsässer, M.; Voon, C.P.; Candeo, A.; Seeliger, I.; Fricker, M.D.; Grefen, C.; Møller, I.M.; et al. ATP sensing in living plant cells reveals tissue gradients and stress dynamics of energy physiology. *eLife* **2017**, *6*, e26770. [CrossRef] [PubMed]

48. Nietzel, T.; Elsässer, M.; Ruberti, C.; Steinbeck, J.; Ugalde, J.M.; Fuchs, P.; Wagner, S.; Ostermann, L.; Moseler, A.; Lemke, P.; et al. The fluorescent protein sensor roGFP2-Orp1 monitors in vivo H₂O₂ and thiol redox integration and elucidates intracellular H₂O₂ dynamics during elicitor-induced oxidative burst in Arabidopsis. *New Phytol.* **2019**, *221*, 1649–1664. [CrossRef]

49. Steinbeck, J.; Fuchs, P.; Negroni, Y.L.; Elsässer, M.; Lichtenauer, S.; Stockdreher, Y.; Feitosa-Araujo, E.; Kroll, J.B.; Niemeier, J.-O.; Humberg, C.; et al. In Vivo NADH/NAD+ Biosensing Reveals the Dynamics of Cytosolic Redox Metabolism in Plants. *Plant Cell* **2020**, *32*, 3324–3345. [CrossRef]

50. Lim, S.-L.; Voon, C.P.; Guan, X.; Yang, Y.; Gardeström, P.; Lim, B.L. In planta study of photosynthesis and photorespiration using NADPH and NADH/NAD+ fluorescent protein sensors. *Nat. Commun.* **2020**, *11*, 3238. [CrossRef]
51. Allwood, J.W.; De Vos, R.C.H.; Moing, A.; Deborde, C.; Erban, A.; Kopka, J.; Goodacre, R.; Hall, R.D. Chapter sixteen—Plant Metabolomics and Its Potential for Systems Biology Research: Background Concepts, Technology, and Methodology. In *Methods in Enzymology*; Jameson, D., Verma, M., Westerhoff, H.V., Eds.; Academic Press: Cambridge, MA, USA, 2011; Volume 500, pp. 299–336.

52. Lu, Q.; Wang, Y.; Xiong, F.; Hao, X.; Zhang, X.; Li, N.; Wang, L.; Zeng, J.; Yang, Y.; Wang, X. Integrated transcriptomic and metabolomic analyses reveal the effects of callose deposition and multihormone signal transduction pathways on the tea plant-Colletotrichum camelliae interaction. *Sci. Rep.* **2020**, *10*, 12858. [CrossRef]

53. Macel, M.; Visschers, I.G.S.; Peters, J.L.; van Dam, N.M.; de Graaf, R.M. High Concentrations of Very Long Chain Leaf Wax Alkanes of Thrips Susceptible Pepper Accessions (*Capsicum* spp.). *J. Chem. Ecol.* **2020**, *46*, 1082–1089. [CrossRef]

54. Lewandowska, M.; Keyl, A.; Feussner, I. Wax biosynthesis in response to danger: Its regulation upon abiotic and biotic stress. *New Phytol.* **2020**, *227*, 698–713. [CrossRef]

55. Hunter, K.; Kimura, S.; Rokka, A.; Tran, H.C.; Toyota, M.; Kukkonen, J.P.; Wrzaczek, M. CRK2 Enhances Salt Tolerance by Regulating Callose Deposition in Connection with PLDα1. *Plant Physiol.* **2019**, *180*, 2004–2021. [CrossRef]

56. Buchberger, A.R.; DeLaney, K.; Johnson, J.; Li, L. Mass Spectrometry Imaging: A Review of Emerging Advancements and Future Insights. *Anal. Chem.* **2018**, *90*, 240–265. [CrossRef]

57. Lima, C.; Muhamadali, H.; Goodacre, R. The Role of Raman Spectroscopy within Quantitative Metabolomics. *Annu. Rev. Anal. Chem.* **2021**, *14*, 323–345. [CrossRef] [PubMed]

58. Shepherd, T.; Robertson, G.W.; Griffiths, D.W.; Birch, A.N.E. Epicuticular wax ester and triacylglycerol composition in relation to aphid infestation and resistance in red raspberry (*Rubus idaeus* L.). *Phytochemistry* **1999**, *52*, 1255–1267. [CrossRef]

59. Wu, X.; Shi, X.; Bai, M.; Chen, Y.; Li, X.; Qi, K.; Cao, P.; Li, M.; Yin, H.; Zhang, S. Transcriptomic and Gas Chromatography–Mass Spectrometry Metabolomic Profiling Analysis of the Epidermis Provides Insights into Cuticular Wax Regulation in Developing ‘Yuluxiang’ Pear Fruit. *J. Agric. Food Chem.* **2019**, *67*, 8319–8331. [CrossRef] [PubMed]

60. Kulkarni, P.; Dost, M.; Bulut, Ö.D.; Welle, A.; Böcker, S.; Boland, W.; Svatoš, A. Secondary ion mass spectrometry imaging and multivariate data analysis reveal co-aggregation patterns of *Populus trichocarpa* leaf surface compounds on a micrometer scale. *Plant J.* **2018**, *93*, 193–206. [CrossRef] [PubMed]

61. Touw, A.J.; Verdecia Mogena, A.; Maedicke, A.; Sontowski, R.; van Dam, N.M.; Tsunoda, T. Both Biosynthesis and Transport Are Involved in Glucosinolate Accumulation During Root-Herbivory in *Brassica rapa*. *Front. Plant Sci.* **2020**, *10*, 1653. [CrossRef]

62. Chowański, S.; Adamski, Z.; Marciniak, P.; Rosiński, G.; Büyükgüzel, E.; Büyükgüzel, K.; Falabella, P.; Scrano, L.; Ventrella, E.; Lelario, F.; et al. A Review of Bioinsecticidal Activity of Solanaceae Alkaloids. *Toxins* **2016**, *8*, 60. [CrossRef] [PubMed]

63. Okamoto, H.; Ducreux, L.J.M.; Allwood, J.W.; Hedley, P.E.; Wright, A.; Gururajan, V.; Terry, M.J.; Taylor, M.A. Light Regulation of Chlorophyll and Glycoalkaloid Biosynthesis During Tuber Greening of Potato, *S. tuberosum*. *Front. Plant Sci.* **2020**, *11*, 753. [CrossRef]

64. Zhou, S.; Richter, A.; Jander, G. Beyond Defense: Multiple Functions of Benzoxazinoids in Maize Metabolism. *Plant Cell Physiol.* **2018**, *59*, 1528–1537. [CrossRef]

65. Bolarinwa, I.F.; Oke, M.O.; Olaniyan, S.A.; Ajala, A.S. A Review of Cyanogenic Glycosides in Edible Plants. In *Toxicology: New Aspects to This Scientific Conundrum*; Soloneski, S., Larramendy, M., Eds.; Intech Open: London, UK, 2016.

66. Muzashvili, T.; Moniuszko-Szajwaj, B.; Pecio, L.; Oleszek, W.; Stochmal, A. Ultraperformance Liquid Chromatography Tandem Mass Spectrometry Determination of Cyanogenic Glucosides in Trifolium Species. *J. Agric. Food Chem.* **2014**, *62*, 1777–1782. [CrossRef]

67. Sun, Z.; Zhang, K.; Chen, C.; Wu, Y.; Tang, Y.; Georgiev, M.I.; Zhang, X.; Lin, M.; Zhou, M. Biosynthesis and regulation of cyanogenic glycoside production in forage plants. *Appl. Microbiol. Biotechnol.* **2018**, *102*, 9–16. [CrossRef] [PubMed]

68. Treutter, D. Significance of flavonoids in plant resistance: A review. *Environ. Chem. Lett.* **2006**, *4*, 147. [CrossRef]

69. Lei, Z.; Sumner, B.W.; Bhatia, A.; Sarma, S.J.; Sumner, L.W. UHPLC-MS Analyses of Plant Flavonoids. *Curr. Protoc. Plant Biol.* **2019**, *4*, e20085. [CrossRef] [PubMed]

70. Wen, W.; Alseekh, S.; Fernie, A.R. Conservation and diversification of flavonoid metabolism in the plant kingdom. *Curr. Opin. Plant Biol.* **2020**, *55*, 100–108. [CrossRef]

71. Bhuiyan, N.H.; Selvaraj, G.; Wei, Y.; King, J. Role of lignification in plant defense. *Plant Signal. Behav.* **2009**, *4*, 158–159. [CrossRef]

72. Ali, M.B.; McNear, D.H. Induced transcriptional profiling of phenylpropanoid pathway genes increased flavonoid and lignin content in Arabidopsis leaves in response to microbial products. *BMC Plant Biol.* **2014**, *14*, 84. [CrossRef] [PubMed]

73. Uarrota, V.G.; Stefen, D.L.V.; Leolato, L.S.; Gindri, D.M.; Nerling, D. Revisiting Carotenoids and Their Role in Plant Stress Responses: From Biosynthesis to Plant Signaling Mechanisms during Stress. In *Antioxidants and Antioxidant Enzymes in Higher Plants*; Gupta, D.K., Palma, J.M., Corpas, F.J., Eds.; Springer International Publishing: Cham, Switzerland, 2018; pp. 207–232.

74. Riedlmeier, M.; Ghirardo, A.; Wenig, M.; Knappe, C.; Koch, K.; Georgii, E.; Dey, S.; Parker, J.E.; Schnitzler, J.-P.; Vlot, A.C. Monoterpenes Support Systemic Acquired Resistance within and between Plants. *Plant Cell* **2017**, *29*, 1440–1459. [CrossRef]

75. Frank, L.; Wenig, M.; Ghirardo, A.; van der Krol, A.; Vlot, A.C.; Schnitzler, J.-P.; Rosenkranz, M. Isoprene and β-caryophyllene confer plant resistance via different plant internal signalling pathways. *Plant Cell Environ.* **2021**, *44*, 1151–1164. [CrossRef]

76. Petry, F.C.; Mercadante, A.Z. Composition by LC-MS/MS of New Carotenoid Esters in Mango and Citrus. *J. Agric. Food Chem.* **2016**, *64*, 8207–8224. [CrossRef] [PubMed]

77. Barros, J.; Dixon, R.A. Plant Phenylalanine/Tyrosine Ammonia-lyases. *Trends Plant Sci.* **2020**, *25*, 66–79. [CrossRef]
78. Mechri, B.; Tekaya, M.; Cheheb, H.; Hammami, M. Determination of Mannitol Sorbitol and Myo-Inositol in Olive Tree Roots and Rhizospheric Soil by Gas Chromatography and Effect of Severe Drought Conditions on Their Profiles. *J. Chromatogr. Sci.* **2015**, *53*, 1631–1638. [CrossRef] [PubMed]

79. Zulfiqar, F.; Akram, N.A.; Ashraf, M. Osmoprotection in plants under abiotic stresses: New insights into a classical phenomenon. *Planta* **2019**, *251*, 3. [CrossRef]

80. William Allwood, J.; Clarke, A.; Goodacre, R.; Mur, L.A.J. Dual metabolomics: A novel approach to understanding plant-pathogen interactions. *Phytochemistry* **2010**, *71*, 590–597. [CrossRef] [PubMed]

81. Lewis, L.A.; Polanski, K.; de Torres-Zabala, M.; Jayaraman, S.; Bowden, L.; Moore, J.; Penfold, C.A.; Jenkins, D.J.; Hill, C.; Baxter, L.; et al. Transcriptional Dynamics Driving MAMP-Triggered Immunity and Pathogen Effector-Mediated Immunosuppression in Arabidopsis Leaves Following Infection with *Pseudomonas syringae* pv tomato DC3000. *Plant Cell* **2015**, *27*, 3038–3064. [CrossRef] [PubMed]

82. De Torres-Zabala, M.; Truman, W.; Bennett, M.H.; Lafforgue, G.; Mansfield, J.W.; Rodriguez Egea, P.; Bögre, L.; Grant, M. *Pseudomonas syringae* pv. tomato hijacks the *Arabidopsis abscisic* acid signalling pathway to cause disease. *EMBO J.* **2007**, *26*, 1434–1443. [CrossRef]

83. De Torres Zabala, M.; Littlejohn, G.; Jayaraman, S.; Studholme, D.; Bailey, T.; Lawson, T.; Tillich, M.; Licht, D.; Bölter, B.; Delfino, L.; et al. Chloroplasts play a central role in plant defence and are targeted by pathogen effectors. *Nat. Plants* **2015**, *1*, 15074. [CrossRef]

84. Lim, E.-K.; Doucet, C.J.; Li, Y.; Elias, L.; Worrall, D.; Spencer, S.P.; Ross, J.; Bowles, D.J. The Activity of Arabidopsis Glycosyltransferases toward Salicylic Acid, 4-Hydroxybenzoic Acid, and Other Benzoates. *J. Biol. Chem.* **2002**, *277*, 586–592. [CrossRef]

85. Dean, J.V.; Delaney, S.P. Metabolism of salicylic acid in wild-type, ugt74f1 and ugt74f2 glucosyltransferase mutants of *Arabidopsis thaliana*. *Physiol. Plant.* **2008**, *132*, 417–425. [CrossRef] [PubMed]

86. Chen, L.; Wang, W.-S.; Wang, T.; Meng, X.-F.; Chen, T.-T.; Huang, X.-X.; Li, Y.-J.; Hou, B.-K. Methyl Salicylate Glucosylation Regulates Plant Defense Signaling and Systemic Acquired Resistance. *Plant Physiol.* **2019**, *180*, 2167–2181. [CrossRef]

87. Návarová, H.; Bernsdorff, F.; Döring, A.-C.; Zeier, J. Pipecolic Acid, an Endogenous Mediator of Defense Amplification and Priming, Is a Critical Regulator of Inducible Plant Immunity. *Plant Cell* **2012**, *24*, 5123–5141. [CrossRef] [PubMed]

88. Hartmann, M.; Zeier, T.; Bernsdorff, F.; Reichel-Deland, V.; Kim, D.; Hohmann, M.; Scholten, N.; Schuck, S.; Bräutigam, A.; Hölzel, T.; et al. Flavin Monooxygenase-Generated N-Hydroxypipecolic Acid Is a Critical Element of Plant Systemic Immunity. *Cell* **2018**, *173*, 456–469. [CrossRef] [PubMed]

89. Chen, Y.-C.; Holmes, E.C.; Rajniak, J.; Kim, J.-G.; Tang, S.; Fischer, C.R.; Mudgett, M.B.; Sattely, E.S. N-hydroxy-pipecolic acid is a mobile metabolite that induces systemic disease resistance in *Arabidopsis*. *Proc. Natl. Acad. Sci. USA* **2018**, *115*, E4920. [CrossRef]

90. Holmes, E.C.; Chen, Y.-C.; Mudgett, M.B.; Sattely, E.S. Arabidopsis UGT76B1 glycosylates N-hydroxy-pipecolic acid and inactivates systemic acquired resistance in tomato. *Plant Cell* **2021**, *33*, 750–765. [CrossRef]

91. Mohnike, L.; Rekhter, D.; Huang, W.; Feussner, K.; Tian, H.; Herrfurth, C.; Zhang, Y.; Feussner, I. The glycosyltransferase UGT76B1 modulates N-hydroxy-pipecolic acid homeostasis and plant immunity. *Plant Cell* **2021**, *33*, 735–749. [CrossRef] [PubMed]

92. Huang, X.-X.; Zhu, G.-Q.; Liu, Q.; Chen, L.; Li, Y.-J.; Hou, B.-K. Modulation of Plant Salicylic Acid-Associated Immune Responses via Glycosylation of Dihydroxybenzoic Acids. *Plant Physiol.* **2018**, *176*, 3103–3119. [CrossRef]

93. Herrfurth, C.; Feussner, I. Quantitative Jasmonate Profiling Using a High-Throughput UPLC-NanoESI-MS/MS Method. In *Jasmonate in Plant Biology: Methods and Protocols*; Champion, A., Laplaze, L., Eds.; Springer: New York, NY, USA, 2020; pp. 169–187.

94. Tripathi, D.; Tanaka, K. A crosstalk between extracellular ATP and jasmonate signaling pathways for plant defense. *Plant Signal. Behav.* **2018**, *13*, e1432229. [CrossRef]

95. Song, C.J.; Steinebrunner, I.; Wang, X.; Stout, S.C.; Roux, S.J. Extracellular ATP Induces the Accumulation of Superoxide via NADPH Oxidases in Arabidopsis. *Plant Physiol.* **2006**, *140*, 1222–1232. [CrossRef]

96. Chen, D.; Cao, Y.; Li, H.; Kim, D.; Ahsan, N.; Thelen, J.; Stacey, G. Extracellular ATP elicits DORN1-mediated RBOHD phosphorylation to regulate stomatal aperture. *Nat. Commun.* **2017**, *8*, 2265. [CrossRef]

97. Chen, D.; Hao, F.; Mu, H.; Ahsan, N.; Thelen, J.J.; Stacey, G. S-acylation of P2K1 mediates extracellular ATP-induced immune signaling in Arabidopsis. *Nat. Commun.* **2021**, *12*, 2750. [CrossRef]

98. Heil, M.; Ibarra-Laclette, E.; Adame-Álvarez, R.M.; Martínez, O.; Ramirez-Chávez, E.; Molina-Torres, J.; Herrera-Estrella, L. How Plants Sense Wounds: Damaged-Self Recognition Is Based on Plant-Derived Elicitors and Induces Octadecanoid Signaling. *PLoS ONE* **2012**, *7*, e30537. [CrossRef] [PubMed]

99. Tripathi, D.; Zhang, T.; Koo, A.J.; Stacey, G.; Tanaka, K. Extracellular ATP Acts on Jasmonate Signaling to Reinforce Plant Defense. *Plant Physiol.* **2018**, *176*, 511–523. [CrossRef] [PubMed]

100. Wang, C.; Zhou, M.; Zhang, X.; Yao, J.; Zhang, Y.; Mou, Z. A lectin receptor kinase as a potential sensor for extracellular nicotinamide adenine dinucleotide in Arabidopsis thaliana. *eLife* **2017**, *6*, e25474. [CrossRef] [PubMed]

101. Wang, C.; Huang, X.; Li, Q.; Zhang, Y.; Li, J.-L.; Mou, Z. Extracellular pyridine nucleotides trigger plant systemic immunity through a lectin receptor kinase/BAK1 complex. *Nat. Commun.* **2019**, *10*, 4810. [CrossRef]

102. Wan, L.; Essuman, K.; Anderson, R.G.; Sasaki, Y.; Monteiro, F.; Chung, E.-H.; Osborne Nishimura, E.; DiAntonio, A.; Milbrandt, J.; Dangl, J.L.; et al. TIR domains of plant immune receptors are NAD⁺-cleaving enzymes that promote cell death. *Science* **2019**, *365*, 799. [CrossRef]
103. Horsefield, S.; Burdett, H.; Zhang, X.; Manik, M.K.; Shi, Y.; Chen, J.; Qi, T.; Gilley, J.; Lai, J.-S.; Rank, M.X.; et al. NAD⁺ cleavage activity by animal and plant TIR domains in cell death pathways. *Science* **2019**, *365*, 793. [CrossRef]

104. Tarkowski, Ł.P.; Signorelli, S.; Höfte, M. γ-Aminobutyric acid and related amino acids in plant immune responses: Emerging mechanisms of action. *Plant Cell Environ.* **2020**, *43*, 1103–1116. [CrossRef]

105. Xian, L.; Yu, G.; Wei, Y.; Rufian, J.S.; Li, Y.; Zhuang, H.; Xue, H.; Morcillo, R.J.L.; Macho, A.P. A Bacterial Effector Protein Hijacks Plant Metabolism to Support Pathogen Nutrition. *Cell Host Microbe* **2020**, *28*, 548–557. [CrossRef]

106. Streubel, J.; Pesce, C.; Hutin, M.; Koebnik, R.; Boch, J.; Szurek, B. Five phylogenetically close rice SWEET genes confer TAL effector-mediated susceptibility to *Xanthomonas oryzae* pv. oryzae. *New Phytol.* **2013**, *200*, 808–819. [CrossRef]

107. David, L.; Kang, J.; Chen, S. Untargeted Metabolomics of Arabidopsis Stomatal Immunity. In *Arabidopsis Protocols*; Sanchez-Serrano, J.J., Salinas, J., Eds.; Springer: New York, NY, USA, 2021; pp. 413–424.

108. Chen, F.; Ma, R.; Chen, X.-L. Advances of Metabolomics in Fungal Pathogen-Plant Interactions. *Metabolites* **2019**, *9*, 169. [CrossRef] [PubMed]

109. Balmer, D.; Flors, V.; Glauser, G.; Mauch-Mani, B. Metabolomics of cereals under biotic stress: Current knowledge and techniques. *Front. Plant Sci.* **2013**, *4*. [CrossRef]

110. Parker, D.; Beckmann, M.; Zubair, H.; Enot, D.P.; Caracuel-Rios, Z.; Overy, D.P.; Snowdon, S.; Talbot, N.J.; Draper, J. Metabolomic analysis reveals a common pattern of metabolic re-programming during invasion of three host plant species by *Magnaporthe grisea*. *Plant J.* **2009**, *59*, 723–737. [CrossRef]

111. Martin-Urdiroz, M.; Oses-Ruiz, M.; Ryder, L.S.; Talbot, N.J. Investigating the biology of plant infection by the rice blast fungus *Magnaporthe oryzae*. *Fungal Genet. Biol.* **2016**, *90*, 61–68. [CrossRef]

112. Yan, X.; Talbot, N.J. Investigating the cell biology of plant infection by the rice blast fungus *Magnaporthe oryzae*. *Curr. Opin. Microbiol.* **2016**, *34*, 147–153. [CrossRef]

113. Leesutthiphonchai, W.; Vu, A.L.; Ah-Fong, A.M.V.; Judelson, H.S. How Does *Phytophthora infestans* Evade Control Efforts? Modern Insight Into the Late Blight Disease. *Phytopathology* **2018**, *108*, 916–924. [CrossRef]

114. Whisson, S.C.; Boevink, P.C.; Wang, S.; Birch, P.R.J. The cell biology of late blight disease. *Curr. Opin. Microbiol.* **2016**, *34*, 127–135. [CrossRef] [PubMed]

115. Van Kan, J.A.L. Licensed to kill: The lifestyle of a necrotrophic plant pathogen. *Trends Plant Sci.* **2006**, *11*, 247–253. [CrossRef]

116. Bai, G.; Shaner, G. Management and resistance in wheat and barley to fusarium head blight. *Annu. Rev. Phytopathol.* **2004**, *42*, 135–161. [CrossRef] [PubMed]

117. Mur, L.A.J.; Kenton, P.; Lloyd, A.J.; Ougham, H.; Prats, E. The hypersensitive response; the centenary is upon us but how much do we know? *J. Exp. Bot.* **2008**, *59*, 501–520. [CrossRef] [PubMed]

118. Liu, X.-H.; Liang, S.; Wei, Y.-Y.; Zhu, X.-M.; Li, L.; Liu, P.-P.; Zheng, Q.-X.; Zhou, H.-N.; Zhang, Y.; Mao, L.-J.; et al. Metabolomics Analysis Identifies Sphingolipids as Key Signaling Moieties in Appressorium Morphogenesis and Function in Magnaporthe oryzae. *mBio* **2019**, *10*, e01467-19. [CrossRef] [PubMed]

119. Jones, O.A.H.; Maguire, M.L.; Griffin, J.L.; Jung, Y.-H.; Shibato, J.; Rakwal, R.; Agrawal, G.K.; Jwa, N.-S. Using metabolic profiling to assess plant-pathogen interactions: An example using rice (*Oryza sativa*) and the blast pathogen *Magnaporthe grisea*. *Eur. J. Plant Pathol.* **2011**, *129*, 539–554. [CrossRef]

120. Tezuka, D.; Kawamata, A.; Kato, H.; Saburi, W.; Mori, H.; Imai, R. The rice ethylene response factor OsERF83 positively regulates disease resistance to *Magnaporthe oryzae*. *Plant Physiol. Biochem.* **2019**, *135*, 263–271. [CrossRef]

121. Rodenburg Sander, Y.A.; Seidl Michael, F.; Judelson Howard, S.; Vu Andrea, L.; Govers, F.; de Ridder, D.; McDowell, J.; Gillian Turgeon, B. Metabolic Model of the Phytophthora infestans-Tomato Interaction Reveals Metabolic Switches during Host Colonization. *mBio* **2019**, *10*, e00454-19. [CrossRef]

122. Rodenburg, S.Y.A.; Seidl, M.F.; de Ridder, D.; Govers, F. Genome-wide characterization of *Phytophthora infestans* metabolism: A systems biology approach. *Mol. Plant Pathol.* **2018**, *19*, 1403–1413. [CrossRef]

123. Fang, C.; Luo, J. Metabolic GWAS-based dissection of genetic bases underlying the diversity of plant metabolism. *Plant J.* **2019**, *97*, 91–100. [CrossRef]

124. Paranidharan, V.; Abu-Nada, Y.; Hamzehzarghani, H.; Kushalappa, A.C.; Mamer, O.; Dion, Y.; Rioux, S.; Comeau, A.; Choiniere, L. Resistance-related metabolites in wheat against *Fusarium graminearum* and the virulence factor deoxynivalenol (DON). *Botany* **2008**, *86*, 1168–1179. [CrossRef]

125. Bollina, V.; Kumaraswamy, G.K.; Kushalappa, A.C.; Choo, T.M.; Dion, Y.; Rioux, S.; Faubert, D.; Hamzehzarghani, H. Mass spectrometry-based metabolomics application to identify quantitative resistance-related metabolites in barley against *Fusarium head blight*. *Mol. Plant Pathol.* **2010**, *11*, 769–782. [CrossRef]

126. Gunnaiah, R.; Kushalappa, A.C.; Duggavathi, R.; Fox, S.; Somers, D.J. Integrated Metabolo-Proteomic Approach to Decipher the Mechanisms by Which Wheat QTL (Fhb1) Contributes to Resistance against *Fusarium graminearum*. *PLoS ONE* **2012**, *7*, e40695. [CrossRef]

127. Cajka, T.; Vaclavikova, M.; Dzuman, Z.; Vaclavik, L.; Ovesna, J.; Hajslova, J. Rapid LC-MS-based metabolomics method to study the *Fusarium infection* of barley. *J. Sep. Sci.* **2014**, *37*, 912–919. [CrossRef] [PubMed]

128. Kumaraswamy, G.K.; Kushalappa, A.C.; Choo, T.M.; Dion, Y.; Rioux, S. Differential metabolic response of barley genotypes, varying in resistance, to trichothecene-producing and -nonproducing (tri5–) isolates of *Fusarium graminearum*. *Plant Pathol. 2012*, *61*, 509–521. [CrossRef]
129. Hu, Z.; Chang, X.; Dai, T.; Li, L.; Liu, P.; Wang, G.; Liu, P.; Huang, Z.; Liu, X. Metabolic Profiling to Identify the Latent Infection of Strawberry by Botrytis cinerea. *Evol. Bioinform. 2019*, *15*, 1176934319838518. [CrossRef]

130. Bhandari, D.R.; Wang, Q.; Friedt, W.; Spengler, B.; Gottwald, S.; Rompp, A. High resolution mass spectrometry imaging of plant tissues: Towards a plant metabolite atlas. *Analyst* **2015**, *140*, 7696–7709. [CrossRef]

131. Heyman, H.M.; Dubery, I.A. The potential of mass spectrometry imaging in plant metabolomics: A review. *Phytochem. Rev.* **2016**, *15*, 297–316. [CrossRef]

132. Cornett, D.S.; Reyzer, M.L.; Chaurand, P.; Caprioli, R.M. MALDI imaging mass spectrometry: Molecular snapshots of biochemical systems. *Nat. Methods* **2007**, *4*, 828–833. [CrossRef]

133. Shroff, R.; Vergara, F.; Muck, A.; Svatoš, A.; Gershenzon, J. Nonuniform distribution of glucosinolates in *Arabidopsis thaliana* leaves has important consequences for plant defense. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 6196. [CrossRef] [PubMed]

134. Etalo, D.W.; De Vos, R.C.H.; Joosten, M.H.A.J.; Hall, R.D. Spatially Resolved Plant Metabolomics: Some Potentials and Limitations of Laser-Ablation Electrospray Ionization Mass Spectrometry Metabolite Imaging. *Plant Physiol.* **2015**, *169*, 1424–1435. [CrossRef]

135. Mitchell, C.; Brennan, R.M.; Graham, J.; Karley, A.J. Plant Defense against Herbivorous Pests: Exploiting Resistance and Tolerance Traits for Sustainable Crop Protection. *Front. Plant Sci.* **2016**, *7*, 1132. [CrossRef] [PubMed]

136. Allema, B.; Hoogendoorn, M.; van Beek, J.; Leendertse, P. Neonicotinoids in European agriculture. Main applications, main crops and scope for alternatives. In *CRa Advice*; CLM Research and Advice: Culemborg, NL, USA, 2017.

137. Visschers, I.G.S.; Peters, J.L.; van de Vondervoort, J.A.H.; Hoogveld, R.H.M.; van Dam, N.M. Thrips Resistance Screening Is Coming of Age: Leaf Position and Ontogeny Are Important Determinants of Leaf-Based Resistance in Pepper. *Front. Plant Sci.* **2019**, *10*, 510. [CrossRef] [PubMed]

138. Omena-Garcia, R.P.; de Ávila Silva, L.; Vital, C.E.; Araújo, W.L.; Ramos, H.J.O.; Nunes-Nesi, A. Identification of metabolite traits from the current metabolomic approaches. *Theor. Exp. Plant Physiol.* **2019**, *31*, 1–19. [CrossRef]

139. Wahyuni, Y.; Stahl-Hermes, V.; Ballester, A.-R.; de Vos, R.C.H.; Voorrips, R.E.; Maharijaya, A.; Molthoff, J.; Zamora, M.V.; Sudarmonowati, E.; Arisi, A.C.M.; et al. Genetic mapping of semi-polar metabolites in pepper fruits (*Capsicum* sp.): Towards unravelling the molecular regulation of flavonoid quantitative trait loci. *Mol. Breed.* **2014**, *33*, 503–518. [CrossRef] [PubMed]

140. Ran, F.A.; Hsu, P.D.; Wright, J.; Agarwala, V.; Scott, D.A.; Zhang, F. Genome engineering using the CRISPR-Cas9 system. *Nat. Protoc.* **2013**, *8*, 2281–2308. [CrossRef]

141. Watson, A.; Ghosh, S.; Williams, M.J.; Cuddy, W.S.; Simmonds, J.; Rey, M.-D.; Asyraf Md Hatta, M.; Hinchliffe, A.; Steed, A.; Reynolds, D.; et al. Speed breeding is a powerful tool to accelerate crop research and breeding. *Nat. Plants* **2018**, *4*, 23–29. [CrossRef]

142. Björkman, M.; Klingen, I.; Birch, A.N.E.; Bones, A.M.; Bruce, T.J.A.; Johansen, T.J.; Meadow, R.; Mølmann, J.; Seljåsen, R.; Smart, L.E.; et al. Phytochemicals of Brassicaceae in plant protection and human health—Influences of climate, environment and agronomic practice. *Phytochemistry* **2011**, *72*, 538–556. [CrossRef] [PubMed]

143. Baaij, B.M.; Kim, H.K.; Grosser, K.; Worrich, A.; de Jong, T.J. Slug herbivory on hybrids of the crop Brassica napus and its wild relative *B. rapa*. *Basic Appl. Ecol.* **2018**, *31*, 52–60. [CrossRef]

144. Xu, S.Y.; Weng, J.-K. Climate change shapes the future evolution of plant metabolism. *Adv. Genet.* **2020**, *1*, e10022. [CrossRef]

145. Alawiye, T.T.; Babalola, O.O. Metagenomic Insight into the Community Structure and Functional Genes in the Sunflower Rhizosphere Microbiome. *Agriculture* **2021**, *11*, 167. [CrossRef]

146. Williams, A.; de Vries, F.T. Plant root exudation under drought: Implications for ecosystem functioning. *New Phytol.* **2020**, *225*, 1899–1905. [CrossRef]

147. Vaughan, M.M.; Huffaker, A.; Schmelz, E.A.; Dafoe, N.J.; Christensen, S.; Sims, J.; Martins, V.F.; Swerbilow, J.A.Y.; Romero, M.; Alborn, H.T.; et al. Effects of elevated [CO₂] on maize defence against mycotoxigenic *Fusarium verticillioides*. *Plant Cell Environ.* **2014**, *37*, 2691–2706. [CrossRef] [PubMed]

148. Mhamdi, A.; Noctor, G. High CO₂ Primes Plant Biotic Stress Defences through Redox-Linked Pathways. *Plant Physiol.* **2016**, *172*, 929–942. [CrossRef]

149. Williams, A.; Pétriacq, P.; Schwarzenbacher, R.E.; Beerling, D.J.; Ton, J. Mechanisms of glacial-to-future atmospheric CO₂ effects on plant immunity. *New Phytol.* **2018**, *218*, 752–761. [CrossRef] [PubMed]

150. Meng, L.; Zhang, T.; Geng, S.; Scott, P.B.; Li, H.; Chen, S. Comparative proteomics and metabolomics of JAZ7-mediated drought tolerance in Arabidopsis. *J. Proteom.* **2019**, *196*, 81–91. [CrossRef]

151. You, J.; Zhang, Y.; Liu, A.; Li, D.; Wang, X.; Dossa, K.; Zhou, R.; Yu, J.; Zhang, Y.; Wang, L.; et al. Transcriptomic and metabolomic profiling of drought-tolerant and susceptible sesame genotypes in response to drought stress. *BMC Plant Biol.* **2019**, *19*, 267. [CrossRef] [PubMed]

152. Lim, C.W.; Baek, W.; Jung, J.; Kim, J.-H.; Lee, S.C. Function of ABA in Stomatal Defense against Biotic and Drought Stresses. *Int. J. Mol. Sci.* **2015**, *16*, 15251–15270. [CrossRef] [PubMed]

153. Vurukonda, S.S.K.P.; Vardharajula, S.; Shrivastava, M.; SkZ, A. Enhancement of drought stress tolerance in crops by plant growth promoting rhizobacteria. *Microbiol. Res.* **2016**, *184*, 13–24. [CrossRef] [PubMed]

154. Prasch, C.M.; Sonnewald, U. Simultaneous Application of Heat, Drought, and Virus to Arabidopsis Plants Reveals Significant Shifts in Signaling Networks. *Plant Physiol.* **2013**, *162*, 1849–1866. [CrossRef]

155. Christensen, S.A.; Santana, E.A.; Alborn, H.T.; Block, A.K.; Chamberlain, C.A. Metabolomics by UHPLC-HRMS reveals the impact of heat stress on pathogen-elicited immunity in maize. *Metabolomics* **2021**, *17*, 6. [CrossRef]
156. Cohen, S.P.; Leach, J.E. High temperature-induced plant disease susceptibility: More than the sum of its parts. *Curr. Opin. Plant Biol.* **2020**, *56*, 235–241. [CrossRef]

157. Williams, A.; Gamir, J.; Gravot, A.; Pétriacq, P. Untangling plant immune responses through metabolomics. In *Plant Metabolomics in Full Swing*; Elsevier: Amsterdam, The Netherlands, 2021; Volume 98, pp. 73–105.

158. Tsugawa, H.; Cajka, T.; Kind, T.; Ma, Y.; Higgins, B.; Ikeda, K.; Kanazawa, M.; VanderGheynst, J.; Fiehn, O.; Arita, M. MS-DIAL: Data-independent MS/MS deconvolution for comprehensive metabolome analysis. *Nat. Methods* **2015**, *12*, 523–526. [CrossRef]

159. Domingo-Almenara, X.; Guijas, C.; Billings, E.; Montenegro-Burke, J.R.; Uritboonthai, W.; Aisporna, A.E.; Chen, E.; Benton, H.P.; Siuzdak, G. The METLIN small molecule dataset for machine learning-based retention time prediction. *Nat. Commun.* **2019**, *10*, 5811. [CrossRef]

160. da Silva, R.R.; Dorrestein, P.C.; Quinn, R.A. Illuminating the dark matter in metabolomics. *Proc. Natl. Acad. Sci. USA* **2015**, *112*, 12549. [CrossRef]

161. Domingo-Almenara, X.; Montenegro-Burke, J.R.; Ivanisevic, J.; Thomas, A.; Sidibé, J.; Teav, T.; Guijas, C.; Aisporna, A.E.; Rinehart, D.; Hoang, L.; et al. XCMS-MRM and METLIN-MRM: A cloud library and public resource for targeted analysis of small molecules. *Nat. Methods* **2018**, *15*, 681–684. [CrossRef] [PubMed]

162. Gomes, P.; Quirós-Guerrero, L.; Silva, C.; Pamplona, S.; Boutin, J.A.; Eberlin, M.; Wolfender, J.-L.; Silva, M. Feature-Based Molecular Network-Guided Dereplication of Natural Bioactive Products from Leaves of *Stryphnodendron pulcherrimum* (Willd.) Hochr. *Metabolites* **2021**, *11*, 281. [CrossRef]

163. Rubio, B.; Fernandez, O.; Cosson, P.; Berton, T.; Caballero, M.; Lion, R.; Roux, F.; Bergelson, J.; Gibon, Y.; Schurdi-Levraud, V. Metabolic Profile Discriminates and Predicts Arabidopsis Susceptibility to Virus under Field Conditions. *Metabolites* **2021**, *11*, 230. [CrossRef] [PubMed]

164. Luna, E.; Flandin, A.; Cassan, C.; Prigent, S.; Chevanne, C.; Kadiri, C.F.; Gibon, Y.; Pétriacq, P. Metabolomics to Exploit the Primed Immune System of Tomato Fruit. *Metabolites* **2020**, *10*, 96. [CrossRef]

165. Walker, H.J. Metabolite imaging by mass spectrometry: A new discovery tool. In *Advances in Botanical Research*; Elsevier: Amsterdam, The Netherlands, 2021; Volume 98, pp. 271–296.

166. Kalmar, J.G.; Oh, Y.; Dean, R.A.; Muddiman, D.C. Investigating host-pathogen meta-metabolic interactions of *Magnaporthe oryzae* infected barley using infrared matrix-assisted laser desorption electrospray ionization mass spectrometry. *Anal. Bioanal. Chem.* **2020**, *412*, 139–147. [CrossRef] [PubMed]

167. Pinzón, A.; Barreto, E.; Bernal, A.; Achenie, L.; González Barrios, A.F.; Isea, R.; Restrepo, S. Computational models in plant-pathogen interactions: The case of Phytophthora infestans. *Theor. Biol. Med. Model.* **2009**, *6*, 24. [CrossRef]

168. Botero, K.; Restrepo, S.; Pinzón, A. A genome-scale metabolic model of potato late blight suggests a photosynthesis suppression mechanism. *BMC Genom.* **2018**, *19*, 863. [CrossRef] [PubMed]

169. Pang, Q.; Zhang, T.; Wang, Y.; Kong, W.; Guan, Q.; Yan, X.; Chen, S. Metabolomics of Early Stage Plant Cell–Microbe Interaction Using Stable Isotope Labeling. *Front. Plant Sci.* **2018**, *9*, 760. [CrossRef] [PubMed]

170. Ceranic, A.; Doppler, M.; Büschl, C.; Parich, A.; Xu, K.; Koutnik, A.; Bürstmayr, H.; Lemmens, M.; Schuhmacher, R. *Preparation of Uniformly Labelled <sup>13</sup>C- and <sup>15</sup>N-Plants Using Customised Growth Chambers*; Research Square: Durham, NC, USA, 2020.

171. Zander, M.; Lewsey, M.G.; Clark, N.M.; Yin, L.; Bartlett, A.; Saldierna Guzmán, J.P.; Hann, E.; Langford, A.E.; Jow, B.; Wise, A.; et al. Integrated multi-omics framework of the plant response to jasmonic acid. *Nat. Plants* **2020**, *6*, 290–302. [CrossRef]

172. Uthe, H.; van Dam, N.M.; Hervé, M.R.; Sorokina, M.; Peters, K.; Weinhold, A. Chapter Six—A practical guide to implementing metabolomics in plant ecology and biodiversity research. In *Advances in Botanical Research*; Pétriacq, P., Bouchereau, A., Eds.; Academic Press: Cambridge, MA, USA, 2021; Volume 98, pp. 163–203.