Abstract

The production of soybean is restricted in sub-Saharan Africa by several stress conditions, including drought because its production is exclusively rain-fed. Identifying drought resistant varieties is of paramount importance. Thus, the objectives of this work were to (i) evaluate the effect of polyethylene glycol 6000 (PEG6000) on soybean at the seedling stage, (ii) determine the root system architecture and physiological characters to water deficit stress, and (iii) establish the correlation among the quantitative variables responsible for drought tolerance in soybean varieties. Twenty soybean accessions (G1 to G20) were subjected to 10% PEG6000 concentration at seedling stages under a controlled environment using a randomized complete block design with 3 replicates. Vegetative growth data were collected. Highly significant differences of proline, carotenoid, chlorophyll a, and chlorophyll b contents were recorded among the 20 accessions in response to PEG application. G16 and G19 had the highest carotenoid, highest chlorophyll a, and chlorophyll b. The highest dry weight was observed in G16 and G10, while the number of leaves was recorded in G19 and G17. G4, G9, G10, and G13 demonstrated the highest dry weight. The PEG-simulated drought stress reduced the average root diameters and the number of lateral roots of all 20 accession plants. G1, G3, G4, G8, G9, and G15 had the longest roots than the control plants as a mechanism to withstand drought stress by seeking water in the deep. Number of leaves was significantly and positively correlated with shoot dry weight, root dry weight, and root diameter but was significantly and negatively correlated with canopy wilting. Proline content was significantly and positively correlated with carotenoid, chlorophyll content, chlorophyll a, and chlorophyll b. G10, G19, G9, G6, G16, G17, G20, G16, and G18 are the tolerant cultivars to drought stress on the basis of growth, physiological, and root system architecture.

1. Introduction

Soybean is an important crop worldwide due to its protein, oil, fatty acids, fibre, and nutrient content essential for health in humans and animals [13]. Processed soybean products include soy cheese, soy cake, soy milk, soy mustard, soy yogurt, and infant foods associated with high protein, which is cheaper when compared with fish and meat [46]. However, the production of soybean is limited in sub-Saharan Africa by numerous stresses, including drought because its production is exclusively rain-fed. Drought conditions often reduce yield at multiple stages of development, including flowering, pollination, and the grain-filling stage [79]. The decline in yields due to drought calls for identification of drought tolerance crops to meet up with the demand of crops including soybeans in the global markets [10]. There are tolerant soybean varieties among the existing genotypes in the world gene banks. Exploiting the existing collection of soybean germplasm will help come up with soybeans resistant to drought conditions. Therefore, there is an imperative need to screen important crops such as soybean for tolerance to drought stress so as to enhance its productivity and yields. Such results would substantially benefit many small farmers, whose life is challenged by constant food insecurity and hunger. IPCC [11, 12] reported that the poorest countries will suffer the most from the consequences of climate change because abiotic stress will reduce yields drastically in the most economically important crops in tropical and subtropical regions, and food production in Africa will severely be compromised. Three main mechanisms are employed by plants to deal with water deficit stress, and these include drought escape and avoidance and tolerance to drought and use various mechanisms to cope with drought stress [13]. In drought escape, plants shorten and complete their life cycle before the commencement of drought stress or complete the critical stages of its development before the onset of water deficit in the soil, which allows the plant to produce some yields instead of total crop failure. Avoidance of drought is the second mechanism, which involves reducing water loss from aerial parts by using cuticle wax or by having the ability to extract moisture from the soil efficiently and the plants continue to keep high water status during periods of stress. The third mechanism is tolerance to drought stress, during periods of drought stress plants continue to maintain turgor by delaying stomatal closure, keeping chloroplast volume, and reducing leaf wilting and plants continue metabolic reactions even at low water potential by synthesizing osmoprotectants, osmolytes, or compatible solutes or by protoplasmic tolerance [14]. This tolerance to internal water deficits allows to prolong the functioning of photosynthesis. Carbon products can then be used for both osmotic adjustment and root growth. Another consequence of maintaining carbon metabolism is a decrease in the frequency of photoinhibition episodes. At the cellular level, osmotic adjustment plays a key role in maintaining turgor at low leaf water potentials.

Methods including mannitol, sorbitol, polyethylene glycol (PEG), and water withholding have been used to impose drought stress in the field or greenhouse or laboratories to understand the mechanisms of plant tolerance to drought stress. PEG-simulated drought stress and water withholding have been the forefront methods of studying drought stress. PEG 6000 is a natural polymer, which is nontoxic and nonionic. PEG 6000 induces osmotic stress in plants by reducing water potential as it is observed during the shortage of rainfall [15, 16].

Some studies have used PEG to simulate drought stress in crops [1722], while the stress study by water withholding is highly reported in the tropics [23, 24]. PEG induces drought in the plant root systems and causes the inaccessibility of available water in the soil to the growing plants. Thus, the present research was conducted to (i) evaluate the effect of PEG on soybean at the seedling stage, (ii) determine the root system architecture and physiological characters due to water deficit, and (iii) establish the correlation among the quantitative variables responsible for drought tolerance in soybean varieties.

2. Materials and Methods

2.1. Study Area and Plant Materials

This study was carried out in the greenhouse at Bowen University, Iwo, Osun State, Nigeria, located between 7°38′Nord and 4°11′East longitude with an altitude of 322 m above sea level.

A total of 20 soybean accessions used in the present study were obtained from the International Institute of Tropical Agriculture (IITA), Ibadan Nigeria (Table 1).

2.2. Experimental Design and Drought Conditions

The study was conducted from June to August 2022 in a randomized complete block design with 3 replicates and twenty soybean accessions. The growing medium for soybean production consisted of topsoil and sawdust in the ratio 2 : 1. Eleven (11) kg of the substrate were filled in each experimental bag. Two factors, namely, accession and PEG6000-simulated drought stress, were studied.

Twenty-one (21) days after sowing, drought treatments were imposed to soybean seedlings thinned to three plants per experimental bag. Two treatments were applied including control and simulation of drought with 10% polyethylene glycol 6000 (10% PEG-6000) for 14 days. 200 ml of 10% PEG was applied to designated seedlings, while 200 ml of water was applied to plants which served as control. Pesticides and fertilizers were not applied during the experiment.

2.3. Data Collection
2.3.1. The Growth Parameters

Measured were plant height (cm), number of leaves, leaf length (cm), leaf width (cm), and above ground wilting, which was recorded using a modified wilting scale of 0–5 [21, 25]. The leaf wilting was rated in stressed plants visually on a scale of 0 to 5. The wilting scale of 0 represented no wilting. 1 represented the unifoliate wilting. 2 represented the 1st wilting of trifoliate leaves. 3 represented the first 2 trifoliolate, 4 represented the first 3 trifoliolate leaves wilting together. 5 represented the whole plant wilting.

2.3.2. Root System Architecture

On the last day of the experiment, the uprooted, treated, and control plants were meticulously washed with tap water to get rid of debris and soil particles. Then, root length, number of lateral roots, and root diameter were measured as morphological root system architecture traits.

2.3.3. Biomass Yield

Biomass yield was measured after 24 hours of drying below and above ground samples in an oven set at 80°C.

2.3.4. Physiological Parameters

On the 14th day of drought imposition, fresh leaves were sampled for the measurement of carotenoid, chlorophyll a, chlorophyll b, and proline content. These photosynthetic pigments were evaluated using a modified Arnon [26] procedure. Carotenoid, chlorophyll a, and chlorophyll b contents were computed using the equation of Porra [27], whereas proline was assessed using Bates et al. [28] procedure. For the measurement of chlorophyll content using SPAD meter, the average of triplicate readings was recorded on the fully expanded leaflet of each plant. SPAD meter is based on the ratio of transmission of near infrared to red wavelengths. The measurement of chlorophyll content using SPAD meter is a nondestructive method unlike in the case of chlorophyll a and chlorophyll b where the destructive method was applied for the determination of chlorophyll a and b using the equations.

2.4. Statistical Analysis

The morphological, root architecture, biomass yield, and physiological data collected were subjected to analyses of variance by using the R statistical package version R-4.0.5. Fischer’s least significant difference (F-LSD) was used for the separation of means at a probability level of 5%. PCA was run using the FactoMineR and factoextra packages, and Pearson correlation was done using corr. Functions in R. A hierarchical cluster analysis was performed using the cluster factoextra package in R.

3. Results

3.1. ANOVA for Genotypes, Peg-Simulated Drought, and Genotypes by Environment Interaction

The Analysis of the Variance indicated that PEG6000-simulated drought stress by variety interaction had a very highly significant effect on root length , carotenoid , proline , and chlorophyll a root diameter (Table 2).

The ANOVA (Table 2) also revealed that variety had a significant effect on root architecture such as the number of lateral roots , root diameter , and root length ; physiological parameters such as carotenoid , proline , chlorophyll a , and chlorophyll b ; and morphological traits including plant height and the number of leaves .

PEG6000-simulated drought stress had each highly significant influence on all the morphological root architecture, physiological, and morphological characters except for the number of leaves where a significant difference was observed. There was no significant difference in chlorophyll content (Table 2).

3.2. Effect of PEG-Simulated Drought Stress on Canopy Wilting

Canopy wilting was recorded in all accessions, although the degree of wilting differed from one accession to the other (Figure 1). The highest wilting was observed in G3 and G10, while the lowest canopy wilting was in G19. No wilting was observed in well-watered plants.

3.3. Effect of PEG-Simulated Drought Stress on Morphological Traits

The analysis of variance of the number of leaves revealed a highly significant difference among varieties (Table 3). Also, there were significant differences in leaf length among the varieties. The average leaf length ranged from 8 to 12 under drought stress, and the highest value was observed in G20. The highest value of plant height was recorded in G6 and G14, and highly significant differences were observed among varieties subjected to simulated drought stress. There was a significant difference between the varieties for leaf width under simulated drought stress. PEG-simulated drought stress inhibited the morphological traits when compared to the control plants.

3.4. Effect of PEG-Simulated Drought Stress on the Morphological Architecture of Root System and Biomass Yield

There were significant differences in the average root length among soybean accessions under normal and drought stress simulated by PEG 6000 (Table 4 and Figure 2). The results showed that there was a reduction of root length in G5, G6, G10, G11, G12, G13, G14, G17, G18, and G19, while G1, G3, G4, G8, G9 and G15 after PEG-6000 simulated drought stress had their root lengths higher than the control plants as a mechanism to seek for water in soil. The PEG-simulated drought stress reduced the average root diameters and the number of lateral roots of all 20 accession plants. Under normal watering conditions and drought stress, there was a significant difference for the above dry matter yield, and the highest value was recorded in G11, G16, and G8 under normal conditions and G16 and G10 under PEG-simulated drought stress (Table 4).

The morphological architecture of root system under (A) control, (B) drought, and (C) PEG-simulated drought stress are revealed in Figure 2. We observed that the main roots of seedling under drought (B) and PEG-simulated drought stress (C) are longer than those of under well-watered plants. The PEG-simulated drought stress has caused the change in the root architecture by extending the main roots and reducing the number of lateral roots.

Figure 3 shows the below dry biomass yield of the 20 soybean accessions. The PEG-simulated drought stress decreased below dry biomass yield in all accessions. The best below dry weight was recorded with G10 followed by G9, G4, G6, and G13.

3.5. Effect of PEG-Simulated Drought Stress on Physiological Traits

Highly significant differences were recorded with accessions regarding carotenoid, chlorophyll a, and chlorophyll b contents in their responses to PEG application (Table 5). There was an increase in carotenoids, chlorophyll a, and chlorophyll B for some varieties and a decrease for others. Increase in carotenoid and chlorophyll a was recorded in G1, G2, G6, G7, G8, G10, G11, G12, G14, G15, G16, G17, and G19 under drought stress. G16 and G19 drought recorded the highest values of 29.46 and 30.13, respectively, for carotenoid; 19.96 and 21.46, respectively, for chlorophyll a; and 10.64 and 9.38, respectively, for chlorophyll b (Table 5). There was no significant difference in the average chlorophyll contents under well-watered plants and water-stressed plants.

Figure 4 presents the concentration of proline under normal and water deficit conditions. A highly significant increase in proline concentration was observed in the leaves of all genotypes except for G1 and G3 under water stress. G6, G2, G16, and G7 accumulated more proline than others.

3.6. Principal Component Analysis

The relationships among physiological, morphological, root architecture system, and biomass yield variables were established using the principal component analysis (Figures 57). A total of 10 axes were generated explaining the total variability (Figure 5). The first two dimensions accounted for 48.14% of the variance culminating in a moderate contribution to the total variation (Figures 5 and 6). But the combination of the first five dimensions explained up to 80.54% of the variance in the data, resulting in a very strong contribution to the total variation (Figures 57). Carotenoid, proline, chlorophyll content, chlorophyll a and chlorophyll b, root diameter, and shoot dry weight significantly contributed to the formation of Dim1, whereas plant height, leaf length, leaf width, root length, and root diameter were significantly and positively correlated with Dim2, and the number of lateral roots was significantly and negatively correlated with Dim2. Canopy wilting and the number of lateral roots showed a significant and positive correlation with Dim3, but the number of leaves was significantly and negatively correlated with Dim3. Root dry weight, shoot dry weight, and the number of leaves were significantly and positively correlated with Dim4 (Figures 5 and 6).

In Figure 6, both variables and accessions were loaded at the same time indicating the relationship among traits and the distances between genotypes. The nearer the vector to the axis, the stronger the correlation. TGm-112, TGm-263, TGm-422, TGm-3972, TGm-4015, and TGm-4502 contributed significantly to Dim1, while TGm-50, TGm-95, TGm-665, TGm-946, TGm-4500, and TGm-4022 significantly contributed to Dim2. TGm-4004 and TGm-4144 are the two significant accessions for the formation of Dim3. Dim4 consisted of TGm-1678, TGm-4400, and TGm-4414, while only TGm-951 significantly contributed to Dim5 (Figure 6).

TGm-1678 (G10) did well for number of leaves, leaf length, root length, above dry weight, below dry weight, chlorophyll a, and chlorophyll b. TGm-4502 (G19) did well for no. of leaves, leaf width, root diameter, carotenoids, chlorophyll a, chlorophyll b, and chlorophyll content. TGm-1326 (G9) did well for leaf length, leaf width, root length, and root diameter. TGm-665 (G6) did well for leaf length, leaf width, root diameter, carotenoids, chlorophyll a, and chlorophyll content. TGm-4015 (G16) did well for root length, number of lateral roots, above dry weight, carotenoids, chlorophyll a, chlorophyll b, and chlorophyll content. TGm-4499 (G17) did well for no. of leaves, above dry weight, carotenoids, chlorophyll a, and chlorophyll b. TGm-4022 (G20) did well for leaf width, plant height, root length, root diameter, carotenoids, and chlorophyll b. G6 did well for leaf length and root length. TGm-951 (G8) did well for leaf length, root length, and chlorophyll a.

3.7. Relationship between Variable Analysis

Proline content was significantly and positively correlated with the following characters: carotenoid (r = 0.67, ), chlorophyll content (r = 0.74, ), chlorophyll a (r = 0.69, ), and chlorophyll b (r = 0.57, ) (Figure 8). There was significant and positive correlation between carotenoid and chlorophyll a (r = 0.97, ), carotenoid and chlorophyll b (r = 0.85, ), and carotenoid and chlorophyll content (r = 0.45, ). Shoot dry weight was positively correlated with proline (r = 0.216, ), root diameter, root length, root dry weight, number of lateral root, leaf length, chlorophyll content, chlorophyll a, chlorophyll b, positively and significantly correlated with (r = 0.448, ). Plant height was significantly and positively correlated with leaf width (r = 0.50, ), root dry weight (r = 0.46, ), and root diameter (r = 0.51, ). The number of leaves was significantly and positively correlated with shoot dry weight (r = 0.45, ), root dry weight (r = 0.47, ), and root diameter (r = 0.51, ) but was significantly and negatively correlated with canopy wilting (r = −0.45, ). Leaf width was significantly and positively correlated with leaf length (r = 0.81, ), root length (r = 0.45, ), and root diameter (r = 0.50, ).

4. Discussion

Soybean is one of the main sources of protein and vegetable oil contributing to the world economy, and its production is faced with the challenges of change in climate, especially being sensitive to drought stress. It is crucial to develop resistant soybean varieties to drought for food security reasons. The first step toward the realization of this is the screening and then the identification of resistant genotypes, which can be used as planting materials or gene reservoirs for soybean breeding programs. There are many mechanisms that plants put in place to lessen the effect of drought stress. Yan et al. [29] suggested that there should be more screening of soybean germplasm for vital characters such as high heritability and stable traits such as drought tolerance in soybean. Thus, this study has been carried out to identify drought-tolerant soybean varieties because it is really a concern for farmers and the world at large to have the yields of soybean compete with other crops such as cereals.

The synthesis of chlorophyll is fundamental to light interception for efficient photosynthesis [30, 31]. Modifying pigmentation could lead to a decrease in photosynthesis and an increase in photooxidation and photoinhibition resulting in more rapid yellowing and death of the leaves [32, 33]. This study showed an increase in carotenoids, chlorophyll a, chlorophyll b, and free proline contents in most of the genotypes, while there was a significant decline in other genotypes under PEG-simulated drought stress. The same observations were made by Wijewardana et al. [34] who reported that carotenoids increased linearly in soybean plants with declining soil moisture content. Basal et al. [22] observed that chlorophyll a and chlorophyll b declined as PEG concentration increased at all stages, but the reduction was insignificant at vegetative stages and significant at reproductive stages. Similar results were obtained by Zhang et al. [35] who observed that soybean plants under water deficit stress showed lower chlorophyll a, chlorophyll b, and carotenoid content. The concentration of free proline decreased in G1 and G3, whereas there was a sharp and remarkable increase in free proline for the other 18 genotypes under water stress. Our results corroborate those obtained with cowpea [36] and beans [37]. The chlorophyll contents were not affected by drought stress after 10 days of drought imposition. Gunes et al. [38] and Silva et al. [40], who worked, respectively, on chickpea cultivars and sugarcane genotypes, considered chlorophyll content, relative water content, and ascorbic acid as secondary indicators for choosing drought-tolerant genotypes and yield as a key characteristic for measuring drought tolerance under water deficit.

PEG-simulated drought stress in this study triggered reductions in plant height, leaf length, leaf width, and the number of leaves. The decline was more pronounced in plant height, leaf length, and leaf width than in the number of leaves. The reasons could be due to their sensitivity to lack of water and the significant decline in metabolic reactions for food production and nutrient transportation, which are essential for plant growth. Similar results were obtained by previous studies [29, 40]. Yan et al. [29] reported that plant height and shoot growth of soybean varieties declined with PEG-simulated drought stress imposition. The abovementioned results are consistent with those of [4143] who showed that significant differences in variety and drought stress for plant height, leaf area index, shoot biomass, root biomass, and total root length.

Morphological root system architecture plays a crucial role in the productivity and survival of plants subjected to either flooding or drought conditions [29, 44, 45]. The root diameters and the number of lateral roots of all 20 accession plants subjected to drought stress were smaller than those of control plants. This is consistent with the results of [29], who recorded a decrease in root diameter but an increase in the number of lateral roots with the PEG-simulated drought stress. We observed variability in root length produced by genotypes as their response to drought stress. The root lengths of G5, G6, G10, G11, G12, G13, G14, G17, G18, and G19 were reduced due to PEG-simulated drought, while under the same drought conditions, G1, G3, G4, G8, G9, and G15 responded differently by increasing their root lengths as a mechanism of tolerance to drought stress. Battisti and Sentelhas [46] reported that soybean varieties that can produce deeper root depth under drought stress are to be considered as resilient to climate change because they can elongate their taproots to search for water deep in the soil.

The PEG-simulated drought stress decreased below dry biomass yield in almost all the accessions. Thus, the below dry biomass yield is more affected than the above dry matter yield by drought stress. This is consistent with the results of [21, 29, 4749] who reported that under water stress, there is a notable decrease in shoot dry biomass and root dry biomass at growth phases.

Plants of the accessions wilted as their response to PEG-simulated drought stress, and the degree of canopy wilting differed from one accession to the other. The same results about soybean wilting under drought stress were observed by previous researchers [21, 29, 50].

Positive and significant correlations among physiological, morphological, and root architecture were observed in this study as shown in the Results section. Liu et al. [51] in their study observed significant and positive correlations between drought resistance and dry weight, root volume, total length, and the number of lateral roots in soybean. We observed the same results with Yan et al. [29] who demonstrated that the wilting index was intensely and negatively correlated with root length, root surface area, number of lateral roots, and root volume.

Under PEG-simulated drought stress, the principle components analysis is used to detect the quantitative characters that are positively and significantly correlated with each component of the plot. In the biplot, the correlation between cultivars and quantitative traits was depicted at each principal component. In addition, PC is for grouping of inherently related cultivars into the same class based on the performance under drought stress. G10 recorded the best no. of leaves, leaf length, root length, above dry weight, below dry weight, chlorophyll a, and chlorophyll b. G19 had the best no. of leaves, leaf width, root diameter, carotenoids, chlorophyll a, chlorophyll b, and chlorophyll content. G9 performed well for leaf length, leaf width, root length, and root diameter. G6 performed well for leaf length, leaf width, root diameter, carotenoids, chlorophyll a, and chlorophyll content. G16 revealed better root length, number of lateral root, above dry weight, carotenoids, chlorophyll a, chlorophyll b, and chlorophyll content. G17 TGm-4499 recorded better no. of leaves, above dry weight, carotenoids, chlorophyll a, and chlorophyll b. G20 TGm-4022 demonstrated better leaf width, plant height, root length, root diameter, carotenoids, and chlorophyll b. G8 showed better leaf length, root length, and chlorophyll a. The soybean cultivars with best physiological, above, and below ground as described above were the best and tolerant cultivars under PEG-simulated drought stress.

5. Conclusion

This study shows remarkable changes in physiology, root system architecture, and dry matter yields of the 20 accessions under PEG-simulated drought stress. G16 and G19 recorded the highest carotenoid, highest chlorophyll a, and chlorophyll b. The highest dry weight was observed with G16 and G10, while the number of leaves was recorded in G19 and G17. G4, G9, G10, and G13 demonstrated the highest dry weight. For chlorophyll content, G6, G8, G16, G19, and G18 were the best. G10, G19, G9, G6, G16, G17, G20, G16, and G18 are the tolerant cultivars to drought stress on the basis of growth, physiological, and root system architecture.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.