Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Metabolite Adjustments in Drought Tolerant and Sensitive Soybean Genotypes in Response to Water Stress

  • Sonia Silvente ,

    silvente@ccg.unam.mx

    Affiliation Centro de Ciencias Genómicas, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, México

  • Anatoly P. Sobolev,

    Affiliation Istituto di Metodologie Chimiche, Laboratorio di Risonanza Magnetica “Annalaura Segre”, CNR, Monterotondo, Rome, Italy

  • Miguel Lara

    Affiliation Centro de Ciencias Genómicas, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, México

Abstract

Soybean (Glycine max L.) is an important source of protein for human and animal nutrition, as well as a major source of vegetable oil. The soybean crop requires adequate water all through its growth period to attain its yield potential, and the lack of soil moisture at critical stages of growth profoundly impacts the productivity. In this study, utilizing 1H NMR-based metabolite analysis combined with the physiological studies we assessed the effects of short-term water stress on overall growth, nitrogen fixation, ureide and proline dynamics, as well as metabolic changes in drought tolerant (NA5009RG) and sensitive (DM50048) genotypes of soybean in order to elucidate metabolite adjustments in relation to the physiological responses in the nitrogen-fixing plants towards water limitation. The results of our analysis demonstrated critical differences in physiological responses between these two genotypes, and identified the metabolic pathways that are affected by short-term water limitation in soybean plants. Metabolic changes in response to drought conditions highlighted pools of metabolites that play a role in the adjustment of metabolism and physiology of the soybean varieties to meet drought effects.

Introduction

Soybean (Glycine max) is one of the most important grain legumes. It represents not only an essential source of protein, oil and micronutrients in human and animal diets, but is also an attractive crop for the production of biodiesel [1]. Soybean growth is affected by unfavorable environmental factors such as extreme temperatures, drought, nutrient deficiency and soil acidity, which form major constraints for soybean crop production.

Soybean plants form root nodule symbiosis with nitrogen-fixing bradyrhizobia, thus rendering the plant independent of N fertilizers. Nodulation and symbiotic nitrogen fixation have long been recognized as being sensitive to environmental stresses, particularly drought [2]. Water stress reduces nitrogen fixation as a result of a decrease in photosynthate supply [3] or a reduction in the O2 flux into the nodule as well as through overloading nodules with nitrogenous compounds [4][7].

Some of the most important responses of a plant against drought stress are associated with the accumulation of minerals [8] and the enhanced synthesis of osmoprotectants, or compatible solutes, which are part of normal metabolism. The accumulation of these compounds helps the stressed cells in water retention [9] and in the maintenance of the structural integrity of the cell membranes [10].

The types of osmoprotectant metabolites and their relative contribution in lowering the osmotic potential differ greatly among plant species. Osmotic adjustment has been reported in legumes with a high tolerance to water stress [11], [12]. Metabolic adjustments in response to the adverse environmental conditions may highlight pools of metabolites that play important roles in metabolism and physiology and may indicate which pathways have been perturbed by the stress.

Nuclear magnetic resonance (NMR) spectroscopy can be used to monitor and quantify the degree of metabolic impact induced by drought or other environmental disturbances [13], [14], since NMR can bring “high-throughput” spectroscopic/structural information on a wide range of metabolites simultaneously with high analytical precision. One of its main advantages is that it avoids biases against various classes of compounds. Molecular identification is easy and straightforward as it can be deduced from the NMR spectrum of the mixture itself by means of 1D and 2D experiments, standard additions and by comparison with database of standard compounds.

In the present investigation, 1H NMR-based metabolic profiling combined with the physiological studies were conducted in two genotypes of soybean differing in their tolerance to drought in order to elucidate metabolite adjustments in relation to the physiological responses in the nitrogen-fixing plants towards water stress. To our knowledge this is the first report on metabolite profiling in soybean under drought stress. NMR based metabolic profiling approach [15], [16] adopted in this study enabled the identification and quantification of a number of metabolites belonging to various classes of compounds from the crude extracts, without involving any separation step.

thumbnail
Figure 1. Responses of drought tolerant and sensitive varieties of soybean towards short-term water stress.

(a–d) trifoliate leaves and (e–h) nodules. (a, e) dry weight (DW); (b) relative water content (RWC), (f) nitrogenase activity, (c, g) ureide levels and (d, h) proline levels in 17-day-old plants subjected to 10 days water sufficient and deficient conditions. TC: Tolerant genotype Control; TD: Tolerant subjected to Drought, SC: Sensitive genotype Control; SD: Sensitive subjected to Drought. Data represent the mean ± standard deviation (SD) of five to 15 replicates. Different letters at the top of each bar indicate significant differences at P<0.01.

https://doi.org/10.1371/journal.pone.0038554.g001

Results

The two soybean cultivars used in the present study were categorized at CIAP-INTA (Centro de Investigaciones Agropecuarias – Instituto Nacional de Tecnología Agropecuaria, Argentina) as tolerant (NA5009RG) and sensitive (DM50048) to drought stress based on their ability to maintain relative water content (RWC) and growth, and withstand oxidative stress through the modulation of cellular malondialdehide (MDA) levels. The drought tolerant line maintained higher RWC and showed greater ability to withstand oxidative damage owing to lower production of MDA, and exhibited sustained growth at reduced soil moisture conditions (Dr. Celina Luna, personal communication). In the present study, a comparison of these two soybean drought tolerant and sensitive genotypes was undertaken to determine the differences in their metabolic profiles/responses during water stress. Drought was imposed on nodulated soybean plants 17 days after inoculation by withholding water during 10 days, and physiological characteristics such as dry weight (DW), relative water content (RWC), chlorophyll and nitrogenase activity (ARA), and metabolite profiles were analyzed in order to establish the effect of water stress on these plant parameters (Figure 1). The results of the study showed that the water stress produced varied effects on leaf and nodule metabolism in both drought tolerant and sensitive soybean varieties. Under the drought condition imposed, both RWC and DW of the leaves showed a reduction, more remarkably in the sensitive variety (a decrease of about 10% in RWC and 42% in DW in the sensitive genotype as compared to 9% and 15%, respectively, in the tolerant one; Figure 1a, b). In addition, water stress affected proline and ureide contents in leaves of the sensitive variety but not in the tolerant genotype (Figure 1c, d). On the other hand, no effect of drought stress was observed on chlorophyll content as well as chlorophyll a/b ratio in both the genotypes (data not shown).

In case of nodules, in both genotypes, water stress conditions even for 10 days resulted in a dramatic reduction in DW, even though the nitrogenase activity per unit DW remained unaffected (Figure 1e, f). Similar to DW, drought conditions also caused fall in both proline and ureide contents in the tolerant variety (Figure 1g, h). In contrast, in the sensitive genotype, only proline declined while ureide content remained unchanged in the nodules as compared to that in well watered plants.

Comparison of the metabolite profiles of leaf and nodule tissues of tolerant and sensitive varieties subjected to drought

Comparison of the 1H NMR spectra revealed no major qualitative differences in the metabolites between leaf and nodule tissues except for minor aromatic compounds (Figure 2 and Table 1). As an example, a typical 1H NMR spectra of leaf and nodule extracts of well watered tolerant soybean plants together with the assignment of the most abundant metabolites (amino acids, sugars, organic acids) are shown in the Figures 2a and 2b, respectively. Among the resonances in the 9.5 – 7 ppm region, characteristic for aromatic and heteroaromatic compounds, only trigonelline was assigned in leaves extracts. On the other hand, asparagine was detected exclusively in nodules.

thumbnail
Figure 2. 1H NMR spectra of (a) leaves and (b) nodules of water-soluble extracts from well-watered tolerant soybean plants.

Assignments: 1, alanine; 2, GABA; 3, glutamine; 4, malic acid; 5, succinic acid; 6, citric acid; 7, aspartate; 8, asparagine; 9, myo-inositol; 10, choline; 11, pinitol; 12, sucrose; 13, fumaric acid; 14, trigonelline. HDO: deuterated water.

https://doi.org/10.1371/journal.pone.0038554.g002

thumbnail
Table 1. List of variables used in statistical analysis for leaves and nodules samples.

https://doi.org/10.1371/journal.pone.0038554.t001

Basically, regardless of the genotype (tolerant or sensitive plants) or watering conditions, all the 1H NMR spectra of extracts from the same type of tissues (leaves or nodules) share the same signals, although their relative intensity is variable. The intensity of selected signals (Table 1 and Materials and Methods section) was used to calculate the relative molecular abundance of about 15 assigned metabolites. On the other hand, the assignment of minor components was hindered by the scarcity of data on these metabolic compounds in literature. Although the number of compounds identified by NMR is limited, the NMR spectra indeed gave a good picture of what really is present in the plant extracts examined. Results of the study demonstrated that water stress induces several changes in various metabolic pathways in both genotypes; the effect being more pronounced in the leaves than in nodules (Figure 3 and Tables 2,3,4).

thumbnail
Figure 3. Schematic representation of the selected metabolic pathways affected by drought in two soybean genotypes contrasting in sensitivity/tolerance to water stress.

Histograms represent relative changes in the level of the metabolites (arbitrary units) in trifoliate leaves and nodules in the plants subjected to water stress. Values are presented as the mean ± standard deviation (SD) of nine independent biological determinations. TC: Tolerant control; TD: Tolerant Drought, SC: Sensitive Control; SD: Sensitive Drought.

https://doi.org/10.1371/journal.pone.0038554.g003

Under drought conditions, levels of the individual sugars varied considerably among the genotypes: for example, sucrose and myo-inositol levels in the leaves decreased drastically in the sensitive genotype, but no significant changes were observed in the tolerant variety. In contrast, in the leaves of the sensitive genotype, pinitol levels increased under drought while it decreased in the tolerant one. In nodules, however, sucrose content decreased in drought in both varieties while pinitol levels increased. Myo-inositol content, on the other hand, did not alter in the nodules of the both varieties when the water stress was imposed.

Individual organic acids that mainly contributed to the differences in total organic acids under drought were 2-oxoglutaric acid, succinic acid and malic acid. Of these three, only succinic acid levels rose while malic acid content decreased in the leaves of drought stressed plants as compared to well-watered plants in both the genotypes, with no significant changes in nodules. 2-oxoglutaric acid, on the other hand, showed downward trend only in the leaves of sensitive variety.

With regard to the free amino acids in the leaves under drought, the contents of alanine and glutamine decreased in both the genotypes. On the other hand, GABA declined only in the tolerant one, whereas aspartate levels increased in the sensitive genotype. In contrast, no significant differences were observed in the amino acid contents in the nodules of both cultivars under control and stress conditions (Figure 3 and Table 4).

thumbnail
Table 2. Two factors ANOVA with a 2×2 between groups design (drought tolerant vs. sensitive plants, well-watered vs. drought stressed plants), leaves samples.

https://doi.org/10.1371/journal.pone.0038554.t002

thumbnail
Table 3. ANOVA on single groups, two types of grouping (control vs stressed, and tollerant vs sensitive) leave samples.

https://doi.org/10.1371/journal.pone.0038554.t003

thumbnail
Table 4. Two factors ANOVA with a 2×2 between groups design (drought tolerant vs. sensitive plants, well-watered vs. drought stressed plants), nodules samples.

https://doi.org/10.1371/journal.pone.0038554.t004

Principal component analysis (PCA) is one of the most popular explorative methods used to reduce multivariate data complexity. This is a method of choice for identifying patterns, and expressing data in ways that highlight similarities and differences between samples [17]. In our study, PCA was applied to 1H NMR spectral data of control and stressed leaves and nodules derived from two soybeans genotypes with varied tolerance to drought, in order to authenticate the differences between the metabolic profiles of the control and stressed tissues statistically and to identify the main metabolites responsible for the differences.

A scores scatter plot of the first two PCs obtained considering all 1H NMR data derived from the leaves shows a good separation of all four groups (LTC: leaf tolerant control, LTD: leaf tolerant drought, LSC: leaf sensitive control and LSD: leaf sensitive drought) along PC1 axis (Figure 4a). It seems that this separation is due to the treatment (control vs stressed) with further separation between stressed sensitive and stressed tolerant genotypes. The greatest separation along PC1 is between LSC and LSD groups whereas the separation between LTC and LTD is less apparent along PC1 axis, but noticeable along PC2 (see arrows on Figure 4a). This behaviour of the data evidences a markedly more profound metabolic impact of drought stress on sensitive plants with respect to tolerant ones. The separation between leaf samples of well watered and stressed plants along PC1 axis seemed to be mainly attributable to aspartate, succinic acid, sucrose, malic acid, alanine, GABA, myo-inositol and 2-oxoglutaric acid as shown in the complementary PCA loading plot (Figure 4b). In case of the tolerant genotype, leaf samples of control and stressed plants are well separated along PC2 due to the metabolites pinitol, citric acid, choline, and allantoin.

A comparison of metabolite mean levels between LSC and LSD samples (sensitive plants) and between LTC and LTD ones (tolerant plants) was performed using ANOVA (Table 3). ANOVA results confirm the observations obtained with PCA. In fact, the levels of 11 out of 15 metabolites was significantly changed in sensitive plants upon the application of drought stress, whereas lesser number of metabolites (7 out of 15) were influenced by the stress in tolerant plants. Considering drought treatment and genotype as two independant factors and possible interaction between them, ANOVA has been applied using a 2×2 between group design (Table 2). This approach was aimed to give a statistical measure of significance for each factor and interaction between them for each variable. The criterion of statistically significant difference between the mean values was p-level less than 0.01. For 8 variables (alanine, aspartate, GABA, 2-oxoglutaric acid, myo-inositol, pinitol, sucrose, allantoin) out of 15 the interaction between two factors was found to be statistically significant, evidencing different responses of sensitive and tolerant plants to the drought stress on the molecular level. In fact, for aspartate, 2-oxoglutaric acid, myo-inositol, pinitol, sucrose and allantoin the trends of changes upon the application of drought stress is opposite in tolerant and sensitive plants. For example, we can see that in tolerant plant samples, the level of pinitol is higher in control than in water stressed plants, while this is reversed when the sensitive plants were subjected to water stress. It is seems that the idea of considering treatment and genotype as independant factors is not adequate, as the stress produces different results in tolerant and sensitive samples.

In nodules, the PCA analysis (Figure 5a) showed that the first two PCs represented 48.2% of the initial variability contained in the original data. The scores plot exhibited separation between all four groups (NSC: Nodule Sensitive Control, NSD: Nodule Sensitive Drought, NTC: Nodule Tolerant Control and NTD: Nodule Tolerant Drought) when PC2 and PC1 were used as variables. It seems that with a few exceptions, the samples of tolerant plants are separated from sensitive ones along PC1, while control samples are separated from stressed ones along PC2. Plot of loadings (Figure 5b) show the variables responsible for this separation. The metabolites sucrose, aspartate, glutamine, GABA, allantoin, and succinic acid play a crucial role in the separation of tolerant from sensitive samples. On the other hand, the separation of controls from the stressed samples is due to the variations in the levels of malic acid, 2-oxoglutaric acid, fumaric acid, and sucrose.

thumbnail
Figure 4. Principal component analysis (PCA) of 15 metabolites in the leaves from the plants grown under water sufficient and deficient conditions.

Score (a) and loading plot (b) of soybean leaf samples. LTC: Leaf Tolerant Control, LTD: Leaf Tolerant Drought, LSC: Leaf Sensitive Control, LSD: Leaf Sensitive Drought.

https://doi.org/10.1371/journal.pone.0038554.g004

thumbnail
Figure 5. Principal component analysis (PCA) of 15 metabolites in nodules from the plants subjected to water sufficient and deficient conditions.

Score (a) and loading plot (b) of soybean nodule samples. NTC: Nodule Tolerant Control, NTD: Nodule Tolerant Drought, NSC: Nodule Sensitive Control, NSD: Nodule Sensitive Drought.

https://doi.org/10.1371/journal.pone.0038554.g005

The ANOVA analysis (Table 4) confirmed the same variables (that were identified by PCA) as statistically significant for the separation of groups. In addition to this, ANOVA revealed that the level of pinitol and citric acid is significantly different in control and stressed nodule tissues. It is notewothy, that only in the case of two metabolites (pinitol and 2-oxoglutaric acid) the interaction between genotype and drought treatment was significant.

Discussion

In plants, the level of tolerance or sensitivity to water stress depends on the species and genotype, length and severity of water loss, as well as on the developmental stage. In Aeluropus lagopoide, Mohsenzadeh S, et al [18] found significant correlation between leaf relative water content (RWC) and relative growth rate, net photosynthesis rate, chlorophyll and proline contents. The leaf RWC directly reflects the water status of plants and it can be used to identify the genotypes tolerant to stress [19]. In our study, the highest level of RWC was found in the cv. NA5009RG (Figure 1b), and according to Rampino et al. [19] this cultivar could be defined as a genotype tolerant to water deficit. In addition, even under water limitation, this variety showed only a marginal reduction in leaf RWC as compared to the sensitive genotype (DM50048) where RWC was significantly affected. Moreover, our results showed a notable reduction in leaf DW only in the sensitive genotype (Figure 1b-a). The decrease in DW under drought in the sensitive soybean variety might be related to the depletion of sucrose in the leaves of this genotype (Figure 3). These results are in agreement with the findings of Reddy et al. [20] who reported that water stress inhibited dry mater production due to limitation of photosynthesis. Differences in the reduction in leaf DW under water stress has also been reported for different legume species: 78% in mungbean, 60 % in cowpea and 37% in peanut compared with the unstressed plants [21].

Despite of the reduction in the RWC and dry weight of the leaves of the sensitive genotype under drought, no parallel decrease in chlorophyll content was observed (data not shown). This result is in agreement with Ashraf and Iram [12] who also observed a lack of effect of drought on chlorophyll content, and suggested that it could be due to the mild moisture stress to which the experimental plants were exposed.

Enhanced tolerance of plants to low water availability is attributed to the accumulation of soluble sugars in water-stressed tissues [22], acting as osmoprotectants [23], [24]. In contrast to this, the results of the present study with the twenty seven-day old plants showed no enhanced accumulation of soluble sugars such as sucrose and myo-inositol in the leaves of both genotypes (Figure 3; indeed sugar content decreased in the sensitive variety), indicating that these sugars do not play an osmoprotectant role at least at the early stages of the plant growth. Similar findings have been reported for other legume species exposed to osmotic stress [25].

In legumes, pinitol is a common sugar alcohol and it has been described as a common osmoprotectant [11], [26]. Our results showed that the tolerant genotype has higher amounts of pinitol even under normal conditions as compared to the sensitive variety. On the other hand, pinitol synthesis was found to be enhanced in the sensitive genotype under water stress. Future work can only establish the exact role of pinitol in osmoprotection of soybean plants.

Accumulation of amino acids was suggested to aid stress tolerance in plants, through osmotic adjustment, detoxification of reactive oxygen species and by intracellular pH regulation [27], [28]. An equivalent role for most of the amino acids, detected in our experiments, seems unlikely because their content, with the exception of aspartic acid, tended to decrease or remain constant even under drought in both genotypes (Figure 3). These results are in agreement with the findings in Phaseolus vulgaris [29].

Proline, an imino acid, is widely regarded as a main osmoprotectant in water stress tolerance in plants. In the present study with the plants at the vegetative stage, we found that water stress does not trigger enhanced proline synthesis (Figure 1d, h). In contrast, when the drought treatment was imposed during the flowering stage, there was a considerable increase in proline levels in the leaves and nodules of both genotypes (data not shown). Increases in proline level were also observed in others varieties of soybean, when drought was imposed at the reproductive stage and the RWC was lower than the values observed in the present work [30], [31]. Fukutoku and Yamato [32] reported that in intact soybean leaves remarkable proline accumulation occurred only when water stress became severe and protein metabolism was disturbed. These results suggest that even under water stress, the stage of the plant and RWC of the leaves seem to be critical in promoting proline synthesis. Role of proline in stress tolerance remains controversial as some authors have reported high proline levels in the susceptible cultivars subjected to stress conditions [33], [34], while the others have observed the opposite trend [35]. It has been suggested that proline functions as an indicator of plant water status but not a measure of level of tolerance [36].

In both soybean cultivars, drought did not trigger the accumulation of organic acids excepting succinate as its concentration was doubled in the leaves of both genotypes after the imposition of drought. Sassi et al. [29] reported a decrease in the total amount of organic acids in the leaves of a sensitive line of bean plants subjected to osmotic stress.

Present study demonstrated that water stress has varied effects on the metabolic processes of leaves and nodules in the soybean cultivars tested; Nodule DW was more affected than the leaf, and in the nodules of both genotypes the decrease in DW mirrored the drought-induced decline in sucrose. However, drought didn't affect the nitrogenase activity. The lack of response of nitrogenase activity to drought in these soybeans cultivars is in contrast to earlier reports [37], [38], [7], and it perhaps reflects the level of intensity of the stress imposed in various studies. In several soybean cultivars, inhibition of nitrogen fixation under drought stress has been attributed to ureide accumulation in leaves [39]. In the present investigation, ureide accumulation in response to drought stress was observed only in the leaves of the sensitive genotype, but not in the tolerant variety (Figure 1). Since the nitrogenase activity was unaffected under drought stress in both the genotypes, it seems that ureide accumulation in the leaves does not have a feedback inhibitory effect on nitrogen fixation. Alamillo et al. [40] suggested that ureide accumulation and nitrogen fixation follow different kinetics and are probably not causally related.

In nodules, malic acid is the most abundant organic acid and is the main carbon substrate for bacteroid respiration and nitrogen fixation activity. It had been suggested that a decrease in nodule malic acid content under certain environmental conditions may lead to the inhibition of nitrogen fixation [3]. Also in the present work, genotype differences in malic acid content or interactions between genotypes and drought were found to be absent (Figure 3), a situation that is consistent with the results of ARA activity under stress condition (Figure 1).

The analysis of metabolites contributes to the understanding of stress biology of plants through the identification of the compounds and the part they play in acclimation or tolerance response. In the present study, metabolite fingerprinting and profiling based on 1H NMR spectra were used to analyze the similarities and differences among leaf and nodule samples obtained from two soybean genotypes with the aim of identifying markers useful for pinpointing water stress response. In this context the results of our study point to six metabolites in leaves (aspartate, 2-oxoglutaric acid, myo-inositol, pinitol, sucrose, allantoin) and two in nodule (2-oxoglutaric acid and pinitol) that were affected differentially in the genotypes when drought was imposed at the vegetative stage in the nodulated soybean plants. These data provide information that may, with further experimentation, allow elucidation of biochemical pathway underlying stress tolerance in soybean.

The results of the study demonstrated that a combination of 1H NMR and multivariate analyses allows comparisons of overall metabolite fingerprints and that this technique can be applied to conclusively identify differences that are due to stress or genotype. The differences under stress conditions between the two genotypes discussed above are reflected in the PCA models of metabolite content as well. PCA of the present study clearly demonstrated that the major variability in metabolites levels (associated with PC1 in PCA) is due to treatment (control vs stressed) in the case of leaves, while in the case of nodules the major variability is due to genetic makeup (tolerant vs sensitive). The phenomena observed in the case of nodules likely depends on the cumulative effects of plant-bacterial genotypes, specificity of interaction and the resultant symbiosis which in turn alters the metabolism outcome.

Materials and Methods

Plant materials and growth conditions

Seeds of the soybean genotypes (Glycine max L. Merr.), namely NA5009RG (drought tolerant) and DM50048 (drought sensitive) were obtained from CIAP-INTA (Centro de Investigación Agropecuaria-Instituto Nacional de Tecnología Agropecuaria), Argentina. Seeds were surface sterilized with 20% (v/v) commercial bleach, washed extensively with sterile-distilled water, and germinated on a sterile moist filter paper at 28°C, in the dark, for 3 days. Subsequently, the seedlings were transferred to vermiculite in pots and inoculated with 1 cm3 of Bradyrhizobium japonicum USDA110, and grown in a greenhouse at 26°C/19°C (day/night temperatures). Plants received nitrogen-free Summerfield nutrient solution [41] twice a week until the stress treatments were imposed.

Drought stress was imposed on 17-day old plants at the vegetative stage by withholding water supply for 10 days until soil water content reached 23% (0.230 g H2O g−1 dry soil). A set of well-watered plants served as a control. Measurements of dry weight, relative water content and nitrogen fixation, as well as sampling of the plant tissues (nodules and first trifoliate leaves) for chlorophyll, proline, ureide and metabolome analysis were performed at the end of the stress period.

The results presented are the means with standard deviations of five to 15 replicates. All data obtained was subjected to one-way analysis of variance (ANOVA), and the mean differences were compared by lowest standard deviations (LSD) test using the STATISTICA package for Windows (version 5.1, 1997). Comparisons with P values<0.01 were considered significantly different.

Leaf relative water content (RWC)

RWC was measured according to Barrs and Weatherly [42]. Briefly, immediately after sampling, the leaves were weighted and then soaked overnight in distilled water at 4°C. After the cold incubation, the leaves were blotted dry and weighed prior to oven-drying at 80°C for 48 h. Subsequently, dry weight of the plant samples was determined. The leaf relative water content was calculated using the following formula: RWC  =  (FW – DW)/(TW-DW)) ×100, where FW is fresh weight, DW is dry weight, and TW is turgid weight (weight after the leaf was kept in distilled water for overnight).

Nitrogen fixation

Nitrogenase activity was determined by acetylene reduction assay (ARA), [43].

Estimation of ureide content

Concentration of ureides present in leaf and nodule cell-free extracts was measured using the colorimetric detection method of Vogels and Van Der Drift [44]. Allantoic acid dissolved at a concentration of 10 mM in water served as a standard for ureide estimation.

Determination of proline

Samples of fresh plant tissues (0.5 g) were homogenized in 5 ml of 3% aqueous sulphosalicyclic acid and supernatant was collected after centrifugation. Two mL extract was reacted with 2 mL acid-nihydrin and 2 mL glacial acetic acid, and incubated for 1 h in a boiling water bath. The reaction was terminated by placing the test tubes on ice after which the reaction mixture was vigorously mixed with 2 ml of toluene. After warming to 25°C, proline present in the upper toluene layer was measured at 520 nm [45].

Metabolite extraction and NMR analysis

Comparative metabolite profiling was performed in leaves and nodules derived from the drought tolerant and sensitive soybean genotypes subjected to water stress. First trifoliate expanded leaves and nodules were collected after 10 days of drought treatment, and immediately frozen in liquid nitrogen, lyophilized and stored at −80°C till they were subjected to NMR analysis.

Water-soluble extracts were derived from 20 mg of lyophilized tissues mixed with 0.9 mL of CH3CN:H2O (1∶1 v/v); Extracts were clarified by centrifugation at 10,000 rpm (5,600 g) for 7 min, and the supernatant obtained was filtered and lyophilized. The dry residue was dissolved in 0.7 mL of 400 mM D2O phosphate buffer (pD = 6.5) containing 1.0 mM of 3-(trimethylsilyl) propionic−2,2,3,3-d4 acid sodium salt (TSPA) and transferred into a standard 5 mm NMR tube. NMR spectra of the extracts were recorded at 300 K on a Bruker AVANCE AQS600 spectrometer operating at the proton frequency of 600.13 MHz and equipped with a Bruker multinuclear z-gradient inverse probe-head capable of producing gradients in the z-direction with the strength of 55.4 G/cm.

Proton spectra were referenced to the signals of TSPA methyl group at δ = 0.00 ppm in D2O phosphate buffer. The 1H spectra of the aqueous extracts were acquired by co-adding 512 transients with a recycle delay of 2.5 s and 32 K data points (acquisition time 40 min). The residual HDO signal was suppressed using a presaturation during the relaxation delay with a long single soft pulse. To avoid possible saturation effects, the experiment was carried out by using a 45° flip angle pulse of 8.0 μs 2D NMR experiments, namely 1H–1H TOCSY and 1H–13C HSQC, were performed using the same experimental conditions as previously reported by Sobolev et al [46].

Fifteen metabolites in leaves and nodules extracts were identified and used for statistical analysis, see Tables 1 and 2. Metabolites were assigned and identified using 2D experiments 1H–1H TOCSY and 1H–13C HSQC and by comparison with the literature data [47].

The signal heights of selected 1H resonances of water-soluble metabolites (Table 1) were measured with respect to the height of TSPA signal used as internal standard. The height of TSPA signal was normalized to 100. The obtained values (relative molecular abundances of selected metabolites) were used in statistical analysis.

Statistical Analysis of NMR data

The statistical treatment of the NMR data was performed using the STATISTICA package for Windows (version 5.1, 1997). Two factors ANOVA has been performed with a 2 × 2 between groups design (drought tolerant vs. sensitive plants, well-watered vs. drought stressed plants). Principal component analysis (PCA) was performed using all 15 variables for leaves and nodules. Before the PCA analysis the variables were mean-centered and each variable was divided by its standard deviation (autoscaling). The effects and interactions represented in bold in the Tables 2,3,4 were statistically significant within the 99% confidence interval.

Acknowledgments

We would like to thank PM Reddy for critically reading the manuscript, Nestor Mariano for his valuable help in statistical analysis and Celina Luna (CIAP-INTA (Centro de Investigaciones Agropecuarias – Instituto Nacional de Tecnología Agropecuaria)) for his kind gift of soybean seeds.

Author Contributions

Conceived and designed the experiments: SS APS. Performed the experiments: SS APS. Analyzed the data: SS APS. Contributed reagents/materials/analysis tools: SS APS. Wrote the paper: SS APS ML.

References

  1. 1. Pimentel D, Patzek TW (2008) Ethanol production using corn, switch grass and wood; biodiesel production using soybean. In Pimentel, D. (ed.), Biofuels, Solar and Wind as Renewable Energy Systems: Benefits and Risks. Springer, Dordrecht, The Netherlands, 375–396:
  2. 2. Sprent JI (1972) The effects of water stress on nitrogen fixing root nodules. IV Effects on whole plants of Vicia faba and Glicine max. New Phytol 71: 603–611.
  3. 3. Arrese-Igor C, Gonzalez EM, Gordon AJ, Minchin FR, Galvez L, et al. (1999) Sucrose synthatse and nodule nitrogen fixation under drought and other environmental stresses. Symbiosis 27: 189–212.
  4. 4. Neo HH, Layzell DB (1997) Phloem glutamine and the regulation of O2 diffusion in legume nodules. Plant Physiol 113: 259–267.
  5. 5. Serraj R, Vadez V, Sinclair TR (2001) Feedback regulation of symbiotic N2 fixation under drought stress. Agronomie 21: 621–626.
  6. 6. King CA, Purcell LC (2005) Inhibition of N2 fixation in soybean is associated with elevated ureides and amino acids. Plant Physiol 137: 1389–1396.
  7. 7. Marino D, Frendo P, Ladrera R, Zabalza A, Puppo A, et al. (2007) Nitrogen fixation control under drought stress. Localized or systemic? Plant Physiol 143: 1968–1974.
  8. 8. Samarah N, Mullen R, Cianzio S (2004) Size distribution and mineral nutrients of soybean seeds in response to drought stress. J Plant Nutr 275: 815–835.
  9. 9. Hare PD, Cress WA, Van Staden J (1998) Dissecting the roles of osmolyte accumulation during stress. Plant Cell Environ 21: 535–553.
  10. 10. Conroy MJ, Goldsberry JR, Hines JE, Stotts DB (1988) Evaluation of aerial transects surveys for wintering American black ducks. Journal of Wildlife Management 52: 694–703.
  11. 11. Ford CW (1984) Accumulation of low molecular weight solutes in water stressed tropical legumes. Phytochem 23: 1007–1015.
  12. 12. Ashraf MY, Iram A (2005) Drought stress induced changes in some organic substances in nodules and other plant parts of two potential legumes differing in salt tolerance. Flora 200: 535–546.
  13. 13. Bligny R, Douce R (2001) NMR and plant metabolism. Curr Opin Plant Biol 4: 191–196.
  14. 14. Charlton AJ, Donarski JA, Harrison M, Jones SA, Godward J, et al. (2008) Responses of the pea (Pisum sativum L.) leaf metabolome to drought stress assessed by nuclear magnetic resonance spectroscopy. Metabolomics 4: 312–327.
  15. 15. Fiehn O (2002) Metabolomics – the link between genotypes and phenotypes. Plant Mol Biol 48: 155–171.
  16. 16. Colquhoun IJ (2007) Use of NMR for metabolic profiling in plant systems. J Pestic Sci 32: 200–212.
  17. 17. Goodacre R, Shann B, Gilbert RJ, Timmins EM, McGovern AC, et al. (2000) Detection of the dipicolinic acid biomarker in Bacillus spores using Curie-point pyrolysis mass spectrometry and Fourier transform infrared spectroscopy. Anal Chem 72: 119–127.
  18. 18. Mohsenzadeh S, Malboobi MA, Razavi K, Farrahi-Aschtiani S Razavi K (2006) Physiological and molecular responses of Aeluropus lagopoides (Poaceae) to water deficit. Environ Exp Bot 56: 314–322.
  19. 19. Rampino P, Pataleo S, Gerardi C, Perotta C (2006) Drought stress responses in wheat: Physiological and molecular analysis of resistant and sensitive genotypes. Plant Cell Envir 29: 2143–2152.
  20. 20. Reddy AR, Chaitanya KV, Vivekanandan M (2004) Drought-induced responses of photosynthesis and antioxidant metabolism in higher plants. J Plant Physiol 161: 1189–1202.
  21. 21. Pandey RK, Herrera WAT, Villepas AN, Pendelton JW (1984) Drought response grain legumes under irrigation gradient. III. Plant Growth Agron J 76: 557–560.
  22. 22. McManus MT, Bieleski RL, Caradus JR, Barker DJ (2000) Pinitol accumulation in mature leaves of white clover in response to a water deficit. Env Exp Botany 43: 11–18.
  23. 23. Ingram J, Bartels D (1996) The molecular basis of dehydration tolerance in plants. Annu Rev Plant Physiol Plant Mol Biol 47: 377–403.
  24. 24. Sanchez FJ, Manzanares M, Andres EF, Tenorio JL, Ayerbe L (1998) Turgor maintenance, osmotic adjustment and soluble sugar and proline accumulation in 49 pea cultivars in response to water stress. Field Crops Res 59: 225–235.
  25. 25. Pinheiro C, Passarinho JA, Ricardo CP (2004) Effect of drought and rewatering on the metabolism of Lupinus albus organs. J Plant Physiol 161: 1203–1210.
  26. 26. Keller R, Ludlow M (1993) Carbohydrate metabolism in drought-stressed leaves of pigeonpea (Cajanus cajan). J Exp Bot 44: 1351–1359.
  27. 27. De Ronde JA, Van Der Mescht A, Steyn HSF (2000) Proline accumulation in response to drought and heat stress in cotton. Afr Crop Sci J 8: 85–92.
  28. 28. Alia MP, Matysik J (2001) Effect of proline on the production of singlet oxygen. Amino Acids 21: 195–200.
  29. 29. Sassi S, Aydi S, Hessini K, Gonzalez EM, Arrese-Igor C, et al. (2010) Long-term mannitol-induced osmotic stress leads to stomatal closure, carbohydrate accumulation and changes in leaf elasticity in Phaseolus vulgaris leaves. Afr J Biotechnol 9: 6061–6069.
  30. 30. Angra S, Kaur S, Singh K, Pathania D, Kaur N, et al. (2010) Water-deficit stress during seed filling in contrasting soybean genotypes: Association of stress sensitivity with profiles of osmolytes and antioxidants. Int Agric Res 5: 328–345.
  31. 31. Masoumi H, Darvish F, Daneshian J, Nourmohammadi G, Habibi D (2011) Chemical and biochemical responses of soybean (Glycine max L.) cultivars to water deficit stress. Aust J Crop Sci 5: 544–553.
  32. 32. Fukutoku Y, Yamada Y (1981) Sources of Proline-nitrogen in Water-stressed Soybean (Glycine max L.) I. Protein Metabolism and Proline Accumulation. Plant Cell Physiol 22: 1387–1404.
  33. 33. Premachandra GS, Hahn DT, Rhodes D, Joly RJ (1995) Leaf water relations and solute accumulations in two grain sorghum lines exhibiting contrasting drought tolerance. J Exp Bot 46: 1833–1841.
  34. 34. Sundaresan S, Sudhakaran PR (1995) Water stress-induced alternations in the proline metabolism of drought-susceptible and tolerant cassava (Manihot esculenta) cultivars. Physiol Plant 94: 635–642.
  35. 35. Hien DT, Jacobs M, Angenon G, Hermans C, Thu TT, et al. (2003) Proline accumulation and Δ1-pyrroline-5-carboxylase syhthetase gene properties in three rice cultivars differing in salinity and drought tolerance. Plant Sci 165: 1059–1068.
  36. 36. Lazcano-Ferrat I, Lovatt CJ (1999) Relationship between relative water content, nitrogen pools, and growth of Phaseolus vulgaris L. and P. acutifolius A. Gray during water deficit. Crop Sci 39: 467–475.
  37. 37. González EM, Aparicio-Tejo PM, Gordon AJ, Minchin FR, Royuela M, et al. (1998) Water-deficit effects on carbon and nitrogen metabolism of pea nodules. J Exp Bot 49: 1705–1714.
  38. 38. Gálvez L, Gónzalez EM, Arrese-Igor C (2005) Evidence for carbon flux shortage and strong carbon/nitrogen interactions in pea nodules at early stages of water stress. J Exp Bot 56: 2551–2561.
  39. 39. Serraj R, Sinclair TR, Purcell LC (1999) Symbiotic N2 fixation response to drought. J Exp Bot 50: 143–155.
  40. 40. Alamillo JM, Díaz-Leal JL, Sánchez-Moran MV, Pineda M (2010) Molecular analysis of ureide accumulation under drought stress in Phaseolus vulgaris L. Plant Cell Environ 33: 1828–1837.
  41. 41. Summerfield RJ, Huxley PA, Minchin F (1977) Plant husbandry and management techniques for growing grain legumes under simulated tropical conditions in controlled environments. Exp Agr 13: 81–92.
  42. 42. Barrs HD, Weatherley PE (1962) A re-examination of the relative turgidity technique for estimating water deficits in leaves. Aust J Biol Sci 15: 413–428.
  43. 43. Hardy RWF, Holsten RD, Jackson EK, Burns RC (1968) The acetylene-ethylene assay for N2 fixation: Laboratory and field evaluation. Plant Physiol 43: 1185–1207.
  44. 44. Vogels GD, Van Der Drift C (1970) Differential analyses of glyoxylate derivatives. Anal Biochem 33: 143–157.
  45. 45. Bates LS, Waldren RP, Teare ID (1973) Rapid determination of free proline for water-stress studies. Plant Soil 39: 205–207.
  46. 46. Sobolev AP, Segre A, Lamanna R (2003) Proton high-field NMR study of tomato juice. Magn Reson Chem 41: 237–245.
  47. 47. Fan TWM (1996) Metabolite profiling by one- and two-dimensional NMR analysis of complex mixtures. Prog Nucl Magn Reson Spectrosc 28: 161–219.