The overall mean for various traits across the broad bean genotypes indicated a considerable range of diversity among the genotypes for each trait (Table 1). The plant height varied from 72.67 to 113.67 cm, Days to 50% flowering with 52.71 to 63.97days, days to maturity varied from 147.00 to 163.50 days. The average number of pods per plant was 17.25. The mean values of 100 seed weight and seed yield per plant are 109.78 g and 64.73 g respectively. The average pod yield per plant is 160gram, protein is 21.53%. The maximum coefficient of variation (10.87) is observed for number of pods per plant. These results highlight that there was substantial variability within the material, making it suitable for further investigations into various traits. Additionally, the coefficient of variation was found to be below 10 percent for all the characters assessed, indicating that the experiment was carried out with precision.
Abbas and Rodany (2009);
Ammar et al., (2015); Arab et al., (2018); Boghara et al., (2016); Dahiya et al., (2021) and
Ertiro et al., (2023) have also reported similar observations for analysis of variance in broad beans.
Genetic divergence analysis
The data collected on different morphological and qualitative characters of broad bean genotypes underwent an analysis for genetic diversity based on Mahalanobis D
2 statistic and principal component analysis. The significant values indicated substantial variability within the studied genotypes, resulting in genetic diversity. The percent contribution of all the 19 characters towards genetic divergence is presented Fig 1 based on the number of times that character appeared in the first rank. The maximum contribution towards diversity was shown by Total soluble solid (44.74%) followed by carbohydrate (41.05%), moisture (4.21%), phenol (2.63%), pod length (2.11%), ascorbic acid (2.11%), 100 seed weight (1.05%), dry matter (1.05%), protein (0.53%) and number of seeds per pod (0.53%) whereas the contribution of other traits was zero per cent.
Jadhav et al., (2023) also reported similar contribution of the traits towards divergence.
Using Tocher’s method, 20 genotypes were categorized into nine clusters based on the relative magnitude of their D
2 estimates. Table 2 illustrates the distribution of various genotypes within each cluster. Among the nine clusters, cluster II was the largest cluster with 4 genotypes followed by the clusters I, III and IV each with three genotypes. Cluster V and VI comprised of 2 genotypes, whereas the clusters VII, VIII and XI are monogenotypic consisting of only one genotype. The diversity within the genotypes is apparent due to the cluster arrangement, indicating significant variability. The findings broadly align with the reports by
Rebaa et al., (2017) and
Dewangan et al., (2022).
The intra and inter-cluster D
2 values among nine clusters were calculated and presented in Table 3. The maximum intra-cluster distance was recorded in cluster IV (3358.2) followed by cluster II (D
2=2899.7), cluster III (D
2=2359.3) and cluster VI (D
2=2298.6), whereas the least intra-cluster distance was observed in cluster I (D
2=2258.9) and cluster V (D
2=1805.8). Cluster VII, VIII and IX had no intra-cluster distance(D
2=0.0) as they contained only one genotype each. Genotypes clustered together are expected to display minimal variation among the measured characteristics. The largest inter-cluster D
2 values were found between cluster I and VII (316.9) and cluster I and III (315.9), indicating a broader genetic diversity among the genotypes in these clusters. The smallest inter-cluster distance (62.7) was observed between cluster VI and cluster XI, as well as between cluster VII and Cluster VIII (64.6), indicating a strong similarity or close relationship among the genotypes in these clusters.
The findings revealed that inter-cluster distances surpassed intra-cluster distances, signifying substantial genetic variations between clusters compared to the limited variations within clusters. The genetic diversity was notably higher in clusters III, V, VII and VIII compared to the remaining clusters, as they displayed the greatest genetic distance. The genotypes within these clusters might be valuable for future breeding programs, with the potential to produce improved offspring in subsequent generations. The results in similar lines were also reported by
Singh et al., (2015) and
Singh et al., (2021). The nearest and farthest clusters from each cluster can also be interpreted based on D
2 values in Table 3.
The cluster means of all the 19 phenotypic and qualitative characters were presented in Table 4. This cluster means value indicates the average performance of all genotypes present in a particular cluster for different traits. The data indicated a wide range of mean values between the clusters.The highest cluster mean value for number of pods per cluster (3.12), number of pods per plant (27.50), seed yield per plant (78.67), pod yield per plant (251.39), protein (24.47) and carbohydrate (52.18) were recorded in cluster III, phenol content (0.19) in cluster IV, number of branches (3.35) in cluster V, for total soluble solids (8.17) in cluster VI, pod length (12.85)in cluster VII, number of seeds per pod (3.25) and ascorbic acid (18.36); similarly for pod width (1.76), 100-seed weight (134.00) and dry matter (88.55) in cluster XI.
The cluster I had the lowest mean values for plant height (74.72). The lowest mean values for days to 50% flowering (54.34) was observed in cluster VIII, days to maturity (150.00) in cluster XI, moisture content (11.25) in cluster XI. Comparative study of cluster mean values suggested that clusters III, V and XI had highest cluster means for pod yield and its contributing traits, therefore, these clusters may be considered superior for selecting promising genotypes for utilizing in future studies. These results largely align with the discoveries made by
Dubey et al., (2022) and
Dewangan et al., (2022).
The top-performing genotypes selected from various clusters, considering multiple characteristics is given in Table 5. The promising broad bean genotypes identified from both the divergence and cluster mean analysis were Vf 2, Vf 8, Vf 1 from cluster III; Vf 9 and Vf 4 from cluster V; Vf 10 and Vf 14 form cluster VI and Vf 19 from cluster XI based on various pod yield and its component traits. These genotypes hold potential for future broad bean breeding programs, enabling the selection, hybridization and recovery of transgressive segregants with the highest yield potential.
Principal component analysis (PCA)
Principal component analysis (PCA) is a commonly employed statistical method for assessing genetic diversity within plant genotypes and identifying the key variables that contribute to this variation
(Price et al., 2006). The variability within the genotypes was distributed across 18 principal components (PCs). Out of these 18 PCs, four had eigenvalues exceeding 1, collectively representing 78.23 per cent of the total variation for yield and its contributing traits among the 20 genotypes. This indicates that the traits identified within these components notably influenced the genotypes. It was in accordance with the findings of
Tiwari and Singh (2019) and
Dewangan et al., (2022).
In a scree plot graph, eigenvalues associated with a factor were plotted in descending order against the number of principal components to illustrate the proportion of variance attributed to these components. The scree plot revealed that, apart from the first four principal components, the rest of the PCs exhibited a negligible amount of variation (Fig 2). The scree plot clearly showed that PC1 held the most significant level of variation, suggesting that choosing genotypes based on traits from PC1 would be advantageous.
The first principal component (PC 1) accounted for approximately 50.29 per cent of the total variation and possessed an eigen value of 9.556. Notably, traits such as plant height (0.697), number of branches per plant (0.745), pod width (0.776), number of pods per cluster (0.690), number of pods per plant (0.842), number of seeds per pod (0.729), 100 seeds weight (0.845), seed yield per plant (0.919), pod yield per plant (0.937), dry matter content (0.711), carbohydrate (0.881) and total soluble solids (0.684) made the most significant contributions to the divergence observed in PC 1. The second principal component (PC 2) possessed an eigen value of 2.884 and explained approximately 15.18 per cent of the total variation. Key contributors to PC 2 included traits like days to 50% flowering (0.720), days to maturity (0.696), dry matter (0.608) and protein content (0.685). The third principal component (PC 3), with an eigen value of 1.370, contributed about 7.21 per cent of the total variation. The major contributor to PC 3 was pod length (0.638). Meanwhile, the fourth principal component (PC 4) contributed roughly 5.54 per cent of the total variation and had an eigen value of 1.053 with pod width (0.413) contributing to the total variation.Similar findings were reported in the studies conducted by
Zayed et al., (2022) and
Dewangan et al., (2022).
The first two principal components were very important and accounted for approximately 65.47% of the total available variability; hence, they were plotted graphically to demonstrate the relationship among the broad bean genotypes. For easy visualization, both genotypes and variables were combined into a single bi-plot graph, which centered on two primary principal components, namely PC1 and PC2 (Fig 3). The length of each trait2 s vector within the graph illustrates its contribution to overall divergence. The trait pod yield per plant displayed the longest vector length, indicating its substantial contribution to overall divergence, which was followed by seed yield per plant, dry matter, carbohydrate, protein and number of pods per plant. These observations concur with those of
Malek et al., (2021).
The angle formed between the trait vectors signifies the direction of association between the traits. It reflected that pod yield per plant, seed yield per plant, pod width and number of pods per plant mostly contributed towards variability among genotypes
viz., Vf 1 and Vf 2. On the other hand, traits such as number of branches per plant and number of pods per cluster contributed significantly to the variability among genotypes Vf 8, Vf 10 and Vf 19. The genotypes Vf 1, Vf 2, Vf 8, Vf 7, Vf 20, Vf 6, Vf 14, Vf 11, Vf 11 and Vf 19 were positioned farthest from the origin in the biplot, signifying that they accounted for more trait variability associated with their respective principal components compared to other genotypes. These outcomes are consistent with the findings of
Tiwari and Singh (2019).
Stepwise multiple regression analysis
The pod yield per plant was regarded as the dependent variable, while the remaining traits served as independent variables. Initially, the model included the seed yield per plant, elucidating 79.1% of the total observed variations, followed by the number of pods per plant. The cumulative variations explained by the combination of the number of pods per plant with the seed yield per plant amount to 85.5% (Table 6). Thus, the stepwise multiple regression analysis identified the crucial economic traits on pod yield of broad beans. These findings could be valuable in effectively selecting promising broad bean genotypes in subsequent generations which agree with those of
Janghel et al., (2020) and
Dewangan et al., (2022).