# Phenotypic Diversity Assessment of Okra (Abelmoschus Esculentus (L.) Moench) Genotypes in Ethiopia Using Multivariate Analysis

Authors: Jemal Mohammed; Wassu Mohammed; Eleni Shiferaw
Journal: Scientifica
 (2022)

---

## Abstract

Okra is a minor crop that has not gained research attention in Ethiopia. Characterization of such underutilized crops has important implications for their utilization. Thus, this study was conducted to assess the genetic diversity of okra genotypes in Ethiopia using agromorphological and biochemical markers. Thirty-six okra genotypes were evaluated for 29 agromorphological and biochemical traits. The results of the analysis of variance showed significant differences among genotypes for most of the traits, except for the number of flower epicalyx and fruit diameter. Results of the principal component analysis indicated that the first eight principal component axes accounted for 3.83 to 30.54% and 82.44% of the total variability. Genetic distances estimated by Euclidean distances from 27 traits ranged from 3.55 to 14.49. The 36 genotypes were grouped into four distinct clusters from the Euclidean distance matrix using the unweighted pair group method with arithmetic mean (UPGMA). The first cluster contained 24 (66.66%) genotypes, and the second cluster contained 10 (27.77%) of the genotypes. This study showed the presence of considerable genetic variation among the genotypes for most of the traits, including fruit yield, seed yield, and nutrient content of seeds, indicating the possibility of using these genotypes to develop okra varieties with high fruit-yielding and good nutritional content.


---

## Body

## 1. Introduction

Today, the world’s food supply is based on a small number of crop species, mostly major cereals (wheat, rice, and maize) leaving an abundance of genetic resources and potentially beneficial traits neglected [1]. In the face of climate change, utilizing the vast pool of minor and underutilized crop species would provide a more varied agricultural system and food sources, ensuring food and nutrition security problems. Underutilized crops play a significant role in ensuring food security, nutrition, and income generation for resource-poor farmers and consumers, especially in the developing world [2].Climate change may increase the relevance of plant species that were previously underutilized or thought to be of minor importance [1]. One approach to maintaining a good fit between crops and environmental challenges because of climate change is to use underutilized (minor, orphan, or neglected) crops and their wild relatives. Okra is among the most underutilized crops cultivated in the southwestern and western parts of Ethiopia [3–5]. The characterization of okra genotypes existing in the country would contribute to developing varieties that could thrive in extreme climatic conditions and would allow further utilization of the crop for enhancing food security.The value of germplasm collections depends on their diversity, and crop improvement prominently relies on existing genetic variation [6, 7]. Shujaat et al. [8] suggested that genetic variation is an important feature to achieve the diversified goals of plant breeding, including higher and quality yield, resistance to diseases, and wider adaptations.The pattern and level of genetic diversity in a given gene pool can be measured in terms of genetic distance, which is a measurement of average genetic divergence between genotypes or populations [9]. Regardless of the dataset (morphological, biochemical, or molecular marker data), multivariate analytical procedures that simultaneously make several measurements on each individual under examination are frequently utilized in genetic diversity studies [10]. This study aimed to assess the phenotypic and biochemical diversity of Ethiopian landrace okra genotypes along with exotic commercial varieties using multivariate analysis for further utilization of the crop and contribute to ensuring food security and alleviating malnutrition.


## 2. Materials and Methods

### 2.1. Description of the Study Site

The field experiment was conducted at the Melkassa Agricultural Research Center, Ethiopia, during the 2018 main cropping season (rainy season). Melkassa is located at 8°24′59.20″ N latitude and 39°19′15.19″ E longitude, with an altitude of 1,548 m above sea level [11]. The biochemical contents were determined at the Ethiopian Biodiversity Institute Nutrition Laboratory (total ash and crude fat), the Debre Ziet Agricultural Research Center (crude fiber), and the Melkassa Agricultural Research Center (total protein).


### 2.2. Experimental Materials and Design

Thirty-six okra genotypes of 24 landrace accessions (collected by the Ethiopian Biodiversity Institute from different okra growing regions of Ethiopia), three genotypes (from the Humera Agricultural Research Center), and nine exotic commercial varieties (eight from India and one from the USA) were used in this study. The 36 genotypes were planted in a 6 × 6 simple lattice design. Three seeds per hill were sown and thinned to one plant per hill when plants reached the 3–4 leaf stage.


### 2.3. Data Collection

Data were collected for phenology traits (days to 50% emergence, days to first flowering, days to 50% flowering, and days to 90% maturity), growth and yield-related traits (plant height, stem diameter, number of primary branches per stem, number of internodes, internodes length, leaf length, leaf width, number flower of epicalyxes, peduncle length, fruit length, fruit diameter, average fruit weight, number of tender fruits per plant, number of mature pod per plant, number of ridges on fruit, fruit yield per plant, fruit yield per hectare, number of seeds per pod, hundred seed weight, seed yield per plant, and seed yield per hectare), and biochemical content of the seed (total ash, total fat, crude fiber, and total protein). Phenology and growth-related traits’ data were recorded according to the IPGRI [12] descriptor list developed for okra.


#### 2.3.1. Total Ash

Total ash was determined following the method of AOAC [13] using the gravimetric method. Crucible was cleaned, dried, and ignited at 550°C for 1 hour and weighed (m1). The flour sample (3 g) weighed (m2) and dried at 120°C for 1 hour. Then, the dried sample was carbonized over a blue flame and ignited in a muffle furnace at 550°C until ashing was complete (over 12 hrs). After being ignited, the sample was cooled to ambient temperature and was weighed (m3). Finally, the total ash content was calculated as follows:(1)Ash%=m3−m1m2−m1∗100,where m1 is the mass of crucible (g), m2 is sample mass with crucible (g), and m3 is the final mass of sample with crucible (g).


#### 2.3.2. Crude Fat

The crude fat content of okra seed was determined by the Soxhlet extraction method according to the AOAC [13]. The flour sample (3 g) was weighed and added into a thimble. The thimble with the sample was placed in a 50 ml beaker and dried in an oven for 2 hours at 110°C. A 150–250 ml dried beaker was weighed and rinsed several times with petroleum ether. The sample contained in the thimble was extracted with petroleum ether in a Soxhlet extraction apparatus for 6–8 hours. After extraction is completed, the extracted fat was transferred into a preweighed beaker (Mi). The beaker with the extracted fat was placed in a fume hood to evaporate the solvent on a steam bath unit no odor of the solvent is detectable. Then, the beaker with contents was removed, cooled in a desiccator, and weighed (Mf). The amount of fat in flour was calculated by using the following formula:(2)Fat%=Mf−MiM∗100,where Mf is the dried mass of the fat with beaker (g), Mi is the mass of beaker (g), and M is the sample mass (g).


#### 2.3.3. Crude Fiber

The crude fiber was determined according to the AOAC [13]. Ground sample (3 g) was weighed (m1), placed in 500 ml beaker, digested with 1.25% sulfuric acid, and washed with water and was further digested with 1.25% sodium hydroxide and filtered in course porous (75 μm) crucible in apparatus at a vacuum of about 25 mm. The residual left after refluxing was washed again with 1.25% sulfuric acid at near boiling point. Then, the residual was dried at 110°C overnight, cooled in a desiccator, and weighed (m2). After being dried, the sample was ashed at 550°C until the ashing was complete, cooled in a desiccator, and weighed again (m3). The total crude fiber was expressed in percentages as follows:(3)Fiber%=m2−m3m1∗100,where m1 is a mass of sample (g), m2 is mass of sample with crucible before ashing (g), and m3 is mass of sample with crucible after ashing (g).


#### 2.3.4. Total Protein

A dried and grounded sample was taken (0.5 g) and added into aKjeldahl digestion flask. One gram of catalyst (Na2SO4 mixed with anhydrous CuSO4 in a ratio of 10 : 1) and 5 ml of concentrated H2SO4 were added into the digestion flask. Then, using a digester, the mixed sample was digested at 350°C for about two hours until the sample was completely digested. Then, the flask was removed from the digester and allowed to cool and the digested sample was diluted by adding 30 ml of distilled water. Then, 25 ml concentrated 40% NaOH was added into the digestion flask to neutralize the acid and make the solution slightly alkaline. The contents were immediately distilled by inserting the digestion tube line into the receiver flask that contained 25 ml of 4% boric acid solution and about 150 ml of distillate collected. Then, the distillate was titrated by a standard acid (0.1 N HCl). The percentage of crude protein was calculated by multiplying the nitrogen percentage by the conversion factor (6.25) [13].(4)Nitrogen%=VHCIA−VHCIB∗NHCI∗0.014Wsample∗100,where V = volume of standard acid used for titration of sample (A) and blank sample (B), N = normality of standard acid used for titration (0.1 N HCl), 0.014 is the molecular weight of nitrogen, and W = weight sample taken for digestion, on a dry basis.(5)Protein%=6.25∗%Nitrogen.


### 2.4. Analysis of Variance

The quantitative field data were subjected to analysis of variance (ANOVA) and computed with R statistical software agricolae package [14]. The biochemical traits were analyzed following the CRD (completely randomized design) procedure. The traits that exhibited significant mean squares in ANOVA were further subjected to multivariate analysis.


### 2.5. Principal Component Analysis

Principal component analysis (PCA) was computed to find out the traits, which accounted more for the total variation. The data were standardized to mean zero and variance of one before computing principal component analysis to avoid differences in measurement scales. The principal component based on the correlation matrix was calculated using the R statistical software FactoMineR package [15].


### 2.6. Euclidean Distance and Clustering of Genotypes

Euclidean distance (ED) was computed from quantitative after subtracting the mean value and dividing it by the standard deviation as established by Sneath and Sokal [16]. R statistical software factoextra package [17] was used for the analysis of distance matrix and constructing dendrogram. The dendrogram was constructed based on the unweighted pair group method with arithmetic mean (UPGMA) from the distance matrix of phenotypic traits.


## 2.1. Description of the Study Site

The field experiment was conducted at the Melkassa Agricultural Research Center, Ethiopia, during the 2018 main cropping season (rainy season). Melkassa is located at 8°24′59.20″ N latitude and 39°19′15.19″ E longitude, with an altitude of 1,548 m above sea level [11]. The biochemical contents were determined at the Ethiopian Biodiversity Institute Nutrition Laboratory (total ash and crude fat), the Debre Ziet Agricultural Research Center (crude fiber), and the Melkassa Agricultural Research Center (total protein).


## 2.2. Experimental Materials and Design

Thirty-six okra genotypes of 24 landrace accessions (collected by the Ethiopian Biodiversity Institute from different okra growing regions of Ethiopia), three genotypes (from the Humera Agricultural Research Center), and nine exotic commercial varieties (eight from India and one from the USA) were used in this study. The 36 genotypes were planted in a 6 × 6 simple lattice design. Three seeds per hill were sown and thinned to one plant per hill when plants reached the 3–4 leaf stage.


## 2.3. Data Collection

Data were collected for phenology traits (days to 50% emergence, days to first flowering, days to 50% flowering, and days to 90% maturity), growth and yield-related traits (plant height, stem diameter, number of primary branches per stem, number of internodes, internodes length, leaf length, leaf width, number flower of epicalyxes, peduncle length, fruit length, fruit diameter, average fruit weight, number of tender fruits per plant, number of mature pod per plant, number of ridges on fruit, fruit yield per plant, fruit yield per hectare, number of seeds per pod, hundred seed weight, seed yield per plant, and seed yield per hectare), and biochemical content of the seed (total ash, total fat, crude fiber, and total protein). Phenology and growth-related traits’ data were recorded according to the IPGRI [12] descriptor list developed for okra.


### 2.3.1. Total Ash

Total ash was determined following the method of AOAC [13] using the gravimetric method. Crucible was cleaned, dried, and ignited at 550°C for 1 hour and weighed (m1). The flour sample (3 g) weighed (m2) and dried at 120°C for 1 hour. Then, the dried sample was carbonized over a blue flame and ignited in a muffle furnace at 550°C until ashing was complete (over 12 hrs). After being ignited, the sample was cooled to ambient temperature and was weighed (m3). Finally, the total ash content was calculated as follows:(1)Ash%=m3−m1m2−m1∗100,where m1 is the mass of crucible (g), m2 is sample mass with crucible (g), and m3 is the final mass of sample with crucible (g).


### 2.3.2. Crude Fat

The crude fat content of okra seed was determined by the Soxhlet extraction method according to the AOAC [13]. The flour sample (3 g) was weighed and added into a thimble. The thimble with the sample was placed in a 50 ml beaker and dried in an oven for 2 hours at 110°C. A 150–250 ml dried beaker was weighed and rinsed several times with petroleum ether. The sample contained in the thimble was extracted with petroleum ether in a Soxhlet extraction apparatus for 6–8 hours. After extraction is completed, the extracted fat was transferred into a preweighed beaker (Mi). The beaker with the extracted fat was placed in a fume hood to evaporate the solvent on a steam bath unit no odor of the solvent is detectable. Then, the beaker with contents was removed, cooled in a desiccator, and weighed (Mf). The amount of fat in flour was calculated by using the following formula:(2)Fat%=Mf−MiM∗100,where Mf is the dried mass of the fat with beaker (g), Mi is the mass of beaker (g), and M is the sample mass (g).


### 2.3.3. Crude Fiber

The crude fiber was determined according to the AOAC [13]. Ground sample (3 g) was weighed (m1), placed in 500 ml beaker, digested with 1.25% sulfuric acid, and washed with water and was further digested with 1.25% sodium hydroxide and filtered in course porous (75 μm) crucible in apparatus at a vacuum of about 25 mm. The residual left after refluxing was washed again with 1.25% sulfuric acid at near boiling point. Then, the residual was dried at 110°C overnight, cooled in a desiccator, and weighed (m2). After being dried, the sample was ashed at 550°C until the ashing was complete, cooled in a desiccator, and weighed again (m3). The total crude fiber was expressed in percentages as follows:(3)Fiber%=m2−m3m1∗100,where m1 is a mass of sample (g), m2 is mass of sample with crucible before ashing (g), and m3 is mass of sample with crucible after ashing (g).


### 2.3.4. Total Protein

A dried and grounded sample was taken (0.5 g) and added into aKjeldahl digestion flask. One gram of catalyst (Na2SO4 mixed with anhydrous CuSO4 in a ratio of 10 : 1) and 5 ml of concentrated H2SO4 were added into the digestion flask. Then, using a digester, the mixed sample was digested at 350°C for about two hours until the sample was completely digested. Then, the flask was removed from the digester and allowed to cool and the digested sample was diluted by adding 30 ml of distilled water. Then, 25 ml concentrated 40% NaOH was added into the digestion flask to neutralize the acid and make the solution slightly alkaline. The contents were immediately distilled by inserting the digestion tube line into the receiver flask that contained 25 ml of 4% boric acid solution and about 150 ml of distillate collected. Then, the distillate was titrated by a standard acid (0.1 N HCl). The percentage of crude protein was calculated by multiplying the nitrogen percentage by the conversion factor (6.25) [13].(4)Nitrogen%=VHCIA−VHCIB∗NHCI∗0.014Wsample∗100,where V = volume of standard acid used for titration of sample (A) and blank sample (B), N = normality of standard acid used for titration (0.1 N HCl), 0.014 is the molecular weight of nitrogen, and W = weight sample taken for digestion, on a dry basis.(5)Protein%=6.25∗%Nitrogen.


## 2.3.1. Total Ash

Total ash was determined following the method of AOAC [13] using the gravimetric method. Crucible was cleaned, dried, and ignited at 550°C for 1 hour and weighed (m1). The flour sample (3 g) weighed (m2) and dried at 120°C for 1 hour. Then, the dried sample was carbonized over a blue flame and ignited in a muffle furnace at 550°C until ashing was complete (over 12 hrs). After being ignited, the sample was cooled to ambient temperature and was weighed (m3). Finally, the total ash content was calculated as follows:(1)Ash%=m3−m1m2−m1∗100,where m1 is the mass of crucible (g), m2 is sample mass with crucible (g), and m3 is the final mass of sample with crucible (g).


## 2.3.2. Crude Fat

The crude fat content of okra seed was determined by the Soxhlet extraction method according to the AOAC [13]. The flour sample (3 g) was weighed and added into a thimble. The thimble with the sample was placed in a 50 ml beaker and dried in an oven for 2 hours at 110°C. A 150–250 ml dried beaker was weighed and rinsed several times with petroleum ether. The sample contained in the thimble was extracted with petroleum ether in a Soxhlet extraction apparatus for 6–8 hours. After extraction is completed, the extracted fat was transferred into a preweighed beaker (Mi). The beaker with the extracted fat was placed in a fume hood to evaporate the solvent on a steam bath unit no odor of the solvent is detectable. Then, the beaker with contents was removed, cooled in a desiccator, and weighed (Mf). The amount of fat in flour was calculated by using the following formula:(2)Fat%=Mf−MiM∗100,where Mf is the dried mass of the fat with beaker (g), Mi is the mass of beaker (g), and M is the sample mass (g).


## 2.3.3. Crude Fiber

The crude fiber was determined according to the AOAC [13]. Ground sample (3 g) was weighed (m1), placed in 500 ml beaker, digested with 1.25% sulfuric acid, and washed with water and was further digested with 1.25% sodium hydroxide and filtered in course porous (75 μm) crucible in apparatus at a vacuum of about 25 mm. The residual left after refluxing was washed again with 1.25% sulfuric acid at near boiling point. Then, the residual was dried at 110°C overnight, cooled in a desiccator, and weighed (m2). After being dried, the sample was ashed at 550°C until the ashing was complete, cooled in a desiccator, and weighed again (m3). The total crude fiber was expressed in percentages as follows:(3)Fiber%=m2−m3m1∗100,where m1 is a mass of sample (g), m2 is mass of sample with crucible before ashing (g), and m3 is mass of sample with crucible after ashing (g).


## 2.3.4. Total Protein

A dried and grounded sample was taken (0.5 g) and added into aKjeldahl digestion flask. One gram of catalyst (Na2SO4 mixed with anhydrous CuSO4 in a ratio of 10 : 1) and 5 ml of concentrated H2SO4 were added into the digestion flask. Then, using a digester, the mixed sample was digested at 350°C for about two hours until the sample was completely digested. Then, the flask was removed from the digester and allowed to cool and the digested sample was diluted by adding 30 ml of distilled water. Then, 25 ml concentrated 40% NaOH was added into the digestion flask to neutralize the acid and make the solution slightly alkaline. The contents were immediately distilled by inserting the digestion tube line into the receiver flask that contained 25 ml of 4% boric acid solution and about 150 ml of distillate collected. Then, the distillate was titrated by a standard acid (0.1 N HCl). The percentage of crude protein was calculated by multiplying the nitrogen percentage by the conversion factor (6.25) [13].(4)Nitrogen%=VHCIA−VHCIB∗NHCI∗0.014Wsample∗100,where V = volume of standard acid used for titration of sample (A) and blank sample (B), N = normality of standard acid used for titration (0.1 N HCl), 0.014 is the molecular weight of nitrogen, and W = weight sample taken for digestion, on a dry basis.(5)Protein%=6.25∗%Nitrogen.


## 2.4. Analysis of Variance

The quantitative field data were subjected to analysis of variance (ANOVA) and computed with R statistical software agricolae package [14]. The biochemical traits were analyzed following the CRD (completely randomized design) procedure. The traits that exhibited significant mean squares in ANOVA were further subjected to multivariate analysis.


## 2.5. Principal Component Analysis

Principal component analysis (PCA) was computed to find out the traits, which accounted more for the total variation. The data were standardized to mean zero and variance of one before computing principal component analysis to avoid differences in measurement scales. The principal component based on the correlation matrix was calculated using the R statistical software FactoMineR package [15].


## 2.6. Euclidean Distance and Clustering of Genotypes

Euclidean distance (ED) was computed from quantitative after subtracting the mean value and dividing it by the standard deviation as established by Sneath and Sokal [16]. R statistical software factoextra package [17] was used for the analysis of distance matrix and constructing dendrogram. The dendrogram was constructed based on the unweighted pair group method with arithmetic mean (UPGMA) from the distance matrix of phenotypic traits.


## 3. Results and Discussion

### 3.1. Analysis of Variance

The results of the analysis of variance for phenology, growth, yield-related traits, and biochemical traits showed a significant (P<0.05) difference. However, the genotypes exhibited a nonsignificant difference for the number of flower epicalyx and fruit diameter (Tables 1 and 2).Table 1
Mean square from analysis variance following simple lattice design for 25 quantitative traits of 36 okra genotypes evaluated at Melkassa in 2018.
TraitsRep. (1)Blocks within rep. (Adj.) (10)Treatments (uadj.) (35)Intrablock error (25)Total (71)CV (%)Days to 50% emergence1.68060.11390.2964∗0.14720.23775.03Days to first flowering0.12502.991722.2679∗∗2.418312.2512.85Days to 50% flowering7.34721.230622.6806∗∗3.193912.5823.03Days to maturity24.5003.800025.1000∗∗4.260014.7532.60Stem diameter (cm)0.11840.15390.1401∗∗0.05340.111213.11Plant height (cm)1773.6459.731153.51∗∗253.55747.6413.94Number of branches0.53220.14841.9838∗∗0.14641.057917.17Number of internodes1.62004.06946.9937∗3.35895.226112.35Internode length (cm)3.57780.95205.0346∗∗0.66702.901215.32Leaf length (cm)0.87120.87128.0159∗3.67436.044313.30Leaf width (cm)13.58916.297121.1354∗∗8.726015.97814.64Number of epicalyx0.03831.86652.4630 ns1.26491.920010.14Peduncle length (cm)0.01680.07720.2807∗∗0.04020.163610.12Fruit length (cm)16.0936.525813.6482∗∗1.67818.46478.38Fruit diameter (cm)0.06210.25270.2067 ns0.11860.180112.08Number of fruit ridge0.52700.64051.6811∗∗0.40921.07049.06Fresh fruit weight (g)118.4932.4883190.64∗∗30.137110.8315.78Number of fruits per plant3.52898.157835.4541∗∗6.656921.0209.42Number of mature pods1.00353.657619.8511∗∗2.791811.29818.38Fruit yield per plant (kg)0.16630.07370.1513∗∗0.04550.103322.22Fruit yield (ton ha−1)84.99028.018863.921∗∗22.09644.43423.43Number of seeds per pod257.4274.1880365.41∗∗87.239224.929.66Hundred seeds’ weight (g)0.48680.19930.6721∗∗0.20040.43686.77Seed yield per plant (g)595.09335.18953.98∗∗170.33585.8522.29Seed yield per hectare (kg)258286145475414055∗∗73925.925426922.29∗∗, ∗ and ns, significant at P<0.01, P<0.05, and nonsignificant, respectively. Rep = replication, Adj = adjusted, Uadj = unadjusted, CV (%) = coefficient of variation in percent. The number in parenthesis in each source of variation represents the degree of freedom.Table 2
Mean square from analysis variance for four seed biochemical contents of 36 okra genotypes evaluated at Melkassa in 2018.
Source of variationDFTotal ash (%)Crude fat (%)Crude fiber (%)Total protein (%)Genotypes350.0948∗∗1.6531∗∗21.52∗∗52.53∗∗Error360.00380.50310.46400.0538CV (%)1.444.622.402.03∗∗, significant at P<0.01, DF = degree of freedom, and CV (%) = coefficient of variation in percent.


### 3.2. Principal Component Analysis

The result of principal component analysis for 27 quantitative traits is presented in Table3. With eigenvalues ranging from 1.033 to 8.247, the principal component analysis resulted in eight principal components (PC1 to PC8). The eight principal components each accounted for a different percentage of the total variance, ranging from 3.83 to 30.54%, for a total variance of 82.44%. The PCs with an eigenvalue of <1 were ignored due to Gutten’s lower bound principle that eigenvalues <1 should be ignored. The first principal component (PC1) contributed to most of the variation (30.54%), followed by PC2, PC3, and PC4, which contributed 14.11%, 10.87%, and 6.98%, of the variation respectively, and the first four PCs accounted for 62.51% of the total variation.Table 3
Factor loading of the first eight principal components for 27 quantitative traits of 36 okra genotypes.
TraitsPrincipal componentsPC1PC2PC3PC4PC5PC6PC7PC8Days to 50% emergence0.270.280.110.26−0.480.130.31−0.41Days to first flowering0.380.80−0.280.02−0.25−0.09−0.020.06Days to 50% flowering0.260.84−0.320.14−0.21−0.110.020.09Days to maturity0.360.75−0.390.10−0.31−0.130.000.07Stem diameter (cm)0.780.25−0.16−0.210.180.14−0.03−0.02Plant height (cm)0.580.12−0.40−0.330.340.080.100.23Number of branches0.470.18−0.26−0.090.470.080.58−0.15Number of internodes0.480.23−0.530.120.430.15−0.19−0.14Internode length (cm)0.46−0.21−0.28−0.33−0.31−0.190.160.31Leaf length (cm)0.77−0.240.120.040.18−0.160.28−0.24Leaf width (cm)0.80−0.270.280.020.040.090.18−0.07Peduncle length (cm)0.750.010.01−0.16−0.300.13−0.03−0.16Fruit length (cm)0.38−0.140.55−0.11−0.40−0.28−0.010.03Number of fruit ridge0.110.510.43−0.01−0.120.39−0.15−0.25Fresh fruit weight (g)0.680.260.550.110.25−0.03−0.160.04Number of fruit per plant0.60−0.48−0.350.01−0.16−0.240.03−0.21Number of mature pod0.59−0.56−0.370.22−0.11−0.10−0.23−0.14Fruit yield per plant (kg)0.810.160.350.160.21−0.17−0.19−0.02Fruit yield (t ha−1)0.810.180.330.190.18−0.16−0.19−0.01Number of seeds per pod0.25−0.050.24−0.15−0.160.800.120.24Hundred seeds’ weight (g)0.550.150.18−0.38−0.09−0.12−0.020.39Seed yield per plant (g)0.75−0.48−0.260.10−0.130.21−0.180.09Seed yield (kg ha−1)0.75−0.48−0.260.10−0.130.21−0.180.09Total ash (%)0.350.090.55−0.530.16−0.18−0.02−0.06Crude fat (%)0.36−0.210.180.24−0.30−0.040.440.17Crude fiber (%)−0.020.010.180.650.29−0.090.310.34Total protein (%)0.250.010.150.64−0.010.10−0.140.25Eigenvalue8.253.812.941.881.831.341.181.03Contribution (%)30.5414.1110.876.986.774.974.363.82Cumulative variance (%)30.5444.6655.5362.5169.2874.2578.6182.44PC1 to PC8 represent the first principle component to the eighth principal component.A similar result on okra was reported by Muluken et al. [18] in which the first three principal components PC1, PC2, and PC3, with values of 32.4%, 16.7%, and 8.2%, respectively, contributed more to the total of 57.3% variation. Amoatey et al. [19] reported the first, second, and third principal components with values of 32.44%, 19.78%, and 9.68% of the total genetic variation, respectively. Ahiakpa [20] also reported that the first principal component (PC1) was (32.44%) the major contributor for variance in okra genotypes.Within the PC1, traits with the largest values closer to one influence the cluster more than traits with lower absolute values closer to zero [21]. Hence, the differentiation of the genotypes into different clusters was because of the cumulative effect of several traits rather than the large contribution of a few traits. In this regard, stem diameter (7.41%), fruit yield per hectare (8.03%), leaf width (7.56%), fruit yield per plant (7.87%), leaf length (7.11%), peduncle length (6.75), seed yield per plant (6.75%), and seed yield per hectare (6.75%) had relatively higher contributions to PC1. This indicates that these traits were responsible for the differentiation of the clusters and had a greater contribution to the total diversity. In PC2, days to 50% flowering (18.58%), days to first flowering (16.87%), date of maturity (14.76%), and the number of mature pods (8.12%) had more contribution, whereas fruit length, fresh fruit weight, and ash content had relatively more contribution in PC3 (Table 3).A biplot was performed based on the first two PCs (Figure1). The genotypes and quantitative traits were shown on a biplot to visualize their associations. The first and the second PC biplots explained 44.66% of the total variability among the genotypes, displaying that stem diameter, leaf length, and leaf width, fruit yield per plant, fruit yield per hectare, seed yield per plant, and seed yield per hectare were considered the most discriminating traits.Figure 1
Biplot (axes PC1 and PC2) of 27 quantitative traits of 36 okra genotypes.The genotypes positioned on the right top quadrant were characterized by late maturity, high fresh fruit weight, much fruit ridges, and high stem diameter. The genotypes depicted in the bottom right quadrant had the highest seed yield, number of fruits per plant, number of mature pods, and longest, and widest leaf. The genotypes distributed around the origin had similar genetic characteristics, while the genotypes that were found far from the origin are considered unrelated genotypes (Figure1). Therefore, these divergent genotypes could be used as potential parents for successful hybridization to develop heterotic groups in the okra-breeding program.PC3 and PC4 biplots are presented in Figure2. These two PCs accounted for 17.85% of the total variability among genotypes, showing that ash content, and total protein content, crude fiber content, fruit length, and fresh fruit weight were the most contributing traits. PC3 and PC4 biplots provided information regarding the similarities and the pattern of differences among the okra genotypes and the association between traits. Genotypes were distributed in all four quadrants on the axes, indicating the presence of wide genetic variability for the traits studied. Overlapped accessions and accessions closer to each other on the axes had similar genetic makeup. However, genotypes that are apart from each other could be considered genetically distinct. Genotypes positioned in the top right quadrant were characterized by high-seed protein and fiber content. The top left quadrant consists of the okra genotypes that are closely related and have a high number of internode and a high number of mature pods. Genotypes found on the right bottom quadrant exhibited the highest seed ash content.Figure 2
Biplot (axes PC3 and PC4) of 27 quantitative traits of 36 okra genotypes.


### 3.3. Cluster Analysis

The optimum number of clusters was determined by the total within sum of square (WSS) (elbow method) using R statistical software version 3.6.3 (Figure3). A dendrogram was constructed based on the unweighted pair group method with arithmetic mean (average) from the distance matrix of phenotypic traits.Figure 3
Determination of the optimum number of clusters K using total within sum of square (WSS) method.The distances of all possible pairs of the 36 okra genotypes from 27 quantitative traits were estimated by Euclidean distance. The distances between okra genotypes ranged from 3.55 to 14.49 with a mean, standard deviation, and coefficient of variation of 7.12, 1.80, and 25.25%, respectively. The highest genetic distance (Euclidean distance) was computed between 29407 and Humera 2 (14.49), whereas the lowest genetic distance was estimated between Dhenu and Mithra (3.55) (Figure4).Figure 4
Scatter plot constructed based on PC1 and PC2 for 27 quantitative traits.Based on PC axes 1 and 2, a scatter plot was constructed for four clusters (Figure4). The plot showed that the genotypes that have similar genetic makeup were grouped in a cluster (near to overlap), and those genotypes that have different genetics were positioned in the opposite corner of the scatter plot.Generally, the Euclidean distances measured among the introduced varieties were lower than the genetic distances among genotypes collected from Ethiopia. This showed that there is a higher chance of improving fruit yield and seed-related traits through the selection and/or hybridization of okra genotypes collected from different okra growing regions of Ethiopia.By characterizing 24 Ethiopian okra genotypes, Fozia [22] reported Euclidean distance that ranged from 1.96 to 11.36 with a mean, standard deviation, and coefficient of variation of 5.85, 1.97, and 33.75%, respectively. The same study also reported that introduced (Indian) varieties had lower (1.96 to 10.01) genetic distances than Ethiopia's okra collection, which ranged from (2.07 to 11.36). Muluken et al. [18] reported that Ethiopian okra collections exhibited a wider genetic distance than exotic varieties. Anteneh [23] estimated the genetic distances of all possible pairs of 25 okra genotypes and reported that the highest genetic distances were observed between okra collections from Ethiopia and introduced commercial varieties from other countries, while the lowest genetic distance was estimated between introduced commercial varieties.The extent of diversity present between genotypes determines the extent of improvement gained through selection and hybridization. The more distant the two genotypes are, the greater the probability of improvement through selection and hybridization. Mihretu et al. [3] also reported the presence of considerable genetic distance among okra collections from Gambela regional state, which is one of the okra-growing regions in Ethiopia.Clustering of genotypes based on Euclidean distances revealed four major clusters. The number and names of genotypes in each cluster along with their collection origin are presented in Table4. Cluster I consists of the majority of the genotypes, which accounted for 24 (66.66%) of the genotypes. Cluster II contains 10 (27.77%) genotypes, while clusters III and IV each contain only single genotypes (Table 4, Figure 5).Table 4
Clusters, number of accessions, name of genotypes, and collection origin of 36 okra genotypes.
ClustersNumber of genotypesGenotypesCollection origin (district)I24Namadhari, SOH 701, 240593, 240581, 240608, Mithra, Dhenu, 240783, 242436, 240613, 242434, 242440, 240782, 240596, 240598, 240614, 240603, 240612, 240594, Kraft, 242589, Humera 1, 240610, Bamaya HumeraSouthwestern Ethiopia (Gog), southwestern Ethiopia (Gambela), southwestern Ethiopia (Itang), northwestern Ethiopia (Dangur), western Ethiopia (Asosa), northwestern Ethiopia (Metema) southwestern Ethiopia (Akobo)II10Arka Anamica, Pocha, 242438, Kiran, Clemson Spinless, 242446, Humera 2, 242442, 242439, 240584Indian commercial variety, American commercial variety, western Ethiopia (Asosa), southwestern Ethiopia (Akobo)III1242450Western Ethiopia (AsosaIV129407Western Ethiopia (Kurmuk)Figure 5
Dendrogram illustrating dissimilarity of 36 okra genotypes by an average method of clustering method from Euclidean distance matrix based on 27 traits.Genotypes clustered in cluster I and cluster II were early maturing, while the two genotypes positioned in cluster III and cluster IV were late-maturing genotypes (Table5). Therefore, genotypes found in cluster I and cluster II could be used for okra production in areas characterized by the low amount of rainfall. These genotypes could also be used by breeders for developing varieties suitable for drought-prone areas. On the contrary, the two genotypes found in clusters III and IV could be used for areas that have a long rainfall season. The highest mean of fruit yield was measured from the genotype (29407), which is found in cluster IV. This genotype also had the highest seed protein content. This genotype could be used for further evaluation to identify genotypes with high fruit yield and desirable nutrient contents.Table 5
Cluster mean value for 27 quantitative traits of 36 okra genotypes.
TraitsClustersIIIIIIIVDays to 50% emergence7.637.607.508.00Days to first flowering54.4653.6059.0061.50Days to 50% flowering58.4858.5065.0065.00Days to maturity79.0878.3085.0087.00Stem diameter (cm)1.841.521.922.27Plant height (cm)120.6593.84162.03113.92Number of branches2.321.754.133.00Number of internodes15.1213.4720.9315.47Internode length (cm)5.694.474.975.85Leaf length (cm)15.0212.4015.1719.00Leaf width (cm)21.2817.3714.8427.10Peduncle length (cm)2.101.641.752.82Fruit length (cm)15.8914.358.2522.87Number of fruit ridge7.146.777.007.91Fresh fruit weight (g)37.4027.0920.8262.41Number of fruit per plant28.5524.0229.5431.46Number of mature pods10.046.559.5011.25Fruit yield per plant (kg)1.040.690.771.87Fruit yield (ton ha−1)21.7414.5815.9438.86Number of seeds per pod100.1192.5273.4077.60Hundred seeds’ weight (g)6.766.285.637.37Seed yield per plant (g)67.3838.5239.2265.22Seed yield per hectare (kg)1403.85802.55817.131358.82Total ash (%)4.354.183.964.63Crude fat (%)15.5814.8813.1416.59Crude fiber (%)29.9431.2928.9932.72Total protein (%)10.7111.376.5620.94


## 3.1. Analysis of Variance

The results of the analysis of variance for phenology, growth, yield-related traits, and biochemical traits showed a significant (P<0.05) difference. However, the genotypes exhibited a nonsignificant difference for the number of flower epicalyx and fruit diameter (Tables 1 and 2).Table 1
Mean square from analysis variance following simple lattice design for 25 quantitative traits of 36 okra genotypes evaluated at Melkassa in 2018.
TraitsRep. (1)Blocks within rep. (Adj.) (10)Treatments (uadj.) (35)Intrablock error (25)Total (71)CV (%)Days to 50% emergence1.68060.11390.2964∗0.14720.23775.03Days to first flowering0.12502.991722.2679∗∗2.418312.2512.85Days to 50% flowering7.34721.230622.6806∗∗3.193912.5823.03Days to maturity24.5003.800025.1000∗∗4.260014.7532.60Stem diameter (cm)0.11840.15390.1401∗∗0.05340.111213.11Plant height (cm)1773.6459.731153.51∗∗253.55747.6413.94Number of branches0.53220.14841.9838∗∗0.14641.057917.17Number of internodes1.62004.06946.9937∗3.35895.226112.35Internode length (cm)3.57780.95205.0346∗∗0.66702.901215.32Leaf length (cm)0.87120.87128.0159∗3.67436.044313.30Leaf width (cm)13.58916.297121.1354∗∗8.726015.97814.64Number of epicalyx0.03831.86652.4630 ns1.26491.920010.14Peduncle length (cm)0.01680.07720.2807∗∗0.04020.163610.12Fruit length (cm)16.0936.525813.6482∗∗1.67818.46478.38Fruit diameter (cm)0.06210.25270.2067 ns0.11860.180112.08Number of fruit ridge0.52700.64051.6811∗∗0.40921.07049.06Fresh fruit weight (g)118.4932.4883190.64∗∗30.137110.8315.78Number of fruits per plant3.52898.157835.4541∗∗6.656921.0209.42Number of mature pods1.00353.657619.8511∗∗2.791811.29818.38Fruit yield per plant (kg)0.16630.07370.1513∗∗0.04550.103322.22Fruit yield (ton ha−1)84.99028.018863.921∗∗22.09644.43423.43Number of seeds per pod257.4274.1880365.41∗∗87.239224.929.66Hundred seeds’ weight (g)0.48680.19930.6721∗∗0.20040.43686.77Seed yield per plant (g)595.09335.18953.98∗∗170.33585.8522.29Seed yield per hectare (kg)258286145475414055∗∗73925.925426922.29∗∗, ∗ and ns, significant at P<0.01, P<0.05, and nonsignificant, respectively. Rep = replication, Adj = adjusted, Uadj = unadjusted, CV (%) = coefficient of variation in percent. The number in parenthesis in each source of variation represents the degree of freedom.Table 2
Mean square from analysis variance for four seed biochemical contents of 36 okra genotypes evaluated at Melkassa in 2018.
Source of variationDFTotal ash (%)Crude fat (%)Crude fiber (%)Total protein (%)Genotypes350.0948∗∗1.6531∗∗21.52∗∗52.53∗∗Error360.00380.50310.46400.0538CV (%)1.444.622.402.03∗∗, significant at P<0.01, DF = degree of freedom, and CV (%) = coefficient of variation in percent.


## 3.2. Principal Component Analysis

The result of principal component analysis for 27 quantitative traits is presented in Table3. With eigenvalues ranging from 1.033 to 8.247, the principal component analysis resulted in eight principal components (PC1 to PC8). The eight principal components each accounted for a different percentage of the total variance, ranging from 3.83 to 30.54%, for a total variance of 82.44%. The PCs with an eigenvalue of <1 were ignored due to Gutten’s lower bound principle that eigenvalues <1 should be ignored. The first principal component (PC1) contributed to most of the variation (30.54%), followed by PC2, PC3, and PC4, which contributed 14.11%, 10.87%, and 6.98%, of the variation respectively, and the first four PCs accounted for 62.51% of the total variation.Table 3
Factor loading of the first eight principal components for 27 quantitative traits of 36 okra genotypes.
TraitsPrincipal componentsPC1PC2PC3PC4PC5PC6PC7PC8Days to 50% emergence0.270.280.110.26−0.480.130.31−0.41Days to first flowering0.380.80−0.280.02−0.25−0.09−0.020.06Days to 50% flowering0.260.84−0.320.14−0.21−0.110.020.09Days to maturity0.360.75−0.390.10−0.31−0.130.000.07Stem diameter (cm)0.780.25−0.16−0.210.180.14−0.03−0.02Plant height (cm)0.580.12−0.40−0.330.340.080.100.23Number of branches0.470.18−0.26−0.090.470.080.58−0.15Number of internodes0.480.23−0.530.120.430.15−0.19−0.14Internode length (cm)0.46−0.21−0.28−0.33−0.31−0.190.160.31Leaf length (cm)0.77−0.240.120.040.18−0.160.28−0.24Leaf width (cm)0.80−0.270.280.020.040.090.18−0.07Peduncle length (cm)0.750.010.01−0.16−0.300.13−0.03−0.16Fruit length (cm)0.38−0.140.55−0.11−0.40−0.28−0.010.03Number of fruit ridge0.110.510.43−0.01−0.120.39−0.15−0.25Fresh fruit weight (g)0.680.260.550.110.25−0.03−0.160.04Number of fruit per plant0.60−0.48−0.350.01−0.16−0.240.03−0.21Number of mature pod0.59−0.56−0.370.22−0.11−0.10−0.23−0.14Fruit yield per plant (kg)0.810.160.350.160.21−0.17−0.19−0.02Fruit yield (t ha−1)0.810.180.330.190.18−0.16−0.19−0.01Number of seeds per pod0.25−0.050.24−0.15−0.160.800.120.24Hundred seeds’ weight (g)0.550.150.18−0.38−0.09−0.12−0.020.39Seed yield per plant (g)0.75−0.48−0.260.10−0.130.21−0.180.09Seed yield (kg ha−1)0.75−0.48−0.260.10−0.130.21−0.180.09Total ash (%)0.350.090.55−0.530.16−0.18−0.02−0.06Crude fat (%)0.36−0.210.180.24−0.30−0.040.440.17Crude fiber (%)−0.020.010.180.650.29−0.090.310.34Total protein (%)0.250.010.150.64−0.010.10−0.140.25Eigenvalue8.253.812.941.881.831.341.181.03Contribution (%)30.5414.1110.876.986.774.974.363.82Cumulative variance (%)30.5444.6655.5362.5169.2874.2578.6182.44PC1 to PC8 represent the first principle component to the eighth principal component.A similar result on okra was reported by Muluken et al. [18] in which the first three principal components PC1, PC2, and PC3, with values of 32.4%, 16.7%, and 8.2%, respectively, contributed more to the total of 57.3% variation. Amoatey et al. [19] reported the first, second, and third principal components with values of 32.44%, 19.78%, and 9.68% of the total genetic variation, respectively. Ahiakpa [20] also reported that the first principal component (PC1) was (32.44%) the major contributor for variance in okra genotypes.Within the PC1, traits with the largest values closer to one influence the cluster more than traits with lower absolute values closer to zero [21]. Hence, the differentiation of the genotypes into different clusters was because of the cumulative effect of several traits rather than the large contribution of a few traits. In this regard, stem diameter (7.41%), fruit yield per hectare (8.03%), leaf width (7.56%), fruit yield per plant (7.87%), leaf length (7.11%), peduncle length (6.75), seed yield per plant (6.75%), and seed yield per hectare (6.75%) had relatively higher contributions to PC1. This indicates that these traits were responsible for the differentiation of the clusters and had a greater contribution to the total diversity. In PC2, days to 50% flowering (18.58%), days to first flowering (16.87%), date of maturity (14.76%), and the number of mature pods (8.12%) had more contribution, whereas fruit length, fresh fruit weight, and ash content had relatively more contribution in PC3 (Table 3).A biplot was performed based on the first two PCs (Figure1). The genotypes and quantitative traits were shown on a biplot to visualize their associations. The first and the second PC biplots explained 44.66% of the total variability among the genotypes, displaying that stem diameter, leaf length, and leaf width, fruit yield per plant, fruit yield per hectare, seed yield per plant, and seed yield per hectare were considered the most discriminating traits.Figure 1
Biplot (axes PC1 and PC2) of 27 quantitative traits of 36 okra genotypes.The genotypes positioned on the right top quadrant were characterized by late maturity, high fresh fruit weight, much fruit ridges, and high stem diameter. The genotypes depicted in the bottom right quadrant had the highest seed yield, number of fruits per plant, number of mature pods, and longest, and widest leaf. The genotypes distributed around the origin had similar genetic characteristics, while the genotypes that were found far from the origin are considered unrelated genotypes (Figure1). Therefore, these divergent genotypes could be used as potential parents for successful hybridization to develop heterotic groups in the okra-breeding program.PC3 and PC4 biplots are presented in Figure2. These two PCs accounted for 17.85% of the total variability among genotypes, showing that ash content, and total protein content, crude fiber content, fruit length, and fresh fruit weight were the most contributing traits. PC3 and PC4 biplots provided information regarding the similarities and the pattern of differences among the okra genotypes and the association between traits. Genotypes were distributed in all four quadrants on the axes, indicating the presence of wide genetic variability for the traits studied. Overlapped accessions and accessions closer to each other on the axes had similar genetic makeup. However, genotypes that are apart from each other could be considered genetically distinct. Genotypes positioned in the top right quadrant were characterized by high-seed protein and fiber content. The top left quadrant consists of the okra genotypes that are closely related and have a high number of internode and a high number of mature pods. Genotypes found on the right bottom quadrant exhibited the highest seed ash content.Figure 2
Biplot (axes PC3 and PC4) of 27 quantitative traits of 36 okra genotypes.


## 3.3. Cluster Analysis

The optimum number of clusters was determined by the total within sum of square (WSS) (elbow method) using R statistical software version 3.6.3 (Figure3). A dendrogram was constructed based on the unweighted pair group method with arithmetic mean (average) from the distance matrix of phenotypic traits.Figure 3
Determination of the optimum number of clusters K using total within sum of square (WSS) method.The distances of all possible pairs of the 36 okra genotypes from 27 quantitative traits were estimated by Euclidean distance. The distances between okra genotypes ranged from 3.55 to 14.49 with a mean, standard deviation, and coefficient of variation of 7.12, 1.80, and 25.25%, respectively. The highest genetic distance (Euclidean distance) was computed between 29407 and Humera 2 (14.49), whereas the lowest genetic distance was estimated between Dhenu and Mithra (3.55) (Figure4).Figure 4
Scatter plot constructed based on PC1 and PC2 for 27 quantitative traits.Based on PC axes 1 and 2, a scatter plot was constructed for four clusters (Figure4). The plot showed that the genotypes that have similar genetic makeup were grouped in a cluster (near to overlap), and those genotypes that have different genetics were positioned in the opposite corner of the scatter plot.Generally, the Euclidean distances measured among the introduced varieties were lower than the genetic distances among genotypes collected from Ethiopia. This showed that there is a higher chance of improving fruit yield and seed-related traits through the selection and/or hybridization of okra genotypes collected from different okra growing regions of Ethiopia.By characterizing 24 Ethiopian okra genotypes, Fozia [22] reported Euclidean distance that ranged from 1.96 to 11.36 with a mean, standard deviation, and coefficient of variation of 5.85, 1.97, and 33.75%, respectively. The same study also reported that introduced (Indian) varieties had lower (1.96 to 10.01) genetic distances than Ethiopia's okra collection, which ranged from (2.07 to 11.36). Muluken et al. [18] reported that Ethiopian okra collections exhibited a wider genetic distance than exotic varieties. Anteneh [23] estimated the genetic distances of all possible pairs of 25 okra genotypes and reported that the highest genetic distances were observed between okra collections from Ethiopia and introduced commercial varieties from other countries, while the lowest genetic distance was estimated between introduced commercial varieties.The extent of diversity present between genotypes determines the extent of improvement gained through selection and hybridization. The more distant the two genotypes are, the greater the probability of improvement through selection and hybridization. Mihretu et al. [3] also reported the presence of considerable genetic distance among okra collections from Gambela regional state, which is one of the okra-growing regions in Ethiopia.Clustering of genotypes based on Euclidean distances revealed four major clusters. The number and names of genotypes in each cluster along with their collection origin are presented in Table4. Cluster I consists of the majority of the genotypes, which accounted for 24 (66.66%) of the genotypes. Cluster II contains 10 (27.77%) genotypes, while clusters III and IV each contain only single genotypes (Table 4, Figure 5).Table 4
Clusters, number of accessions, name of genotypes, and collection origin of 36 okra genotypes.
ClustersNumber of genotypesGenotypesCollection origin (district)I24Namadhari, SOH 701, 240593, 240581, 240608, Mithra, Dhenu, 240783, 242436, 240613, 242434, 242440, 240782, 240596, 240598, 240614, 240603, 240612, 240594, Kraft, 242589, Humera 1, 240610, Bamaya HumeraSouthwestern Ethiopia (Gog), southwestern Ethiopia (Gambela), southwestern Ethiopia (Itang), northwestern Ethiopia (Dangur), western Ethiopia (Asosa), northwestern Ethiopia (Metema) southwestern Ethiopia (Akobo)II10Arka Anamica, Pocha, 242438, Kiran, Clemson Spinless, 242446, Humera 2, 242442, 242439, 240584Indian commercial variety, American commercial variety, western Ethiopia (Asosa), southwestern Ethiopia (Akobo)III1242450Western Ethiopia (AsosaIV129407Western Ethiopia (Kurmuk)Figure 5
Dendrogram illustrating dissimilarity of 36 okra genotypes by an average method of clustering method from Euclidean distance matrix based on 27 traits.Genotypes clustered in cluster I and cluster II were early maturing, while the two genotypes positioned in cluster III and cluster IV were late-maturing genotypes (Table5). Therefore, genotypes found in cluster I and cluster II could be used for okra production in areas characterized by the low amount of rainfall. These genotypes could also be used by breeders for developing varieties suitable for drought-prone areas. On the contrary, the two genotypes found in clusters III and IV could be used for areas that have a long rainfall season. The highest mean of fruit yield was measured from the genotype (29407), which is found in cluster IV. This genotype also had the highest seed protein content. This genotype could be used for further evaluation to identify genotypes with high fruit yield and desirable nutrient contents.Table 5
Cluster mean value for 27 quantitative traits of 36 okra genotypes.
TraitsClustersIIIIIIIVDays to 50% emergence7.637.607.508.00Days to first flowering54.4653.6059.0061.50Days to 50% flowering58.4858.5065.0065.00Days to maturity79.0878.3085.0087.00Stem diameter (cm)1.841.521.922.27Plant height (cm)120.6593.84162.03113.92Number of branches2.321.754.133.00Number of internodes15.1213.4720.9315.47Internode length (cm)5.694.474.975.85Leaf length (cm)15.0212.4015.1719.00Leaf width (cm)21.2817.3714.8427.10Peduncle length (cm)2.101.641.752.82Fruit length (cm)15.8914.358.2522.87Number of fruit ridge7.146.777.007.91Fresh fruit weight (g)37.4027.0920.8262.41Number of fruit per plant28.5524.0229.5431.46Number of mature pods10.046.559.5011.25Fruit yield per plant (kg)1.040.690.771.87Fruit yield (ton ha−1)21.7414.5815.9438.86Number of seeds per pod100.1192.5273.4077.60Hundred seeds’ weight (g)6.766.285.637.37Seed yield per plant (g)67.3838.5239.2265.22Seed yield per hectare (kg)1403.85802.55817.131358.82Total ash (%)4.354.183.964.63Crude fat (%)15.5814.8813.1416.59Crude fiber (%)29.9431.2928.9932.72Total protein (%)10.7111.376.5620.94


## 4. Conclusions

The results showed the presence of considerable genetic diversity for the studied morphological and biochemical traits. This variation could be exploited to develop varieties with different desirable agronomic traits like early maturing, high yield, and good nutrient content through either selection and/or hybridization using the okra genotypes collected and conserved in Ethiopia.In addition, the study revealed the potential of the landrace okra genotypes as sources of nutrients. This indicates the importance of neglected crops that could be utilized for ensuring food security and alleviating malnutrition in developing countries like Ethiopia, where malnutrition is a widespread problem. It is also recommended to extend the research in okra to include micronutrient content analysis, molecular diversity study, and sequencing okra genotypes to identify important agronomic and biochemical traits and to characterize the genes responsible for the traits.



---

*Source: 3306793-2022-03-12.xml*