Effects of Drought Stress on Some Agronomic and Morpho-Physiological Traits in Durum Wheat Genotypes

Durum wheat performance in the Mediterranean climate is limited when water scarcity occurs before and during anthesis. The present research was performed to determine the effect of drought stress on several physiological and agro-morphological traits in 17 durum wheat genotypes under two conditions (control and drought) over two years. The results of analysis of variance indicated that the various durum wheat genotypes responded differently to drought stress. Drought stress significantly reduced the grain filling period, plant height, peduncle length, number of spikes per plot, number of grains per spike, thousand grains weight, grain yield, biomass, and harvest index in all genotypes compared to the control condition. The heatmap-based correlation analysis indicated that grain yield was positively and significantly associated with phenological characters (days to heading, days to physiological maturity, and grain filling period), as well as number of spikes per plant, biomass, and harvest index under drought conditions. The yield-based drought and susceptible indices revealed that stress tolerance index (STI), geometric mean productivity (GMP), mean productivity (MP), and harmonic mean (HM) were positively and significantly correlated with grain yields in both conditions. Based on the average of the sum of ranks across all indices and a three-dimensional plot, two genotypes (G9 and G12) along with the control variety (G1) were identified as the most tolerant genotypes. Among the investigated genotypes, the new breeding genotype G12 showed a high drought tolerance and yield performance under both conditions. Hence, this genotype can be a candidate for further multi-years and locations test as recommended for cultivation under rainfed conditions in arid and semi-arid regions.


Introduction
Among cereal crops, durum wheat (Triticum durum Desf.) is widely cultivated in the Mediterranean region and other semi-arid areas of the world [1]. It represents an essential food source, as it provides carbohydrates, dietary proteins, fiber, calcium, zinc, fats, and energy [2]. Durum wheat is usually grown under rainfed conditions where fluctuations in precipitation have caused water scarcity to act as a major limiting factor for crop productivity, particularly when water deficit stress occurs during the flowering and grain filling period stages [3]. It is clear that plant survival and development in stressful environments-especially drought conditions-depend on its adaptive mechanisms and defense systems, such as nutrient and water uptake from the deep layers of soil, osmotic adjustment by the active accumulation of free amino acids, ions, sugars, and proline, biochemical activities, stomatal control of water loss, and phenological plasticity [4,5].
The pattern of biomass allocation is a key growth parameter of an adaptive strategy in plants, especially in drought conditions. It is strictly associated with the size and morphological construction of crop organs [6]. Studies of agro-morphological, biochemical, and physiological responses under rainfed conditions can advance our understanding of the ability of crops to respond and adapt to drought-prone environments. Several agro-morphological traits such as plant height, peduncle length, number of fertile tillers per plant, spike length, number of grains per spike, thousand grains weight, and even awn length not only affect stress tolerance to limited soil moisture in wheat, they also reveal how adaptive genotypes cope with drought via morphological changes [7]. Similar to other crops, the important physio-chemical responses and parameters for evaluating drought response in wheat include chlorophyll content, stomatal conductance, leaf water potential, photosynthetic rate, and transpiration rate [8].
Despite intensive durum wheat breeding programs in the last century, in-depth comprehension of the agronomic and morpho-physiological responses of wheat to drought stress is not available. Hence, agro-morphological assessments might be useful to improve the successive operations of wheat processing and transformation. In particular, the selection of cultivars resistant to drought stress and the adoption of correct agronomic treatments [9] might prevent problems and complications during essential operations of wheat milling and flour processing [10]. Moreover, the development, selection, and cultivation of durum wheat cultivars resistant to drought stress which do not need significant input of fertilizers might decrease carbon and water footprints, reducing environmental impacts [11]. Over a long period of selection and breeding of wheat by plant breeders, complex strategies involving agro-morphological traits along with physiological and biochemical mechanisms have been incorporated to cope with drought stress [12]. Hence, knowledge of the association between yield performance and these traits, as well as accurate estimations of stress tolerance in plant genetic material, can be used to improve the efficiency of breeding programs in diverse environments [13].
The objectives of the current study were to evaluate the 16 agronomic and morpho-physiological responses of 17 durum wheat genotypes exposed to drought stress and identify effective traits to facilitate a screening process in future durum wheat breeding programs that will lead to improved yield performance under rainfed conditions.

Plant Materials and Experimental Layout
The experiments were carried out in the Dryland Agricultural Research Sub-Institute (DARSI) of Iran (Sararood Station, Kermanshah, Iran, 34 • 19 N; 47 • 17 E; 1351 m a.s.l.) A set of 16 durum wheat genotypes along with a bread wheat local check (Table 1) were investigated under irrigated (supplemental irrigation; two times irrigation from heading to mid-grain filling stages) and rainfed (hereon referred to as drought stress) conditions during two cropping seasons (2011-2012 and 2012-2013). The soil texture at the station was clay loam, and the average long-term annual precipitation is an estimated 425 mm. For each experiment, a randomized complete block design with three replications was used. The plot size was 7.2 m 2 , which included 6 rows, each 6 m long, spaced 20 cm apart. The soil was supplemented with 50 kg N ha −1 and 50 kg P 2 O 5 ha −1 .

Physiological Assay
The potential quantum efficiency of photosystem II (Fv/Fm), canopy temperature, relative chlorophyll content (SPAD reading), and stomatal conductance (SC; mmol m −2 s −1 ) were taken on the flag leaf. The Fv/Fm measurements and canopy temperature were taken between 11:00 and 14:00, and between 12:00 and 14:00 using a portable Fluorometer (Optic-Science OS-30p, Hudson, NY, USA) and infrared thermometer (Model: E200IR, Germany), respectively. The SPAD index was measured for three flag leaves in each plot using the SPAD chlorophyll meter (Minolta Co. Ltd., Tokyo, Japan). SC was estimated from three random flag leaves in each plot using the leaf porometer (Decagon Devices, Inc., Pullman, WA, USA).

Agro-Morphological Assay
The following agronomic traits were evaluated: days to heading (DH), days to maturity (DM), grain filling period (GFP), plant height (PH), peduncle length (PL; cm), flag leaf length (FL; cm), number of spikes per plot (NSPP), spike length (SL; cm), number of grains per spike (NGPS), thousand grains weight (TGW; g), grain yield (GY; g plot −1 ), biological biomass (Bio; g plot −1 ), and harvest index (HI; %). All traits were measured based on five randomly selected samples of each genotype during different stages from development stage to maturity. Days to heading and maturity were recorded as number of days from planting until 50% of the heads and physiological maturity in a plot. Plant height was recorded from ground to the tip of the main spike at maturity. Length of flag leaf was measured from leaf sheath to the tip of leaf. Peduncle length was measured from the last node from the top of the plant to the rachis. The main spike length was measured excluding awns. Grain yield and biomass were calculated according to one meter square (g m −2 ) for each genotype, and then converted to g plot −1 . Biomass was estimated from the above-ground tissues including the tillers per plant and spikes. Harvest index was calculated by dividing the grain yield by the biological yield.

Statistical Analysis
The combined analysis of variance (ANOVA) was performed using SAS software ver. 9.1 [14]. Pearson's correlation coefficients were calculated for each pair of traits per growth condition using XLSTAT software (Addisonsoft XLSTAT, Paris, French). Several susceptibility and tolerance indices, along with their ranks, were calculated using the iPASTIC toolkit [15]. Furthermore, the ranking method was used to identify the best genotypes according to the different measured indices following the methodology proposed by Ketata et al. [16]. Based on this approach, the average rank sum (ASR) for all indices was used as an indicator to select the best genotypes. Accordingly, the best performance of each trait received the lowest rank, hence the genotypes with the lowest ASR and standard deviation (SD) values were identified as the best genotypes.

Analysis of Variance, Range of Data, and Means Comparison
The results of the combined analysis of variance (ANOVA) on grain yield and other measured traits indicated significant effects for year (except for PL, TGW, and Fv/Fm), stress (except for SL, FL and SC), and genotypes (TGW, SC, CT, Fv/Fm, and Bio). The interaction of year and stress effects was significant for all traits except SL, PE, FL, TGW, Fv/Fm, SC, and Bio. Moreover, the genotype and year interaction was significant for DH, DM, GFP, NSPP, SC, CT, GY, HI, and Bio, whereas the genotype and stress interaction was significant for DH, DM, GFP, PL, NGPS, NSPP, GY, and Bio. On the other hand, the three-way interaction between year, stress and genotype was only significant for DH, DM, and GFP traits ( Table 2). The coefficient of variance (CV) ranged from 0.48 to 43.06%; the highest value was observed for Fv/Fm (43.06%), followed by PE (32.69%) and GY (26.56%).  Table 2). The response of each genotype for the measured traits under control and drought conditions is shown in Table 3. Under the control condition, DH ranged from 214.83 to 219 days (mean~218), and genotypes G3, G6, and G8 showed a short time of spike appearance. Under the drought condition, DH varied between 213.67 and 217.67 days (mean~215), and most of the tested genotypes showed a similar pattern; G5, G10, and G15 had the highest values. There was low variability in days to maturity under both conditions; this trait varied between 170 and 176.67 days (average~173) under the control condition, and between 173.83 and 179.33 days (average~176) under the drought conditions. Under the control condition, genotypes G7, G10, and G12 were identified as the early-maturing genotypes; under drought condition, the early-maturing genotypes were G6, G7, and G8. The GFP indicated a low variability in both conditions. In the control, GFP ranged from 41.50 to 47.83 days (average~45), and the best genotypes were G4, G5, and G8. In the drought condition, this phonological trait varied between 35.17 and 41.50 days (average~39), and G6, G8, and G10 had the shortest grain filling period. The PL ranged from 24.47 to 33.28 cm (average 29.52 cm) under the control condition, whereas in the drought condition it ranged from 16.86 to 26.53 cm (average 23.85 cm). Under control conditions, G7, G8, and G16 showed the highest values for plant height; under drought, the highest plant heights were seen in G6, G12, and G14.
In the control condition, there was no significant difference among the tested genotypes in terms of PE, while under drought conditions this trait varied between 3.97 and 11.83 cm (average 8.58 cm). Therefore, the ideal genotypes with the highest peduncle length were G8, G16, and G17. The range of variability across all 17 tested durum genotypes for FL was limited in both growth conditions; it ranged from 12.95 to 17.97 cm (average 15.13 cm) and from 12.69 to 16.94 cm (average 15.34 cm) under control and drought conditions, respectively. The highest values for FL under the control condition were recorded in genotypes G2, G5, and G15; under drought condition, the highest values were seen in genotypes G2, G5, and G13 compared to other genotypes. Likewise, the range of variability for SL was limited, and it varied between 6.08 and 9.45 cm (average 6.94 cm) and between 6.03 and 8 cm (average 6.93 cm) under control and drought conditions, respectively. The longest spikes were observed in G3 and G8 ( Table 2).
The NSPP indicated a high variability, ranging from 431.67 to 895.83 (average 576.23) under control conditions, and from 297.50 to 627.83 (average 457.52) under the drought condition. G3, G14, and G16 in the control and G2, G3, and G8 in the drought condition had the highest number of spikes per plants. The range of variability for NGPS was considerable--it spanned from 23.32 to 49.07 (average 42.27) and from 19.50 to 35.13 (average 31.19) under control and drought conditions, respectively. Genotypes G12, G16, and G17 in the control condition, and G1, G13, and G14 in the drought condition indicated the highest number of grains per spike. In the control condition, there was no significant difference among genotypes in terms of TGW, whereas under the drought condition genotypes responded to drought in different ways. Indeed, under the drought condition, TGW varied between 19.06 and 27.61 (average 24.05), and the highest values were recorded in G2, G4, and G6. The tested durum wheat genotypes in the present study revealed a high variability for GY. In the control condition, GY ranged from 463.77 to 888.57 g plot −1 (average 721.18 g plot −1 ); genotypes G1, G10, and G16 indicated the highest grain yield. Under the drought condition, GY ranged from 356.17 to 546.63 g plot −1 (average 458 g plot −1 ); G1, G9, and G12 were identified as the highest yield genotypes (Table 2). There was a high rate of variation in Bio across the 17 durum wheat genotypes in the control condition. Three genotypes (G1, G14, and G16) showed the highest amount of dry biomass per plot. Furthermore, HI ranged from 27.75 to 37.36% (average 33.39%), and G1, G10, and G17 were identified as the best genotypes. However, there were no significant differences in Bio and HI among the tested genotypes under the drought condition. Furthermore, under both growth conditions there were no significant differences among the tested genotypes for Fv/Fm and CT traits. Nevertheless, SC displayed a range of variability, which spanned from 25.79 to 58.19 mmol m −2 s −1 (average 40.02 mmol m −2 s −1 ). Genotypes G8, G13, and G14 were distinct from other genotypes as the best lines, with a high stomatal conductance under drought conditions.

Association Among Measured Traits
To visualize the effects of drought stress on the relationships among the measured traits, we used two heatmaps based on Pearson's correlation coefficients among grain yield and physiological and agro-morphological traits in both control and drought conditions. Under the control condition, GY showed a positive and significant correlation with DH, GFP, CT, NGPS, and Bio. The association pattern between Bio and other traits was similar to GY. NGPS showed a positive and significant correlation with DH and GFP. Among the physiological traits, only CT showed a positive and significant correlation with some traits such as GY, Bio, NGPS, DH, and GFP. Other correlations between measured traits are shown in Figure 1. Under the drought condition, positive and significant correlations were found between the following parameters: GY and DH, DM, GFP, NSPP, Bio and HI; PL and PE; CT and SC, FL and SL; SC and NGPS; TGW and NGPS (Figure 1).

Association Among Measured Traits
To visualize the effects of drought stress on the relationships among the measured traits, we used two heatmaps based on Pearson's correlation coefficients among grain yield and physiological and agro-morphological traits in both control and drought conditions. Under the control condition, GY showed a positive and significant correlation with DH, GFP, CT, NGPS, and Bio. The association pattern between Bio and other traits was similar to GY. NGPS showed a positive and significant correlation with DH and GFP. Among the physiological traits, only CT showed a positive and significant correlation with some traits such as GY, Bio, NGPS, DH, and GFP. Other correlations between measured traits are shown in Figure 1. Under the drought condition, positive and significant correlations were found between the following parameters: GY and DH, DM, GFP, NSPP, Bio and HI; PL and PE; CT and SC, FL and SL; SC and NGPS; TGW and NGPS (Figure 1).

Screening Drought-Tolerant Genotypes due to Yield-Based Indices
To identify the most tolerant genotypes, several drought tolerance and susceptibility indices were estimated. The actual values and ranking pattern for each genotype, along with the average of sum of ranks (ASR) across all indices for each genotype, are shown in Tables S1-S3. In the first year (2012-2013), genotypes G8, G9, and G10 had the lowest ASR values and highest grain yield and hence were distinguished from other genotypes (Table S1). In the second year (2013-2014), G1, G2, and G3

Screening Drought-Tolerant Genotypes Due to Yield-Based Indices
To identify the most tolerant genotypes, several drought tolerance and susceptibility indices were estimated. The actual values and ranking pattern for each genotype, along with the average of sum of ranks (ASR) across all indices for each genotype, are shown in Tables S1-S3. In the first year (2012-2013), genotypes G8, G9, and G10 had the lowest ASR values and highest grain yield and hence were distinguished from other genotypes (Table S1). In the second year (2013-2014), G1, G2, and G3 revealed the lowest ASR values and the highest grain yield under drought conditions (Table S2). Regarding the average data across both years, G1, G9, and G12 were identified as the most tolerant genotypes compared to the others (Table S3). The correlation coefficients among grain yields and the seven indices are shown in Figure 2. In all analyses, the MP, GMP, HM, and STI indices indicated a positive and significant association with Yp and Ys. Thus, these indices can be used as alternative indicators for selecting the high-yielding genotypes under drought conditions. Based on Fernandez's theory (1992), all tested genotypes were classified into four groups (A-D). Accordingly, the genotypes belonging to group A-with high yield performance under control and drought conditions-can be selected as the highest yielding and most tolerant genotypes. To illustrate the separation of the tested genotypes into each group, a three-dimensional plot based on stress tolerance index (STI) and grain yield under control and drought conditions (Yp and Ys, respectively) was rendered. As shown in Figure 3, the following genotypes fell in group A: G1, G7, G10, G13, G14, G15, G16, and G17 in the first year; G1 and G12 in the second year; G1, G9, G10, G11, and G12 across both years.  (Table S2).
Regarding the average data across both years, G1, G9, and G12 were identified as the most tolerant genotypes compared to the others (Table S3). The correlation coefficients among grain yields and the seven indices are shown in Figure 2. In all analyses, the MP, GMP, HM, and STI indices indicated a positive and significant association with Yp and Ys. Thus, these indices can be used as alternative indicators for selecting the high-yielding genotypes under drought conditions. Based on Fernandez's theory (1992), all tested genotypes were classified into four groups (A-D). Accordingly, the genotypes belonging to group A-with high yield performance under control and drought conditions-can be selected as the highest yielding and most tolerant genotypes. To illustrate the separation of the tested genotypes into each group, a three-dimensional plot based on stress tolerance index (STI) and grain yield under control and drought conditions (Yp and Ys, respectively) was rendered. As shown in Figure 3, the following genotypes fell in group A: G1, G7, G10, G13, G14, G15, G16, and G17 in the first year; G1 and G12 in the second year; G1, G9, G10, G11, and G12 across both years.

Discussion
Among the various environmental stresses, drought is the major detrimental factor for growth and productivity [17]. The development of high-yielding cultivars with acceptable stability and adaptability across different environments is one promising strategy for improving wheat yield under drought-prone environments [18]. In general, drought stress has a negative effect on agronomic and physiological traits in wheat, regardless of the genotype and the time of stress application [19]. Among different agronomic traits, yield and its related traits are considered to be the most important components for improving wheat productivity. Indeed, grain yield can be increased through manipulation of yield-related traits such as spike length, number of spikelets per spike, and grains per spike [20]. However, several factors including pollen abortion, increased number of sterile tillers and reduction in food reserves can negatively influence the yield performance during drought stress [21]. In the present study, genotypic variation in response to drought stress was investigated in 17 durum wheat genotypes based on several phenological, morpho-physiological and agronomic traits. A significant interaction effect between stress and genotype revealed differences in how the various genotypes reacted under control and drought conditions ( Table 2). Drought stress negatively affected the grain filling period, plant height, peduncle length, number of spikes per plot, number of grains per spike, thousand grains weight, grain yield, biomass and harvest index across all tested genotypes. In particular, drought stress reduced the thousand grains weight, number of grains per spike and peduncle length more than the other traits. The significant reduction in grain yield due to drought stress has also been shown in previous studies [7,19,22].
Grain filling period is an important phenological trait that is associated with current photosynthesis and relocation of assimilates from reserve pools in vegetative tissues. Importantly, drought shortens the life cycle and duration of grain filling. Under such conditions, the grain filling period decreases due to accelerated leaf senescence, reduced photosynthesis, and sink limitations.

Discussion
Among the various environmental stresses, drought is the major detrimental factor for growth and productivity [17]. The development of high-yielding cultivars with acceptable stability and adaptability across different environments is one promising strategy for improving wheat yield under drought-prone environments [18]. In general, drought stress has a negative effect on agronomic and physiological traits in wheat, regardless of the genotype and the time of stress application [19]. Among different agronomic traits, yield and its related traits are considered to be the most important components for improving wheat productivity. Indeed, grain yield can be increased through manipulation of yield-related traits such as spike length, number of spikelets per spike, and grains per spike [20]. However, several factors including pollen abortion, increased number of sterile tillers and reduction in food reserves can negatively influence the yield performance during drought stress [21]. In the present study, genotypic variation in response to drought stress was investigated in 17 durum wheat genotypes based on several phenological, morpho-physiological and agronomic traits. A significant interaction effect between stress and genotype revealed differences in how the various genotypes reacted under control and drought conditions ( Table 2). Drought stress negatively affected the grain filling period, plant height, peduncle length, number of spikes per plot, number of grains per spike, thousand grains weight, grain yield, biomass and harvest index across all tested genotypes. In particular, drought stress reduced the thousand grains weight, number of grains per spike and peduncle length more than the other traits. The significant reduction in grain yield due to drought stress has also been shown in previous studies [7,19,22].
Grain filling period is an important phenological trait that is associated with current photosynthesis and relocation of assimilates from reserve pools in vegetative tissues. Importantly, drought shortens the life cycle and duration of grain filling. Under such conditions, the grain filling period decreases due to accelerated leaf senescence, reduced photosynthesis, and sink limitations. Furthermore, the shortened grain filling period directly influences grain number and grain size, which largely accounts for the decrease in wheat yields [23]. Our results indicated that drought stress decreased the grain filling period by~14% as compared to the control condition (Table 2), suggesting substantial decrease in grain yield in the 17 durum wheat genotypes. Drought stress remarkably increased canopy temperature (~21%) compared with the control conditions (Table 2). During drought periods, plants are often subjected to high temperatures. Increasing leaf temperature limits the activity of photosynthetic enzymes (such as Rubisco) leading to a reduction of carbon dioxide fixation [24]. Limitation of carbon dioxide fixation increases the rate of active oxygen formation in chloroplasts, which in turn would either manifest as oxidative damage to the plant or result in the activation of defense systems that could prevent such damage from occurring [25]. In this situation, most of the plants can adapt through the protection of their photosystem apparatus II (PSII) from photodamaging effects by dissipating the excess excitation energy thermally through down-regulating PSII [26]. The Fv/Fm parameter can be applied to determine the potential efficiency of PSII [27,28]. As shown in Table 3, drought stress increased this parameter (3.5%) compared to the control condition (Table 1). This result indicates that some of the tested genotypes (especially G9) may have a strong capacity to protect themselves through maintenance of high levels of PSII activity. Although several studies reported that drought stress reduces Fv/Fm [1,17,29], our results revealed that the process of energy displacement from the antenna chlorophyll complex to the PSII may not be affected by drought stress in some of the tested genotypes. The stomatal conductance (Gs) as a physiological trait has two important roles: (1) maintenance of water in leaf tissue by control of transpiration, and (2) determining the rate of gas exchange by movement of gases between the inner and outer environments of the leaf tissue. Drought stress reduces Gs and the reaction centers in PSII, as well as exposure of excess energy in the chloroplasts. All of these events restrict photosynthesis [30]. Our results revealed that although Gs was reduced by drought stress (~3% compared to the control conditions) across all tested genotypes, some of the genotypes nevertheless showed high stomatal conductivity compared to the control. Indeed, this result confirms that some of the genotypes have a protected photosynthesis apparatus and can maintain the health of photosynthesis mechanisms under drought conditions.
The efficiency of breeding programs in diverse environments can be improved by gaining an understanding of the associations between grain yield and different selection indices. Although many agronomic, morphological, and physiological traits have been studied for their use in breeding programs for drought tolerance, only a few traits are currently recommended for application in practical programs. For instance, several traits such as stay-green character [31], canopy temperature and spike fertility [32], number of grains per plot [33], stomatal conductance and relative chlorophyll content [17,29], water use efficiency [34], root architecture [35], chlorophyll fluorescence [17,29,36] were used as selection indicators due to their correlations with grain yield under drought conditions. In the present study, a positive correlation between grain yield and other phenological characteristics (number of days to spike heading, physiological maturity, and grain filling period), as well as number of spikes per plant, biomass, and harvest index under drought conditions were found. Furthermore, negative significant correlations were observed between grain yield and canopy temperature, flag leaf, and plant height. Altogether this suggests that it is relatively easy to develop high-yielding durum genotypes with increasing numbers of spike and phonological characters-especially grain filling period-as well as decreased plant stature, canopy temperature, and flag leaf length. There was no significant correlation between yield and Fv/Fm, thousand grains weight, and peduncle length, which indicates that grain yield can potentially be improved without changing these traits under drought stress conditions (Figure 1). Following Fernandez's theory [37], the genotypes formed four distinct groups according to their yield response to stress conditions. The first group (A) contains the genotypes with relatively equal yield in both control and stress conditions. The second group (B) includes the genotypes that show a high yield performance only in the control condition. The genotypes that have a high yield performance in stress conditions are classified as the third group (C), and finally the fourth group (D) includes the genotypes with low yield performance in both conditions. Based on these results, several yield-based stress tolerance and susceptibility indices have been applied to detect the tolerant genotypes. In the present study, seven drought tolerance and susceptibility indices were calculated based on grain yields under control and drought conditions. According to all analyses for each year separately as well as the average of both years, significant and positive associations between grain yields with STI, HM, MP, and GMP were found ( Figure S1). This result reveals the capacity of these indices to identify genotypes with high potential yield and tolerance to drought conditions (Group A). Likewise, the ability to separate Group A genotypes from the others using these indices is consistent with the findings reported for canola [38], chickpea [39], maize [40], and durum wheat [22]. Among the yield-based indices, STI is a suitable index for characterization of the most tolerant to drought stress. Indeed, its applicability has been confirmed in numerous studies [1,17,22,29,35,38,41,42]. Based on the three-dimensional plots, several different genotypes were distributed in Group A each year ( Figure S2). However, the genotype G1 showed the same response in all plots. To better understand the responses of the tested durum wheat genotypes and identify the most drought-tolerant genotypes, the average of the sum of ranks (ASR) across all indices was used. Based on this data, three genotypes (G1, G9, and G12) with high yield performance (546.63, 518.60, and 513.08 g plot −1 , respectively) indicated low ASR values; hence, these genotypes can be recommended as the drought-tolerant genotypes. In contrast, G5, G13, and G14 with low yield performance (356.17, 387.62, and 384.60, g plot −1 , respectively) revealed high ASR values, which indicated their susceptibility to drought stress. Among all genotypes, G12 as a new breeding line reveals good performance under both conditions (817.67 g plot −1 under control, and 513.08 g plot −1 under drought conditions), similar to the check variety (G1). Therefore, it can be recommended for cultivation under rainfed conditions in the arid and semi-arid regions of Iran and similar climates.

Conclusions
To summarize, our results revealed that durum wheat genotypes were significantly affected by drought stress for grain yield and several agronomic and morpho-physiological traits. This resulted in considerable variation in the measured traits and drought tolerance in these genotypes that should be considered for durum wheat breeding programs. As seen in this study, positive associations of number of grains per spike and phenological characters with grain yield make it possible to evaluate genotypes with high yield performance by using these character attributes in rainfed conditions. Among the tested genotypes, G12 as a new breeding line responded well to drought stress. It can thus be used for cultivation in drought-prone environments and to develop drought-tolerant varieties.
Supplementary Materials: The following are available online at http://www.mdpi.com/2071-1050/12/14/5610/s1, Table S1. Yield performance of 17 durum wheat genotypes under control (Yp) and drought stress (Ys) conditions, along with the yield-based drought tolerance and susceptibility indices in the first year (2011-2012); Table S2. Yield performance of 17 durum wheat genotypes under control (Yp) and drought stress (Ys) conditions, along with the yield-based drought tolerance and susceptibility indices in the second year (2012-2013); Table S3. Yield performance of 17 durum wheat genotypes under control (Yp) and drought stress (Ys) conditions, along with the yield-based drought tolerance and susceptibility indices over both years (2011)(2012)(2013).