SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

Drought yield index to select high yielding rice lines under different drought stress severities

Anitha Raman1, Satish Verulkar2, Nimai Mandal3, Mukund Variar3, V Shukla3, J Dwivedi4, B Singh5, O Singh6, Padmini Swain6, Ashutosh Mall1, S Robin7, R Chandrababu7, Abhinav Jain8, Tilatoo Ram9, Shailaja Hittalmani10, Stephan Haefele1, Hans-Peter Piepho11 and Arvind Kumar1*

Author Affiliations

1 International Rice Research Institute (IRRI), DAPO Box 7777, Metro Manila, Philippines

2 Indira Gandhi Krishi Vishwavidyalaya (IGKV), Raipur, India

3 Central Rainfed Upland Rice Research Station (CRURRS), Hazaribag, India

4 Narendra Dev University of Agriculture and Technology (NDUAT), Faizabad, India

5 Birsa Agricultural University (BAU), Ranchi, India

6 Central Rice Research Institute (CRRI), Cuttack, India

7 Tamil Nadu Agricultural University (TNAU), Coimbatore, India

8 Barwale Foundation (BF), Hyderabad, India

9 Directorate of Rice Research, Hyderabad, India

10 University of Agricultural Sciences, Bangalore, India

11 Universitaet Hohenheim, Bioinformatics Unit, 70593, Stuttgart, Germany

For all author emails, please log on.

Rice 2012, 5:31  doi:10.1186/1939-8433-5-31


The electronic version of this article is the complete one and can be found online at: http://www.thericejournal.com/content/5/1/31


Received:3 September 2012
Accepted:24 September 2012
Published:4 October 2012

© 2012 2012 Raman et al.

licensee Springer. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Drought is the most severe abiotic stress reducing rice yield in rainfed drought prone ecosystems. Variation in intensity and severity of drought from season to season and place to place requires cultivation of rice varieties with different level of drought tolerance in different areas. Multi environment evaluation of breeding lines helps breeder to identify appropriate genotypes for areas prone to similar level of drought stress. From a set of 129 advanced rice (Oryza sativa L.) breeding lines evaluated under rainfed drought-prone situations at three locations in eastern India from 2005 to 2007, a subset of 39 genotypes that were tested for two or more years was selected to develop a drought yield index (DYI) and mean yield index (MYI) based on yield under irrigated, moderate and severe reproductive-stage drought stress to help breeders select appropriate genotypes for different environments.

Results

ARB 8 and IR55419-04 recorded the highest drought yield index (DYI) and are identified as the best drought-tolerant lines. The proposed DYI provides a more effective assessment as it is calculated after accounting for a significant genotype x stress-level interaction across environments. For rainfed areas with variable frequency of drought occurrence, Mean yield index (MYI) along with deviation in performance of genotypes from currently cultivated popular varieties in all situations helps to select genotypes with a superior performance across irrigated, moderate and severe reproductive-stage drought situations. IR74371-70-1-1 and DGI 75 are the two genotypes identified to have shown a superior performance over IR64 and MTU1010 under all situations.

Conclusion

For highly drought-prone areas, a combination of DYI with deviation in performance of genotypes under irrigated situations can enable breeders to select genotypes with no reduction in yield under favorable environments compared with currently cultivated varieties. For rainfed areas with variable frequency of drought stress, use of MYI together with deviation in performance of genotypes under different situations as compared to presently cultivated varieties will help breeders to select genotypes with superior performance under all situations.

Keywords:
Rice; Drought; Genotype x environment interaction; Drought yield index; Mixed model

Background

Rainfed lowland rice ecosystems are highly variable and unpredictable in nature ( Yoshida 1977 ). Multiple abiotic stresses such as unfavorable soil conditions, regional weather patterns, topography, pests and weeds all contribute to the complexity of the ecosystem. The worldwide harvested area of rainfed lowland rice is estimated to be 46 to 48 million hectares. Of this, 90% is in South and Southeast Asia. Rice farming in these rainfed areas is risk-prone. Yields remain low, about 1.5 to 2.5 t ha-1 in most areas. The income of most farmers is low and they are challenged by erratic yields. In Asia, about 50% of all the rice land is rainfed and although rice yields in irrigated ecosystems have doubled and tripled over the past 30 years, only modest gains have occurred in rainfed rice systems ( Fischer et al. 2003 ). The water supply in rainfed areas principally comes from rainfall. Uncertainty in the timing of rainfall and variability in its intensity and its distribution cause either flood or drought stress in rainfed lowland rice production ecosystem.

Drought is the major constraint to productivity and the cause of yield instability in rainfed lowlands. At least 23 million ha of rice area in Asia are estimated to be drought-prone ( Pandey et al. 2005 ). Severe and regular droughts affect mainly eastern India, northeast Thailand and parts of Myanmar and Laos. In India, the rainfed lowland rice area is about 20.4 million hectares, which accounts for 32.4% of the total area under rice in the country. Out of the total of 20.4 million ha of rainfed rice area, approximately 7.3 million ha of lowland area are drought-prone ( Pandey and Bhandari 2008 ). The timing of drought, early season, mid-season or terminal stage, has a major influence on how much yield loss occurs ( Fischer et al. 2003 ). Drought patterns also differ among locations and among years, and in some years no drought occurs during the growing season. Therefore, poverty reduction strategies in rainfed areas must focus on stabilizing yields, that is, on breeding varieties with improved yield under drought stress as well as good response to irrigated conditions. To a farmer’s eyes, a drought-resistant cultivar is one that yields better than any other available cultivar particularly under water-limited conditions ( Blum 2006 ). The objective of a drought-tolerance breeding program is to select varieties that outperform currently available varieties in the target population of environments (TPE). The TPE is the future set of drought-prone environments in which the varieties developed by the breeding program will be grown. The environments in the TPE vary in predictable ways such as annual rainfall patterns, toposequence, soil type and farmers practices and in unpredictable ways such as random drought or disease incidence ( Fischer et al. 2003 ). Variability among environments within TPE (locations and seasons) is particularly common in rice.

The ability of crop cultivars to perform reasonably well in drought-stressed environments is paramount for stability of production. The relative yield performance of genotypes in drought-stressed and non-stressed environments can be used as an indicator to identify drought-resistant varieties in breeding for drought-prone environments. Several drought indices have been suggested on the basis of a mathematical relationship between yield under drought conditions and non-stressed conditions. These indices are based on either drought resistance or drought susceptibility of genotypes.

Let (Yi)S denote the yield of the ith genotype under stress, (Yi)NS the yield of the ith genotype under non-stress (i.e., irrigated) conditions and yS and yNS the mean yields of all genotypes evaluated under stress and non-stress conditions, respectively. Rosielle and Hamblin (1981 ) defined stress tolerance (TOL) as the differences in yield between the stress and non-stress environments, i.e., T O L = Y i NS Y i S . Higher values of TOL indicate susceptibility of a given cultivar. Hossain et al. (1990 ) defined mean productivity index (MPI) as the average of (Yi)NS and (Yi)S. This index has an upward bias when the differences between non-stress and stress conditions are large and it favors genotypes with higher yield potential and lower stress tolerance. Higher values mean a higher rate of productivity. Mean relative performance is calculated as M R P = Y i S Y S + Y i NS Y NS and relative efficiency is given by R E I = Y i S Y S * Y i NS Y NS , Ramirez Vallejo and Kelly (1998 ) computed the geometric mean of productivity (GMP), which is the square root of the product of yield under stress and yield under non-stress: G M P = Y i S * Y i NS . This index is suitable when the breeding objective is directed toward testing performance under favorable and stress conditions, taking into consideration variability in drought intensity over environment and years. Fernandez (1992 ) defined a stress tolerance index as S T I = Y i NS * Y i S Y NS 2 , which can be used to identify genotypes that produce high yield under both stress and non-stress conditions. A high value of STI implies higher tolerance of stress. A stress susceptibility index (SSI) that assesses the reduction in yield caused by unfavorable compared with favorable environments was suggested by Fischer and Maurer (1978 ). SSI is expressed by S S I = 1 Y i S Y i NS S I , SI, the stress intensity is estimated as S I = 1 Y S Y NS . Lower SSI values indicate lower differences in yield across stress levels, in other words, more resistance to drought. A modified formula for Schneider’s stress severity index ( Schneider et al. 1997 ) is defined by Singh et al. (2011 ) as S S S I = 1 Y i S Y i NS 1 Y S Y NS . The SSSI estimates the relative tolerance for yield reduction of a genotype relative to the population mean reduction in grain yield response due to stress. Selections based on these indices were carried out by many authors ( Pantuwan et al. 2002 , Ouk et al. 2006 , Golabadi et al. 2006 , Sio-Se Mardeh et al. 2006 , Talebi et al. 2009 , Khayatnezhad et al. 2010 , Nouri et al. 2011 , Singh et al. 2011 ).

Phenology is important in determining grain yield response also because quick maturing cultivars often escape from severe stress while late maturing cultivars may be affected by terminal stress ( Singh et al. 1996 ). Research on genetic variation in grain yield of pearl millet under post flowering indicated that as much as 50% of the total variation in grain yield under stress is explained by yield potential and time of flowering. Bidinger et al. (1987 ) calculated a drought response index based on a regression model to quantify the remaining part of the variation associated with tolerance / susceptibility and to identify traits linked to tolerance. The following quadratic equations were used to develop stress indices for lines in advanced screening. Mid-season stress is defined by the regression equation Y ˆ s = a + b 1 Y c + b 2 b l + b 3 b l 2 and terminal stress by Y ˆ s = a + b 1 Y c + b 2 b l , where Y s ^ is the regression estimate of the stress yield, Ys is the measured stress yield, Yc is the non-stress control yield, bl the days to flowering under non-stress, a the intercept and b1, b2 and b3 the regression coefficients. The drought index is then given by Y s Y ˆ s S . E . Y ˆ s , where S.E. indicates the standard error. Attention is focused only on those cultivars with indices of less than −1.3 or more than +1.3. These represent the upper and lower 10% of the normal distribution of the indices. Following Bidinger et al. (1987 ), Slim and Saxena (1993 ) used a DRI to describe the response of individual chickpea genotypes to drought conditions and fitted multiple regression of stressed grain yield on unstressed grain yield and days to flowering: Y ˆ 0 = a b F + C Y i . Here, Y 0 ^ is the regression estimate of yield under drought, Yi is the yield potential, F is the number of days to flowering under non-stress, b and C are the regression coefficients and a is the intercept. The drought response index is then given by D R I = Y 0 Y ˆ 0 S . E . Y ˆ 0 , where S.E indicates the standard error. ( Garrity and O’Toole 1995 ) proposed a field screening method for reproductive phase drought resistance in rice based on the drought index suggested by Bidinger et al. (1982 ) to adjust for genotypic variation in phenlogy before genotypic difference in drought tolerance is estimated.

Under severe stress, yield reduction in rice is 65-85% compared with that in non-stress conditions ( Kumar et al. 2008 ). Rice is particularly sensitive to drought stress during reproductive growth, even under moderate drought stress ( Hsiao 1982 , O’Toole 1982 ). In rice, moderate stress can be broadly characterized by a 31 to 64% loss in grain yield as compared with non-stress conditions ( Kumar et al. 2008 ). The objective of our study is to develop a drought yield index (DYI) that takes into account yield under both moderate and severe drought stress for the identification of breeding lines with superior performance over current cultivated varieties. The mean yield index (MYI) is proposed for rainfed areas where, in different years or during different growth stages of the crop in the same year, water availability fluctuates between a normal supply of water due to favorable rain and the occurrence of mild and moderate to severe drought stress based on the duration of days without rain during the cropping season. The MYI enables breeders to identify breeding lines with superior performance over current varieties under all situations.

Methods

The genotypes evaluated in this study were 134 advanced breeding lines of 105–120 days of maturity duration generated from crosses of popular high-yielding varieties with a diverse array of donors for drought tolerance, and traditional drought-tolerant landraces. The entries were tested under a diverse set of conditions that ranged from favorable irrigated conditions (non-stress) to conditions with moderate to severe reproductive-stage drought in drought-prone eastern India. The lines were tested at three locations, Raipur, Hazaribag and Faizabad, in relatively similar shallow rainfed drought-prone ecosystems, from 2005 to 2007. At all sites, the experiments were planted under irrigated control and reproductive-stage drought stress. Experiments were planted in alpha lattice designs with three replicates. From the set of 134 genotypes, a subset of 39 genotypes (Table 1) that were evaluated in two or more years under non-stress, moderate stress and severe stress conditions in 16 environments at three sites was selected for analysis. A site-year-stress-level cross combination is referred to as a trial. The details of trials are presented in Table 2. For each drought-stress experiment, an irrigated control (non-stress) trial was planted at each of the locations. In non-stress experiments, standing water was maintained at each site from transplanting to 10 days before maturity by providing water by rain or by supplementary irrigation through a pump as and when required. The reproductive-stage drought-stress experiments were irrigated like the non-stress experiments by keeping standing water up to 28 days after transplanting. Thereafter, the stress fields were drained to allow them to dry and for stress to develop. The stress experiments were not provided with any supplemental irrigation after drainage even if the stress was very severe. This has been reflected in the difference in mean yield of the stress trials at different sites (Table 2). Some trials did not achieve the targeted stress levels due to rains during the wet season and were not exposed to the desired severity of drought stress. Therefore, after harvest, each stress trial was classified for observed stress intensity as moderately stressed if the yield reduction compared to the irrigated control trial was 31–65%, and severely drought-stressed if the yield reduction surpassed 65%.

Table 1. Genotype mean yield and index values of breeding lines for the drought yield index (DYI)

Table 2. Description of shallow rainfed lowland experimental sites

Statistical analysis

Yield is the final trait used to assess drought resistance of genotypes. The definition of drought resistance in terms of yield is linked to the drought stress scenarios presented by the TPE ( Blum 2006 ). To study the yield performance of genotypes in a TPE, breeders evaluate their germplasm in multi-environment trials (METs). Testing under a range of environments, in particular drought stress levels, produces genotype x stress-level interactions. Multi-environment trials are often analyzed using mixed models ( Annichiarico 2002 , Smith et al. 2005 ). Data generated from METs are highly unbalanced due to variation in the entries tested from year to year, variation in the number of replicates at each site and variation in the sites that are included from year to year. In addition, drought-susceptible varieties may contribute to missing values when exposed to severe stress. Missing values, parameter estimation and prediction of genotype performance are effectively handled by mixed models ( Piepho and Möhring 2006 ). When the mixed model and its error structure are complex it is worthwhile exploring different data transformations such that the transformation leads to a simpler analysis of the data ( Piepho 2009 ).

Line means for each trial were first estimated fitting a mixed model by the REML algorithm taking lines as fixed and replicates and blocks within replicates as random. A conventional combined analysis of variance across years and sites was then done on log-transformed line means with all effects except stress level and genotype set to random. The performance Ylijk of the ith genotype at the jth site within the kth year at the sth stress level is modeled as:

Y lijk = μ + s i + g i + s g li + b r jk + s b r ljk + g b r ijk + e lijk (1)

where Yijk is the log-transformed yield, μ is the overall mean, si is the effect of the sth stress level, gi is the effect of the ith genotype, sgli is the interaction of the ith genotype at the lth stress level, brjk is the interaction of the jth site with the kth year, sbrljk is the interaction of the lth stress level, jth site and kth year, gbrijk is the interaction of the ith genotype, jth site and kth year and elijk is a residual error. A separate residual error variance was fitted to each stress level. The least squares means for stress x genotype interactions were computed and pair wise comparisons involving the same genotype were made. The model was fitted using the MIXED procedure in SAS ( Littell et al., 2006 ).

We propose the following drought yield index:

D Y I = Y i NS Y i S G NS G S

where the subscripts NS and S refer to control and stress, respectively, Yi is the mean yield of a genotype on the untransformed scale, while G is the geometric mean across genotypes. Since geometric means are used, the index is commensurate with a multiplicative model on the original scale, which transforms to an additive model on the logarithmic scale. As GNS and Gs are the same for all genotypes, DYI is the ratio Y i NS Y i S that determines the ranking of genotypes. The ratios were calculated after back transforming the log-transformed means to the original scale. This ratio translates into a simple difference on the log-scale. The index to compute then just boils down to the mean difference for a genotype between a control and stressed environment. The mean yield index (MYI) is defined as the average of predicted means across the non-stress, moderate and severe stress levels. In addition the deviation of genotype means from that of the means of popular varieties IR64 and MTU1010 under all three stress levels were calculated.

Based on the mean grain yield across trials under non-stress, moderate and severe stress conditions, conventional drought tolerance indices, namely, the stress susceptibility index (SSI), stress tolerance (TOL), stress tolerance index (STI) and geometric mean productivity (GMP), were calculated for both moderate and severe stress conditions. Correlations between these indices and with the drought yield index (DYI) were evaluated using the CORR procedure of SAS.

Results

ANOVA

Results of a conventional analysis of variance of the log-transformed data for grain yield indicated highly significant genotype x stress-level interaction (p<0.0001). The REML estimates of variance components and fixed effects are presented in Table 3. For each category, based on the AIC values (not presented), it is concluded that the model with separate error variances for each stress level was better than the one with constant residual variance. At each stress level, i.e., moderate and severe, those entries with non-significant differences between means of the non-stress and the respective stress level are the drought-resistant genotypes.

Table 3. REML estimates of variance and fixed effects parameters

Grain yield performance under non-stress, moderate stress and severe stress

Genotype mean yields range from 3.9 t ha-1 to 5.0 t ha-1 under non-stress, from 1.78 t ha-1 to 3.40 t ha-1 under moderate stress and from 0.45 t ha-1 to 2.19 t ha-1 under severe stress (Table 1). There was a high positive correlation (0.81) between the moderate stress and severe stress means and a low positive correlation between the non-stress yields with both moderate stress (0.07) and severe stress (0.004) yields. IR74371-70-1-1, DGI75, Khiradhan, IR64 and IR74371-78-1-1 were the best performers under non-stress. Of these, IR74371-70-1-1 and DGI75 yielded well at all stress levels. IR36 and PMK2 yielded low at all stress levels. IR74371-78-1-1 yielded well under non-stress and moderate stress and yield was reasonable under severe stress. ARB8 was the top performer under both moderate and severe stress, followed by IR55419-04 and ARB7. Khiradhan yielded very well under non-stress, but its yield declined by almost 50% under moderate stress and by 85% under severe stress. MTU 1010 was a moderate performer at all stress levels.

Proposed drought yield index (DYI)

The drought yield index values for the moderate and severe stress cases appear in Table 1. In each case, the ratios of the non-stress to the stress means were relatively low for those genotypes with a non-significant difference between the means of the non-stress and the respective stress levels. Annada, ARB5, ARB7, ARB8 and NDR 1098–6 were tolerant of stress. Significant differences were noted between non-stress and severe stress category means for all entries under severe stress, which implies that the performance under non-stress and severe stress was considerably different. However, the entries that had relatively closer performance under non-stress and severe stress were ARB7, ARB8 and IR55419-04.

Other drought tolerance indices

Moderate stress

ARB8 ranked first for most indices. With respect to STI, GMP, MRP and REI, the following seven lines were usually among the top seven ARB8, ARB3, DGI75, DGI307, DSL104-1, IR74371-70-1-1 and IR74371-78-1-1. The best performers for SSI and TOL were Annada, ARB2, ARB5, ARB7 and IR55419-04.

Severe stress

ARB8 ranked first for most indices. With respect to MP, STI, GMP, MRP and REI, the following lines were usually found within the top seven: DGI75, DGI307 and IR74371-70-1-1. The best performers for SSI and TOL were ARB5, ARB7, CB0-15-24 and IR55419-04 (Table 4).

Table 4. Conventional drought response indices of genotypes under moderate and severe stress conditions

Correlation between indices

There is a strong positive correlation between (a) REI, MRP, STI, GMP and MPI and among (b) DYI, TOL, SSI and SSSI whereas the correlation between indices in groups (a) and (b) is highly negative. Table 5 presents the correlations of the severe stress levels. The correlation between predicted means under moderate and severe stress levels was fairly high (0.69). A low positive correlation existed between predicted means of non-stress with that of moderate and severe stress levels respectively.

Table 5. Correlation coefficients between drought yield index and common drought indices under severe drought stress

Mean yield index

ARB 8 recorded the highest MYI value of 3.34, followed by IR74371-70-1-1, DGI 75, DGI307, IR74371-78-1-1, and IR55419-04 with MYI values of 3.30, 3.25, 3.21, 3.20 and 3.17, respectively.

Deviations from IR64 and MTU1010 means

The deviations of genotype means from means of the cultivated high-yielding check varieties (IR64 and MTU1010) are presented in Table 6. Two entries, IR74371-70-1-1 and DGI 75, showed positive deviations from IR64 under all three situations, whereas five entries, ARB3, DGI 75, DGI 307, IR74371-70-1-1 and IR74371-78-1-1, showed positive deviations from MTU 1010 under all three situations. However, many genotypes showed positive deviation in performance over both IR64 and MTU1010 under moderate and severe stress but only two genotypes, IR74371-70-1-1 and DGI 75, showed positive deviations from both IR64 and MTU1010 under all situations. ARB8, with the highest MYI however, did not show positive deviation from the performance of IR64 under irrigated situations because of its lower yield than IR64 under irrigated control situation. This indicates that the higher MYI of ARB8 is due to its significantly higher performance than IR64 under moderate and severe drought stress conditions.

Table 6. Mean yield index, deviations from IR64 means and MTU 1010 means across non-stress moderate and severe stress levels

Discussion

In our study, the significance of genotype x stress-level interaction was first investigated using an appropriate mixed model. The effect of stress level and genotypes was considered fixed because these were the only levels and genotypes of interest. The objective is to test the significance of genotype x stress interaction across all potential drought-prone sites and years. Interest is not specifically in those sites and years but the target population from which they were drawn. Therefore, years and sites were considered as random. When treatments are evaluated at different stress levels, they tend to have different levels of error variation. Therefore, heterogeneity of error variances at different stress levels was incorporated. The least squares means for the stress x genotype interaction were computed after appropriately accounting for all potential sources of variation. Then, pair wise stress-level comparison of means (non-stress vs moderate stress and non-stress vs severe stress) involving the same entry was carried out to look for those entries with statistically non-significant differences between means across the two stress levels.

At the reproductive stage, yield reduction in rice is significant even with a moderate stress ( Verulkar et al. 2010 ). However, this is lower than that under severe drought at the reproductive stage. Most studies point out a negative relationship between yield potential (yield under non-stress) and yield under severe stress, but varieties with high yield potential will generally have an advantage over varieties with lower yield potential under moderate drought stress. Hence, our study aims at calculating the response index under both moderate and severe stress to identify varieties that are tolerant of both moderate and severe drought. The proposed drought yield index (DYI) is similar to ones proposed by previous authors ( Fernandez 1992 , Araghi and Assad 1998 ) based on the ratio of means under non-stress to the respective stress level (moderate or severe). However, the main advantage of the current method is that DYI is calculated as a linear function of parameter estimates within the framework of a mixed model after accounting for a significant genotype x stress-level interaction across environments and it therefore provides more precise statistical inference, which was not the case with the earlier ratio-based measures.

A high positive correlation between the moderate and severe stress means indicates that the entries that performed well under moderate stress also performed well under severe stress. STI and GMP identified DGI75 (5.13, 2.76, 1.86 t/ha) and IR74371-70-1-1 (5.10, 2.92, 1.87 t/ha) among the top three stable performers. These two, along with ARB8 (4.47, 3.35, 2.19 t/ha) and IR55419-04 (4.39, 2.96, 2.15 t/ha), were the top five. This implies that STI and GMP may be useful in identifying entries that yield well under non-stress and yield reasonably well under severe stress. The proposed DYI, TOL and SSI favor genotypes with good yield under stress. For these indices, DGI75, DGI307 or IR74371-70-1-1 do not stand out among the top performers. These indices are therefore more useful for identifying stress-tolerant genotypes that perform well in stress environments.

ARB8 is again the top performer in the case of MRP and REI as well. This is because these indices are respectively the sum and product of two ratios, (i) genotype control mean/overall control mean and (ii) genotype severe-stress mean/overall severe-stress mean. The index values are inflated if either (i) or (ii) is higher. If (ii) is high, the entry enters the set of top performers though the performance under non-stress is not the topmost and vice versa. So, these indices are not very effective in distinctively discriminating entries that perform well under both non-stress and stress.

For areas where severe stress is a recurrent phenomenon, selection of genotypes with high DYI can be useful. However, selection based on DYI may lead to the identification of genotypes with high yield in moderate or severe drought stress but not very high yield or yield equivalent to that of current cultivated varieties under normal irrigated situations. MYI that considers the performance of genotypes under different situations ranging from irrigated control to moderate drought stress and severe drought stress provides better opportunities for breeders to select genotypes with good performance under all situations. The careful study of the deviation in performance of 39 genotypes from two popular varieties (IR64 and MTU1010) revealed that the performance of only two genotypes (IR74371-70-1-1 and DGI 75) showed a positive deviation over the performance of both popular varieties under the three conditions. Even though genotypes such as ARB 8, DGI307, IR55419-04, and IR74371-78-1-1 showed a high MYI, they did not show a positive deviation in their performance under all situations due to their lower yield than IR64 under irrigated control situations. For areas with a variable rainfall pattern and with equal probability of occurrence of normal rainfall, moderate stress or severe stress situations, breeders should combine MYI with a positive deviation in the performance of genotypes under all situations in order to select a variety with better performance under all situations over current varieties. In fact, even for highly drought-prone areas, combinations of DYI and deviations in the performance of genotypes from that of current varieties will enable breeders to select genotypes with yield very near that of current cultivated varieties under irrigated normal situations in addition to their higher yield under moderate and severe drought situations.

Conclusions

Many authors have reported that selection of stable genotypes should be based on a combination of indices. For regions with severe drought problems and higher frequency of drought occurrence, our study indicated that selection based on DYI will result in the identification of lines with significantly higher performance over current cultivated varieties under moderate to severe drought at the expense of slightly lower yield under normal irrigated situations. For regions where moderate to severe drought occurs once in three to four years and normal rainfall persists over other years, selection of genotypes based on a combination of MYI and a positive response over existing high-yielding varieties under all situations will enable breeders to select genotypes with higher yield over current varieties under all situations. From among 39 genotypes, the combination of MYI and deviation identified IR74371-70-1-1 and DGI 75 as two promising high-yielding drought-tolerant cultivars.

Competing interests

Authors declare that they have no competing interests

Authors’ contributions

AK, SH oversee the evaluation of breeding lines in multi environment trials and contributed to manuscript preparation, AR, HPP contributed to data analysis and manuscript preparation, SBV, MPN, MV, VDS, JLD, BNS, ONS, PS, AKM, SR, RC, AJ, TR, SH contributed to evaluation of breeding lines in multi environment trials and preparation of the manuscript. All authors read and approved the final manuscript.

Acknowledgements

Authors thank the Generation Challenge Program (GCP), Mexico and Bill and Melinda Gates Foundation (BMGF), USA for financial support to conduct this study.

References

  • Annichiarico P (2002) Genotype-environment interactions: challenges and opportunities for plant breeding and cultivar recommendations. FAO Plant Production and Protection Paper No. 174 . Food and Agriculture Organization, Rome. OpenURL

  • Araghi SG, Assad MT (1998) Evaluation of four screening techniques for drought resistance and their relationship to yield reduction ratio in wheat. Euphytica 103:293-299 Publisher Full Text OpenURL

  • Bidinger FR, Mahalakshmi V, Talukdar BS, Alagarswamy G (1982) Improvement of drought resistance in pearl millet . In: Drought resistance in crops with emphasis on rice . International Rice Research Institute, Los Banos, Laguna, Philippines. pp 357-376 OpenURL

  • Bidinger FR, Mahalakshmi V, Rao GDP (1987) Assessment of drought resistance in pearl millet [Pennisetum americanum (L.) Leeke]. I. Factors affecting yield under stress. Aust J Agric Res 38:37-48 Publisher Full Text OpenURL

  • Blum A (2006) Drought adaptation in cereal crops: A prologue . In: Ribaut JM (ed) Drought adaptation in cereals , Haworth Food Products Press, New York, USA. pp 301-333 OpenURL

  • Fernandez GCJ (1992) Effective selection criteria for assessing plant stress tolerance . In: Kus EG (ed) Adaptation of Food Crop Temperature and Water Stress , Proceeding of 4th International Symposium, Asian Vegetable and Research and Development Center, Shantana, Taiwan. pp 257-270 OpenURL

  • Fischer KS, Lafitte R, Fukai S, Atlin G, Hardy B (eds) (2003) Breeding rice for drought-prone environments International Rice Research Institute, Los Baños, Philippines. OpenURL

  • Fischer RA, Maurer R (1978) Drought resistance in spring wheat cultivars. I. Grain yield responses. Aust J Agric Res 29:892-912 OpenURL

  • Garrity DP, O’Toole JC (1995) Selection for reproductive stage avoidance in rice, using infrared thermometry. Agron J 87:773-779 Publisher Full Text OpenURL

  • Golabadi M, Arzani A, Mirmohammadi Maibody SAM (2006) Assessment of drought tolerance in segregating populations in durum wheat. Afr J Agric Res 1:162-171 OpenURL

  • Hossain ABS, Sears AG, Cox TS, Paulsen GM (1990) Desiccation tolerance and its relationship to assimilate partitioning in winter wheat. Crop Sci 30:622-627 Publisher Full Text OpenURL

  • Hsiao TC (1982) The soil-plant-atmosphere continuum in relation to drought and crop production . In: Drought resistance in crops with emphasis on rice . IRRI, Los Baños, Philippines. pp 39-52 OpenURL

  • Khayatnezhad M, Gholamin R, Jamaati-e-Somarin S, Zaibhi-e-Mahmoodabad R (2010) Study of drought tolerance of maize genotypes using the stress tolerance index. Am-Eurasian J Agric Environ Sci 9:359-363 OpenURL

  • Kumar A, Bernier J, Verulkar S, Lafitte HR, Atlin GN (2008) Breeding for drought tolerance: direct selection for yield, response to selection and use of drought-tolerant donors in upland and lowland-adapted populations. Field Crops Res 107:221-231 Publisher Full Text OpenURL

  • Littell RC, Milliken GA, Stroup WW, Wolfinger RD (2006) SAS System for Mixed Models (second edition) . SAS Institute Inc., Cary, NC. OpenURL

  • Nouri A, Etminan A, Jaime A, da Silva T, Mohammadi R (2011) Assessment of yield, yield-related traits and drought tolerance of durum wheat genotypes (Triticum turjidum var. durum Desf.). Aust J Crop Sci 5:8-16 OpenURL

  • Ouk M, Basnayake J, Tsubo M, Fukai S, Fischer K, Cooper M, Nesbitt H (2006) Use of drought response index for identification of drought tolerant genotypes in rainfed lowland rice. Field Crops Res 99:48-58 Publisher Full Text OpenURL

  • O’Toole JC (1982) Adaptation of rice to drought prone environments . In: Drought resistance in crops with emphasis on rice . IRRI, Los Baños, Philippines. pp 195-213 OpenURL

  • Pandey S, Bhandari H (2008) Drought: economics costs and research implications . In: Serraj R, Bennett J, Hardy B (eds) Drought Frontiers in Rice: Crop Improvement for Increased Rainfed Production , World Scientific Publishing and International Rice Research Institute, Singapore. pp 3-17 OpenURL

  • Pandey S, Bhandari H, Sharan R, Naik D, Taunk SK, Sastri ASRAS (2005) Economic costs of drought and rainfed rice farmers’ coping mechanisms in eastern India, Final report . International Rice Research Institute, Los Baños, Philippines. OpenURL

  • Pantuwan G, Fukai S, Cooper M, Rajatasereekul S, O’Toole JC (2002) Yield response of rice (Oryza sativa L.) genotypes to different types of drought under rainfed lowlands. I. Grain yield and yield components. Field Crops Res 73:153-168 Publisher Full Text OpenURL

  • Piepho HP, Möhring J (2006) Selection in cultivar trials – is it ignorable? Crop Sci 46:192-201 Publisher Full Text OpenURL

  • Piepho HP (2009) Data transformation in statistical analysis of field trials with changing treatment variance. Agron J 101:865-869 Publisher Full Text OpenURL

  • Ramirez Vallejo P, Kelly JD (1998) Traits related to drought resistance in common bean. Euphytica 99:127-136 Publisher Full Text OpenURL

  • Rosielle AA, Hamblin J (1981) Theoretical aspects of selection for yield in stress and non-stress environment. Crop Sci 21:943-946 Publisher Full Text OpenURL

  • Schneider KA, Rosales-Serna R, Ibarra-Perez F, Cazares-Enriquez B, Acosta-Gallegos JA, Ramirez-Vallejo P, Wassini N, Kelly JD (1997) Improving common bean performance under drought stress. Crop Sci 37:43-50 Publisher Full Text OpenURL

  • Singh BU, Rao KV, Sharma HC (2011) Comparison of selection indices to identify sorghum genotypes resistant to the spotted stemborer Chilo partellus (Lepidoptera: Noctuidae). Int J Trop Insect Sci 31:38-51 Publisher Full Text OpenURL

  • Sio-Se Mardeh A, Ahmadi A, Poustini K, Mohammadi V (2006) Evaluation of drought resistance indices under various environmental conditions. Field Crops Res 98:222-229 Publisher Full Text OpenURL

  • Singh VP, et al. (1996) Physiology of stress tolerance in rice . Proceedings of the International Conference on Stress Physiology of Rice, 28 Feb – 5 March 1994, Lucknow, U.P. India. p 239 OpenURL

  • Slim SN, Saxena MC (1993) Adaptation of spring sown chickpea to the Mediterranean basin. I. Response to moisture supply. Field Crops Res 34:121-136 Publisher Full Text OpenURL

  • Smith AB, Cullis BR, Thompson R (2005) The analysis of crop cultivar breeding and evaluation trials: an overview of current mixed model approaches. J Agric Sci 143:1-14 OpenURL

  • Talebi R, Fayazl F, Naji AM (2009) Effective selection criteria for assessing drought stress tolerance in durum wheat (Triticum durum Desf.). Gen Appl Plant Physiol 35:64-74 OpenURL

  • Verulkar SB, Mandal NP, Dwivedi JL, Singh BN, Sinha PK, Mahato RN, Dongre P, Singh ON, Bose LK, Swain P, Robin S, Chandrababu R, Senthil S, Jain A, Shashidhar HE, Hittalmani S, Vera Cruz S, Paris T, Raman A, Haefele S, Serraj R, Atlin G, Kumar A (2010) Breeding resilient and productive genotypes adapted to drought-prone rainfed ecosystem of India. Field Crops Res 117:197-208 Publisher Full Text OpenURL

  • Yoshida S (1977) Rice . In: Alvim PdT, Kozlowski TT (eds) Ecophysiology of tropical crops , Academic Press, New York. OpenURL