Effect of Genotype, Year, and Their Interaction on the Accumulation of Bioactive Compounds and the Antioxidant Activity in Industrial Hemp (Cannabis sativa L.) Inflorescences.

The phytochemical content and the antioxidant activity in the inflorescences of six industrial hemp (Cannabis sativa L.) genotypes, four monoecious (Codimono, Carmaleonte, Futura 75, and Santhica 27), and two dioecious (Fibrante and Carmagnola Selezionata), were assessed for three consecutive years from 2018 to 2020. The total phenolic content, total flavonoid content, and antioxidant activity were determined by spectrophotometric measurements, whereas HPLC and GC/MS were used to identify and quantify the phenolic compounds, terpenes, cannabinoids, tocopherols, and phytosterols. All the measured traits were significantly affected by genotype (G), cropping year (Y), and their interaction (G × Y), although the Y effect prevailed as a source of variation, ranging from 50.1% to 88.5% for all the metabolites except cannabinoids, which were equally affected by G, Y, and G × Y interaction (33.9%, 36.5%, and 21.4%, respectively). The dioecious genotypes presented a more constant performance over the three years compared to the monoecious genotypes, with the highest and most stable phytochemical content observed in the inflorescences of Fibrante, which was characterized by the highest levels of cannabidiol, α-humulene and ß-caryophyllene, which may confer on the inflorescences of this genotype a great economic value due to the important pharmacological properties of these metabolites. Conversely, the inflorescences of Santhica 27 were characterized by the lowest accumulation of phytochemicals over the cropping years, with the notable exception of cannabigerol, a cannabinoid that exhibits a wide range of biological activities, which was found at its highest level in this genotype. Overall, these findings can be used by breeders in future programs aimed at the selection of new hemp genotypes with improved levels of phytochemicals in their inflorescences, which can provide better health and industrial benefits.

Enrichment of provitamin A content in durum wheat grain by suppressing ß-carotene hydroxylase 1 genes with a TILLING approach.

KEY MESSAGE: The suppression of the HYD-1 gene by a TILLING approach increases the amount of ß-carotene in durum wheat kernel. Vitamin A deficiency is a major public health problem that affects numerous countries in the world. As humans are not able to synthesize vitamin A, it must be daily assimilated along with other micro- and macronutrients through the diet. Durum wheat is an important crop for Mediterranean countries and provides a discrete amount of nutrients, such as carbohydrates and proteins, but it is deficient in some essential micronutrients, including provitamin A. In the present work, a targeting induced local lesions in genomes strategy has been undertaken to obtain durum wheat genotypes biofortified in provitamin A. In detail, we focused on the suppression of the ß-carotene hydroxylase 1 (HYD1) genes, encoding enzymes involved in the redirection of ß-carotene toward the synthesis of the downstream xanthophylls (neoxanthin, violaxanthin and zeaxanthin). Expression analysis of genes involved in carotenoid biosynthesis revealed a reduction of the abundance of HYD1 transcripts greater than 50% in mutant grain compared to the control. The biochemical profiling of carotenoid in the wheat mutant genotypes highlighted a significant increase of more than 70% of ß-carotene compared to the wild-type sibling lines, with no change in lutein, α-carotene and zeaxanthin content. This study sheds new light on the molecular mechanism governing carotenoid biosynthesis in durum wheat and provides new genotypes that represent a good genetic resource for future breeding programs focused on the provitamin A biofortification through non-transgenic approaches.

Comparative Analysis Based on Transcriptomics and Metabolomics Data Reveal Differences between Emmer and Durum Wheat in Response to Nitrogen Starvation.

Mounting evidence indicates the key role of nitrogen (N) on diverse processes in plant, including development and defense. Using a combined transcriptomics and metabolomics approach, we studied the response of seedlings to N starvation of two different tetraploid wheat genotypes from the two main domesticated subspecies: emmer and durum wheat. We found that durum wheat exhibits broader and stronger response in comparison to emmer as seen from the expression pattern of both genes and metabolites and gene enrichment analysis. They showed major differences in the responses to N starvation for transcription factor families, emmer showed differential reduction in the levels of primary metabolites while durum wheat exhibited increased levels of most of them to N starvation. The correlation-based networks, including the differentially expressed genes and metabolites, revealed tighter regulation of metabolism in durum wheat in comparison to emmer. We also found that glutamate and γ-aminobutyric acid (GABA) had highest values of centrality in the metabolic correlation network, suggesting their critical role in the genotype-specific response to N starvation of emmer and durum wheat, respectively. Moreover, this finding indicates that there might be contrasting strategies associated to GABA and glutamate signaling modulating shoot vs. root growth in the two different wheat subspecies.

Provitamin A Biofortification of Durum Wheat through a TILLING Approach.

Macro- and micronutrients, essential for the maintenance of human metabolism, are assimilated daily through the diet. Wheat and other major cereals are a good source of nutrients, such as carbohydrates and proteins, but cannot supply a sufficient amount of essential micronutrients, including provitamin A. As vitamin A deficiency (VAD) leads to several serious diseases throughout the world, the biofortification of a major staple crop, such as wheat, represents an effective way to preserve human health in developing countries. In the present work, a key enzyme involved in the branch of carotenoids pathway producing ß-carotene, lycopene epsilon cyclase, has been targeted by a Targeting Induced Local Lesions in Genomes (TILLING) approach in a "block strategy" perspective. The null mutant genotype showed a strong reduction in the expression of the lcyE gene and also interesting pleiotropic effects on an enzyme (ß-ring hydroxylase) acting downstream in the pathway. Biochemical profiling of carotenoids in the wheat mutant lines showed an increase of roughly 75% in ß-carotene in the grains of the complete mutant line compared with the control. In conclusion, we describe here the production and characterization of a new wheat line biofortified with provitamin A obtained through a nontransgenic approach, which also sheds new light on the molecular mechanism governing carotenoid biosynthesis in durum wheat.

Evolutionary Metabolomics Reveals Domestication-Associated Changes in Tetraploid Wheat Kernels.

Domestication and breeding have influenced the genetic structure of plant populations due to selection for adaptation from natural habitats to agro-ecosystems. Here, we investigate the effects of selection on the contents of 51 primary kernel metabolites and their relationships in three Triticum turgidum L. subspecies (i.e., wild emmer, emmer, durum wheat) that represent the major steps of tetraploid wheat domestication. We present a methodological pipeline to identify the signature of selection for molecular phenotypic traits (e.g., metabolites and transcripts). Following the approach, we show that a reduction in unsaturated fatty acids was associated with selection during domestication of emmer (primary domestication). We also show that changes in the amino acid content due to selection mark the domestication of durum wheat (secondary domestication). These effects were found to be partially independent of the associations that unsaturated fatty acids and amino acids have with other domestication-related kernel traits. Changes in contents of metabolites were also highlighted by alterations in the metabolic correlation networks, indicating wide metabolic restructuring due to domestication. Finally, evidence is provided that wild and exotic germplasm can have a relevant role for improvement of wheat quality and nutritional traits.

Metabolomic analysis can detect the composition of pasta enriched with fibre after cooking.

BACKGROUND: Several studies have demonstrated that metabolomics has a definite place in food quality, nutritional value, and safety issues. The aim of the present study was to determine and compare the metabolites in different pasta samples with fibre, and to investigate the modifications induced in these different kinds of pasta during cooking, using a gas chromatography-mass spectrometry-based metabolomics approach. RESULTS: Differences were seen for some of the amino acids, which were absent in control pasta, while were present both in the commercially available high-fibre pasta (samples A-C) and the enriched pasta (samples D-F). The highest content in reducing sugars was observed in enriched samples in comparison with high-fibre pasta. The presence of stigmasterol in samples enriched with wheat bran was relevant. Cooking decreased all of the metabolites: the high-fibre pasta (A-C) and Control showed losses of amino acids and tocopherols, while for sugars and organic acids, the decrease depended on the pasta sample. The enriched pasta samples (D-F) showed the same decreases with the exception of phytosterols, and in pasta with barley the decrease of saturated fatty acids was not significant as for tocopherols in pasta with oat. Principal component analysis of the metabolites and the pasta discrimination was effective in differentiating the enriched pasta from the commercial pasta, both uncooked and cooked. CONCLUSIONS: The study has established that such metabolomic analyses provide useful tools in the evaluation of the changes in nutritional compounds in high-fibre and enriched pasta, both before and after cooking. © 2015 Society of Chemical Industry.

Impact of domestication on the phenotypic architecture of durum wheat under contrasting nitrogen fertilization.

The process of domestication has led to dramatic morphological and physiological changes in crop species due to adaptation to cultivation and to the needs of farmers. To investigate the phenotypic architecture of shoot- and root-related traits and quantify the impact of primary and secondary domestication, we examined a collection of 36 wheat genotypes under optimal and nitrogen-starvation conditions. These represented three taxa that correspond to key steps in the recent evolution of tetraploid wheat (i.e. wild emmer, emmer, and durum wheat). Overall, nitrogen starvation reduced the shoot growth of all genotypes, while it induced the opposite effect on root traits, quantified using the automated phenotyping platform GROWSCREEN-Rhizo. We observed an overall increase in all of the shoot and root growth traits from wild emmer to durum wheat, while emmer was generally very similar to wild emmer but intermediate between these two subspecies. While the differences in phenotypic diversity due to the effects of primary domestication were not significant, the secondary domestication transition from emmer to durum wheat was marked by a large and significant decrease in the coefficient of additive genetic variation. In particular, this reduction was very strong under the optimal condition and less intense under nitrogen starvation. Moreover, although under the optimal condition both root and shoot traits showed significantly reduced diversity due to secondary domestication, under nitrogen starvation the reduced diversity was significant only for shoot traits. Overall, a considerable amount of phenotypic variation was observed in wild emmer and emmer, which could be exploited for the development of pre-breeding strategies.

Metabolic profiling of a mapping population exposes new insights in the regulation of seed metabolism and seed, fruit, and plant relations.

To investigate the regulation of seed metabolism and to estimate the degree of metabolic natural variability, metabolite profiling and network analysis were applied to a collection of 76 different homozygous tomato introgression lines (ILs) grown in the field in two consecutive harvest seasons. Factorial ANOVA confirmed the presence of 30 metabolite quantitative trait loci (mQTL). Amino acid contents displayed a high degree of variability across the population, with similar patterns across the two seasons, while sugars exhibited significant seasonal fluctuations. Upon integration of data for tomato pericarp metabolite profiling, factorial ANOVA identified the main factor for metabolic polymorphism to be the genotypic background rather than the environment or the tissue. Analysis of the coefficient of variance indicated greater phenotypic plasticity in the ILs than in the M82 tomato cultivar. Broad-sense estimate of heritability suggested that the mode of inheritance of metabolite traits in the seed differed from that in the fruit. Correlation-based metabolic network analysis comparing metabolite data for the seed with that for the pericarp showed that the seed network displayed tighter interdependence of metabolic processes than the fruit. Amino acids in the seed metabolic network were shown to play a central hub-like role in the topology of the network, maintaining high interactions with other metabolite categories, i.e., sugars and organic acids. Network analysis identified six exceptionally highly co-regulated amino acids, Gly, Ser, Thr, Ile, Val, and Pro. The strong interdependence of this group was confirmed by the mQTL mapping. Taken together these results (i) reflect the extensive redundancy of the regulation underlying seed metabolism, (ii) demonstrate the tight co-ordination of seed metabolism with respect to fruit metabolism, and (iii) emphasize the centrality of the amino acid module in the seed metabolic network. Finally, the study highlights the added value of integrating metabolic network analysis with mQTL mapping.

Effects of Heat Stress on Metabolite Accumulation and Composition, and Nutritional Properties of Durum Wheat Grain.

Durum wheat (Triticum turgidum (L.) subsp. turgidum (L.) convar. durum (Desf.)) is momentous for human nutrition, and environmental stresses can strongly limit the expression of yield potential and affect the qualitative characteristics of the grain. The aim of this study was to determine how heat stress (five days at 37 °C) applied five days after flowering affects the nutritional composition, antioxidant capacity and metabolic profile of the grain of two durum wheat genotypes: "Primadur", an elite cultivar with high yellow index, and "T1303", an anthocyanin-rich purple cultivar. Qualitative traits and metabolite evaluation (by gas chromatography linked to mass spectrometry) were carried out on immature (14 days after flowering) and mature seeds. The effects of heat stress were genotype-dependent. Although some metabolites (e.g., sucrose, glycerol) increased in response to heat stress in both genotypes, clear differences were observed. Following the heat stress, there was a general increase in most of the analyzed metabolites in "Primadur", with a general decrease in "T1303". Heat shock applied early during seed development produced changes that were observed in immature seeds and also long-term effects that changed the qualitative and quantitative parameters of the mature grain. Therefore, short heat-stress treatments can affect the nutritional value of grain of different genotypes of durum wheat in different ways.

Impact of Chitosan-Based Foliar Application on the Phytochemical Content and the Antioxidant Activity in Hemp (Cannabis sativa L.) Inflorescences.

In the present study, the phytochemical content and the antioxidant activity in the inflorescences of the monoecious hemp cultivar Codimono grown in southern Italy were assessed, and their elicitation was induced by foliar spray application of 50 mg/L and 250 mg/L of chitosan (CHT) at three different molecular weights (low, CHT L; medium, CHT M; high CHT H). The analysis of the phytochemical profile confirmed that cannabinoids were the most abundant class (54.2%), followed by flavonoids (40.3%), tocopherols (2.2%), phenolic acids (1.9%), and carotenoids (1.4%). Cannabinoids were represented almost exclusively by cannabidiol, whereas cannabigerol and Δ9-tetrahydrocannabinol were detected at very low levels (the latter was below the legal limit of 0.3%). The most abundant flavonoids were orientin and vitexin, whereas tocopherols were mainly represented by α-tocopherol. The antioxidant activity was found to be positively correlated with flavonoids and tocopherols. Statistical analysis revealed that the CHT treatments significantly affected the phytochemical content and the antioxidant activity of hemp inflorescences. Notably, a significant increase in the total phenolic content (from +36% to +69%), the α-tocopherol (from +45% to +75%) and ß+γ-tocopherol (from +35% to +82%) contents, and the ABTS radical scavenging activity (from +12% to +28%) was induced by all the CHT treatments. In addition, treatments with CHT 50 solutions induced an increase in the total flavonoid content (from +12% to +27%), as well as in the vitexin (from +17% to +20%) and orientin (from +20% to +30%) contents. Treatment with CHT 50 L almost always resulted in the greatest increases. Overall, our findings indicated that CHT could be used as a low-cost and environmentally safe elicitor to improve the health benefits and the economic value of hemp inflorescences, thus promoting their employment in the food, pharmaceutical, nutraceutical, and cosmetic supply chains.

Can the use of magnetized water affect the seedling development and the metabolite profiles of two different species: Lentil and durum wheat?

Seedlings of durum wheat and lentil were utilized to investigate the efficiency of magnetic water on growth and metabolic epicotyl profile. Tap water was passed through a magnetic device with a flow rate of max. 12900 - 13200 Gauss (G). Seeds and plantlets were grown on sand-free paper soaked by magnetized water, with unmagnetized tap water used in a control group. The growth parameters were collected at three time points (48, 96, and 144 hours after treatment), the same times at which metabolomic analysis was conducted on seeds, roots, and epicotyls. Although the effects varied with the species, tissues, and time point considered, compared with tap water (TW), the use of magnetized water treatment (MWT) led to higher root elongation in both genotypes. On the contrary, epicotyl length was not affected by treatment both in durum wheat and lentil. The results indicate that the use of magnetized water in agriculture can be considered a sustainable technology to promote plant development and quality with reduced and more efficient water usage, leading to cost-saving and environmental protection.

Chitosan-Induced Activation of the Antioxidant Defense System Counteracts the Adverse Effects of Salinity in Durum Wheat.

The present study was carried out with the aim of (i) evaluating the effect of chitosan (CTS) on the growth of durum wheat under salinity and (ii) examining CTS-regulated mechanisms of salinity tolerance associated with the antioxidant defense system. To achieve these goals, durum wheat seedlings were treated with CTS at different molecular weight, low (L-CTS, 50-190 kDa), medium (M-CTS, 190-310 kDa) and high (H-CTS, 310-375 kDa). The results obtained show that exposure to 200 mM NaCl reduced the shoot and the root dried biomass by 38% and 59%, respectively. The growth impairment induced by salinity was strongly correlated with an increase in the superoxide anion production (5-fold), hydrogen peroxide content (2-fold) and malondialdehyde (MDA) content (4-fold). Seedlings responded to the oxidative stress triggered by salinity with an increase in the total phenolic content (TPC), total flavonoid content (TFC) and total antioxidant activity (TAA) by 67%, 51% and 32%, respectively. A salt-induced increase in the activity of the antioxidant enzymes superoxide dismutase and catalase (CAT) of 89% and 86%, respectively, was also observed. Treatment of salt-stressed seedlings with exogenous CTS significantly promoted seedling growth, with the strongest effects observed for L-CTS and M-CTS, which increased the shoot biomass of stressed seedlings by 32% and 44%, respectively, whereas the root dried biomass increased by 87% and 64%, respectively. L-CTS and M-CTS treatments also decreased the superoxide anion production (57% and 59%, respectively), the hydrogen peroxide content (35% and 38%, respectively) and the MDA content (48% and 56%, respectively) and increased the TPC (23% and 14%, respectively), the TFC (19% and 10%, respectively), the TAA (up to 10% and 7%, respectively) and the CAT activity (29% and 20%, respectively). Overall, our findings indicate that CTS exerts its protective role against the oxidative damages induced by salinity by enhancing the antioxidant defense system. L-CTS and M-CTS were the most effective in alleviating the adverse effect of NaCl, thus demonstrating that the CTS action is strictly related to its molecular weight.

Domestication of Crop Metabolomes: Desired and Unintended Consequences.

The majority of the crops and vegetables of today were domesticated from their wild progenitors within the past 12 000 years. Considerable research effort has been expended on characterizing the genes undergoing positive and negative selection during the processes of crop domestication and improvement. Many studies have also documented how the contents of a handful of metabolites have been altered during human selection, but we are only beginning to unravel the true extent of the metabolic consequences of breeding. We highlight how crop metabolomes have been wittingly or unwittingly shaped by the processes of domestication, and highlight how we can identify new targets for metabolite engineering for the purpose of de novo domestication of crop wild relatives.

Characterization of Nutritional Quality Traits of a Common Bean Germplasm Collection.

Food legumes are at the crossroads of many societal challenges that involve agriculture, such as climate change and food sustainability and security. In this context, pulses have a crucial role in the development of plant-based diets, as they represent a very good source of nutritional components and improve soil fertility, such as by nitrogen fixation through symbiosis with rhizobia. The main contribution to promotion of food legumes in agroecosystems will come from plant breeding, which is guaranteed by the availability of well-characterized genetic resources. Here, we analyze seeds of 25 American and European common bean purified accessions (i.e., lines of single seed descent) for different morphological and compositional quality traits. Significant differences among the accessions and superior genotypes for important nutritional traits are identified, with some lines showing extreme values for more than one trait. Heritability estimates indicate the importance of considering the effects of environmental growth conditions on seed compositional traits. They suggest the need for more phenotypic characterization in different environments over different years to better characterize combined effects of environment and genotype on nutritional trait variations. Finally, adaptation following the introduction and spread of common bean in Europe seems to have affected its nutritional profile. This finding further suggests the relevance of evolutionary studies to guide breeders in the choice of plant genetic resources.

Relationship between Seed Morphological Traits and Ash and Mineral Distribution along the Kernel Using Debranning in Durum Wheats from Different Geographic Sites.

Debranning was applied to durum wheat to the study the relationship between kernel shape and size, and ash and mineral distribution having implications for semolina yield. To this aim four durum wheat genotypes carried out over three environments were selected to determine the morphological and yield traits as well as the distribution along the kernel of the ash, macro- (Na, K, P, Ca, and Mg), and micro-elements (Mn, Fe, Cu, Zn, and Mo). A descendent ash gradient within the kernel reflects the decreases in the minerals that occurred during debranning. Perciasacchi with high seed weight (TKW) and greater thickness followed by Cappelli showed a more uniform distribution of ash content along the kernels. High r Pearson coefficient (p < 0.01) showed an inverse relationship between thickness and ash decay. Since thickness was strongly affected by the genotype, it could represent a useful trait for breeding programs to predict the milling quality.

Metabolomics Provides Valuable Insight for the Study of Durum Wheat: A Review.

Metabolomics is increasingly being applied in various fields offering a highly informative tool for high-throughput diagnostics. However, in plant sciences, metabolomics is underused, even though plant studies are relatively easy and cheap when compared to those on humans and animals. Despite their importance for human nutrition, cereals, and especially wheat, remain understudied from a metabolomics point of view. The metabolomics of durum wheat has been essentially neglected, although its genetic structure allows the inference of common mechanisms that can be extended to other wheat and cereal species. This review covers the present achievements in durum wheat metabolomics highlighting the connections with the metabolomics of other cereal species (especially bread wheat). We discuss the metabolomics data from various studies and their relationships to other "-omics" sciences, in terms of wheat genetics, abiotic and biotic stresses, beneficial microbes, and the characterization and use of durum wheat as feed, food, and food ingredient.

Analysis of metabolic and mineral changes in response to salt stress in durum wheat (Triticum turgidum ssp. durum) genotypes, which differ in salinity tolerance.

The key mechanisms of salinity tolerance (ST) in durum wheat were investigated, with five genotypes used to determine changes in morpho-physiological traits and mineral and metabolite contents after exposure to 50, 100 and 200 mM NaCl. Plant growth impairment was evident at the highest salt level. Under this condition, a wide range of shoot Na+ contents and ST were observed within genotypes. However, no significant correlation was seen between ST and Na+ exclusion from the shoots, which indicates that tissue tolerance also has a role. Consistent with this, there was significant correlation between ST and the Na+:K+ ratio in the shoots. Indeed, the maintenance of the shoot Na+ and K+ homeostasis was found to be essential to achieve osmotic adjustment, which relied substantially on inorganic osmolytes, and to avoid toxicity symptoms, such as chlorophyll loss, which appeared only at the highest salinity level. Consistently, the metabolite changes occurred mainly in the shoots, with a dual response to salinity: (i) a conserved response that was common to all the genotypes and resulted in the accumulation of proline and in the depletion of organic acids, including some intermediates of the Krebs cycle; and (ii) a genotype-specific response that involved the accumulation of GABA, threonine, leucine, glutamic acid, glycine, mannose and fructose and appeared related to the different tolerance of genotypes to salinity. The lower magnitude of response to salinity detected in the roots confirmed the major role of the shoots in the determination of ST of durum wheat.

Mineral composition of durum wheat grain and pasta under increasing atmospheric CO2 concentrations.

The concentrations of 10 minerals were investigated in the grain of 12 durum wheat genotypes grown under free air CO2 enrichment conditions, and in four of their derived pasta samples, using inductively coupled plasma mass spectrometry. Compared to ambient CO2 (400ppm; AMB), under elevated CO2 (570ppm; ELE), the micro-element and macro-element contents showed strong and significant decreases in the grain: Mn, -28.3%; Fe, -26.7%; Zn, -21.9%; Mg, -22.7%; Mo, -40.4%; K, -22.4%; and Ca, -19.5%. These variations defined the 12 genotypes as sensitive or non-sensitive to ELE. The pasta samples under AMB and ELE showed decreased mineral contents compared to the grain. Nevertheless, the contributions of the pasta to the recommended daily allowances remained relevant, also for the micro-elements under ELE conditions (range, from 18% of the recommended daily allowance for Zn, to 70% for Mn and Mo).

Stoichiometric Correlation Analysis: Principles of Metabolic Functionality from Metabolomics Data.

Recent advances in metabolomics technologies have resulted in high-quality (time-resolved) metabolic profiles with an increasing coverage of metabolic pathways. These data profiles represent read-outs from often non-linear dynamics of metabolic networks. Yet, metabolic profiles have largely been explored with regression-based approaches that only capture linear relationships, rendering it difficult to determine the extent to which the data reflect the underlying reaction rates and their couplings. Here we propose an approach termed Stoichiometric Correlation Analysis (SCA) based on correlation between positive linear combinations of log-transformed metabolic profiles. The log-transformation is due to the evidence that metabolic networks can be modeled by mass action law and kinetics derived from it. Unlike the existing approaches which establish a relation between pairs of metabolites, SCA facilitates the discovery of higher-order dependence between more than two metabolites. By using a paradigmatic model of the tricarboxylic acid cycle we show that the higher-order dependence reflects the coupling of concentration of reactant complexes, capturing the subtle difference between the employed enzyme kinetics. Using time-resolved metabolic profiles from Arabidopsis thaliana and Escherichia coli, we show that SCA can be used to quantify the difference in coupling of reactant complexes, and hence, reaction rates, underlying the stringent response in these model organisms. By using SCA with data from natural variation of wild and domesticated wheat and tomato accession, we demonstrate that the domestication is accompanied by loss of such couplings, in these species. Therefore, application of SCA to metabolomics data from natural variation in wild and domesticated populations provides a mechanistic way to understanding domestication and its relation to metabolic networks.