The present study aimed to investigate the genetic diversity of foxtail millet germplasm and select its superior genotypes. To this end, 130 foxtail millets were studied under an augmented complete block design at Shahid Bahonar University Research Farm. Analysis of Variance (ANOVA) was also carried out on the genotypes to determine the homogeneity of experimental field. The results proved that there were significant differences among the control genotypes regarding all their traits except for the plant loading trait. Based on the results of stepwise regression for grain yield as a dependent variable in foxtail millet genotypes, the number of grains per panicle, the number of panicles, 1000-seed weight, and height respectively had the best efficiency as index traits in selecting genotypes. Based on the cluster analysis, the foxtail millet genotypes were classified into 5 groups. Comparison of the genotypes demonstrated that the genotypes S27 and S7 and cultivar Boston respectively, had the highest biological yield and Bastan and genotype S27 had the highest grain yield among others. The current study revealed that genotypes S7 and S27 had a higher yield than Iran’s only breed cultivar, namely Bastan and therefore they can be used in breeding studies in order to produce superior cultivars.
Keywords: Biological yield, Cluster analysis, Foxtail millet, Grain yield, Stepwise regression
Foxtail millet is widely planted in northern China and India. The amount of nutrients this plant contains is as much as or even more than that of wheat or rice. Every single grain has in it 12.3% of protein, 4.7% of fat and 60.6% of carbohydrate. The protein in this grain is as much as that of wheat, but more than that of rice. The amount of minerals, vitamins and amino acid it contains is very high. It also has a high rate of dietary fibers and non-starchy polysaccharide. Furthermore, it has a low glycemic index and releases sugars slowly [1-3]. Foxtail millet is a diploid (2n=2x=18) and self-pollinating plant belonging to Setaria italica, which is mainly planted in India, Russia, China and Middle Eastern countries . Foxtail millet has so far been under-researched but identifying its genotypes, particularly local genotypes and analyzing its genetic diversity and relations have a key role in improving its yield and in protecting and using this plant . Regarding the low amount of genetic reserves and considering the fact that many useful genes have been lost and also exposition of crops to unhealthy environmental conditions and biotic and abiotic stress, has made studying and identifying genetic diversity in plants a priority which is regarded as the first and foremost step to identifying and maintaining genetic resources which in turn are the basis for carrying out genetic researches and plant breeding programs. Breeding and producing new cultivars are dependent upon the exact selection among plants which is itself dependent on identifying cultivars and their diversity [5,6]. The most important criterion in identifying superior genotypes is studying morphological markers and yield components. One of the ways to identify the morphological traits relating to yield is identifying the correlation between the traits and the yield [7,8]. Research findings prove that there is positive correlation between grain yield and number of panicles . In their study on 20 foxtail millet genotypes, Channappagoudar et al.  demonstrated that there is a significantly positive correlation between grain yield and biomass and plant height, photosynthesis rate and number of tillers.
The current research aimed to investigate the genetic diversity in 130 foxtail millet genotypes. It also tried to analyze the relations between morphological traits and grain yield. Clustering the genotypes according to different traits using cluster analysis and selecting superior genotypes was among the final purposes of the present study.
MATERIALS AND METHODOLOGY
In the present study, 130 Foxtail millet genotypes (Table 1) were collected from different parts of Iran. They were then planted under an augmented complete block design at Shahid Bahonar University of Kerman Research Farm at a proper time and in the same manner. Each one of the genotypes was planted in a 2 m long line. The distance between the lines was 40 cm and each plant was at a distance of 8 cm to the next. The processes of planting, growing and harvesting were the common ways as regards those of the foxtail millet. Ten plants were randomly chosen from each plot to measure morphological traits, namely the plant height (cm), panicle length (cm), number of panicles per plant, number of grains per panicle, dry weight per plant (g), grain weight per plant (g), 1000-seed weight, length and width of the leaf (cm), number of tillers per plant, number of leaves, plant loading, bristle color, grain color, grain yield (kg/ha), biological yield (kg/ha) and harvest index. Due to the qualitative nature of the traits like bristle color and grain color, the data were codified in the following way:
Green bristle (1), purple bristle (2), yellow grain (1), brown grain (2), red grain (3).
To perform variance analysis for the measured traits, based on the type of the augment, randomized complete block design (CRD) was used for control cultivars and the experimental error in this analysis was applied as standard error in different comparisons. To compare the means, LSD was calculated at 1% level and the statistical analysis of the data was carried out according to descriptive statistics like minimum, maximum, mean, standard deviation, diversity coefficient of the traits. Phenotype correlation of the traits and their significance was also measured. To identify the influence of each trait on the yield, linear stepwise regression was used for the grain yield as the dependent variable and other traits as independent variables. Cluster analysis was also carried out applying a Euclidean distance matrix through UPGMA. Statistical software like SAS 9.2, Microsoft Excel 2017, XlStat 2017 were used in this study.
RESULTS AND DISCUSSION
In the present research, analysis of the control genotypes was performed according to the randomized complete block design and the results are as presented in Tables 2 and 3. Regarding the implementation of the augment design and the results of the variance analysis, the amount of each of the traits in the genotypes under study was corrected based on the control genotypes. The results showed a significant difference among the control genotypes regarding all the traits except for plant loading. Among the different traits, plant height and leaf length were significant at 5% level. The maximum height was related to genotype S69 (130.60 cm) and the minimum height belonged to genotype S83 (31.60 cm). The longest leaf belonged to genotype S72 (42.30 cm) and shortest was related to genotype S41 (12.30 cm). The traits number of tillers, number of leaves, width of the leaf, length of the panicle, number of the panicles, number of grains per panicle, bristle color, grain color, dry weight per plant, grain weight per plant, 1000-seed weight, harvest index, grain yield and biological yield were also at 1% level of significance. The findings of the research showed that genotypes S27, S7 and Bastan cultivar had the highest biological yield, respectively. Bastan cultivar and genotype S27 also had the highest grain yield among other genotypes.
Based on the corrected traits, the descriptive statistics including the maximum, minimum, mean, standard deviation, diversity coefficient of the traits was calculated. As presented in Table 3, the traits number of grains per panicle, biological yield, grain weight per plant, Dry weight per plant and grain yield had the highest diversity coefficient and the 1000-seed weight the lowest. On the whole, most of the traits had high phenotype diversity. The range of variation revealed a wide spectrum for most of the traits which indicated a high diversity among the genotypes under study. As an instance, grain yield range of variation was between 340.63 and 2931.04 and 1769.04 to 29408.45 kg/ha for biological yield. The high range of variation for most of the traits under study demonstrates a significant variation among Iranian millets. Therefore, selecting a special genotype in Foxtail millet germplasm in order to reach different breeding goals, including improving grain yield and biological yield can be of many benefits . In a similar way observed a significant deviation among these traits in their study on different millets. In order to better investigate the genotypes regarding the traits under question, a comparison of means was made between the genotypes using an LSD method as shown in Table 4.
Phenotype correlation was measured through Pearson correlation coefficient as demonstrated in Table 5. The results show that the correlation of grain yield with all traits except for 1000-seed weight and number of tillers was positively significant. Among these traits, number of grains per panicle and biological yield had the highest correlation (80%) with grain yield. Channappagoudar et al.  also reported a correlation between grain yield and height and biological yield. Santhakumar  observed a positive correlation between height, panicle length, biological yield and grain yield. Grain yield and panicle height and length were reported to have had a significant correlation .
The present research showed a significant relation between most of the traits except for the 1000-seed weight and biological yield. The correlation is negative (-0.23) with harvest index, but positive with other traits. The harvest index had a significantly positive correlation with grain yield, but due to the inverse relationship between biological yields, it had a negative correlation. The harvest index showed the highest negative correlation (-0.34) with the number of leaves. There is also a significantly positive relationship between plant height and leaf length, leaf width, panicle length, number of grains per panicle, grain yield and biological yield. The research also proved a significantly positive correlation between leaf length and height, leaf width, panicle length, number of grains per panicle, grain yield and biological yield. Among the traits under study, the highest correlation was related to the number of tillers and number of panicles (0.86), grain yield with number of grains per panicle (0.80) and grain yield with biological yield (0.80). Based on these findings, it can be concluded that the number of grains per panicle and biological yield are the most important traits in improving grain yield in Foxtail millet. Azari Nasrabad and Mirzaee  reported similar results and findings in their studies.
Table 6 represents the results of the stepwise regression for grain yield (kg/ha) as a dependent variable and other measured traits as independent variables. Biological yield, harvest index, grain weight per plant and dry weight per plant were controlled and not considered as dependent variables due to their close relationship with grain yield. In the regression model and considering the regression coefficient obtained, the traits with higher coefficients have a more important role in improving the yield. Among the measured traits, four traits got into the regression model including the number of grains per panicle, number of panicles, 1000-seed weight and height. The first trait in the model was the number of grains per panicle which had 64.5 of grain yield variation. The R2 related to the regression model shows that almost 81% of grain yield variance change is contingent upon the traits in the model. Dadashi et al.  reported that the number of grains per panicle, the number of fertile tillers, and 1000-seed weight got into the model and according to R2, these traits were reported to have justified 96% of the changes in the grain yield.
The final model obtained for grain yield is as follows:
Y, X1, X2, X3 and X4 are grain yield (kg/ha), number of grains per panicle, number of panicles, 1000-seed weight, and height respectively, which are regarded as the main parameters in grain yield. Therefore, in order to select genotypes with high ratios of yield panicles, number of grains per, these traits ought to be taken into consideration. Generally in cereals, the number of fertile panicle, and the average grain weight are regarded as the major components of the yield .
In order to cluster the genotypes under question, cluster analysis was used through UPGMA applying a Euclidean distance matrix. The tree diagram obtained from cluster analysis (Figure 1) shows five genotype groups containing 65, 55, 4, 5 and 1 genotypes. The genotypes in group 1 were almost average in most of their traits. Comparison of the characteristics of the genotypes of these groups showed that the mean of the most of the traits related to group 2 genotypes was lower than the total mean and the mean of the traits related to group 3 was higher than that of the total. The genotypes of the group 3 have a great potential regarding their grain yield, biological yield and height. The Bastan cultivar with high grain and biological yields also belong this group. This is while the genotypes belonging to this group are somehow week regarding these traits. Therefore, if programs are about to improve grain and biological yields, selecting from group 3 may be a proper choice. Group 4 genotypes are very close to those of group 3 considering most of their traits and had the highest rates regarding 1000-seed weight and harvest index traits. Group 5 contained only 1 genotype, which had a high rate regarding most of the traits. It was reported that cluster analysis using Ward method put the different barley genotypes in 4 groups at normal conditions and 3 at drought stress.
The results of the cluster analysis showed that there is a high genetic diversity in the germplasm under study. Based on the clustering and the estimated means of the traits for the genotypes in each cluster, proper parents can be selected in the improvement programs to carry out purposeful crossings. Since there is closer affinity between the genotypes of these clusters compared to other clusters, crossings between clusters which have more differences can create a higher genetic diversity and heterosis.
To put everything in a nutshell, the present study showed that there was a sufficient diversity in the germplasm considering all the traits investigated. According to the results of the stepwise regression model for grain yield as a dependent variable in Foxtail millet genotypes, the number of grains per panicle, number of panicles, 1000-seed weight, height and number of tillers as indexing traits in selecting genotypes had the highest efficiency in selecting genotypes. Also, the present study revealed genotypes (S27 and S7) that had a higher rate of yield compared to the only Iranian cultivar, namely Bastan. These genotypes can be used in breeding and improving studies to produce superior cultivars.
1. Kumari R, Dikshit N, Sharma D, Bhat KV (2011) Analysis of molecular genetic diversity in a representative collection of foxtail millet (Setaria italica L.) from different agro-ecological regions of India. Physiol Mol Biol Plants 17: 363-374.
2. Li P, Brutnell TP (2011) Setaria viridis and Setaria italica, model genetic systems for the Panicoid grasses. J Exp Botany 62: 3031-3037.
3. Sapkota M, Pandey M, Thapa D (2016) Agromorphological characterisation of foxtail millet (Setaria italica L.) at Rampur, Chitwan, Nepal. Int J Appl Sci Biotechnol 4: 298-307.
4. Siyadat A, Modhaj A, Esfahani M (2013) Cereals. University of Ferdoosi Mashhad Publications: Mashhad, Iran, pp: 317-332.
5. Aghaee-Sarbarzeh M, Amini A (2012) Genetic variability for agronomy traits in bread wheat genotype collection of Iran. Seed Plant Improv J 1: 581-599.
6. Mohammadi SA, Prasanna BM (2003) Analysis of genetic diversity in crop plants, salient statistical tools and considerations. Crop Sci 43: 1235-1248.
7. Aydin N, Sermet C, Zeki Mut H, Bayramoglu O, Ozcan H (2010). Path analyses of yield and some agronomic and quality traits of bread wheat (Triticum aestivum L.) under different environments. Afr J Biotechnol 9: 5131-5134.
8. Sabaghpour SH, Safikhani M, Sarker A, Ghaffari A, Ketata H (2004). Present status and future prospects of lentil cultivation in Iran. 5th European Conference on Grain Legumes, AEP Publ: Dijon, France, p: 23.
9. Shooshi Dezfooli A, Mehrani A (2011) Investigation of correlation between yield and component yield in foxtail millet species. J Agronomy Sci 41: 413-421.
10. Channappagoudar BB, Hiremath SM, Biradar NR, Koti RV, Bharamagoudar TD (2008) Physiological basis of yield variation in foxtail millet. Karnataka J Agric Sci 20: 481-486.
11. Santhakumar G (1999) Correlation and path analysis in foxtail millet. J Maharashtra Agric Univ 24.
12. Siles MM, Ken Russell W, Baltensperger DD, Nelson LA, Johnson B, et al. (2004) Heterosis for grain yield and other agronomic traits in foxtail millet. Crop Sci 44: 1960-1965.
13. Azari Nasrabad M, Mirzaee R (2012). Effect of sowing date on grain yield and yield components of foxtail millet (Setaria italica) promising lines. Seed Plant Prod J 28: 95-105.
14. Dadashi MR, Nourinia A, Asgar M, Azizi SH (2010) Correlation assessment of a number of physiological and morphological characteristics of barley hybrids with grain yield. J Ecophysiol Crops Weeds 15: 29-40.
15. Koocheki A, Banayan-Aval A, Rezvani P, Mahdavi-Damghani A, Jamiolahmadi M, et al. (2005) The plant ecophysiology. University of Ferdoosi, Mashhad Publications: Mashhad, Iran, p: 271.
- Proteomics and Bioinformatics (ISSN:2641-7561)
- Journal of Genomic Medicine and Pharmacogenomics (ISSN:2474-4670)
- Food and Nutrition-Current Research (ISSN:2638-1095)
- Journal of Microbiology and Microbial Infections (ISSN: 2689-7660)
- Journal of Veterinary and Marine Sciences (ISSN: 2689-7830)
- Journal of Agriculture and Forest Meteorology Research (ISSN:2642-0449)
- Journal of Biochemistry and Molecular Medicine (ISSN:2641-6948)