Harnessing Genetic Variability, Correlation and Path Analysis among Garden Pea (Pisum sativum L.) Genotypes for Food Security and Sustainability

1Department of Horticulture, Lovely Professional University, Phagwara-144 411, Punjab, India.
2Faculty of Agriculture, Guru Kashi University, Talwandi Sabo, Bathinda, Punjab-151302, India.
3Department of Horticulture (Vegetable Science), Bihar Agricultural University, Sabour-813 210, Bhagalpur, Bihar, India.
4Department of Agronomy, Lovely Professional University, Phagwara-144 411, Punjab, India.
5Department of Genetics and Plant Breeding, Lovely Professional University, Phagwara-144 411, Punjab, India.

Background: A field experiment was conducted at Lovely Professional University, Punjab (2022-23) to evaluate fourteen garden pea (Pisum sativum L.) genotypes for genetic variability, trait associations and direct and indirect effects of yield-related traits through path coefficient analysis. Understanding genetic variation and interrelationships among yield-contributing characters is essential for improving selection efficiency and enhancing productivity in garden pea breeding programs.

Methods: Genotypes were evaluated in a randomized complete block design (RCBD) with three replications. Observations were recorded on morphological and yield-related parameters including plant height, number of nodes to first flower, number of primary branches/plants, pod length, number of pods/plants, 100-seed weight and pod yield/plant. Data were statistically analysed for analysis of variance (ANOVA), phenotypic coefficient of variation (PCV), genotypic coefficient of variation (GCV), heritability, genotypic and phenotypic correlations and path coefficient analysis to determine direct and indirect effects of various traits on pod yield.

Result: Analysis of variance revealed significant genetic variability among genotypes for all studied traits. PCV values exceeded GCV for most traits, indicating environmental influence, except for nodes/plant. High GCV, PCV and heritability were observed for traits like number of primary branches/plants, pod yield/plant, 100-seed weight and numbers of pods/plant. Pod yield showed strong positive correlations with number of pods, 100-seed weight and node to first pod appearance. Path analysis indicated that number of primary branches, pod length, plant height and nodes to first flower appearance had the strongest direct effects on yield, suggesting these as key selection criteria in breeding programs.

Garden pea (Pisum sativum L. sub sp. hortense Asch.; 2n=2x=14) is a member of Leguminosae family and its primary center of origin are Central Asia, Abyssinia and Mediterranean region, while Near East is considered its secondary origin (Zohary and Hopf, 1973: Srivastava, 2025). Peas are cultivated for their immature pods, dry seeds, canned, frozen or dehydrated form (Singh et al., 2022). It has a high nutritional content, particularly proteins and several other health-building components such as carbohydrates, minerals, vitamins and amino acids, like lysine (Ghobary, 2010; Singh et al., 2019). Inclusion of pea in crop rotation is important agronomically due to its high nitrogen fixation efficiency and better potential to utilize minerals (Adgo and Schulze, 2002). A crop performance is the outcome of interactions between the crop’s genetic makeup and environment (Bharti et al., 2024: Neupane et al., 2023). Germplasm screening is vital for identifying promising genotypes, while understanding genetic variability enables the discovery of key genes for crop improvement (Azmat et al., 2011; Nwangburuka et al., 2011). Path coefficient analysis further clarifies trait associations by partitioning correlation coefficients into direct and indirect effects, thereby highlighting the relative importance of each factor (Bhatt, 1973; Patil et al., 2024). Such insights are essential for developing effective breeding programs.  Thus, this study was designed to assess the extent of genetic variability and the relationships between traits that contribute to green pod yield/plant in peas.
Experimental site
 
This particular experimental trial was conducted at the Research Farm of Lovely Professional University, Punjab, India, during Rabi season-2022-2023.
 
Experimental material and design
 
There were fourteen distinct genotypes of garden pea selected for this study including SATAYA-1010 (Kulak Seeds Company), Peas Vasundhara (Tycoon Seeds), Zanna + (Aggarwal Seeds), Imported Peas OS-10 (Omaxe Seeds), Samag (Samag Seeds), Sweet Ruby (Shriram Seeds), AP-3 (CSAUA and T, Kanpur), GS-10 (Advanta Golden Seeds), Niarli (Agro Seeds), Peas A-1 (Hara Seeds), Tako-10 (Takoyama  Seeds), 11 Danno Wali (tiger) (Kalash Seed Pvt. Ltd.), BKS peas (Bundelkhand Seeds) and PB-89 (PAU, Ludhiana, Punjab used as a check variety) procured from PAU, Ludhiana, Punjab and other private companies. The selected genotypes were arranged in a randomized block design (RBD) with three replications.
 
Seed sowing
 
Seed sowing was done at 30 cm × 10 cm spacing and all essential cultural-practices were followed as suggested by PAU, Ludhiana, Punjab for raising a healthy crop.
 
Analysis of variance (ANOVA)
 
The collected data for the various studied parameters were analyzed using analysis of variance (ANOVA) to assess significance levels at P = 0.01 and 0.05 (Panse and Sukhatme, 1961).

Variability, heritability, path and correlation analysis
 
Variability coefficients, genetic advance (GA) and heritability in broad-sense (h²b) were analyzed by following the methods of Johnson et al., (1955) and Burton and De Vane (1953). Subsequently, Correlation and path analyses at both phenotypic and genotypic levels were carried out as described by Wright (1934), Al-Jibouri et al. (1958) and Dewey and Lu (1959).
Analysis of variance (ANOVA)
 
The analysis of variance (ANOVA) results demonstrated a significant variation among genotypes for the considered traits in our study.
 
Variability analysis
 
The magnitude of PCV was greater than GCV for most parameters, with the exception of nodes/plant, which demonstrated an additive gene effect driven by the environment (Table 1). Characters like number of pods/plants, 100-seed weight, number of primary branches/plant and pod yield/plant exhibited the highest GCV and PCV values, both reaching 20%. These results suggest that certain traits exhibit higher variation genetically, that might be beneficial for crop improvement under selective breeding. However, environmental factors must also be considered when selecting traits with high variability. These results were in conformity with of Pal and Singh (2012); Katoch et al., (2016); Pandey et al., (2017); Kumar et al., (2019) and Gupta et al., (2020).

Table 1: Estimation of genetic variability and its components for various traits among garden pea genotypes.


       
Moderate GCV and PCV values (ranging 10% to 20%) were recorded for traits like number of seeds/100 g of pods, number of seeds/pods, plant height, number of pods/100 g and TSS. These findings are consistent with those of Thakur et al., (2016) and Barcchiya et al., (2018). The lowest values (less than 10%) of PCV and GCV were noticed for various characters including days to 50% flowering, node to 1st pod appears, internodal length (cm), node to 1st flower appears and nodes/plant. Similar results were presented by Katoch et al., (2016); Kumar et al., (2019) and Jaiswal et al., (2013).
 
Heritability in broad sense (%)
 
The heritability values were found to range from 68.90% for the trait node to 1st pod appear to a remarkably high 99.25% for the trait pod yield/plant (Table 1). Among traits, the highest heritability was observed in pod yield/plant (99.25%), followed closely by plant height (97.65%), number of seeds/100 g pods (96.72%), total soluble solids (96.18%), 100-seed weight (95.94%), number of pods/plant (95.25%), number of primary branches/plant (93.46%), number of pods/100 g (92.00%), number of seeds/pod (89.76%), number of nodes/plant (89.31%), pod length (88.77%), days to 50% flowering (86.33%), internodal length (71.77%), node to 1st flower appear (69.78%) and node to 1st pod appears (68.90%). These findings suggest that there is substantial genetic control over these traits in the plant species under study. These results were in confirmation with the findings of Pal and Singh (2012); Katoch et al., (2016); Kumar et al., (2019); Gupta et al., (2020); Thakur et al., (2016); Barcchiya et al., (2018); Jaiswal et al., (2013); Sharma and Bora, (2013); Georgieva et al., (2016); Gudadinni et al., (2017); Kumar et al., (2018); Singh et al., (2019) and Kumar et al., (2023).
 
Genetic advance (GA) as percentage of mean
 
Results presented in Table 1 revealed that higher heritability infused with high GA in percent of mean was observed for 100-seed weight, followed by number of pods/plants, number of primary branches/plants, number of seeds/100 g pods, number of pods/100 g, pod yield/plant and plant height. Lower genetic advance along with high heritability as % of mean was observed for several characters viz., internodal length, days to 50% flowering, node to 1st pod appears, pod length, node to 1st flower appears and numbers of nodes/plant. These outcomes align with the conclusions of Pandey et al., (2017); Kumar et al., (2019); Thakur et al., 2016); Georgieva et al., (2016); Jagadeesh et al., (2023) and Aziz-ur-Rahman et al. (2021).
 
Correlation analysis
 
The results of correlation analysis among traits revealed a significant positive relationship between pod yield/plant (Table 2). Various traits such as node to 1st pod appearance, number of pods/plants and 100-seed weight may be considered as key influencers of pod yield in peas. Analogous outcome was recorded formerly by Mohanty et al., (2020) and Thapa et al., (2020). The trait plant height showed a high negative significant value with total soluble solids (TSS), while it was negatively significant with node to 1st pod appears and number of primary branches/plants. Based on the results, traits such as number of pods/plants, 100-seed weight and nodes to 1st pod appearance should be prioritized when selecting high-yielding genotypes for garden peas, as they show a strong correlation with pod yield.

Table 2: Genotypic and phenotypic correlation coefficient for different traits among garden pea genotypes.


 
Path analysis
 
Path analysis at genotypic and phenotypic both levels for various traits among genotypes is presented in Table 3. At genotypic level, the strongest direct positive effect on pod yield/plant was observed for the number of pods/plant (0.8576), followed by pod length (0.4974) and the number of pods/100 g (0.2559). Similar kind of findings were reported by Pal and Singh (2012) and Khan et al., (2017). The positive correlation of these traits with pod yield/plant can be attributed to the positive indirect effects of other related traits (Fig 1). Whereas, negative direct effect on pod yield/plant were found for traits like plant height (-0.4177), number of seeds/pod (-0.3873), number of seeds/100 g pods (-0.3159), days to 50% flowering (-0.1656), number of nodes/plant (-0.0307) and TSS (-0.4070).

Table 3: Direct and indirect effect of different studied traits on pod yield/plant (g) at phenotypic and genotypic level among garden pea genotypes.



Fig 1: Genotypic path analysis illustrating the direct and indirect effects of yield-contributing traits on pod yield per plant in garden pea genotypes.


       
Whereas, the direct highest positive effect on pod yield/plant was refereed number of primary branches/plant (4.4041), followed by node to 1st flower appears (2.1220), pod length (1.997 cm) and plant height (1.968 cm) at the phenotypic level (Fig 2). Similarly, findings were recorded by Saxena et al., (2014); Gupta et al. (2020) and Mehandi and Mishra, (2023) for pod length and plant height and by Mohanty et al., (2020) for plant height. Hence, these traits should be considered in the future selection procedures for getting higher pod yield/plant (Khan et al., 2017). The negative direct effect on pod yield/plant were found by numerous traits such as number of nodes/plant (-1.410), internodal length (-1.093), number of pods/plant (-1.670), days to 50% flowering (-1.323), number of seeds/pod (-1.344) and number of seeds/100 g pods (-1.339) at phenotypic level. These outcomes align with the findings of Bijalwan et al. (2018). The maximum positive indirect effect was reported by number of primary branches/plants on characters like number of pods/plant (2.800), number of nodes/plant (2.422), number of seeds/pod (2.138) and node to 1st pod appears (2.030). In the present research, the residual effect was calculated to be 0.08, indicating that roughly 99.92% variation in pod yield may explained by the investigated traits. Nevertheless, it’s evident that other factors, not accounted for in this analysis, must be included to comprehensively address the remaining variability in yield.

Fig 2: Phenotypic path analysis illustrating the direct and indirect effects of yield-contributing traits on pod yield per plant in garden pea genotypes.

The character like number of primary branches/plants was found highly genetically variable and exhibits high heritability, indicating its potential for selective breeding to enhance yield and productivity. Similarly, the number of pods/plants showed the high genetic variability, moderate heritability and associated with high genetic advance (GA), making it a key determinant of pod yield/plant. The trait pod yield/plant had higher GA and heritability thus, it emerges as a critical target for breeding programs aimed at improving overall yield and productivity. Correlation analysis revealed a strong positive relationship among pod yield/plant and variables such as number of pods/plants, 100-seed-weight and node to 1st pod appears. These findings suggested that these characteristics may play an important role in determining pod output in garden peas. Path analysis revealed that at the genotypic level, the trait number of pods/plants had the greatest direct beneficial effect on pod yield/plant, followed by pod length and number of pods/100 g. Similarly, at the phenotypic level, the number of primary branches/plants had the maximum favourable direct effect on pod yield/plant, followed by plant height, node to I flower appears and pod length. Thus, by identifying and harnessing genetic traits that enhance yield, resilience and nutritional value, we can contribute to the sustainable developmental goals (SDGs) by improving crop productivity, reducing environmental impact and ensuring a reliable food supply for growing populations. The obtained findings not only advance our scientific understanding but also support global efforts towards sustainable development and the well-being of communities worldwide.
Authors are thankful to Lovely Professional University, Jalandhar, Punjab, India for providing the infrastructural support.
 
Authors’ contribution
 
Vishal Tripathi designed the technical programme, edited and reviewed the work. Arkit Neupane conducted the research and collected the data. Suhel Mehandi did the statistical analysis. Yaman Kumar, Etalesh Goutam, Amrita Kumari, Rupinder Singh and Rajeev prepared the manuscript. All the authors approved the final version of the manuscript.
 
Disclaimers
 
The views and conclusions expressed in this article are solely those of the authors and do not necessarily represent the views of their affiliated institutions. The authors are responsible for the accuracy and completeness of the information provided, but do not accept any liability for any direct or indirect losses resulting from the use of this content.
 
Informed consent
 
All animal procedures for experiments were approved by the Committee of Experimental Animal care and handling techniques were approved by the University of Animal Care Committee.
The authors declare that there are no conflicts of interest regarding the publication of this article. No funding or sponsorship influenced the design of the study, data collection, analysis, decision to publish, or preparation of the manuscript. 

  1. Adgo, E. and Schulze, J. (2002). Nitrogen fixation and assimilation efficiency in Ethiopian and German pea varieties. Plant and Soil. 239: 291-299. doi: 10.1023/A:1015048331366.

  2. Al-Jibouri, H.A., Miller, P.A. and Robinson, H.F. (1958). Genotypic and environmental variances and covariances in an upland cotton cross of interspecific origin. Agronomy Journal. 50: 633-636.

  3. Aziz-ur-Rahman, Katoch, V., Rathour, R., Sharma, S., Rana, S.S. and Sharma, A. (2021). Studies on genetic parameters, correlation and path coefficient analysis in er2 introgressed garden pea genotypes. Legume Research. 44(6): 621-626. doi: 10.18805/LR-4142.

  4. Azmat, M.A., Nawab, N.N., Khan, A.A., Ashraf, M., Niaz, S. and Mahmood, K. (2011). Characterization of pea germplasm. International Journal of Vegetable Science. 17: 246- 258. doi: 10.1080/19315260.2010.544380.

  5. Barcchiya, J., Naidu, A.K., Mehta, A.K. and Upadhyay, A. (2018). Genetic variability, heritability and genetic advance for yield and yield components in pea (Pisum sativum L.). International Journal of Chemical Studies. 6: 3324-3327.

  6. Bharti, Mal, D., Goutam, E. and Kumar, L. (2024). Performance of french bean (Phaseolus vulgaris L.) genotypes under trans-gangetic plains region. Indian Journal of Agricultural Research. 58(1): 101-105. doi: 10.18805/IJARe.A-6148.

  7. Bhatt, G.M. (1973). Significance of path coefficient analysis in determining the nature of character association. Euphytica. 22: 338-343. doi: 10.1007/BF00022643.

  8. Bijalwan, P., Raturi, A. and Mishra, A.C. (2018). Character association and path analysis studies in garden pea (Pisum sativum L.) for yield and yield attributes. International Journal of Current Microbiology and Applied Sciences. 7: 3491- 3495. doi: 10.20546/ijcmas.2018.703.401.

  9. Burton, C.W. and De Vane, E.H. (1953). Estimating heritability in tall fescue (Festuca arundinacae) from replicated clonal material. Agronomy Journal. 45: 478-481.

  10. Dewey, D.R. and Lu, K. (1959). A correlation and path-coefficient analysis of components of crested wheatgrass seed production. Agronomy Journal. 51: 515-518.

  11. Georgieva, N., Nikolova, I. and Kosev, V. (2016). Evaluation of genetic divergence and heritability in pea (Pisum sativum L.). Journal of BioScience and Biotechnology. 5: 61-67.

  12. Ghobary, H.M.M. (2010). Study of relationship between yield and some yield components in garden pea (Pisum sativum L.) by using correlation and path analysis. Journal of Agricultural Research. 36: 351-360.

  13. Gudadinni, P., Bahadur, V., Ligade, P., Topno, S.E. and Prasad, V.M. (2017). Study on genetic variability, heritability and genetic advance in garden pea (Pisum sativum var. hortense L.). International Journal of Current Microbiology and Applied Sciences. 6: 2384-2391. doi: 10.20546/ ijcmas.2017.603.282.

  14. Gupta, A., Singh, B., Kumar, M., Chand, P., Vaishali and Sharma, V.R. (2020a). Studies on genetic variability, heritability and genetic advance in table pea (Pisum sativum var. hortense L.). International Journal of Current Microbiology and Applied Sciences. 9: 3449-3454. doi: 10.20546/ ijcmas.2020.909.428.

  15. Gupta, A., Singh, B., Prakash, S., Rout, S. and Singh, J. (2020b). Correlation coefficient and path analysis in table pea (Pisum sativum var. hortense L.). Current Journal of Applied Science and Technology. 39: 56-60. doi: 10.9734/ CJAST/2020/v39i3130989.

  16. Jagadeesh, K., Mahto, C.S. and Kumar, N. (2023). Genetic variability studies in field pea (Pisum sativum L.) for yield and associated characters. Environment Conservation Journal. 24: 244-249. doi: 10.36953/ECJ.13172379.

  17. Jaiswal, N.K., Gupta, A.K., Dewangan, H. and Lavanya, G.R. (2013). Genetic variability analysis in field pea (Pisum sativum L.). International Journal of Scientific Research. 4: 2006- 2007.

  18. Johnson, H.W., Robinson, H.F. and Comstock, R.E. (1955). Estimates of genetic and environmental variability in soyabean. Agronomy Journal. 47: 314-318.

  19. Katoch, V., Singh, P., Devi, M.B., Sharma, A., Sharma, G.D. and Sharma, J.K. (2016). Study of genetic variability, character association, path analysis and selection parameters for heterotic recombinant inbred lines of garden peas (Pisum sativum var. hortense L.) under mid-hill conditions of Himachal Pradesh, India. Legume Research. 39(2): 163- 169. doi: 10.18805/lr.v0iOF.6775.

  20. Khan, M.R.A., Mahmud, F., Reza, M.A., Mahbub, M.M., Shirazy, B.J. and Rahman, M.M. (2017). Genetic diversity, correlation and path analysis for yield and yield components of pea (Pisum sativum L.). World Journal of Agricultural Sciences. 13: 11-16. doi: 10.5829/idosi.wjas.2017.11.16.

  21. Kumar, R., Singh, B., Tomar, H., Kumar, V., Chaudhary, A.K., Marwah, A. and Yashraj. (2023). Genetic heritability, variability, genetic advance, correlation and path analysis assessment in garden pea (Pisum sativum L.). International Journal of Plant and Soil Science. 35: 69-79. doi: 10.9734/ijpss/ 2023/v35i143022.

  22. Kumar, T.N.V., Alloli, T.B., Hadimani, H.P., Ajjappalavar, P.S., Satish, D., Abdul, K. and Hanchinamani, C.N. (2019). Studies on genetic variability, heritability and genetic advance in garden pea (Pisum sativum L.) varieties. International Journal of Current Microbiology and Applied Sciences. 8: 3032-3038. doi: 10.20546/ijcmas.2019.812.353.

  23. Kumar, V., Singh, J., Srivastava, C.P. and Lal, C. (2018). Genetic variability, heritability and genetic advance for seed yield and related traits in diverse genotypes of pea (Pisum sativum L.). International Journal of Chemical Studies. 6: 1174-1178.

  24. Mehandi, S. and Mishra, S.P. (2023). Genetic diversity and population structure analyses in mungbean (Vigna radiata L. Wilczek). Journal of Food Chemistry and Nanotechnology. 9(S1): S316-S322. doi: 10.17756/jfcn.2023-s1-041.

  25. Mohanty, A.T., Verma, A. and Pandey, B. (2020). Evaluation of correlation and path coefficients analysis for yield attributing traits in garden pea (Pisum sativum L.) under Tarai tract of Uttarakhand. International Journal of Current Microbiology and Applied Sciences. S11: 1315- 1322.

  26. Neupane, A., Tripathi, V., Swain, A.R., Sharma, J.K., Kumar, Y. and Goutam, E. (2023). Varietal mean performance of garden pea (Pisum sativum var. hortense L.) under subtropical conditions of Punjab, India. International Journal of Environment and Climate Change. 13: 1746-1755. doi: 10.9734/IJECC/2023/v13i113558.

  27. Nwangburuka, C.C., Kehinde, O.B., Ojo, D.K., Denton, O.A. and Popoola, A.R. (2011). Morphological classification of genetic diversity in cultivated okra [Abelmoschus esculentus (L.) Moench] using PCA and SLCA. African Journal of Biotechnology. 10: 11165-11172. doi: 10.5897/ AJB11.285.

  28. Pal, A.K. and Singh, S. (2012). Correlation and path analysis in garden pea (Pisum sativum L. var. hortense). Asian Journal of Horticulture. 7: 569-573.

  29. Pandey, M., Singh, V.B., Yadav, G.C., Tyagi, N., Vishen, G.S., Bhargav, K.K. and Pandey, P. (2017). Correlation and path coefficient analysis among different characters in genotypes of vegetable pea. Vegetable Science. 44: 139-141.

  30. Panse, V.G. and Sukhatme, P.V. (1961). Statistical Methods for Agricultural Workers. ICAR, New Delhi. pp. 1-328.

  31. Patil, P., Shrivastav, S.P., Patil, K., Landge, R. and Gurjar, D. (2024). Genetic variability, heritability, genetic advance and divergence analysis in wheat (Triticum aestivum L.). Indian Journal of Agricultural Research. 58(2): 209-214. doi: 10.18805/IJARe.A-6036.

  32. Saxena, R.R., Vidayakar, V., Vishkarma, M.V., Yadav, P.S., Meena, M.L. and Lal, G.M. (2014). Genetic variability and heritability analysis for some quantitative traits in field pea (Pisum sativum L.). Bioscan. 9: 895-898.

  33. Sharma, V.K. and Bora, L. (2013). Studies on genetic variability and heterosis in vegetable pea (Pisum sativum L.) under high hills condition of Uttarakhand, India. African Journal of Agricultural Research. 8: 1891-1895.

  34. Singh, S., Singh, B., Sharma, V.R., Kumar, M. and Sirohi, U. (2022). Assessment of genetic diversity and population structure in pea (Pisum sativum L.) germplasm based on morphological traits and SSR markers. Legume Research. 45(6): 683- 688. doi: 10.18805/LR-4751.

  35. Singh, S., Verma, V., Singh, B., Sharma, V.R. and Kumar, M. (2019). Genetic variability, heritability and genetic advance studies in pea (Pisum sativum L.) for quantitative characters. Indian Journal of Agricultural Research. 53(5): 542-547. doi: 10.18805/IJARe.A-5245.

  36. Srivastava, A., Nath, S., Singh, P. and Sharma, P. (2025). Assessment of yield-contributing traits and genetic variability in field pea (Pisum sativum L.) under heat stress using half- diallel derived genotypes. Legume Research. 48(10): 1627-1634. doi: 10.18805/LR-5526.

  37. Thakur, S., Thakur, R. and Mehta, K.D. (2016). Genetic variability and association studies for green pod yield and component horticultural traits in garden pea under high hill dry temperate conditions of Tabo valley of Spiti district of Himachal Pradesh. International Journal of Science Environment and Technology. 5: 1987-1992.

  38. Thapa, U., Nandi, S. and Gurung, D. (2020). Performance of garden pea (Pisum sativum var. hortense L.) genotypes and their genetic variation, genetic advance, character association and path analysis for pod, seed and quality characters. International Journal of Current Microbiology and Applied Sciences. 9: 3572-3586. doi: 10.20546/ ijcmas.2020.909.442.

  39. Wright, S. (1934). The method of path coefficients. Annals of Mathematical Statistics. 5: 161-215.

  40. Zohary, D. and Hopf, M. (1973). Domestication of pulses in the old world: Legumes were companions of wheat and barley when agriculture began in the near East. Science. 182: 887-894. doi: 10.1126/science.182.4115.887.

Harnessing Genetic Variability, Correlation and Path Analysis among Garden Pea (Pisum sativum L.) Genotypes for Food Security and Sustainability

1Department of Horticulture, Lovely Professional University, Phagwara-144 411, Punjab, India.
2Faculty of Agriculture, Guru Kashi University, Talwandi Sabo, Bathinda, Punjab-151302, India.
3Department of Horticulture (Vegetable Science), Bihar Agricultural University, Sabour-813 210, Bhagalpur, Bihar, India.
4Department of Agronomy, Lovely Professional University, Phagwara-144 411, Punjab, India.
5Department of Genetics and Plant Breeding, Lovely Professional University, Phagwara-144 411, Punjab, India.

Background: A field experiment was conducted at Lovely Professional University, Punjab (2022-23) to evaluate fourteen garden pea (Pisum sativum L.) genotypes for genetic variability, trait associations and direct and indirect effects of yield-related traits through path coefficient analysis. Understanding genetic variation and interrelationships among yield-contributing characters is essential for improving selection efficiency and enhancing productivity in garden pea breeding programs.

Methods: Genotypes were evaluated in a randomized complete block design (RCBD) with three replications. Observations were recorded on morphological and yield-related parameters including plant height, number of nodes to first flower, number of primary branches/plants, pod length, number of pods/plants, 100-seed weight and pod yield/plant. Data were statistically analysed for analysis of variance (ANOVA), phenotypic coefficient of variation (PCV), genotypic coefficient of variation (GCV), heritability, genotypic and phenotypic correlations and path coefficient analysis to determine direct and indirect effects of various traits on pod yield.

Result: Analysis of variance revealed significant genetic variability among genotypes for all studied traits. PCV values exceeded GCV for most traits, indicating environmental influence, except for nodes/plant. High GCV, PCV and heritability were observed for traits like number of primary branches/plants, pod yield/plant, 100-seed weight and numbers of pods/plant. Pod yield showed strong positive correlations with number of pods, 100-seed weight and node to first pod appearance. Path analysis indicated that number of primary branches, pod length, plant height and nodes to first flower appearance had the strongest direct effects on yield, suggesting these as key selection criteria in breeding programs.

Garden pea (Pisum sativum L. sub sp. hortense Asch.; 2n=2x=14) is a member of Leguminosae family and its primary center of origin are Central Asia, Abyssinia and Mediterranean region, while Near East is considered its secondary origin (Zohary and Hopf, 1973: Srivastava, 2025). Peas are cultivated for their immature pods, dry seeds, canned, frozen or dehydrated form (Singh et al., 2022). It has a high nutritional content, particularly proteins and several other health-building components such as carbohydrates, minerals, vitamins and amino acids, like lysine (Ghobary, 2010; Singh et al., 2019). Inclusion of pea in crop rotation is important agronomically due to its high nitrogen fixation efficiency and better potential to utilize minerals (Adgo and Schulze, 2002). A crop performance is the outcome of interactions between the crop’s genetic makeup and environment (Bharti et al., 2024: Neupane et al., 2023). Germplasm screening is vital for identifying promising genotypes, while understanding genetic variability enables the discovery of key genes for crop improvement (Azmat et al., 2011; Nwangburuka et al., 2011). Path coefficient analysis further clarifies trait associations by partitioning correlation coefficients into direct and indirect effects, thereby highlighting the relative importance of each factor (Bhatt, 1973; Patil et al., 2024). Such insights are essential for developing effective breeding programs.  Thus, this study was designed to assess the extent of genetic variability and the relationships between traits that contribute to green pod yield/plant in peas.
Experimental site
 
This particular experimental trial was conducted at the Research Farm of Lovely Professional University, Punjab, India, during Rabi season-2022-2023.
 
Experimental material and design
 
There were fourteen distinct genotypes of garden pea selected for this study including SATAYA-1010 (Kulak Seeds Company), Peas Vasundhara (Tycoon Seeds), Zanna + (Aggarwal Seeds), Imported Peas OS-10 (Omaxe Seeds), Samag (Samag Seeds), Sweet Ruby (Shriram Seeds), AP-3 (CSAUA and T, Kanpur), GS-10 (Advanta Golden Seeds), Niarli (Agro Seeds), Peas A-1 (Hara Seeds), Tako-10 (Takoyama  Seeds), 11 Danno Wali (tiger) (Kalash Seed Pvt. Ltd.), BKS peas (Bundelkhand Seeds) and PB-89 (PAU, Ludhiana, Punjab used as a check variety) procured from PAU, Ludhiana, Punjab and other private companies. The selected genotypes were arranged in a randomized block design (RBD) with three replications.
 
Seed sowing
 
Seed sowing was done at 30 cm × 10 cm spacing and all essential cultural-practices were followed as suggested by PAU, Ludhiana, Punjab for raising a healthy crop.
 
Analysis of variance (ANOVA)
 
The collected data for the various studied parameters were analyzed using analysis of variance (ANOVA) to assess significance levels at P = 0.01 and 0.05 (Panse and Sukhatme, 1961).

Variability, heritability, path and correlation analysis
 
Variability coefficients, genetic advance (GA) and heritability in broad-sense (h²b) were analyzed by following the methods of Johnson et al., (1955) and Burton and De Vane (1953). Subsequently, Correlation and path analyses at both phenotypic and genotypic levels were carried out as described by Wright (1934), Al-Jibouri et al. (1958) and Dewey and Lu (1959).
Analysis of variance (ANOVA)
 
The analysis of variance (ANOVA) results demonstrated a significant variation among genotypes for the considered traits in our study.
 
Variability analysis
 
The magnitude of PCV was greater than GCV for most parameters, with the exception of nodes/plant, which demonstrated an additive gene effect driven by the environment (Table 1). Characters like number of pods/plants, 100-seed weight, number of primary branches/plant and pod yield/plant exhibited the highest GCV and PCV values, both reaching 20%. These results suggest that certain traits exhibit higher variation genetically, that might be beneficial for crop improvement under selective breeding. However, environmental factors must also be considered when selecting traits with high variability. These results were in conformity with of Pal and Singh (2012); Katoch et al., (2016); Pandey et al., (2017); Kumar et al., (2019) and Gupta et al., (2020).

Table 1: Estimation of genetic variability and its components for various traits among garden pea genotypes.


       
Moderate GCV and PCV values (ranging 10% to 20%) were recorded for traits like number of seeds/100 g of pods, number of seeds/pods, plant height, number of pods/100 g and TSS. These findings are consistent with those of Thakur et al., (2016) and Barcchiya et al., (2018). The lowest values (less than 10%) of PCV and GCV were noticed for various characters including days to 50% flowering, node to 1st pod appears, internodal length (cm), node to 1st flower appears and nodes/plant. Similar results were presented by Katoch et al., (2016); Kumar et al., (2019) and Jaiswal et al., (2013).
 
Heritability in broad sense (%)
 
The heritability values were found to range from 68.90% for the trait node to 1st pod appear to a remarkably high 99.25% for the trait pod yield/plant (Table 1). Among traits, the highest heritability was observed in pod yield/plant (99.25%), followed closely by plant height (97.65%), number of seeds/100 g pods (96.72%), total soluble solids (96.18%), 100-seed weight (95.94%), number of pods/plant (95.25%), number of primary branches/plant (93.46%), number of pods/100 g (92.00%), number of seeds/pod (89.76%), number of nodes/plant (89.31%), pod length (88.77%), days to 50% flowering (86.33%), internodal length (71.77%), node to 1st flower appear (69.78%) and node to 1st pod appears (68.90%). These findings suggest that there is substantial genetic control over these traits in the plant species under study. These results were in confirmation with the findings of Pal and Singh (2012); Katoch et al., (2016); Kumar et al., (2019); Gupta et al., (2020); Thakur et al., (2016); Barcchiya et al., (2018); Jaiswal et al., (2013); Sharma and Bora, (2013); Georgieva et al., (2016); Gudadinni et al., (2017); Kumar et al., (2018); Singh et al., (2019) and Kumar et al., (2023).
 
Genetic advance (GA) as percentage of mean
 
Results presented in Table 1 revealed that higher heritability infused with high GA in percent of mean was observed for 100-seed weight, followed by number of pods/plants, number of primary branches/plants, number of seeds/100 g pods, number of pods/100 g, pod yield/plant and plant height. Lower genetic advance along with high heritability as % of mean was observed for several characters viz., internodal length, days to 50% flowering, node to 1st pod appears, pod length, node to 1st flower appears and numbers of nodes/plant. These outcomes align with the conclusions of Pandey et al., (2017); Kumar et al., (2019); Thakur et al., 2016); Georgieva et al., (2016); Jagadeesh et al., (2023) and Aziz-ur-Rahman et al. (2021).
 
Correlation analysis
 
The results of correlation analysis among traits revealed a significant positive relationship between pod yield/plant (Table 2). Various traits such as node to 1st pod appearance, number of pods/plants and 100-seed weight may be considered as key influencers of pod yield in peas. Analogous outcome was recorded formerly by Mohanty et al., (2020) and Thapa et al., (2020). The trait plant height showed a high negative significant value with total soluble solids (TSS), while it was negatively significant with node to 1st pod appears and number of primary branches/plants. Based on the results, traits such as number of pods/plants, 100-seed weight and nodes to 1st pod appearance should be prioritized when selecting high-yielding genotypes for garden peas, as they show a strong correlation with pod yield.

Table 2: Genotypic and phenotypic correlation coefficient for different traits among garden pea genotypes.


 
Path analysis
 
Path analysis at genotypic and phenotypic both levels for various traits among genotypes is presented in Table 3. At genotypic level, the strongest direct positive effect on pod yield/plant was observed for the number of pods/plant (0.8576), followed by pod length (0.4974) and the number of pods/100 g (0.2559). Similar kind of findings were reported by Pal and Singh (2012) and Khan et al., (2017). The positive correlation of these traits with pod yield/plant can be attributed to the positive indirect effects of other related traits (Fig 1). Whereas, negative direct effect on pod yield/plant were found for traits like plant height (-0.4177), number of seeds/pod (-0.3873), number of seeds/100 g pods (-0.3159), days to 50% flowering (-0.1656), number of nodes/plant (-0.0307) and TSS (-0.4070).

Table 3: Direct and indirect effect of different studied traits on pod yield/plant (g) at phenotypic and genotypic level among garden pea genotypes.



Fig 1: Genotypic path analysis illustrating the direct and indirect effects of yield-contributing traits on pod yield per plant in garden pea genotypes.


       
Whereas, the direct highest positive effect on pod yield/plant was refereed number of primary branches/plant (4.4041), followed by node to 1st flower appears (2.1220), pod length (1.997 cm) and plant height (1.968 cm) at the phenotypic level (Fig 2). Similarly, findings were recorded by Saxena et al., (2014); Gupta et al. (2020) and Mehandi and Mishra, (2023) for pod length and plant height and by Mohanty et al., (2020) for plant height. Hence, these traits should be considered in the future selection procedures for getting higher pod yield/plant (Khan et al., 2017). The negative direct effect on pod yield/plant were found by numerous traits such as number of nodes/plant (-1.410), internodal length (-1.093), number of pods/plant (-1.670), days to 50% flowering (-1.323), number of seeds/pod (-1.344) and number of seeds/100 g pods (-1.339) at phenotypic level. These outcomes align with the findings of Bijalwan et al. (2018). The maximum positive indirect effect was reported by number of primary branches/plants on characters like number of pods/plant (2.800), number of nodes/plant (2.422), number of seeds/pod (2.138) and node to 1st pod appears (2.030). In the present research, the residual effect was calculated to be 0.08, indicating that roughly 99.92% variation in pod yield may explained by the investigated traits. Nevertheless, it’s evident that other factors, not accounted for in this analysis, must be included to comprehensively address the remaining variability in yield.

Fig 2: Phenotypic path analysis illustrating the direct and indirect effects of yield-contributing traits on pod yield per plant in garden pea genotypes.

The character like number of primary branches/plants was found highly genetically variable and exhibits high heritability, indicating its potential for selective breeding to enhance yield and productivity. Similarly, the number of pods/plants showed the high genetic variability, moderate heritability and associated with high genetic advance (GA), making it a key determinant of pod yield/plant. The trait pod yield/plant had higher GA and heritability thus, it emerges as a critical target for breeding programs aimed at improving overall yield and productivity. Correlation analysis revealed a strong positive relationship among pod yield/plant and variables such as number of pods/plants, 100-seed-weight and node to 1st pod appears. These findings suggested that these characteristics may play an important role in determining pod output in garden peas. Path analysis revealed that at the genotypic level, the trait number of pods/plants had the greatest direct beneficial effect on pod yield/plant, followed by pod length and number of pods/100 g. Similarly, at the phenotypic level, the number of primary branches/plants had the maximum favourable direct effect on pod yield/plant, followed by plant height, node to I flower appears and pod length. Thus, by identifying and harnessing genetic traits that enhance yield, resilience and nutritional value, we can contribute to the sustainable developmental goals (SDGs) by improving crop productivity, reducing environmental impact and ensuring a reliable food supply for growing populations. The obtained findings not only advance our scientific understanding but also support global efforts towards sustainable development and the well-being of communities worldwide.
Authors are thankful to Lovely Professional University, Jalandhar, Punjab, India for providing the infrastructural support.
 
Authors’ contribution
 
Vishal Tripathi designed the technical programme, edited and reviewed the work. Arkit Neupane conducted the research and collected the data. Suhel Mehandi did the statistical analysis. Yaman Kumar, Etalesh Goutam, Amrita Kumari, Rupinder Singh and Rajeev prepared the manuscript. All the authors approved the final version of the manuscript.
 
Disclaimers
 
The views and conclusions expressed in this article are solely those of the authors and do not necessarily represent the views of their affiliated institutions. The authors are responsible for the accuracy and completeness of the information provided, but do not accept any liability for any direct or indirect losses resulting from the use of this content.
 
Informed consent
 
All animal procedures for experiments were approved by the Committee of Experimental Animal care and handling techniques were approved by the University of Animal Care Committee.
The authors declare that there are no conflicts of interest regarding the publication of this article. No funding or sponsorship influenced the design of the study, data collection, analysis, decision to publish, or preparation of the manuscript. 

  1. Adgo, E. and Schulze, J. (2002). Nitrogen fixation and assimilation efficiency in Ethiopian and German pea varieties. Plant and Soil. 239: 291-299. doi: 10.1023/A:1015048331366.

  2. Al-Jibouri, H.A., Miller, P.A. and Robinson, H.F. (1958). Genotypic and environmental variances and covariances in an upland cotton cross of interspecific origin. Agronomy Journal. 50: 633-636.

  3. Aziz-ur-Rahman, Katoch, V., Rathour, R., Sharma, S., Rana, S.S. and Sharma, A. (2021). Studies on genetic parameters, correlation and path coefficient analysis in er2 introgressed garden pea genotypes. Legume Research. 44(6): 621-626. doi: 10.18805/LR-4142.

  4. Azmat, M.A., Nawab, N.N., Khan, A.A., Ashraf, M., Niaz, S. and Mahmood, K. (2011). Characterization of pea germplasm. International Journal of Vegetable Science. 17: 246- 258. doi: 10.1080/19315260.2010.544380.

  5. Barcchiya, J., Naidu, A.K., Mehta, A.K. and Upadhyay, A. (2018). Genetic variability, heritability and genetic advance for yield and yield components in pea (Pisum sativum L.). International Journal of Chemical Studies. 6: 3324-3327.

  6. Bharti, Mal, D., Goutam, E. and Kumar, L. (2024). Performance of french bean (Phaseolus vulgaris L.) genotypes under trans-gangetic plains region. Indian Journal of Agricultural Research. 58(1): 101-105. doi: 10.18805/IJARe.A-6148.

  7. Bhatt, G.M. (1973). Significance of path coefficient analysis in determining the nature of character association. Euphytica. 22: 338-343. doi: 10.1007/BF00022643.

  8. Bijalwan, P., Raturi, A. and Mishra, A.C. (2018). Character association and path analysis studies in garden pea (Pisum sativum L.) for yield and yield attributes. International Journal of Current Microbiology and Applied Sciences. 7: 3491- 3495. doi: 10.20546/ijcmas.2018.703.401.

  9. Burton, C.W. and De Vane, E.H. (1953). Estimating heritability in tall fescue (Festuca arundinacae) from replicated clonal material. Agronomy Journal. 45: 478-481.

  10. Dewey, D.R. and Lu, K. (1959). A correlation and path-coefficient analysis of components of crested wheatgrass seed production. Agronomy Journal. 51: 515-518.

  11. Georgieva, N., Nikolova, I. and Kosev, V. (2016). Evaluation of genetic divergence and heritability in pea (Pisum sativum L.). Journal of BioScience and Biotechnology. 5: 61-67.

  12. Ghobary, H.M.M. (2010). Study of relationship between yield and some yield components in garden pea (Pisum sativum L.) by using correlation and path analysis. Journal of Agricultural Research. 36: 351-360.

  13. Gudadinni, P., Bahadur, V., Ligade, P., Topno, S.E. and Prasad, V.M. (2017). Study on genetic variability, heritability and genetic advance in garden pea (Pisum sativum var. hortense L.). International Journal of Current Microbiology and Applied Sciences. 6: 2384-2391. doi: 10.20546/ ijcmas.2017.603.282.

  14. Gupta, A., Singh, B., Kumar, M., Chand, P., Vaishali and Sharma, V.R. (2020a). Studies on genetic variability, heritability and genetic advance in table pea (Pisum sativum var. hortense L.). International Journal of Current Microbiology and Applied Sciences. 9: 3449-3454. doi: 10.20546/ ijcmas.2020.909.428.

  15. Gupta, A., Singh, B., Prakash, S., Rout, S. and Singh, J. (2020b). Correlation coefficient and path analysis in table pea (Pisum sativum var. hortense L.). Current Journal of Applied Science and Technology. 39: 56-60. doi: 10.9734/ CJAST/2020/v39i3130989.

  16. Jagadeesh, K., Mahto, C.S. and Kumar, N. (2023). Genetic variability studies in field pea (Pisum sativum L.) for yield and associated characters. Environment Conservation Journal. 24: 244-249. doi: 10.36953/ECJ.13172379.

  17. Jaiswal, N.K., Gupta, A.K., Dewangan, H. and Lavanya, G.R. (2013). Genetic variability analysis in field pea (Pisum sativum L.). International Journal of Scientific Research. 4: 2006- 2007.

  18. Johnson, H.W., Robinson, H.F. and Comstock, R.E. (1955). Estimates of genetic and environmental variability in soyabean. Agronomy Journal. 47: 314-318.

  19. Katoch, V., Singh, P., Devi, M.B., Sharma, A., Sharma, G.D. and Sharma, J.K. (2016). Study of genetic variability, character association, path analysis and selection parameters for heterotic recombinant inbred lines of garden peas (Pisum sativum var. hortense L.) under mid-hill conditions of Himachal Pradesh, India. Legume Research. 39(2): 163- 169. doi: 10.18805/lr.v0iOF.6775.

  20. Khan, M.R.A., Mahmud, F., Reza, M.A., Mahbub, M.M., Shirazy, B.J. and Rahman, M.M. (2017). Genetic diversity, correlation and path analysis for yield and yield components of pea (Pisum sativum L.). World Journal of Agricultural Sciences. 13: 11-16. doi: 10.5829/idosi.wjas.2017.11.16.

  21. Kumar, R., Singh, B., Tomar, H., Kumar, V., Chaudhary, A.K., Marwah, A. and Yashraj. (2023). Genetic heritability, variability, genetic advance, correlation and path analysis assessment in garden pea (Pisum sativum L.). International Journal of Plant and Soil Science. 35: 69-79. doi: 10.9734/ijpss/ 2023/v35i143022.

  22. Kumar, T.N.V., Alloli, T.B., Hadimani, H.P., Ajjappalavar, P.S., Satish, D., Abdul, K. and Hanchinamani, C.N. (2019). Studies on genetic variability, heritability and genetic advance in garden pea (Pisum sativum L.) varieties. International Journal of Current Microbiology and Applied Sciences. 8: 3032-3038. doi: 10.20546/ijcmas.2019.812.353.

  23. Kumar, V., Singh, J., Srivastava, C.P. and Lal, C. (2018). Genetic variability, heritability and genetic advance for seed yield and related traits in diverse genotypes of pea (Pisum sativum L.). International Journal of Chemical Studies. 6: 1174-1178.

  24. Mehandi, S. and Mishra, S.P. (2023). Genetic diversity and population structure analyses in mungbean (Vigna radiata L. Wilczek). Journal of Food Chemistry and Nanotechnology. 9(S1): S316-S322. doi: 10.17756/jfcn.2023-s1-041.

  25. Mohanty, A.T., Verma, A. and Pandey, B. (2020). Evaluation of correlation and path coefficients analysis for yield attributing traits in garden pea (Pisum sativum L.) under Tarai tract of Uttarakhand. International Journal of Current Microbiology and Applied Sciences. S11: 1315- 1322.

  26. Neupane, A., Tripathi, V., Swain, A.R., Sharma, J.K., Kumar, Y. and Goutam, E. (2023). Varietal mean performance of garden pea (Pisum sativum var. hortense L.) under subtropical conditions of Punjab, India. International Journal of Environment and Climate Change. 13: 1746-1755. doi: 10.9734/IJECC/2023/v13i113558.

  27. Nwangburuka, C.C., Kehinde, O.B., Ojo, D.K., Denton, O.A. and Popoola, A.R. (2011). Morphological classification of genetic diversity in cultivated okra [Abelmoschus esculentus (L.) Moench] using PCA and SLCA. African Journal of Biotechnology. 10: 11165-11172. doi: 10.5897/ AJB11.285.

  28. Pal, A.K. and Singh, S. (2012). Correlation and path analysis in garden pea (Pisum sativum L. var. hortense). Asian Journal of Horticulture. 7: 569-573.

  29. Pandey, M., Singh, V.B., Yadav, G.C., Tyagi, N., Vishen, G.S., Bhargav, K.K. and Pandey, P. (2017). Correlation and path coefficient analysis among different characters in genotypes of vegetable pea. Vegetable Science. 44: 139-141.

  30. Panse, V.G. and Sukhatme, P.V. (1961). Statistical Methods for Agricultural Workers. ICAR, New Delhi. pp. 1-328.

  31. Patil, P., Shrivastav, S.P., Patil, K., Landge, R. and Gurjar, D. (2024). Genetic variability, heritability, genetic advance and divergence analysis in wheat (Triticum aestivum L.). Indian Journal of Agricultural Research. 58(2): 209-214. doi: 10.18805/IJARe.A-6036.

  32. Saxena, R.R., Vidayakar, V., Vishkarma, M.V., Yadav, P.S., Meena, M.L. and Lal, G.M. (2014). Genetic variability and heritability analysis for some quantitative traits in field pea (Pisum sativum L.). Bioscan. 9: 895-898.

  33. Sharma, V.K. and Bora, L. (2013). Studies on genetic variability and heterosis in vegetable pea (Pisum sativum L.) under high hills condition of Uttarakhand, India. African Journal of Agricultural Research. 8: 1891-1895.

  34. Singh, S., Singh, B., Sharma, V.R., Kumar, M. and Sirohi, U. (2022). Assessment of genetic diversity and population structure in pea (Pisum sativum L.) germplasm based on morphological traits and SSR markers. Legume Research. 45(6): 683- 688. doi: 10.18805/LR-4751.

  35. Singh, S., Verma, V., Singh, B., Sharma, V.R. and Kumar, M. (2019). Genetic variability, heritability and genetic advance studies in pea (Pisum sativum L.) for quantitative characters. Indian Journal of Agricultural Research. 53(5): 542-547. doi: 10.18805/IJARe.A-5245.

  36. Srivastava, A., Nath, S., Singh, P. and Sharma, P. (2025). Assessment of yield-contributing traits and genetic variability in field pea (Pisum sativum L.) under heat stress using half- diallel derived genotypes. Legume Research. 48(10): 1627-1634. doi: 10.18805/LR-5526.

  37. Thakur, S., Thakur, R. and Mehta, K.D. (2016). Genetic variability and association studies for green pod yield and component horticultural traits in garden pea under high hill dry temperate conditions of Tabo valley of Spiti district of Himachal Pradesh. International Journal of Science Environment and Technology. 5: 1987-1992.

  38. Thapa, U., Nandi, S. and Gurung, D. (2020). Performance of garden pea (Pisum sativum var. hortense L.) genotypes and their genetic variation, genetic advance, character association and path analysis for pod, seed and quality characters. International Journal of Current Microbiology and Applied Sciences. 9: 3572-3586. doi: 10.20546/ ijcmas.2020.909.442.

  39. Wright, S. (1934). The method of path coefficients. Annals of Mathematical Statistics. 5: 161-215.

  40. Zohary, D. and Hopf, M. (1973). Domestication of pulses in the old world: Legumes were companions of wheat and barley when agriculture began in the near East. Science. 182: 887-894. doi: 10.1126/science.182.4115.887.
In this Article
Published In
Agricultural Science Digest

Editorial Board

View all (0)