Metabolomics for Biomarker Discovery: Key Signatory Metabolic Profiles for the Identification and Discrimination of Oat Cultivars

Omics-driven exploration and mining of key functional genes for the improvement of food and fiber crops

The deployment of omics technologies has obtained an incredible boost over the past few decades with the advances in next-generation sequencing (NGS) technologies, innovative bioinformatics tools, and the deluge of available biological information. The major omics technologies in the limelight are genomics, transcriptomics, proteomics, metabolomics, and phenomics. These biotechnological advances have modernized crop breeding and opened new horizons for developing crop varieties with improved traits. The genomes of several crop species are sequenced, and a huge number of genes associated with crucial economic traits have been identified. These identified genes not only provide insights into the understanding of regulatory mechanisms of crop traits but also decipher practical grounds to assist in the molecular breeding of crops. This review discusses the potential of omics technologies for the acquisition of biological information and mining of the genes associated with important agronomic traits in important food and fiber crops, such as wheat, rice, maize, potato, tomato, cassava, and cotton. Different functional genomics approaches for the validation of these important genes are also highlighted. Furthermore, a list of genes discovered by employing omics approaches is being represented as potential targets for genetic modifications by the latest genome engineering methods for the development of climate-resilient crops that would in turn provide great impetus to secure global food security.

Vetiver aerial parts and roots ameliorate rheumatoid arthritis in complete Freund's adjuvant rat model, a phytochemical profiling and mechanistic study

ETHNOPHARMACOLOGICAL RELEVANCE

Vetiver (Chrysopogon zizanioides) is indigenous to India where it is traditionally used to relief rheumatisms, lumbagos and sprains. Vetiver anti-inflammatory activity has not been previously investigated, and its specific interactions with body inflammation cascade remain largely unknown.

AIM OF THE STUDY

The present work was performed to validate the ethnobotanical use of the plant and compare the anti-inflammatory activities of the ethanolic extracts of the most traditionally used part (aerial part) to that of the root. Furthermore, we attempt to reveal the molecular mechanism of this anti-inflammatory activity in correlation to the chemical composition of C. zizanioides aerial (CA) and root parts (CR).

MATERIALS AND METHODS

Ultraperformance liquid chromatography coupled to high resolution mass spectrometry (UHPLC/HRMS) was used for comprehensive analysis of both CA and CR. The anti-inflammatory effect of both extracts was evaluated in complete Freund's adjuvant (CFA)-induced RA model in Wistar rats.

RESULTS

Phenolic metabolites were predominant in CA and 42 were identified for the first time, while only 13 were identified in CR. Meanwhile, triterpenes and sesquiterpenes were confined to the root extract. In CFA arthritis model, CA showed better anti-inflammatory activity than CR marked by an increase in serum level of IL-10 with simultaneous decrease in pro-inflammatory markers; IL-6, ACPA and TNF-α and was evident in histopathological examination. This anti-inflammatory effect was accompanied by down-regulation of JAK2/STAT3/SOCs3, ERK1/ERK2, TRAF6/c-FOS/NFATC1, TRAF6/NF-κB/NFATC1 and RANKL pathways which were all upregulated after CFA injection. These pathways were modulated to larger extent by CA, except for ERK1/ERK2 which was downregulated more effectively by CR. This differential effect between CA and CR can be explained by the variability in their phytoconstituents profile.

CONCLUSION

In agreement with the ethnobotanical preference, CA extract was more effective than CR extract in reducing the symptoms of RA probably due to its enrichment with flavonoids, lignans, and flavolignans. Both CA and CR reduced the production of inflammatory cytokines through modulating various biological signaling pathways. These findings support the traditional use of vetiver leaves as a remedy for RA and suggest that the use of the whole plant may offer advantage by synergistically affecting more inflammatory pathways.

Untargeted Metabolomics of Alternaria solani-Challenged Wild Tomato Species Solanum cheesmaniae Revealed Key Metabolite Biomarkers and Insight into Altered Metabolic Pathways

Untargeted metabolomics of moderately resistant wild tomato species Solanum cheesmaniae revealed an altered metabolite profile in plant leaves in response to Alternaria solani pathogen. Leaf metabolites were significantly differentiated in non-stressed versus stressed plants. The samples were discriminated not only by the presence/absence of specific metabolites as distinguished markers of infection, but also on the basis of their relative abundance as important concluding factors. Annotation of metabolite features using the Arabidopsis thaliana (KEGG) database revealed 3371 compounds with KEGG identifiers belonging to biosynthetic pathways including secondary metabolites, cofactors, steroids, brassinosteroids, terpernoids, and fatty acids. Annotation using the Solanum lycopersicum database in PLANTCYC PMN revealed significantly upregulated (541) and downregulated (485) features distributed in metabolite classes that appeared to play a crucial role in defense, infection prevention, signaling, plant growth, and plant homeostasis to survive under stress conditions. The orthogonal partial least squares discriminant analysis (OPLS-DA), comprising a significant fold change (≥2.0) with VIP score (≥1.0), showed 34 upregulated biomarker metabolites including 5-phosphoribosylamine, kaur-16-en-18-oic acid, pantothenate, and O-acetyl-L-homoserine, along with 41 downregulated biomarkers. Downregulated metabolite biomarkers were mapped with pathways specifically known for plant defense, suggesting their prominent role in pathogen resistance. These results hold promise for identifying key biomarker metabolites that contribute to disease resistive metabolic traits/biosynthetic routes. This approach can assist in mQTL development for the stress breeding program in tomato against pathogen interactions.

Avenanthramides, Distinctive Hydroxycinnamoyl Conjugates of Oat, Avena sativa L.: An Update on the Biosynthesis, Chemistry, and Bioactivities

Avenanthramides are a group of N-cinnamoylanthranilic acids (phenolic alkaloid compounds) that are produced in oat plants as phytoalexins, in response to pathogen attack and elicitation. The enzyme catalysing the cinnamamide-generating reaction is hydroxycinnamoyl-CoA: hydroxyanthranilate N-hydroxycinnamoyltransferase (HHT, a member of the super family of BAHD acyltransferases). HHT from oat appears to have a narrow range of substrate usage, with preferred use of 5-hydroxyanthranilic acid (and to a lesser extent, other hydroxylated and methoxylated derivatives) as acceptor molecules, but is able to use both substituted cinnamoyl-CoA and avenalumoyl-CoA thioesters as donor molecules. Avenanthramides thus combine carbon skeletons from both the stress-inducible shikimic acid and phenylpropanoid pathways. These features contribute to the chemical characteristics of avenanthramides as multifunctional plant defence compounds, as antimicrobial agents and anti-oxidants. Although avenanthramides are naturally and uniquely synthesised in oat plants, these molecules also exhibit medicinal and pharmaceutical uses important for human health, prompting research into utilisation of biotechnology to enhance agriculture and value-added production.

Small molecule metabolites: discovery of biomarkers and therapeutic targets

Metabolic abnormalities lead to the dysfunction of metabolic pathways and metabolite accumulation or deficiency which is well-recognized hallmarks of diseases. Metabolite signatures that have close proximity to subject's phenotypic informative dimension, are useful for predicting diagnosis and prognosis of diseases as well as monitoring treatments. The lack of early biomarkers could lead to poor diagnosis and serious outcomes. Therefore, noninvasive diagnosis and monitoring methods with high specificity and selectivity are desperately needed. Small molecule metabolites-based metabolomics has become a specialized tool for metabolic biomarker and pathway analysis, for revealing possible mechanisms of human various diseases and deciphering therapeutic potentials. It could help identify functional biomarkers related to phenotypic variation and delineate biochemical pathways changes as early indicators of pathological dysfunction and damage prior to disease development. Recently, scientists have established a large number of metabolic profiles to reveal the underlying mechanisms and metabolic networks for therapeutic target exploration in biomedicine. This review summarized the metabolic analysis on the potential value of small-molecule candidate metabolites as biomarkers with clinical events, which may lead to better diagnosis, prognosis, drug screening and treatment. We also discuss challenges that need to be addressed to fuel the next wave of breakthroughs.

Metabolic Variations among Three New Tea Varieties Cultivated in Shandong, China

Cultivar identification is a necessary step in tea breeding programs. Rapid identification methods would greatly improve these breeding processes. To preliminarily identify the three new Lucha tea varieties (LC6, LC7, and LC17) cultivated in Shandong, we measured their main agronomic characters and biochemical components. Then, we analyzed the metabolic profiles of these tea varieties and Fuding Dabaicha (FD) using a UPLC-ESI-MS/MS system. Their biochemical components indicated that the Lucha varieties had excellent varietal characteristics, with higher amino acid contents. Furthermore, secondary metabolism changed a lot in the Lucha tea varieties compared with that in the FD, with their accumulations of flavonoids and phenolic acids showing significant differences. These differential flavonoids were dominated by flavones and flavanone, flavonols, flavonoid carbonosides, and flavanols monomer. Flavanols especially, including epicatechin glucoside, epicatechin-3-(3″-O-methyl)gallate, epigallocatechin-3-O-(3,5-O-dimethyl)gallate, and epitheaflavic acid-3-O-Gallate, showed higher levels in the Lucha varieties. The phenolic acids containing caffeoyl groups showed higher levels in the Lucha varieties than those in the FD, while those containing galloyl groups showed a reverse pattern. Nitrogen metabolism, including amino acids, also showed obvious differences between the Lucha varieties and FD. The differential amino acids were mainly higher in the Lucha varieties, including 5-L-glutamyl-L-amino acid, N-monomethyl-L-arginine, and N-α-acetyl-L-ornithine. By using these approaches, we found that LC6, LC7, and LC17 were excellent varieties with a high yield and high quality for making green teas in Shandong.

A Metabolomic Study of the Analgesic Effect of Lappaconitine Hydrobromide (LAH) on Inflammatory Pain

Lappaconitine (LA) is a C-18 diterpene alkaloid isolated from Aconitum sinomontanum Nakai that has been shown to relieve mild to moderate discomfort. Various researchers have tried to explain the underlying mechanism of LA's effects on chronic pain. This article uses metabolomics technology to investigate the metabolite alterations in the dorsal root ganglion (DRG) when lappaconitine hydrobromide (LAH) was injected in an inflammatory pain model, to explain the molecular mechanism of its analgesia from a metabolomics perspective. The pain model used in this study was a complete Freund's adjuvant (CFA)-induced inflammatory pain model in rats. There were two treatment groups receiving different dosages of LAH (4 mg/kg LAH and 8 mg/kg LAH). The analgesic mechanism of LAH was investigated with an analgesic behavioral test, tissue sections, and metabolomics. The results of the analgesic behavioral experiment showed that both 4 mg/kg LAH and 8 mg/kg LAH could significantly improve the paw withdrawal latency (PWL) of rats. The tissue section results showed that LAH could reduce the inflammatory response and enlargement of the paw and ankle of rats and that there was no significant difference in the tissue sections of the DRG. The metabolomics results showed that retinol metabolism and glycerophospholipid metabolism in the CFA-induced inflammatory pain model were significantly affected and may exacerbate the inflammatory reactions and initiate persistent pain; in addition, the linoleic acid metabolism, arachidonic acid metabolism, and alanine, aspartate, and glutamate metabolism were also slightly affected. Among them, the alpha-linolenic acid metabolism was up-regulated after LAH treatment, while the retinol metabolism was down-regulated. These results suggest that LAH could effectively reduce inflammatory pain and might achieve this by regulating the lipid metabolism in the rat DRG.

Mass Spectrometry-Based Nontargeted and Targeted Analytical Approaches in Fingerprinting and Metabolomics of Food and Agricultural Research

Mass spectrometry (MS)-based techniques have been extensively applied in food and agricultural research. This review aims to address the advances and applications of MS-based analytical strategies in nontargeted and targeted analysis and summarizes the recent publications of MS-based techniques, including flow injection MS fingerprinting, chromatography-tandem MS metabolomics, direct analysis using ambient mass spectrometry, as well as development in MS data deconvolution software packages and databases for metabolomic studies. Various nontargeted and targeted approaches are employed in marker compounds identification, material adulteration detection, and the analysis of specific classes of secondary metabolites. In the newly emerged applications, the recent advances in computer tools for the fast deconvolution of MS data in targeted secondary metabolite analysis are highlighted.

Untargeted metabolite profiling to elucidate rhizosphere and leaf metabolome changes of wheat cultivars (Triticum aestivum L.) treated with the plant growth-promoting rhizobacteria Paenibacillus alvei (T22) and Bacillus subtilis

The rhizosphere is a highly complex and biochemically diverse environment that facilitates plant–microbe and microbe–microbe interactions, and this region is found between plant roots and the bulk soil. Several studies have reported plant root exudation and metabolite secretion by rhizosphere-inhabiting microbes, suggesting that these metabolites play a vital role in plant–microbe interactions. However, the biochemical constellation of the rhizosphere soil is yet to be fully elucidated and thus remains extremely elusive. In this regard, the effects of plant growth-promoting rhizobacteria (PGPR)–plant interactions on the rhizosphere chemistry and above ground tissues are not fully understood. The current study applies an untargeted metabolomics approach to profile the rhizosphere exo-metabolome of wheat cultivars generated from seed inoculated (bio-primed) with Paenibacillus (T22) and Bacillus subtilis strains and to elucidate the effects of PGPR treatment on the metabolism of above-ground tissues. Chemometrics and molecular networking tools were used to process, mine and interpret the acquired mass spectrometry (MS) data. Global metabolome profiling of the rhizosphere soil of PGPR-bio-primed plants revealed differential accumulation of compounds from several classes of metabolites including phenylpropanoids, organic acids, lipids, organoheterocyclic compounds, and benzenoids. Of these, some have been reported to function in plant–microbe interactions, chemotaxis, biocontrol, and plant growth promotion. Metabolic perturbations associated with the primary and secondary metabolism were observed from the profiled leaf tissue of PGPR-bio-primed plants, suggesting a distal metabolic reprograming induced by PGPR seed bio-priming. These observations gave insights into the hypothetical framework which suggests that PGPR seed bio-priming can induce metabolic changes in plants leading to induced systemic response for adaptation to biotic and abiotic stress. Thus, this study contributes knowledge to ongoing efforts to decipher the rhizosphere metabolome and mechanistic nature of biochemical plant–microbe interactions, which could lead to metabolome engineering strategies for improved plant growth, priming for defense and sustainable agriculture.

Metabolomics and Chemoinformatics in Agricultural Biotechnology Research: Complementary Probes in Unravelling New Metabolites for Crop Improvement

Simple Summary

The world is facing an overarching threat to food security, particularly in developing nations. The issue is further exacerbated by the apparent impacts of biotic and abiotic stresses driving down crop yields and productivity. Conventional strategies to improve yields and sustain productivity have been employed, including plant breeding for favourable and resilient agronomic traits. However, the efficacy and success rates of these methods are declining, partly due to the rapid changes in climate variability and the emergence of new and resistant phytopathogens. Additionally, the process of creating new and improved transgenic varieties of crops is long and can be expensive. Thus, new and innovative technologies are required for crop improvement. This review explores recent advances in the science of metabolomics and chemoinformatics, which have presented an avenue for rapid and robust analysis; moreover, it explores the elucidation of the complex plant metabolome, providing the opportunity to decipher the reactionary mechanisms of plants to the surrounding environment through their metabolic activity. As such, specific metabolites can, thus, be selected as biomarkers for crop improvement based on their functional characteristics under varying environmental conditions (growth, development, and defence). This new knowledge can enhance breeding practices through rapid and robust metabolic engineering techniques for sustainable agriculture.

Abstract

The United Nations (UN) estimate that the global population will reach 10 billion people by 2050. These projections have placed the agroeconomic industry under immense pressure to meet the growing demand for food and maintain global food security. However, factors associated with climate variability and the emergence of virulent plant pathogens and pests pose a considerable threat to meeting these demands. Advanced crop improvement strategies are required to circumvent the deleterious effects of biotic and abiotic stress and improve yields. Metabolomics is an emerging field in the omics pipeline and systems biology concerned with the quantitative and qualitative analysis of metabolites from a biological specimen under specified conditions. In the past few decades, metabolomics techniques have been extensively used to decipher and describe the metabolic networks associated with plant growth and development and the response and adaptation to biotic and abiotic stress. In recent years, metabolomics technologies, particularly plant metabolomics, have expanded to screening metabolic biomarkers for enhanced performance in yield and stress tolerance for metabolomics-assisted breeding. This review explores the recent advances in the application of metabolomics in agricultural biotechnology for biomarker discovery and the identification of new metabolites for crop improvement. We describe the basic plant metabolomics workflow, the essential analytical techniques, and the power of these combined analytical techniques with chemometrics and chemoinformatics tools. Furthermore, there are mentions of integrated omics systems for metabolomics-assisted breeding and of current applications.

Characterization of Ikaria Heather Honey by Untargeted Ultrahigh-Performance Liquid Chromatography-High Resolution Mass Spectrometry Metabolomics and Melissopalynological Analysis

Honey represents a valuable food commodity, known since ancient times for its delicate taste and health benefits due to its specific compositional characteristics, mainly the phenolic compound content. “Anama” honey is a monofloral honey produced from the nectar of Erica manipuliflora plant, a heather bush of the Greek island of Ikaria, one of the Mediterranean’s longevity regions. “Anama” is characterized by a unique aroma and taste, with a growing demand for consumption and the potential to be included in the list of products with a protected designation of origin. The aim of this study was to determine the chemical and botanical profile of authentic Anama honey samples and find similarities and differences with honey samples of a different botanical origin from the same geographical area. Untargeted Ultrahigh-Performance Liquid Chromatography-Hybrid Quadrupole-Orbitrap High-Resolution Mass Spectrometry (UHPLC-HRMS) metabolomics study was conducted on authentic heather, pine, and thyme honey samples from Ikaria and neighboring islands. The Principal Component Analysis (PCA), Orthogonal Projections to Latent Structures Discriminant Analysis (OPLS-DA), and differential analysis were performed using the entire metabolic profile of the samples and allowed the identification of chemical markers for sample discrimination. Thirty-two characteristic secondary metabolites (cinnamic acids, phenolic acids, flavonoids, terpenes) and other bioactive phenolic compounds, some of them not previously reported in a heather honey (aucubin, catalpol, domesticoside, leonuriside A, picein among others), emerged as potential chemical indicators of Anama honey. Melissopalynological analysis was also carried out to decipher the botanical and geographical origin of Anama honey. The relative frequency of the pollen of dominant plants of the Ericaceae family and a multitude of nectariferous and nectarless plants contributing to the botanical profile of Anama was evaluated. The identification of the pollen sources enabled a potential correlation of differentially increased secondary metabolites and chemicals with their botanical origin. The physicochemical profile of Anama was also determined, including the parameters of pH, color, electrical conductivity, diastase, moisture, as well as sugars, supporting the high quality of this heather honey.

UHPLC/MS-Based Serum Metabolomics Reveals the Mechanism of Radiation-Induced Thrombocytopenia in Mice

Radiation-induced thrombocytopenia is a common and life-threatening side effect of ionizing radiation (IR) therapy. However, the underlying pathological mechanisms remain unclear. In the present study, irradiation was demonstrated to significantly reduce platelet levels, inhibit megakaryocyte differentiation, and promote the apoptosis of bone marrow (BM) cells. A metabolomics approach and a UHPLC-QTOF MS system were subsequently employed for the comprehensive analysis of serum metabolic profiles of normal and irradiated mice. A total of 66 metabolites were significantly altered, of which 56 were up-regulated and 10 were down-regulated in irradiated mice compared to normal mice on day 11 after irradiation. Pathway analysis revealed that disorders in glycerophospholipid metabolism, nicotinate and nicotinamide metabolism, sphingolipid metabolism, inositol phosphate metabolism, and tryptophan metabolism were involved in radiation-induced thrombocytopenia. In addition, three important differential metabolites, namely L-tryptophan, LysoPC (17:0), and D-sphinganine, which were up-regulated in irradiated mice, significantly induced the apoptosis of K562 cells. L-tryptophan inhibited megakaryocyte differentiation of K562 cells. Finally, serum metabolomics was performed on day 30 (i.e., when the platelet levels in irradiated mice recovered to normal levels). The contents of L-tryptophan, LysoPC (17:0), and D-sphinganine in normal and irradiated mice did not significantly differ on day 30 after irradiation. In conclusion, radiation can cause metabolic disorders, which are highly correlated with the apoptosis of hematopoietic cells and inhibition of megakaryocyte differentiation, ultimately resulting in thrombocytopenia. Further, the metabolites, L-tryptophan, LysoPC (17:0), and D-sphinganine can serve as biomarkers for radiation-induced thrombocytopenia.

Metabolomic Characterisation of Discriminatory Metabolites Involved in Halo Blight Disease in Oat Cultivars Caused by Pseudomonas syringae pv. coronafaciens

The metabolome is the underlying biochemical layer of the phenotype and offers a functional readout of the cellular mechanisms involved in a biological system. Since metabolites are considered end-products of regulatory processes at a cellular level, their levels are considered the definitive response of the biological system to genetic or environmental variations. The metabolome thus serves as a metabolic fingerprint of the biochemical events that occur in a biological system under specific conditions. In this study, an untargeted metabolomics approach was applied to elucidate biochemical processes implicated in oat plant responses to Pseudomonas syringae pv. coronafaciens (Ps-c) infection, and to identify signatory markers related to defence responses and disease resistance against halo blight. Metabolic changes in two oat cultivars (“Dunnart” and “SWK001”) responding to Ps-c, were examined at the three-leaf growth stage and metabolome changes monitored over a four-day post-inoculation period. Hydromethanolic extracts were analysed using an ultra-high-performance liquid chromatography (UHPLC) system coupled to a high-definition mass spectrometer (MS) analytical platform. The acquired multi-dimensional data were processed using multivariate statistical analysis and chemometric modelling. The validated chemometric models indicated time- and cultivar-related metabolic changes, defining the host response to the bacterial inoculation. Further multivariate analyses of the data were performed to profile differential signatory markers, putatively associated with the type of launched defence response. These included amino acids, phenolics, phenolic amides, fatty acids, flavonoids, alkaloids, terpenoids, lipids, saponins and plant hormones. Based on the results, metabolic alterations involved in oat defence responses to Ps-c were elucidated and key signatory metabolic markers defining the defence metabolome were identified. The study thus contributes toward a more holistic understanding of the oat metabolism under biotic stress.

Comparative Metabolite Profiling of Wheat Cultivars (Triticum aestivum) Reveals Signatory Markers for Resistance and Susceptibility to Stripe Rust and Aluminium (Al3+) Toxicity

Plants continuously produce essential metabolites that regulate their growth and development. The enrichment of specific metabolites determines plant interactions with the immediate environment, and some metabolites become critical in defence responses against biotic and abiotic stresses. Here, an untargeted UHPLC-qTOF-MS approach was employed to profile metabolites of wheat cultivars resistant or susceptible to the pathogen Puccinia striiformis f. sp. tritici (Pst) and Aluminium (Al³⁺) toxicity. Multivariate statistical analysis (MVDA) tools, viz. principal component analysis (PCA) and hierarchical cluster analysis (HiCA) were used to qualify the correlation between the identified metabolites and the designated traits. A total of 100 metabolites were identified from primary and secondary metabolisms, including phenolic compounds, such as flavonoid glycosides and hydroxycinnamic acid (HCA) derivatives, fatty acids, amino acids, and organic acids. All metabolites were significantly variable among the five wheat cultivars. The Pst susceptible cultivars demonstrated elevated concentrations of HCAs compared to their resistant counterparts. In contrast, ‘Koonap’ displayed higher levels of flavonoid glycosides, which could point to its resistant phenotype to Pst and Al³⁺ toxicity. The data provides an insight into the metabolomic profiles and thus the genetic background of Pst- and Al³⁺-resistant and susceptible wheat varieties. This study demonstrates the prospects of applied metabolomics for chemotaxonomic classification, phenotyping, and potential use in plant breeding and crop improvement.

A Metabolomics Approach and Chemometric Tools for Differentiation of Barley Cultivars and Biomarker Discovery

One of the ultimate goals of plant breeding is the development of new crop cultivars capable of withstanding increasing environmental stresses, to sustain the constantly growing population and economic demands. Investigating the chemical composition of the above and underground tissues of cultivars is crucial for the understanding of common and specific traits thereof. Using an untargeted metabolomics approach together with appropriate chemometrics tools, the differential metabolite profiles of leaf and root extracts from five cultivars of barley (‘Erica’, ‘Elim’, ‘Hessekwa’, ‘S16’ and ‘Agulhas’) were explored and potential signatory biomarkers were revealed. The study was conducted on seedlings grown for 21 days under identical controlled conditions. An ultra-high performance liquid chromatography coupled to high-resolution mass spectrometry (UHPLC-HRMS) was employed to analyse hydromethanolic leaf and root extracts of barley cultivars. Furthermore, unsupervised and supervised learning algorithms were applied to mine the generated data and to pinpoint cultivar-specific metabolites. Among all the classes of metabolites annotated, phenolic acids and derivatives formed the largest group and also represented the most discriminatory metabolites. In roots, saponarin, an important allelochemical differentially distributed across cultivars, was the only flavonoid annotated. The application of an untargeted metabolomics approach in phenotyping grain crops such as barley was demonstrated, and the metabolites responsible for differentiating between the selected cultivars were revealed. The study provides insights into the chemical architecture of barley, an agro-economically relevant cereal crop; and reiterates the importance of metabolomics tools in plant breeding practices for crop improvement.