- 1Department of Genetics and Plant Breeding, Banda University of Agriculture and Technology, Banda, India
- 2Crop Improvement Division, Indian Council of Agricultural Research (ICAR)-Indian Institute of Wheat and Barley Research, Karnal, India
- 3Krishi Vigyan Kendra, Jhansi, Banda University of Agriculture and Technology, Banda, India
- 4Uttar Pradesh (UP) Council of Agricultural Research, Lucknow, India
- 5Department of Plant Science and Landscape Architecture, University of Connecticut, Storrs, CT, United States
Climate change has resulted in extreme temperature and drought around the globe, which poses serious threat to food security. Both heat and drought stress affects the production and productivity of wheat crop. The present study was undertaken to evaluate 34 landraces and elite cultivars of Triticum spp. for phenological and yield-related traits under optimum, heat, and combined heat–drought stress environments during 2020–2021 and 2021–2022. The pooled analysis of variance showed significant genotype × environment interaction, suggesting an influence of stress on trait expression. The trait performance of genotypes exhibited significant reduction under combined heat–drought stress as compared to optimum and heat stress environments. The maximum seed yield penalty was observed under combined heat–drought stress environment as compared to heat stress alone. Regression analysis indicated significant contribution of number of grains per spike towards stress tolerance. Based on Stress Tolerance Index (STI), genotypes Local-17, PDW 274, HI-8802, and HI-8713 were identified to be tolerant to both heat and combined heat and drought stress at Banda, whereas genotypes DBW 187, HI-8777, Raj 4120, and PDW 274 were tolerant at Jhansi location. The genotype PDW 274 showed stress tolerance under all treatments at both the locations. The genotypes PDW 233 and PDW 291 showed highest stress susceptibility index (SSI) across the environments. The number of grains per spike and test kernel weight were positively associated with seed yield across the environments and locations. The selected genotypes Local-17, HI 8802, and PDW 274 were identified to be the potential sources of heat and combined heat–drought tolerance, which may be utilized in hybridization to develop tolerant wheat genotypes and also for mapping of underlying genes/quantitative trait loci (QTLs).
Introduction
Wheat is the second most important cereal crop after rice, cultivated on approximately 30 million hectares with a production of 108 million tons in India (Annonymous, 2020). Globally, it is a vital staple food crop, feeding 2.5 billion people in more than 85 countries and a source of 20% calories in the human diet (Yashavanthakumar et al., 2021). Wheat productivity was reported to be affected up to 21% due to drought on a global scale from 1980 to 2015 (Daryanto et al., 2016). Global food security is threatened by the ever-increasing world population and by climate change (Lesk et al., 2016). The reduced rainfall, changing rainfall pattern, and distribution coupled with shorter winter seasons result in extreme temperature due to climate change (Lobell et al., 2011). Among abiotic stresses, drought and heat stress are the most important devastating factors for the growth and development of crop plants, affecting yield potential and end-use quality of food products (Delahunty et al., 2015). Furthermore, unpredictable yield losses by temperature and drought stress depends on annual precipitation patterns, surface water evapo-transpiration, and soil water holding capacity. The response of plants to water stress depends on several factors, such as developmental stage, intensity and duration of stress, and cultivar genetics (Wahid et al., 2007; Eskandari and Kazemi, 2010). Under water-deficit conditions, various physiological activities, particularly water and nutrient uptake, photosynthesis, partitioning of photosynthates, and metabolic activities, are affected, which ultimately results in significant yield loss (Farooq et al., 2009; Praba et al., 2009). Heat stress also impacts various metabolic activities like protein synthesis, inactivation of enzymatic activities, and physiological activities of the cell including cell membrane damage. The cell division is also adversely affected by heat stress (Smertenko et al., 1997). The damages from heat shocks seriously limit plant growth and favor oxidative damage. The terminal heat stress accelerates the assimilation of photosynthates due to early senescence, which results in shrunken grains and poor quality (Porter, 2005). Although the increasing temperature is also beneficial for crop production in some cooler regions of the world, the overall impact on global food security is still negative (Challinor et al., 2014).
Climate change has a negative impact on wheat production through extreme drought and temperature stress conditions during crop cycle. Hence, the identification of heat- and drought-tolerant genotypes is of paramount importance to the breeders. Phenotyping remains a key criterion for screening breeding materials based on drought-adaptive and constitutive morpho-physiological characteristics including yield and its components (Monneveux et al., 2012; Passioura, 2012). Seed yield and component traits-based physiological indices are widely used in wheat breeding programs for the selection of tolerant genotypes. Attempts to measure the degree of tolerance with a single parameter provide results of a limited value because of the multiplicity of the factors and their interactions contributing to drought tolerance under field conditions. Thus, there is a need to select genotypes with a good combination of agronomical and physiological important traits, cumulatively contributing to improved yields under target drought and heat environments (Tardieu, 2012). In combination, heat and drought stress harm plant growth and productivity more than any other environmental stress, although heat and drought stress individually have been extensively studied in wheat crop. A few studies were conducted to understand the plant response to high temperature during grain filling stage and in combination with drought stress conditions. The experimental sites in the Bundelkhand region are characterized by high temperature and water deficit during rabi season (rabi season is when the crops are sown around mid-November, preferably after the monsoon rains are over, and harvesting begins in April/May), and the studies on evaluation of wheat germplasm for heat and drought tolerance in Bundelkhand region are very limited. Generally, the crop is exposed to heat stress by planting late, and comparison is made with the normal planting, whereas the drought stress tolerance is estimated by comparing the performance of wheat genotypes under normal and water-deficit conditions (Liu et al., 2015; Gupta et al., 2020; Shokat et al., 2020). In the present study, the heat and combination of heat–drought stress tolerance was assessed to capture the agro-morphological variability for tolerance that could be utilized for developing tolerant wheat genotypes suiting the stress environments.
Materials and methods
Plant material and experimental site
The present investigation was carried out with 34 diverse wheat genotypes (Table 1), of which 18 were landraces collected across the Bundelkhand region (Uttar Pradesh); seven commercial durum cultivars were obtained from Indian Institute of Agricultural Research (IARI), Regional Research Station, Indore (MP), which are recommended for cultivation in Central Zone (CZ); seven durum cultivars were obtained from Punjab Agricultural University (PAU), Ludhiana (Punjab) representing Northwest Plane Zone (NWPZ), and two bread wheat cultivars were from IIWBR Karnal (Haryana). The landraces were collected from six districts (Banda, Mahoba, Chitrakoot, Jhansi, Hamirpur, and Lalitpur) of the Bundelkhand region, conserved by farmers for traditional cultivation. The local wheat collections possessed good adaptability to rainfed farming and are preferred by farmers of this region. The experiments were conducted for two consecutive crop seasons 2020–2021 and 2021–2022 at two locations (Table 2):
Table 2 Details of the experiments conducted at BUAT, Banda and KVK, Jhansi during 2020–2021 and 2021–2022 crop seasons.
a. Crop Research Farm, Banda University of Agriculture and Technology Banda, UP (25.5269°N latitude and 80.3418°E longitude)
b. Experimental Farm, Krishi Vigyan Kendra, Jhansi (25.5347° N latitude, 78.5742° E longitude)
Bundelkhand region in Uttar Pradesh state of India is known for extreme heat and drought stress, which severely affect the growth, development, and yield of wheat crop. This region frequently experiences extreme temperature, usually 50°C during summer and 5°C in the winter season. The summer season (April–June) characterized by strong heat wave and dusty, gusty, and dry wind blowing over the entire region, which often lead to fatal heat strokes. The average rainfall of this region is 800–900 mm per annum; however, most of rainwater is lost as runoff due to undulated topography. During winter season (November–March), few rain showers are received occasionally, which supplies inadequate moisture to crops. The whole Bundelkhand region is known for predictable drought and heat stress during the anthesis to grain filling stage of winter season crops.
Experimental design
The experiments were conducted under three environments, viz., the optimum environment (OE), heat stress environment (HSE), and the heat–drought environment (HDSE) (Table 2). The timely sown (OE) and late sown experiments (HSE) were given irrigation (flood irrigation) at the required field capacity throughout the crop period, whereas in the late sown restricted irrigation, no irrigation was given during the reproductive phase (HDSE). The genotypes were planted in a two-row plot of 2 m length with two replicates following Alpha lattice design for each condition. The distance of 5 m was kept between the trials to check the moisture percolations in the drought experiment. The sowing was done with the modified precision numeric seed planter. The sowing and harvesting dates for each environment for both the years is given in Table 2. The crop was supplemented with 120 kg/ha nitrogen and 60 kg/ha phosphorus. Total N was split into three doses of 40 kg/ha each and applied at the three-leaf stage, tillering, and again at heading stage. The recommended agronomic practices of weeding and crop protection were applied to raise healthy crop plants.
Measurement of phenological and yield-related traits
The data on yield and yield-contributing traits were recorded for each genotype in all three environments during 2020–2021 and 2021–2022.
• Days 50% heading (DH): number of days taken from date of sowing to appearance of anthers on 50% plants in a plot
• Days to maturity (DM): number of days taken from date of sowing until physiological maturity completed
• Numbers of tillers/plant (NTP): actual count of the productive numbers of tillers of five plants (spike bearing) per plant
• Spike length (SL): measured from base of spike to the end tip of spike excluding awns in centimeters.
• Number of spikelets/spike: counting number of spikelets per spike on main spike
• Plant height (PH): taken from base of soil surface to terminal of spike excluding awns in centimeters of 10 plants
• Number of grain/spike (NGS): counting grain numbers in five randomly selected spike and divide by 5
• Test kernel weight (TKW): 100 seeds were randomly taken for each genotype, and weight was recorded in grams
• Biological yield (BY): measured by whole plant weight including shoot and leaves with kernel after drying
• Seed yield per plot (SY): measured in grams using a sensitive balance after the moisture of the seed is adjusted to 12%. Total dry weight of grains harvested from each plot was taken as grain yield per plot
• Harvest index (%): estimated by dividing grain yield per plot to biological yield per plot
Calculation of stress tolerance indices
The experiments were conducted under three different environments, viz., the optimum environment (OE), heat stress environment (HSE), and the heat–drought environment (HDSE). The selection indices for heat stress and combined heat–drought stress were calculated based on the relationship among yield under normal and stress environment. The stress tolerance indices were calculated as follows:
Stress tolerance index (STI): STI = (Fernandez, 1992)
Stress susceptibility index: SSI = (Fischer and Maurer, 1978)
Mean productivity (MP): MP = (Rosielle and Hamblin, 1981)
Tolerance (TOL): TOL = Yp−Ys (Rosielle and Hamblin, 1981)
Modified stress tolerance index (MSTI): MSTI = (Farshadfar and Sutka, 2002)
where Yp is the yield under optimal conditions, Ys is the yield under stress, is the average yield of all entries under optimal conditions, and is the average yield of all entries under stress.
MP favors higher yield potential and lower stress tolerance, and the selections based on MP generally increases the average performance of genotypes in both stress and non-stress environments and fails to distinguish between stress tolerant and high yielding genotypes. Larger values of TOL represents more sensitivity to stress; thus, smaller values are favored for selecting tolerant genotypes. The smaller the values of SSI and TOL, the greater the value of tolerance but fails to distinguish between genotypes performing better under both environments and performing better only under stress. Selection based on TOL and SSI favors those genotypes having lower yield potential under normal and high yield potential under stress environments. Thus, TOL also fails to distinguish between tolerant and high yielding genotypes. The STI is used for the identification of wheat genotypes that produce higher yields under non-stressed and stressed environments. The genotypes selected based on STI will have higher stress tolerance and grain yield. Genotypes having higher values of STI are tolerant and have high yield potential, as it involves the stress tolerance intensity and geometric mean productivity (GMP). Thus, STI is the best one to distinguish between genotypes performing superior/poor either under stress or non-stress environments or under both environments.
Statistical analysis
The analysis of variance (ANOVA) was performed, first taking locations, years, and environments as random factor. In the second stage, ANOVA was performed location wise separately for all three environments. The analysis of alpha lattice experimental design for all agro-morphological and physiological traits was computed using Agricolae function of R package. The location, year, environment, genotypes, and their interaction, standard errors and least significant difference test (LSD), and Tukey’s honest significant difference (HSD) test were performed. The basic statistical parameters and graphical representation through box plot of data matrix for all studied traits were performed with IBM SPSS software and PAST3 (version 1.0.0.0) to know about the variability in the traits among germplasm. The correlation analysis was conducted for each environment (pooled across years and locations) using correlation plot function of “metan” package (Olivoto and Lúcio, 2020) in R package.
Results
Environment assessment
The Banda and Jhansi locations are typically represented by hot and dry environment of Central Zone (CZ) in India. The weather data were recorded at the Meteorological Observatory Unit (MOU) of the University at Banda and KVK, Jhansi for two years (Figure 1). The maximum temperature was recorded more than 35°C at the grain filling stage during 2020–2021 and 2021–2022 at the Banda and Jhansi locations. At Banda location, the cropping season during 2020–2021 was extended due to 10 mm rainfall in the month of March, while during 2021–2022, the temperature in the last week of February increased rapidly exposing the crop to terminal heat. During 2020–2021, maximum precipitation (958 mm) was received from June to October during 2020–2021, while 929 mm of rainfall was recorded in the monsoon season during 2021–2022. During the month of February 2021, very less rainfall of 10.57 mm was received in the month of February 2021, whereas rainfall of 104 mm was received at early vegetative stage of the crop in the first week of January 2022. During 2022, the rainfall was received only in the early vegetative stage, and the crop was exposed to drought stress during reproductive phase. The OE trials were irrigated at both locations to avoid the effect of heat and drought stress. The relative humidity (RH) at the grain filling stage during 2020–2021 was high due to rain in the month of March, while RH was low during 2021–2022 at the time of grain filling.
Figure 1 Weather parameter rainfall (mm), temperature °C (maximum and minimum) relative humidity (%) in two growing season 2020–2021 and 2021–2022 at Banda and Jhansi locations.
Analysis of variance
Combined analysis of variance using linear mixed model (LML) indicated significant effects of environment, location, year, and their interaction on most of the measured traits (Table 3A). The genotypes exhibited highly significant variation for all traits except TKW and HI. The pooled analysis of variance indicated significant differences in all interactions; thus, a separate analysis of variance was performed for each location. The individual location ANOVA analysis also indicated significant differences between genotypes, and Gen × Env interaction was also significant for most of the traits (Tables 3B, C).
Table 3A Combined analysis of variance of agronomic traits across environment and locations in two growing periods (2020–2021 and 2021–2022).
Effects of terminal heat and combined heat–drought stress on yield and yield components
Phenotypic trait variation was observed across 34 genotypes under OE, HSE, and HDSE, indicating considerable genetic variation in performance across environments. Some genotypes exhibited higher reduction, while other genotypes managed to grow and remain relatively healthy. The performance of genotypes under combined heat–drought stress was severely affected compared to heat stress alone at Banda and Jhansi. The mean of 11 agro-morphological traits of the wheat lines under optimal and stress conditions is presented in Table 4 and Figure 2. The highest heritability was reported for PH (98%) followed by TKW (95%), whereas least was reported for NSS (45%).
Table 4 Basic statistics of agronomic traits under optimum, heat stress, and heat–drought stress environments at Banda and Jhansi for two consecutive growing seasons.
Figure 2 Performance of wheat genotypes for different agronomic traits DFF, day to 50% flowering (number of days); DM, days to maturity (number of days); PH, plant height (cm); NT, number of tillers per plant (number); SL, spike length (cm); NSS, number of spikelets per spike; NGS, number of grains per spike; TKW, test kernel weight (g); BY, biological yield (g); HI, harvest index; SY g/plot, seed yield (g/plot) under all three environments at Banda and Jhansi locations.
Impact of heat stress on agro-morphological traits
Heat stress significantly affected the performance of all the traits. At Banda location, highest reduction was reported in NT (24%) followed by SY (20%), DM (18%), and BY (16%) as compared to optimum environment, whereas least reduction of 4% was reported in PH and NSS (Table 5). Similarly at Jhansi, highest reduction was reported in SY (44%) followed by (34%), and DM (22%) as compared to optimum environment, whereas least reduction of<5% was reported in PH and NSS (Table 5). The maximum yield reduction of 48% was recorded in genotype PDW-233 at Banda, whereas yield reduction of 66% was reported in Local-5 at Jhansi location. The genotype Local-4 recorded minimum SY reduction of 5.6% at Banda, while a minimum 15% SY reduction was recorded in the HI 8627 genotype at Jhansi.
Table 5 Percent reduction in yield components of bread wheat under heat stress and combined heat–drought stress environment at Banda and Jhansi location.
Impact of combined heat–drought stress on agro-morphological traits
Combined heat–drought stress also affected agro-morphological traits across years at both locations. At Banda, highest reduction was recorded for BY (37%) followed by SY (36%) and NT (27%), whereas the least affected trait was NSS, which showed a reduction of 8% under combined heat–drought stress. Spike yield (59%) showed the highest reduction at Jhansi followed by (57%) and DM (25%), whereas NSS recorded the minimum reduction of 7% (Table 5). The PH and NSS traits were the least affected traits under heat–drought stress at both locations. Highest SY was reported for genotype Local-17 (985 g) at Banda, whereas the genotype PDW-274 recorded with a maximum SY of 490 g at Jhansi. Genotype Local-3 and HI8627 recorded minimum SY reduction of 19% and 35% at Banda and Jhansi, respectively. Under HDSE, DFF significantly reduced from 88 days to 74 days. Genotypes Local-17 and Local-13 were earliest in heading at Banda, whereas genotypes Local-3, Local-9, HI-8777, Local-13, Local-14, and DBW 187 were earliest in heading at the Jhansi.
Among all genotypes, Local-17 and HI 8802 performed well under HSE and HDSE with higher yield at the Banda location. At Jhansi, genotypes DBW 187 and HI 8777 performed well under HSE with higher yield, while under HDSE genotype, PDW 274 and HI 8777 were better performing genotypes.
Correlations analysis
Correlation analysis was undertaken to assess the association of traits under normal and stress environments using pooled environment wise data (Tables 6A, B). The SY showed significant positive correlation with NGS (0.73**) and HI (0.78**), whereas it showed significantly negative correlation with NT and PH under optimum environment. Under heat stress environment, SY showed significantly positive correlation with HI (0.89**), NGS (0.80**), NSS (0.49**), and TKW (0.66**). In contrast to this, in the combined heat–drought environment, SY indicated significant correlation with BY (0.37**), NGS (0.68**), NSS (0.37**), HI (0.42*), and TKW (0.38**). There was a significantly negative correlation of PH (−0.44*) and NT (−0.61) under HSE, whereas SY had a significantly negative correlation with DFF (−0.53**) and DM (−0.61**). NGS, NSS, and TKW are the main traits contributing towards yield under stress environments, whereas BY is a key trait to be considered under combined heat and drought stress.
Table 6 Correlation analysis between different agro-morphological traits based on pooled data across different environments.
The regression analysis was also undertaken to assess the contribution of different traits towards final grain yield. Under optimum environment, NGS (R2 = 0.53) followed by TKW (R2 = 0.20), PH (R2 = 0.06), and BY (R2 = 0.03) contributed significantly towards SY (Figure 3A). Under HSE, NGS (R2 = 0.63) contributed more towards SY followed by TKW (R2 = 0.09) and NSS (R2 = 0.04) (Figure 3B), whereas under HDSE, NGS (R2 = 0.46) contributed highest followed by DFF (R2 = 0.19) and DM (R2 = 0.07) (Figure 3C).
Figure 3 Regression analysis of seed yield (SY) with number of grains/spike (NGS) under (A) optimum, (B) heat stress, and (C) combined heat–drought stress environments.
Identification of heat-tolerant genotypes
Stress indices (STI, SSI, MP, TOL, and MSTI) were estimated based on seed yield under stress environments (HSE and HDSE) and OE conditions. Based on the STI and SSI, the genotypes were grouped into tolerant and sensitive to heat and combined drought-heat stress conditions (Supplementary Tables S1–S4).
Heat-tolerant genotypes
Based on STI, five genotypes local-17 (2.16) and PDW 274 (1.88), HI-8802 (1.83), HI-8713 (1.54), and Local-5 (1.32) performed better under heat stress environment at Banda (Table 7). Similarly, five genotypes DBW 187 (1.77), HI-8777 (1.69), Local-4 (1.53), Raj 4120 (1.48), and PDW 274 (1.40) were tolerant at Jhansi (Table 7). Thus, these genotypes are considered as heat stress tolerant. Based on SSI, genotypes PDW 233 (2.05) and Local-16 (1.96) were identified as heat susceptible at Banda (Supplementary Table S1), whereas at Jhansi, the maximum value of SSI was exhibited in Local-5 (1.34) and PDW 233 (1.33) (Supplementary Table S3); thus, these genotypes are considered as heat stress susceptible. The maximum value of TOL was recorded in PDW 233 (437) and Local-16 (365) in the HSE at Banda (Supplementary Table S1), whereas genotypes Local-5 (522) and PDW 233 (516) had a maximum value of TOL under HSE at Jhansi (Supplementary Table S3).
Table 7 Seed yield and tolerance indices of five highly heat tolerant genotypes of wheat at Banda and Jhansi. based on mean of two years data.
Heat–drought-tolerant genotypes
Under HSDE, highest values for STI were recorded in genotypes Local-17, PDW 274, HI-8802, HI-8713, and HI-8759 at Banda (Table 8). At Jhansi, genotypes PDW 274, DBW 187, HI-8777, Raj 4120, and Local-4 had the highest value of STI under HDSE. Hence, these genotypes can be considered as heat–drought stress tolerant. Based on SSI, the higher value was recorded in genotypes PDW 233 and PDW 291 at Banda (Supplementary Table S2), whereas at Jhansi, genotypes PDW 233 and DWL 5023 recorded the highest values (Supplementary Table S4). The maximum values of TOL recorded were in genotypes PDW 233 and PDW 291 at Banda, whereas genotypes PDW 233 and Local-5 recorded maximum values at Jhansi.
Table 8 Seed yield and stress tolerance indices of five highly heat–drought-tolerant genotypes of wheat at Banda and Jhansi based on mean of two years data.
In summary, genotypes Local-17 and PDW 274 were selected to be heat and combined heat and drought tolerant along with yield at Banda. Similarly, genotypes DBW 187, HI 8777, and PDW 274 performed better at Jhansi under both environments compared to the rest of the genotypes.
Discussion
The confounding effect of heat along with drought stress are more detrimental to plant performance than any individual stress (Barnabas et al., 2008; Li et al., 2013). In the present study, 34 wheat genotypes were evaluated at unexplored field screening sites Banda and Jhansi for heat and combined heat–drought stress during 2020–2021 and 2021–2022. The study aims to quantify the effects of stress on phenological traits development and yield loss. The delayed sowing and restricted irrigation are efficient strategies to identify heat- and drought-tolerant genotypes in a number of crops like bread wheat (Qaseem et al., 2019; Hussain et al., 2021; Yashavanthakumar et al., 2021), Durum wheat (Aberkane et al., 2021; Ayed et al., 2021; Laddomada et al., 2021), mustard (Limbalkar et al., 2021), lentil (Sellami et al., 2021; El haddad et al., 2020; Choudhury et al., 2012), and chickpea (Krishnamurthy et al., 2011; Devasirvatham et al., 2015). The results demonstrated that high temperature coupled with water stress is more destructive to the plant development process and, finally, reflected in yield reduction. The heat and drought stress in combination is more pronounced at Jhansi as compared to that at Banda location.
Effect of stress on plant phenology and grain yield
The weather conditions of the Bundelkhand region are hot and dry during the reproductive phase for the cool season crops like wheat. The temperature and rainfall pattern are consistently favoring experiments focused on screening germplasm for abiotic stresses, more specifically heat and drought stress. The meteorological data clearly demonstrate that more than 95% rainfall was received from June to September in both cropping seasons at Banda. Similarly, at the Jhansi location during monsoon season, the major part of the total rainfall was received. Therefore, at both locations, weather conditions were ideal for challenging plants to stress conditions (Figure 1). During crop season, very less amount of precipitation was received in the reproductive phase, while crop water requirement in the vegetative phase was maintained by irrigation. The temperature in HSE and HDSE experiments was >35°C during the anthesis period, and it reached up to 40°C during the grain filling stage, although experiments under optimum environments were raised under non-stressed conditions.
Our result showed that combined stress penalty on seed yield at Jhansi was more severe (59% under HDSE as compared to 44% under HSE). The maximum yield penalty (36%) was recorded in HDSE, whereas yield penalty of 20% was estimated under heat stress at Banda location (Table 4). The maximum yield for each genotype was harvested at Banda due to the long duration of crop compared to Jhansi. The climatic conditions particularly rainfall and temperature patterns and recorded yield penalty due to stress environments at both locations reflect that Bundelkhand is an appropriate location for terminal heat and drought stress evaluation.
The combined stress is more detrimental on plant phenological traits than heat stress alone at both the locations (Table 5). Several reports have shown that the combination of heat and drought stress is common in dry and semi-dry regions worldwide and reported extensive yield losses. The results of the present study are in agreement with the study of Yashavanthakumar et al. (2021), which reported the confounding effect of heat and drought resulting in yield loss (56%), followed by drought (41.1%) and least affected by heat stress alone (4.8%) over the control environment. Results of the present study also showed that heat and heat–drought stress adversely affected all measured trait significantly, which is also reported in other wheat crop experiments (Hussain et al., 2021; Yashavanthakumar et al., 2021) and Durum wheat (Aberkane et al., 2021; Ayed et al., 2021; Laddomada et al., 2021). Drought stress at pre- and post-anthesis stage significantly reduced grain yield (Y), spike length (SL), number of grains spikes-1 (NGS) and thousand kernel weight (TKW), while kernel abortion (KA) was increased (Shokat et al., 2020). Our study is in agreement with Qaseem et al. (2019), who reported that the combined heat and drought stress is a more devastating factor to plant development as compared to any stress alone. However, there was less impact of stress conditions on PH and NSS compared to optimum environments. Although, both traits are highly stable and were less affected by the environmental factors. This finding is contrary to previous studies, which have suggested that stress has a significant impact on plant height reduction (Sattar et al., 2020). In conclusion, the genotypes Local-17, PDW 274, HI 8802, and HI 8713 have the highest yield potential under OE at the Banda location, whereas DBW 187, HI 8777, Local-4, and PDW 274 were superior at Jhansi location.
Correlation among traits
The association among all measured agro-morphological traits with seed yield was studied using pooled data over environments under OE, HSE, and HSDE (Tables 6A, B). The NGS and TKW directly positively associated with seed yield across the normal and stress environments, indicating that the genotypes were able to maintain both the traits under stress conditions. The selection of the genotypes under stress conditions is most reliable taking NGS and TKW as targeted traits (Dodig et al., 2012; Sareen et al., 2014). There was a significantly negative correlation of PH (−0.44*) and NT (−0.61) with SY under HSE, whereas SY had a significantly negative correlation with DFF (−0.53**) and DM (−0.61**). The traits DFF and DM indicated that heat and drought escape mechanisms as short duration genotype are preferred for drought and high-temperature tolerance (Reynolds et al., 2007; Punia et al., 2011; Mondal et al., 2013; Gupta et al., 2020). The days to maturity and early heading escapes the stress during the critical reproductive phase and increases productivity before the onset of heat and drought stress (Blum, 2010). NGS, NSS, and TKW are the main traits contributing towards yield under stress environments, whereas BY is a key trait to be considered under combined heat and drought stress.
The regression analysis indicated that NGS is the key trait that is contributing towards grain yield under stress environments (Figure 3). High temperature reduces grain weight (Wardlaw et al., 1980) due to reduction in both duration and rate of grain filling and high respiration rate (Tashiro and Wardlaw, 1990); however, yield reduction is primarily due to reduction in grain weight (Guttieri et al., 2001). Tyagi et al. (2003) also reported that grain weight per spike could be used as a measure of heat tolerance. Similarly in the present study, TKW is found to be significantly contributing towards SY under heat stress.
Selecting stress-tolerant genotypes
In the present study, analysis of variance revealed the significant difference among 34 wheat genotypes for all measured traits. The significant variation in genotype, environment, year, and location indicated variable performance of genotypes across environments. The genotype × environment interaction for SY and its associated trait shows the possibilities of selecting genotypes with less yield penalty under stress environments. The significant effect of different environments on the performance of genotypes was reported by Mohammadi et al. (2010) and Talebi et al. (2010) in bread wheat.
Stress selection indices STI, SSI, MP, TOL and MSTI, are widely used for the identification of tolerant genotypes under stress environments. The selection of heat- and drought-stress tolerant genotypes was done using stress tolerance index in wheat was adopted in several studies (Puri et al., 2015; Aberkane et al., 2021a; Poudel et al., 2021). The STI is a more appropriate index to identify tolerant genotypes in the panel of germplasm. Earlier many researchers used STI to screen germplasm in chickpeas (Erdemci, 2018; Shabani et al., 2018). The genotype having higher values of STI, MP, and MSTI indicates their superior performance under both normal and stress environments. Additionally, the genotypes selected based on STI will have higher stress tolerance and grain yield, as it involves the stress tolerance intensity and GMP.
MP favors higher yield potential and lower stress tolerance, and the selections based on MP generally increase the average performance of genotypes in both stress and non-stress environments and fails to distinguish between stress-tolerant and high-yielding genotypes.
Larger values of TOL represent more sensitivity to stress; thus, smaller values are favored for selecting tolerant genotypes. The smaller the values of SSI and TOL, the greater the value of tolerance but fails to distinguish between genotypes performing better under both environments and performing better only in under stress. Selection based on TOL and SSI favors those genotypes having lower yield potential under normal and high yield potential under stress environments. Thus, TOL also fails to distinguish between tolerant and high-yielding genotypes.
The STI is used for the identification of wheat genotypes that produce higher yields under non-stressed and stressed environments. Genotypes having higher values of STI are tolerant and have high yield potential, as it involves the stress tolerance intensity and GMP. Thus, STI is the best one to distinguish between genotypes performing superior/poor either under stress or non-stress environments or under both environments.
In our screening experiments under HSE and HDSE at Banda, genotypes Local-17, PDW 274, HI 8802, and HI 8713 have the highest values of STI, MP, and MSTI; therefore, these are considered as tolerant to heat and combined heat–drought stress (Table 7). Similarly at Jhansi, DBW 187, HI 8777, Local-4, Raj 4120, and PDW 274 showed the highest value for STI, MP, MSI, and MSTI under HSE and HDSE. Our result is in support of the finding of Kamrani et al. (2017) for the selection of higher-yielding and heat tolerant genotypes based on MP and STI. The genotype having a higher value of TOL and SSI indicates more reduction in yield under HSE and HDSE, thus indicating stress susceptibility (Tables 7, 8) (Nouri et al., 2011; Poudel et al., 2021). The smaller the values of SSI and TOL, the greater the value of tolerance but fails to distinguish between genotypes performing better under both environments and performing better only in under stress. The selection based on TOL and SSI favors those genotypes having lower yield potential under normal and stress environments. The susceptibility indices TOL and SSI were less effective than tolerance indices STI, MP, and MSTI to select suitable lines for drought stress (Kamrani et al., 2017).
In our study, genotype PDW 274 showed higher values of stress indices STI, MP, and MSTI and ranked in the top 5 genotypes across the environment at both locations (Tables 6, 7). The other genotypes PDW 233 and PDW 291 showed the highest stress susceptibility index (SSI) across the environments at both locations. The reliability of STI for the selection of stress-tolerant genotypes with high seed yield under normal and stressed conditions further confirms the consistency of this index in selecting for high productivity under stress conditions (Fernandez, 1992).
Conclusion
In conclusion, high genetic variation was found for agro-morphological traits in a collection of wheat lines cultivated across the Bundelkhand region of Uttar Pradesh. This means that the landraces and recently released cultivars for commercial cultivation had good potential to harness in breeding for heat and drought stress. The selection based on stress indices identified tolerant lines, which could be used to improve stress tolerance in wheat. The stable traits NGS and TKW have shown positive association with SY across environments. Genotypes Local-17, PDW 274, HI-8802 and HI-8713, DBW 187, HI-8777, and Raj 4120 are identified as stress tolerant (heat and drought) at both the locations. The combination of selection indices and consistent association with physiological traits showed the potential to increase the selection efficiency of superior genotypes.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author/s.
Author contributions
Conceptualization: HK, CS and MK. Methodology: HK. Investigation: ShK, SK and AK performed multi-year-location field experiments. Data curation: ShK. Statistical analysis: RK, VG, ShK, HK, and CS. Validation: HK and VG. Resources: MK and GP. Writing—original draft preparation: HK. Review and editing: HK, VG, SK, AS and RK. Visualization: RK. All authors contributed to the article and approved the submitted version.
Acknowledgments
The authors acknowledge BUAT, Banda, and KVK, Jhansi for providing experimental land and all necessary facilities for conducting the multi-year-location trials. Authors are equally thankful to all the Bundelkhand farmers who have preserved the seed of landraces and to Regional Research Station IARI, Indore, PAU Ludhiana, and IIWBR, Karnal for providing seeds of recent varieties for evaluation in trials.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1136455/full#supplementary-material
References
Aberkane, H., Amri, A., Belkadi, B., Filali-Maltouf, A., Valkoun, J., Kehel, Z. (2021). Contribution of wild relatives to durum wheat (Triticum turgidum subsp. durum) yield stability across contrasted environments. Agronomy. 11, 1992. doi: 10.3390/agronomy11101992
Annonymous. (2020). Annual report, IIWBR. Available at: https://iiwbr.icar.gov.in/wp-content/uploads/2022/04/Annual-Report-Hindi-2020.
Ayed, S., Othmani, A., Bouhaouel, I., Teixeira da Silva, J. A. (2021). Multi-environment screening of durum wheat genotypes for drought tolerance in changing climatic events. Agronomy. 11, 875. doi: 10.3390/agronomy11050875
Barnabas, B., Jager, K., Fehér, A. (2008). The effect of drought and heat stress on reproductive processes in cereals. Plant Cell Environ. 31, 11–38. doi: 10.1111/j.1365-3040.2007.01727.x
Challinor, A. J., Watson, J., Lobell, D. B., Howden, S. M., Smith, D. R., Chhetri, N. A. (2014). Meta-analysis of crop yield under climate change and adaptation.Nat. Clim. Change. 4, 287–291. doi: 10.1038/nclimate2153
Choudhury, D. R., Tarafdar, S., Das, M., Kundagrami, S. (2012). Screening lentil (Lens culinaris medik.) germplasms for heat tolerance. TrendsBiosci. 5, 143–146.
Daryanto, S., Wang, L., Jacinthe, P. A. (2016). Global synthesis of drought effects on maize and wheat production. PloS One 11, e0156362. doi: 10.1371/journal.pone.0156362
Delahunty, A., Nuttall, J., Nicolas, M., Brand, J. (2015). “Genotypic heat tolerance in lentil,” in Proceedings of the 17th ASA Conference (Hobart), 20–24.
Devasirvatham, V., Gaur, P. M., Raju, T. N., Trethowan, R. M., Tan, D. K. Y. (2015). Field response of chickpea (Cicer arietinum l.) to high temperature. Field Crops Res. 172, 59–71. doi: 10.1016/j.fcr.2014.11.017
Dodig, D., Zoriæ, M., Kandiæ, V., Peroviæ, D., Šurlan–Momiroviæ, G. (2012). Comparison of responses to drought stress of 100 wheat accessions and landraces to identify opportunities for improving. Plant Breed. 131, 369–379. doi: 10.1111/j.1439-0523.2011.01941.x
El haddad, N., Rajendran, K., Smouni, A., Es-Safi, N. E., Benbrahim, N., Mentag, R., et al. (2020). Screening the FIGS set of lentil (Lens culinaris medikus) germplasm for tolerance to terminal heat and combined drought-heat stress. Agronomy 10 (7), 1036. doi: 10.3390/agronomy10071036
Erdemci, I. (2018). Investigation of genotype× environment interaction in chickpea genotypes using AMMI and GGE biplot analysis. Turkish J. Field Crop 231, 20–26. doi: 10.17557/tjfc.414846
Eskandari, H., Kazemi, K. K. (2010). Response of different bread wheat (Triticumaestivum l.) genotypes to post-anthesis water deficit. Not. Sci. Biol. 2, 49–52. doi: 10.15835/nsb245002
Farooq, M., Wahid, A., Basra, S. M. A., Din, I. U. (2009). Improving water relations and gas exchange with brassinosteroids in rice under drought stress. J. Agron. Crop Sci. 195, 262–269. doi: 10.1111/j.1439-037X.2009.00368.x
Farshadfar, E., Sutka, J. (2002). Screening drought tolerance criteria in maize. Acta Biol. Hung. 50, 411–416. doi: 10.1556/AAgr.50.2002.4.3
Fernandez, G. C. (1992). “Effective selection criteria for assessing plant stress tolerance,” in Proceeding of the International Symposium on Adaptation of Vegetables and other Food Crops in Temperature and Water Stress, Aug. 13-16, Shanhua, Taiwan, pp.257–270.
Fischer, R., Maurer, R. (1978). Drought resistance in spring wheat cultivars. 1. grain yield responses. Aust. J. Agric. Res. 29, 897. doi: 10.1071/AR9780897
Gupta, P. K., Balyan, H. S., Sharma, S., Kumar, R. (2020). Genetics of yield, abiotic stress tolerance, and bio fortification in wheat (Triticum aestivum l.). Theor. Appl. Genet. 133, 1569–1602. doi: 10.1007/s00122-020-03583-3
Guttieri, M. J., Stark, J. C., De Souza, B. K. O. (2001). Relative sensitivity of spring wheat grain yield and quality parameters to moisture deficit. Crop Sci. 41, 327–335. doi: 10.2135/cropsci2001.412327x
Hussain, T., Akram, Z., Shabbir, G., Manaf, A., Ahmed, M. (2021). Identification of drought tolerant chickpea genotypes through multi trait stability index. Saudi J. Biol. Sci. 28, 6818–6828. doi: 10.1016/j.sjbs.2021.07.056
Kamrani, M., Hoseini, Y., Ebadollahi, E. (2017). Evaluation for heat stress tolerance in durum wheat genotypes using stress tolerance indices. Arch. Agron. Soil Sci. 64, 38–45. doi: 10.1080/03650340.2017.1326104
Krishnamurthy, L., Gaur, P. M., Basu, P. S., Chaturvedi, S. K., Tripathi, S., Vadez, V., et al. (2011). Large Genetic variation for heat tolerance in the reference collection of chickpea (CicerarietinumL.) germplasm. Plant Genet.Resour. 9, 59–69. doi: 10.1017/S1479262110000407
Laddomada, B., Blanco, A., Mita, G., D’Amico, L., Singh, R. P., Ammar, K., et al. (2021). Drought and heat stress impacts on phenolic acids accumulation in durum wheat cultivars. Foods. 10 (9), 2142. doi: 10.3390/foods10092142
Lesk, C., Rowhani, P., Ramankutty, N. (2016). Influence of extreme weather disasters on global crop production. Nature 529, 84–87.https://doi.org/10.1038/nature16467. doi: 10.1038/nature16467
Li, Y. F., Wu, Y., Hernandez-Espinosa, N., Peña, R. J. (2013). Heat and drought stress on durum wheat: responses of genotypes, yield, and quality parameters. J. Cereal Sci. 57, 398–404. doi: 10.1016/j.jcs.2013.01.005
Limbalkar, O. M., Singh, R., Kumar, P., Nanjundan, J., Parihar, C. M., Vasisth, P., et al. (2021). Deployment of Brassica carinata a. Braun derived Brassica juncea (L.) czern. lines for improving heterosis and water use efficiency under water deficit stress conditions. Front. Plant Sci. 12, 765645. doi: 10.3389/fpls.2021.765645
Liu, H., Searle, I. R., Mather, D. E., Able, A. J., Able, J. A. (2015). Morphological, physiological and yield responses of durum wheat to preanthesis water-deficit stress are genotype-dependent. Crop Pasture Sci 66, 1024–1038. doi: 10.1071/CP15013
Lobell, D. B., Schlenker, W., Costa-Roberts, J. (2011). Climate trends and global crop production since 1980. Science 333, 616–620. doi: 10.1126/science.1204531
Mohammadi, R., Armion, M., Kahrizi, D., Amri, A. (2010). Efficiency of screening techniques for evaluating durum wheat genotypes under mild drought conditions. Int. J. Plant Prod. 4, 11–24. doi: 10.22069/IJPP.2012.677
Mondal, S., Singh, R. P., Crossa, J., Huerta-Espino, J., Sharma, I., Chatrath, R., et al. (2013). Earliness in wheat: a key to adaptation under terminal and continual high-temperature stress in south Asia. Field Crops Res. 151, 19–26. doi: 10.1016/j.fcr.2013.06.015
Monneveux, P., Jing, R., Misra, S. (2012). Phenotyping for drought adaptation in wheat using physiological traits. Front. Physiol. 3, 429. doi: 10.3389/fphys.2012.00429
Nouri, A., Etminan, A., Teixeira da Silva, J. A., Mohammadi, R. (2011). Assessment of yield, yield-related traits and drought tolerance of durum wheat genotypes (Triticumturjidum var. durum desf.). Aust. J. Crop Sci. 5 (1), 8–16. doi: 10.3316/informit.834329535176781
Olivoto, T., Lúcio, A. D. C. (2020). Metan: an r package for multi-environment trial analysis. Methods Ecol. Evol. 11 (6), 783–789. doi: 10.1111/2041-210X.13384
Passioura, J. B. (2012). Phenotyping for drought tolerance in grain crops: when is it useful to breeders. Funct. Plant Biol. 39, 851–859. doi: 10.1071/FP12079
Porter, J. R. (2005). Rising temperatures are likely to reduce crop yields. Nature. 436, 174. doi: 10.1038/436174b
Poudel, P. B., Poudel, M. R., Puri, R. R. (2021). Evaluation of heat stress tolerance in spring wheat (Triticum aestivum l.) genotypes using stress tolerance indices in western region of Nepal. J. Agric. Res. 5, 100179. doi: 10.1016/j.jafr.2021.100179
Praba, M. L., Cairns, J. E., Babu, R. C., Lafitte, H. R. (2009). Identification of physiological traits underlying cultivar differences in drought tolerance in rice and wheat. J. Agron. Crop Sci. .195, 30–46. doi: 10.1111/j.1439-037X.2008.00341.x
Punia, S. S., Shah, A. M., Ranwha, B. R. (2011). Genetic analysis for high-temperature tolerance in bread wheat. Afr. Crop Sci. J. 19, 149–163.
Puri, R. R., Gautam, N. R., Joshi, A. K. (2015). Exploring stress tolerance indices to identify terminal heat tolerance in spring wheat in Nepal. J. Wheat Res. 7,1, 13–17.
Qaseem, M. F., Qureshi, R., Shaheen, H. (2019). Effects of pre-anthesis drought, heat and their combination on the growth, yield and physiology of diverse wheat (Triticum aestivum l.) genotypes varying in sensitivity to heat and drought stress. Sci. Rep. 9 (1), 1–12. doi: 10.1038/s41598-019-43477-z
Reynolds, M. P., Saint Pierre, C., Saad, A. S. I., Vargas, M., Condon, A. G. (2007). Evaluating potential genetic gains in wheat associated with stress-adaptive trait expression in elite genetic resources under drought and heat stress. Crop Sci. 47, 172–189. doi: 10.2135/cropsci2007.10.0022IPBS
Rosielle, A. A., Hamblin, J. (1981). Theoretical aspects of selection for yield in stress and non-stress environment. Crop Sci. 21 (6), 943. doi: 10.2135/cropsci1981.0011183X002100060033x
Sareen, S., Tyagi, B. S., Sarial, A. K., Tiwari, V., Sharma, I. (2014). Trait analysis, diversity, and genotype x environment interaction in some wheat landraces evaluated under drought and heat stress conditions. Chil. J. Agric. Res. 74, 135–142. doi: 10.4067/S0718-58392014000200002
Sattar, A., Sher, A., Ijaz, M., Ul-Allah, S., Rizwan, M. S., Hussain, M., et al. (2020). Terminal drought and heat stress alter physiological and biochemical attributes in flag leaf of bread wheat. PloS One 15 (5), e0232974. doi: 10.1371/journal.pone.0232974
Sellami, M. H., Pulvento, C., Lavini, A. (2021). Selection of suitable genotypes of lentil (lens culinarismedik.) under rainfed conditions in south italy using multi-trait stability index (MTSI). Agronomy 11 (9), 1807. doi: 10.3390/agronomy11091807
Shabani, A., Zebarjadi, A., Mostafaei, A., Saeidi, M., Poordad, S. S. (2018). Evaluation of drought stress tolerance in promising lines of chickpea (Cicer arietinum l.) using drought resistance indices. Env. Stresses Crop Sci. 11 (2), 289–299. doi: 10.22077/escs.2018.420.1079
Shokat, S., Sehgal, D., Vikram, P., Liu, F., Singh, S. (2020). Molecular markers associated with agro-physiological traits under terminal drought conditions in bread wheat. Int. J. Mol. Sci. 21 (9), 3156. doi: 10.3390/ijms21093156
Smertenko, A., Draber, P., Viklicky, V., Opatrny, Z. (1997). Heat stress affects the organization of microtubules and cell division in Nicotiana tabacum cells. Plant Cell Environ. 20, 1534–1542. doi: 10.1046/j.1365-3040.1997.d01-44.x
Talebi, R., Fayyaz, F., Naji, A. M. (2010). Genetic variation and interrelationships of agronomic characteristics in durum wheat under two constructing water regimes.Braz. Arch. Biol. Technol. 53, 785–791. doi: 10.1590/S1516-89132010000400006
Tardieu, F. (2012). Any trait or trait-related allele can confer drought tolerance: just design the right drought scenario. J. Exp. Bot. 63, 25–31. doi: 10.1093/jxb/err269
Tashiro, T., Wardlaw, I. F. (1990). The effect of high temperature at different stages of ripening on grain set, grain weight and grain dimensions in the semi-dwarf wheat ‘Banks.’. Ann. Bot. (Lond) 65, 51–61. doi: 10.1093/oxfordjournals.aob.a087908
Tyagi, P. K., Pannu, R. K., Sharma, K. D., Chaudhary, B. D., Singh, D. P. (2003). Post-anthesis dry-matter accumulation and its partitioning in different wheat (Triticum aestivum) genotypes under varying growing environments. Ind. J. Agron. 49, 163–167.
Wahid, A., Gelani, S., Ashraf, M., Foolad, M. R. (2007). Heat tolerance in plants: an overview. Environ. Exp. Bot. 61, 199–223. doi: 10.1016/j.envexpbot.2007.05.011
Wardlaw, I. F., Sofield, I., Cartwright, P. M. (1980). limiting the rate of dry matter accumulation in the grain of wheat grown at high temperature. Aust. J. Plant Physiol. 7, 387–400.
Keywords: drought stress, elite cultivars, landraces, terminal heat stress, wheat
Citation: Kumar S, Kumar H, Gupta V, Kumar A, Singh CM, Kumar M, Singh AK, Panwar GS, Kumar S, Singh AK and Kumar R (2023) Capturing agro-morphological variability for tolerance to terminal heat and combined heat–drought stress in landraces and elite cultivar collection of wheat. Front. Plant Sci. 14:1136455. doi: 10.3389/fpls.2023.1136455
Received: 03 January 2023; Accepted: 05 April 2023;
Published: 09 May 2023.
Edited by:
Ilaria Marcotuli, University of Bari Aldo Moro, ItalyReviewed by:
Hugh James Earl, University of Guelph, CanadaSajid Shokat, Nuclear Institute for Agriculture and Biology (NIAB), Pakistan
Copyright © 2023 Kumar, Kumar, Gupta, Kumar, Singh, Kumar, Singh, Panwar, Kumar, Singh and Kumar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hitesh Kumar, hiteshkmr25@gmail.com; Rahul Kumar, rahul.kumar@uconn.edu; Vikas Gupta, vikas.gupta@icar.gov.in