Next Article in Journal
Recovery of 2R.2Sk Triticale-Aegilops kotschyi Robertsonian Chromosome Translocations
Previous Article in Journal
Crop Rotation Compared with Continuous Rye Cropping for Weed Biodiversity and Rye Yield. A Case Study of a Long-Term Experiment in Poland
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of Different Multivariate Statistical Methods for Screening the Drought Tolerant Genotypes of Pearl Millet (Pennisetum americanum L.) and Sorghum (Sorghum bicolor L.)

by
Mohammadreza Negarestani
1,
Enayatollah Tohidi-Nejad
1,2,*,
Gholamreze Khajoei-Nejad
1,2,
Babak Nakhoda
3 and
Ghasem Mohammadi-Nejad
1,2
1
Department of Agronomy and Plant Breeding, Shahid Bahonar University of Kerman, P. O. Box 76169-133 Kerman, Iran
2
Research and Technology Institute of Plant Production (RTIPP), Shahid Bahonar University of Kerman, P.C. 7616913439 Kerman, Iran
3
Department of molecular physicology, Agricultural Biotechnology Research Institute of Iran, Agricultural Research Education and Extension Organization (AREEO), P. O. Box 19395-1113 Karaj, Iran
*
Author to whom correspondence should be addressed.
Agronomy 2019, 9(10), 645; https://doi.org/10.3390/agronomy9100645
Submission received: 4 September 2019 / Revised: 7 October 2019 / Accepted: 11 October 2019 / Published: 17 October 2019

Abstract

:
Drought is the main limiting factor of crops production in major regions of the world. Forage plants such as pearl millet and sorghum are drought tolerant and suitable for arid regions to grow. In this study, for selecting and introducing the best drought tolerant genotypes, seven pearl millet and five sorghum genotypes in three locations (Kerman, Jiroft, and Bardsir in Kerman Province) were studied with different climatic conditions. The experiments were conducted for three consecutive years of cultivation (2016, 2017, and 2018) under regularly irrigated conditions and the 100% (I100, full irrigation) plant water requirement and under water deficiency, (50%) plant water requirement (I50) in two randomized complete block designs in triplicate separately at each location. Eight drought tolerance/susceptibility indices were used. Multivariate factor analysis (FA) and stress tolerance score (STS) methods were employed to compare the most drought tolerant genotypes in each region. The STS method was more efficient and effective than the FA method for detecting the most drought tolerant genotype with any number of used genotypes. Based on the results of drought tolerance indices and STS, the IP13150 and IP13151 genotypes of pearl millet and speed feed genotype of sorghum in Kerman, Jiroft, and Bardsir respectively, were introduced as the most drought tolerant genotypes for three consecutive years of cultivation.

1. Introduction

Drought stress is the most important factor limiting the growth of plants all over the world; hence, reduced growth due to drought stress is more than other environmental stresses significantly [1].
Pearl millet (Pennisetum americanum L.) and forage sorghum (Sorghum bicolor L.) are adapted to hot arid climates [2] and drought stress, which are appropriate for regions with irrigation water shortage is the main cause of limitation in crop production [3,4]. Due to its low water requirement, low dependence on the soil type as well as high yield and quality of forage, pearl millet can be considered as a suitable strategy for feeding livestock [5]. Different mechanisms have been proposed for drought tolerance in pearl millet and sorghum, including strategies against drought stress, improving drought tolerance, and premature plants. It is important to know how each genotype tolerates drought stress. Researches have revealed that millet can escape from the impacts of drought stress through early flowering and lack of synchronicity in tillers development [6,7].
Bruck et al. (2000) [8] reported the low crop harvest index in pearl millet due to drought stress. Different indexes have been reported for identifying the tolerance and response of cultivars to stress based on yield. Some researchers have suggested the tolerance (TOL) and the mean productivity (MP) indexes. High values of TOL indicate the relative sensitivity of genotypes to stress [9]. MP index is also defined as the total yield of a genotype under stress and non-stress conditions. The stress tolerance index (STI) is considered as a criterion for selecting drought tolerant cultivars; high values of STI indicate high tolerance and high yield potential [10].
Geometric mean productivity (GMP) is another index having more potential than MP for separation of cultivars [10]. The stress susceptibility index (SSI) is another selection index. Researchers believe that genotypes with lower SSI than unit tolerate drought more. Hence, reduction in their yield under drought conditions is less than the reduced average yield of total genotypes [11]. The yield stability index (YSI) assesses the yield of a cultivar under stress conditions compared to its yield under non-stress conditions and it can be a suitable index for identifying varieties tolerant of drought stress. Consequently, it is believed that cultivars with higher YSI have better yield under both conditions [12]. The yield index (YI) only categorizes the cultivars based on yield and under stress conditions, thus, it does not detect cultivars with high yield under both stress and non-stress conditions [13]. Ouk et al. [14] have reported that the yield under drought stress is associated with three factors, including yield potential, drought tolerance, and appropriate phenology. Drought tolerance is one of the traits controlled by many genes and it is one of the quantitative traits. The main problem is to identify a standard assessment test for selecting drought tolerant genotypes [15]. The aim of this study was selecting drought tolerant genotypes in regions with different climates using all drought tolerance/susceptibility indexes under stress and non-stress conditions simultaneously.
Ghollamin et al. (2010) [16] in a study stated that on durum wheat properties under drought stress and without drought stress 82.58% of its total yield variability was justified by factor analysis [16]. Abdulshahi et al. (2012) selected drought resistant genotypes from 40 different wheat genotypes using the STS method [17].

2. Materials and Methods

2.1. Treatments

The treatments were performed to determine drought resistant genotypes, among the pearl millet cultivars including HHVBC, IP22269, HYB-1, HYB-2, ICMVO-O5222, IP13150, and IP13151 and different sorghum genotypes including S35, ICSV-112, SPEED FEED, ICSV-14001, and Pegah were selected in three regions almost different alternate in climate, including: Kerman, Jiroft, and Bardsir, whose geographical features and other characteristics are given in Table 1.
Water treatments including 100% (I100, full irrigation) plant water requirement and 50% plant water requirement (I50) were used respectively, as regularly irrigated and drought stress condition (I100, I50).

2.2. The Characteristics of Experimental Farms

The first field experiment site was conducted at the Research Farm in Faculty of Agriculture, Shahid Bahonar University of Kerman and had a hot arid climate. The second field experiment site for the study was located in the Research Farm, Faculty of Agriculture, Islamic Azad University of Jiroft and had a hot arid climate. The third field experiment site located in the Research Farm of Agriculture Organization Office of Bardsir had a temperate semi-arid. The temperature range of the three regions is given in Figure 1. This experiment was conducted in the mentioned regions during the 2016, 2017, and 2018 crop years. The soil properties of the tested sites are also given in Table 2.

2.3. The Stages of Implementation

At the first step, the studied farms were plowed and leveled using a land leveler so that the slope of the ground became almost nearly zero and the runoff problem was eliminated. Fifty kg ha−1 of P2O5 fertilizer and 70 kg ha−1 of urea fertilizer were added to the soils in three farms before implantation. In this experiment, two randomized complete block designs were used in each area (each with three replications). Each replicate consisted of 12 experimental units (plots) assigned to each genotype. Therefore, each randomized complete block design consisted of 36 experimental units and a total of 72 experimental units in total.
Each block design was irrigated with water treatment. The amount of irrigation water was determined by two separate weighing lysimeters embedded on the periphery of the farm (one for pearl millet and another for sorghum). The irrigation water was transferred to the test site through a poly vinyl chloride (PVC) pipe (75 mm) using an electric pump from the water source and the required water was given to each plot by a water volume meter. Table 3 shows the number of irrigation in the growing season and the amount of water consumed in each area.
Water treatments were applied when the plant was well established (four-leaf stage). Five lines were implanted in each experiment unit. Plants were cut from 10 cm height of soil surface for harvesting. After harvesting, the fresh forage was weighed. The growing season was considered 135 days from planting date to final harvesting with the same period for all three years. In some areas, such as Jiroft with high temperatures, three crops were harvested (first crop: 72 days after planting, the second: 120 days, and the third: 135 days after planting), and in Kerman and Bardsir regions only two crops were harvested (first crop: 84 days after planting and the second crop: 135 days after planting).

2.4. Calculation of Drought Tolerance/Susceptibility Indices

I100 (full irrigation) and I50 treatments were considered as water treatments under regularly irrigated and water deficiency conditions. Fresh yield in regularly irrigated conditions was (Yp) and in water stress conditions, it was (Ys). The tolerance/susceptibility indices for each plant genotype were obtained using the following equations [16].
S S I = 1 Y s Y p 1 Y ¯ s Y ¯ p .
where Ys is the wet forage yield under water deficiency and Yp is the wet forage yield under regularly irrigated condition, Y ¯ s and Y ¯ p are the mean yield (t ha−1) of all genotypes under water deficiency and regularly irrigated conditions, respectively.
M P = Y s + Y p 2 ( t   h a 1 )
T O L = Y P Y s ( t   h a 1 )
G M P = Y P × Y s ( t   h a 1 )
S T I = Y p + Y s ( Y ¯ p ) 2
Y I = Y s Y ¯ s
Y S I = Y s Y P
(β)
β is the coefficient linear regression of yield of each genotype in each condition of the environmental index.
The new stress tolerance score (STS) method was used to apply these indices. This method is easier than the factor analysis (FA) method and could be applied to any number of genotypes, which is calculated by the following equation.
S T S = M P + S T I + G M P + Y I + Y S I S S I T O L β
As it is deduced from the tolerance/susceptibility indices equations, high values are related to MP, STI and GMP, YI, and YSI indices and low values are related to SSI, Tol, and β indices, that are indicated higher resistance to drought stress. Therefore, MP, STI, GMP, YI, and YSI have a positive sign and coefficient, while SSI, Tol, and β have a negative coefficient [17]. Raw data cannot be used for calculating STS; hence, these indices must be first standardized and then used; the data are standardized using the following equation.
Z i j = X i j X ¯ i S i ,
where Z i j is the standard score for jth genotype ith index. X i j is the raw data of jth genotype in ith index, and Si is the standard deviation of ith index.

2.5. Statistical Analysis

ANOVA, simple correlation, and mean comparisons with LSD were performed by SAS 9.4 statistical software (SAS Institute Inc. Cary, NC, USA). The FA was used by means of MINITAB 19 software (State College, PA, USA).

3. Results

The results of ANOVA are shown that Yp and Ys were significant in all indices at all three sites in the tested years. Comparison of means indicated that IP13151 and IP13150 genotypes of pearl millet in Kerman and Jiroft, respectively and the speed feed genotype of sorghum in the Bardsir area had the highest yield under stress and non-stress conditions. According to Table 4, Table 5 and Table 6, there was a significant positive correlation between the yields (Yp) and (Ys), so that in Kerman (r = 0.76**, r = 0.98**, and r = 0.99**) and in Jiroft (r = 0.98**, r = 0.98**, and r = 0.99**).
This correlation implies that by using normal conditions it could perform a direct selection for stress conditions, and this role indicates the yield potential under drought stress conditions [18]. However, Si-o-seh Mardeh et al. (2006) [19] reported a negative correlation between yields in stress and non-stress conditions, which was inconsistent with the results of the present study. Eight drought tolerance/susceptibility indices (SSI, MP, TOl, GMP, STI, YI, YSI, and β) were used for selecting the favorable genotypes with acceptable yield in different climatic conditions.
According to the Table 4, Table 5 and Table 6, four indices MP, GMP, STI, and YI had positive and strong correlations with Yp and Ys. This was also reported by Ilker et al. (2011) [20] that high-yielding genotypes could be assessed by MP, GMP, and STI indices under stress and non-stress conditions; while Jafari (2009) [21] introduced the STI index for this selection. According to the results of the present study, MP, GMP, STI, and YI were the most appropriate and effective indices. On the other hand, three indices of SSI, TOL, and β correlated negatively with yield under stress conditions.
According to Table 7, the first two factors were calculated using the FA method with varimax rotation in 2016, 2017, and 2018, and there was a high factor loading among MP, TOL, GMP, and YI in Kerman, Bardsir, and Jiroft regarding the first factor (FA1) for every three years. These indices had positive coefficients for FA1. There was a high factor loading between the SSI and YSI indices for the second factor (FA2) with a positive coefficient. However, these two indices had negative coefficients in Jiroft and Bardsir in the second year. However, FA1 and FA2 were considered as drought tolerant factors. Higher FA1 + FA2 score shows higher drought tolerance and lower genotypic susceptibility scores indicate the susceptibility to drought.

4. Discussion

The strong and positive correlation among MP, GMP, STI, and YI indices were shown that they could produce identical results. Selections based on these indices resulted in genotypes with high yield.
In contrast, there was a strong and negative correlation between β and TOL indices with Ys (Table 4, Table 5 and Table 6). Hence, according to Fernandez (1992) [10], low-yield potential genotypes can be selected based on TOL under non-stress conditions and high-yield potential genotypes can be selected under stress conditions. Thus, the idea of selection based on STI, MP, GMP, and YI was employed to select favorable low-yield potential genotypes.
In this study, the IP13151 and TP13150 pearl millet genotypes and speed feed sorghum genotype (with the highest yield under stress conditions) were introduced as the most drought tolerant genotypes in Kerman, Jiroft, and Bardsir in the tested years, respectively (Table 8). Comparison of this result with the result of multivariate statistical methods that can be used the indices simultaneously, such as FA and STS methods are given in this table.
Based on the first two factors (FA1 + FA2), for example, in Kerman, in the three consecutive crop years, the ICSV-112 genotype with a score of 7.71 in the first year, and a score of 4.20 in the second year was the most resistant genotype, respectively. However, in the third year, the highest FA score was for S35 genotype of pearl millet with the value of 0.86, which was contradictory to previous result, which introduced the IP13151 genotype as the highest potential yield. In Jiroft, genotypes IP13150 (15.96), HHVBC (1.64), and ICSV-112 (1.07) were also, respectively, introduced based on FA rankings. This result was contradictory to the previous result. HYB-2 (0.90), speed feed (1.84), and S35 (0.78) genotypes were also introduced as the drought tolerant genotypes in Bardsir for three consecutive years, respectively.
The other statistical multivariate method that can use all the indices simultaneously, is the STS method; it is very easy to use and is calculated according to Equation (9). It is observed from the STS result (Table 8), in Kerman, that the STS scores were 4.60, 8.72, and 6.53 for the IP13151 genotype in the three consecutive crop years, respectively, that was introduced similar to results previously obtained based on the indices. Furthermore, in Jiroft, the obtained STS results of IP13150 genotype, with the sores 8.86, 3.36, and 4.90, respectively, were presented as the genotype with an acceptable yield in these climatic conditions under drought stress conditions. In Bardsir, the highest STS scores were related to the speed feed genotype of sorghum that were 6.95, 4.07, and 5.74, respectively. STS biplot graphs indicating the resistant genotypes in each region and year are presented in Figure 2.

5. Conclusions

The multivariate statistics such as FA and STS that simultaneously applies all drought tolerance/susceptibility indices can be helpful for selecting drought tolerant genotypes. The multivariate STS method is a simpler, more understandable, and more efficient method than the FA method. Contrary to the FA method, STS could be employed to select and introduce the most tolerant and susceptible genotype even with a small number of genotypes. FA is a complex method and it usually leads to errors in selecting genotypes with a small number of genotypes. In this study, it was tried to select the most tolerant genotype to drought stress using STS in different climatic conditions and the results were compared with the results of the FA.
Finally, based on the new stress tolerance score (STS) method, the best genotype was identified, hence this method can be recommended to use by other researchers.

Author Contributions

Collected the data and wrote the first draft of manuscript, M.N.; Planned, supervised the project and reviewed the final manuscript, E.T.-N. and G.K.-N.; Provided the plant materials, selected the experimental locations and genetic consultant, B.N.; Conduct data analysis, G.M.-N.

Acknowledgments

We appreciate the support of Research and Technology Institute of Plant Production (RTIPP), Shahid Bahonar University of Kerman, Kerman, Iran. Special thanks to Sonia Aghighi for editing the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Blum, A. Drought resistance–is it really a complex trait? Funct. Plant Biol. 2011, 38, 753–757. [Google Scholar] [CrossRef]
  2. Yadav, O.; Bhatnagar, S. Evaluation of indices for identification of pearl millet cultivars adapted to stress and non-stress conditions. Field Crop. Res. 2001, 70, 201–208. [Google Scholar] [CrossRef]
  3. Zegada-Lizarazu, W.; Iijima, M. Deep root water uptake ability and water use efficiency of pearl millet in comparison to other millet species. Plant Prod. Sci. 2005, 8, 454–460. [Google Scholar] [CrossRef]
  4. Diouf, O.; Brou, Y. Response of pearl millet to nitrogen as affected by water deficit. Agronomie 2004, 24, 77–84. [Google Scholar] [CrossRef]
  5. Gholamhoseini, M.; Aghaalikhani, M.; Malakouti, M. Effect of natural zeolite and nitrogen rates on canola forage quality and quantity. J. Sci. Technol. Agric. Nat. Resour. 2008, 12, 537–548. [Google Scholar]
  6. Bidinger, F.; Mahalakshmi, V.; Rao, G.D.P. Assessment of drought resistance in pearl millet (Pennisetum americanum (L.) Leeke). II. Estimation of genotype response to stress. Aust. J. Agric. Res. 1987, 38, 49–59. [Google Scholar] [CrossRef]
  7. Mahalakshmi, V.; Bidinger, F. Water deficit during panicle development in pearl millet: Yield compensation by tillers. J. Agric. Sci. 1986, 106, 113–119. [Google Scholar] [CrossRef]
  8. Bruck, H.; Payne, V.; Sattelmacher, B. Effects of phosphorus and water supply on yield, transpirational water-use efficiency, and carbon isotope discrimination of pearl millet. Crop Sci. 2000, 40, 120–125. [Google Scholar] [CrossRef]
  9. Rosielle, A.; Hamblin, J. Theoretical aspects of selection for yield in stress and non-stress environment 1. Crop Sci. 1981, 21, 943–946. [Google Scholar] [CrossRef]
  10. Fernandez, G.C. Effective selection criteria for assessing plant stress tolerance. In Proceedings of the International Symposium on Adaptation of Vegetables and other Food Crops in Temperature and Water Stress, Shanhua, Taiwan, 13–16 August 1992. [Google Scholar]
  11. Fischer, R.; Maurer, R. Drought resistance in spring wheat cultivars. I. Grain yield responses. Aust. J. Agric. Res. 1978, 29, 897–912. [Google Scholar] [CrossRef]
  12. Bouslama, M.; Schapaugh, W.T. Evaluation of three screening techniques for heat and drought tolerance. Crop Sci. 1984, 24, 933–937. [Google Scholar] [CrossRef]
  13. Gavuzzi, P.; Rizza, F.; Palumbo, M.; Campanile, R.G.; Ricciardi, G.L.; Borghi, B. Evaluation of field and laboratory predictors of drought and heat tolerance in winter cereals. Can. J. Plant Sci. 1997, 77, 523–531. [Google Scholar] [CrossRef]
  14. Ouk, M.; Basnayake, J.; Tsubo, M.; Fukai, S.; Fischer, K.S.; Cooper, M.; Nesbitt, H. Use of drought response index for identification of drought tolerant genotypes in rainfed lowland rice. Field Crop. Res. 2006, 99, 48–58. [Google Scholar] [CrossRef]
  15. Hao, Z.F.; Li, X.H.; Su, Z.J.; Xie, C.X.; Li, M.S.; Liang, X.L.; Weng, J.F.; Zhang, D.G.; Li, L.; Zhang, S.H. A proposed selection criterion for drought resistance across multiple environments in maize. Breed. Sci. 2011, 61, 101–108. [Google Scholar] [CrossRef] [Green Version]
  16. Gholamin, R.; Zaeifizadeh, M.; Khayatnzhad, M. Factoranalysis for performance and other characteristics in durum wheat under drought stress and without stress. Middle East J. Sci. Res. 2010, 6, 599–603. [Google Scholar]
  17. Abdolshahi, R.; Safarian, A.; Nazari, M.; Pourseyedi, S.; Mohamadi-Nejad, G. Screening drought-tolerant genotypes in bread wheat (Triticum aestivum L.) using different multivariate methods. Arch. Agron. Soil Sci. 2013, 59, 685–704. [Google Scholar] [CrossRef]
  18. Dadbakhsh, A.; Yazdan-Sepas, A. Evaluation of drought tolerance indices for screening bread wheat genotypes in end-season drought stress conditions. Adv. Environ. Biol. 2011, 5, 1040–1045. [Google Scholar]
  19. Mardeh, A.S.S.; Ahmadi, A.; Poustini, K.; Mohammadi, V. Evaluation of drought resistance indices under various environmental conditions. Field Crop. Res. 2006, 98, 222–229. [Google Scholar] [CrossRef]
  20. Lker, E.; Tatar, Ö.; Tonk, F.A.; Tosun, M. Determination of tolerance level of some wheat genotypes to post-anthesis drought. Turk. J. Field Crop. 2011, 16, 59–63. [Google Scholar]
  21. Jafari, A.; Paknejad, F.; Jami, A.-A. Evaluation of selection indices for drought tolerance of corn (Zea mays L.) hybrids. Int. J. Plant Prod. 2012, 3, 33–38. [Google Scholar]
Figure 1. The temperature In Kerman, Jiroft, and Bardsir in three years.
Figure 1. The temperature In Kerman, Jiroft, and Bardsir in three years.
Agronomy 09 00645 g001
Figure 2. Bi-plot based on the stress tolerance score (STS) in genotypes pearl millet and sorghum in three years. (a) 2016, (b) 2017, and (c) 2018.
Figure 2. Bi-plot based on the stress tolerance score (STS) in genotypes pearl millet and sorghum in three years. (a) 2016, (b) 2017, and (c) 2018.
Agronomy 09 00645 g002
Table 1. Geographical characteristics of the study sites in different regions.
Table 1. Geographical characteristics of the study sites in different regions.
KermanJiroftBardsir
Above sea level1756 m680 m2047 m
Longitude30°29′28°40′29°55′
Latitude57°6′57°44′56°34′
Average rainfall in 10 year122.7 mm180 mm92.5 mm
Soil typeL-SLL-S
Growth season20 June–31 October20 March–2 August30 April–11 September
Table 2. Soil properties of the study sites in Kerman, Jiroft, and Bardsir during the growing season.
Table 2. Soil properties of the study sites in Kerman, Jiroft, and Bardsir during the growing season.
N%P
(mg kg−1)
K
(mg kg−1)
pHEC
(dsm−1)
OC%Fe
(mg kg−1)
Zn
(mg kg−1)
Kerman0.494707.54.21.281
Jiroft0.5114407.54.61.16.21.2
Bardsir0.38.84607.54.11.46.80.8
Table 3. The number of cut and irrigation and amount of water consumed per plant and irrigation treatment in Kerman, Jiroft, and Bardsir.
Table 3. The number of cut and irrigation and amount of water consumed per plant and irrigation treatment in Kerman, Jiroft, and Bardsir.
PlantCutNumber of IrrigationTotal Amount of Applied Water (m3 ha−1)
Kerman
I100
Pearl millet215.286744.51
Sorghum215.287495.91
I50
Pearl millet215.283372.26
Sorghum215.283747.95
Jiroft
I100
Pearl millet327.757726.8
Sorghum327.758731.28
I50
Pearl millet327.753863.4
Sorghum327.754365.64
Bardsir
I100
Pearl millet215.286455.44
Sorghum215.287423.75
I50
Pearl millet215.283227.72
Sorghum215.283711.72
Table 4. Simple correlation coefficients between Yp, Ys, and drought tolerance/susceptibility indices of genotypes in Kerman in three years.
Table 4. Simple correlation coefficients between Yp, Ys, and drought tolerance/susceptibility indices of genotypes in Kerman in three years.
2016YpYsSSIMPTOLSTIGMPYIYSIβ
Yp1
Ys0.76 **1
SSI0.71 *0.081
MP0.98 **0.84 **0.591
TOL0.97 **−0.58 **0.85 **0.92 **1
STI0.96 **0.89 **0.510.99 **0.88 **1
GMP0.96 **0.90 **0.500.99 **0.87 **0.99 **1
YI0.76 **0.99 **0.090.85 **0.59 *0.89 **0.90 **1
YSI−0.68−0.05−0.99 **−0.57−0.83 **−0.48−0.47−0.061
β−0.69 **0.08−0.97 **−0.58 **−0.82 **−0.50 *0.480.09−0.97 *1
2017YpYsSSIMPTOLSTIGMPYIYSIβ
Yp1
Ys0.98 **1
SSI−0.11−0.261
MP0.99 **0.99 **−0.84 **1
TOL0.99 **0.39−0.82 **0.99 **1
STI0.99 **0.99 **−0.85 **0.99 **0.99 **1
GMP0.99 **0.99 **−0.85 **0.99 **0.99 **0.99 **1
YI0.98 **0.99 **−0.86 **0.99 **0.99 **0.99 **0.99 **1
YSI0.23 **0.38 **−0.78 **0.80 **0.78 **0.81 **0.80 **0.81 **1
β0.12−0.030.48−0.06−0.01−0.09−0.07−0.1−0.411
2018YpYsSSIMPTOLSTIGMPYIYSIβ
Yp1
Ys0.99 **1
SSI−0.22 **−0.29 **1
MP0.99 **0.99 **−0.23 **1
TOL0.99 **-0.43 *−0.180.99 **1
STI0.99 **0.99 **−0.250.99 **0.99 **1
GMP0.99 **0.99 **−0.250.99 **0.99 **0.99 **1
YI0.99 **0.99 **−0.290.99 **0.99 **0.99 **0.99 **1
YSI−0.08 **−0.01 **−0.76 **−0.060.110.06−0.04−0.0091
β0.160.08−0.87 **0.140.200.120.120.08−0.8 **1
Significant at the * 0.05 and ** 0.01 probability level.
Table 5. Simple correlation coefficients between Yp, Ys, and drought tolerance/susceptibility indices of genotypes in Jiroft in three years.
Table 5. Simple correlation coefficients between Yp, Ys, and drought tolerance/susceptibility indices of genotypes in Jiroft in three years.
2016YpYsSSIMPTOLSTIGMPYIYSIβ
Yp1
Ys0.98 **1
SSI0.008−0.181
MP0.99 **0.98 **−0.031
TOL0.99 **0.350.100.99 **1
STI0.99 **0.99 **−0.080.99 **0.97 **1
GMP0.99 **0.99 **−0.080.99 **0.98 **0.99 **1
YI0.98 **0.99 **−0.170.98 **0.95 **0.99 **0.99 **1
YSI−0.090.09−0.97 **−0.04−0.180.0060.0030.091
β0.96 **−0.91 **0.190.96 **0.98 **0.93 **0.94 **0.91 **−0.271
2017YpYsSSIMPTOLSTIGMPYIYSIβ
Yp1
Ys0.98 **1
SSI0.02−0.151
MP0.99 **0.98 **−0.011
TOL0.99 **−0.96 **0.110.99 **1
STI0.99 **0.99 **−0.070.99 **0.97 **1
GMP0.99 **0.99 **−0.060.99 **0.98 **0.99 **1
YI0.98 **0.99 **−0.150.98 **0.96 **0.99 **0.99 **1
YSI−0.20−0.02−0.94 **−0.16−0.26−0.11−0.11−0.031
β0.96 **−0.92 **0.200.96 **0.98 **0.93 **0.95 **0.92 **−0.371
2018YpYsSSIMPTOLSTIGMPYIYSIβ
Yp 1
Ys0.99 **1
SSI−0.13−0.601
MP0.99 **0.99 **−0.131
TOL0.99 **−0.76 **−0.110.99 **1
STI0.99 **0.99 **−0.110.99 **0.99 **1
GMP0.99 **0.99 **−0.140.99 **0.99 **0.99 **1
YI0.99 **0.99 **−0.160.99 **0.99 **0.99 **0.99 **1
YSI−0.24−0.19−0.50 **−0.22−0.26−0.19−0.21−0.191
β0.98 **−0.97 **−0.150.98 **0.98 **0.95 **0.97 **0.97 **−0.391
Significant at the ** 0.01 probability level.
Table 6. Simple correlation coefficients between Yp, Ys, and drought tolerance/susceptibility indices of genotypes in Bardsir in three years.
Table 6. Simple correlation coefficients between Yp, Ys, and drought tolerance/susceptibility indices of genotypes in Bardsir in three years.
2016YpYsSSIMPTOLSTIGMPYIYSIβ
Yp1
Ys0.99 **1
SSI0.250.131
MP0.99 **0.99 **0.221
TOL0.99 **−0.41 *0.320.99 **1
STI0.99 **0.99 **0.170.99 **0.98 **1
GMP0.99 **0.99 **0.20.99 **0.99 **0.99 **1
YI0.99 **0.99 **0.130.99 **0.97 **0.99 **0.99 **1
YSI−0.19−0.08−0.86−0.16−0.25−0.1−0.14−0.071
β0.47−0.350.95 **0.440.63 **0.380.410.35−0.861
2017YpYsSSIMPTOLSTIGMPYIYSIβ
Yp1
Ys0.99 **1
SSI−0.52−0.591
MP0.99 **0.99 **−0.541
TOL0.99 **0.38−0.470.99 **1
STI0.99 **0.99 **−0.560.99 **0.99 **1
GMP0.99 **0.99 **−0.560.99 **0.99 **0.99 **1
YI0.85 **0.86 **−0.60.85 **0.84 **0.85 **0.86 **1
YSI0.520.59−10.540.470.560.560.60 *1
β0.98 **−0.96 **−0.410.97 **0.99 **0.97 **0.97 **0.80 **0.411
2018YpYsSSIMPTOLSTIGMPYIYSIβ
Yp1
Ys0.99 **1
SSI−0.46−0.411
MP0.99 **0.99 **−0.451
TOL0.98 **−0.48 *−0.480.99 **1
STI0.99 **0.99 **−0.430.99 **0.99 **1
GMP0.99 **0.99 **−0.430.99 **0.99 **0.99 **1
YI0.99 **0.99 **−0.410.99 **0.98 **0.99 **0.99 **1
YSI−0.19−0.110.23−0.17−0.24−0.16−0.15−0.11
β0.26−0.16−0.570.240.320.220.210.16−0.81
Significant at the * 0.05 and ** 0.01 probability level.
Table 7. Factor analysis (FA) for drought tolerance/susceptibility indices of genotypes in Kerman, Jiroft, and Bardsir in three years.
Table 7. Factor analysis (FA) for drought tolerance/susceptibility indices of genotypes in Kerman, Jiroft, and Bardsir in three years.
Kerman 2016 2017 2018
Factor loading Factor loading Factor loading
IndexFA1FA2Com.FA1FA2Com.FA1FA2Com.
SSI−0.2930.9550.996−0.1640.9860.997−0.0540.9980.99
MP0.960−0.27610.999−0.03810.9990.0121
TOL0.992−0.1230.9960.9910.12110.9920.1170.999
STI0.952−0.3010.9980.992−0.1010.9940.997−0.0260.998
GMP0.946−0.32110.995−0.09310.999−0.0261
YI0.918−0.3940.9970.981−0.19110.994−0.1011
YSI0.2900.9560.9960.1670.9850.9970.0540.9980.997
B−0.2210.9740.970.1520.9870.9950.0950.9940.998
%Var0.5170.4770.9940.710.2870.9980.5740.4240.998
Jiroft 2016 2017 2018
Factor loading Factor loading Factor loading
IndexFA1FA2Com.FA1FA2Com.FA1FA2Com.
SSI−0.207−0.9770.9480.0440.99810.2300.9311
MP0.9180.3950.9980.9990.03910.998−0.0561
TOL0.9840.1740.9920.9850.1700.9990.994−0.0951
STI0.8530.5120.9920.995−0.0190.9910.993−0.0450.989
GMP0.8880.4580.9990.999−0.00710.999−0.0411
YI0.8240.56410.994−0.10110.999−0.0121
YSI0.2070.9770.976−0.046−0.99810.1020.9601
B0.990−0.0400.9790.9550.2570.9790.976−0.1030.973
%Var0.5430.4420.9860.7330.2630.9960.7420.2530.995
Bardsir 2016 2017 2018
Factor loading Factor loading Factor loading
IndexFA1FA2Com.FA1FA2Com.FA1FA2Com.
SSI0.0850.9961−0.275−0.9560.998−0.4180.6581
MP0.9850.16710.9490.3070.9910.992−0.0101
TOL0.9630.26710.9710.2270.9720.980−0.1850.999
STI0.9920.1180.9980.9350.34310.995−0.0610.998
GMP0.9890.14010.9410.33110.995−0.0800.981
YI0.9960.0760.9970.7980.4280.9680.998−0.0291
YSI−0.0820.9960.9960.261−0.9600.9470.1610.9660.997
B0.3120.9490.970.9810.14410.1260.9710.998
%Var0.6210.3780.9940.6690.30810.630.3690.369
Note: Com. = Communality.
Table 8. Mean yield in regularly irrigated (Yp) and water deficiency (Ys) conditions and drought tolerance/susceptibility indices for genotypes in three years.
Table 8. Mean yield in regularly irrigated (Yp) and water deficiency (Ys) conditions and drought tolerance/susceptibility indices for genotypes in three years.
Kerman
2016 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC67.6927.321.2947.5140.370.27430.750.42.23.21−1.92
IP2226978.8634.511.2156.6944.350.452.150.950.442.070.26−0.43
HYB-184.8636.141.2460.548.720.4555.3610.432.14−0.021.59
HYB-290.9237.161.2864.0453.760.558.121.020.412.24−0.372.23
IP1315099.6539.81.369.7359.850.5962.981.10.42.17−0.524.34
ICM-0522282.8133.121.357.9749.690.4152.370.910.42.271.22−1.19
IP1315199.9640.81.2870.3859.160.6163.861.130.412.27−0.894.6
S3571.3436.181.0653.7635.160.3850.7610.511.763.89−3.88
ICS-11268.735.341.0352.0233.360.3649.160.970.531.727.71−6.35
SPEEDFEED92.9237.681.2965.355.240.5259.171.040.412.1−0.273.01
ICS-1400181.7536.951.1859.3544.80.4554.951.020.452.01−0.351.39
PEGAH66.0233.861.0549.9432.160.3347.280.930.511.73−2.430.52
LSD2.0730.093.686.560.053.380.080.050.21
2017 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC85.5427.541.0456.2858.530.2948.070.920.322.920.82−6.31
IP2226993.3430.331.0362.84650.3653.771.030.322.920.37−3.49
HYB-196.4631.561.0256.4464.890.3855.171.070.332.84−1.441.14
HYB-2103.4433.771.0260.2369.680.4359.11.150.332.86−1.282.99
IP13150111.7437.031.0165.174.710.5164.321.260.332.83−2.285.95
ICM-0522297.5230.571.0455.0866.950.3754.61.040.312.971.3−0.62
IP13151115.7138.221.0365.0877.50.5566.51.30.332.85−2.058.72
S3572.4122.581.0448.0149.830.240.420.770.312.911.61−3.17
ICS-11274.6322.321.0646.4352.310.2140.810.760.32.044.2−4.12
SPEEDFEED94.2431.61.0161.7162.640.3754.271.080.342.77−3.065.19
ICS-1400183.6227.841.0154.5855.790.2948.250.950.332.76−2.59−0.15
PEGAH63.3421.38143.7641.970.1736.80.730.342.64−3.50−1.28
LSD6.013.020.093.556.320.053.320.080.050.2
2018 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC84.526.631.0155.5757.870.2747.440.930.323.09−1.15−0.68
IP2226994.2429.84162.0464.40.3453.021.040.323.08−1.511.71
HYB-196.0330.311.0163.1765.720.3553.951.060.323.09−1.261.67
HYB-2102.8432.041.0167.4470.80.457.41.120.313.14−0.351.39
IP13150111.5636.321.0173.9475.250.4963.651.270.333.01−0.644.6
ICM-0522297.1930.291.0163.7466.910.3654.251.060.313.14−0.380.3
IP13151113.7736.02174.977.750.564.011.260.323.1−2.116.53
S3571.9822.041.0247.0149.950.1939.820.770.313.170.86−6.29
ICS-11273.9522.691.0248.3251.260.240.960.790.313.160.66−5.65
SPEEDFEED95.0230.32162.6764.70.3553.671.060.323.06−2.032.78
ICS-1400183.1826.61154.956.580.2747.040.930.323.04−2.200.85
PEGAH63.0720.28141.6842.780.1535.760.710.323−2.62−1.98
LSD6.353.090.093.796.510.063.50.080.050.21
Jiroft
2016 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC85.6925.311.0255.560.380.2646.570.720.31.46−4.99−2.92
IP2226995.6627.741.0361.767.920.3251.510.790.291.7−1.59−2.92
HYB-196.7229.091.0162.9167.630.3453.040.830.31.650.75−1.06
HYB-2103.132.910.9968.0170.190.4158.230.940.321.666.952.73
IP13150109.9838.460.9474.2271.520.5165.011.090.351.615.968.86
ICM-0522297.5728.351.0362.9669.220.3352.590.810.291.73−0.66−2.65
IP13151116.5237.80.9877.1678.720.5366.361.070.321.814.395.45
S3570.8221.21.0246.0149.620.1838.740.60.31.01−10.64−3.17
ICS-11273.8421.411.0347.6352.430.1939.760.610.291.16−10.75−4.39
SPEEDFEED94.3531.020.9762.6963.330.3554.10.880.331.474.513.21
ICS-1400183.6625.471.0154.5758.180.2646.160.720.31.38−4.36−1.66
PEGAH63.519.75141.6343.750.1535.40.560.310.76−11.91−1.82
LSD8.313.260.095.546.060.064.970.090.020.24
2017 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC85.7126.021.0355.8759.690.2747.230.890.31.421.64−4.03
IP2226995.6830.091.0162.8965.590.3553.651.030.311.570.53−0.74
HYB-196.7531.60.9964.1765.150.3755.291.080.331.52−0.671.68
HYB-2103.1433.01168.0770.130.4158.351.120.321.6501.4
IP13150117.1837.29177.2479.90.5366.11.270.321.840.183.36
ICM-0522297.7730.121.0263.9567.650.3554.271.030.311.631.2−1.70
IP13151108.5835.250.9971.9273.340.4661.961.20.321.7−0.473.21
S3572.3522.98147.6749.370.240.750.780.321.030.05−2.65
ICS-11273.9322.851.0248.3951.080.2141.090.780.311.091.17−4.27
SPEEDFEED94.3531.570.9862.9662.410.3654.581.080.331.45−1.482.86
ICS-1400183.6628.180.9855.9255.480.2848.550.960.341.24−1.671.69
PEGAH63.6221.40.9842.5142.220.1636.890.730.340.7−1.63−0.52
LSD8.311.950.085.076.580.043.990.060.020.25
2018 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC85.3226.57155.9558.740.2847.620.760.311.38−0.73−1.15
IP2226995.2830.09162.6965.190.3553.550.860.321.550.561.54
HYB-196.130.26163.1865.840.3553.930.870.311.570.31.37
HYB-2102.932.11167.570.790.457.480.920.311.69−0.551.48
IP13150117.0536.31177.1879.740.5366.081.070.321.84−0.814.9
ICM-0522297.8830.4164.1467.480.3654.550.870.311.62−0.990.26
IP13151113.8135.58174.778.240.4963.631.020.311.83−0.403.43
S3571.1322.5146.8148.630.2400.640.320.970.9−0.98
ICS-11274.0222.961.0148.4951.050.2141.230.660.311.091.07−2.71
SPEEDFEED94.0330.150.9962.0963.880.3453.240.860.321.510.092.87
ICS-1400183.3326.50.9954.9256.830.2746.990.760.321.31−1.180.57
PEGAH63.3720141.6943.370.1535.60.570.320.740.49−2.05
LSD7.572.290.014.895.410.044.120.070.010.22
Bardsir
2016 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC75.6927.151.0151.4248.540.2945.330.890.362.690.26−4.00
IP2226979.2128.54153.8850.670.3247.550.940.362.68−0.07−2.24
HYB-181.7529.56155.6652.180.3449.160.970.362.67−0.33−0.89
HYB-292.4932.891.0162.6959.60.4355.161.080.362.730.9−0.36
IP1315095.6734.371.0165.0261.30.4757.341.130.362.710.142.39
ICM-0522281.429.181.0155.2952.220.3448.730.960.362.70.32−2.44
IP1315196.0234.76165.3961.260.4857.771.140.362.68−0.403.79
S3577.1328.830.9852.9848.30.3247.150.950.372.58−2.713
ICS-11277.8128.15152.9849.650.3146.80.930.362.67−0.37−1.98
SPEEDFEED97.3335.810.9966.5761.510.559.031.180.372.64−1.606.95
ICS-1400182.730.370.9956.5352.330.3650.1110.372.63−1.441.89
PEGAH67.8925.050.9946.4742.850.2441.240.820.372.6−1.78−1.79
LSD4.392.10.023.192.580.043.010.080.010.07
2017 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC75.6827.461.0151.5748.220.345.580.750.361.27−2.25−2.66
IP2226979.0828.791.0153.9450.30.3347.710.560.361.33−1.86−3.35
HYB-181.5729.651.0155.6151.920.3549.180.810.361.37−2.05−1.76
HYB-292.3233.491.0162.9158.830.4455.60.910.361.54−2.36−0.27
IP1315095.4935.51165.559.980.4858.230.970.371.540.62.68
ICM-0522281.2230.01155.6151.210.3549.370.820.371.35−0.17−0.36
IP1315195.9935.76165.8760.230.4958.590.970.371.540.822.95
S3576.9428.54152.7448.410.3146.860.780.371.260.28−0.56
ICS-11277.6228.78153.248.840.3247.270.780.371.280.26−0.54
SPEEDFEED97.436.60.996760.80.5159.7110.381.551.844.07
ICS-1400182.7730.091.0156.4352.680.3649.90.820.361.39−2.02−1.54
PEGAH68.1724.791.0146.4843.380.2441.10.680.361.1−2.01−3.20
LSD4.481.920.023.073.12.072.810.190.010.09
2018 GenotypeYp (t ha−1)Ys (t ha−1)SSIMP (t ha−1)TOL (t ha−1)STIGMP (t ha−1)YIYSIβFA1 + FA2STS
HHVBC75.6726.47151.0749.20.2944.750.90.352.76−1.51−0.98
IP2226979.0727.361.0153.2251.720.3146.510.930.352.79−2.84−0.80
HYB-181.5627.961.0154.7653.610.3347.750.950.342.82−3.97−2.39
HYB-291.9432.44162.1959.50.4354.621.10.352.75−1.184.38
IP1315095.1133.59164.3561.520.4656.521.140.352.75−1.185.51
ICM-0522280.8728.28154.5752.590.3347.820.960.352.76−1.70−1.02
IP1315194.4733.480.9964.8861.830.4657.041.150.352.74−0.785.08
S3576.727.380.9952.0449.310.345.830.930.362.70.780.64
ICS-11278.1127.850.9952.9850.250.3146.640.950.362.710.570.26
SPEEDFEED9834.40.866.263.60.4858.071.170.352.77−1.955.74
ICS-1400183.2529.510.9956.3853.750.3549.5610.352.73−0.361.91
PEGAH68.0524.110.9946.0843.950.2340.50.820.352.710.41−4.01
LSD4.781.910.013.33.070.042.980.060.010.07

Share and Cite

MDPI and ACS Style

Negarestani, M.; Tohidi-Nejad, E.; Khajoei-Nejad, G.; Nakhoda, B.; Mohammadi-Nejad, G. Comparison of Different Multivariate Statistical Methods for Screening the Drought Tolerant Genotypes of Pearl Millet (Pennisetum americanum L.) and Sorghum (Sorghum bicolor L.). Agronomy 2019, 9, 645. https://doi.org/10.3390/agronomy9100645

AMA Style

Negarestani M, Tohidi-Nejad E, Khajoei-Nejad G, Nakhoda B, Mohammadi-Nejad G. Comparison of Different Multivariate Statistical Methods for Screening the Drought Tolerant Genotypes of Pearl Millet (Pennisetum americanum L.) and Sorghum (Sorghum bicolor L.). Agronomy. 2019; 9(10):645. https://doi.org/10.3390/agronomy9100645

Chicago/Turabian Style

Negarestani, Mohammadreza, Enayatollah Tohidi-Nejad, Gholamreze Khajoei-Nejad, Babak Nakhoda, and Ghasem Mohammadi-Nejad. 2019. "Comparison of Different Multivariate Statistical Methods for Screening the Drought Tolerant Genotypes of Pearl Millet (Pennisetum americanum L.) and Sorghum (Sorghum bicolor L.)" Agronomy 9, no. 10: 645. https://doi.org/10.3390/agronomy9100645

APA Style

Negarestani, M., Tohidi-Nejad, E., Khajoei-Nejad, G., Nakhoda, B., & Mohammadi-Nejad, G. (2019). Comparison of Different Multivariate Statistical Methods for Screening the Drought Tolerant Genotypes of Pearl Millet (Pennisetum americanum L.) and Sorghum (Sorghum bicolor L.). Agronomy, 9(10), 645. https://doi.org/10.3390/agronomy9100645

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop