Functional Metabolomics—A Useful Tool to Characterize Stress-Induced Metabolome Alterations Opening New Avenues towards Tailoring Food Crop Quality

: The breeding of stress-tolerant cultivated plants that would allow for a reduction in harvest losses and undesirable decrease in quality attributes requires a new quality of knowledge on molecular markers associated with relevant agronomic traits, on quantitative metabolic responses of plants to stress challenges, and on the mechanisms controlling the biosynthesis of these molecules. By combining metabolomics with genomics, transcriptomics and proteomics datasets a more comprehensive knowledge of the composition of crop plants used for food or animal feed is possible. In order to optimize crop trait developments, to enhance crop yields and quality, as well as to guarantee nutritional and health factors that provide the possibility to create functional food or feedstuffs, knowledge about the plants’ metabolome is crucial. Next to classical metabolomics studies, this review focuses on several metabolomics-based working techniques, such as sensomics, lipidomics, hormonomics and phytometabolomics, which were used to characterize metabolome alterations during abiotic and biotic stress in order to ﬁnd resistant food crops with a preferred quality or at least to produce functional food crops.


Importance of Metabolomics for Agricultural Research
In the environment, plants are often exposed to an enormous number of biotic and abiotic stress factors, such as pathogen and insect infestation as well as extreme temperatures, drought, salinity, pollutants, heavy metals or nutritional deficiencies. That leads to harvest or quality losses, such as the formation of a pronounced bitter off-flavor and sometimes causes toxicological problems as well as huge global economic losses [1]. The breeding of stress-tolerant cultivated plants that would allow for a reduction in harvest losses and undesirable decrease in quality attributes requires a new quality of knowledge on molecular markers associated with relevant agronomic traits, on quantitative metabolic responses of plants on stress challenges, and on the mechanisms controlling the biosynthesis of these molecules [2]. Thereby, knowledge of the biologically active metabolites, especially secondary metabolites as well as metabolic networks of primary and secondary metabolites affected by plant stress conditions, is essential.
The metabolome of plants represents the complete set of low-molecular weight metabolites (such as primary metabolites including all necessary metabolic intermediates, hormones and other signaling molecules, as well as secondary metabolites) in a given organism, a biological cell, tissue, Figure 1. Functional metabolomics-a useful tool to discover metabolome alterations during abiotic and biotic stress to find the perfect functional pheno-/genotype to produce functional food.
Next to the genetic code, environmental stress especially influences the metabolome composition of plants. Stress conditions are sensed by the plant, which then activates a network of signaling pathways, including the participation of phytohormons, and leads to changes, on the one hand, in the primary metabolism and, on the other hand, in the up-regulation of phytochemicals [20].
For example, during fungi infestations, antifungal phytochemicals, especially metabolic pathways leading to isoprenoids, phenylpropanoids, alkaloids, fatty acids, and polyketides are upregulated [21]. Many phytoalexins, which partly have been described as flavor-active, nutritional components and as a source for development of health-promoting food products, have been well documented in crop plants in the field of plant defense [21]. In contrast, primary metabolism is documented to be mainly influenced by abiotic stress challenges [5,[20][21][22]. For instance, in Arabidopsis leaves, next to the amino acid profiles, oligosaccharides, γ-amino butyrate and the metabolites from the tricarboxylic acid cycle are influenced during drought conditions [23].
Analysis of those changes in the molecular composition of plant material by means of mass spectrometry (MS) or nuclear magnetic resonance spectroscopy (NMR) based techniques is called metabolomics, which can profile the impact of time, stress, nutritional status and environmental perturbation on hundreds of metabolites simultaneously. It is an essential technology for functional genomics and systems biology which can visualize and answer questions about biological systems Figure 1. Functional metabolomics-a useful tool to discover metabolome alterations during abiotic and biotic stress to find the perfect functional pheno-/genotype to produce functional food. Therefore, an important aspect of the use of crop plants for human consumption is the characterization and determination of their nutritional value as well as their individual flavor-active, health-promoting, non-toxic, functional constituents. The aim of crop plant research should be to find an optimal functional pheno-/genotype as a golden standard with high resistance and best possible quality criteria for future breeding experiments (cf. Figure 1). Analytical marker compounds can help to support molecular breeding efforts to obtain phenotype knowledge-based plant populations, no matter which analysis method is used, bulk segregant analysis, mapping by sequencing or quantitative trait loci (QTL) analysis [10,11]. In sum, metabolomics, especially the characterization of functional relevant metabolites, also called functional metabolomics, can be used as a tool for functional genomics [2,9,12,13]. For example, integration of metabolomic data with genomics and transcriptomics, especially that obtained from transcriptome coexpression network analysis, helped to decode several functions of genes and to gain a deep insight into the biological processes of the model plant Arabaidopsis [2,9]. Although more studies are desirable, at the moment only manuscripts are available in which secondary, functional metabolites of non-model plants are characterized [2,9,12,13]. Next to studies using potatoes, strawberries, carrots or Brassica napus [9,[14][15][16][17][18], Rambla et al. [18] performed a QTL analysis using a cross between Solanum lycopersicum and a red-fruited wild tomato species Solanum pimpinellifolium to relate 102 QTLs to 39 corresponding volatile organic compounds, which partly contribute to the typical flavor impression of tomatoes. In addition, by means of genome-wide association studies 323 associations among 143 single nucleotide polymorphisms and 89 metabolites were identified, including non-volatile, functional, secondary metabolites, which can be used for new human drug production [19].
Next to the genetic code, environmental stress especially influences the metabolome composition of plants. Stress conditions are sensed by the plant, which then activates a network of signaling pathways, including the participation of phytohormons, and leads to changes, on the one hand, in the primary metabolism and, on the other hand, in the up-regulation of phytochemicals [20].
For example, during fungi infestations, antifungal phytochemicals, especially metabolic pathways leading to isoprenoids, phenylpropanoids, alkaloids, fatty acids, and polyketides are up-regulated [21]. Many phytoalexins, which partly have been described as flavor-active, nutritional components and as a source for development of health-promoting food products, have been well documented in crop plants in the field of plant defense [21]. In contrast, primary metabolism is documented to be mainly influenced by abiotic stress challenges [5,[20][21][22]. For instance, in Arabidopsis leaves, next to the amino acid profiles, oligosaccharides, γ-amino butyrate and the metabolites from the tricarboxylic acid cycle are influenced during drought conditions [23].
Analysis of those changes in the molecular composition of plant material by means of mass spectrometry (MS) or nuclear magnetic resonance spectroscopy (NMR) based techniques is called metabolomics, which can profile the impact of time, stress, nutritional status and environmental perturbation on hundreds of metabolites simultaneously. It is an essential technology for functional genomics and systems biology which can visualize and answer questions about biological systems [3,24]. Kushalappa and Gunnaiah [25] already predicted that while in the past, genetic tools were mainly used for crop improvement on yield, flagged as the green revolution, in the next few decades functional genomics, which allows the identification of genes that enhance crop yields without losses in nutritional value while minimizing loss due to stress, will be the main topic. Metabolomics is especially helpful to discover metabolites up-regulated during stress conditions [2]. These metabolites then can be tested for their biological functions. Thereby, the combination and understanding of information received from genomics, transcriptomics, proteomics and metabolomics studies is a complementary tool for functional genomics and systems biology investigations. Although nowadays, next to the human genome [26], the genomes of several plants such as Arabidopsis thaliana [27], rice [28,29], tomato [30], and barley [31] have been decoded by means of automated nucleotide sequencing, knowledge about the regulation of gene expression, about the action of these genetic products, and especially about the metabolic networks resulting from catalytic proteins, is rather fragmentary [3].
Besides aiming to avoid quality or yield losses induced by abiotic or biotic stress conditions, there are exciting new strategies of using plant stress as a way of producing phytoalexin-enriched functional foods [32] or to enhance the flavor quality of crop plants by stress or elicitor challenges to produce new functional foodstuffs with high quality [33]. Dixon et al. [3] additionally suggest that the long-term goal should be the application of new knowledge of functional metabolites for "rational custom-designed breeding by classical methods as well the application of genetic engineering techniques to improve and develop new aromatic plants". Therefore, publicly open databases, containing all spectroscopic, spectrometric and bioactivity data as well as the biosynthetic pathways of worthy metabolites are necessary in the future [3].
Targeted as well as non-targeted metabolomic-based working techniques, such as phytometabolomics, sensomics, lipidomics or hormonomics can help to unravel metabolic pathways, signal transmissions as well as metabolites with different bioactivity which can be simultaneously up-or down-regulated by stress conditions. After decoding marker compounds by means of these non-targeted approaches, marker metabolites can be isolated from the plant material and characterized by means of 1/2D-NMR or MS based structure identification techniques as well as by biological assays (e.g., using antifungal susceptibility test systems). Although the last-named structure identification techniques and the testing of the metabolite's bioactivity do not belong to the metabolomics approach, its application is mandatory to understand biological processes.

Analytical Techniques Used to Characterize Stress-Induced Metabolome Alterations in Plants
Enormous technological progress has been achieved since the early days of plant metabolome characterization. Thin-layer chromatography (TLC) together with gas chromatography coupled to a flame ionization detector (GC-FID) have so far been applied to characterize primary and secondary metabolites in plant materials. Nowadays, besides NMR, high-throughput, sensitive techniques using gas chromatography, liquid chromatography, ultra-high performance liquid chromatography/time-of-flight as well as capillary electrophoresis coupled to mass spectrometry (GC-MS, LC-MS/MS, UHPLC-TOF-MS or CE-MS), together with bioinformatics tools, are state-of the art working techniques [3,13]. With the development of mass spectrometry imaging approaches, including matrix-assisted laser desorption (MALDI) or desorption electrospray (DESI) ionization procedures coupled to high resolution mass spectrometers, it is possible to perform metabolomics analysis in situ, as well [34]. Although NMR-based techniques are used in plant metabolomics, it is less common in plant-stress metabolomics. However, new working techniques including NMR-technology, such as differential off-line LC-NMR (DOLC-NMR), show promise for characterizing metabolites upor down-regulated during biotic or abiotic stress influences on plants [35].
Although the term plant metabolomics is defined as the identification and quantification of all low-molecular weight substances in an organism, at a defined point of harvest time and development stage and in a given organ, tissue or cell type, there is no single work-up protocol or standard technique available to detect all metabolites using only one platform. Therefore, several analytical separation and detection techniques are usually combined to visualize abiotic and biotic stress induced metabolome alterations of different plants induced by abiotic and biotic stress [4].
In the last four decades, different methods have been established for the targeted and non-targeted characterization of plant metabolites (cf. the following reviews: Shualev et al. [20]; Urano et al. [22]; Obata and Fernie, [5]; Arbona et al. [4]; Kumar et al. [2]). Nowadays, we differentiate between the non-targeted, sometimes also called untargeted, and the targeted metabolomics approach. While the non-targeted approach represents a comprehensive analysis of all measurable compounds, including unknown metabolites, targeted metabolomics is the measurement of defined groups of chemically characterized and biochemically annotated substances [36]. For both, targeted and non-targeted analyses GC-MS, UHPLC-TOF-MS, as well as LC-MS/MS are the most common techniques to characterize stress-induced metabolome alterations in plants [5].
In particular, primary metabolites, such as carbohydrates, amino acids or fatty acids and all volatile compounds are mainly characterized by means of targeted GC-MS with or without derivatization [37]. Due to the fact that standard compounds are commercially available and spectrometric data is published in several public libraries, such as NIST (National Institute of Standards and Technology), the identification and exact quantification of primary metabolites by means of GC-MS analysis is straightforward. In contrast, structure identification of secondary metabolites exclusively by means of MS data is rather challenging [4,38]. Metabolomics approaches mainly identify secondary non-volatile compounds by means of non-targeted UHPLC-TOF-MS techniques. In contrast to GC-MS, LC-MS libraries for structure identification are less developed, as instrument-type dependent mass spectra and MS fragmentation patterns, as well as retention time shifts depending on used LC columns complicate the comparison of structure identification results using different instruments [39]. In the future, a worldwide concept of raw data processing, extensive mass spectral libraries, and powerful database management systems that can store and provide both raw and meta data is desirable.
At the moment, biological active compounds are unequivocally identified by means of MS and 1/2D NMR experiments after their isolation by means of preparative HPLC techniques. But the MS spectra obtained has mostly not been shared with the metabolomics community due to the absence of a standardized data-handling concept. Once identified, sensitive quantitation methods using stable isotope dilution assays (SIDA) are carried out by means of LC-MS/MS or GC-MS techniques. Therefore, nowadays, determination of the exact concentration of the isotopically labelled internal standards can best be achieved by means of quantitative NMR spectroscopy (q NMR) [40].
Although the aforementioned analytical metabolomics techniques nowadays are reliable and very sensitive, they still have their limitations regarding small sample volumes or low quantities of metabolites (partly in the nanomolar range). This is especially so when metabolomics is used to answer specific questions in the field of functional genomics or systems biology to get a comprehensive understanding of plant growth, development, defense or productivity, and knowledge about metabolic pathways and networks at cellular levels is required [9,12,[41][42][43]. In addition to single cell-cell genomics, transcriptomics or proteomics, several groups developed plant single-cell and single-cell type metabolomics during the last decade [43][44][45][46]. While using single-cell metabolomics, the metabolome of one individual cell, such as a single epidermal cell, is analyzed, while in single-cell-type metabolomics the metabolome of a population of cells, such as guard cells or trichomes, is characterized [43]. Whereas in studies using single animal cells, more than 300 compounds could, nowadays, be visualized by means of capillary electrophoresis-electrospray ionization-time-of-flight-mass spectrometry (CE-ESI-TOF-MS) methods. The amounts of metabolites characterized in smaller plant-based single-cells is still much lower [42,43]. To gain a better understanding of biological systems and characterize metabolic pathways in single cells or in situ, special attention should, in the future, be given to the development of more sensitive MS apparatuses. This is especially the case in the development of more sensitive visualizing MS techniques, such as MALDI or DESI imaging experiments, and in the improvement of extraction protocols [13,43]. Newly developed derivatisation strategies are also promising strategies to improve the MS sensitivities in future. Next to limitations in measuring low sample amounts, data interpretation sometimes is challenging.
A major challenge for combining all omics-techniques is the integration of detailed metabolomics data into transcriptomic and proteomic profiling data sets [3,9,12,13,20]. Although Arbona et al. [4] only recently reviewed the role of metabolomics as a tool to investigate abiotic plant stress in the context of systems biology, and how metabolomics can be integrated into other omics data sets, there is still a lack of knowledge. Furthermore, Urano et al. [22] summarized arguments that holistic omics analysis and their data interpretation are essential to identify the broad function of metabolite regulatory networks during the responses to plant stress. Therefore, bioinformatics knowledge and tools should be used more often in the future to combine the information from both research fields. Although databases and pathway viewers (e.g., an overview is given by Bino et al. [12] or de Souoza et al. [47]) have been further developed in recent years, inter alia Sumner et al. [13] also emphasized that standardized databases integrating metabolomics with other global-omics data are required.

Phytometabolomics-From Plant Stress to Metabolic Response
Plants continuously encounter various biotic and abiotic environmental stresses during their growth and development phases. While biotic stress is caused by pathogens, parasites, predators, and other competing organisms, abiotic stress arises from inappropriate levels of physical components in the environment, such as temperature or water extremes [21]. Both stress factors lead likewise to yield and quality losses in plant crops used for human nutrition. The aim is, on the one hand, to avoid those immense quality and yield losses to ensure feeding the world's growing population and on the other hand, to develop new functional foods and to reduce insecticides and pesticides used in agriculture. Therefore, knowledge of molecular networks and pathway activation of plants during their stress response is essential.
Only recently Mithöfer and Boland [48] asked the question "Do you speak chemistry?". They highlighted that plants always respond by emitting chemical compounds to signal their environment. Moreover, they pointed out that plants use metabolites as chemical sensing and communication systems. On the one hand, signaling molecules and phytoalexins are formed, on the other hand it is suggested that plants "cry for help" by producing secondary metabolites as indirect defense system. Those cry-for-help compounds, usually belonging to the group of volatile organic compounds, attract predators or herbivores helping the plant to get rid of other stress triggers [48,49]. Therefore, it is important to notice that stress-specific as well as common metabolites, so-called generalist and specialist, can either be formed or their amounts up-regulated [5].
The understanding of metabolic stress pathways is just possible if comparative metabolite characterization by means of metabolomics analysis of stressed and non-stressed plants followed by structure identification experiments is performed. In the following three sub-sections an initial overview of the literature about biotic and abiotic stress studies is undertaken.

From Abiotic Plant Stress to Metabolic Response
It is well expected that climatic changes due to the global warming will lead to massive abiotic stress factors influencing the metabolome spectra of all plants in the future. Next to water stress, including soil flooding and drought stress, salt stress, temperature stress, light stress, sulfur and phosphorus stress, oxidative stress, as well as heavy metal stress will play important roles for food crop producers. To understand and avoid quality and yield losses induced by abiotic stress factors, several metabolomics studies have been carried out (cf. the following reviews: Shulaev et al. [20]; Kaplan et al. [21]; Urano et al. [22]; Obata and Fernie [5]; Arbona et al. [4]).
Among them, Arbona et al. [4] already summarized in their review "Metabolomics as a tool to investigate abiotic stress tolerance in plants" that the plants' metabolic responses to abiotic stress can mainly be found in the primary metabolism. Next to metabolic responses in mono-and disaccharides, sugar alcohols, amino acids, especially proline, polyamines and members of the tricarboxylic acid cycle were found to be influenced by abiotic factors. Thereby, the carbohydrate metabolism in particular is directly influenced by stress conditions, because during stress instead of glucose, plants use fructanes and starch as an energy source [50].
In their review, Obata and Fernie [5] compared several metabolic fingerprints of Arabidopsis leaves influenced by dehydration, salt, heat and cold, high light and sulfur limitation, ultraviolet (UV) light quality change, low nitrogen amounts, and potassium limitation with each other. Thereby, they highlighted that specific compounds and compound classes are generally accumulated during abiotic stress factors (except of light treatment), such as sucrose, raffinose, proline, other branched chain amino acids or γ-aminobutyric acid (GABA). However, the amounts of such up-regulated "generalists" varies between stress factors. They concluded, that the stress-specific plant response is the result of an inhibition or activation of a defined metabolic pathway. In particular, enzymes activities are known to be influenced by temperature or ion concentrations in plants. By contrast, other metabolites are accumulated only under specific conditions, such as trehalose.
Next to the carbohydrates, the lipid composition of plants is additionally influenced by abiotic stress factors. Lipids in plants are a crucial and diverse class of biomolecules that forms the so called plant lipidome. Lipidomics is defined as the identification and quantification of all lipids within a biological system at a defined stage of development [51]. Within the field of omics-analysis lipidomics represents a sub-unit of targeted and non-targeted metabolomics approaches, mainly using MS-and NMR-based techniques. The development of lipidomics and its increased promise in systems biology over the past two decades has been the subject of several reviews (e.g., Blanksby and Mitchell [34]). Changes in the lipidome, induced by modifications of its biosynthesis, regulation, adaption, remodeling, function, role, and interaction, as well as membrane lipid remodeling are relevant responses from plant cells to counteract biotic and abiotic stress challenges [51]. In particular, the oxidation of polyunsaturated fatty acids is one of the most fundamental reactions in lipid chemistry, which for example leads to off-flavor (perceived rancidity) or to phytohormon precursor formations [52]. Different groups of lipids, such as fatty acids, phosphatidic acids, inositol phosphates, diacylglycerols, oxylipins, shingolipids, and N-acetylethanolamine are involved in signaling systems during stress conditions as well [53][54][55][56][57][58]. Those compounds are shown to be directly synthesized after abiotic stress influences by a wide range of enzymes, such as fatty acid amide hydrolases, phospholipases, acyl hydrolases, diacylglycero kinases or phytoshingosine kinases [58]. In addition, it is known that in response to abiotic stress factors, lipids migrate into cell walls to repair damage through membrane remodeling [58][59][60]. Moreover, a combination of transcriptomics and lipidomic profiling showed that cold stress induces the prokaryotic pathway and suppresses the eukaryotic pathway for glycerolipid biosynthesis [61]. Next to cold stress, light and temperature stress, as well as nutrient starvation as typical abiotic stress factors were studied by means of lipidomic approaches in Arabidopsis (for details cf. the following review: Tenenboim et al. [51]). Only recently, lipidomics studies led to the clarification of the role of leaf lipids in thyme plant response to drought stress [58].
In the future, more stress-related lipidomic studies are expected, due to the fact that in nutritional research lipid profiling also plays a crucial role [51]. The food industry, for example, shows great interest in the use of lipidomics to characterize defined crop plant genotypes. As a replacement of palm kernel fat by a natural alternative, new plant genotypes having the same triglyceride pattern are needed. Furthermore, although most genes annotated as lipid-related are functionally characterized, their corresponding metabolites are not assigned. That would be an interesting field of research in the future, as well.
Next to numerous papers and reviews which highlight metabolome alterations induced by one single stress factor, Karin Köhl's review [62] gives a very nice overview of the metabolites influenced by combined abiotic stress effects in crops. Due to the fact that plants often have to deal with more than one stress factor, further metabolomics studies in crops which take several stress factors into account would be useful in the future [62].
In sum, although integrated omics data of model plants such as Arabidopsis have markedly increased our knowledge of understanding the mechanisms involved in plants' response to various abiotic stress factors [22,63,64], knowledge about secondary metabolome alterations in abiotic stressed food crops is rather fragmentary. Only few studies indicate abiotic stress influences the amounts of phenolic compounds, glucosinolates, carotinoides terpene derivatives and phytohormons in plants [4,51]. Urano et al. [22] emphasize again and again that phytohormon levels, especially those of abscisic acid, vary a lot in different plant compartments. Next to knowledge gaps in the stress influence on the secondary metabolome, more studies are needed to characterize stress combinations as well [5].

From Biotic Stress Metabolomics to Metabolic Response
During biotic stress exposure, plants use qualitative and quantitative measures to resist pathogen attacks [65]. While, in the past, the qualitative resistance, based on monogenetic inheritance, has been successfully used to elite cultivars to improve resistance, the quantitative resistance, which is presumed to be durable, non-race-specific, and effective against a large bouquet of pathogens, is mainly unknown [25]. However, a huge number of quantitative trait loci (QTLs) associated with qualitative resistance against pathogen-associated molecular patterns (PAMs), found for example in soy beans (Glycine max), and with quantitative resistance, found for blast in barley (Hordeum vulgare), or powdery mildew in wheat, have already been identified for biotic stress challenges [25,[66][67][68][69][70]. However, both the mechanisms of resistance and the metabolites up-regulated during plant stress controlled by the QTLs are mainly unexplored [66]. In particular, targeted metabolomics approaches enable the possibility to monitor marker compounds up-regulated during plant stress. They can then be isolated, identified and biologically characterized, for example by application of an antifungal test.
Non-targeted metabolomics in comparison with bioinformatics and database research can also give an idea of which metabolic pathways are activated during a plant's stress response. Moreover, targeted analysis of known bioactive, e.g., anti-microbial or anti-fungal, compounds can help to identify QTLs, by using targeted metabolomics techniques for phenotyping experiments in breeding crossing lines. Thereby, alongside breeding improvement, metabolomics helps in the understanding of biological functions and the study of host-pathogen interactions [25].
It is already known from plant-pathogen interaction studies that signaling molecules, such as ethylene, salicilic acid, jasmonic acid or inositol, as well as antifungal/antimicrobial phytochemicals and cell wall compartments are formed or their concentrations up-regulated during plant stress [25]. Numerous studies already indicate that a complex network of secondary metabolites is influenced by biotic stress challenges. Next to phenylpropanoids, including flavonoids [34,50,51,71], alkaloids [72,73], terpenoids [33], and fatty acids [74] are also known to be up-regulated [25]. In addition, lipids act under pathogen or herbivore attack as mechanical barriers. The plants' wax layers form the first line of defense against pathogens and herbivores [51,75]. Additionally, some lipids act as signaling molecules, namely oxylipins and jasmonates, which participate in immune response cascades. Phytoalexins, a group of fatty acid degradation products, have antimicrobial and antifungal activities [51]. In addition, only recently Caroline Gutjahr's lab could show that fatty acids, like carbohydrates, are also transported from the plant host to fungi [76].
Moreover, as plants produce a diverse array of more than 100,000 low molecular weight natural products [32], countless studies are available using targeted quantifications of phytoalaxins in plants.
In 2009, Boue et al. [32] published a forward-looking paper called "phytoalexin-enriched functional foods" in which they proposed to use stress conditions to enhance the amounts of phytoalexins in food plants to produce functional foods. Several phytoalaxins are known for their health beneficial properties, including antioxidant activities, anti-inflammatory activities, cholesterol-lowering abilities, and anticancer activities [32]. In particular, flavonoids, which are ubiquitous in many food plants, have been linked to important health-promoting activities. For example, consumption of legumes, especially soy, have been linked to the reduction of cancer risks and coronary heart diseases simply due to its high flavonoid content [32]. Therefore, besides resistance and yield stability, functional food ingredients from the phytoalexin family could also be a possible target for plant breeders. In addition, mild stress conditions could possibly help to enhance the amounts of bioactive secondary metabolites, thereby producing functional foods [32].
Next to aspects in nutrition, diet and health, food and environmental safety can especially be monitored by means of metabolomics [3]. Toxic compounds can, for example, be formed as a stress response by plants, such as furanocoumarins in celery and glycoalcaloides in potatoes, or by fungi contamination, which can produce human toxic mycotoxins [77]. Therefore, versatile, sensitive, reliable and fast-targeted LC-MS/MS and GC-MS methods using stable isotope dilution techniques have been developed in the past to monitor food quality [78,79].
In sum, in the future more versatile and sensitive multiple-targeted MS quantification methods to evaluate food and crop quality are required.

Functional Phytometabolomics-Characterization Approach of Plant Stress Metabolites
In conclusion, the science-driven breeding of stress-tolerant cultivated plants that would allow for a reduction in harvest losses and undesirable decrease in quality attributes requires a new quality of knowledge on molecular markers associated with relevant agronomic traits, on quantitative metabolic stress responses of plants, and the mechanisms controlling their biosynthesis. The field of "functional phytometabolomics", using targeted and non-targeted MS or NMR techniques to quantitatively assess key metabolome alterations in plant-derived crops and foods induced by biotic stress challenges as well as abiotic stress conditions, is, therefore, a promising field of research.
In the phytometabolomics approach, metabolites up-regulated during stress challenges are visualized by means of working techniques used in the field of metabolomics. Markers previously not published in literature could e.g., be visualized by retention time and mass to charge ratio in UHPLC-TOF-MS analysis. They are isolated in purities higher than 98% from the plant material using medium pressure liquid chromatography (MPLC) and preparative HPLC techniques and are identified by means of LC-MS, LC-MS/MS, UHPLC-TOF-MS and 1/2D NMR experiments. Testing of different biological activities (anti-fungal, anti-bacterial, anti-oxidant activities etc.) of those compounds allows initial insights into their biological functions. To translate the knowledge on how stress-resistant traits master their successful defense against stress conditions into breeding programs, genotype-specific metabolome alterations have to be characterized. Subsequently, the gene clusters controlling the biosynthetic pathways of key stress metabolites have to be identified by means of genome-wide association and QTL mapping studies. Just to highlight two examples: Firstly, Matsuda et al. [19] characterized several bioactive flavon glycosides, which can be used for new human drug discovery in rice by means of a metabolom genome-wide association study. Thereby, they highlighted "that one plant species produces more diverse phytochemicals than previously expected, and plants still contain many useful compounds for human applications". Secondly, only recently Rambla et al. [18] used a QTL analysis, to show that 102 QTLs correspond to 39 different volatile organic compounds, including flavor active key metabolites, in tomatoes. This research will help to navigate breeding programs and to optimize post-harvest treatment of plant-derived food products from producer to consumer/processor towards the production of high-quality food products.

Phytohormone Profiling by Means of Plant Hormonomics
Phytohormones are a class of low molecular weight, structurally diverse, but highly bioactive compounds in plants. They act as chemical messengers, triggering and controlling physiological processes during plant growth and development (e.g., cell elongation, regulation of apical dominance, vascular differentiation, fruit development, latal and adventitious root formation) as well as in response to abiotic and biotic stress conditions [2,80,81]. Next to ethylene, auxins, cytokinins, brassinosteroids, gibberellins, jasmonates, salicylates, polyamines, abscisates and signal peptides, strigolactones are part of the phytohormone family [81,82]. During stress exposure those phytohormone classes interact with each other by means of synergistic or antagonistic cross-talks, resulting in each other's biosynthesis or up-regulation response [83]). In the past, several studies provided evidence that plant hormones are necessary for plants to adapt during stress conditions, especially, during abiotic stress factors by mediating a wide range of adaptive responses [83]. They play a key role in the plant's intricate signal networks, often immediately altering gene expression by inducing or preventing the degradation of transcriptional regulators via the ubiquitin-proteasome system [83,84]. Thanks to Kumar et al. [2], metabolic engineering of phytohormons can be used to improve quality and stress tolerance of crops. Although the analysis of plant hormones, such as auxins, especially their quantification by means of SIDA-LC-MS/MS revealed significant insights in their tissue-and cell-type-specific analysis, their distribution profiles in plant organs, tissues, and cells still remains elusive [81]. A wide variety of targeted GC-MS and LC-MS/MS methods have been published dealing with the quantification of exact amounts of single members of the phytohormone family (e.g., cf.: Porifiro et al. [85]; Novák et al. [81]). However, in the last two decades, only one method was described to characterize those physiologically important molecules in their network interactions in one single run [86]. Only recently, Ondřej Novák and his team developed a new, versatile and sensitive targeted UHPLC-MS/MS method, which enables the simultaneous quantification of 101 phytohormone-related metabolites (phytohormones and their precursors) in less than 20 mg plant material [86]. This allows the characterization of the majority of known phytohormones as well as their biosynthetic precursors, that regulate diverse processes in plants by intricate signaling networks. With this newly developed, targeted metabolomics approach, so called "plant hormonomics", Šimura et al. [86] were able to detect 45 and quantify a total of 43 endogenous compounds out of the 101 phytohormones in both root and shoot samples, in salt-stressed and non-stressed, 12-day-old Arabidopsis thaliana seedlings. Subsequent multivariate statistical analysis cross-compared with data obtained from transcriptomic studies enabled the identification of the main phytohormones involved in the adaption of Arabidopsis thaliana to salt stress. The multivariate statistical analysis revealed that 23 of the quantified 43 metabolites significantly differed between the salt-stressed and non-stressed roots. In contrast, the shoots and the roots differed in the concentrations of 15 compounds. Those hormone profiles obtained were cross-compared with transciptomic data. Well in line with findings from Rhy and Cho [87], who found that jasmoic acid and abscisic acid promote salt tolerance in Arabidopsis, Šimura et al. [86] observed increased levels of abscisic acid, its oxidation products phaseic acid, dihydrophaseic acid, and jasmonic acid in salt-stressed root samples. In addition, different responses of gibberillinic acid derivatives, especially its active form GA4, and its transcriptomic data, were found by Šimura et al. [86]. In sum, this newly developed multiple parallel analysis of phytohormones, called "plant hormonomics" enables the real-time profiling of hormone networks of large collections of phytohormones, their precursors, transport forms and degradation products in single stressed or non-stressed samples [86,88].
Although, Šimura et al. [86] were already able to quantify the main phytohormones in less than 20 mg of salt stressed Arabidopsis samples, special derivatisation reactions of the analytes prior to their analysis will conceivably lead to an increase in sensitivity and selectivity during MS analysis in the future [85,[88][89][90]. For example, the response/sensitivity during LC-ESI-MS/MS analysis of metabolites only present in very low concentrations in plants, like indole-3-acetic acid, a member of the auxin family, can be tremendously increased up to 200-fold after methylation [89]. In addition, several aldehyde trapping derivatisation reagents are known to enhance the sensitivity during ESI-MS analysis of biological mixtures [91]. In future, it could be an extremely important and advantageous step in hormonomics to analyze a combination of derivatised and non-derivatised phytohormone classes, simultaneously in one method.

Sensomics-A Phenotyping Tool to Characterize Crops Flavor Impression
In general, more than 83% of consumers say that flavor influences their decision the most when purchasing any kind of food products or beverages [92]. Therefore, it is no surprise that, although aided by visual inspection, the final recognition and quality evaluation of food crops made by consumers are mainly mediated by its flavor perception. Human flavor perception is induced by the interaction of volatile odor-active and non-volatile taste-active molecules with~390 odorant receptor proteins located in the regio olfactoria in the nose and~40 taste receptor proteins on the tongue [93]. To meet consumers' demand for continuously available, fresh foods with a premium quality, the "flavor blueprint" of a golden standard, which is a combinatorial code of the entire set of odor-and taste-active metabolites in their natural concentrations, has to be known [94]. To decode all flavor-active molecules of a crop, the so-called sensometabolome, several high-end working techniques, including combinations of state-of-the-art metabolomics and chromatography approaches in combination with human sensory science experiments, such as aroma or taste dilution analysis, have been developed in the past. This so called molecular sensory science or sensomics approach has especially been shaped by Peter Schieberle's and Thomas Hofmann's working groups during the last three decades [93,94]. Today, it is well accepted that the presence of certain structural elements, so-called olfactophores and gustophores, as well as specific concentrations exceeding the sensory thresholds are important prerequisites of low molecular weight metabolites to become flavor-active [94]. Although a total of more than 10,000 volatiles occur in crop plants, the use of the so-called sensomics approach gave evidence that the typical aroma impression of a crop-based foodstuff is caused by a limited number of aroma-active volatiles [93]. In conclusion, only a surprisingly small group of so called "key odorants in food" contribute to the typical aroma profile of a crop plant. For example, only two compounds, namely ethyl (2S)-2-methylbutanoate and 1-(ethylsulfanyl)ethane-1-thiol, are necessary to explain the typical aroma impression of durian fruits [95], and to mimic the typical aroma impression of mangos only eight compounds are required [96]. While only a small number of aroma-active compounds contribute to the key odorants in food and interact with a huge number of odorant receptors, several thousands of non-volatile taste-active compounds are already known. In particular, the application of taste dilution analysis followed by dose/activity considerations led to the discovery of many bitter, sweet, umami, pungent, astringent, salty, or sour sensing molecules in several plants, such as carrots [97,98], cocoa [99], asparagus [100][101][102], pepper [103], red currants [104], tea [105], stevia [106] or spinach [107] in the past. It is nowadays well accepted that not only a single flavor-impact molecule, but a combinatorial code of multiple odor-and taste-active key compounds, each in its specific concentration, reflect the chemosensory phenotype and trigger the typical flavor profile of food products. In particular, for flavor improvement of crops, the analysis of aroma-and taste-active compounds presents a major challenge for flavor improvement of crops [94]. It has been shown in the past, that the biosynthesis of several key odorants in food and tastants is controlled by genes, the expression of which is altered or even induced by biotic or abiotic stress challenges. But the primary target of crop production then was field performance, yield, and storage characteristics, while ignoring quality traits, such as the flavor code [33,108].
On the one hand, induced by several biotic and abiotic stress factors during growth in the field as well as during post-harvest storage, the attractive sensory quality of miscellaneous crop plants is hindered by sporadic off-flavors, which is often the reason for consumer complaints and therefore a major problem for plant processors. On the other hand, mild stress factors can lead to an increase in in the concentration of flavor-active constituents and hence to a more intensive desirable aroma or taste impressions [33,108].
Moreover, abiotic factors, such as mechanical stress, are reported to increase the amounts of the key bitter tastants, members of the C17-polycetylenes, present in native carrots (Daucus carota L.). The increase causes a perceived bitter off-flavor, occurring especially often during the production process for infant diet carrot products [97,98,[109][110][111][112][113][114][115]. In addition, a decrease in flavor quality accompanied by an increase in bitter taste has also been reported in raw hazelnuts (Corylus avellana L.) upon biotic stress challenges, such as upon infection by bugs, belonging to the hemipteran family, like Gonocerus acuteangulatus and Coreus marginatus [116].
Besides the non-volatile, taste-active metabolites, aroma-active compounds are also known to be influenced by individual stress factors [33]. In particular, studies dealing with tea (Camellia sinensis), which is enjoyed as freshly brewed green, black, oolong, or decaffeinated tea infusion, and its quality changes induced by stress are available in literature. For example, the tea green leafhopper (Empoasca (Matsumurasca) onukii Matsuda) attacks, at pre-harvest stages, can decisively influence the unique aroma quality of tea leaves as a result of the upregulation of the linalool synthase (CsLIS1 and CsLIS2) causing higher concentrations of the key odorant, linalool [33,117]. Besides linalool, the formation of key odorants, such as (E)-nerolidol, can also be influenced by a combination of low-temperature stress and mechanical damage [33,118]. Therefore, Wüst [33] concluded that the tea-related findings illustrate how the use of the stress response of plants within the sensometabolome can lead to an improvement of flavor of agricultural products. In addition, he pointed out, that the plants' contact with stress elicitors, such as methyl jasmonate instead of the actual biotic or abiotic stress factors could also lead to "stress induced" flavor improvement [33,118].
In order to gain a more comprehensive knowledge on the chemical mechanisms involved in quality changes of cultivated crop plants in response to biotic or abiotic stress challenges or to improve the flavor quality by the application of moderate, well controlled stress, numerous volatile, aroma-active as well as non-volatile taste-active key metabolites, the so-called sensometabolites, of stressed and non-stressed plant genotypes should be comparatively characterized by means of a fast and robust high-throughput GC-MS systems with high peak separation capacity and sensitivity, such as comprehensive two-dimensional gas chromatography/time-of-flight mass spectrometry (GC × GC/TOF-MS) or UPLC-TOF-MS metabolic profiling analysis in the future. This strategy aims at reducing the flavor deficiencies in modern commercial varieties as a "green" alternative to genetic engineering. The workflow for a successful implementation of this approach-from the identification of key odorants by molecular science techniques to the investigation of mechanisms controlling their biosynthesis-is complex and calls for interdisciplinary research [33,108].

Conclusions
Plants continuously encounter various biotic and abiotic environmental stresses during their growth and development phases, which leads likewise to yield and quality losses. To avoid those previously listed economical losses, to reduce insecticides and pesticides used in agriculture, to ensure feeding the world's growing population, and to develop new functional foods knowledge about molecular networks and the pathway activation of plants during their stress response is essential.
This review has emphasized the importance of metabolomics-based working techniques to discover metabolome alterations during abiotic and biotic stress conditions. Therefore, metabolomics is a promising tool for knowledge-based targeted breeding programs. It also shows that due to missing databases and non-standardized LC-MS conditions, basic metabolomics, lipidomics and phytohormonics strategies, without isolation and unequivocal structure identification experiments, are sometimes insufficient. Therefore, techniques using biological and molecular structural characterizations of marker metabolites in combination with metabolomics techniques, such as phytometabolomics or sensomics approaches, are useful solutions to produce high-quality phytoalexin-enriched functional foods in the future.