Genetic Diversity in Local and Exotic Avena sativa L. (Oat) Germplasm Using Multivariate Analysis

: Avena sativa L., also known as Oat belongs to the Poaceae family, is one of the most signiﬁcant crops that is grown for its seeds, fodder as well as for human consumption as oatmeal. In the current study, 236 genotypes of A. sativa were analysed for genetic diversity through agro-morphological and SDS-PAGE analysis. Cluster analysis based on agro-morphological characteristics grouped all the genotypes into nine clusters, whereas genotype numbers 537 and 728 were highly different from others. The seed yield production of cluster 9 genotypes was the highest per plant (38.2 ± 0.20 g), while cluster 2 genotypes produced maximum biomass per plant (122.5 ± 9.55 g) as compared to other clusters. In a principal component analysis where four variables were studied, and the observed total variations were 57.60%. Among the genotypes, a maximum grain yield of 38.2 g (each) was recorded for genotypes 22,350 and 728, followed by genotypes 737 and 22,390 (with 36.4 g and 35.6 g of seed productions, respectively). The SDS-PAGE analysis resulted in 13 bands and all the genotypes were grouped into seventeen clusters. At the extreme periphery of the dendrogram, genotype 537 and 22,332 were considered to be the most diverse genotypes. Our ﬁndings have implications for both understanding the diversity and relationships among these diverse genotypes of A. sativa and will provide a basis for obtaining the elite germplasm optimally adapted to local conditions. The selected genotypes based on agronomic performance may be potential breeding material to raise successful future cultivars.


Introduction
Avena sativa L. (oat) belongs to the family Poaceae is a highly nutritive plant cultivated globally [1]. Oat is ranked sixth in the world after wheat, rice, maize, barley and sorghum as a cereal producing plant [2]. Oat is a major component of infant foods due to its high nutritional profile, lack of allergenicity, palatable flavour, good shelf-life, stability and low cost. Oat based food includes oat bran, oatmeal, oat flour and oat flakes which are mainly used as breakfast cereals. Porridge, hot cereals, bread, biscuits, infant food, muesli and granola bars are a few examples of food products produced from oats [3]. Oat is also a decent source of lipids and unsaturated fatty acids. The high lipid content in the seed makes it a potential oilseed crop [4]. It is a major winter forage crop and is cultivated total seed storage proteins using SDS-PAGE and (iii) to find out promising genotypes that can fulfil fodder demands of farmers in terms of yield-related traits.

Experimentation
The present research work was conducted in the Botanical Garden and Herbarium, (34 • 40 318 N and 72 • 0 753 E with 726 m altitude using Altimeter eTrex®30x (Garmin, Schaffhausen, New Taipei City, 221, Taiwan) in the University of Malakand, Pakistan (with an average of annual temperature of 28 • C and annul precipitation of 80 mm). The data was recorded for two consecutive years (November 2018 to May 2020). The seeds of 236 genotypes were acquired from the Gene-bank of Plant Genetic Resources Institute, Islamabad, Pakistan. Out of which, 143 genotypes were indigenous (Pakistani origin) and 92 were exotic (USA origin) ( Table S1). The most prominent variety locally grown (Dir Lower, Khyber Pakhtunkhwa) was used as a control variety.

Germination of Plant Material
The experiments were performed in triplicate applying randomised complete block design. About 6 g of dried seeds were sown in a given row. The dimensions of plots where the seeds were grown have 10 m 2 area with 2.5 m length and row to row distance of 0.3 m. The experiment was performed in a loamy soil: 3.60% clay, 54% silt, sand 44.40%, nitrogen 0.04%, phosphorus 39.5 mg/kg, potassium 128 mg/kg and pH 8. The recommended practices were applied for soil preparation, ploughing, levelling, irrigation (after 2 days of sowing) and weeding. The herbicides (Dual-gold 960 EC, Syngenta, Basel, Switzerland) was applied before sowing to avoid unwanted plants growth.

Morphological Evaluation
A total of 21 morphological traits (7 qualitative and 14 quantitative) were monitored following guidelines of IPBGR-1985 devised for Oat (Avena sativa L.). Out of these the seven qualitative traits were observable-stem thickness, plant vigour, plant lodging, leaf rigidity, rigidity of flag leaf, erectness of spikelet and glume colour of spikelet which were recorded in percentages while the quantitative traits were days to flowering, days to maturity, number of tiller per plant, number of leaves per plant, plant height, number of spikelet, peduncle length, spike length, leaf length, leaf width, biomass, yield per plant, harvest index and 1000 seeds weight which were recorded for two consecutive years specified above. Days to flowering were measured from the sowing date till the first flower initiate. Days to maturity were noted when the colour of the plant became yellow. The total number of tillers per plant and leaves were counted for 5 plants. The plant height was the length recorded from base to top of the plant with meter rod after maturation. The spikelet was enumerated for five plants in each cultivar, whereas spike length, peduncle length, leaf length and leaf width were measured with a metre rod. The weight of resultant biomass was recorded where the whole plant was cutting into pieces and digital balance was used to record its mass. Seed yield per plant was determined after harvesting and threshing of all the spikelet. The average seed weight was calculated by taking 1000 seeds of each genotype and weighing it's through digital balance. Harvest index was calculated as the ratio between the grain yield and biological yield.

Total Seed Protein Analysis
The seeds proteins of all 236 genotypes were extracted and subjected to SDS-PAGE analysis to determine genetic variation among these germplasms. Using mortar and pestle, a total of 10 healthy and mature seeds were ground into fine powder. About 0.02 g of the seeds powder were transferred to 1.5 mL centrifuge tubes already containing 400 µL protein extraction buffer (Distilled water 70 mL, Tris 0.60 g, SDS 0.2 g, Urea 30.3 g, 2-Mercaptoethanol 1 mL and pH 8). The powder was homogenised at room temperature for 1 min before centrifugation at 12,000 rpm for 10 min. About 14% polyacrylamide gel was used to make SDS-PAGE slabs, consisting of two phase's separation gel and stacking gel. The composition of separation gel was 20% acrylamide with 0.135 percent N.N-methylene-acrylamide in 0.15 M Tris-HCl buffer with a pH of 8.8 and SDS 0.27%. For polymerization of the gel, 10% APS (Ammonium per sulphate) and fifteen microliters Tetramethylethylene-diamine (TEMED) was added to the solution. For making of stacking gel 0.8% N.N-methylene-bis-acrylamide, 30% acrylamide 0.25 Tris (pH of 6.8), 0.2% SDS, 70 µL APS and 17 µL TEMED were mixed together. About 10 µL of the sample solution was loaded into the well of the stacking gel through a micropipette. The gel plates were mounted on gel assembly and a 100 voltage current was applied to get the electrophoretogram. When the blue line of bromophenol blue reached the bottom of the gel the electric current supply was disconnected. The gel plates were taken out from the assembly and transfer to a plate filled with staining solution and gently rotated on a shaker for 2 h. After that for de-staining, gels were kept in the destining buffer for several hours to overnight till the disappearance of the blue colour for the gel and protein subunits were clearly observed [20].

Data Analysis
For the estimation of multivariate diversity and selection of candidate lines for high yield and fodder quality, 5 plants were randomly selected from each line (genotype) and statistically analysed. The distribution of qualitative traits among the germplasms, data frequency in percent was calculated using Microsoft Excel (Redmond, WA, USA). The quantitative data obtained for cluster analysis, principal component analysis and Pearson correlation analysis were computed using Statistica version 7 (Palo Alto, CA, USA), PC ORD version 5 (Wagga Wagga, Australia) and NTSys PC version 2.1 (Bangalore, India). Similarly, the binary data matrices of seed storage protein profile were analysed via Unweighted Pair Group Arithmetic Average (UPGMA), genetic similarity analysis and genetic linkage analysis, where 0 has been used for the absence of a given band while 1 for the presence of a band. A two-way cluster dendrogram has been constructed using PC-ORD software version 5 (Wagga Wagga, Australia).

Diversity in Qualitative Traits
Significant variations were observed for stem thickness: thin (10.25%), intermediate (40.25%) and thick (49.15%). Three types of plant growth vigour: good (58.90%), medium (29.66%) and poor (11.44%) were recorded as well. Lodging is one of the most important traits: 4.24% of genotypes had lodging while 95.76% of genotypes had no lodging. Three types of rigidity of flag leaf were noted: bent 35.59%, slightly bent 34.74% and stiff 29.66%. Three types of rigidity among leaves of the germplasms were noted: in 9.32% leaves were bent, in 48.31% slightly bent while in 46.61% germplasm they were stiff. Three types of glume colour were observed: green 85.32%, silver 13.98% and light green 1.69%. Erectness of spikelet observed were 3 types: droop 16.10%, semi-erect 35.59% and erect 48.31%.

Diversity in Quantitative Traits
The descriptive statistics (range, mean standard deviation and CV %) of different traits recorded in the current study are presented in Table 1. Overall, significant variations were observed in all the studied traits and the coefficient of variation (CV) are ranged from 2.56 to 30.73% with the highest CV recorded for harvest index with an increasing trend from 1st to 2nd year (CV = 29.50 and 30.73% correspondingly). An increase in CV% for 1st to 2nd year in the number of leaves plant −1 from 19.95 to 25.29% and plant height (17.25 to 18.26%) were observed. A similar trend in leaf length and leaf width from 1st to 2nd year was reordered (leaf length means correspondingly for 1st and 2nd year = 33.102 ± 6.10 and 33.63 ± 7.18 while for leaf width these are 1.72 ± 0.27 and 1.74 ± 0.31). A decreasing trend in CVs for the number of tiller/plant (24.13% and 20.18%) and plant biomass (23.43 and 20.17%) were observed correspondingly from the 1st to 2nd year. The CVs were recorded for days to maturity (2.76 and 2.56%), days to flowering (3.06% and 3.02%) with decreasing trends from 1st to 2nd year. The early maturity of 172 days was observed for indigenous accessions; 537 and 538 while for exotic varieties it was observed for 727, 753 and 758 accessions. For first-year plants, all remaining accessions matured till 196 days whereas in the 2nd year the maturity period was from 174 to 190 days. A significant variation in yield per plant −1 were observed, ranged from 8.80 to 38.20 g with maximum yield pant −1 38.20 g recorded for the genotypes 22,350 (Pakistani origin) and 728 (USA origin) whereas minimum yield recorded was 8.80 g for the genotype 545 (Pakistani origin).

Correlation Analysis
The correlation analysis were performed considering the 14 quantitative traits that were monitored for 2 consecutive years ( Table 2). DM was found to strongly correlate with DF, in particular in 2020 (r = 0.86 **). In 2019, NT/P showed a negative correlation with DF and DM with the value r = 0.23 * for each, while no significant differences in the correlations were registered for the year 2020. There is a significant correlation between PH and NT/P in 2019 (r = 0.50 **), which is reduced (r = 0.25 *) in the 2nd year (2020). LW and Y/P have a low correlation coefficient to NT/P for 2019 (r = 0.16 * and r = 0.21 **) and 2020 (r = 0.29 * and r = 0.20 *), respectively. Plant biomass and LL also showed the lowest correlation coefficient (r = 0.11 * and r = 0.39 **) with LL. Similarly, for the harvest index, the data were strongly significant with Y/P (0.57 **), while there was a negative correlation with BM (0.58 **), which for the 2nd, it was strongly positively correlated with Y/P (0.62 **), and a negative correlation with PH (−0.40 **) and BM (−0.64 **). The parameter 1000 seed weight was found to be strongly significant with Y/P (0.53 **) and HI (0.28 **), while there was a negative correlation with LW (−0.23 **), which for the 2nd year was found to be significant with HI (0.16 *) and negative with NSPK (−0.20 *), respectively.

Principal Component Analysis (PCA)
The principal component analysis (PCA) helps select the best cultivars among a given population. Table 3 shows the PCA of the 14 studied quantitative traits. Altogether, PC1, PC2, PC3 and PC4 explained 57.6% (2019) and 56.4 % (2020) of the total variation in the data set.
In the 1st PC component, the total variations were 20.66% and they were associated positively with DF and DM, negatively with NT/P, PH, NSPKL and Y/P; whereas for the 2nd year the total variation was 19.83%, with a positive value of NT/P, PH, LW and BM while negative in HI, NSPKL, SPKL, LL, Y/P and HI. The 2nd principal component accounted for 33.79% of the total variation and was positively related to DF, DM, NSPKL, SPKL, LL, LW and BM, whereas, in the 2nd year, the total variation was 35.225% with a positive association with DF, DM and NSPKL, and negatively with LW, Y/P, HI and 1000 seed weight. Furthermore, DM, Y/P, HI and 1000 seed weight contributed positively to PC3, which accounted for 46.35% of the overall variations, whereas PH and LL were found negative. While for the 2nd year plants the total variation was 47.77%. In addition, the overall variation in PC4 was 57.59% for the 1st year and 56.43% for 2nd year the association of SPK L, LW and HI were positive while 1000 seed weight was negative, whereas LL, BM and 1000 seed weight were associated positive while NT/P, SPK L, LW, NSPKL and HI were recorded negative. Among the PCS (PC1, PC2, PC3 and PC4) DM was ranked high followed by Y/P, HI and LW. Biplot analysis was used to determine the associations between the traits. As a result, the plot defined by the two PCS, PC2 and PC3 showed a total variation of 80.15% for 2019 and were placed against each other to find the level of genetic diversity and geographical association among the genotypes. The separation based on PC2 and PC3 revealed that the genotypes were scattered in all the quarters, which showed the high level of genetic diversity in the studied genotypes. Among the studied genotypes of both local (Pakistan) and exotic (USA) was found diversely scattered and most of the genotypes was found condensed toward the centre of the plot. According to the biplot, the genotype-728, 22 (Figure 1). For the 2nd year genotypes, the total variation of PC1 and 2 was 55.05% and were plotted against each other to get the scatter diagram. Genetic diversity based on origin-wise shows that the genotypes of Pakistan and the USA were found scattered.

Cluster Analysis
Wards method was followed while performing the cluster analysis. All the genotypes were divided into three main lineages (lineage 1, lineage 2 and lineage 3) at a distance of 25%, which were further subdivided into 9 sub-clusters at a distance of 62.5% for the year 2019. Lineage 1 was further divided into 5 sub-clusters whereas lineage 2 was divided into 1 cluster Lineage 3 into 3 sub-clusters, with genotypes 537 and 728 as the most diverse among the studied genotypes. In a given group there was less variation among the genotypes, but they were considerably different from members of the other groups ( Figure 3). Cluster analysis of the 2nd year (2020) plants demonstrated 100 genotypes which were placed into two main Lineages (lineage1 and lineage 2) at a distance of 25%, which were further subdivided into 5 sub-clusters. Lineage 1 was further divided into 2 sub-clusters while lineage 2 into 3 sub-clusters. The number of genotypes in both lineage 1 and 2 are presented in Figure 4. The passport information of 236 genotypes along with their accessions detail and origin (from Pakistan and the USA) used in the current study has been shown in Figure S1 Supplementary Materials). The mean performance of all traits were calculated for both the years where cluster 6 was found to have the lowest DF and DM value (2019) while in (2020) while the lowest DM were recorded in cluster 4. The high grain yield genotypes were from cluster 1 (2019) with a mean value of 38.2 ± 0.2 g, followed by cluster 3 genotypes with a mean value of 22.16 ± 5.71, whereas in 2020, the highest grain yield was found in cluster 2, which may be selected in future breeding practices of the oat plant. In cluster 2, the highest BM mean value = 122.5 ± 9.55 was recorded (Tables 4 and 5).

Total Seed Storage Protein
All the genotypes were investigated for total seed storage proteins using SDS-PAGE, where 13 polymorphic bands were observed with the molecular weight ranging from 25 to 180 kDA ( Figure 5). The highest genetic diversity was observed in Band 7 (19.39%) with the frequency distribution of 128, followed by Band 10 (17.95%) with a frequency distribution of 136 while the minimum genetic variance was recorded in Band 1 (0.72%) with a frequency distribution of 232. The genetic diversity and frequency distribution of 13 Bands are presented in Table S2 Supplementary Materials).
Based on the cluster dendrogram all the 236 genotypes were grouped into 3 main lineages (lineage 1, lineage 2 and lineage 3) at a distance of 25% which were further subdivided into subgroups at a distance of 50%. The 1st lineage is divided into 3 subclusters, 2nd into 1 sub-cluster while the 3rd into 3 sub-clusters. The maximum number of genotypes 33.23% were found in cluster 1 followed by cluster 5 with 19.49% genotypes while a minimum of 4.23% genotypes were observed in cluster 2. The 1st sub-cluster of lineage 1 was consist of 80 genotypes in which 7 genotypes were of USA origin and 73 from Pakistan origin. The 2nd sub-cluster of the same lineage was consist of 10 genotypes of USA origin. The 3rd sub-cluster of 1st lineage is composed of 34 genotypes of USA origin. The 4th cluster of lineage 2 consists of a total of 22 genotypes from which 8 genotypes were of Pakistan origin and 14 from the USA. Sub-cluster 5 comprises 46 genotypes of which 10 were of USA origin and 36 from Pakistan. Sub-cluster 6 was consist of 34 genotypes of which 10 genotypes were of USA origin and 24 from Pakistan. Similarly, cluster 7 was comprised of 10 genotypes of which 2 genotypes were from the USA and 8 from Pakistan origin ( Figure 6).

Discussion
For improvement in any breeding program, there is always need to evaluate the material available to find out the good lines for breeding purposes. Different agronomic traits have been used as a prime objective to evolve new breeding materials. Determination of correlation helps in taking the decision related to selection strategies. In breeding programs, the studies on diversity are the key factors to bring genetic improvement in a given population. This path analysis further helps in determining the direct and indirect effect on the traits of interest. All these traits are affected by genetic and environmental factors [21]. Genetic variation within and between population crops species is of major interest for plant breeders. An estimate of the extent of variation within and between populations of a species is useful in analysing the genetic structure of germplasm. Ahmad et al. [22] studied clustering patterns in 75 germplasms of oat to get high yield varieties with desired traits.
The extent of genetic diversity in genetic materials paves the way for any program of genetic improvement and genetic heterogeneity is measured in different ways, and the most commonly used approach is multivariate analysis [23]. Multivariate analyses like morphological and biochemical differences along with the response of a plant towards molecular markers are employed to determine genetic differences among the population of a given plant. Based on morphological differences, one can pick and preserve genotypes for future use [18]. Unfortunately, morphological parameters are strongly influenced by the climate, which may positively or negatively affect a given plant and data for every locality around the globe is necessary for successful breeding of the given plant in a given locality [24].
For most of the traits studied significant differences between the genotypes were observed such as days to flowering, days to maturity, plant height and the number of tiller plant 1 , yield plant 1 , total biomass and 1000 seed weight. Early flowering initiation was observed in genotypes 543, 657, 661, 663 and 736 with days to flowering less than 150 days which is a good sign for breeders. The other important factor in breeding practices is the time taken to maturity which for 4 genotypes was found to be less than 175 days. Similar results have been reported previously by other investigators as well [25][26][27][28] where they observed significant variation for days to flowering and days to maturity in various oat germplasm. Green fodder yield is highly dependent on the number of tillers plant-1, which is certainly the most important factor. The maximum number of tillers per plant 14.2 was found for accession 727 while the minimum of 3.40 was recorded for accession 602. Our result is in close agreement with that of Wagh et al., [25] and Naeem et al., [27] where the authors found significant variations for the number of tillers per plant. A maximum plant height of 145.20 cm was recorded for the accession 804 while the minimum was 59.60 for 22,355. Our results are in close agreement with that of Bibi et al. [29] where the researchers found significant differences for days to 50% flowering, days to 50% maturity and plant height. The number of leaves plays an important role in plant development and growth which directly affect the yield both for cereal and fodder plants. The number of leaves per plant has a direct effect on the yield of fodder crops. Accession 22,389 with a maximum number of leaves 10 was the most successful genotype while the minimum number of leaves 4 were recorded for accession 687 and 746. The same character with the same trend although for different accession have been studied by other investigators as well [30,31]. Ahmad et al. [32] suggest in their study that diversity in leaf and other traits in various germplasm at different areas may be also recognised to have variation in genetic makeup, soil nutrients and environmental conditions. Higher fodder yields in oat genotypes can be due to their bigger leaf area, which results in higher photosynthetic activities and a high ability to retain assimilative photosynthesis products [33]. The most important characteristic of a seed-producing plant is its seed yield per plant as the seed is the edible part with high nutritive values. Significant variations were found with respect to seed yield in the genotypes studied whereas the mean value calculated was 20.78 ± 5.17. Maximum Y/P 38.20 g were recorded for the accession 22,350 (Pakistan) and 728, for 723 (USA)-22,324-22,347 the grain yield were (34 g), while for accession 22,390 (Pakistan) the seed yield was (35.6 g) and for 737 (USA) the yield of (36.4g) was recorded. Similar results have been obtained by other investigators as well [27,34]. The knowledge of the relationship between yield and its component characteristics is paramount to make the best use of these relationships in the selection of promising plant types.
In the studied qualitative traits, a number of useful variations were observed. Among the genotypes, with regard to stem thickness, the plant was categorised into three groups: thin (10.25%), intermediate (40.25%) and thick (49.15%) stem plants. Previously, significant variation in stem thickness and other important agro morphological traits have been reported by Kaur et al. and are in close agreement with our study [35]. Chakraborty [36] have evaluated 180 genotypes of Oat for different morphological traits and have found significant variation in growth habit, foliage, plant stature panicle shape and awn. Sumathi and Balamurugan [37] have evaluated 11 Oat varieties and have found significant variation in various traits.
Correlation analysis is another criterion for selecting a plant for successful breeding practices. The findings of this current study are comparable with the reported studies [30,[38][39][40] where the contributors found a significant correlation for plant height, 1000 seed weight and filling period. Ahmad et al. [32] reported in their study that plant height were found a positive association with culm diameter, leaf stem ratio, number of tillers and 1000 seed weight and seed length. It can be assumed that agro-morphological screening is extremely necessary and that the selection of possible future lines must take into account both positive and negative associations of useful traits to prevent yield depression or any other quality trait [41].
Moreover, the clustering revealed that genetic diversity and the geographical area were not clustered together since the accessions collected from the same geographical region were divided into separate clusters and those collected from different geographical regions tended to be clustered together and were directly related. The studied genotypes of A. sativa were placed into three main lineages (lineage 1, lineage 2 and lineage 3) and 9 sub-clusters after applying cluster analysis. In a given group variations among the plants were minimum but considerably high when compared to other group plants. In the hybridization program, these varieties would be used to achieve desirable characteristic's plants at least in the locality where the study has been performed. Our results are comparable with the reported study in literature [35] where the authors have found significant divergence in 96 oat genotypes and have divided them into 6 clusters. The maximum number of germplasm in a single group shows that these genotypes are more nearly related to one another and having less divergence among them. It can be inferred from the study that the heterogeneity observed within the 236 genotypes could not only be due to their geographical origin or environment but also due to their genetic makeup as predicted by cluster and biochemical analysis. Such clustering helps in selecting desirable individuals with specific traits for the crossing program. Days to maturity, plant height, number of spikelets, spike length, leaf length and width, biomass, yield/plant, harvest index and seed weight has contributed more towards divergence, while peduncle length, number of leaves and teller per plant shows less contribution toward divergence, so the direct selection from these traits would be helpful. Similar results have been reported by Ahmed et al. [22].
The biochemical analysis using SDS-PAGE has a great contribution to the estimation of genetic diversity. For the past few decades, biochemical analysis has been getting popular, while determining the genetic diversities among the given plant germplasms and are more beneficial as compared to classical morphological analysis [42]. For protein profiling, SDS-PAGE is one of the most helpful techniques, which is widely used to study genetic differences among different crops [43]. The seed storage protein profiling were made and all the genotypes were grouped into three main lineages which were then further divided into 7 sub-clusters. Among the studied genotypes, the 1st group consist of 4 accessions, the 2nd consist of 16 genotypes while the 3rd is comprised of 34 genotypes. Similarly, a greater degree of polymorphism has been observed by Mirza et al., [44] in Avena fatua where the authors found different polymorphic bands. Similarly, a barely greater degree of polymorphism has been recorded [45] where the contributor found a total of 12 bands.

Conclusions
The diversity reported in a considerable number of A. sativa germplasm originating from diverse agro-ecologies are exploited by the plant breeders for genetic improvement of this plant. To get an improved variety of given germplasm the differences in morphological and biochemical parameters are investigated which are the basic criterion for successful breeders. In this study, 236 genotypes were studied for differences in morphological and biochemical traits. Based on seed yield as a nutritive parameter, accession 22,350, 728, 737, 22,390, 22,347, 22,324, 723 and 22,395 was the best genotypes among the studied genotypes. The selected accessions can be used as breeding material to improve the nutritional value of high-yielding oat cultivars. Accession 22,348, 722, 686, 22,335, 843, 664 and 788 as high yield fodder in terms of biomass was the most successful accession. The studied parameters in terms of correlation showed significantly positive associations with each other and thus have the potential to be used in future for the breeding program while selecting the high candidates' lines. Based on cluster analysis accessions 537 early maturing and 728 based on high yielding were the most appropriate genotypes to be grown on a large scale which would increase per acre production of the selected plant.