Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 10 October 2022
Sec. Crop and Product Physiology
This article is part of the Research Topic Improving crop nutritional security for sustainable agriculture in the era of climate change View all 12 articles

Single and interactive effects of variables associated with climate change on wheat metabolome

Kristýna Ve
e&#x;ovKristýna Večeřová1Michal OravecMichal Oravec1Swati PuranikSwati Puranik1Hana Findurov,Hana Findurová1,2Barbora VeselBarbora Veselá1Emmanuel Opoku,Emmanuel Opoku1,2Kojo Kwakye Ofori-Amanfo,Kojo Kwakye Ofori-Amanfo1,3Karel KlemKarel Klem1Otmar UrbanOtmar Urban1Pranav Pankaj Sahu*Pranav Pankaj Sahu1*
  • 1Laboratory of Ecological Plant Physiology, Global Change Research Institute of the Czech Academy of Sciences, Brno, Czechia
  • 2Department of Agrosystems and Bioclimatology, Faculty of AgriSciences, Mendel University in Brno, Brno, Czechia
  • 3Department of Forest Ecology, Faculty of Forestry and Wood Technology, Mendel University in Brno, Brno, Czechia

One of the key challenges linked with future food and nutritional security is to evaluate the interactive effect of climate variables on plants’ growth, fitness, and yield parameters. These interactions may lead to unique shifts in the morphological, physiological, gene expression, or metabolite accumulation patterns, leading to an adaptation response that is specific to future climate scenarios. To understand such changes, we exposed spring wheat to 7 regimes (3 single and 4 combined climate treatments) composed of elevated temperature, the enhanced concentration of CO2, and progressive drought stress corresponding to the predicted climate of the year 2100. The physiological and metabolic responses were then compared with the current climate represented by the year 2020. We found that the elevated CO2 (eC) mitigated some of the effects of elevated temperature (eT) on physiological performance and metabolism. The metabolite profiling of leaves revealed 44 key metabolites, including saccharides, amino acids, and phenolics, accumulating contrastingly under individual regimes. These metabolites belong to the central metabolic pathways that are essential for cellular energy, production of biosynthetic pathways precursors, and oxidative balance. The interaction of eC alleviated the negative effect of eT possibly by maintaining the rate of carbon fixation and accumulation of key metabolites and intermediates linked with the Krebs cycle and synthesis of phenolics. Our study for the first time revealed the influence of a specific climate factor on the accumulation of metabolic compounds in wheat. The current work could assist in the understanding and development of climate resilient wheat by utilizing the identified metabolites as breeding targets for food and nutritional security.

Introduction

Plants are sensitive to changes in temperature, water availability, and atmospheric CO2 concentration (Thornton et al., 2014). As reported by the Intergovernmental Panel on Climate Change (IPCC), extreme climate events like heat waves and drought periods will become more intense and frequent, which could have devastating consequences on plant growth, development, and production (IPCC, 2018). The combinations of these climate factors may have even a more severe impact on food availability to the increasing world population. The strategy to improve the resilience of economically important crops has mostly relied on the screening of genotypes by exposing them to a single climate factor (e.g., drought, temperature, CO2, etc.). However, the adaptive response under simultaneous exposure to multiple climate factors may vary widely, resulting in an undesired or failed crop selection process. Therefore, to prevent such crop failure in future, an assessment of the impact of combined treatments on plant’s performance and adaptability is required.

Plants respond to stress by adjusting their morpho-physiological status such as by inhibiting photosynthesis (Jagtap et al., 1998), and/or altering vegetative growth and biomass allocation (Pandey et al., 2015). The molecular mechanisms of acclimation constitute the activation of genes for adaptive responses, such as the heat shock proteins (HSPs), and/or enzymes and additional transcripts facilitating detoxification and signaling of reactive oxygen species (ROS) (Rizhsky et al., 2004). In addition, under severe conditions, membrane, and protein damage lead to the accumulation of lipid peroxidation products (Jiang and Huang, 2001). This induces altered redox homeostasis and increase in ROS and consequently enhanced cellular oxidative damage (Krasensky and Jonak, 2012; Munné-Bosch et al., 2012). In addition, a wide range of different metabolites such as primary metabolites (carbohydrates, tricarboxylic acid cycle intermediates, and amino acids) essential for the growth and development, secondary metabolites (phenolics, terpenes, and nitrogen-containing compounds) along with the hormones (such as abscisic acid, jasmonic acid, salicylic acid, and ethylene) are also produced to regulate plant climate response (Fang et al., 2019; Kumar et al., 2020; Kumari et al., 2020). The primary metabolites, such as amino acids, play a pivotal role in enhancing fitness under various environmental cues. For example, the accumulation of amino acid proline is considered an indicator of adaptive behavior of plants exposed to various abiotic factors (Hayat et al., 2012; Fàbregas and Fernie, 2019; Girija et al., 2022). It helps in the maintenance of cell turgidity and reduction in electrolyte leakage to prevent cellular oxidative bursts (Hayat et al., 2012). Besides being the major metabolic resource and structural constituent of cells, carbohydrates such as sucrose and hexoses, also regulate stress-signaling and control the expression of growth as well as stress-related genes (Rosa et al., 2009). The TCA cycle intermediates are also involved in providing environmental stress tolerance attributes (Tahjib-Ul-Arif et al., 2021) and their accumulation under drought and heat stress has been reported to assist in the maintenance of physiological performance by regulating the photosynthetic rates and cellular redox status (Fernie et al., 2004; Das et al., 2017). On the other hand, secondary metabolites mitigate the plant’s survival by protecting it from adverse environments (Agostini-Costa et al., 2012; Kurepin et al., 2017). For example, the antioxidants nature of flavonoids (a major class of phenylpropanoids) regulates various functions such as maintaining water homeostasis (Naing and Kim, 2021), and stomatal dynamics (Li et al., 2021) to improve the crop stress adaptability.

Identification of changes in the composition of such key compounds and derivate through the metabolomic approach serves as a powerful strategy to advance our understanding of the adaptation responses in plants (Obata and Fernie, 2012; Sardans et al., 2020; Girija et al., 2022). Various studies have provided novel insights into the metabolic signature among crop species (e.g., rice, maize, and wheat) wide-ranging across growth stages and genotypes (Peng et al., 2017; Xu et al., 2019; Chen et al., 2020). However, investigating the effects of climate change is complex because multiple interactions between various environmental factors could prompt unique shifts in the biochemical composition. A recent study in Arabidopsis revealed the difference in the metabolite accumulations under single and multi-climate factors (Zinta et al., 2018). While elevated CO2 was found to largely enhance the cellular concentration of saccharides and amino acids, the combination with other factors caused the opposite effect on the accumulation of these compounds. Very few studies have been done to identify the shifts in metabolite pattern under combined heat and drought stress, particularly in crops (Obata et al., 2015; Templer et al., 2017). These studies revealed that most of the metabolic changes are additive, i.e., represent the sum of responses to the individual factors, and drought contributed most to such changes.

Bread wheat (Triticum aestivum) is a major cereal crop and widely cultivated for its grains (Shewry, 2009). The climate variability risks the livelihood and food security of about ~2 billion people who rely on wheat globally. The decline in the climate resilience of most of the wheat cultivars in the last 15 years was already observed in Europe, due to the repeated selection for few desirable traits (Kahiluoto et al., 2019). As most traits are complex and multigenic in nature, impact assessment of environmental stress on wheat crop relies not only on the type of factor exposed but also on its phenological growth stage. For example, the water availability at tillering is a critical phase in terms of grain yield potential (Curtis et al., 2002). Depending upon the severity and duration of drought stress, flowering and grain-filling phases are also considered as key stages which may result in the substantial yield loss (Farooq et al., 2014). The flowering stage is more sensitive to temperature stress which regulates the seed-setting and grain quality (Aiqing et al., 2018). The sensitivity of wheat to any of these climate variables can result in altered metabolic processes coupled with lower biomass accumulation and grain yield (Ihsan et al., 2016). Despite the advancement in knowledge, the impact of climate factor combinations on metabolite accumulation in wheat remains ambiguous. Therefore, considering its economic importance, we aimed to examine its acclimation responses to current as well as future climate scenarios by using an experimental model of single and combined climate treatments. This work allowed us to identify the most important compounds contributing to metabolome variations in wheat leaves under multifactorial treatments and their potential to the risks associated with climate change. The key metabolites and pathways identified in this study could be the potential targets to understand their future climate-induced regulation, and development of nutritionally enriched crops.

Materials and methods

Experimental regimes and plant growth set-up

The experiment started with setting climate conditions representing the typical growth environment of Central Europe for spring wheat. Three major climate variables representing key factors associated with an ongoing climate change, i.e., air temperature, water availability and atmospheric CO2 concentration were used. The data for air temperature (daily maximum, minimum and average) representing the years 2020 (current) and 2100 (future) were retrieved from the World Bank portal (https://climateknowledgeportal.worldbank.org/country/czech-republic/climate-data-projections). For the control set-up, an hourly dynamic of temperature (ambient temperature; aT; maximum/minimum= 24/10°C), constant CO2 concentration (ambient CO2; aC; 400 ± 50 ppm) and well-watered soil (water holding capacity; WHC = 90 ± 5%) conditions were maintained to mimic the environment for the year 2020. To understand the influence of individual climate variables and their combinations on acclimation responses in the year 2100, seven experimental setups were used. These included single factor treatments comprising of elevated temperature (eT), elevated CO2 (eC) and drought (D), and combinations of two-factors (eT+D, eC+D, eC+eT) and three-factors (eC+eT+D). For these setups, the chambers were enabled with the daily air temperature of maximum/minimum= 28/12°C (eT) and CO2 concentration of 700 ppm (eC). The applied 700 ppm CO2 concentration as an elevated condition was based on the projections for year 2100 (Meehl et al., 2007) as well as previous studies (Reyenga et al., 2001; Jauregui et al., 2015; Soba et al., 2019). Further, a progressive drought was applied by with-holding the water supply after 7 days of optimization until WHC reached 25 ± 5% and kept constant thereafter throughout the experiment. As each setup had different evaporation rates, a minimum of 10-days continuous drought was applied depending upon the onset of 25% WHC. The pattern of daily changes in the relative humidity (RH; maximum day/night, 65/90%), maximum day light intensity (800 µmol/m2/s) and light duration (15/9 h) were kept as constant for all experimental set-ups. The selection of light intensity (800 µmol/m2/s) was based on previous studies conducted on wheat to study the impact of single or multiple climate factors in a growth cabinet conditions (Nuttall et al., 2018; Chavan et al., 2019).

The spring wheat cultivar Cadenza was selected for our study as it the most widely used research variety (Fernández-Gómez et al., 2020) and genomic resources are also available for further study (Walkowiak et al., 2020). Seeds after stratification (4°C/48 h) and germination (24°C/48 h) were planted in the standard substrate (TS2; Klasmann-Deilmann, Geeste, Germany). Pots (dimension: 11×11×25 cm) containing the seedlings were later transferred to climate chambers (FytoScope FS-SI 3400; PSI, Drásov, Czech Republic) for respective treatments. Seedlings were initially allowed to optimize for seven days under temperature (daily maximum/minimum= 15/10°C), constant CO2 concentration (400 ± 50 ppm), RH (daily minimum/maximum, 65/90%), maximum/minimum light intensity (800/0 µmol/m2/s), light duration (15/9 h) and water saturated conditions. Afterwards, they were subjected to the control and treatment/experimental set-ups (detailed above) with each setup containing at least three plants per treatment (n=3). All physiological measurements and sampling for metabolite profiling were done on the youngest fully developed leaves (second from the top) of the primary shoot at end of the tillering stage (Feekes-6) as it is an agronomically important developmental stage that contributes to grain productivity (Khadka et al., 2020). All samplings and measurements were conducted between 11.00 and 14.00 hours Central European Time (CET), when the temperature and light intensity reached maxima and relative air humidity reached a minimum for the day.

Physiological measurements

The measurements of photosynthesis and photosynthesis-related processes were done using Li-6800 gas-exchange system (LI-COR Biosciences, Lincoln, Nebraska, USA). To determine light-saturated CO2 assimilation rate (A), stomatal conductance (Gsw), transpiration rate (Tr), and water use efficiency (WUE = A/Tr) the leaves investigated were exposed to the light intensity of 1200 µmol/m2/s, while other microclimatic conditions inside the assimilation chamber were kept constant corresponding to daily maxima of individual treatments: CO2 concentration (400 ppm or 700 ppm for aC and eC regimes, respectively), leaf temperature (24°C or 28°C) and VPD (1.04 kPa and 1.32 kPa for aT and eT regimes, respectively). The measurements were conducted on three biological replicates. All data after testing for normal distribution were subject to one- or two-way analysis of variance (ANOVA) followed by post-hoc comparison using Tukey’s multiple range test (P ≤ 0.05). Principal component analysis (PCA) was used to analyze relations between variables. All statistical analyses were performed using the STATISTICA package version 12 (Stat Soft, Tulsa, Oklahoma, USA).

Metabolite profiling of wheat leaf at vegetative stage

The metabolomic analysis was performed using the established protocol for the metabolic studies in cereals in the laboratory of metabolomics and isotopic analyses at the CzechGlobe (Klem et al., 2019). Briefly, detached leaves at 25-days post-treatments were immediately stored in liquid nitrogen. Lyophilized samples (100 mg dry weight) were subsequently homogenized and extracted using a methanol: chloroform: H2O solution (1:2:2). Aliquots of the polar phase were used to analyze saccharides, phenolic compounds, amino acids, phytohormones, and Krebs cycle intermediates employing an UltiMate 3000 system of high-performance liquid chromatography (HPLC) coupled with an LTQ Orbitrap XL high-resolution mass spectrometer (HRMS) (ThermoFisher Scientific, Waltham, Massachusetts, USA). The HRMS is equipped with a heated electrospray ionization source and was operated in full scan mode with a resolution of 60000. For the separation of individual compounds, a Hypersil GOLD chromatographic column (ThermoFisher Scientific, Waltham, Massachusetts, USA) of 150 mm (length), 2.1 mm (ID), and 3 μm (film thickness) was used. Moreover, gas chromatography coupled with mass spectrometry (GC-MS) was employed to analyze a spectrum of fatty acids, sterols, saccharides, and other non-polar compounds. Analyses were performed with a TSQ Quantum XLS triple Quadrupole (ThermoFisher Scientific, Waltham, Massachusetts, USA) on a Rxi-5Sil MS capillary column (Restek, Bellefonte, Pennsylvania, USA) of 30 m (length) with 5 m of integra-Guard, 0.25 mm (ID) and 0.25 µm (film thickness). The compounds detected by the chromatography were identified based on the comparison with the own mass library which had been created from the measurement of standards using HPLC-HRMS and GC-MS in full scan mode and containing about 300 compounds. The statistical analysis such as partial least squares discriminant analysis (PLS-DA), mean decrease accuracy (MDA), one-way ANOVA of the identified metabolites was done using MetaboAnalyst software (Pang et al., 2021). Correlation analyses were performed using the STATISTICA package version 12 (Stat Soft, Tulsa, Oklahoma, USA).

Results

Setting the regimes to study the impact of climate

As predicted by RCP 8.5 scenarios of IPCC, temperature, drought, and atmospheric CO2 concentration were used as the three major factors to study the impact of climate change (Figure 1A). Current (2020) and expected future (2100) growing conditions in Central Europe corresponding to June were simulated as described in the materials and methods. The growth of plants was monitored regularly until 25 day-post treatment (DPT). The tiller and leaf number were counted after 25 DPT, however, there was no significant difference at this stage (Figure 1B).

FIGURE 1
www.frontiersin.org

Figure 1 A schematic representation of experimental pipeline and developmental changes monitored in wheat cv. Cadenza. (A) Climate setup. Experimental setup to mimic the natural daily temperature dynamics (Tmax/Tmin) in June of years 2020 and 2100, atmospheric CO2 concentration [ambient (400 ppm) and elevated (700 ppm)] and water availability [water holding capacity (WHC) at 25 ± 5% (Drought, D) and 90 ± 5% (water saturated pots)]. The symbol of ± indicates fluctuations in the parameters. (B) Phenotypic differences observed at 25 DPT showing the variable effects of treatments. The bar scale represents 5 cm. Ctr, control treatment; eC, elevated CO2; eT, elevated temperature; D, drought; eC+eT, elevated CO2 and elevated temperature; eC+D, elevated CO2 and drought; eT+D, elevated temperature and drought; eC+eT+D, elevated CO2, elevated temperature, and drought.

Concentration of CO2 determines the physiological responses

It was revealed that the pattern of change in the physiological parameters was mostly treatment-specific (Figures 2A–E). The variance (PC1 ~86%) was clearly observable due to the effect of eT (Figure 2A). Transpiration rate (Tr) and stomatal conductance (Gsw) significantly declined (P ≤ 0.05) under eC, eC+eT and eC+D treatments in comparison to the control condition (Figures 2B, C). There was no clear trend in the pattern of the photosynthesis rate (A) except for an increase under eC and eC+eT+D conditions (Figure 2D). On the contrary, eC+eT treatment led to a significant reduction of A values. Noticeably, eC alone or in combinations with other factors significantly (P ≤ 0.05) increased the water use efficiency (WUE), although this effect was less pronounced under eC+eT+D condition. These findings indicate a critical role of CO2 concentration in the adjusting of physiological state under future conditions.

FIGURE 2
www.frontiersin.org

Figure 2 Impact of single and multiple climate variables on the physiology of wheat cv. Cadenza in tillering stage. (A) Principal component analysis plot of the physiological traits measured at 25 DPT. An account of (B) Transpiration rate, Tr; (C) Stomatal conductance, Gsw; (D) Rate of photosynthesis, A; (E) Water use efficiency, WUE, under single and combined climate variable treatments. Three biological replicates (n=3) per treatment were used to evaluate the physiological traits. The box plot shows the inter-quartile range with the mean. Values marked with the same letter do not differ according to Tukey’s test multiple range tests (P ≤ 0.05). Ctr, control treatment; eC, elevated CO2; eT, elevated temperature; D, drought; eC+eT, elevated CO2 and elevated temperature; eC+D, elevated CO2, and drought; eT+D, elevated temperature, and drought; eC+eT+D, elevated CO2, elevated temperature, and drought.

Metabolites accumulate differentially under various regimes

Overall, 77 metabolites were identified by GC-MS and HPLC-HRMS analysis among the leaf samples exposed to different climate regimes (Supplementary Table 1). Out of these, 44 metabolites significantly (P ≤ 0.05) varied among the treatments (Supplementary Table 2). The most influenced metabolite groups included amino acids (16), Krebs cycle intermediates (8), phenolics (10), saccharides (4), fatty acids (3) and nitrogenous bases (3) (Supplementary Table 2). A distinct separation of metabolites between the regimes with aT and those with eT was observed, accounting for ~20% of the variation (PC1; Figure 3A). This finding indicates that temperature has a strong effect on the metabolome. The main reasons for the metabolite differences in different regimes were identified by Partial least square discriminant analysis (PLS-DA). Based on their VIP scores, which is the measure of a variable’s importance, the metabolites such as tyrosin (Tyr), ferulic acid (Feru), pyruvic acid (Pyru), 3-hydroxybenzoic acid (3hyd), epigallocatechin gallate (Epgg), guanine (G), uracil (U), and thymidine monophosphate (TMP) showed largest differences among all the regimes (Figure 3B). Metabolites having VIP score >1.5 were considered significant (Figure 3C). The accumulation patterns of individual metabolite categories under the subjected regimes are described in detail below. Overall, the results show a treatment-specific accumulation of metabolites during the vegetative growth stage in spring wheat.

FIGURE 3
www.frontiersin.org

Figure 3 A supervised Partial Least Squares Discriminant Analysis (PLS-DA) of all the metabolites accumulated in the leaves of wheat cv. Cadenza after exposure to climate regimes. (A) Scores plot showing the distribution of metabolites based on the regimes. PC1 separated the metabolites based on the ambient (aT, gray circle) and elevated temperature (eT, orange circle). (B) Loading plot showing the metabolites discriminated based on PLS-DA. The metabolites such as phenolics (green), amino acids (dark blue), Krebs cycle intermediates (red), lipids (yellow) and nitrogenous bases (brown) have been shown. (C) The metabolites ranked by variable importance in projection (VIP). Ctr, control treatment; eC, elevated CO2; eT, elevated temperature; D, drought; eC+eT, elevated CO2 and elevated temperature; eC+ D, elevated CO2 and drought; eT+D, elevated temperature and drought; eC+eT+D, elevated CO2, elevated temperature and drought. The color code represents the specific metabolite families: dark blue, amino acid; light blue, saccharides; green, phenolics; yellow, fatty acids; red, Kreb’s cycle acids; brown, nucleotides. Compounds in grey represent non-significant (P ≤ 0.05) metabolites among the treatments identified by one-way ANOVA. Three biological replicates (n=3) per treatment were used in the study.

Out of 19 amino acids detected in the metabolite profiling, 16 amino acids were significantly (P ≤ 0.05) accumulated among the treatments (Supplementary Table 2). These included branched chain amino acids (BCAAs) such as Leucine (Leu), aromatic (Tryptophan; Try, Phenylalanine; Phe, and Tyrosine; Tyr), hydrophilic-basic (Lysine; Lys, Arginine; Arg and Histidine; His) and hydrophobic (Glycine; Gly, Proline; Pro, Hydroxy-proline; Hpro, Methionine; Met) amino acids (Supplementary Table 2). The first component (46.8%) separated the eT and eC regimes from other treatments (Figure 4A). Further, PLS-DA analysis discriminated Tyr and hydrophilic-basic amino acids from the group of all amino acids (Figure 4B). The VIP score (>1.5; P ≤ 0.05) suggested Tyr and Lys as the key amino acids that tend to accumulate under eT and eT+D conditions in comparison to Ctr and D regimes, respectively (Figure 4C). Proline and BCAAs which accumulated primarily under eC, eC+D and eC+eT+D regimes have been reported to have osmoregulatory roles (Girija et al., 2022). On the other hand, eC alleviated the effect of eT on the accumulation of amino acids, as the accumulation of Tyr and Lys decreased under eC+eT and eC+eT+D regimes when compared to aC counterparts (Figure 4C).

FIGURE 4
www.frontiersin.org

Figure 4 A supervised Partial Least Squares Discriminant Analysis (PLS-DA) of key metabolite categories accumulated in the leaves of wheat cv. Cadenza after exposure to climate regimes. (A) Scores plot showing the distribution of metabolites (amino acids; phenolics, Krebs cycle intermediates; saccharides) based on the regimes. Component 1 separated the metabolites based on the ambient temperature (aT, gray circle), elevated temperature (eT, orange circle). Ambient CO2 (aC, purple circle), and elevated CO2 concentrations (eC, blue circle). (B) Loading plot showing the metabolites (amino acids; phenolics, Krebs cycle intermediates; saccharides) discriminated based on PLS-DA. The metabolites in the bigger font significantly (P ≤ 0.05) varied among the treatments. (C) The metabolites ranked by variable importance in projection (VIP). Scale bar represents high and low VIP scores corresponding to the metabolites. Ctr, control treatment; eC, elevated CO2; eT, elevated temperature; D, drought; eC+eT, elevated CO2 and elevated temperature; eC+ D, elevated CO2 and drought; eT+D, elevated temperature and drought; eC+eT+D, elevated CO2, elevated temperature and drought. Abbreviations of the metabolites are provided in Supplementary Table 1.

The metabolite profiling identified 10 phenolic compounds significantly (P ≤ 0.05) accumulated among the treatments (Supplementary Table 2). The scores plot between the first and second components identified two distinct groups related to the aT and eT regimes (Figure 4A). Based on the PLS-DA and VIP scores, 3hyd, Feru, and Epgg were identified as the key metabolites varying under treatment conditions (Figure 4B). The accumulation of these phenolics was temperature-dependent and unlike amino acids, CO2 exposure did not affect their abundance (Figure 4C). Compared with eT and eC, the effect of D on metabolome was minor or negligible.

One-way ANOVA revealed that all the Krebs cycle-related intermediates were significantly (P ≤ 0.05) accumulated among the climate variable treatments (Supplementary Table 2). Similarly, to the phenolic compounds, the scores plot between the first and second components identified two distinct groups associated with the aT and eT regimes (Figure 4A). Acids, such as α-ketoglutaric (Ket), Pyru and malic acid (Mal), were discriminated from other Krebs cycle intermediates based on the PLS-DA and VIP score (Figure 4B). An eT treatment increased contents of Pyru and Ket, while the content of Mal decreased (Figure 4C).

Fructose (Fruc), glucose (Gluc), sucrose (Suc), and galactinol (Galc) were the most significantly (P ≤ 0.05) accumulated saccharides among the climate variable treatments (Supplementary Table 2). The pattern of saccharide abundance showed two distinct groups associated with the aC and eC regimes (Figure 4A). Based on the PLS-DA and VIP score (> 1.0), Fruc and Gluc were identified as the most important saccharides whose abundance was mainly associated with the concentration of CO2 (Figures 4B, C). Apart from the above metabolic compounds, nitrogenous bases such as uracil (U), guanine (G) and thymidine monophosphate (TMP) were found to be significantly (P ≤ 0.05) accumulated (Supplementary Table 2). The most accumulated fatty acids at different regimes were identified to be linoleic acid (Lino), oleic acid (Olei), and palmitoleic acid (Palm) (Supplementary Table 2).

The correlation between the 44 significant metabolites and physiological traits was also evaluated. The analysis identified a significant negative correlation between the physiological traits (Gsw and Tr) and the accumulation of Fruc, Gluc, Arg, and shikimic acid (Shi) (Supplementary Table 3). Contrastingly, WUE showed positive correlation with these compounds.

Discussion

Changing environment causes profound physiological and metabolic changes in the plants (Gianoli and Molina-Montenegro, 2021; Zandalinas et al., 2022). While the recent years have witnessed substantial advances in understanding the acclimation responses of wheat genotypes to various environmental drivers, most of the studies so far have mainly focused on single- or two-factor treatments, with limited possibility to assess the impacts of their mutual interactions (Zandalinas et al., 2018). Hence, the question of how plants will acclimate to a comprehensive future change in growth conditions is still far from complete. Therefore, understanding the underlying molecular mechanisms of adaptation response is crucial not only for the basic understanding of wheat stress response but also for applied wheat breeding. Through quantification of key metabolites, we reveal a differential metabolomic response against eight climate regimes (single and combined treatment) comprised of increased temperature, drought, and elevated CO2 in spring wheat.

The general mechanism to cope with changes in the environment involves an adjustment in the physiological and molecular traits (Pandey et al., 2017; Sharma et al., 2021; Alamri et al., 2022; Sharma et al., 2022). The shift in physiological responses can be attributed to specific abiotic factors and their nature i.e., individual and combined (Zandalinas and Mittler, 2022). Gas exchange measurements suggested that at the single-factor level, the effect of eC was predominant and led to the significant changes in Tr, Gsw and WUE, while eT and D had only minor effects. Similar patterns of physiological traits were also observed under double-factor regimes combining elevated CO2 with eT or D. As no prominent effect of eT and D on Tr and Gsw under aC was observed, this indicated that the response of these traits was mainly controlled by CO2 enrichment either at higher temperatures or decreased soil water availability. However, the physiological pattern exhibited in single and double factor combinations inclusive of an eC regime (eC+D and eC+eT) were not replicated in the triple factor regime (eC+eT+D). This suggests that the effect of eC could have been compromised by the simultaneous counteraction of the other two factors (eT and D). Higher A and Tr were observed during eC+eT+D treatment. In general, eC and eT have opposing effects on Tr (Kirschbaum and McMillan, 2018). An increase in the temperature enhances the Tr, however, due to partial stomatal closure under high CO2 concentration environment, a reduction in the Tr could be expected (Kirschbaum and McMillan, 2018). The increase in Tr under eC+eT+D in our study could be due to the influence of D, which could have helped to curtail the damage to the photosynthesis process (indicated by high A) through an evaporative cooling mechanism (Šigut et al., 2015).

The next focus was on investigating the metabolite accumulation to understand their role in plant acclimation under various regimes. The metabolic change patterns were highly treatment specific. For example, eC dominantly controlled the accumulation of soluble monosaccharides (Gluc and Fruc). A concomitant increase in the photosynthetic activity (A) was also observed in eC and the triple factor combinations eC+eT+D (Figure 2D). As monosaccharides are critical molecules for the primary carbon metabolism and the first substrates for starch biosynthesis, these results indicate that the accumulation of Gluc corresponded to the supplementary carbon fixation in these samples. The effect of higher CO2 concentration on the accumulation of compounds related to the sugar metabolites has also been reported in Arabidopsis thaliana (Ainsworth and Rogers, 2007; Zinta et al., 2018), and crops like wheat (Levine et al., 2008) and rice (Ainsworth, 2008). Moreover, saccharides are well known to regulate the rate of transpiration, particularly via stomatal dynamics (Li et al., 2018). The significant correlation identified in this study between the Gluc and Fruc accumulation with Tr and WUE of wheat leaves under various regimes could be attributed to the stomatal function.

On the other hand, metabolites associated with the Krebs cycle and phenolics were found to be mostly dependent on the exposure of temperature (Figure 4). High-temperature stress has been found to affect the levels of pyruvic acid, fumarate, malate, and citrate in poplar and soybean (Sicher, 2013; Ren et al., 2019). The spring wheat may have an inherent capacity to withstand the projected temperature increase (4°C) in the future climate by limiting the conversion of soluble sugars to organic acids of the Krebs cycle and secondary metabolites during elevated temperatures. The increase of carbon flux under eC conditions may then be utilized for synthesizing the compounds needed for acclimation responses and defense mechanisms as has been reported previously (Kiani-Pouya et al., 2017). Moreover, at the molecular level, the expression of genes linked with saccharides, Krebs cycle intermediates, and phenolics changes dynamically upon exposure to climate change-related factors (Krasensky and Jonak, 2012; Glaubitz et al., 2015; Xalxo et al., 2020). This may consequently facilitate modifications in the levels of many saccharides and phenolic compounds imparting specific acclimation responses. Accumulation of Feru (important for cell wall rigidity and strength), 3hyd (potential antioxidant), and Epgg (flavonoid) was dominantly controlled by the high temperature (Figure 4), suggesting that these compounds may help to maintain the cellular integrity and functionality during interaction with the environment.

Amino acids are an important molecular form of organic nitrogen in plants acting as intermediates for many metabolic reactions (Trovato et al., 2021). The pattern of amino acid accumulation exhibited a unique trend where eT played a pivotal role in the accumulation of amino acids like Tyr and Lys, but their enhanced accumulation pattern diminished when combined with eCO2 (Figure 4). CO2 enrichment has been previously shown to reverse the effects of elevated growth temperatures on many soluble amino acids in soybean leaves (Sicher, 2013). Noticeably, the level of amino acid and Krebs cycle intermediates showed contrasting patterns in eT and eC regimes, suggesting that a subtle balance may exist between these pathways.

In general, phenolics play a pivotal role in plants as antioxidants to scavenge the ROS upon exposure to abiotic stresses (Rice-Evans et al., 1997). The increase in their cellular concentration suggests the activation of redox homeostasis to protect against the possible damage to proteins and lipids (Awasthi et al., 2015). Similar to our study, an increase in the production of phenolics, and particularly flavonoids, under heat stress has been reported previously (Wu et al., 2016). Overall, these results show that the differential accumulation of metabolites regulates plastic response under future climate scenarios.

Our study represents a comprehensive analysis of metabolite homeostasis under different climate scenarios. This study highlighted that the metabolomic acclimation of spring wheat to a single and combination of climate factors is mainly associated with temperature and CO2 concentration. An environment with elevated CO2 mitigated some of the effects of elevated temperature on physiological performance and metabolism possibly by maintaining the rate of carbon fixation and accumulation of key metabolites and intermediates linked with sugar metabolism, Krebs cycle, and synthesis of phenolics. Usually, drought stress is manifested by substantial changes in plant metabolome including changes in simple saccharides that have a role of osmolytes. However, in this study, the effect of D on physiology and metabolome was minor/negligible. It is possible that the drought exposure applied to the targetted stage (Feekes-6) was likely mild, or for a short duration. In conclusion, the key metabolites and pathways identified in this study and their interactions with climate variables will be helpful in understanding their future climate-induced regulation and the manipulations of potential targets for the development of climate resilient and nutritionally enriched wheat crops (Irfan and Datta, 2017). Further, a better understanding of these strategies can provide an impetus to research on gene function discovery and biochemical evolution, which is foundational for improved metabolic engineering.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

PPS and KK conceptualized the research. KV, SP, MO, HF, BV, EO, KKO-A, and PPS were involved in the methodology and analysis. SP, KV, PPS, OU and KK wrote and edited the manuscript. All authors have read and agreed to the final version of the manuscript.

Funding

This research was funded by Czech Science Foundation, grant no. 20-25845Y (awarded to PPS). A part of the research was also supported by the project SustES “Adaptation strategies for sustainable ecosystem services and food security under adverse environmental conditions” (CZ.02.1.01/0.0/0.0/16_019/0000797) to SP, BV, OU and KK. Infrastructure facilities used in the experiments were supported by the Ministry of Education, Youth and Sports of CR within the CzeCOS program, grant number LM2018123HF.

Acknowledgments

We express our gratitude to the Prof. RNDr. Ing. Michal V. Marek, Director of the Global Change Research Institute CAS, Brno for providing access to the necessary facilities and support.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.1002561/full#supplementary-material

References

Agostini-Costa, S., da, T., F., R., R., H., Silveira, D. (2012). “Secondary metabolites,” in Chromatography and its applications (London, UK: InTech). XXXA., M. doi: 10.5772/35705

CrossRef Full Text | Google Scholar

Ainsworth, E. A. (2008). Rice production in a changing climate: A meta-analysis of responses to elevated carbon dioxide and elevated ozone concentration. Global Change Biol. 14, 1642–1650. doi: 10.1111/j.1365-2486.2008.01594.x

CrossRef Full Text | Google Scholar

Ainsworth, E. A., Rogers, A. (2007). The response of photosynthesis and stomatal conductance to rising [CO2]: Mechanisms and environmental interactions. Plant Cell Environ. 30, 258–270. doi: 10.1111/j.1365-3040.2007.01641.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Aiqing, S., Somayanda, I., Sebastian, S. V., Singh, K., Gill, K., Prasad, P. V. V., et al. (2018). Heat stress during flowering affects time of day of flowering, seed set, and grain quality in spring wheat. Crop Sci. 58, 380–392. doi: 10.2135/cropsci2017.04.0221

CrossRef Full Text | Google Scholar

Alamri, S., Siddiqui, M. H., Mukherjee, S., Kumar, R., Kalaji, H. M., Irfan, M., et al. (2022). Molybdenum-induced endogenous nitric oxide (NO) signaling coordinately enhances resilience through chlorophyll metabolism, osmolyte accumulation and antioxidant system in arsenate stressed-wheat (Triticum aestivum l.) seedlings, environmental pollution, Vol. 292. 269–7491. doi: 10.1016/j.envpol.2021.118268

CrossRef Full Text | Google Scholar

Awasthi, R., Bhandari, K., Nayyar, H. (2015). Temperature stress and redox homeostasis in agricultural crops. Front. Environ. Sci. 3. doi: 10.3389/fenvs.2015.00011

CrossRef Full Text | Google Scholar

Chavan, S. G., Duursma, R. A., Tausz, M., Ghannoum, O. (2019). Elevated CO2 alleviates the negative impact of heat stress on wheat physiology but not on grain yield. J. Exp. Bot. 70, 6447–6459. doi: 10.1093/jxb/erz386

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J., Hu, X., Shi, T., Yin, H., Sun, D., Hao, Y., et al. (2020). Metabolite-based genome-wide association study enables dissection of the flavonoid decoration pathway of wheat kernels. Plant Biotechnol. J. 18, 1722–1735. doi: 10.1111/pbi.13335

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, A., Rushton, P., Rohila, J. (2017). Metabolomic profiling of soybeans (Glycine max L) reveals the importance of sugar and nitrogen metabolism under drought and heat stress. Plants. 6, 21. doi: 10.3390/plants6020021

CrossRef Full Text | Google Scholar

Fàbregas, N., Fernie, A. R. (2019). The metabolic response to drought. J. Exp. Bot. 70, 1077–1085. doi: 10.1093/jxb/ery437

PubMed Abstract | CrossRef Full Text | Google Scholar

Fang, C., Fernie, A. R., Luo, J. (2019). Exploring the diversity of plant metabolism. Trends Plant Sci. 24, 83–98. doi: 10.1016/j.tplants.2018.09.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Curtis, B. C., Rajaram, S. (2002). Bread wheat: Improvement and production,”. (Rome: Food and Agriculture Organization of the United Nations), 554.

Google Scholar

Farooq, M., Hussain, M., Siddique, K. H. M. (2014). Drought stress in wheat during flowering and grain-filling periods. Crit. Rev. Plant Sci. 33, 331–349. doi: 10.1080/07352689.2014.875291

CrossRef Full Text | Google Scholar

Fernández-Gómez, J., Talle, B., Tidy, A. C., Wilson, Z. A. (2020). Accurate staging of reproduction development in cadenza wheat by non-destructive spike analysis. J. Exp. Bot. 71, 3475–3484. doi: 10.1093/jxb/eraa156

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernie, A. R., Carrari, F., Sweetlove, L. J. (2004). Respiratory metabolism: glycolysis, the TCA cycle and mitochondrial electron transport. Curr. Opin. Plant Biol. 7, 254–261. doi: 10.1016/j.pbi.2004.03.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Gianoli, E., Molina-Montenegro, M. A. (2021). Evolution of physiological performance in invasive plants under climate change. Evolution 75, 3181–3190. doi: 10.1111/evo.14314

PubMed Abstract | CrossRef Full Text | Google Scholar

Girija, A., Han, J., Corke, F., Brook, J., Doonan, J., Yadav, R., et al. (2022). Elucidating drought responsive networks in tef (Eragrostis tef) using phenomic and metabolomic approaches. Physiologia Plantarum. 174, e13597. doi: 10.1111/ppl.13597

PubMed Abstract | CrossRef Full Text | Google Scholar

Glaubitz, U., Erban, A., Kopka, J., Hincha, D.K., Zuther, E.. (2015). High night temperature strongly impacts TCA cycle, amino acid and polyamine biosynthetic pathways in rice in a sensitivity-dependent manner. Journal of Experimental Botany. 66, 6385–6397. doi: 10.1093/jxb/erv352

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayat, S., Hayat, Q., Alyemeni, M. N., Wani, A. S., Pichtel, J., Ahmad, A. (2012). Role of proline under changing environments. Plant Signaling Behav. 7, 1456–1466. doi: 10.4161/psb.21949

CrossRef Full Text | Google Scholar

Ihsan, M. Z., El-Nakhlawy, F. S., Ismail, S. M., Fahad, S., Daur, I. (2016). Wheat phenological development and growth studies as affected by drought and late season high temperature stress under arid environment. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.00795

CrossRef Full Text | Google Scholar

IPCC (2018). Summary for policymakers- global warming of 1.5°C. An IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty. Masson-Delmotte, V., Zhai, P., Pörtner, H.-O., Roberts, D., Skea, J., Shukla, P.R., et al (Geneva, Switzerland: World Meteorological Organization)

Google Scholar

Irfan, M., Datta, A. (2017). Improving food nutritional quality and productivity through genetic engineering. Int. J. Cell Sci. Mol. Biol. 2, 555576. doi: 10.19080/IJCSMB.2017.02.555576

CrossRef Full Text | Google Scholar

Jagtap, V., Bhargava, S., Streb, P., Feierabend, J. (1998). Comparative effect of water, heat and light stresses on photosynthetic reactions in Sorghum bicolor (L.) moench. J. Exp. Bot. 49, 1715–1721. doi: 10.1093/jxb/49.327.1715

CrossRef Full Text | Google Scholar

Jauregui, I., Aroca, R., Garnica, M., Zamarreño, Á.M., García-Mina, J. M., Serret, M. D., et al. (2015). Nitrogen assimilation and transpiration: Key processes conditioning responsiveness of wheat to elevated [CO2] and temperature. Physiologia Plantarum. 155, 338–354. doi: 10.1111/ppl.12345

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Huang, B. (2001). Drought and heat stress injury to two cool-season turfgrasses in relation to antioxidant metabolism and lipid peroxidation. Crop Sci. 41, 436–442. doi: 10.2135/cropsci2001.412436x

CrossRef Full Text | Google Scholar

Kahiluoto, H., Kaseva, J., Balek, J., Olesen, J. E., Ruiz-Ramos, M., Gobin, A., et al. (2019). Decline in climate resilience of european wheat. Proc. Natl. Acad. Sci. United States America. 116, 123–128. doi: 10.1073/pnas.1804387115

CrossRef Full Text | Google Scholar

Khadka, K., Earl, H. J., Raizada, M. N., Navabi, A. (2020). A physio-morphological trait-based approach for breeding drought tolerant wheat. Front. Plant Sci. 11. doi: 10.3389/fpls.2020.00715

CrossRef Full Text | Google Scholar

Kiani-Pouya, A., Roessner, U., Jayasinghe, N. S., Lutz, A., Rupasinghe, T., Bazihizina, N., et al. (2017). Epidermal bladder cells confer salinity stress tolerance in the halophyte quinoa and atriplex species. Plant Cell Environ. 40, 1900–1915. doi: 10.1111/pce.12995

PubMed Abstract | CrossRef Full Text | Google Scholar

Kirschbaum, M. U. F., McMillan, A. M. S. (2018). Warming and elevated CO2 have opposing influences on transpiration. Which is more important? Curr. Forestry Rep. 4, 51–71. doi: 10.1007/s40725-018-0073-8

CrossRef Full Text | Google Scholar

Klem, K., Gargallo-Garriga, A., Rattanapichai, W., Oravec, M., Holub, P., Veselá, B., et al. (2019). Distinct morphological, physiological, and biochemical responses to light quality in barley leaves and roots. Front. Plant Sci. 10. doi: 10.3389/fpls.2019.01026

PubMed Abstract | CrossRef Full Text | Google Scholar

Krasensky, J., Jonak, C. (2012). Drought, salt, and temperature stress-induced metabolic rearrangements and regulatory networks. J. Exp. Bot. 63, 1593–1608. doi: 10.1093/jxb/err460

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumari, M., Joshi, R., Kumar, R. (2020). Metabolic signatures provide novel insights to picrorhiza kurroa adaptation along the altitude in Himalayan region. Metabolomics. 16, 77. doi: 10.1007/s11306-020-01698-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, R., Joshi, R., Kumari, M., Thakur, R., Kumar, D., Kumar, S. (2020). Elevated CO2 and temperature influence key proteins and metabolites associated with photosynthesis, antioxidant and carbon metabolism in Picrorhiza kurroa. J. Proteomics. 219, 103755. doi: 10.1016/j.jprot.2020.103755

PubMed Abstract | CrossRef Full Text | Google Scholar

Kurepin, L. V., Ivanov, A. G., Zaman, M., Pharis, R. P., Hurry, V., Hüner, N. P. A. (2017). “Interaction of glycine betaine and plant hormones: protection of the photosynthetic apparatus during abiotic stress,” in Photosynthesis: Structures, mechanisms, and applications eds. Hou, H., Najafpour, M., Moore, G, Allakhverdiev, S. (Cham: Springer International Publishing), 185–202. doi: 10.1007/978-3-319-48873-8_9

CrossRef Full Text | Google Scholar

Levine, L. H., Kasahara, H., Kopka, J., Erban, A., Fehrl, I., Kaplan, F., et al. (2008). Physiologic and metabolic responses of wheat seedlings to elevated and super-elevated carbon dioxide. Adv. Space Res. 42, 1917–1928. doi: 10.1016/j.asr.2008.07.014

CrossRef Full Text | Google Scholar

Li, B., Fan, R., Sun, G., Sun, T., Fan, Y., Bai, S., et al. (2021). Flavonoids improve drought tolerance of maize seedlings by regulating the homeostasis of reactive oxygen species. Plant Soil. 461, 389–405. doi: 10.1007/s11104-020-04814-8

CrossRef Full Text | Google Scholar

Li, Y., Xu, S., Wang, Z., He, L., Xu, K., Wang, G. (2018). Glucose triggers stomatal closure mediated by basal signaling through HXK1 and PYR/RCAR receptors in arabidopsis. J. Exp. Bot. 69, 1471–1484. doi: 10.1093/jxb/ery024

PubMed Abstract | CrossRef Full Text | Google Scholar

Meehl, G. A., Stocker, T. F., Collins, W. D., Friedlingstein, P., Gaye, A. T., Gregory, J. M., et al (2007). "Global climate projections," in Climate change 2007: The physical science basis. contribution of working group I to the fourth assessment report of the intergovernmental panel on climate change, eds.Solomon, S., Qin, H. L. M. D., Manning, M., Chen, Z., Marquis, M., Averyt, K. B. (Cambridge, UK and New York, USA: Cambridge University Press).

Google Scholar

Munné-Bosch, S., Queval, G., Foyer, C. H. (2012). The impact of global change factors on redox signaling underpinning stress tolerance. Plant Physiol. 161, 5–19. doi: 10.1104/pp.112.205690

PubMed Abstract | CrossRef Full Text | Google Scholar

Naing, A. H., Kim, C. K. (2021). Abiotic stress-induced anthocyanins in plants: Their role in tolerance to abiotic stresses. Physiologia Plantarum. 172, 1711–1723. doi: 10.1111/ppl.13373

PubMed Abstract | CrossRef Full Text | Google Scholar

Nuttall, J. G., Barlow, K. M., Delahunty, A. J., Christy, B. P., O’Leary, G. J. (2018). Acute high temperature response in wheat. Agron. J. 110, 1296–1308. doi: 10.2134/agronj2017.07.0392

CrossRef Full Text | Google Scholar

Obata, T., Fernie, A. R. (2012). The use of metabolomics to dissect plant responses to abiotic stresses. Cell. Mol. Life Sci. 69, 3225–3243. doi: 10.1007/s00018-012-1091-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Obata, T., Witt, S., Lisec, J., Palacios-Rojas, N., Florez-Sarasa, I., Araus, J. L., et al. (2015). Metabolite profiles of maize leaves in drought, heat and combined stress field trials reveal the relationship between metabolism and grain yield. Plant Physiol. 169, 2665–2683. doi: 10.1104/pp.15.01164

PubMed Abstract | CrossRef Full Text | Google Scholar

Pandey, P., Irulappan, V., Bagavathiannan, M. V., Senthil-Kumar, M. (2017). Impact of combined abiotic and biotic stresses on plant growth and avenues for crop improvement by exploiting physio-morphological traits. Front. Plant Sci. 8. doi: 10.3389/fpls.2017.00537

CrossRef Full Text | Google Scholar

Pandey, P., Ramegowda, V., Senthil-Kumar, M. (2015). Shared and unique responses of plants to multiple individual stresses and stress combinations: Physiological and molecular mechanisms. Front. Plant Sci. 6. doi: 10.3389/fpls.2015.00723

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, Z., Chong, J., Zhou, G., de Lima Morais, D. A., Chang, L., Barrette, M., et al. (2021). MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 49, W388–W396. doi: 10.1093/nar/gkab382

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, M., Shahzad, R., Gul, A., Subthain, H., Shen, S., Lei, L., et al. (2017). Differentially evolved glucosyltransferases determine natural variation of rice flavone accumulation and UV-tolerance. Nature Communications. 8, 1975. doi: 10.1038/s41467-017-02168-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, S., Ma, K., Lu, Z., Chen, G., Cui, J., Tong, P., et al. (2019). Transcriptomic and metabolomic analysis of the heat-stress response of Populus tomentosa Carr. Forests. 10, 383. doi: 10.3390/f10050383

CrossRef Full Text | Google Scholar

Reyenga, P. J., Howden, S. M., Meinke, H., Hall, W. B. (2001). Global change impacts on wheat production along an environmental gradient in south Australia. Environ. Int. 27, 195–200. doi: 10.1016/S0160-4120(01)00082-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Rice-Evans, C., Miller, N., Paganga, G. (1997). Antioxidant properties of phenolic compounds. Trends Plant Sci. 2, 152–159. doi: 10.1016/S1360-1385(97)01018-2

CrossRef Full Text | Google Scholar

Rizhsky, L., Liang, H., Shuman, J., Shulaev, V., Davletova, S., Mittler, R. (2004). When defense pathways collide. The response of arabidopsis to a combination of drought and heat stress. Plant Physiol. 134, 1683–1696. doi: 10.1104/pp.103.033431

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosa, M., Prado, C., Podazza, G., Interdonato, R., González, J. A., Hilal, M., et al. (2009). Soluble sugars. Plant Signaling Behav. 4, 388–393. doi: 10.4161/psb.4.5.8294

CrossRef Full Text | Google Scholar

Sardans, J., Gargallo-Garriga, A., Urban, O., Klem, K., Walker, T. W. N., Holub, P., et al. (2020). Ecometabolomics for a better understanding of plant responses and acclimation to abiotic factors linked to global change. Metabolites. 10, 239. doi: 10.3390/metabo10060239

CrossRef Full Text | Google Scholar

Sharma, M., Irfan, M., Kumar, A., Kumar, P., Datta, A. (2021). “Recent insights into plant circadian clock response against abiotic stress,” Journal of Plant Growth Regulation. 221, 112403. doi: 10.1007/s00344-021-10531-y

CrossRef Full Text | Google Scholar

Sharma, M., Kumar, P., Verma, V., Sharma, R., Bhargava, B., Irfan, M. (2022). Understanding plant stress memory response for abiotic stress resilience: Molecular insights and prospects. Plant Physiol. Biochem. 179, 10–24. doi: 10.1016/j.plaphy.2022.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Shewry, P. R. (2009). Wheat. J. Exp. Bot. 60, 1537–1553. doi: 10.1093/jxb/erp058

PubMed Abstract | CrossRef Full Text | Google Scholar

Sicher, R. (2013). Combined effects of CO2 enrichment and elevated growth temperatures on metabolites in soybean leaflets: Evidence for dynamic changes of TCA cycle intermediates. Planta. 238, 369–380. doi: 10.1007/s00425-013-1899-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Šigut, L., Holisova, P., Klem, K., Prtova, M., Calfapietra, C., Marek, M. V., et al(2015). Does long-term cultivation of saplings under elevated CO2 concentrationinfluence their photosynthetic response to temperature? Ann. Bot. 116, 929–939. doi: 10.1093/aob/mcv043

PubMed Abstract | CrossRef Full Text | Google Scholar

Soba, D., Ben Mariem, S., Fuertes-Mendizábal, T., Méndez-Espinoza, A. M., Gilard, F., González-Murua, C., et al. (2019). Metabolic effects of elevated CO2 on wheat grain development and composition. J. Agric. Food Chem. 67, 8441–8451. doi: 10.1021/acs.jafc.9b01594

PubMed Abstract | CrossRef Full Text | Google Scholar

Tahjib-Ul-Arif, M., Zahan, M. I., Karim, M. M., Imran, S., Hunter, C. T., Islam, M. S., et al. (2021). Citric acid-mediated abiotic stress tolerance in plants. Int. J. Mol. Sci. 7235. doi: 10.3390/ijms22137235

PubMed Abstract | CrossRef Full Text | Google Scholar

Templer, S. E., Ammon, A., Pscheidt, D., Ciobotea, O., Schuy, C., McCollum, C., et al. (2017). Metabolite profiling of barley flag leaves under drought and combined heat and drought stress reveals metabolic QTLs for metabolites associated with antioxidant defense. J. Exp. Bot. 68, 1697–1713. doi: 10.1093/jxb/erx038

PubMed Abstract | CrossRef Full Text | Google Scholar

Thornton, P. K., Ericksen, P. J., Herrero, M., Challinor, A. J. (2014). Climate variability and vulnerability to climate change: A review. Global Change Biol. 20, 3313–3328. doi: 10.1111/gcb.12581

CrossRef Full Text | Google Scholar

Trovato, M., Funck, D., Forlani, G., Okumoto, S., Amir, R. (2021). Editorial: Amino acids in plants: Regulation and functions in development and stress defense. Front. Plant Sci. 12. doi: 10.3389/fpls.2021.772810

PubMed Abstract | CrossRef Full Text | Google Scholar

Walkowiak, S., Gao, L., Monat, C., Haberer, G., Kassa, M. T., Brinton, J., et al. (2020). Multiple wheat genomes reveal global variation in modern breeding. Nature. 588, 277–283. doi: 10.1038/s41586-020-2961-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, G., Johnson, S. K., Bornman, J. F., Bennett, S. J., Clarke, M. W., Singh, V., et al. (2016). Growth temperature and genotype both play important roles in sorghum grain phenolic composition. Sci. Rep. 6, 21835. doi: 10.1038/srep21835

PubMed Abstract | CrossRef Full Text | Google Scholar

Xalxo, R., Yadu, B., Chandra, J., Chandrakar, V., Keshavkant, S. (2020). “Alteration in carbohydrate metabolism modulates thermotolerance of plant under heat stress,” in Heat stress tolerance in plants eds. Wani, S.H., Kumar, V. (New Jersey, USA: Wiley), 77–115. doi: 10.1002/9781119432401.ch5

CrossRef Full Text | Google Scholar

Xu, G., Cao, J., Wang, X., Chen, Q., Jin, W., Li, Z., et al. (2019). Evolutionary metabolomics identifies substantial metabolic divergence between maize and its wild ancestor, teosinte. Plant Cell. 31, 1990–2009. doi: 10.1105/tpc.19.00111

PubMed Abstract | CrossRef Full Text | Google Scholar

Zandalinas, S. I., Balfagón, D., Gómez-Cadenas, A., Mittler, R. (2022). Plant responses to climate change: Metabolic changes under combined abiotic stresses. J. Exp. Bot. 73, 3339–3354. doi: 10.1093/jxb/erac073

PubMed Abstract | CrossRef Full Text | Google Scholar

Zandalinas, S. I., Mittler, R. (2022). Plant responses to multifactorial stress combination. New Phytol. 234, 1161–1167. doi: 10.1111/nph.18087

PubMed Abstract | CrossRef Full Text | Google Scholar

Zandalinas, S. I., Mittler, R., Balfagón, D., Arbona, V., Gómez-Cadenas, A. (2018). Plant adaptations to the combination of drought and high temperatures. Physiologia Plantarum. 162, 2–12. doi: 10.1111/ppl.12540

PubMed Abstract | CrossRef Full Text | Google Scholar

Zinta, G., AbdElgawad, H., Peshev, D., Weedon, J. T., Van den Ende, W., Nijs, I., et al. (2018). Dynamics of metabolic responses to periods of combined heat and drought in arabidopsis thaliana under ambient and elevated atmospheric CO2. J. Exp. Bot. 69, 2159–2170. doi: 10.1093/jxb/ery055

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: climate change, wheat, metabolomics, physiology, elevated CO2, temperature, drought

Citation: Večeřová K, Oravec M, Puranik S, Findurová H, Veselá B, Opoku E, Ofori-Amanfo KK, Klem K, Urban O and Sahu PP (2022) Single and interactive effects of variables associated with climate change on wheat metabolome. Front. Plant Sci. 13:1002561. doi: 10.3389/fpls.2022.1002561

Received: 25 July 2022; Accepted: 08 September 2022;
Published: 10 October 2022.

Edited by:

Mohammad Irfan, Cornell University, United States

Reviewed by:

Aiswarya Girija, Aberystwyth University, United Kingdom
Rajiv Kumar, Institute of Himalayan Bioresource Technology Council of Scientific and Industrial Research, (CSIR), India

Copyright © 2022 Večeřová, Oravec, Puranik, Findurová, Veselá, Opoku, Ofori-Amanfo, Klem, Urban and Sahu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Pranav Pankaj Sahu, sahu.p@czechglobe.cz

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.