Genotype-by-Environment Interaction E ﬀ ects under Heat Stress in Tropical Maize

: Spring maize area has emerged as a niche market in South Asia. Production of maize during this post-rainy season is often challenged due to heat stress. Therefore, incorporating heat stress resilience is an important trait for incorporation in maize hybrids selected for deployment in this season. However, due to the signiﬁcant genotype × environment interaction (GEI) e ﬀ ects under heat stress, the major challenge lies in identifying maize genotypes with improved stable performance across locations and years. In the present study, we attempted to identify the key weather variables responsible for signiﬁcant GEI e ﬀ ects, and identify maize hybrids with stable performance under heat stress across locations / years. The study details the evaluation of a set of prereleased advanced maize hybrids across heat stress vulnerable locations in South Asia during the spring seasons of 2015, 2016 and 2017. Using factorial regression, we identiﬁed that relative humidity ( RH ) and vapor pressure deﬁcit ( VPD ) as the two most important environmental covariates contributing to the large GEI observed on grain yield under heat stress. The study also identiﬁed reproductive stage, starting from tassel emergence to early grain-ﬁlling stage, as the most critical crop stage highly susceptible to heat stress. Across-site / year evaluation resulted in identiﬁcation of six high yielding heat stress resilient hybrids.


Introduction
Maize acreage in South Asia is increasing due to its higher profitability and increased demand in feed utilization [1]. Post-rainy spring season (February-May) maize is emerging as a potential segment for expanding the maize production niche in the region. However, the production capacity of maize during this season is often limited by the high daily temperatures [2]. Given the diverse topography of the target population of environments (TPE) for spring maize cultivation, the disparities in other climatic, biotic, and abiotic factors augment the difficulties in developing and identifying suitable maize cultivars with high and stable performance across locations/year during spring season. These disparities contribute to the differential responses of cultivars and are partitioned as genotype × environment (G × E) interaction component [3][4][5]. A clear understanding of TPE is essential for systematic testing and deployment of cultivars [6]. In general, the spring maize segment in South Asia comprises high temperatures (>37 • C) zones that can be grouped into two subregions with distinct agroecological conditions: (1) hot-dry region with relative humidity <20% and (2) warm-humid region with relative humidity >60%. In addition, in certain regions within this TPE, farmers frequently provide supplemental irrigations to the crops to manage the microclimate of the crop by increasing humidity within crop canopy and obtain better yields [7]. Frequent irrigations can increase relative humidity and reduce the canopy temperature by increasing the evapotranspiration, and thereby realize good yields even under high temperatures [8][9][10][11]. However, the cost of frequent irrigations is too high to be economically sustainable for maize cultivation during spring season.
Response of maize to drought stress and heat stress has been well documented [12,13], and water deficit in combination with high temperature can cause drastic reduction in grain yield [8,[14][15][16][17]. Though both stresses have adverse effects on maize yield, a low correspondence in maize cultivars identified individually for tolerance to drought stress and heat stress, perhaps indicates a complete or partially independent mechanism of tolerance for combination of stresses [18]. Heat stress alone, or in combination with drought, is likely to become an increasing constraint to maize production in the region of maize-dependent countries [19]. A study has shown that an increase in temperature of 2 • C would result in a greater reduction in maize yields than a decrease in precipitation of 20% [20]. Similarly, a recent study in Tanzania also indicated that increasing temperatures would result in a greater reduction in maize yields than increased intraseasonal variability in precipitation [21]. Under heat stress, the crop usually simultaneously experiences (physiological) drought at high temperature. This explains relatively high yield penalties under heat stress in comparison to drought. Unless frequently irrigated for maintaining high humidity to offset the effect of physiological drought, crops face compound effects of heat plus drought, which results in relatively higher yield losses. The intensity of heat stress is defined by the level of vapor pressure deficit (VPD), which is a function of maximum temperature (Tmax) and relative humidity at Tmax. Therefore, the maize crop exposed to heat stress at different locations may respond differently depending on the level of VPD at Tmax. These studies highlight the need to incorporate tolerance to heat plus physiological drought into maize germplasm to offset predicted yield losses due to heat stress. In the present study, we attempted to assess the extent of G × E effects under heat stress during spring season maize, characterize the key weather parameters contributing to significant G × E effects, and identify suitable maize hybrids with high and stable performance under heat stress across locations and years.

Materials
Materials involved in the present study were various precommercial hybrids developed from the advanced stage lines of the heat stress breeding pipeline of the International Maize and Wheat Improvement Center (CIMMYT)-Asia program. Parental lines of these experimental hybrids were Asia-adapted lines having founder parents derived from various pools and populations developed by the CIMMYT stress breeding program, including drought-tolerant population (DTP), La Posta Sequia, Pool 16 BN Sequia, Pool 18 Sequia, and Pool 26. They were constituted during the mid-1980s using putative drought-tolerant sources, including Tuxpeno Sequia C 8 , Latente, Michoacán 21, Suwan 1, crosses of CIMMYT populations 22, 32, 62, 64, and 66, landraces, Corn Belt hybrids, and germplasm from Thailand, Brazil, and South Africa. The details of the selection and improvement procedure of these populations are described in previous studies [22][23][24]. Parental lines of these advanced-stage Agronomy 2020, 10, 1998 3 of 18 hybrids were selected from the extensive heat screens of the program that were conducted across various sites of Asia during 2013 and 2014 spring season evaluations. Details of these parental lines involved are presented in Supplementary Table S1.
The present study was based on the performance of prerelease hybrids, selected from breeding pipeline, tested across three years, 2015-2017, in locations across South Asia during spring season (Table 1). Experimental hybrids across all the years were not common. For each year, some poor-performing entries were dropped from the trial and new advanced stage experimental hybrids were included in the next trial from the breeding pipeline at the final stage of evaluations. In the year 2015, 19 advance stage hybrids were evaluated across 8 locations, followed by 32 hybrids during 2016 across 16 locations, and 27 hybrids in 2018 across 11 locations. Similar to the test entries, locations were also not the same each year and were chosen from a pool of locations based on the target population of the environment. However, the entries were balanced across years and locations, and their performance was evaluated against the best available commercial checks for the location.

Field Trials across Locations under Heat Stress
Field trials were conducted during spring season at multiple locations in India, Nepal, Pakistan, and Bangladesh (Table 1). Trials were laid-out following an alpha lattice design in 4 m row plots with two replications using a recommended plant spacing of 75 cm × 20 cm, which represents an optimum Agronomy 2020, 10, 1998 4 of 18 plant population of 66,666 plants ha −1 . Recommended agronomic practices were followed to maintain a good plant stand, and to keep trials free from weeds and other biotic pressure, including local diseases and insect pests. Entries across locations within a year were balanced and their performance was evaluated against popular commercial hybrids for spring season. The number of entries tested across years and the locations are presented in Table 1. Entries were evaluated for days to 50% anthesis, days to 50% silking, and grain yield under heat stress. Grain yield was recorded on per plot basis and converted to tones per hectare after adjusting for kernel moisture content at 12.5%.

Characterization of Locations and Weather Parameters
All the trials were conducted during the spring season (planted during March), at selected locations that normally record average maximum temperatures above 35 • C during most part of the crop cycle. Daily weather variables were collected during the trial period and used for characterization of the heat stress environments at different locations. Daily weather data were recorded during the trial period (from planting to maturity) on (i) temperature ( • C), including maximum (Tmax) and minimum temperature (Tmin), and (ii) relative humidity (%) at Tmax and vapor pressure deficit (VPD), which was calculated using Tmax and relative humidity at Tmax as the absolute difference between saturation vapor pressure and actual atmospheric vapor pressure [25], as follows: where Tmax represents the maximum temperature (in • C) and RH represents relative humidity (%). The climatic variables were then averaged over fortnightly intervals (designated at quarters, Q) starting from the date of planting until dough stage of the crop, for further analysis and environmental characterization. Daily averages of the weather variable from each fortnight were designated as Q1 to Q7 representing the prevailing weather condition at major crop growth stages, starting from seedling emergence (VE) and ending at dough stage (R4). Q1 to Q3 represented the vegetative phase of the crop, Q4 to Q6 represented the reproductive or flowering phase, and Q7 represented the early grain-filling phase. Hence, in total, 28 environmental variables spanning across seven fortnights of crop growth stages were used in the current study for factorial regression analysis. In addition, based on the mean yield of the trials, the environments were grouped as high yielding (>6 t ha −1 ), moderate yielding (3-5.9 t ha −1 ), and low yielding (<2.9 t ha −1 ) and across-environment analysis was performed within each type of environment. Relative efficiency of indirect selection across these three major environmental classifications were also followed using the procedure outlined by Vanessa et al. (2012).

Statistical Analysis
The residual maximum likelihood (REML) approach was used to analyze grain yield from a single site trial dataset, treating incomplete blocks, replications as random factors, and entry as fixed. Single location repeatability of the trials was computed using the genotypic variance estimates (σ 2 g) and single location residual (σ 2 ε ) as Trials with good repeatability (ω 2 > 0.50) were then subjected to across-location analysis and variance components estimated within each year using the model where µ denotes the overall mean; g i the genetic effect of genotype i; e j the effect of location j; g eij the interaction between genotype i and location j; r k (e j ) the effect of the replication k nested in the location j; b l [re] kj the incomplete block l nested in the replication k and location j; and ε jklm the residual effect of the plot m nested in block l, replication k, and location j. All factors were considered as random effects for estimating the variance components in the combined analysis for each environmental group. In each type of the environmental groups, heritability (H) of grain yield was calculated as where e denotes the number of locations and r the number of replicates, σ 2 g is the genotypic variance, σ 2 ge is the variance due to genotype × environment, and σ 2 ε is the residual variance. To assess the gains made in each environmental class over the period of three years, the average performance of the best 10 percent hybrids in each environment class (low, moderate, and high yielding environment) was compared to the mean of the commercial checks using Welch-t test for unequal sample size and unequal variances in order to determine if the differences in the means were significant.
To measure stability of the entries in each year, the related linear regression genotypic stability of Tai [26] was estimated using the software GEA-R [27]. Two parameters, α (determining the hybrid response to environmental effect) and λ (deviation from linear response), were used for determining the stability of the hybrid. A hybrid with average stability will have ∝= 0 and λ = 1, and a perfectly stable hybrid will have ∝= −1 and λ = 1 [28]. These parameters can also be represented in two orthogonal axes to depict the entries, which would help in determining the stability of the entries [27,28].
To understand the effect of environmental variables on the genotype × environment interaction, a factorial regression was performed on the sampled trial locations, where data on required weather variables were available. Daily averaged weather variables from each fortnight during the trial period were considered as the environmental covariate and grain yield as the dependent variable. Factorial regression approximates the interaction effects, and the model has been defined in previous studies [29][30][31] to interpret the interaction. The general factorial regression model for which the GEI includes environmental covariates [27] is as follows: where Y ij is the yield of the ith genotype (I = 1 to I) in the jth environment (j = 1 to J); µ is the grand mean; g i and e j are the genotype and environment deviations from the grand mean, respectively; z in are the environmental covariates; ζ jn are the genotype factors; N is the number of environmental covariates; and ε ij is the error term. A forward regression was used in the study for identifying the most important environmental covariable.

Heat Stress at Different Testing Locations
Selected locations were sampled for climatic variables during the trialing season across the test environments. Daily weather data were averaged over fortnightly intervals representing vegetative (Q1-Q3), reproductive or flowering (Q4-Q6), and early grain-filling (Q7) stages (Figure 1).

Performance of Hybrids across Locations
Significant variation was observed among the heat stress locations in terms of performance of test hybrids (Table 2). Mean yield of the trial was used for classifying the locations as low yielding (<3 t ha −1 ), moderate yielding (3-6 t ha −1 ), and high yielding (>6 t ha −1 ).
Only few of the locations tested (E12 in year 2015; E2, E9, E18 in 2016; and E22, E23 in 2017) had a mean trial yield of less than 3 t ha −1 (Table 1). Combined analysis of grain yield for each tested genotype across the years showed significant variation among genotypes (p < 0.01) ( Table 2). Genotypic variance accounted for 7, 21, and 14% of the total phenotypic variance for 2015, 2016 and 2017 year trials, respectively. Genotype × location variance was highly significant (p < 0.0001) and accounted for more than 30 percent of the total phenotypic variance for grain yield across the locations tested in each year. Consequently, the performance of genotypes was not similar across the location, with few genotypes showing clear crossover interactions. Performance of genotypes tested across all the locations are presented in Figure 2.
The trial repeatability for grain yield in individual locations tested was high (>50%) ( Table 1) except for one location in India (Bejjanki) that was dropped from further analysis. The heritability estimate across locations was 40% in 2015 and more than 70% across locations in 2016 and 2017. Correlation coefficient estimate between trial repeatability and trial grain yield was nonsignificant but positive (r = 0.30). Considering high G × E interaction effects and a large variation in the trial yields obtained in the locations tested (Table 2), a stratified ranking approach was used to identify the best cultivars in each year. The entries, in each location were ranked based on the grain yield, and the entry that performed consistently in top one-third of the total entries tested was considered best suited for next stage of evaluation ( Figure 3). Tai stability parameters was estimated for each genotype to determine the stability of these entries (Figure 4). In the year 2015, three hybrids (CAH1516, ZH141592, and VH112887) were identified as the best entries; these hybrids consistently ranked in top one-third fraction of the total entries tested in at least 3 or 4 of the 7 locations tested (43-57% of locations tested). Grain yield of the test entries ranged from 2.5 t ha −1 -11.0 t ha −1 across the tested locations. Among these, ZH141592 was identified with average stability. In addition, another entry, CAH153, was identified as relatively more stable in this year, though not top ranking across all locations, but this entry also had good yields across the locations (Table 3). While none of the identified entries surpassed the best commercial check in the test location, specific entries were identified in each of the test locations that had yields comparable to the best commercial check (Figure 3). In the year 2016, three hybrids consistently ranked in the top one-third fraction of the tested entries. CAH1432, CAH1719, and ZH15440 ranked top in more than 60% of the tested locations (9 of 15 locations tested). Amongst them, CAH1719 was identified as the most stable entry, whereas CAH1432 and ZH15440 were identified as unstable entries. Average yield of the entry CAH1719 across the test location ranged between 1.6 and 1.9 t ha −1 in low yielding locations and 7.3 t ha −1 in high yielding locations, respectively (Table 3). This hybrid also performed significantly better or on par with the commercial hybrids in most of the moderate and high yielding locations. However, in low yielding locations, the yield of CAH1719 was significantly lower than the commercial checks ( Figure 3). LSD-Least significant difference.
Agronomy 2020, 10, 1998 9 of 18 0.0001) and accounted for more than 30 percent of the total phenotypic variance for grain yield across the locations tested in each year. Consequently, the performance of genotypes was not similar across the location, with few genotypes showing clear crossover interactions. Performance of genotypes tested across all the locations are presented in Figure 2.   Table 1.

of 18
Agronomy 2020, 10, x FOR PEER REVIEW 11 of 20 † CAH1714 1.60 ± 0. 30 5.24 ± 1.36 7.56 ± 0.39 † Entries identified as stable based on Tai stability parameters.   In addition to these, entries CAH153, CAH1715, ZH15434, and ZH141593 were identified with more than average stability. However, the performance of CAH153 and CAH1715 surpassed the other two entries as they ranked in top one-third of the fraction in over 8 of the 15 locations tested. In 2017 trials, two entries (ZH15400 and CAH1714) were identified as adaptable entries and they were at least on par with the best commercial check in the locations tested (Figure 3). These  In addition to these, entries CAH153, CAH1715, ZH15434, and ZH141593 were identified with more than average stability. However, the performance of CAH153 and CAH1715 surpassed the other two entries as they ranked in top one-third of the fraction in over 8 of the 15 locations tested. In 2017 trials, two entries (ZH15400 and CAH1714) were identified as adaptable entries and they were at least on par with the best commercial check in the locations tested ( Figure 3). These entries ranked in top one-third of the entries tested in 9 of the 11 locations tested (82 percent). However, Tai stability parameters indicate ZH15400 to be highly unstable, while CAH1714 was identified as having average stability. Another entry, ZH15333, was identified as the most stable and ranked in top one-third of the entries tested in 6 of the 11 locations tested. These identified entries out yielded or were on par to the best commercial check in the trial (Figure 4).
To assess the progress being made in terms of grain yield within the heat stress breeding programs, a post hoc t-test was used to compare the average yield of the top 10% entries with the performance of the checks. The locations were clustered based on the trial mean yield at an interval difference of 1 t ha −1 and the mean yield of the top ranking (top 10%) entries in these locations were compared to the mean performance of the checks in these locations. The average performance of these top entries in comparison to the checks performance along with the significance is depicted in Figure 5. The average yield of the best test entries often out yielded the performance of checks, particularly in all moderate yielding locations. However, these differences were not substantial in high yielding locations.
Agronomy 2020, 10, x FOR PEER REVIEW 13 of 20 entries ranked in top one-third of the entries tested in 9 of the 11 locations tested (82 percent). However, Tai stability parameters indicate ZH15400 to be highly unstable, while CAH1714 was identified as having average stability. Another entry, ZH15333, was identified as the most stable and ranked in top one-third of the entries tested in 6 of the 11 locations tested. These identified entries out yielded or were on par to the best commercial check in the trial (Figure 4). To assess the progress being made in terms of grain yield within the heat stress breeding programs, a post hoc t-test was used to compare the average yield of the top 10% entries with the performance of the checks. The locations were clustered based on the trial mean yield at an interval difference of 1 t ha −1 and the mean yield of the top ranking (top 10%) entries in these locations were compared to the mean performance of the checks in these locations. The average performance of these top entries in comparison to the checks performance along with the significance is depicted in Figure 5. The average yield of the best test entries often out yielded the performance of checks, particularly in all moderate yielding locations. However, these differences were not substantial in high yielding locations.

Contribution of Environmental Factor in G × E Interaction Effects
A sample of test locations (Table 1) was selected from the pool of locations in each of the years tested and the daily averages of the environmental factors were determined over an interval of two weeks throughout the trialing season in that location. These environmental factors were then treated as environmental covariates for fitting the

Contribution of Environmental Factor in G × E Interaction Effects
A sample of test locations (Table 1) was selected from the pool of locations in each of the years tested and the daily averages of the environmental factors were determined over an interval of two weeks throughout the trialing season in that location. These environmental factors were then treated as environmental covariates for fitting the factorial regression model onto the grain yields of the tested entries in these locations. In the year 2015, environmental factors from 7 locations were used in the analysis, while from the 2016 and 2017 year trials, 6 and 5 locations, respectively, were used for the analysis. The genotype × location variance in each year for the sample location accounted for 11 to 33 percent of the total phenotypic variation for grain yield (Table 4) across the three years. Table 4. Sums of squares and percent variation estimates using environmental covariates to partition genotype × environment interaction variance for grain yield into genotype × environmental covariate interactions, using factorial regression.

Source of Variation
Year Factorial interaction revealed several significant GE interaction components for grain yield in each year. The model for interaction component in the year 2015 consisted of the five interactions with genotype: Genotype × Q5_Tmin, Genotype × Q6_VPD, Genotype × Q2_VPD, Genotype × Q7_Tmax, and Genotype × Q4_Tmax, cumulatively accounting for 88 percent of the GE interaction variation with 70 degrees of freedom (p-value ranging from 0.002 to 2.58 × 10 −9 ). In 2016, the genotype × environmental interaction model consisted of genotype interaction with Q7_Tmax, Q4_Tmax, Q3_RH, and Q4_RH, these four factors accounted for 86 percent variation with 136 degrees of freedom. In 2017 Genotype × Q7_RH, Genotype × Q6_Tmax, and Genotype × Q6_RH cumulatively accounted for 67 percent of variation with 87 degrees of freedom (Table 4). Across the years, Tmax along with reproductive stage, RH, and VPD accounted due to most of the variation for genotype × environmental interaction ( Figure 6).
consisted of genotype interaction with Q7_Tmax, Q4_Tmax, Q3_RH, and Q4_RH, these four factors accounted for 86 percent variation with 136 degrees of freedom. In 2017 Genotype × Q7_RH, Genotype × Q6_Tmax, and Genotype × Q6_RH cumulatively accounted for 67 percent of variation with 87 degrees of freedom (Table 4). Across the years, Tmax along with reproductive stage, RH, and VPD accounted due to most of the variation for genotype × environmental interaction (Figure 6).

Discussion
The locations used in the current study are representative of spring maize cultivation zones in South Asia and are characterized as heat stress prone with recorded average daily Tmax > 35 °C and Tmin >23.0 °C during most part of spring season. The purpose of the study was to understand the weather variables in the spring environment, their contribution in heat stress intensity, and identify hybrids best suited for cultivation in this region. A comparison of weather parameters across the years and locations, for the duration of the crop period indicated that in general, daily maximum temperature (Tmax) was highest during most part of the flowering stage. This stage is also reported as the most susceptible growth stage for maize to high temperature [32][33][34], causing significant losses to grain yield, mostly by impacting pollen viability [35,36]. Distinct differences were observed in mean grain yield of the trial, within each year ranging between low (<3 t ha −1 ), moderate (3-5 t ha −1 ), and high yielding (> 6 t ha −1 ) locations in spring season. Broad-sense heritability of the trials was high (>0.50) across all these managed stress locations. Significant effects of stress on repeatability and grain yield has been reported in earlier studies [37,38], suggesting that as stress increases (causing lower grain yield), the repeatability of the trial is also affected. Contrary to these findings, in this study the correlation coefficient estimates between trial grain yield and repeatability were positive but nonsignificant across all the three groups of yielding locations (low, moderate, and

Discussion
The locations used in the current study are representative of spring maize cultivation zones in South Asia and are characterized as heat stress prone with recorded average daily Tmax > 35 • C and Tmin >23.0 • C during most part of spring season. The purpose of the study was to understand the weather variables in the spring environment, their contribution in heat stress intensity, and identify hybrids best suited for cultivation in this region. A comparison of weather parameters across the years and locations, for the duration of the crop period indicated that in general, daily maximum temperature (Tmax) was highest during most part of the flowering stage. This stage is also reported as the most susceptible growth stage for maize to high temperature [32][33][34], causing significant losses to grain yield, mostly by impacting pollen viability [35,36]. Distinct differences were observed in mean grain yield of the trial, within each year ranging between low (<3 t ha −1 ), moderate (3-5 t ha −1 ), and high yielding (> 6 t ha −1 ) locations in spring season. Broad-sense heritability of the trials was high (>0.50) across all these managed stress locations. Significant effects of stress on repeatability and grain yield has been reported in earlier studies [37,38], suggesting that as stress increases (causing lower grain yield), the repeatability of the trial is also affected. Contrary to these findings, in this study the correlation coefficient estimates between trial grain yield and repeatability were positive but nonsignificant across all the three groups of yielding locations (low, moderate, and high). This could be related to the lower sample size and/or low variation among the estimated repeatability parameters across the managed stress trials evaluated. Across-location heritability estimates within each year were moderate (0.40 in 2015) to high (>0.70 in 2016 and 2017); however, a few strong crossover performances were also observed among entries, particularly between locations in Bangladesh in comparison to locations in India and Pakistan during the three years.
While the trials were evaluated in spring season across South Asia, higher grain yields of test entries were recorded at some of the locations, such as Jessore (Bangladesh) in 2015 (>10 t ha −1 ). It was also observed that, in general, locations in Bangladesh had higher trial mean yields across years as compared to other locations. This could possibly be attributed to the warm humid (high relative humidity) climatic conditions prevailing in this region. Higher trial yields were also observed in some of the locations of India (such as Sabour) that has similar climatic conditions and topography to locations tested in Bangladesh. In maize, the advantage of humid over drier atmosphere, resulting in higher dry matter production and eventually higher yields, have been reported [39]. Frequent irrigations even during high temperatures reduces the canopy temperature and thereby cuts yield losses due to high atmospheric temperature [9,40]. Though resource intensive and unsustainable in long run, this is also a common practice currently followed in several parts of spring maize production areas of South Asia.
Studies have shown that environmental covariables other than the one primarily studied could also play a significant role in the performance of entries [41], and often breeding designs do not take into account these factors for evaluating entries and selecting test locations. Factorial regression models could help delineate environmental variables that could influence crop traits [29,31]. As required, this study involved testing of entries in a wide range of test locations (5-7) sampled from a larger set of test locations commonly having high maximum temperatures during the cropping season (spring) across South Asia. As expected, a strong G × E interaction was observed in the test locations with some entries exhibiting significant crossover performance for grain yield across the three years. Among environmental variables that influence the crop yield, relative humidity (RH) and vapor pressure deficit (VPD), apart from maximum temperature, are the two major parameters that were found to contribute to the significant variation for grain yield among the test entries, particularly during flowering and the early grain-filling stage of the crop across all three years.
Reproductive stage has been reported as the most critical stage of the maize crop to high temperature stress [19,[42][43][44] and premonsoon season maize particularly in parts of South Asia is most vulnerable to severe heat stress during anthesis and early grain-filling stages [45]. The effect of heat stress is further exacerbated by the low levels of moisture in the atmosphere (low and high VPD), causing pollen to dry down faster and lose its viability, thereby directly impacting the yields [36,46]. The analysis only identified minimum temperature as an important variable for the year 2015 dataset, while in 2016 and 2017 this weather parameter was nonsignificant. It was also observed that the effects of minimum temperature (>30 • C) could potentially reduce biomass and crop yields, as night-time temperature could directly impact night-time respiration rates of plants [32]. The inability of this analysis to identify minimum temperature across years as an important variable might possibly be because the environments sampled for the current study did not pose enough variation for minimum temperature to be expressed in the observed yield losses, particularly in 2016 and 2017. Locations with higher relative humidity (such as parts of Bangladesh and locations in Bihar, India) recorded higher grain yields in comparison to parts where the atmosphere was relatively much drier due to low humidity, such as Hyderabad (India) and Sahiwal (Pakistan), even under a comparable temperature regime across the locations during the crop growth. Therefore, it is pertinent to consider these additional environmental factors before selection of test locations and entries.
Spring maize cultivation area within South Asia consists of locations spread across varied topographies, and therefore the grain yield of the entries tested was not consistent across these locations, with few entries exhibiting substantial crossover performance. The variance estimates of genotype × location interaction was also high and at least more than twice the estimated genotypic variance across the three years, suggesting a high influence of the environment on the genotypic performance. Overall, the variance due to genotype × environment in the current study was estimated at 13% of the total phenotypic variance. High genotype × environment interaction, although suggesting more test locations are needed to identify more stable hybrids, also presents an opportunity to select entries with specific adaptations [47]. The current study identified CAH153 in 2015; CAH153, CAH1719, and CAH1715 in 2016; and ZH15333 in 2017 as the most stable entries. A combination of parametric (Tai stability) and nonparametric (stratified ranking) approach identified entry ZH141592 in the year 2015 as a potential entry for advancement to the next stage of testing, as it had good stability and also performed well in low (3.48 t ha −1 ), moderate (5.65 t ha −1 ), and high (8.23 t ha −1 ) yielding locations. In 2016, entries CAH153, CAH1719, and CAH1715 performed well in many locations tested and had good stability and yield levels in both moderate and high yielding locations. In 2017, entry ZH15333 performed well (6 of the 11 locations tested) and had more than average stability, while CAH1714 in 2017 performed well in most of the test locations (9 of 11 locations) but had only average stability. A comparison of the average grain yield performance of the top 10 percent entries in each location and the commercial checks clearly indicates that the products from the heat stress breeding pipeline performed better than or were at least on par to the commercial checks used in the study. Interestingly, the performance of the top entries in comparison to commercial checks were clearly higher in moderate yielding locations (3 to 7 t ha −1 grain yield) in comparison to high yielding locations (optimal with >7 t ha −1 grain yield), and low yielding locations (severe stressed <3 t ha −1 grain yield). These results are comparable to a study by [48], where 41 hybrids from a stress breeding program showed distinct advantage when compared to 42 commercial hybrids from a private seed company at different yielding locations and levels, suggesting no yield drags of the stress breeding pipeline products, even under high yielding environment if an appropriate selection strategy is chosen. In fact, with adequate weighing of managed stress-prone location performance of entries in the breeding pipeline, the yield levels of entries could be marginally increased even in severely stressed low yielding environments [49].

Conclusions
Findings of the present study clearly showed that spring season environments for maize in South Asia have varied weather conditions and there are strong genotype × environment interaction effects. The factorial regression analysis provided evidence for the effects of relative humidity and vapor pressure deficit as key determinants for yields under heat stress. This suggests the need to delineate the test locations based on these key weather parameters (relative humidity and vapor pressure deficit) while selecting heat stress resilient genotypes for spring season cultivation. Further, this study also provides evidence that effective selection of entries under heat stress in the breeding pipeline may not cause any yield drags under a high yielding condition, but rather contribute in stable performance across variable heat stress regimes.