Original Research ARTICLE
Metabolite Profiling of Barley Grains Subjected to Water Stress: To Explain the Genotypic Difference in Drought-Induced Impacts on Malting Quality
- 1Agronomy Department, Zhejiang University, Hangzhou, China
- 2Zhejiang Academy of Agricultural Sciences, Hangzhou, China
Grain weight and protein content will be reduced and increased, respectively, when barley is subjected to water stress after anthesis, consequently deteriorating the malt quality. However, such adverse impact of water stress differs greatly among barley genotypes. In this study, two Tibetan wild barley accessions and two cultivated varieties differing in water stress tolerance were used to investigate the genotypic difference in metabolic profiles during grain-filling stage under drought condition. Totally, 71 differently accumulated metabolites were identified, including organic acids, amino acids/amines, and sugars/sugar alcohols. Their relative contents were significantly affected by water stress for all genotypes and differed distinctly between the wild and cultivated barleys. The principal component analysis of metabolites indicated that the Tibetan wild barley XZ147 possessed a unique response to water stress. When subjected to water stress, the wild barley XZ147 showed the most increase of β-amylase activity among the four genotypes, as a result of its higher lysine content, less indole-3-acetic acid (IAA) biosynthesis, more stable H2O2 homeostasis, and more up-regulation of BMY1 gene. On the other hand, XZ147 had the most reduction of β-glucan content under water stress than the other genotypes, which could be explained by the faster grain filling process and the less expression of β-glucan synthase gene GSL7. All these results indicated a great potential for XZ147 in barley breeding for improving water stress tolerance.
Global warming and climate change have become a primary concern worldwide (Intergovernmental Panel on Climate Change [IPCC], 2012). Since the middle of the 20th century, there have been considerable changes in the nature of droughts, extreme weather events, and floods in many regions of the world, which caused marked damage to crop production and great threat to global food security (Wheeler and Von Braun, 2013; Lesk et al., 2016). It has been estimated that, during 1964–2007, droughts significantly reduced cereal production by 10% on average, and this percentage was increasing annually, due to the rising drought severity, increasing vulnerability, and exposure to drought (Lesk et al., 2016). In recent two decades, the effect of drought stress on crop growth, yield, and quality was increasingly becoming a major issue of scientific concerns (reviewed by Kang et al., 2009; Alqudah et al., 2010). And, breeding programs with the aim to enhance crop productivity under drought stress are a top priority in the era of climate change (Jagadish et al., 2015). The occurrence of water stress at the reproductive stage is the most critical, as it strongly impacts yield and seed quality. Understanding the response mechanism of drought stress at crop reproductive stage will help to partially address some of concerns for improving crop tolerance to drought stress and to minimize consequent impacts.
Barley is the fourth most important cereal crop worldwide in terms of planting area, and is primarily used for food, brewing, and animal feed (Druka et al., 2011). In last few decades, malt quality of barley grains was heavily addressed and found to be associated with many chemical constituents and enzymes, such as β-amylase activity and β-glucan content (Wei et al., 2009a). β-Amylase, which is a key factor affecting the capacity of starch degradation during grain germination (Beck and Ziegler, 1989), is closely related to diastatic power (DP) (Liu et al., 2005). There have two genes encoding β-amylase been identified: BMY1 and BMY2. Of them, BMY1 is mainly expressed during the gain filling stages, and plays a vital role in regulating the gain β-amylase activity and malting quality (Li et al., 2002). And, low β-glucan content is required for brewing as it favors high wort filtration rate and malt extract. In barley grains, the β-glucan content was found to be largely regulated by three genes, GSL1, GSL4, and GSL7 (Schober et al., 2009). Furthermore, it has been documented that both β-glucan content and β-amylase activity are not only genetically controlled, but also greatly affected by environmental factors (Zhang et al., 2001, 2006; Djukić and Knežević, 2014; Rakszegi et al., 2014). It was reported that β-glucan content in barley grains was highly reduced when plants suffered from heat or water stress during grain-filling process (Macnicol et al., 1993; Christensen and Scheller, 2012). On the other hand, total protein content in barley grains was generally increased under the drought or heat stress condition, consequently resulting in higher β-amylase activity (Todaka et al., 2000; Yin et al., 2002; Qi et al., 2006). Such reduction of β-glucan content and increase of β-amylase activity were also observed in our previous study on water stress at the grain-filling stage, and they were closely associated with the reduction of grain yield and increase of protein content, which would consequently cause deterioration and instability of grain yield and malt quality (Wu et al., 2015). Furthermore, our previous study also found such effect of water stress on the β-glucan content and β-amylase activity varied dramatically among barley genotypes, especially between the Tibetan wild barley and cultivated ones (Wu et al., 2015). Therefore, to achieve the stability of grain yield and malt quality under water stress condition, it is imperative to reveal the mechanisms in the adverse impact of water stress on malt quality, including β-glucan content and β-amylase activity, and the difference among barley genotypes in the response to water stress.
Metabolite profiling may provide a comprehensive approach to assess a broader spectrum of constituent analysis and has been proven to be a suitable tool for investigating the metabolites changes caused by genetic modification and environmental condition (Röhlig et al., 2009; Frank et al., 2011; Cañas et al., 2015). Currently, metabolomic analysis has been extensively conducted on the changes of metabolite profiles under water stress in many species, such as tomato (Rivero et al., 2015), soybean (Silvente et al., 2012), maize (Sun et al., 2015), and rice (Li et al., 2015). In barley, several metabolite profiling studies have also been performed to determine the influence of water stress on free amino acids (Lanzinger et al., 2015), abscisic acid, and the oxidative status (Thameur et al., 2014). However, to our knowledge, no such report has been addressed to the effect of water stress on malt quality. In the present study, four wild and cultivated barley genotypes contrasting in drought tolerance were used to compare metabolic changes in the response to water stress using gas chromatography-mass spectrometry (GC-MS). The objectives of the current work are to determine the mechanisms in the influence of water stress on β-amylase activity and β-glucan content, and to reveal the reasons why barley genotypes differ in their responses to water stress.
Materials and Methods
Plant Materials and Sampling
Four barley genotypes contrasting in drought tolerance were used in this study based on our previous research (Wu et al., 2015). Of them, XZ5 (sensitive to drought) and XZ147 (tolerant to drought) are the Tibetan wild barley accessions, while Triumph (Tr) is a cultivar and TL43 is an ABA-insensitive mutant derived from Tr (Romagosa et al., 2001). Ten seeds of each barley genotype were sown in a pot with 7.5 kg (7 l) sandy-clay soil in mid-November 2015 at Zijin’gang Campus, Zhejiang University, Hangzhou, China. Three weeks later, only five uniform healthy barley seedlings were preserved for the experiment. During the growth of barley seedlings, all the pots were well irrigated to the water content of 40% (equaling to water potential of -0.15 MPa). When half the plants of each barley genotype were at anthesis stage in March 2016, the spikes at the same stage were tagged for further measurements, and the water stress treatment (drought, abbreviation of T) was conducted subsequently by stopping water supply to make the soil water content drop to be around 14% (equaling to water potential of -0.75 MPa, taking approximately 3 days, Supplementary Figure S1) and maintained for another 14 days. Thereafter, the water supply was resumed to control level. For control condition (normal water level, abbreviation of CK), the water content in soil was continuously maintained around 40% (equaling to water potential of -0.15 MPa). During the water stress treatment, all the pots were irrigated everyday by a repeated watering-testing procedure to get the required soil water content (40% for the control and 14% for drought). For each time, 20 mL water was added to each pot without disturbing the plants and 10 min later soil water content was tested using HH2 Moisture Meter (Delta-T Devices, Cambridge, United Kingdom). Once the soil water content reached to the required value, the watering was stopped.
At the 3rd (3d) and 7th (7d) of water stress treatment, the tagged spikes of each genotype were collected from both control and water stress treatment, and immediately stored at -80°C for the subsequent analysis. At maturity, the spikes on main stems were harvested and oven-dried at 37°C. Oven-dried grains were ground with a cyclone mill (SCINO CT410, FOSS) equipped with a 0.5 mm sieve and the flour samples were stored at -20°C until to be analyzed.
Measurements of β-Amylase Activity and β-Glucan Content
β-Amylase activity was measured using a Betamyl Assay Kit (Megazyme International, Ireland Ltd.) according to McCleary and Codd’s (1989) method. Total β-glucan content was assayed using a commercial kit (Megazyme International, Ireland Ltd.) according to McCleary and Codd’s (1991) method.
Metabolites of barley grains were extracted according to Lisec et al. (2006) with some modification. In brief, finely grinded barley powder (approx. 100 mg) was mixed well with 1,400 μl of 100% methanol (pre-cooled at -20°C) and 60 μl of ribitol (0.2 mg/ml stock in dH2O, as an internal quantitative standard) in 2 ml centrifuge tube, and then placed in a shaking bath at 70°C for 10 min and centrifuged at 11,000 × g for 10 min. The supernatant was collected into a 10 ml centrifuge tube containing 750 μl chloroform (pre-cooled at -20°C) and 1,500 μl deionized water (4°C), vortexed for 30 s, and then centrifuged at 2,200 × g for 15 min. One hundred and fifty microliters of supernatant was collected into a new 1.5 ml centrifuge tube and vacuum-dried for 1 h. Thereafter, 40 μl methoxyamine pyridine solution (20 mg/ml) was added into the tube, incubated in a shaker at 37°C for 2 h, then 70 μl MSTFA reagent was added, and shaked for 30 min at 37°C. Metabolites contents were subsequently determined using Agilent 6890N GC/5975B MSD (Agilent, United States). The program of temperature rise was set as: initial temperature of 70°C for 2 min, 10°C/min rate up to 140°C, 4°C/min rate up to 240°C, 10°C/min rate up to 300°C, and then staying for 8 min.
The raw signals were imported into software AMDIS (Version 2.71) to search for metabolites from its default universal database. The total mass of signal integration area was normalized for each sample, with the total integral area of each sample being normalized to 10,00,000. Finally, the normalized data were imported into MetaboAnalyst online analysis software1, employing PLS-DA model and the first principal component of VIP (variable importance in the projection) values (VIP.1) combined with Student’s t-test (T-test) (p, 0.01), to find differentially accumulated metabolites (Lisec et al., 2006).
Quantitative Real-time PCR
The relative transcript level of the genes encoding β-amylase and β-glucan synthases was determined through quantitative real-time PCR (qRT-PCR). Total RNAs in barley grains were extracted using a RNA plant Plus Reagent Kit [Tiangen Biotech (Beijing) Co. Ltd.] according to the instruction. RNA was reverse-transcripted to cDNAs using a PrimeScript RTF Reagent Kit with gDNA Eraser (Takara Bio Inc.), and then stored at -20°C for subsequent PCR analysis. The qRT-PCR was performed on CFX96 TouchTM Real-Time PCR Detection System (Bio-Rad Laboratories, Inc.) with an iTaq Universal SYBR Green Supermix (Bio-Rad Laboratories, Inc.). Primer sequences for qRT-PCR analysis were listed in Supplementary Table S1.
Prior to data analysis, quantitative normalization within replicates was transformed by logarithmic base of 2. MetaboAnalyst online analysis software2 was used to build heatmap diagram (Xia and Wishart, 2011). Meanwhile, hierarchical cluster analysis (HCA) and principal component analysis (PCA) models were performed for all samples. The different significance of metabolites between treatment and control or among genotypes was tested using T-test and ANOVA analysis on SPSS 20.0 software.
Effect of Water Stress on Grain Weight, β-Amylase Activity, and β-Glucan Content
Drought stress dramatically reduced grain weight of barley genotypes Tr, TL43, and XZ5 with XZ5 showing the most reduction, whereas no significant effect of drought on grain weight was observed for XZ147 (Figure 1A). Similarly, drought stress caused a significant decrease in β-glucan content of barley grains, in particular for XZ147 (37.38%) (Figure 1B). On the other hand, β-amylase activity was remarkably increased in XZ147, TL43, and Tr by 93.76, 37.57, and 32.55% under water stress compared with the control, but only a tiny increase (11.52%) in β-amylase activity was seen for XZ5 (Figure 1C).
FIGURE 1. The effect of water stress on grain weight (A), β-amylase activity (B), and β-glucan content (C) of different barley genotypes. Different letters indicate the significant difference between each genotype × treatment combination at 95% probability. Data are mean ± SD.
Changes of Metabolites in the Response to Water Stress
In the present study, the metabolic profiles in barley grains at filling stage were investigated. Totally, 71 differently accumulated metabolites were successfully identified, and their relative contents were significantly affected by water stress for all genotypes (Figure 2). According to the chemical profile, these metabolites could be classified into three fractions: organic acids (fraction I), amino acids and amines (fraction II), and sugars and alcohols (fraction III). In order to determine the effect of water stress on the metabolites, PCA was conducted separately on each metabolite fraction. The two treatments (water stress and control) and four genotypes were clearly separated by two principal components (Figure 3). For organic acid fraction, samples from the control and water stress for XZ5, Tr, and TL43 were clearly separated by PC1, which could explain 51.3% of the total variation, whereas the two treatments for XZ147 were separated by PC2, which only explained 18.7% of the total variation (Figure 3A). For fraction of amino acids and amines, PC1 clearly separated the samples of the two water treatments for all barley genotypes, explaining 79.2% of the total variation; and PC2 obviously separated the wild (XZ5 and XZ147) and cultivated genotypes (Tr and TL43) under the condition of water stress, explaining 10.3% of the total variation (Figure 3B). Furthermore, PCAs based on the metabolite data of sugars and alcohols (fraction III) showed a strong influence of water stress, and could explain 38.2 and 26.7% of the total variation, respectively (Figure 3C).
FIGURE 2. Hierarchical cluster analysis of 71 metabolites and 8 genotype × treatment combinations (distance measure using Euclidean and clustering algorithm using Ward.D). CK, control condition with soil water content of 40%; T, water stress treatment with soil water content of 14% and Tr, Triumph. The scale –3 (dark blue, the lowest) to 3 (dark red, the highest) indicated the relative content of each metabolite for each genotype × treatment combination.
FIGURE 3. Principal component analysis (PCA) of metabolic profiles in grains of four barley genotypes under the two water treatments. (A) PCA of fraction I, organic acid; (B) PCA of fraction II, amino acids and amines; and (C) PCA of fraction III, sugars and sugar alcohols. CK, control condition with soil water content of 40%; T, water stress treatment with soil water content of 14%; Tr, Triumph; PC1, the first principal component; and PC2, the second principal component.
The loading plots revealed that constituents from all identified substance classes (organic acids, amino acids and amines, and sugars and alcohols) were responsible for the variation (Figure 4). Fraction I (acids) and fraction II (amino acids and amines) were the major prominent factors for shifting on PC1 and PC2 (Figures 4A,B). The loading plots of fraction III (sugars and alcohols) showed that monosaccharides, including fructose, sucrose, and maltose, as well as alcohol and mannitol were obviously increased under water stress (Figure 4C and Table 1). Furthermore, the majority of metabolites in fractions I and II were also affected by water stress (Figures 4A,B and Table 1).
FIGURE 4. The corresponding loading plots of metabolite profiling data from fraction I (organic acid, A), fraction II (amino acids and amines, B), and fraction III (sugars and sugar alcohols, C) between PC1 and PC2.
TABLE 1. Relative concentration and the fold changes of the major metabolites in the grains of the different barley genotypes under water stress.
The relative contents and fold changes [calculated using formula: log2 (drought/control)] of the 25 dominant metabolites in grains were listed in Table 1. Among them, the relative contents of 19 metabolites were increased by water stress for all barley genotypes, including several osmoprotectants like mannitol and L-proline (Table 1). However, such water stress-induced increase in these metabolites showed an obvious genotypic difference, with the wild barley XZ5 showing the most increase in these metabolites except mannitol. On the other hand, the relative content of 2-ketoglutaric acid was reduced in all genotypes except TL43 under water stress, and the reduced extent differed dramatically between the two wild barleys (-1.44-fold in XZ5 and -1.32-fold in XZ147) and the cultivar Tr (-0.30-fold) (Table 1). Furthermore, there were five metabolites showing different changes in response to water stress between the four barley genotypes. For instance, water stress significantly reduced the relative content of maltose in XZ5 by -0.48-fold, but no significant change in it was seen for the other three genotypes. The effect of water stress on myo-inositol content differed among barley genotypes, with XZ147 being reduced (-0.40-fold); and XZ5, TL43, and Tr being increased (0.27-, 0.34-, and 0.38-fold, respectively). Similarly, water stress reduced xylitol content in wild barleys (-1.51-fold in XZ5 and -3.75-fold in XZ147) but increased it in the cultivated genotype TL43 (1.22-fold), and caused little change in Tr (0.11-fold).
Effect of Water Stress on Metabolic Pathway
The metabolites with significant changes under drought stress were illustrated on the metabolic pathway (Figure 5). The elevated level of sucrose and TCA cycle components (citric acid, succinic acid, fumaric acid, malic acid) probably indicated the enhanced energy metabolism (Warth et al., 2015) and elevated respiratory rates (Bolton, 2009) for plant defense actions, such as drought stress protective program (Figures 2, 5). Along with the increase of intermediates in TCA cycle, several amino acids displayed increased level at different extent in response to drought stress, which is always associated with different plant defense mechanisms (Figures 2, 5; Warth et al., 2015). Quinic acid – the precursor of the shikimate pathway showed higher accumulation under drought of all genotypes, especially in wild genotypes. This pathway provides important aromatic secondary metabolites such as phenylpropanoids as well as the plant hormone auxins (Figure 5; Sugawara et al., 2009; Zhao, 2010; Mashiguchi et al., 2011). In mitochondrion of plants, arginine is transferred to ornithine and urea by arginase, which is crucial for the mobilization of nitrogen (Witte, 2011). In our study, surprisingly, only in wild barley XZ147 the urea biosynthesis was inhibited, indicating that XZ147 processes a unique nitrogen metabolic way different from the other genotypes (Figure 5). Furthermore, the enhanced glutamate recycling could significantly maintain GSH levels under abiotic stress (Figures 2, 5; Paulose et al., 2013). These results revealed that drought stress was able to modify both the primary carbohydrate metabolism and the primary nitrogen metabolism, and induce several defense responses as well.
FIGURE 5. Schematic overview of metabolic pathways in response to water treatments for different barley genotypes. CK, control condition with soil water content of 40%; T, water stress treatment with soil water content of 14%; and Tr, Triumph. Black arrows indicate specific metabolic steps.
Effect of Water Stress on Transcriptional Levels of the Genes Encoding β-Amylase and β-Glucan Synthases
The effect of water stress on the expression of genes encoding β-amylase and β-glucan synthesis was presented in Figure 6. It could be seen that the response of the expression of four examined genes to water stress was highly duration- and genotype-dependent. At 3d of water stress, transcriptional level of BMY1 gene (encoding β-amylase) was little affected. However, at 7d, expression of BMY1 was distinctly increased in XZ5 (2.3-fold), XZ147 (3.3-fold), and TL43 (2.1-fold), but no significant change was observed for Tr (Figure 6A) in comparison with control. The transcriptional levels of the three genes encoding β-glucan synthase, GSL1, GSL4, and GSL7, differed greatly among genotypes under water stress. The effect of water stress on expressional patterns of GSL1 and GSL4 was quite similar (Figures 6B,C). At 3d of water stress, only Tr increased the expression of GSL1 by 2.1-fold, while at 7d, the expression of both GSL1 and GSL4 in XZ5 was increased by 2.5- and 1.9-fold, respectively. The effect of water stress on GSL7 expression varied greatly over the time of water stress. At 3d of water stress treatment, the expression of GSL7 was significantly increased in XZ5 (2.9-fold), TL43 (4.3-fold), and Tr (6.0-fold) but dramatically decreased in XZ147 (only 3.8% of the control) (Figure 6D). With the exposure time increased to 7d, the water stress caused up-regulation of GSL7 in XZ5, TL43, and Tr and the down-regulation of GSL7 in XZ147 were both reduced (Figure 6D). In addition, it was easily found that the expression level of GSL7 was much lower in wild barley than in cultivated ones (Figure 6D), being coincide with the results of the β-glucan content (Figure 1B).
FIGURE 6. Changes in transcriptional levels of BMY1 (A), GSL1 (B), GSL4 (C), and GSL7 (D) in caryopsis of different barley genotypes after onset of water stress for 3d and 7d. GSLs, β-glucan synthase genes and BMY1, β-amylase gene 1. Different letters indicate the significant difference at 95% probability. Data are mean ± SD.
It is well documented that metabolic changes happen in the plants exposed to water stress, resulting in yield loss (Boyer, 1982; Bray et al., 2000). In this study, water stress induced a significant reduction of grain weight for all genotypes, compared with the control (Figure 1), being consistent with the results obtained in our previous work (Wu et al., 2015). However, such reduction of grain weight differed greatly between genotypes, with XZ147 being the least and XZ5 being the most. It was reported that carbohydrates and starch account for 78–83 (MacGregor and Fincher, 1993) and 50–70% (Henry, 1988) of barley grain weight, respectively. Photosynthetic product (glucose) is transported to grains in the form of sucrose, and used for synthesis of starch as well as β-glucan. The inhibited synthesis of starch by water stress may be a major reason for the reduction of grain weight (Chaves, 1991). In this study, the available sucrose transported to grains was dramatically increased under water stress for all barley genotypes, with XZ147 having the most increase (Figures 2, 5). Thanks to its much quicker grain-filling process under water stress (Wu et al., 2017), XZ147 remained relatively smaller change in grain weight than other genotypes (Figure 1A).
In the present study, the PCA revealed that the clustering between barley genotypes strongly depended on fraction II (amino acids and amines) and fraction III (sugars and sugar alcohols) (Figures 4B,C). So, it can be assumed that these low molecular compounds (amino acids, amines, sugars, and sugar alcohols) played an important role in barley’s tolerance to water stress. Indeed, both our and numerous previous studies have found that the content of amino acids in plant vegetative tissues and reproductive grains changed remarkably when subjected to water stress (Widodo et al., 2009; Sicher et al., 2012; Nam et al., 2014; Table 1). It is well known that osmotic adjustment by the accumulation of compatible solutes is a key mechanism for maintaining cell turgor under water stress (Serraj and Sinclair, 2002; Hummel et al., 2010). Proline is one of the main osmoprotectant in plants when subjected to osmotic stress, such as drought or salinity (Delauney and Verma, 1993; Liu and Zhu, 1997). It has been also reported that the accumulation of mannitol and inositol increased in plants when exposed to osmotic stress like salinity (Abebe et al., 2003; Sanchez et al., 2008). In this study, the contents of both proline and mannitol were increased by water stress in all used barley genotypes, whereas the drought tolerant wild barley XZ147 showed the least increase among all genotypes (Table 1). Moreover, XZ147 showed the reduction in inositol content under water stress relative to the control, while other three genotypes showed the increase of this metabolite. Very similar results were also obtained in the study on the metabolic responses to salt stress of barley (Widodo et al., 2009). These results suggested that the increasing accumulation of osmoprotectants like proline and mannitol is an important adaptive strategy of barley to survive under the terminate drought stress, especially for the sensitive genotypes. Likewise, synthesis of some amino acids and organic acids was significantly enhanced under water stress (Table 1), which could be in favor of osmotic adjustment and membrane stability (Ueda et al., 2004; Chen et al., 2007; Widodo et al., 2009). It has been previously documented that poly-amines (PAs) associated with some related amino acids play important roles in the adaptation to abiotic and biotic stresses by regulating carbon/nitrogen homeostasis or acting as signaling molecules and compatible solutes under drought stress (Moschou et al., 2008; Liu et al., 2010; Hussain et al., 2011; Zeier, 2013). However, the biosynthesis of these compatible solutes by plants is always at the high cost of photo-assimilates and energy, which consequently sacrifice the grain yield and quality (Bolton, 2009). In this study, we found a much higher increased accumulation of osmoprotectants (such as proline, mannitol, fructose, sucrose, and citric acid) in drought-sensitive wild XZ5 than drought-tolerant wild XZ147 (Table 1), which consequently severely affected the yield formation and the two malt quality traits (β-amylase activity and β-glucan content) (Figure 1), completely coinciding with the findings of Serraj and Sinclair (2002).
β-Amylase activity is an important malt quality parameter, which is positively correlated with DP (Ovesná et al., 2012). Furthermore, a number of studies have demonstrated that β-amylase induction could help plants cope with unfavorable growing conditions (Dreier et al., 1995; Nielsen et al., 1997; Todaka et al., 2000; Kaplan and Guy, 2004), possibly by hydrolyzing more starch to maltose to function as a compatible-solute stabilizing factor (Lu and Sharkey, 2004; Kaplan and Guy, 2005). Hejgaard and Boisen (1980) found that the genotypes rich in lysine had higher β-amylase activity because of more serine protease inhibitor (Z protein) and chymotrypsin inhibitors (CI-1 and CI2). In this study, the grain lysine content was dramatically increased by water stress especially for wild genotypes (Figures 2, 5), being consistent with the results reported previously (Macnicol et al., 1993; MacGregor et al., 1994; Wei et al., 2009b). So, it could be assumed that the increased grain lysine content is an important reason for the increased β-amylase activity under water stress. On the other hand, during grain filling stage of cereal crops or ripening stage of fruits, the increased levels of auxins could be observed along with the inhibition of β-amylase activity for starch accumulation (Eeuwens and Schwabe, 1975; Bangerth et al., 1985; Aufhammer et al., 1986; Purgatto et al., 2001). In the present study, we found that the content of L-Tryptophan, the precursor for auxin biosynthesis (Mano and Nemoto, 2012) was extensively increased by water stress for all used barley genotypes except XZ147. This might be one other explanation for the higher β-amylase activity in XZ147 than in the other genotypes under water stress. Abiotic stress like drought stress always causes excess reactive oxygen species (ROS) accumulation, which will result in cellular damage and consequently yield losses. Fortunately, plants have obtained numerous antioxidants to scavenge the excess ROS, for instance GSH and AsA, which play a key role in intracellular ROS level regulation (Apel and Hirt, 2004; Sairam and Tyagi, 2004; Møller et al., 2007). On the other hand, Wei et al. (2009b) reported that exogenous H2O2 treatment could elevate barley β-amylase activity, suggesting that the enhanced grain H2O2 concentration would increase β-amylase activity under osmotic stress. However, the increase in GSH concentration might extinguish this effect. In the present study, we found that the precursor of GSH biosynthesis, L-glutamine, was increased by water stress in all used barley genotypes, with XZ5 being significantly higher than the other genotypes. Such over accumulation of GSH in XZ5 would be helpful for the scavenging of ROS, but not benefit for the induction of β-amylase activity (Figure 1C). All these results indicated that the induction of barley grain β-amylase activity was highly correlated with the lysine content, indole-3-acetic acid (IAA) biosynthesis, and H2O2 homeostasis. Furthermore, our results revealed that the expression level of BMY1, which encodes β-amylase in barley grains, was dramatically induced in XZ147 after 7 days of water stress being accompanied by higher β-amylase activity (Figures 1C, 6A). Therefore, it can be hypothesized that the reason for the higher grain β-amylase activity in XZ147 than the other genotypes under water stress might be attributed to its higher lysine content, less IAA biosynthesis, more stable H2O2 homeostasis, and also the higher up-regulation of BMY1 gene.
β-Glucan content is a critical parameter for malt brewing, with high β-glucan content always resulting in slower wort filtration rate and beer haze (Bamforth, 1982). It has been reported that β-glucan accumulated in barley grains throughout the whole filling stage (Christensen and Scheller, 2012). Thus, it may be understandable that the faster grain filling process causes the lower β-glucan content. Indeed, our results revealed that the grain β-glucan content in wild barley (XZ147 and XZ5) was decreased much more than that in cultivated barley (Tr and TL43) when subjected to water stress (Figure 1B), completely coinciding with our previous results that the grain filling process of wild barley (XZ147 and XZ5) was much faster than the cultivated barley (Tr and TL43) after onset of drought stress (Wu et al., 2017). Furthermore, it has been also reported that the synthesis of β-glucan in barley grains was largely controlled by the genes GSL1, GSL4, and GSL7 (Schober et al., 2009). In the present study, we found that, in comparison with the cultivated barley, the expression level of GSL7 in wild barley was much lower and it was reduced much more with the prolongation of water stress (Figure 6D), which could account for the lower β-glucan content in the grains of the wild barley (Figure 1B).
In summary, water stress remarkably altered grain weight, β-amylase activity, and β-glucan content of the four genotypes used in this study. When subjected to water stress, the wild barley XZ147 showed the most increase of β-amylase activity among the four genotypes which might be attributed to its higher lysine content, less IAA biosynthesis, more stable H2O2 homeostasis, and more up-regulation of BMY1 gene. XZ147 also had the most reduction of β-glucan content under water stress than the other genotypes, which could be explained by the faster grain filling process and the less expression of β-glucan synthase gene GSL7. Obviously, XZ147 is interesting for malt barley breeders to develop the new cultivars with high water stress tolerance and stable malt quality.
FZ and GZ designed the project; XW and KC carried out the experimental work; XW and FZ wrote the paper with contributions from GZ; and all authors read and approved the paper.
This research was supported by the National Natural Science Foundation of China (31201166, 31330055), the China Agriculture Research System (CARS-05), and Jiangsu Collaborative Innovation Center for Modern Crop Production (JCIC-MCP).
Conflict of Interest Statement
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.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.01547/full#supplementary-material
Abebe, T., Guenzi, A. C., Martin, B., and Cushman, J. C. (2003). Tolerance of mannitol- accumulating transgenic wheat to water stress and salinity. Plant Physiol. 131, 1748–1755. doi: 10.1104/pp.102.003616
Alqudah, A. M., Samarah, N. H., and Mullen, R. E. (2010). “Drought stress effect on crop pollination, seed set, yield and quality,” in Alternative Farming Systems, Biotechnology, Drought Stress and Ecological Fertilization, ed. E. Lichtfouse (Berlin: Springer), 193–213.
Aufhammer, W., Zinsmaier, P., and Bangerth, F. (1986). Variation of dry matter accumulation at definite positions within wheat ears and levels of indole-3yl acetic acid (IAA). Plant Growth Regul. 4, 305–310. doi: 10.1007/BF00024930
Bray, E. A., Bailey-Serres, J., and Weretilnyk, E. (2000). “Responses to abiotic stresses,” in Biochemistry and Molecular Biology of Plants, eds W. Gruissem, B. Buchannan, and R. Jones (Rockville, MD: American Society of Plant Physiologists), 1158–1249.
Cañas, R. A., Canales, J., Muñoz-Hernández, C., Granados, J. M., Ávila, C., García-Martín, M. L., et al. (2015). Understanding developmental and adaptive cues in pine through metabolite profiling and co-expression network analysis. J. Exp. Bot. 66, 3113–3127. doi: 10.1093/jxb/erv118
Chen, Z., Pottosin, I. I., Cuin, T. A., Fuglsang, A. T., Tester, M., Jha, D., et al. (2007). Root plasma membrane transporters controlling K+/Na+ homeostasis in salt-stressed barley. Plant Physiol. 145, 1714–1725. doi: 10.1104/pp.107.110262
Djukić, N. H., and Knežević, D. S. (2014). Molecular characterization and genetic diversity analysis β-glucan content variability in grain of oat (Avena sativa L.). Genetika 46, 529–536. doi: 10.2298/GENSR1402529D
Dreier, W., Schnarrenberger, C., and Börner, T. (1995). Light-and stress-dependent enhancement of amylolytic activities in white and green barley leaves: β-amylases are stress-induced proteins. J. Plant Physiol. 145, 342–348. doi: 10.1016/S0176-1617(11)81900-0
Druka, A., Sato, K., and Muehlbauer, G. J. (2011). “Genome analysis: the state of knowledge of barley genes,” in Barley: Production, Improvement, and Uses, ed. S. E. Ullrich (Hoboken, NJ: Blackwell Publishing), 85–111.
Hejgaard, J., and Boisen, S. (1980). High-lysine proteins in Hiproly barley breeding: identification, nutritional significance and new screening methods. Hereditas 93, 311–320. doi: 10.1111/j.1601-5223.1980.tb01369.x
Hummel, I., Pantin, F., Sulpice, R., Piques, M., Rolland, G., Dauzat, M., et al. (2010). Arabidopsis plants acclimate to water deficit at low cost through changes of carbon usage: an integrated perspective using growth, metabolite, enzyme, and gene expression analysis. Plant Physiol. 154, 357–372. doi: 10.1104/pp.110.157008
Hussain, S. S., Ali, M., Ahmad, M., and Siddique, K. H. (2011). Polyamines: natural and engineered abiotic and biotic stress tolerance in plants. Biotechnol. Adv. 29, 300–311. doi: 10.1016/j.biotechadv.2011.01.003
Jagadish, K. S., Kishor, P. B. K., Bahuguna, R. N., von Wirén, N., and Sreenivasulu, N. (2015). Staying alive or going to die during terminal senescence-An enigma surrounding yield stability. Front. Plant Sci. 6:1070. doi: 10.3389/fpls.2015.01070
Kaplan, F., and Guy, C. L. (2005). RNA interference of Arabidopsis beta-amylase8 prevents maltose accumulation upon cold shock and increases sensitivity of PSII photochemical efficiency to freezing stress. Plant J. 44, 730–743. doi: 10.1111/j.1365-313X.2005.02565.x
Lanzinger, A., Frank, T., Reichenberger, G., Herz, M., and Engel, K. H. (2015). Metabolite profiling of barley grain subjected to induced drought stress: responses of free amino acids in differently adapted cultivars. J. Agric. Food Chem. 63, 4252–4261. doi: 10.1021/acs.jafc.5b01114
Li, C. D., Langridge, P., Zhang, X. Q., Eckstein, P. E., Rossnagel, B. G., Lance, R. C. M., et al. (2002). Mapping of barley (Hordeum vulgare L.) β-amylase alleles in which an amino acid substitution determines β-amylase isoenzyme type and the level of free β-amylase. J. Cereal Sci. 35, 39–50. doi: 10.1006/jcrs.2001.0398
Li, X., Lawas, L. M. F., Malo, R., Glaubitz, U., Erban, A., Mauleon, R., et al. (2015). Metabolic and transcriptomic signatures of rice floral organs reveal sugar starvation as a factor in reproductive failure under heat and drought stress. Plant Cell Environ. 38, 2171–2192. doi: 10.1111/pce.12545
Lisec, J., Schauer, N., Kopka, J., Willmitzer, L., and Fernie, A. R. (2006). Gas chromatography mass spectrometry–based metabolite profiling in plants. Nat. Protoc. 1, 387–396. doi: 10.1038/nprot.2006.59
Liu, G., Ji, Y., Bhuiyan, N. H., Pilot, G., Selvaraj, G., Zou, J., et al. (2010). Amino acid homeostasis modulates salicylic acid–associated redox status and defense responses in Arabidopsis. Plant Cell 22, 3845–3863. doi: 10.1105/tpc.110.079392
Liu, Z., Cheng, F., and Zhang, G. (2005). Grain phytic acid content in japonica rice as affected by cultivar and environment and its relation to protein content. Food Chem. 89, 49–52. doi: 10.1016/j.foodchem.2004.01.081
MacGregor, A. W., Dushnicky, L. G., Schroeder, S. W., and Ballance, G. M. (1994). Changes in barley endosperms during early stages of germination. J. Inst. Brew. 100, 85–90. doi: 10.1002/j.2050-0416.1994.tb00811.x
MacGregor, A. W., and Fincher, G. B. (1993). “Carbohydrates of the barley grain,” in Barley Chemistry and Technology, eds A. W. MacGregor and R. S. Bhatty (St Paul, MN: American Association of Cereal Chemists), 73–130.
Macnicol, P. K., Jacobsen, J. V., Keys, M. M., and Stuart, I. M. (1993). Effects of heat and water stress on malt quality and grain parameters of Schooner barley grown in cabinets. J. Cereal Sci. 18, 61–68. doi: 10.1006/jcrs.1993.1034
Mashiguchi, K., Tanaka, K., Sakai, T., Sugawara, S., Kawaide, H., Natsume, M., et al. (2011). The main auxin biosynthesis pathway in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 108, 18512–18517. doi: 10.1073/pnas.1108434108
Moschou, P. N., Paschalidis, K. A., Delis, I. D., Andriopoulou, A. H., Lagiotis, G. D., Yakoumakis, D. I., et al. (2008). Spermidine exodus and oxidation in the apoplast induced by abiotic stress is responsible for H2O2 signatures that direct tolerance responses in tobacco. Plant Cell 20, 1708–1724. doi: 10.1105/tpc.108.059733
Nam, K. H., Kim, D. Y., Shin, H. J., Nam, K. J., An, J. H., Pack, I. S., et al. (2014). Drought stress-induced compositional changes in tolerant transgenic rice and its wild type. Food Chem. 153, 145–150. doi: 10.1016/j.foodchem.2013.12.051
Ovesná, J., Kučera, L., Vaculová, K., Štrymplová, K., Svobodova, I., and Milella, L. (2012). Validation of the β-amy1 transcription profiling assay and selection of reference genes suited for a RT-qPCR assay in developing barley caryopsis. PLOS ONE 7:e41886. doi: 10.1371/journal.pone.0041886
Paulose, B., Chhikara, S., Coomey, J., Jung, H., Vatamaniuk, O., and Dhankher, O. P. (2013). A γ-glutamyl cyclotransferase protects Arabidopsis plants from heavy metal toxicity by recycling glutamate to maintain glutathione homeostasis. Plant Cell 25, 4580–4595. doi: 10.1105/tpc.113.111815
Purgatto, E., Lajolo, F. M., Oliveira Do Nascimento, J. R., and Cordenunsi, B. R. (2001). Inhibition of β-amylase activity, starch degradation and sucrose formation by indole-3-acetic acid during banana ripening. Planta 212, 823–828. doi: 10.1007/s004250000441
Qi, J. C., Zhang, G. P., and Zhou, M. X. (2006). Protein and hordein content in barley seeds as affected by nitrogen level and their relationship to beta-amylase activity. J. Cereal Sci. 43, 102–107. doi: 10.1016/j.jcs.2005.08.005
Rakszegi, M., Lovegrove, A., Balla, K., Láng, L., Bedő, Z., Veisz, O., et al. (2014). Effect of heat and drought stress on the structure and composition of arabinoxylan and β-glucan in wheat grain. Carbohydr. Polym. 102, 557–565. doi: 10.1016/j.carbpol.2013.12.005
Romagosa, I., Prada, D., Moralejo, M. A., Sopena, A., Muñoz, P., Casas, A. M., et al. (2001). Dormancy, ABA content and sensitivity of a barley mutant to ABA application during seed development and after ripening. J. Exp. Bot. 52, 1499–1506. doi: 10.1093/jexbot/52.360.1499
Sanchez, D. H., Siahpoosh, M. R., Roessner, U., Udvardi, M., and Kopka, J. (2008). Plant metabolomics reveals conserved and divergent metabolic responses to salinity. Physiol. Plant. 132, 209–219. doi: 10.1111/j.1399-3054.2007.00993.x
Schober, M. S., Burton, R. A., Shirley, N. J., Jacobs, A. K., and Fincher, G. B. (2009). Analysis of the (1,3)-β-glucan synthase gene family of barley. Phytochemistry 70, 713–720. doi: 10.1016/j.phytochem.2009.04.002
Sicher, R. C., Timlin, D., and Bailey, B. (2012). Responses of growth and primary metabolism of water-stressed barley roots to rehydration. J. Plant Physiol. 169, 686–695. doi: 10.1016/j.jplph.2012.01.002
Silvente, S., Sobolev, A. P., and Lara, M. (2012). Metabolite adjustments in drought tolerant and sensitive soybean genotypes in response to water stress. PLOS ONE 7:e38554. doi: 10.1371/journal.pone.0038554
Sugawara, S., Hishiyama, S., Jikumaru, Y., Hanada, A., Nishimura, T., Koshiba, T., et al. (2009). Biochemical analyses of indole-3-acetaldoxime-dependent auxin biosynthesis in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 106, 5430–5435. doi: 10.1073/pnas.0811226106
Thameur, A., Ferchichi, A., and López-Carbonell, M. (2014). Involvement of abscisic acid metabolites and the oxidative status of barley genotypes in response to drought. Can. J. Plant Sci. 94, 1481–1490. doi: 10.4141/cjps-2014-069
Ueda, A., Kathiresan, A., Inada, M., Narita, Y., Nakamura, T., Shi, W., et al. (2004). Osmotic stress in barley regulates expression of a different set of genes than salt stress does. J. Exp. Bot. 55, 2213–2218. doi: 10.1093/jxb/erh242
Warth, B., Parich, A., Bueschl, C., Schoefbeck, D., Neumann, N. K. N., Kluger, B., et al. (2015). GC–MS based targeted metabolic profiling identifies changes in the wheat metabolome following deoxynivalenol treatment. Metabolomics 11, 722–738. doi: 10.1007/s11306-014-0731-1
Wei, K., Dai, F., Wu, F., and Zhang, G. (2009a). The variation of β-amylase activity and protein fractions in barley grains as affected by genotypes and post-anthesis temperatures. J. Inst. Brew. 115, 208–213. doi: 10.1002/j.2050-0416.2009.tb00370.x
Wei, K., Jin, X., Chen, X., Wu, F., Zhou, W., Qiu, B., et al. (2009b). The effect of H2O2 and abscisic acid (ABA) interaction on β-amylase activity under osmotic stress during grain development in barley. Plant Physiol. Biochem. 47, 778–784. doi: 10.1016/j.plaphy.2009.05.005
Widodo, Patterson, J. H., Newbigin, E., Tester, M., Bacic, A., and Roessner, U. (2009). Metabolic responses to salt stress of barley (Hordeum vulgare L.) cultivars, Sahara and Clipper, which differ in salinity tolerance. J. Exp. Bot. 60, 4089–4103. doi: 10.1093/jxb/erp243
Wu, X. J., Chen, X., Zeng, F. R., and Zhang, G. P. (2015). The genotypic difference in the effect of water stress after anthesis on the malt quality parameters in barley. J. Cereal Sci. 65, 209–214. doi: 10.1016/j.jcs.2015.07.012
Wu, X. J., Zeng, F. R., and Zhang, G. P. (2017). PEG-simulated drought stress and spike in vitro culture are used to study the impact of water stress on barley malt quality. Plant Growth Regul. 81, 243–252. doi: 10.1007/s10725-016-0201-z
Yin, C., Zhang, G. P., Wang, J. M., and Chen, J. X. (2002). Variation of beta-amylase activity in barley as affected by cultivar and environment and its relation to protein content and grain weight. J. Cereal Sci. 36, 307–312. doi: 10.1006/jcrs.2002.0467
Zhang, G. P., Chen, J. X., Dai, F., Wang, J. M., and Wu, F. B. (2006). The effect of cultivar and environment on β-amylase activity is associated with the change of protein content in barley grains. J. Agron. Crop. Sci. 192, 43–49. doi: 10.1111/j.1439-037X.2006.00181.x
Zhang, G. P., Chen, J. X., Wang, J. M., and Ding, S. (2001). Cultivar and environmental effects on (1→3, 1→4)-β-D-glucan and protein content in malting barley. J. Cereal Sci. 34, 295–301. doi: 10.1006/jcrs.2001.0414
Keywords: barley (Hordeum vulgare L.), metabolite profiling, β-amylase activity, β-glucan content, water stress
Citation: Wu X, Cai K, Zhang G and Zeng F (2017) Metabolite Profiling of Barley Grains Subjected to Water Stress: To Explain the Genotypic Difference in Drought-Induced Impacts on Malting Quality. Front. Plant Sci. 8:1547. doi: 10.3389/fpls.2017.01547
Received: 28 April 2017; Accepted: 23 August 2017;
Published: 07 September 2017.
Edited by:Luis A. J. Mur, Aberystwyth University, United Kingdom
Reviewed by:Daniela Romano, University of Catania, Italy
Hui Liu, University of Western Australia, Australia
Copyright © 2017 Wu, Cai, Zhang and Zeng. 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) or licensor 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: Fanrong Zeng, firstname.lastname@example.org