Metabolome profiling dissects the oat (
Plants are sessile organisms that must adapt and survive recurring exposure events to a multitude of environmental stressors (abiotic and biotic). Biotic stress includes attacks by numerous pathogens like bacteria, fungi, nematodes, herbivores and oomycetes . Plants have an elaborate and sophisticated detection and signalling system that allows them to recognise pathogen attacks quickly and initiate dynamic defence responses at a molecular level via a number of genes and transcription factors, including those from the phytohormonal pathway . These interactions and how plants respond to various stress combinations has been widely studied and has led to the development of the ‘Stress Combinations and their Interactions in Plants Database’ (SCIPDb; http://www.nipgr.ac.in/scipdb.php ) that provides data on the morpho-physio-biochemical (phenome) and molecular (transcriptome and metabolome) responses . Plants utilise a multitude of pre-formed barriers such as the cell wall and cuticle to ward off potential pathogens. To counteract plant defence mechanisms, phytopathogens have evolved specialised mechanisms to infiltrate and overthrow host plant defences. In response, plants readily produce antimicrobial compounds that aid in preventing pathogen entry, maintaining cellular integrity and provide additional structural support. Pathogens that overcome these barriers trigger signalling pathways in the potential host that lead to the expression of defence response genes. Plants therefore rely on receptors to recognise pathogen invasion and activate signal transduction pathways that involve a multitude of genes and their products . These signalling pathways activate calcium-dependent protein kinases, MAPKs, transcription factors, G-proteins, ubiquitin, and hormones that leads to various responses such as the hypersensitive response (HR), cell wall modification, stomatal closure, production of reactive oxygen species (ROS), specific proteins targeted at the invader (e.g., ß-glucananses and chitinases, defensins and protease inhibitors), or specific metabolites (e.g., phytoalexins) to protect against further infection . The Pseudomonas syringae ( Ps ) species complex contains some of the best studied plant pathogens, existing as more than 60 pathovars able to infect a wide range of plants . P . syringae utilises two types of secretion systems when infecting plants namely, the type II secretion system (T2SS) to transport proteins and the type III secretion system (T3SS) to deliver effectors into plant cells along with toxins that play an essential role in causing disease . Numerous phytotoxins have been discovered in Ps strains, including coronatine, syringolin A, syringomycin, syringopeptin, phaseolotoxin, and tabtoxin. Coronatine mimics jasmonoyl-isoleucine (JA-Ile) and is one of the most well-studied Ps toxins . Coronatine inhibits the activity of salicylic acid during stomatal closure, causing the stomata to reopen and allow pathogen invasion . Once effectors are delivered and recognised by resistance (R) proteins, effector-triggered immunity (ETI) is triggered, and an HR induced  . According to this scenario, effectors are ‘pathogen-derived Avr proteins that trigger resistance via activation of specific cognate host R proteins’. The latter are defined as ‘proteins that confers resistance by mediating direct or indirect recognition of a pathogen Avr protein’ . Nucleotide-binding domain and leucine-rich repeat-containing proteins (NLRs) are ubiquitous in plants and a constitute a major class of immune receptors. NLRs detect specific pathogen effectors through diverse mechanisms and different models have been proposed to describe the mechanism of detection. These include the direct -, guard -, decoy -, integrated decoy—and NLR-like models . Successful pathogens are able to overcome host plant defences by releasing/injecting effectors, however these are insufficiently potent in nonhost plant species. Effectors interfere with signal perception and transduction events within the host cell and may target transcription of specific defence-related genes . Very little is known about the genes that control whether effectors can inhibit basal defence or not, despite significant progress in the understanding of the molecular features of nonhost resistance to plant diseases. Effectors from virulent pathogens are successful due to their ability to overthrow basal resistance compared to nonhosts. According to recent research, nonhost resistance should be regarded as a form of basal resistance that is polygenetically inherited and shares a high resemblance and correlation with innate plant resistance to adaptive pathogens . Nonhost resistance can also be divided into type I and type II, based on the presence or absence of visual phenotypic changes. While type II nonhost resistance is linked to visual necrosis/cell death brought on by a HR, type I nonhost resistance does not cause any visual symptoms. In response to generic pathogen-derived elicitors like PAMPs, type I nonhost resistance often utilises passive or preformed barriers to activate defence mechanisms. Type I nonhost resistance is comparable to pattern-triggered immunity (PTI) and type II nonhost resistance to ETI. Despite some similarities between host and nonhost resistance, nonhost resistance is more complicated, and the mechanism of resistance might change depending on the pathogen and the plant species . Phytoanticipins also play an important role in preventing infections by nonadapted pathogens. The most well-known examples of these preformed chemical defences are glucosinolates in Arabidopsis and in oat, the avenacins . In view of the above, alterations in plant metabolism therefore have a significant impact on the fate of attempted infections. Accordingly, studying plant-microbe interactions such as those between healthy and infected or resistant and susceptible plants can reveal metabolomic changes in signalling pathways that are crucial for elucidating and determining the outcome of these interactions . Since metabolites directly reflect biochemical processes (coordinated by genes, mRNA transcripts and proteins) the examination, at both a qualitative and quantitative level, provides new understanding of the biochemical systems behind the phenotype while also creating a thorough metabolic profile of the specific system . Metabolites play a variety of roles in plant-pathogen interactions, including pathogen surveillance, signal transduction, enzyme control, cell-to-cell transmission, and antimicrobial activity . Plants produce primary and secondary (or specialised) metabolites with most of the metabolome consisting of primary metabolites that are crucial for plant growth and development. Secondary metabolites, on the other hand, promote interaction with the environment and aids in plant survival under threatening conditions such as exposure to biotic and abiotic stressors. The synthesis of these specialised metabolites requires precursors that branch off from primary metabolite pathways . The shikimate pathway is the first step in the biosynthesis of aromatic amino acids; when activated due to stress it produces tryptophan, tyrosine, and phenylalanine, which supports and boosts secondary metabolite production. Depending on the stress and environment, different metabolites accumulate in distinct plant organs. Phytoalexins, for example, exhibit antibacterial properties against phytopathogens and accumulate in high concentrations in leaves. In addition to their antibacterial effects, several of these metabolites help to build polymeric barriers against disease penetration  Although Ps has been widely studied, limited studies have been done on oat under biotic stress. Here, these model Ps bacteria were used to study host and nonhost defence mechanism as well as the PTI and ETI triggered by the respective pathovars. An untargeted metabolomics approach was therefore applied to qualitatively profile and identify as many metabolites involved in the P . syringae–Avena sativa L. interaction as possible. Monitoring metabolite levels and changes therein, can complement and corroborate transcriptome and/or proteomic data on plant-pathogen interactions, thus revealing pathogen attack mechanisms. Transcriptomic and proteomic analyses are used to investigate changes in messenger RNA and proteins, respectively. Metabolites can also influence the outcomes of these gene-centered pathways. Thus, combining metabolomics and other omics data adds new layers of information to studies of plant-pathogen interactions, such as identifying metabolites with antimicrobial properties, metabolomic profile differences between infected and non-infected plants, and pathogenic requirements for infection and colonisation . Therefore, the integration of other omics approaches or a multi-omics approach with metabolomics plays an important role in untangling the mechanisms underlying pathogen attacks . In this report, an untargeted metabolomics approach was used to study the underlying mechanisms of A . sativa L. treated with various pathovars of Ps . While some oat cultivars exhibit a higher level of disease resistance/tolerance, the molecular mechanisms underlying these interactions are still poorly understood . In addition to P . syringae pv. coronafaciens (pathogenic on oat), P . syringae pv. tabaci and P . syringae pv. tomato DC3000 were chosen as putative nonhost pathogens and the hrcC mutant of P . syringae pv. tomato DC3000 (lacking the T3SS) was chosen to represent only the PTI component of the oat response to infiltration with the Ps pathovars. A few studies have reported halo blight occurrence (caused by Ps-c ), highlighting the economic and crop loss implications caused by the pathogen across several geographically dispersed countries . By analysing the cellular and molecular responses between the plant and pathogen, sustainable means of combating disease could be developed and be particularly useful in breeding programs by identifying biomarkers and pathways associated with susceptibility, tolerance or resistance .   Oat plant cultivation Seeds of the Dunnart oat cultivar was procured (Agricol, Pretoria, South Africa) and chosen for further infection studies using an initial screening procedure, as described below. Seedlings were cultivated in 10 cm, 250 mL pots (±15 seeds per pot) with a pasteurized (90˚C) germination mixture soil (Culterra, Muldersdrift, South Africa) and watered twice weekly with a water-soluble chemical fertiliser (Multisol ‘N’, Culterra, Muldersdrift, South Africa). The were cultivated in greenhouse conditions with a 12 h/12 h (light/dark) cycle, a light intensity of approximately 80 μmol/m 2 /s, and a temperature of 25°C. The study was designed to track the cultivar’s response to bacterial infection over time (2–6 days post inoculation—d.p.i). Triplicate pots were grown for each time point and were cultivated under identical conditions. When the seedlings achieved 3-week maturity (three-leaf seedling stage), they were thinned out and selected for uniformity in developmental stage before treatment with the bacterial suspensions described in the following sections. The experimental design included three independent biological replicates. Preparation of Pseudomonas syringae pathovars Dr. W. Kriel (Starke Ayres Seeds, Bredell, South Africa) provided an isolate pathogenic on oat namely P . syringae , pv. coronafaciens ( Ps-c ) with sequence information available on GenBank (PP940110), P . syringae pv. tabaci ( Ps-t ) was obtained from Prof. Y. Ichinose, Okayama University, Japan, and P . syringae pv. tomato DC3000 and its hrcC mutant from Dr. B. Kemmerling, University of Tuebingen, Germany. Inqaba Biotechnical Industries (Pretoria, South Africa) performed 16S rRNA sequencing to validate the pathovar identities. Ps isolates were then cultured and maintained on nutrient agar (Merck, Modderfontein, South Africa). A colony was selected under sterile circumstances in a laminar flow cabinet and cultivated overnight at 28°C in nutrient broth (Merck, Modderfontein, South Africa) on an orbital shaking incubator. The overnight culture’s OD 600 was measured and diluted with 0.1% Tween 20 and phosphate buffered saline (PBS) to a value of around 0.3. The same dilution was used for non-inoculated nutrient broth as a control (applied to vehicle control plants). The pathogenicity of P . syringae strains on several oat cultivars was examined by drop-inoculation  on leaf segments (OD600 ≈0.1, 0.2, 0.3). The inoculated leaf segments were stored in a high-humidity container for 5 days in a regulated growth environment with a 12-hour light-dark cycle at 25°C and visual symptoms scored daily. ‘Dunnart’ was chosen for further metabolomic analysis as a cultivar with a tolerant reaction based on initial visual observation tests that demonstrated its ability to tolerate infection during the 5-day period. Inoculation of oat seedlings At the three-leaf growth stage (approximately 21 days post emergence), the plants were treated by spraying with the Ps-c , Ps-t , DC3000 and hrcC − bacterial suspensions (prepared in PBS with 0.1% Tween 20), diluted to OD 600 ≈0.3. The vehicle control (VC) plants were sprayed with a solution free of the bacteria and the healthy control (HC) groups were untreated (i.e., not sprayed with either solution), all grown under normal growth conditions. The 50 mL of either the inoculum (containing the bacteria) or the control (0.1% Tween 20 in PBS) solution was evenly sprayed onto the leaves of the treated and VC groups respectively. The plants were then incubated in darkness in an incubator for 1 h to provide 100% relative humidity. Following the 1 h incubation, the plants were removed and another 50 mL of either inoculum or control solution was applied to the treated and VC groups, respectively, and further incubated for 6 h. After incubation, the plants were then exposed again to the same initial conditions: with cycles of 12 h light/dark, a light intensity of 80 μmol/m 2 /s and temperature of 25°C. Post-treatment harvesting of plants was done for treated, VC and HC groups at 2, 4 and 6 d.p.i. by cutting the leaves and immediately snap freezing with liquid nitrogen to quench metabolic activity associated with possible wounding and handling of the tissue. Leaves were stored at −80°C until metabolite extraction. Luminol-based reactive oxygen species (ROS) assay ROS production was measured using a modified version of the luminometry method . Briefly, leaves were pressure infiltrated with the respective bacterial suspensions (OD 600 ≈0.3) or with H 2 O (control). The infiltrated sites were then excised and washed for 1 h in sterile water with agitation. The single leaf discs (4 mm in diameter) were then rinsed 3 times and individually placed into a 96-well microtitre plate with 20 μg/mL horseradish peroxidase (HRP) and 34 μg/mL of luminol (Sigma-Aldrich, St. Louis, MO, USA). The luminescence was then measured over a period of 30 min using a Synergy HT Biotek microplate reader (Biotek Instruments, Winooski, VT, USA). Peroxidase (POX) assay POX activity associated with the innate immune response was measured using a modified version of the method described by . Leaf discs (4 mm in diameter) were punched from leaves of the oat plants and washed with agitation for 1 h in 1 mL of a 1 x MS (Murashige and Skoog) salt solution to remove any background activity caused by cutting the leaves. The discs were then carefully transferred to individual wells in a 96-well plate using a small spatula to minimise damage. Each well then received 50 μL of 1 x MS solution with either H 2 O (control) or one of the respective overnight bacterial suspensions (OD 600 ≈0.3). The plates were then sealed with parafilm and incubated for 20 h with agitation at room temperature. The leaf discs were removed and 50 μL of 5-aminosalicylic acid (Sigma-Aldrich, St. Louis, MO, USA) (1 mg/mL) at pH 6.0 with 0.01% hydrogen peroxide (Barrs Pharmaceuticals, Cape Town, RSA) was added to each well using a multichannel pipette in order to minimise timing differences. The reaction was allowed to proceed for 3 min and then stopped by adding 20 μM NaOH with a multichannel pipette in the same order as before, again to minimise timing differences and allow equal time for the enzymatic reactions to proceed. The plates were then analysed using a Synergy HT Biotek microplate reader (Biotek Instruments, Winooski, VT, USA) in absorbance mode at 600 nm (A 600 nm). Metabolite extraction and sample preparation The harvested leaf material was quenched with liquid nitrogen before being crushed into powder with a pre-cooled mortar and pestle. The samples were weighed (1 g) into 50 mL Falcon tubes with addition of 80% cold (4°C) aqueous analytical grade methanol (Romil Chemistry, Cambridge, UK).at a 1:10 m/v ratio. The suspensions were then homogenised with a probe sonicator (Bandelin Sonopuls, Berlin, Germany) at 55% power for 10 sec per sample. To avoid cross-contamination, equipment was cleansed between each sample. The homogenates were centrifuged at 5100 x g for 20 min at 4°C in a benchtop centrifuge after which the supernatants were kept and concentrated by evaporating the methanol under vacuum to approximately 1 mL using a rotary evaporator set to 55°C. The concentrated samples were transferred to 2 mL microcentrifuge tubes and dried in a centrifugal evaporator under vacuum. The dried extracts were then reconstituted by dissolving in 500 μL of 50% aqueous methanol (MilliQ deionised water and LC-grade methanol (Romil, Cambridge, UK). The samples were subsequently filtered through nylon syringe filters (0.22 μm) into chromatography vials fitted with 500 μL inserts, capped, and kept at 4°C until analysis. Sample analyses using ultra-high-performance liquid chromatography (UHPLC) and quadrupole time-of-flight mass spectrometry (qTOF-MS) An Acquity UHPLC system (Waters Corporation, Manchester, UK) was used to analyse 2 μL of each sample, separated into its respective components using a binary solvent on an HSS T3 reverse-phase column (Waters Corporation, Billerica, MA, USA; 2.1 × 150 mm × 1.8 μm), thermostatted at 60°C. The solvents used were MilliQ water and acetonitrile (Romil Chemistry, Cambridge, UK), both containing 0.1% formic acid (Sigma-Aldrich, Munich, Germany) and 2.5% isopropanol (IPA, Romil, Cambridge, UK). The run was set to 30 min per injection with an elution gradient carried out via a binary solvent system consisting of 0.1% aqueous formic acid with 2.5% isopropanol (solvent A) and 0.1% formic acid and 2.5% isopropanol in acetonitrile (Romil, Cambridge, UK; solvent B) at a flow rate of 0.4 mL/min. The initial conditions were 95% A and 5% B and held for 1 min. A gradient was applied to change the chromatographic conditions to 10% A and 90% B at 25 min; and changed to 5% A and 95% B at 25.10 min. These conditions were held for 2 min and then changed to the initial conditions at 28 min. The analytical column was allowed to equilibrate for 2 min before each subsequent injection. Pooled quality control (QC) samples were also prepared to monitor the stability of the LC-MS system and assess the reliability and reproducibility of each analysis . Additionally, blank samples (50% MeOH) were also randomly included in the run to monitor the background noise and potential carry-over of analytes. Each sample was analysed in triplicate (analytical/technical replicates), and together with the three biological replicates, this generated n = 9, to have the minimum required number of replicates for metabolomic studies that involve multivariate analyses. A high definition SYNAPT G1 quadrupole time-of-flight (qTOF) mass spectrometry system (Waters Corporation, Manchester, UK) was coupled to the UHPLC chromatography system to detect metabolites and acquire data in both positive and negative electrospray ionisation (ESI) operation modes. The controlling software was MassLynx XS TM (Waters, Manchester, UK). A reference calibrant, leucine encephalin (554.2615 Da) was used as the ‘lockmass’ calibrant and allowed for typical mass accuracies between 1 to 3 mDa. The respective capillary and sampling cone voltages were set as 2.5 kV and 30 V. The desolvation temperature used was 450°C, with the source temperature set to 120°C, cone gas flow was set to 50 L/h, and the desolvation gas flow set to 550 L/h. An m/z range of 50–1200 was set with a scan time of 0.1 s. The desolvation-, collision- and cone gas used at a flow rate of 700 L/h was high-purity nitrogen. Data was acquired using five different collision energies (MS E ), ramping from 0–50 eV to cause fragmentation of the initial ions to ensure that information regarding the fragmentation of the respective compounds could be obtained for downstream structural elucidation and metabolite annotation . Data analysis The data sets obtained were explored and processed using the applications manager, MarkerLynx XS TM (Waters Corporation, Manchester, UK). The software makes use of a patented ApexTrack algorithm. The following parameters were used for processing: 2–23 min retention time (Rt) range and m/z range 50–1200 Da. The Rt window was set to 0.20 min and the mass window to 0.05 Da. The mass tolerance was 0.05 Da and the intensity threshold was set to 100 counts. The corrected data matrices were exported to ‘soft independent modelling of class analogy’ (SIMCA-version 16) software (Sartorius, Umeå, Sweden) for multivariate data analysis (MVDA). Unsupervised models, namely principal component analysis (PCA) and hierarchical clustering analysis (HCA) were used to reduce the dimensionality of the data sets and to explore the underlying structures and characteristics of the data. Supervised orthogonal projection to latent structures discriminant analysis (OPLS-DA) was used to compare the response of the oat plants to the different pathovars and identify discriminatory ions among the respective treatments. The OPLS-DA models were validated using rigorous validation methods that included cross-validated analysis of variance (CV-ANOVA) and receiver operator characteristic (ROC) analysis . Metabolite annotation and semi-quantitative comparison Metabolites were identified based on their respective accurate masses, fragmentation patterns, and possible elemental compositions generated from the MarkerLynx TM software. Each putatively suggested empirical formula was exported and searched for in various databases such as MetaCyc ( https://metacyc.org/ ), Plant metabolic network (PMN) ( https://plantcyc.org/ ), ChemSpider, Mass bank of North America ( https://mona.fiehnlab.ucdavis.edu/ ), Dictionary of Natural Products ( www.dnp.chemnetbase.com ), and the Kyoto Encyclopedia of Genes and Genomes (KEGG, www.genome.jp/kegg/ ). Metabolites were putatively identified to level 2 of the Metabolomics Standards Initiative (MSI)  unless specified otherwise. The generated matrixes of the annotated metabolites were exported from MarkerLynx XS™ as.csv files. MetaboAnalyst 5.0 ( https://www.metaboanalyst.ca/ ), an online platform for statistical, functional and integrative analysis of metabolomics data , was then utilised for visualisation of the MS files which contained the m/z , Rt and peak intensities of metabolites separated by chromatography and detected as ions by mass spectrometry. Data processing, integrity, missing values, filtering and normalisation were performed on MetaboAnalyst 5.0 followed by Pareto-scaling before statistical analyses to reduce variance within the features. A comparison of the presence and relative concentration/intensity of the identified metabolites among the various treatments were performed via dendrogram heatmap analyses using a Pearson distance measure and the Ward clustering algorithm in MetaboAnalyst 5.0 to visualise the distribution of metabolites. Additionally, the software’s metabolomics pathway analysis (MetPA) tool was employed to identify the crucial metabolic pathways induced in oat treated with the respective P . syringae pathovars. Each annotated metabolite’s KEGG (Kyoto Encyclopedia of Genes and Genomes; www.genome.jp/kegg/pathway.html ) identifier was employed as an input, and the KEGG metabolic pathways ( Arabidopsis thaliana ) served as the knowledge base for the construction of the pathways.   Seeds of the Dunnart oat cultivar was procured (Agricol, Pretoria, South Africa) and chosen for further infection studies using an initial screening procedure, as described below. Seedlings were cultivated in 10 cm, 250 mL pots (±15 seeds per pot) with a pasteurized (90˚C) germination mixture soil (Culterra, Muldersdrift, South Africa) and watered twice weekly with a water-soluble chemical fertiliser (Multisol ‘N’, Culterra, Muldersdrift, South Africa). The were cultivated in greenhouse conditions with a 12 h/12 h (light/dark) cycle, a light intensity of approximately 80 μmol/m 2 /s, and a temperature of 25°C. The study was designed to track the cultivar’s response to bacterial infection over time (2–6 days post inoculation—d.p.i). Triplicate pots were grown for each time point and were cultivated under identical conditions. When the seedlings achieved 3-week maturity (three-leaf seedling stage), they were thinned out and selected for uniformity in developmental stage before treatment with the bacterial suspensions described in the following sections. The experimental design included three independent biological replicates.   Pseudomonas syringae pathovars Dr. W. Kriel (Starke Ayres Seeds, Bredell, South Africa) provided an isolate pathogenic on oat namely P . syringae , pv. coronafaciens ( Ps-c ) with sequence information available on GenBank (PP940110), P . syringae pv. tabaci ( Ps-t ) was obtained from Prof. Y. Ichinose, Okayama University, Japan, and P . syringae pv. tomato DC3000 and its hrcC mutant from Dr. B. Kemmerling, University of Tuebingen, Germany. Inqaba Biotechnical Industries (Pretoria, South Africa) performed 16S rRNA sequencing to validate the pathovar identities. Ps isolates were then cultured and maintained on nutrient agar (Merck, Modderfontein, South Africa). A colony was selected under sterile circumstances in a laminar flow cabinet and cultivated overnight at 28°C in nutrient broth (Merck, Modderfontein, South Africa) on an orbital shaking incubator. The overnight culture’s OD 600 was measured and diluted with 0.1% Tween 20 and phosphate buffered saline (PBS) to a value of around 0.3. The same dilution was used for non-inoculated nutrient broth as a control (applied to vehicle control plants). The pathogenicity of P . syringae strains on several oat cultivars was examined by drop-inoculation  on leaf segments (OD600 ≈0.1, 0.2, 0.3). The inoculated leaf segments were stored in a high-humidity container for 5 days in a regulated growth environment with a 12-hour light-dark cycle at 25°C and visual symptoms scored daily. ‘Dunnart’ was chosen for further metabolomic analysis as a cultivar with a tolerant reaction based on initial visual observation tests that demonstrated its ability to tolerate infection during the 5-day period.   At the three-leaf growth stage (approximately 21 days post emergence), the plants were treated by spraying with the Ps-c , Ps-t , DC3000 and hrcC − bacterial suspensions (prepared in PBS with 0.1% Tween 20), diluted to OD 600 ≈0.3. The vehicle control (VC) plants were sprayed with a solution free of the bacteria and the healthy control (HC) groups were untreated (i.e., not sprayed with either solution), all grown under normal growth conditions. The 50 mL of either the inoculum (containing the bacteria) or the control (0.1% Tween 20 in PBS) solution was evenly sprayed onto the leaves of the treated and VC groups respectively. The plants were then incubated in darkness in an incubator for 1 h to provide 100% relative humidity. Following the 1 h incubation, the plants were removed and another 50 mL of either inoculum or control solution was applied to the treated and VC groups, respectively, and further incubated for 6 h. After incubation, the plants were then exposed again to the same initial conditions: with cycles of 12 h light/dark, a light intensity of 80 μmol/m 2 /s and temperature of 25°C. Post-treatment harvesting of plants was done for treated, VC and HC groups at 2, 4 and 6 d.p.i. by cutting the leaves and immediately snap freezing with liquid nitrogen to quench metabolic activity associated with possible wounding and handling of the tissue. Leaves were stored at −80°C until metabolite extraction.   ROS production was measured using a modified version of the luminometry method . Briefly, leaves were pressure infiltrated with the respective bacterial suspensions (OD 600 ≈0.3) or with H 2 O (control). The infiltrated sites were then excised and washed for 1 h in sterile water with agitation. The single leaf discs (4 mm in diameter) were then rinsed 3 times and individually placed into a 96-well microtitre plate with 20 μg/mL horseradish peroxidase (HRP) and 34 μg/mL of luminol (Sigma-Aldrich, St. Louis, MO, USA). The luminescence was then measured over a period of 30 min using a Synergy HT Biotek microplate reader (Biotek Instruments, Winooski, VT, USA).   POX activity associated with the innate immune response was measured using a modified version of the method described by . Leaf discs (4 mm in diameter) were punched from leaves of the oat plants and washed with agitation for 1 h in 1 mL of a 1 x MS (Murashige and Skoog) salt solution to remove any background activity caused by cutting the leaves. The discs were then carefully transferred to individual wells in a 96-well plate using a small spatula to minimise damage. Each well then received 50 μL of 1 x MS solution with either H 2 O (control) or one of the respective overnight bacterial suspensions (OD 600 ≈0.3). The plates were then sealed with parafilm and incubated for 20 h with agitation at room temperature. The leaf discs were removed and 50 μL of 5-aminosalicylic acid (Sigma-Aldrich, St. Louis, MO, USA) (1 mg/mL) at pH 6.0 with 0.01% hydrogen peroxide (Barrs Pharmaceuticals, Cape Town, RSA) was added to each well using a multichannel pipette in order to minimise timing differences. The reaction was allowed to proceed for 3 min and then stopped by adding 20 μM NaOH with a multichannel pipette in the same order as before, again to minimise timing differences and allow equal time for the enzymatic reactions to proceed. The plates were then analysed using a Synergy HT Biotek microplate reader (Biotek Instruments, Winooski, VT, USA) in absorbance mode at 600 nm (A 600 nm).   The harvested leaf material was quenched with liquid nitrogen before being crushed into powder with a pre-cooled mortar and pestle. The samples were weighed (1 g) into 50 mL Falcon tubes with addition of 80% cold (4°C) aqueous analytical grade methanol (Romil Chemistry, Cambridge, UK).at a 1:10 m/v ratio. The suspensions were then homogenised with a probe sonicator (Bandelin Sonopuls, Berlin, Germany) at 55% power for 10 sec per sample. To avoid cross-contamination, equipment was cleansed between each sample. The homogenates were centrifuged at 5100 x g for 20 min at 4°C in a benchtop centrifuge after which the supernatants were kept and concentrated by evaporating the methanol under vacuum to approximately 1 mL using a rotary evaporator set to 55°C. The concentrated samples were transferred to 2 mL microcentrifuge tubes and dried in a centrifugal evaporator under vacuum. The dried extracts were then reconstituted by dissolving in 500 μL of 50% aqueous methanol (MilliQ deionised water and LC-grade methanol (Romil, Cambridge, UK). The samples were subsequently filtered through nylon syringe filters (0.22 μm) into chromatography vials fitted with 500 μL inserts, capped, and kept at 4°C until analysis.   An Acquity UHPLC system (Waters Corporation, Manchester, UK) was used to analyse 2 μL of each sample, separated into its respective components using a binary solvent on an HSS T3 reverse-phase column (Waters Corporation, Billerica, MA, USA; 2.1 × 150 mm × 1.8 μm), thermostatted at 60°C. The solvents used were MilliQ water and acetonitrile (Romil Chemistry, Cambridge, UK), both containing 0.1% formic acid (Sigma-Aldrich, Munich, Germany) and 2.5% isopropanol (IPA, Romil, Cambridge, UK). The run was set to 30 min per injection with an elution gradient carried out via a binary solvent system consisting of 0.1% aqueous formic acid with 2.5% isopropanol (solvent A) and 0.1% formic acid and 2.5% isopropanol in acetonitrile (Romil, Cambridge, UK; solvent B) at a flow rate of 0.4 mL/min. The initial conditions were 95% A and 5% B and held for 1 min. A gradient was applied to change the chromatographic conditions to 10% A and 90% B at 25 min; and changed to 5% A and 95% B at 25.10 min. These conditions were held for 2 min and then changed to the initial conditions at 28 min. The analytical column was allowed to equilibrate for 2 min before each subsequent injection. Pooled quality control (QC) samples were also prepared to monitor the stability of the LC-MS system and assess the reliability and reproducibility of each analysis . Additionally, blank samples (50% MeOH) were also randomly included in the run to monitor the background noise and potential carry-over of analytes. Each sample was analysed in triplicate (analytical/technical replicates), and together with the three biological replicates, this generated n = 9, to have the minimum required number of replicates for metabolomic studies that involve multivariate analyses. A high definition SYNAPT G1 quadrupole time-of-flight (qTOF) mass spectrometry system (Waters Corporation, Manchester, UK) was coupled to the UHPLC chromatography system to detect metabolites and acquire data in both positive and negative electrospray ionisation (ESI) operation modes. The controlling software was MassLynx XS TM (Waters, Manchester, UK). A reference calibrant, leucine encephalin (554.2615 Da) was used as the ‘lockmass’ calibrant and allowed for typical mass accuracies between 1 to 3 mDa. The respective capillary and sampling cone voltages were set as 2.5 kV and 30 V. The desolvation temperature used was 450°C, with the source temperature set to 120°C, cone gas flow was set to 50 L/h, and the desolvation gas flow set to 550 L/h. An m/z range of 50–1200 was set with a scan time of 0.1 s. The desolvation-, collision- and cone gas used at a flow rate of 700 L/h was high-purity nitrogen. Data was acquired using five different collision energies (MS E ), ramping from 0–50 eV to cause fragmentation of the initial ions to ensure that information regarding the fragmentation of the respective compounds could be obtained for downstream structural elucidation and metabolite annotation .   The data sets obtained were explored and processed using the applications manager, MarkerLynx XS TM (Waters Corporation, Manchester, UK). The software makes use of a patented ApexTrack algorithm. The following parameters were used for processing: 2–23 min retention time (Rt) range and m/z range 50–1200 Da. The Rt window was set to 0.20 min and the mass window to 0.05 Da. The mass tolerance was 0.05 Da and the intensity threshold was set to 100 counts. The corrected data matrices were exported to ‘soft independent modelling of class analogy’ (SIMCA-version 16) software (Sartorius, Umeå, Sweden) for multivariate data analysis (MVDA). Unsupervised models, namely principal component analysis (PCA) and hierarchical clustering analysis (HCA) were used to reduce the dimensionality of the data sets and to explore the underlying structures and characteristics of the data. Supervised orthogonal projection to latent structures discriminant analysis (OPLS-DA) was used to compare the response of the oat plants to the different pathovars and identify discriminatory ions among the respective treatments. The OPLS-DA models were validated using rigorous validation methods that included cross-validated analysis of variance (CV-ANOVA) and receiver operator characteristic (ROC) analysis .   Metabolites were identified based on their respective accurate masses, fragmentation patterns, and possible elemental compositions generated from the MarkerLynx TM software. Each putatively suggested empirical formula was exported and searched for in various databases such as MetaCyc ( https://metacyc.org/ ), Plant metabolic network (PMN) ( https://plantcyc.org/ ), ChemSpider, Mass bank of North America ( https://mona.fiehnlab.ucdavis.edu/ ), Dictionary of Natural Products ( www.dnp.chemnetbase.com ), and the Kyoto Encyclopedia of Genes and Genomes (KEGG, www.genome.jp/kegg/ ). Metabolites were putatively identified to level 2 of the Metabolomics Standards Initiative (MSI)  unless specified otherwise. The generated matrixes of the annotated metabolites were exported from MarkerLynx XS™ as.csv files. MetaboAnalyst 5.0 ( https://www.metaboanalyst.ca/ ), an online platform for statistical, functional and integrative analysis of metabolomics data , was then utilised for visualisation of the MS files which contained the m/z , Rt and peak intensities of metabolites separated by chromatography and detected as ions by mass spectrometry. Data processing, integrity, missing values, filtering and normalisation were performed on MetaboAnalyst 5.0 followed by Pareto-scaling before statistical analyses to reduce variance within the features. A comparison of the presence and relative concentration/intensity of the identified metabolites among the various treatments were performed via dendrogram heatmap analyses using a Pearson distance measure and the Ward clustering algorithm in MetaboAnalyst 5.0 to visualise the distribution of metabolites. Additionally, the software’s metabolomics pathway analysis (MetPA) tool was employed to identify the crucial metabolic pathways induced in oat treated with the respective P . syringae pathovars. Each annotated metabolite’s KEGG (Kyoto Encyclopedia of Genes and Genomes; www.genome.jp/kegg/pathway.html ) identifier was employed as an input, and the KEGG metabolic pathways ( Arabidopsis thaliana ) served as the knowledge base for the construction of the pathways.   Evaluation of the oat leaf symptoms in response to the Pseudomonas syringae pathovars The disease severity of the oat cultivar (Dunnart) to the respective pathovars was assessed using visual observation, with scoring ending 6 d after first infection. Visual observation scoring was used with a scale of 0 to 5, with 0 indicating no symptoms, 1–3 indicating minor symptoms, 3–5 indicating moderate to severe disease symptoms . The typical symptoms associated with Ps-c infections include a HR that is characterised by necrosis at the inoculation site where a water-soaked lesion presented followed by the appearance of a characteristic yellow halo . Typically, Ps-t causes brown spots that spreads across the surface of the leaf and later develops yellow halos, however when constrained by veins, the spots turn angular. Tabtoxin is produced by both of these pathovars and is responsible for the characteristic yellow halo formation. DC3000 infection is associated with chlorosis of the leaves, as can be seen in , which is mainly due to the presence of the phytotoxin coronatine (COR). The hrcC mutant lacks the T3SS, therefore its growth in the plant is restricted. Some studies have found that tomato plants inoculated with hrcC − still developed some chlorotic patches that are suggestive of COR formation, however there was a lack of necrotic lesions that the wild type (DC3000) typically produces . Perception assays ROS luminescence assay The amount of light released owing to the oxidation of luminol by HRP is measured and represented in relative light units (RLU) to determine the ROS burst. The data provided is portrayed as the amount of RLU generated every minute over a 28 min period to track ROS production as an early defence response . Data is also represented as total RLU generated over the same period. The luminol based perception assays elucidates the capacity of the plants to detect and respond to infection with the respective Ps pathovars. The luminescence curves and total RLU indicate the production of ROS by the oat leaf discs in response to exposure with Ps-c , Ps-t , DC3000 and hrcC − compared to the control (dH 2 O). Twelve independent experiments were carried out with n = 36 samples (3 biological replicates for each treatment condition). The total RLU graphs for the respective treatments show that the Ps-c treatment produced the greatest response followed by Ps-t , DC3000 and the lowest production seen in the hrcC − treatment. The ROS burst for the Ps-c treatment lasted approximately 23 min before gradually decreasing back to basal levels. Treatment with Ps-t reached a maximum at 7 min and lasted for about 21 min. DC3000 treatment reached a maximum between 7–9 min and returned to basal levels at around 21 min post-treatment. The hrcC − treated group showed a shorter burst with lower RLU compared to the other treatments. POX activity A 96-well microtiter plate-based assay was used for analysing plant PTI by evaluating the activity of POX enzymes produced in response to bacterial treatment. The experiment was carried out using a fixed concentration of bacteria (OD 600 ≈0.3) to assess the kinetics of the plant POX response to the various strains of Ps . Unlike the ROS burst as detected by the luminol-based assay, which occurs over minutes, detectable POX activity is not found until 16 h after treatment, with increasing levels seen after 22 h . POX activity increased 9-, 9-, 7-, and 6-fold in the treated ( Ps-c , Ps-t , DC3000 and hrcC − ) plants, respectively compared to the control . The highest POX activity was observed for the Ps-c -treated plants. UHPLC-qTOF-MS Extracts from leaf tissues were analysed on an UHPLC-qTOF-MS system equipped with an ESI source for the identification of individual metabolites originating from metabolite classes of varied polarity. The data was collected in both ESI  modes, which was important for the investigation of particular metabolite classes with distinct chemical characteristics. After pathogen treatment, visual inspection of the base peak intensity (BPI) MS chromatograms  in negative ionisation mode revealed distinct variations in peak intensities as well as the presence/absence of peaks across all samples. The chromatograms demonstrate that bacterial infiltration with the different pathovars causes distinct metabolic changes. Initial optimisation tests revealed that the majority of extractable metabolites ionised better in the ESI  mode; thus, only these data sets are provided and illustrated further. This section describes an untargeted approach that was used to elucidate and identify as many statistically significant metabolites involved in oat plant defence as possible. Chemometric analyses for the elucidation of metabolic changes in response to the different Pseudomonas syringae pathovars The unique chromatographic profiles offered a visual evaluation confirming the occurrence of metabolic reprogramming resulting from pathogen treatment. This therefore prompted obtaining more information regarding these underlying changes by data mining and comparative chemometric analyses to reveal more significant structures within the datasets that distinguish the control from the treated conditions. The unsupervised PCA reduces the multi-dimensionality of the dataset and subsequently enables biological interpretability by projecting the data in a lower-dimensional plane, thereby exposing underlying structures, groupings and trends in the data sets. The principal components (PCs) of each model illustrates distinct treatment-related groupings when compared to the control for the Dunnart cultivar at 2-, 4- and 6 d.p.i.  respectively. Based on the respective PCA plots, HCA plots were constructed as dendrograms that outline the similarities and differences between the individual groups/clusters in a hierarchical format. The dendrograms at the respective time points illustrate that the data clusters into two main branches separating the control and treated groups . Further separation is apparent with subsequent branches forming within the treated groups at 2 d.p.i.  showing clusters of DC3000 and Ps-t vs Ps-c and hrcC − treated groups. At 4 d.p.i. overlapping clustering is observed  for the control and hrcC − treated groups suggesting similar metabolite content vs the other branch showing clustering of DC3000, Ps-t and Ps-c . A similar trend and grouping is observed at 6 d.p.i. . The control and treated groups were then further compared using a supervised discriminant analysis method with an example illustrated in where the OPLS-DA score plot demonstrates differential sample clustering between the control and hrcC − treated group. To choose discriminating ions between treatments the corresponding loadings S-plot was constructed with an example of the control vs hrcC − treatment shown in . The S-plot provides a visual explanation of the OPLS-DA model by displaying the properties ( m/z ions) that contribute to class differentiation. The loadings S-plot assisted in identifying metabolite signatures that appeared statistically significant for the respective treatments. The covariance (magnitude) and correlation (reliability) of the samples in the model are respectively represented on the axes as p and p(corr). MS spectral-based metabolite identification was used to identify discriminant ions with a |p(corr)| of ≥0.3, ≤0.3 and a co-variance value of |(p1)| ≥0.05, ≤0.05. The reliability of each model was assessed using cross-validation analysis of variance testing (CV-ANOVA) as a diagnostic tool, with models of significance having p -values of <0.05 . The permutation test ( n = 100) revealed that the OPLS-DA models had higher R 2 and Q 2 values than the 100 permuted models, implying that the produced models were statistically superior to the permuted models . The respective computed OPLS-DA models were used as binary classifiers and showed perfect discrimination, with the ROC curve passing through the top left corner indicating 100% sensitivity and specificity . OPLS-DA models with corresponding loadings S-plots were generated for all treatments along with ROC curves and permutation plots to test the validity of each model. Due to the large number of generated models only one example is shown, however the R 2 X(cum), R 2 Y, Q 2 (cum), and CV-ANOVA determined p -values for each generated model are captured in . A list of annotated (putatively identified) metabolites were generated after significant ions were selected from the respective loadings S-plots and presented in . The experimental section ‘Metabolite annotation and semi-quantitative comparison’ describes how the statistically significant variables were annotated. The potential chemical structures were investigated by examining the fragmentation patterns that were produced under various collision energies MS E . A total of 51 metabolites were identified for the different treatments and time points. In the table, metabolites are shown that presented as discriminatory for a particular treatment (with a variable importance in projection score, VIP > 1.0) as well as metabolites that were readily identified and had a fold change ≥ 1.5 (indicated with *). The annotated metabolites were also categorised based on the respective metabolite classes which included: fatty acids, amino acids, phenolic acids, phenolic acid amides, steroidal saponins, flavonoids, and alkaloids. Time-related profiling of oat plants upon treatment with pathovars of Pseudomonas syringae In order to gain more insight into the metabolic reprogramming that occurs in the leaf tissue of the oat seedlings in response to treatment with the P . syringae pathovars, biochemical interpretations were drawn from the putatively annotated metabolites . In addition to the amino acids and fatty acids which are classified as primary metabolites, phenolic acids and phenolic amides, saponins, flavonoids, and alkaloids were among the metabolite groups identified as secondary or specialised metabolites. Changes in these metabolite classes were analysed from a metabolome perspective and revealed alterations involved due to the bacterial treatments. The presence and amount of the respective metabolites in the various treatment and control groups were considered using data visualisation tools available on MetaboAnalyst 5.0. A heatmap was constructed using the average integrated peak areas of the different metabolites . The infographic clearly distinguishes between the treated and control groups. The distinguishing characteristics are evidently illustrated among the respective treatments. The profile could be beneficial in identifying metabolic markers linked with disease resistance in oat plants against the respective pathovars of P . syringae . By profiling these host and nonhost responses deeper insight can be gained into oat plant defence mechanisms. Venn diagrams were constructed to highlight overlapping and unique metabolites. represents Venn diagrams constructed at the respective time points 2- , 4–  and 6 d.p.i. . Unique metabolites presented for the respective groups and time-specific metabolites were noted for certain treatments. Starting with 2 d.p.i., 7 metabolites were found as unique for the control group (feruloylserotonin, syringin, 3-O-feruloylquinic acid, sinapaldehyde glucoside, isovitexin-7-O-glucosde, coumaroylquinic acid and dirhamnosyl linolenic acid). At 4- and 6 d.p.i. only 3 (namely isoamoritin, a triprenylated flavanone, sinapaldehyde glucoside and syringin) and 6 (isoamoritin, sinapaldehyde glucoside, coumaroylquinic acid, syringin, 3-O-feruloylquinic acid and 9-oxo-12,13-dihydroxy-10,15-octadecadienoic acid) metabolites presented as unique to the control at these time points respectively. The only common metabolites that presented across the three time points were syringin and sinapaldehyde glucoside. The Ps-c -treated group revealed hordenine, protocatechuic acid hexose, dihydroferulic acid 4-O-glucuronide and gentisic acid glucoside as signatory metabolites at 2 d.p.i. At 4- and 6 d.p.i. only avenanthramide (Avn) A presented as unique for this treatment. Interestingly, for the Ps-t treatment each time point presented different metabolites, at 2 d.p.i. (palmitoleic-linoleic glucoside and tricin ether glucopyranoside), 4 d.p.i. (feruloylserotonin) and 6 d.p.i. (hordenine). The DC3000 treatment showed linoleic acid and ononin (a isoflavone glycoside) as signatory features at 2 d.p.i. and ethyl 7-epi-12-hydroxyjasmonate glucoside at 4- and 6 d.p.i.. Treatment with the hrcC mutant revealed the prenylated isoflavonoid kanzonol I as signatory at 2 d.p.i., whereas 3 metabolites presented at 4 d.p.i. (26-desglucoavenacoside A, sinapic acid glucose and isovolubilin, an isoflavone rhamnoside) and 6 d.p.i. (sinapic acid glucose and isovolubilin). The only shared metabolites for this treatment being sinapic acid glucose and isovolubilin at 4- and 6 d.p.i.. Among the control and all the treated groups, only one common overlapping metabolite was noted at 2 d.p.i. (avenacoside A). Among the four treated groups, 3 metabolites overlapped for 2- and 4 d.p.i. namely 26-desglucoavenacoside A, naringenin, hydroxylinolenic acid for 2 d.p.i. and Avn E, quercetin dimethyl ether methylbutyrate and naringenin for 4 d.p.i.. At 6 d.p.i., naringenin and tryptophan were present as unique for the four treated groups with naringenin (as the precursor of flavonoid metabolites) being the common overlapping metabolite for all the treatments across all time points. Other inferences that can be drawn from the Venn diagrams without mentioning all the specific overlapping metabolites, is that among the respective groups the greatest overlapping responses are seen among the Ps-c and hrcC − treated groups. Some of the key metabolites that overlap for these two treatments are Avns A, B and O/R at 2 d.p.i. and Avn O/R being commonly present across all time points. These defence-related metabolites interact with one another through various metabolic pathways rather than acting independently. Metabolic pathway mapping was thus carried out to elucidate the most pertinent pathways implicated in oat responses to inoculation with the respective pathovars of P . syringae , which will aid in the biochemical interpretation of the post-infection metabolic perturbations. MetaboAnalyst 5.0 was used to carry out the Metabolomics Pathway Analysis (MetPA). This extremely sensitive web-based application is helpful when it comes to analysing and visualising metabolomic data and can pick up small variations between various compounds. Relative concentration variations of compounds with known KEGG or HMDB identifiers can be used to construct these biological pathways . The computed metabolic pathways are presented according to significance or pathway impact as shown in . The most significant pathways (displayed on the y-axis) were the phenylpropanoid—and phenylalanine, tyrosine and tryptophan metabolism pathways, whereas the most impactful (displayed on the x-axis) were the linoleic acid pathway and (in general) pathways involved in the biosynthesis of secondary metabolites. Some of the most frequently occurring secondary metabolites that are involved in plant growth and defence against abiotic and biotic stressors includes phenolics, flavonoids, coumarins, and lignins that are produced via the phenylpropanoid pathway . The phenylpropanoid—and general secondary metabolite biosynthesis pathways show overlap where the conversion of phenylalanine to p -coumaroyl-CoA can be seen. Coumaroyl-CoA (KEGG ID = C00223) is an important intermediate in the synthesis of a multitude of secondary metabolites such as the respective flavonoids and the phenolic amides. The presence and distribution of these substances, from both pathways, at the cellular, tissue, and organ levels across the plant kingdom emphasises the myriads of biological and biochemical processes that are vital to the survival of plants . It is well known that these phenolic acids and derivatives such as phenolic amide conjugates can act as naturally occurring antibiotic molecules (pre-existing phytoanticipins or inducible phytoalexins) in plant-pathogen interactions . Unsaturated fatty acids (linoleic acids, C18:2) are abundant in plant membranes, which makes them crucial for plant cell structure and function and supportive of adaptive responses through membrane remodelling . In addition, they play a role in the formation oxylipins and signalling molecules such as 12-oxo-phytodienoic acid (OPDA). OPDA in turn acts as precursor of jasmonic acid (JA) and derivatives, which are produced in response to tissue damage brought on by insects, pathogens, herbivores, or mechanical stress  .   Pseudomonas syringae pathovars The disease severity of the oat cultivar (Dunnart) to the respective pathovars was assessed using visual observation, with scoring ending 6 d after first infection. Visual observation scoring was used with a scale of 0 to 5, with 0 indicating no symptoms, 1–3 indicating minor symptoms, 3–5 indicating moderate to severe disease symptoms . The typical symptoms associated with Ps-c infections include a HR that is characterised by necrosis at the inoculation site where a water-soaked lesion presented followed by the appearance of a characteristic yellow halo . Typically, Ps-t causes brown spots that spreads across the surface of the leaf and later develops yellow halos, however when constrained by veins, the spots turn angular. Tabtoxin is produced by both of these pathovars and is responsible for the characteristic yellow halo formation. DC3000 infection is associated with chlorosis of the leaves, as can be seen in , which is mainly due to the presence of the phytotoxin coronatine (COR). The hrcC mutant lacks the T3SS, therefore its growth in the plant is restricted. Some studies have found that tomato plants inoculated with hrcC − still developed some chlorotic patches that are suggestive of COR formation, however there was a lack of necrotic lesions that the wild type (DC3000) typically produces .   ROS luminescence assay The amount of light released owing to the oxidation of luminol by HRP is measured and represented in relative light units (RLU) to determine the ROS burst. The data provided is portrayed as the amount of RLU generated every minute over a 28 min period to track ROS production as an early defence response . Data is also represented as total RLU generated over the same period. The luminol based perception assays elucidates the capacity of the plants to detect and respond to infection with the respective Ps pathovars. The luminescence curves and total RLU indicate the production of ROS by the oat leaf discs in response to exposure with Ps-c , Ps-t , DC3000 and hrcC − compared to the control (dH 2 O). Twelve independent experiments were carried out with n = 36 samples (3 biological replicates for each treatment condition). The total RLU graphs for the respective treatments show that the Ps-c treatment produced the greatest response followed by Ps-t , DC3000 and the lowest production seen in the hrcC − treatment. The ROS burst for the Ps-c treatment lasted approximately 23 min before gradually decreasing back to basal levels. Treatment with Ps-t reached a maximum at 7 min and lasted for about 21 min. DC3000 treatment reached a maximum between 7–9 min and returned to basal levels at around 21 min post-treatment. The hrcC − treated group showed a shorter burst with lower RLU compared to the other treatments. POX activity A 96-well microtiter plate-based assay was used for analysing plant PTI by evaluating the activity of POX enzymes produced in response to bacterial treatment. The experiment was carried out using a fixed concentration of bacteria (OD 600 ≈0.3) to assess the kinetics of the plant POX response to the various strains of Ps . Unlike the ROS burst as detected by the luminol-based assay, which occurs over minutes, detectable POX activity is not found until 16 h after treatment, with increasing levels seen after 22 h . POX activity increased 9-, 9-, 7-, and 6-fold in the treated ( Ps-c , Ps-t , DC3000 and hrcC − ) plants, respectively compared to the control . The highest POX activity was observed for the Ps-c -treated plants.   The amount of light released owing to the oxidation of luminol by HRP is measured and represented in relative light units (RLU) to determine the ROS burst. The data provided is portrayed as the amount of RLU generated every minute over a 28 min period to track ROS production as an early defence response . Data is also represented as total RLU generated over the same period. The luminol based perception assays elucidates the capacity of the plants to detect and respond to infection with the respective Ps pathovars. The luminescence curves and total RLU indicate the production of ROS by the oat leaf discs in response to exposure with Ps-c , Ps-t , DC3000 and hrcC − compared to the control (dH 2 O). Twelve independent experiments were carried out with n = 36 samples (3 biological replicates for each treatment condition). The total RLU graphs for the respective treatments show that the Ps-c treatment produced the greatest response followed by Ps-t , DC3000 and the lowest production seen in the hrcC − treatment. The ROS burst for the Ps-c treatment lasted approximately 23 min before gradually decreasing back to basal levels. Treatment with Ps-t reached a maximum at 7 min and lasted for about 21 min. DC3000 treatment reached a maximum between 7–9 min and returned to basal levels at around 21 min post-treatment. The hrcC − treated group showed a shorter burst with lower RLU compared to the other treatments.   A 96-well microtiter plate-based assay was used for analysing plant PTI by evaluating the activity of POX enzymes produced in response to bacterial treatment. The experiment was carried out using a fixed concentration of bacteria (OD 600 ≈0.3) to assess the kinetics of the plant POX response to the various strains of Ps . Unlike the ROS burst as detected by the luminol-based assay, which occurs over minutes, detectable POX activity is not found until 16 h after treatment, with increasing levels seen after 22 h . POX activity increased 9-, 9-, 7-, and 6-fold in the treated ( Ps-c , Ps-t , DC3000 and hrcC − ) plants, respectively compared to the control . The highest POX activity was observed for the Ps-c -treated plants.   Extracts from leaf tissues were analysed on an UHPLC-qTOF-MS system equipped with an ESI source for the identification of individual metabolites originating from metabolite classes of varied polarity. The data was collected in both ESI  modes, which was important for the investigation of particular metabolite classes with distinct chemical characteristics. After pathogen treatment, visual inspection of the base peak intensity (BPI) MS chromatograms  in negative ionisation mode revealed distinct variations in peak intensities as well as the presence/absence of peaks across all samples. The chromatograms demonstrate that bacterial infiltration with the different pathovars causes distinct metabolic changes. Initial optimisation tests revealed that the majority of extractable metabolites ionised better in the ESI  mode; thus, only these data sets are provided and illustrated further. This section describes an untargeted approach that was used to elucidate and identify as many statistically significant metabolites involved in oat plant defence as possible.   Pseudomonas syringae pathovars The unique chromatographic profiles offered a visual evaluation confirming the occurrence of metabolic reprogramming resulting from pathogen treatment. This therefore prompted obtaining more information regarding these underlying changes by data mining and comparative chemometric analyses to reveal more significant structures within the datasets that distinguish the control from the treated conditions. The unsupervised PCA reduces the multi-dimensionality of the dataset and subsequently enables biological interpretability by projecting the data in a lower-dimensional plane, thereby exposing underlying structures, groupings and trends in the data sets. The principal components (PCs) of each model illustrates distinct treatment-related groupings when compared to the control for the Dunnart cultivar at 2-, 4- and 6 d.p.i.  respectively. Based on the respective PCA plots, HCA plots were constructed as dendrograms that outline the similarities and differences between the individual groups/clusters in a hierarchical format. The dendrograms at the respective time points illustrate that the data clusters into two main branches separating the control and treated groups . Further separation is apparent with subsequent branches forming within the treated groups at 2 d.p.i.  showing clusters of DC3000 and Ps-t vs Ps-c and hrcC − treated groups. At 4 d.p.i. overlapping clustering is observed  for the control and hrcC − treated groups suggesting similar metabolite content vs the other branch showing clustering of DC3000, Ps-t and Ps-c . A similar trend and grouping is observed at 6 d.p.i. . The control and treated groups were then further compared using a supervised discriminant analysis method with an example illustrated in where the OPLS-DA score plot demonstrates differential sample clustering between the control and hrcC − treated group. To choose discriminating ions between treatments the corresponding loadings S-plot was constructed with an example of the control vs hrcC − treatment shown in . The S-plot provides a visual explanation of the OPLS-DA model by displaying the properties ( m/z ions) that contribute to class differentiation. The loadings S-plot assisted in identifying metabolite signatures that appeared statistically significant for the respective treatments. The covariance (magnitude) and correlation (reliability) of the samples in the model are respectively represented on the axes as p and p(corr). MS spectral-based metabolite identification was used to identify discriminant ions with a |p(corr)| of ≥0.3, ≤0.3 and a co-variance value of |(p1)| ≥0.05, ≤0.05. The reliability of each model was assessed using cross-validation analysis of variance testing (CV-ANOVA) as a diagnostic tool, with models of significance having p -values of <0.05 . The permutation test ( n = 100) revealed that the OPLS-DA models had higher R 2 and Q 2 values than the 100 permuted models, implying that the produced models were statistically superior to the permuted models . The respective computed OPLS-DA models were used as binary classifiers and showed perfect discrimination, with the ROC curve passing through the top left corner indicating 100% sensitivity and specificity . OPLS-DA models with corresponding loadings S-plots were generated for all treatments along with ROC curves and permutation plots to test the validity of each model. Due to the large number of generated models only one example is shown, however the R 2 X(cum), R 2 Y, Q 2 (cum), and CV-ANOVA determined p -values for each generated model are captured in . A list of annotated (putatively identified) metabolites were generated after significant ions were selected from the respective loadings S-plots and presented in . The experimental section ‘Metabolite annotation and semi-quantitative comparison’ describes how the statistically significant variables were annotated. The potential chemical structures were investigated by examining the fragmentation patterns that were produced under various collision energies MS E . A total of 51 metabolites were identified for the different treatments and time points. In the table, metabolites are shown that presented as discriminatory for a particular treatment (with a variable importance in projection score, VIP > 1.0) as well as metabolites that were readily identified and had a fold change ≥ 1.5 (indicated with *). The annotated metabolites were also categorised based on the respective metabolite classes which included: fatty acids, amino acids, phenolic acids, phenolic acid amides, steroidal saponins, flavonoids, and alkaloids.   Pseudomonas syringae In order to gain more insight into the metabolic reprogramming that occurs in the leaf tissue of the oat seedlings in response to treatment with the P . syringae pathovars, biochemical interpretations were drawn from the putatively annotated metabolites . In addition to the amino acids and fatty acids which are classified as primary metabolites, phenolic acids and phenolic amides, saponins, flavonoids, and alkaloids were among the metabolite groups identified as secondary or specialised metabolites. Changes in these metabolite classes were analysed from a metabolome perspective and revealed alterations involved due to the bacterial treatments. The presence and amount of the respective metabolites in the various treatment and control groups were considered using data visualisation tools available on MetaboAnalyst 5.0. A heatmap was constructed using the average integrated peak areas of the different metabolites . The infographic clearly distinguishes between the treated and control groups. The distinguishing characteristics are evidently illustrated among the respective treatments. The profile could be beneficial in identifying metabolic markers linked with disease resistance in oat plants against the respective pathovars of P . syringae . By profiling these host and nonhost responses deeper insight can be gained into oat plant defence mechanisms. Venn diagrams were constructed to highlight overlapping and unique metabolites. represents Venn diagrams constructed at the respective time points 2- , 4–  and 6 d.p.i. . Unique metabolites presented for the respective groups and time-specific metabolites were noted for certain treatments. Starting with 2 d.p.i., 7 metabolites were found as unique for the control group (feruloylserotonin, syringin, 3-O-feruloylquinic acid, sinapaldehyde glucoside, isovitexin-7-O-glucosde, coumaroylquinic acid and dirhamnosyl linolenic acid). At 4- and 6 d.p.i. only 3 (namely isoamoritin, a triprenylated flavanone, sinapaldehyde glucoside and syringin) and 6 (isoamoritin, sinapaldehyde glucoside, coumaroylquinic acid, syringin, 3-O-feruloylquinic acid and 9-oxo-12,13-dihydroxy-10,15-octadecadienoic acid) metabolites presented as unique to the control at these time points respectively. The only common metabolites that presented across the three time points were syringin and sinapaldehyde glucoside. The Ps-c -treated group revealed hordenine, protocatechuic acid hexose, dihydroferulic acid 4-O-glucuronide and gentisic acid glucoside as signatory metabolites at 2 d.p.i. At 4- and 6 d.p.i. only avenanthramide (Avn) A presented as unique for this treatment. Interestingly, for the Ps-t treatment each time point presented different metabolites, at 2 d.p.i. (palmitoleic-linoleic glucoside and tricin ether glucopyranoside), 4 d.p.i. (feruloylserotonin) and 6 d.p.i. (hordenine). The DC3000 treatment showed linoleic acid and ononin (a isoflavone glycoside) as signatory features at 2 d.p.i. and ethyl 7-epi-12-hydroxyjasmonate glucoside at 4- and 6 d.p.i.. Treatment with the hrcC mutant revealed the prenylated isoflavonoid kanzonol I as signatory at 2 d.p.i., whereas 3 metabolites presented at 4 d.p.i. (26-desglucoavenacoside A, sinapic acid glucose and isovolubilin, an isoflavone rhamnoside) and 6 d.p.i. (sinapic acid glucose and isovolubilin). The only shared metabolites for this treatment being sinapic acid glucose and isovolubilin at 4- and 6 d.p.i.. Among the control and all the treated groups, only one common overlapping metabolite was noted at 2 d.p.i. (avenacoside A). Among the four treated groups, 3 metabolites overlapped for 2- and 4 d.p.i. namely 26-desglucoavenacoside A, naringenin, hydroxylinolenic acid for 2 d.p.i. and Avn E, quercetin dimethyl ether methylbutyrate and naringenin for 4 d.p.i.. At 6 d.p.i., naringenin and tryptophan were present as unique for the four treated groups with naringenin (as the precursor of flavonoid metabolites) being the common overlapping metabolite for all the treatments across all time points. Other inferences that can be drawn from the Venn diagrams without mentioning all the specific overlapping metabolites, is that among the respective groups the greatest overlapping responses are seen among the Ps-c and hrcC − treated groups. Some of the key metabolites that overlap for these two treatments are Avns A, B and O/R at 2 d.p.i. and Avn O/R being commonly present across all time points. These defence-related metabolites interact with one another through various metabolic pathways rather than acting independently. Metabolic pathway mapping was thus carried out to elucidate the most pertinent pathways implicated in oat responses to inoculation with the respective pathovars of P . syringae , which will aid in the biochemical interpretation of the post-infection metabolic perturbations. MetaboAnalyst 5.0 was used to carry out the Metabolomics Pathway Analysis (MetPA). This extremely sensitive web-based application is helpful when it comes to analysing and visualising metabolomic data and can pick up small variations between various compounds. Relative concentration variations of compounds with known KEGG or HMDB identifiers can be used to construct these biological pathways . The computed metabolic pathways are presented according to significance or pathway impact as shown in . The most significant pathways (displayed on the y-axis) were the phenylpropanoid—and phenylalanine, tyrosine and tryptophan metabolism pathways, whereas the most impactful (displayed on the x-axis) were the linoleic acid pathway and (in general) pathways involved in the biosynthesis of secondary metabolites. Some of the most frequently occurring secondary metabolites that are involved in plant growth and defence against abiotic and biotic stressors includes phenolics, flavonoids, coumarins, and lignins that are produced via the phenylpropanoid pathway . The phenylpropanoid—and general secondary metabolite biosynthesis pathways show overlap where the conversion of phenylalanine to p -coumaroyl-CoA can be seen. Coumaroyl-CoA (KEGG ID = C00223) is an important intermediate in the synthesis of a multitude of secondary metabolites such as the respective flavonoids and the phenolic amides. The presence and distribution of these substances, from both pathways, at the cellular, tissue, and organ levels across the plant kingdom emphasises the myriads of biological and biochemical processes that are vital to the survival of plants . It is well known that these phenolic acids and derivatives such as phenolic amide conjugates can act as naturally occurring antibiotic molecules (pre-existing phytoanticipins or inducible phytoalexins) in plant-pathogen interactions . Unsaturated fatty acids (linoleic acids, C18:2) are abundant in plant membranes, which makes them crucial for plant cell structure and function and supportive of adaptive responses through membrane remodelling . In addition, they play a role in the formation oxylipins and signalling molecules such as 12-oxo-phytodienoic acid (OPDA). OPDA in turn acts as precursor of jasmonic acid (JA) and derivatives, which are produced in response to tissue damage brought on by insects, pathogens, herbivores, or mechanical stress  .   Surface and intracellular recognition are all components of the host’s immune surveillance with significant convergence and reciprocal interplay between PTI and ETI. Effector-triggered susceptibility (ETS) is the outcome of pathogens being able to inject suppressors into plant cells to overcome activated plant defences . In this molecular arms race, both the pathogen and the host go through natural selection to diversify their effectors and R proteins . ETI is the cornerstone of traditional gene-for-gene resistance and is initiated by the recognition of effector or Avr proteins by corresponding plant R proteins. Upon recognition the majority of R proteins participate in the initiation of the HR by producing ROS . As such, the boundary between ‘susceptibility’ or ‘resistance’ (i.e., non-effective or effective nonhost resistance) is frequently an equilibrium between these two states that might be perturbed by the gain or loss of individual genes in the host or the microbe . In general, the phenotypic outcome would be classified as resistant, tolerant, or susceptible based on the sum of PTI and ETI, minus inhibitory effects from ETS (also considering the timing and intensity of responses). The development and severity of symptoms in oat crops vary depending on the plant-pathogen interaction and environmental factors . Since the oat plants were infected in a controlled environment, the symptomatic differences reflect responses to the infection with the respective pathovars. It is therefore critical to determine the precise variety or cultivar, as significant differences can occur within cultivars of the same species in terms of susceptibility or resistance/tolerance to various diseases . To resist changing environmental and pathogenic threats, plants rely on an inherent sophisticated and multi-layered innate immune system . As a result, identifying the metabolic phenotypes associated with oat defence responses to infection with the respective Ps pathovars would reveal more information about the cellular and metabolic pathways involved in the plant-pathogen relationship . A variety of studies have reported important factors involved in plant nonhost responses such as ROS , one of the earliest observable defence strategies in plants. One of the first studies on the role of ROS was done on type II nonhost responses  and showed an accumulation of H 2 O 2 as a strong form of resistance against bacteria. ROS are known to have two essential roles in plant defence to infection, one being the accumulation leading to a HR at the infected site and inhibiting further pathogen growth. The second role is the activation of resistance genes . In a study by Smith and Heese  the early production of ROS in Arabidopsis plants treated with the DC3000 and hrcC mutant strains of P . syringae pv. tomato was evaluated. The ROS production showed similar responses in terms of timing, intensity and duration. Similarly, the timing, intensity and duration of ROS production for both these treatments in this study also showed no significant differences. The triggered ROS production clearly functions independent of the T3SS, as evidenced by the lack of statistically significant differences between the two strains. The ROS response observed in this investigation is thus caused by PTI-dependent processes, which are consistent with effectors being transported into host cells at a much later time post infection. Peroxidases are another important component of plant defence and are activated in host plant tissues by pathogen infection to prevent cellular spread of infection by creating structural barriers (e.g., lignification using mono-lignols as substrates) or by releasing large amounts of ROS to create toxic conditions . The POX assay used measures activity of POX enzymes produced in response to activation of plant innate immunity and is purportedly a marker of PTI . The results in this study showed increases in POX activity in the treated plants as compared to controls. Such increases in POX activity signify increases in levels of lignin formation, suberisation and the HR . Due to POXs being able to operate as both ROS producers and catalytic enzymes depending on a variety of factors such as the availability of substrates and a range of other reaction conditions, it becomes difficult to identify the enzyme responsible for the observed POX activity. However, due to the initial washing phase of the assay, it is unlikely that the POX activity is associated with the ROS burst since any activity would be lost during this step, therefore it is more likely to be involved in cell-wall reinforcement. The result/outcome of the host response to an attempted infection may be affected by qualitative and quantitative variations, in particular metabolites or classes of metabolites within wider metabolomic profiles. Since the metabolites were seen to accumulate in the leaves in various amounts and with distinct accumulation patterns, it implies that differential reprogramming has occurred over time due to the respective treatments. High or low accumulation at particular time points signify early-, late-, or oscillatory reactions. According to the time-dependent reprogramming, infected plants modify their metabolomes toward inducible defence responses to restrict pathogen entry and further multiplication. The extensively researched host resistance is frequently referred to as gene-for-gene resistance and is typically extremely specific to a particular genotype or cultivar of plant and against a certain pathogen. Nonhost resistance, on the other hand, is a broad-spectrum resistance displayed by the entire plant species against a given disease and is not pathogen specific. Nonhost resistance is multi-tiered with many barriers reliant on a specific host to prevent pathogen colonisation . These barriers range from preformed barriers and phytoanticipins, to induced defence responses like the HR, lignin accumulation, production of antimicrobials (phytoalexins), and induction of pathogenesis-related (PR) proteins . Relevant to the evaluation of the results of this study, nonhost resistance has been described ‘as a gradual transitional phenomenon that is influenced by various exogenous factors and that there is ‘a continuum between nonhost and host plants, with several possible intermediate forms’ (e.g., ‘marginal/near/intermediate host’ and ‘near/intermediate nonhost’ or ‘apparent nonhost’) . In this study, an untargeted metabolomics approach has been applied to highlight the metabolic reprogramming that occurs in the oat seedlings (Dunnart cultivar, tolerant to Ps-c ) under treatment with different pathovars of P . syringae . Among these, Ps-c is known to cause a host response while Ps-t , DC3000 and the hrcC mutant elicit a nonhost response ( hrcC encodes a putative outer membrane protein that is conserved in all T3SS). Among the nonhost responses, Ps-t can further be described as showing a type II response since symptoms indicative of an HR appeared on leaves treated with this particular pathovar. DC3000 and its hrcC mutant on the other hand showed a response that can be categorised as a type I nonhost response. A mRNA profiling study of compatible and incompatible Interactions of Arabidopsis with P . syringae strains demonstrated that a big portion of the difference between incompatible and compatible interactions can be explained quantitatively with a saturating response curve model, where the plant response in an incompatible interaction was strong but that of a compatible interaction was not . Accordingly, by comparing the responses triggered by DC3000 and the hrcC mutant, a metabolic profile can be elucidated that reflects the outcome of effector proteins on the plant defence response (i.e., PTI vs PTI and ETI). The metabolites identified (phytoanticipins and phytoalexins) were categorised into classes of phenolic acids and—amides, saponins, flavonoids, and alkaloids. The summary figures below show the distribution of these metabolite classes among the different treatments . The study revealed phenolic amides as discriminatory metabolites, unique in the treated groups compared to the controls. These metabolites, also sometimes referred to as hydroxycinnamic acid amides (HCAAs), are secondary metabolites that are widely distributed throughout the plant kingdom. The structure of these metabolites is often characterised by the association of at least one hydroxycinnamic acid derivative linked via an amide bond to an aromatic mono- or an aliphatic polyamine . The most commonly found hydroxycinnamic acids in nature include p -coumaric -, caffeic -, ferulic—and sinapic acids along with some glycosylated forms. In oat crops, these phenolic amides consist of an anthranilic acid bound to a hydroxycinnamic acid and named as avenanthramides (Avns) when dimerised . The two precursor metabolites involved in the synthesis of Avns are p -coumaric acid that forms p -coumaroyl-CoA (early phenylpropanoid pathway) and caffeoyl-CoA that is converted to feruloyl-CoA in the general secondary metabolite biosynthesis pathway . The enzymatic synthesis of p -coumaric acid by the stress-responsive phenylalanine ammonia lyase (PAL) from phenylalanine and cinnamate-4-hydroxylase (C4′H) or directly from tyrosine by tyrosine ammonia lyase (TAL), initiates the Avn biosynthesis. The enzyme hydroxyanthranilate N-hydroxycinnamoyl transferase (HHT) catalyses the condensation of 5-hydroxyanthranilic acid with p -coumaric acid to form Avn A after (4CL) converts p -coumaric acid into its activated CoA thioester. On the other hand, the p -coumaroyl-CoA is frequently transformed to p -coumaroyl shikimate or quinate before being hydroxylated to generate caffeoyl-CoA by the p -coumaroyl-CoA ester 3′-hydroxylase. In the presence of HHT, the caffeoyl-CoA is then condensed with 5-hydroxyanthranilic acid to create Avn C. The caffeoyl-CoA O-methyltransferase (CCoAOMT) enzyme then methylates Avn C to produce Avn B . To date, over 40 different types have been reported in leaves and grains of oat plants based on their structure . The most abundant Avns are A, B and C. These metabolites were identified as discriminatory in the Ps-c , hrcC − and DC3000 treated groups . Other less abundant Avns were also identified as discriminatory and included Avns O/R, L and E ( Ps-c and hrcC − ). Avns are produced in response to pathogen infection or when oat leaves are treated with various elicitors. In oat, these Avns function as phytoalexins that are induced to act as chemical defence barriers and substrates for cell wall reinforcement upon exposure to pathogens. Additional to the antimicrobial activity, Avns are also potent antioxidant and radical scavenging compounds due to the cinnamic acid structure and the hydroxyanthranilic acid moiety . A greater variety of Avns were detected as discriminatory in the host response to Ps-c compared to the nonhost response to Ps-t . When infected with the hrcC mutant and the DC3000 strain (type I nonhost response), the cultivar presented a multitude of Avns. Since these metabolites are produced in all treatments (albeit with differences in apparent quantitative levels), it is suggested that the initial PAMP recognition is enough to trigger the production of these metabolites as the majority of the identified Avns were present as discriminatory ions for the hrcC mutant treated group. Avenacoside A and B and their biologically active counterparts (26-desglucoavenacoside A and B) were detected as discriminatory between the treated and control groups. Avenacosides are saponin phytoanticipins that are physiologically inactive until they are converted by avenacosidase into biologically active 26-desglucoavenacosides in response to tissue damage or pathogen infection . Here, avenacoside A and B can be seen as showing a decrease in level from the control to the infected groups, with 26-desglucoavenacoside A and B increasing in the treated groups. Since the conversion to the physiologically active form was noted in all the treatments, it is evidence that the respective pathovars elicited a clear defence response. These compounds are important in plant defence and have been greatly studied due to their antifungal properties where they are able to bind with sterols in the pathogen membrane and disrupt its integrity, which is the primary mechanism of action against the pathogen. As a result of saponin aggregation with sterol groups, this mechanism is thought to lead to the formation of transmembrane pores, which ultimately leads to cell death . It was proposed that similar saponins may also interfere with the bacterial outer membrane’s permeability . Lipopolysaccharides covers over 90% of the surface of Gram-negative bacteria’s cell walls. The study therefore proposed that the lipid A component of lipopolysaccharides may interact with the saponin thereby increasing the permeability of the bacterial cell wall as a result. Phenolic acids are produced and accumulate in plant tissues in response to stress and/or pathogen attack, where they act as protective agents against invading pathogens. These phenolics are often present as simple glycosides, conjugates to esters or amides or linked to the cell wall; and exhibit insecticidal, antimicrobial, antioxidant, and free radical scavenger activities . Depending on their chemical composition, they can be classified as simple phenolics, tannins, coumarins, flavonoids, chromones, xanthones, stilbenes, or lignans. These chemicals are mostly synthesised from shikimic acid as a precursor compound by the phenylpropanoid- and/or malonate pathway. The shikimic acid pathway is used to produce phenylalanine, leading to (hydroxy)cinnamic acids, and their derivatives, including coumarins, lignans, and simple phenols through deamination, hydroxylation, and methylation . Here a greater number of phenolic acids were present in the Ps-c and Ps-t -treated groups. More flavonoids were however present in the hrcC − and DC3000-treated groups. The role of flavonoids in the tomato-DC3000 interaction was evaluated and revealed that these compounds may operate as the plant’s first line of defence against pathogen invasion by preventing the production and assembly of a functioning T3SS and reducing the swimming and swarming ability of DC300 due to the loss of flagella . However, the study also showed that DC3000 has developed efflux pumps like MexAB-OprM that prevent the intracellular build-up of flavonoids. Thanks to these efflux pumps further colonisation by the bacteria is promoted by the subsequent secretion of T3 effector (T3E) into plant cells thereby inhibiting continued generation of these antimicrobials. The identified C-glycosylated flavones in oat are also found in other cereals like wheat, rice and maize. Vitexin, isovitexin, orientin and isoorientin are among the most commonly found C-glycosides and are derived from naringenin as flavanone precursor . Known derivatives of vitexin include isovitexin, vitexin-2-O-rhamnoside, methylvitexin, rhamnopyranosyl‐ and vitexin‐2‐O‐xyloside . These compounds are known as radical scavengers and antioxidants, while the antimicrobial activity of apigenin, luteolin, and their C-glucosides were tested and found to generally be more potent against Gram-negative bacteria than Gram-positive ones . Among the triggered alterations, linoleic—and linolenic acid fatty acids were identified. Previous metabolomic studies have highlighted the presence of oxygenated lipids in plant defence in response to biotic stress . These fatty acids are responsible for the biosynthesis of JA and derivatives via the octadecanoic pathway. Jasmonates are signalling hormones produced in response to pathogen attack and causes the plant to adjust its metabolism to produce potent defensive secondary metabolites like phenolics, flavonoids, alkaloids and terpenes . In this study, ethyl 7-epi-12-hydroxyjasmonate glucoside was identified among the discriminatory metabolites in the treated groups of Ps-c and DC3000. This metabolite has been linked to the JA metabolic pathway, and its production was found to be induced when plants detect a microorganism . Alkaloid metabolites are essential for plant defence, particularly as antibacterial substances. A range of amino acids namely aspartate, lysine, tyrosine, and tryptophan (identified as discriminatory, ) are precursors in the synthesis of alkaloids . Due to the presence of nitrogen atoms (proton accepting) and amine hydrogen group(s) (proton donating), these compounds exhibit outstanding biological activities that are mostly attributed to their capacity to form hydrogen bonds with proteins, enzymes, and receptors . Two alkaloids were identified as discriminatory in the treated groups namely hordenine and feruloylserotonin that possibly contributed to the antimicrobial activity against the respective pathovars. Hordenine was only identified as discriminatory for the Ps-c and Ps-t treated groups. Hordenine is a phenethylamine alkaloid with numerous bioactivities, including antibiotic action against microbes, that was first discovered in barley and is now also identified in a variety of grasses. Hordenine was reported to act as a plant allelochemical that can restrict the growth of weeds or guard against pathogen attacks. Additionally, it has also been found to be involved in plant defence responses through the activation of JA-dependent pathways . Feruloylserotonin is the other alkaloid (also a hydroxycinnamic acid amide) that has been identified in Ps-t treated groups. Relatedly, a study on rice showed that the tryptophan pathway is activated after treatment with Bipolaris oryzae , producing serotonin and its amide-conjugates . Moreover, it was reported that serotonin is absorbed into the cell walls of the rice plant tissues and that mutants lacking the elements required for serotonin production and deposition shows an increased susceptibility towards pathogen infection.   In order to improve on defence strategies against pathogens that are constantly evolving, the development of a systems biology approach to elucidate the biochemical and molecular mechanisms underlying plant immune responses has become essential. Oat plants have received very little attention, with minimal studies on how these plants respond at a metabolic level to pathogenic threats. The adaptive metabolic alterations and mechanisms involved in the oat response to pathogen infection with a range of P . syringe pathovars have therefore been covered in this work. The untargeted metabolomics approach was successfully implemented to acquire a deeper understanding of how oat plants defend themselves at a metabolome level against biotic stress. A broad-based chemical defence response was revealed by multivariate data analysis, which also highlighted signature metabolites and discriminatory markers from a range of metabolite classes. The metabolic markers identified in this study show variations in activated oat defences and shed light on the processes and interactions that underlie this plant pathosystem. The study also revealed metabolic pathways that might be responsible for the metabolic changes brought on by inoculation. The observed trend in the metabolite profiles at various time points provided reasonable support that the respective treatments caused defence activation involving the linoleic acid pathway and secondary metabolite biosynthesis pathways. Phenylalanine, tyrosine, and tryptophan metabolism, the phenylpropanoid pathway and the isoquinoline alkaloid biosynthesis pathway, were highlighted as associated with the synthesis of metabolites that belong to the fatty acid, amino acid, phenolic acid, phenolic amide, saponin, flavonoid, and alkaloid classes. Combined, the different biological functions of the secondary metabolites are essential in providing protection against pathogen infections and maintaining it under a variety of environmental conditions. In terms of nonhost resistance where the plant is able to restrict the pathogen from causing severe infection due to the pathogen not being well adapted in terms of effectors to overthrow plant defences, it can be seen that the Dunnart cultivar showed great overlapping responses to inoculation with the nonhost pathogens ( Ps-t , DC3000 and hrcC − ). Overlapping responses can also be noted among the host and nonhost treated plants as can be seen with the activation of similar metabolic pathways. This might be the result of converging signalling pathways but also the relative size of the precursor flux through contributing pathways as well as the timing and extent of the initiating elicitation events. However, some metabolites were found to be present in the host response ( Ps-c ) and absent in the nonhost response ( Ps-t ) (e.g. avenanthramide A, B and O/R) and vice versa in Ps-t but not Ps-c (e.g. feruloylserotonin). By comparing the responses of the hrcC − treated group vs DC3000, it is noted that even though the mutant did not cause visual symptoms, PTI was sufficient to initiate many of the same responses as the host response to Ps-c where great overlap was seen between these treated groups. This suggests that the pathways involved in the synthesis of key defence metabolites (like avenanthramides) are induced as an early response to initial PAMP recognition, even in the absence of subsequent injection of effectors. The greatest difference between the hrcC − and DC3000 treated groups is seen in the presence of phenolics and flavonoids with the hrcC − treatment showing a greater abundance of these metabolite classes compared to that of DC3000. This could be due to the DC3000 secreting effectors into the plant that are able to reduce the mobilisation of small molecule defences beyond the initial PTI response. Understanding what prevents DC3000 and the hrcC mutant from being virulent in oat plants and other non-host plants could be beneficial for future research. Although oat is an economically important food crop, limited in-depth research is currently being conducted. Thus, a disease model including the DC3000 and hrcC mutants’ effect on oat would be a valuable reference for crop breeding and resistance studies. This could provide a greater understanding of how the host manages infection by P . syringae pathovars and could lead to insight into whether effector virulence is target specific in different plants, requiring the use of specific effectors for different plants within a strain’s host range. Disease resistance and the study thereof is crucial for crop defense, and more so the study of non-host resistance since it is considered to be more durable and outlast host-specific resistance. The data presented in this study could also add to growing databases like the SCIPDb  where the molecular and phenotypical response of plants to a combination of stressors are highlighted. Climate change may increase disease risks by altering pathogen evolution and facilitating pathogen spread to new areas. This type of research will therefore become increasingly important when it comes to defining ‘what it takes to be a pathogen’  and as plant disease outbreaks continue to threaten global food security. This study will ultimately provide a basis for insights into the dynamics of the oat metabolome under diverse biotic stresses which, in turn, can be applied in future studies to improve crop resistance and contribute to the development of new strategies for crop improvement endeavours.   S1 Fig Pseudomonas syringae interaction with the plant cell and the respective PAMPs commonly recognised by plant cell receptors. (DOCX) S2 Fig The use of mass spectral fragmentation data for the confirmation of elemental composition and potential structural elucidation. (DOCX) S1 Table Statistical validation of the generated OPLS-DA models from the ESI data of oat seedlings treated with the respective Pseudomonas syringae pathovars ( Ps-c , Ps-t , DC3000 and hrcC mutant). (DOCX)   