Legume Research

  • Chief EditorJ. S. Sandhu

  • Print ISSN 0250-5371

  • Online ISSN 0976-0571

  • NAAS Rating 6.80

  • SJR 0.391

  • Impact Factor 0.8 (2024)

Frequency :
Monthly (January, February, March, April, May, June, July, August, September, October, November and December)
Indexing Services :
BIOSIS Preview, ISI Citation Index, Biological Abstracts, Elsevier (Scopus and Embase), AGRICOLA, Google Scholar, CrossRef, CAB Abstracting Journals, Chemical Abstracts, Indian Science Abstracts, EBSCO Indexing Services, Index Copernicus
Legume Research, volume 47 issue 8 (august 2024) : 1274-1281

Diversity Analysis for Phenotypic and Qualitative Traits of Broad Bean (Vicia faba L.): an Underutilized Vegetable Crop of India by Multivariate Analysis

S. P. Thamaraiselvi 1, A. Ajay Raja2, S. Geethanjali2, P. Raja3, S. Karthikeyan4,*
  • ORCID ID: 0000-0001-7416-9982
1Department of Horticulture, Tamil Nadu Agricultural University, Coimbatore-641 003, Tamil Nadu, India.
2Department of Plant Breeding and Genetics, Tamil Nadu Agricultural University, Coimbatore-641 003, Tamil Nadu, India.
3Department of Microbiology, Tamil Nadu Agricultural University, Coimbatore-641 003, Tamil Nadu, India.
4Department of Floriculture and Landscape Architecture, Horticultural College and Research Institute, Tamil Nadu Agricultural University, Coimbatore-641 003, Tamil Nadu, India.
  • Submitted09-05-2024|

  • Accepted21-06-2024|

  • First Online 22-06-2024|

  • doi 10.18805/LR-5347

Cite article:- Thamaraiselvi P. S., Raja Ajay A., Geethanjali S., Raja P., Karthikeyan S. (2024). Diversity Analysis for Phenotypic and Qualitative Traits of Broad Bean (Vicia faba L.): an Underutilized Vegetable Crop of India by Multivariate Analysis . Legume Research. 47(8): 1274-1281. doi: 10.18805/LR-5347.

Background: Broad bean (Vicia faba L.) is an important protein rich underutilized vegetable crop cultivated for its tender green pods and dry seeds. The present study was carried out to assess the genetic divergence among the broad bean genotypes collected for developing high-yielding cultivars.

Methods: The experimental material consisted of 20 broad bean genotypes, raised in randomized block design with three replications at Horticultural Research Station, Nilgiris, Tamil Nadu, India during Kharif 2023. Data concerning 12 quantitative traits and 7 qualitative traits underwent examination using Mahalanobis D2 statistic, principal component analysis and stepwise multiple regression analysis.

Result: The cluster analysis categorized the genotypes into nine clusters, with the highest number of genotypes in Cluster II (4) and the lowest in Clusters VII, VIII and IX (1 each). Cluster IV displayed the highest intra-cluster distance (D2=3358.2), signifying substantial variation among the genotypes within this cluster. The largest inter-cluster distance (D2=316.9) was observed between Clusters I and VII, indicating significant divergence from other clusters, suggesting that genotypes within these clusters could be valuable as diverse parents for hybridization and selection. The PCA combined with stepwise multiple regression analysis highlighted the primary variations and traits influencing pod yield. 

The broad bean (Vicia faba L.; Family: Fabaceae) is an important underutilized winter sown legume crop and serves as a valuable source of protein (20-30%) for both human and animal consumption (Qahtan et al., 2021). The cultivation of broad beans traces its roots to the Stone Age, marking it as one of the earliest crops to be domesticated (Duc et al., 2010).It is known by several other names, bakla bean (Hindi name), faba bean, horse bean, winter bean, English bean, field bean and Windsor bean (Singh et al., 2017), whereas in Tamil Nadu it is recognized as “double beans”. They are well-suited to a wide range of climatic conditions and contribute to the overall sustainability of cropping systems. Both green and dry cotyledons of broad bean are consumed as vegetable and pulse crop respectively in the hilly regions of India (Gasim and Link, 2007). In the Nilgiris District of Tamil Nadu, broad beans are used in culinary preparations during important occasions and festival times.
       
Out of over 50 countries that produce broad beans, approximately 90% of the production is centred in the regions of Asia, European Union and Africa (FAO, 2020). In India, despite its agricultural potential for producing high yields and protein, it is regarded as an underutilized crop and cultivation in India is confined to limited acreage in Bihar, Madhya Pradesh and some part of Uttar Pradesh (Singh et al., 2012). Lack of a good cultivar along with proper package of practises may be hindering the expansion of this crop (Raja et al., 2023). By adopting improved, nutritionally rich and stable genotypes in cultivation, it is possible to transform the status of this crop from being underutilized to a potentially high-value crop.Therefore, there is a need to evaluate and identify the genotype that consistently delivers both high yield and quality.
       
Harnessing the genetic diversity within the gene pool is of utmost significance for developing high-yielding broad bean cultivars. Multivariate methods such as D2 analysis, principal component analysis and multiple regression analysiscan serve as valuable tools for assessing genetic diversity. Hence, the current research was planned to evaluate various broad bean genotypes based on specific traits identified through multivariate techniques and effective selection criteria to identify superior genotypes among them, which can be used as potential parents in future breeding programs.
The experimental material consisted of twenty broad bean genotypes collected across different geographical locations of The Nilgiris district of Tamil Nadu ranging from 918 to 2192 m above MSL were collected. The crop was raised in randomized block design (RBD) with three replications at the experimental farm (Woodhouse farmat an elevation of 2535 MSL; latitudes of 11.4°N and 76.4°E) of Horticultural Research Station, TNAU, Udhagamandalam,  Nilgiris, Tamil Nadu, India during Kharif 2023. The experimental plot consisted of three rows each separated by 4 m. The spacing between rows and individual plants was maintained at 45 cm and 15 cm, respectively.
       
The observations were recorded on five randomly selected plants in each replications for 12 quantitative and 7 qualitative traits viz., plant height (cm), number of branches per plant, days to 50% flowering, days to maturity, pod length (cm), pod width (cm), number of pods per cluster, number of pods per plant, number of seeds per pod, 100 seed weight (g), seed yield per plant (g), pod yield per plant (g), moisture (%), dry matter (%), ascorbic acid (mg/100 g), protein content (%), total phenol (%), carbohydrate (%) and total soluble solids (°Brix) to study the genetic divergence. Multivariate analysis using Mahalanobis D2 statistic (Mahalanobis, 1936) based on minimum genetic distance using Tocher2 s method as described by Rao (1952) was performed using TNAUSTAT statistical software. The software R version 3.4.3 was utilized to conduct Principal Component Analysis by Kaiser, 1958 and Jeffers, 1967.Stepwise multiple regression analysis described by Lush (1940) was carried out using SPSS software 2.0.
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.
 

Table 1: Variations in morphological and qualitative traits of broad bean genotypes.


 
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 D2 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.
 

Fig 1: Per cent contribution of genetic diversity in broad bean genotypes.


       
Using Tocher’s method, 20 genotypes were categorized into nine clusters based on the relative magnitude of their D2 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).
 

Table 2: Mahalanobis D2 clustering pattern of broad bean genotypes.


       
The intra and inter-cluster D2 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 (D2=2899.7), cluster III (D2=2359.3) and cluster VI (D2=2298.6), whereas the least intra-cluster distance was observed in cluster I (D2=2258.9) and cluster V (D2=1805.8). Cluster VII, VIII and IX had no intra-cluster distance(D2=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 D2 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.
 

Table 3: Average intra (diagonal) and inter cluster (above diagonal) distance (D2) values of clusters in broad bean genotypes.


       
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 D2 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.
 

Table 4: Cluster mean values for various traits among 20 broad bean genotypes.


       
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.
 

Table 5: Superior and diverse broad bean genotypes selected from different clusters for various phenotypic and qualitative traits.


 
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.
 

Fig 2: Scree plot showing contribution of various principal components towards divergence.


       
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).
 

Fig 3: Biplot illustrating the relationship between PC1 and PC2.


       
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).
 

Table 6: Summary of stepwise multiple regression analysis for pod yield and its component traits.

The broad bean genotypes Vf 2, Vf 8, Vf 1, Vf 9, Vf 4 and Vf 19 could be utilized in hybridization programmes to enhance the heterotic response for increased yields. The findings from this study offer essential guidance for the selection of diverse broad bean genotypes and will serve as basis for research in hybridization programmes as well as contributing to germplasm conservation.
The authors acknowledge the support rendered by Tamil Nadu Agricultural University, Coimbatore in the taking up the experiments.
All authors declared that there is no conflict of interest.

  1. Abbas, E. and El-Rodany, (2009). Estimation of genetic variability in faba bean (Vicia faba L.) genotypes under new reclaimed soil in north delta. Journal of Plant Production. 34(10): 10057-10066.

  2. Ammar, M.H., Anwar, F., El Harty, E.H., Migdadi, H.M., Abdel Khalik, S.M., Al Faifi, S., Muhammad Farooq, A. and Alghamdi, S.S. (2015). Assessment of genetic diversity among faba bean genotypes using agro-morphological and molecular markers. Journal of Agronomy and Crop Science. 22(3): 340-350.

  3. Arab, S.A., Azza, F., El-Sayed and Mohamed, M.K.A. (2018). Genetic diversity in some faba bean landraces using morphological characters and yield components. Journal of Plant Production. 9(12): 975-980.

  4. Boghara, M.C., Dhaduk, H.L., Kumar, S., Parekh, M.J., Patel, N.J. and Sharma, R. (2016). Genetic divergence, path analysis and molecular diversity analysis in cluster bean (Cyamopsis tetragonoloba L. Taub.). Industrial Crops and Products. 89: 468-477.

  5. Dahiya, G.S., Kumar, R., Arya, R.K. and Kumari, P. (2021). Evaluation of faba bean elite genotypes for seed yield and quality parameters under Haryana conditions. Forage Res. 47: 271-276.

  6. Dewangan, N.K., Dahiya, G.S., Janghel, D.K. and Dohare, S. (2022). Diversity analysis for seed yield and its component traits among faba bean (Vicia faba L.) germplasm lines. Legume Research: An International Journal. 45(6): 689-694. doi: 10.18805/LR-4301.

  7. Dubey, P., Shweta and Mishra, A. (2022). Assessment of genetic diversity in exotic collections of fava bean (Vicia faba L.). The Pharma Innovation Journal. 11(7): 1445-1450.

  8. Duc, G., Bao, S., Baum, M., Redden, B., Sadiki, M., Suso, M.J., Vishniakova, M. and Zong, X. (2010). Diversity maintenance and use of Vicia faba L. genetic resources. Field Crops Research. 115(3): 270-278.

  9. Ertiro, Abo, T., Kebede, G.Y., Assen, K.Y., Haile, G.A. and Gutu, D.T. (2023). Variability and association of some morpho- agronomic traits in advanced faba bean (Vicia faba L.) genotypes at potential areas of South Eastern Ethiopia. Asian J. Res. Crop Sci. 8(4): 198-210.

  10. FAO (Food and Agriculture Organization), (2020). Crop Production and Trade Data. http://www.fao.org/faostat/en/#data (accessed December 30, 2020).

  11. Gasim, S. and Link, W. (2007). Agronomic performance and the effect of self-fertilization on german winter faba beans. Journal of Central European Agriculture. 8(1): 121-128.

  12. Jadhav, R.A., Mehtre, S.P., Patil, D.K. and Gite, V.K. (2023). Multivariate analysis using D2 and principal component analysis in mung bean [Vigna radiata (L.) Wilczek] for study of genetic diversity. Legume Research-An International Journal. 46(1): 10-17. doi: 10.18805/LR-4508

  13. Janghel, D.K., Kumar, K., Verma, S.S. and Chhabra, A.K. (2020). Genetic relationships and principal component analysis in elite chickpea (Cicer arietinum L.) genotypes for seed yield and its component traits. Legume Research-An International Journal. 43(6): 770-778. doi: 10.18805/LR- 4036.

  14. Jeffers, J.N.R. (1967). Two case studies in the application of principal component analysis. Applied Statistics. 16(3): 225-236.

  15. Kaiser, H.F. (1958). The varimax criterion for analytical rotation in factor analysis. Psychometrika. 23: 187-200.

  16. Lush, J.L. (1940). Intra-sire correlation and regression of offspring on dams as a method of estimating heritability of characters. Journal of Animal Science. 1: 293-301.

  17. Mahalanobis, P.C. (1936). On the generalized distance in statistics. Proceedings of the Indian National Science Academy. 2: 49-55.

  18. Malek, N., Aci, M.M., Khamassi, K., Lupini, A., Rouissi, M. and Hanifi-Mekliche, L. (2021). Agro-morphological and molecular variability among algerian faba bean (Vicia faba L.) accessions. Agronomy. 11(8): 1456. https://doi.org/10.3390/agronomy11081456.

  19. Price, A.L., Patterson, N.J., Plenge, R.M., Weinblatt, M.E., Shadick, N.A. and Reich, D. (2006). Principal components analysis corrects for stratification in genome-wide association studies. Nature Genetics. 38(8): 904-909.

  20. Qahtan, A.A., Al-Atar, A., Abdel-Salam, E.M., El-Sheikh, M.A., Gaafar, A.R.Z. and Faisal, M. (2021). Genetic diversity and structure analysis of a worldwide collection of faba bean (Vicia faba) genotypes using ISSR markers. International Journal of Agriculture and Biology. 25: 683- 691.

  21. Raja, A.A., Karthikeyan, S., Geethanjali, S., Raja, P. and Thamaraiselvi, S.P. (2023). Morphological evaluation of the local genotypes of broad bean (Vicia faba L.) for yield attributes under the nilgiris condition. International Journal of Environment and Climate Change. 13(10): 1916-1924.

  22. Rao, C.R. (1952). Advanced Statistical Methods in Biometrical Research. John Wiley and sons. New York. 236-272.

  23. Rebaa, F., Abid, G., Aouida, M., Abdelkarim, S., Aroua, I., Muhovski, Y., Baudoin, J.P., Mhamdi, M., Sassi, K. and Jebara, M. (2017). Genetic variability in tunisian populations of fababean (Vicia faba L.) assessed by morphological and SSR markers. Physiology and Molecular Biology of Plants. 23(2): 397-409.

  24. Singh, A.K., Bhatt, B.P., Upadhyaya, A., Kumar, S., Sundaram, P.K., Singh, B.K. and Bharati, R.C. (2012). Improvement of faba bean (Vicia faba L.) yield and quality through biotechnological approach: A review. African Journal of Biotechnology. 11(87): 15264-15271.

  25. Singh, D.V., Nath, S., Singh, S.P., Mishra, U. and Singh, S. (2021). Genetic diversity association of different quantitative traits in Faba bean (Vicia faba L.). The Pharma Innovation Journal. 10(4): 514-516.

  26. Singh, S.K., Gautam, S.C., Yadav, C.B. and Nivas, R. (2015). Studies on association of yield and quality contributing parameters in faba bean (Vicia faba L.). Journal of AgriSearch. 2(4): 257-262.

  27. Singh, Y., Sharma, S., Sekhon, B.S., Verma, A. and Sharma, S. (2017). Performance of broad bean (Vicia faba L.) genotypes for pod yield and yield contributing morpho- physiological components under mid hill conditions of North Western Himalayas. Biotechnology Journal International. 19(3): 1-10.

  28. Tiwari, J.K. and Singh, A.K. (2019). Principal component analysis for yield and yield traits in faba bean (Vicia faba L.). Journal of Food Legumes. 32(1): 13-15.

  29. Zayed, E.M., Zeinab, E.G. and Saad, K.I. (2022). Genetic diversity and principle component analysis (PCA) of faba bean landraces based on yield-traits and protein SDS-PAGE. Journal of Global Agriculture and Ecology. 13(4): 1-16.

Editorial Board

View all (0)