Molecular Characterization of Cowpea Genotypes for Genetic Diversity using RAPD and ISSR Techniques

M
I
Ishu Kumar Khute1
R
Roshan Parihar2,*
S
Sunil Kumar Nag1
S
Shubha Banerjee3
N
Nirmodh Prabha4
1Department of Genetics and Plant Breeding, College of Agriculture, Indira Gandhi Krishi Vishwavidyalaya, Raipur-492 012,  Chhattisgarh, India.
2Zonal Agricultural Research Station, (Indira Gandhi Krishi Vishwavidyalaya, Raipur), Kumhrawand, Jagdalpur-494 005, Chhattisgarh, India.
3Department of Plant Molecular Biology and Biotechnology, College of Agriculture, Indira Gandhi Krishi Vishwavidyalaya, Raipur-492 012, Chhattisgarh, India.
4College of Agriculture and Research Station, (Indira Gandhi Krishi Vishwavidyalaya, Raipur), Kawardha, Raipur-492 012,  Chhattisgarh, India.
  • Submitted04-12-2025|

  • Accepted19-02-2026|

  • First Online 28-03-2026|

  • doi 10.18805/LR-5619

Background: Cowpea (Vigna unguiculata (L.) Walp) is an important pulse crop with wide nutritional and agronomic value. Genetic diversity is essential for its improvement and molecular markers such as RAPD and ISSR provide valuable information for assessing phylogenetic relationships and species distinctiveness.

Methods: In the present investigation, 78 cowpea germplasm lines, including three check varieties, were evaluated in a randomized complete block design with three replications. Out of these, 48 diverse genotypes were selected for molecular characterization using 22 primers (10 RAPD and 12 ISSR). DNA was extracted from young leaves and subjected to PCR amplification. Nine RAPD primers generated 1180 scoreable bands, while eight ISSR primers amplified 1215 bands across the 48 genotypes.

Result: The polymorphic information content (PIC) values for RAPD primers ranged from 0.27 (RAPD-1) to 0.41 (RAPD-4), with a mean of 0.35. ISSR primers showed PIC values between 0.23 (UBC 815) and 0.42 (UBC 835), with a mean of 0.35. Similarity coefficients among genotypes based on RAPD markers ranged from 0.58 to 1.00, while ISSR-based similarity values ranged from 0.55 to 0.97 with a mean of 0.729. Dendrogram analysis grouped the genotypes into two major clusters. Combined RAPD and ISSR data produced similarity coefficients between 0.58 and 0.94 and the integrated dendrogram revealed 78–83% similarity with individual marker systems. The results indicate a high degree of genetic variation among the cowpea accessions, which can be effectively utilized in pre-breeding programs to broaden the genetic base and facilitate selection of superior lines.

Cowpea (Vigna unguiculata (L.) Walp.) (2n = 22) is an important pulse crop belonging to the often cross-pollinated category. It plays a vital role in global food security and human nutrition due to its rich nutritional and nutraceutical properties (Olabanji et al., 2018). Cowpea is an annual herb with a strong tap root system and spreading growth habit. It is a multipurpose crop used as food for humans, feed and fodder for livestock and as green manure. The crop also provides livelihood opportunities for farmers and small-scale entrepreneurs (Moolendra et al., 2018). Cowpea is popularly known as southern pea or black-eyed pea.
       
Globally, cowpea was cultivated on about 15 million hectares in 2022, with more than 90 percent of the area located in Africa (FAOSTAT, 2023). World production was approximately 9.8 million tonnes, with Niger and Nigeria being the major producers, occupying nearly 6 and 4.8 million hectares, respectively (FAO, 2023). In India, cowpea is grown on about 1.5 million hectares with an annual production of 2.21 million tonnes and an average productivity of 683 kg ha{ ¹ (Kaushik et al., 2016). The crop is widely cultivated in states such as Maharashtra, Gujarat, Rajasthan, Punjab, Haryana, Madhya Pradesh, Kerala, Karnataka, Andhra Pradesh and Tamil Nadu. Chhattisgarh is one of the important cowpea-growing states, contributing about 0.81 per cent of national area and 0.25 per cent of production, with approximately 43.13 thousand hectares under cultivation and an average productivity of 314 kg ha-1 (Agriculture Department, Government of Chhattisgarh; FAOSTAT, 2023).
       
Genetic diversity is the major factor influencing population improvement and selection efficiency. Assessment of genetic variability in cowpea is essential for developing improved varieties with higher yield and stress tolerance. However, cultivated cowpea possesses a narrow genetic base due to its predominantly self-pollinating nature and domestication bottlenecks (Govindaraj et al., 2015). Therefore, evaluation of genetic diversity is crucial for effective conservation and utilization of genetic resources (Khan et al., 2012).
       
Molecular marker techniques such as RAPD (Random Amplified Polymorphic DNA) and ISSR (Inter Simple Sequence Repeat) are widely used for assessing genetic diversity and establishing phylogenetic relationships. RAPD markers are simple, rapid and cost-effective tools for detecting genetic variation among cultivars (Paiva et al., 1994; Subramanian et al., 2000; Dwivedi et al., 2001). ISSR markers target microsatellite regions and provide highly reproducible and polymorphic multilocus profiles. DNA polymorphism-based molecular techniques have become essential for characterizing germplasm and identifying genetic variation in breeding materials. The development of molecular tools that directly detect DNA-level variation has significantly improved the efficiency and accuracy of genetic diversity analysis in crops lacking extensive classical genetic information (Nath et al., 2017; Asadova et al., 2024, Ofem et al., 2025).
Genotype collection
 
78 germplasm lines in all (Table 1), including three check varieties, were assessed using a randomized block design with three replicates at Research cum Instructional Farm, Department of Genetics and Plant Breeding, College of Agriculture, Indira Gandhi Krishi Vishwavidalaya, Raipur, Chhattisgarh. Out of these only 48 selected diverse genotypes of cowpea were analyzed to assess the molecular diversity.

Table 1: List of 78 cowpea genotypes collected from different locations of the Chhattisgarh state.


 
DNA extraction and electrophoresis
 
From each of the 48 genotypes, young and healthy leaves were harvested independently and the modified mini prep CTAB procedure was used to extract the genomic DNA. Place the chopped leaf pieces in a 2.0 ml Eppendorf tube, add 700 μl of CTAB buffer and then leave it at 4°C for three to four hours. Next, Put the leaf in a water bath at 65°C for 20 minutes, grind it if necessary, add extra CTAB buffer, add 700 μl of chloroform: iso amyl alcohol (24:1), vortex the sample and centrifuge it for 20 minutes at 14,000 rpm. Next, move the top layer of clear material into a fresh 1.5 ml Eppendorf tube, add 700 μl of ice-cold isopropanol in and 175 μl of potassium acetate and incubate for two hours at 4°C. After 20 minutes of centrifuging at 14,000 rpm and discarding the supernatant, wash the pellet with 50 μl of 70% ethanol and centrifuge it once more for 10 minutes after the centrifugation air dry the pellet and add 100 μl of TE buffer to dissolve the pellet. The Nano Drop Spectrophotometer was used to quantify the DNA samples. 2.0 μl of sample was taken and placed on the Nano Drop followed by closing the lid, enter 1st blank (TE) whose range should be 0.0 to 0.1. After quantification, the final DNA concentration for PCR amplification was 50 ηg/μl after the DNA was diluted with TE or sterile water (nuclease free water). 1 μl of diluted template DNA of each genotype was dispensed at the bottom of PCR plate (AXYGEN) for both RAPD and ISSR primers. About 19 μl to each sample and the PCR was performed.
       
Initial denaturation at 93ºC for 3 min and again denaturation step is carried out for 1 min. followed by annealing for 1 min at 36ºC. The amplification was completed with 2 min final extension for 72ºC. The electrophoresis of the amplification products was performed on a 1.2% agarose gel that was produced in 1X TAE buffer and stained with ethidium bromide. After electrophoresis, the products are visualized under UV transilluminator and recorded under GEL-DOC system.
 
Molecular marker RAPD and ISSR assay
 
Nine RAPD primer and eight ISSR primers were used for PCR reaction (Table 2).

Table 2: Total number of bands and PIC value of RAPD and ISSR markers.

The PIC (Polymorphic information content) index indicated extent of polymorphic bands generated by a primer has been used extensively in many genetic diversity studies. It depends on the number of detectable alleles and distribution of their frequency and is equivalent to gene diversity. The abundant allelic variation is desirable for rapid genotyping in large populations with limited primers, since it favors the attainment of high PIC, especially if the markers are associated with gene of interest. The moderate values of PIC for the RAPD primers and ISSR primers could be attributed to the diverse nature of the cowpea genotypes and/or highly informative RAPD and ISSR markers used in this study. A total number of 22 primers (RAPD 10 primers, ISSR 12 primers) were used to estimate genetic variability between cowpea genotypes (Table 2). Out of 10 RAPD markers 9 RAPD marker and 08 out of 12 ISSR markers are amplified in all genotypes. A total 1180 bands were observed among the 48 cowpea genotypes based on RAPD analysis with 9 primers. There were 1215 scoreable bands in all from ISSR markers. The likelihood that polymorphism would occur between two randomly chosen genotypes at a given locus is determined by using the PIC value as a measure of variability at that locus. With an average of 0.35, the estimated PIC value varied from 0.27 for primer RAPD 1 to 0.41 for RAPD 4 on the other hand, the average PIC for ISSR primers was 0.35, with a range of 0.23 for primer UBC 815 to 0.42 for primer UBC 835. With an average of 131.11 bands per primer, the number of scoreable markers generated varied from 41 to 226.
       
The similarity coefficient analysis was performed unit UPGMA coefficient (NT SYSPC 2014, software) based on 9 RAPD markers ranged from 0.58 to 1.00 (Fig 1). The genotype RCC63 and RCC64 showed the highest similarity index (1.00) and CP1 and CP41 showed the lowest similarity index (0.58). The mean similarity index was 0.79. RAPD markers to generate a dendrogram of cowpea genotypes. The dendrogram separated cowpea genotypes and showing two clusters; the cluster I (24 genotypes) divided into the two sub cluster A (4 genotypes) and B (20 genotypes). Sub cluster A further divided two cluster a‘ and a‘‘ with cluster a‘ containing only one genotypes RCC47 and Cluster a‘‘ having 3 genotype RCC35, RCC51, RCC40. Sub cluster B included two cluster b‘ and b‘‘. Cluster b‘ have 3 genotype Vellyani-1, RCC5 and RCC6 whereas 17 genotypes containing in cluster b‘‘. Among the 48 genotypes, eight ISSR markers yielded 1215 scoreable bands. With an average of 151.88 bands per primer, the number of scoreable markers generated varied from 57 to 330 (Similar results also reported by Nounagnon et al., 2024 and Munyao 2023). Whereas for ISSR marker ISSR primer the similarity coefficients for the cowpea genotypes based on 1215 ISSR markers ranged between 0.55 and 0.97 (Fig 2). The mean similarity index was 0.729. The genetic similarity matrix was used to obtain dendrogram through UPGMA based cluster analysis showing two clusters; the cluster II (24 genotypes) further divided sub cluster C (1 genotype) and D (23 genotypes). Sub cluster C have only one genotype RCC66. Sub cluster D included two cluster d‘ and d‘‘. Cluster d‘ contained the 10 genotypes and the cluster d‘‘ included 13 genotype. Cluster I (24 genotypes) further divided into two sub cluster A (3 genotypes) and B (21 genotypes). The sub cluster A classified into two cluster a‘ and a‘‘. Cluster a‘ included only one genotype RCC51 and cluster a‘‘ have two genotype RCC47, RCC49. Sub cluster B also further divided into two cluster b‘ and b‘‘. The cluster b‘ containing 15 genotypes, whereas b‘‘ containing 6 genotypes MFC-09-3, BL1, BL2, RCC40, RCC47 and RCC40. Highest genetic distance between genotype MFC-8-14 and TSFC-12-15 and the least distance UPC-5286 with MFC-09-3, RCC7 with RCC8 and RCC61 with RCC68 (0.97) In general, the genetic similarity between the majorities of the genotypes analyzed was quite high. The similarity coefficient analysis was performed using UPGMA coefficient (NT SYSPC 2014, software) based on 9 RAPD markers ranged from 0.58 to 1.00 which indicates there was high variability among the cowpea cultivar under study. Similar results also reported by Raza et al. (2019) and Giachino (2019) for many studies, including population differentiation and fingerprinting. To determine the genetic similarity and dissimilarity between all the genotypes, the dendogram divided the cowpea genotypes into two clusters. Similarity coefficients for the cowpea genotypes based on 1215 bands of ISSR markers ranged between 0.55 and 0.97. Similar results also reported by Vyas et al., 2016, Diallo et al., 2024, Saxena and Rukam 2020). The mean similarity index was 0.729. Comparative analysis of RAPD and ISSR markers was done. RAPD and ISSR primers amplified an average of 131.11 and 151.88 fragments, respectively (Plate 1-8). The range of the similarity coefficient was 0.58 to 0.94. The combined cluster of ISSR and RAPD was 78–83% similar to each of the individual cluster analyses of ISSR or RAPD. Similar findings also reported by Moolendra et al. (2018); Ba et al., (2004); Gajera et al. (2014) and Sharma et al. (2014) and Khare et al(2024).

Fig 1: UPGMA based clustering of 48 cowpea genotypes using genetic distance matrix for RAPD markers.



Fig 2: UPGMA based clustering of 48 cowpea genotypes using genetic distance matrix for ISSR markers.



Plate 1: Polymorphic RAPD banding pattern for the locus RAPD 2, L=Molecular size ladder 100 bp (CP1 to CP24).



Plate 2: Polymorphic RAPD banding pattern for the locus RAPD 2, L=Molecular size ladder 100 bp. (CP25 to CP48).



Plate 3: Polymorphic RAPD banding pattern for the locus RAPD 9, L=Molecular size ladder 100 bp (CP1 to CP24).



Plate 4: Polymorphic RAPD banding pattern for the locus RAPD 9, L= Molecular size ladder 100 bp (CP25 to CP48).



Plate 5: Polymorphic ISSR banding pattern for the locus UBC 809, L=Molecular size ladder 100 bp (CP1 to CP24).



Plate 6: Polymorphic ISSR banding pattern for the locus UBC 809, L=Molecular size ladder 100 bp (CP25 to CP48).



Plate 7: Polymorphic ISSR banding pattern for the locus UBC 810, L=Molecular size ladder 100 bp (CP1 to CP24).



Plate 8: Polymorphic ISSR banding pattern for the locus UBC 810, L=Molecular size ladder 100 bp (CP25 to CP 48).

The likelihood that polymorphism would occur between two randomly chosen genotypes at a given locus is determined by using the PIC value as a measure of variability at that locus. Cowpea genotypes were divided into two clusters by the dendogram and the combined cluster of RAPD and ISSR was 78–83% similar to each of the individual cluster analyses (RAPD or ISSR). The estimated PIC value ranged from 0.27 for primer RAPD 1 to 0.41 for RAPD 4 with an average of 0.35. Whereas for ISSR primers PIC ranged from 0.23 for primer UBC 815 to 0.42 for UBC 835 with average of 0.35. The number of scorable markers produced per primer ranged from 41 to 226 with an average of 131.11 bands per primer. Eight ISSR markers produced a total of 1215 scorable bands among the 48 genotypes. The number of scorable markers produced per primer ranged from 57 to 330 with an average of 151.88 bands per primer. Based on the current findings, it is clear that the cowpea accessions exhibit a high degree of genetic diversity. Therefor it will be helpful in pre-breeding for generating variability for selection of superior lines.
       
Based on the overall results, both RAPD and ISSR markers were found to be effective in detecting genetic polymorphism among the cowpea genotypes. However, ISSR markers generated a higher average number of scorable bands and exhibited slightly higher polymorphism levels and similarity resolution compared to RAPD markers. The relatively higher PIC values and greater number of amplified fragments obtained with ISSR primers indicate that ISSR markers are comparatively more informative and reliable for assessing genetic diversity in cowpea. Identification of greater genetic distance and higher levels of polymorphism among genotypes can significantly facilitate effective hybridization. Utilization of genetically diverse parents in breeding programs enhances the chances of obtaining superior recombinants and broadens the genetic base, thereby contributing to the development of improved cowpea varieties.
The authors sincerely acknowledge the facilities and support provided by the Research cum Instructional Farm and Molecular biology Laboratory of Department of Genetics and Plant Breeding, College of Agriculture, Indira Gandhi Krishi Vishwavidyalaya, Raipur, Chhattisgarh, for successful conduct of the field experiments and for providing the laboratory facilities for carrying out RAPD and ISSR marker analysis. The authors express their heartfelt gratitude to the Head of the Department and all faculty members for their constant encouragement, guidance and valuable suggestions throughout the course of this research work. We also extend our sincere thanks to the technical and supporting staff for their cooperation and assistance in completion of the study.
 
Disclaimers
 
The views and conclusions expressed in this article are solely those of the authors and do not necessarily represent the views of other research 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.
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. This study did not involve any experiments on humans or animals. Therefore, ethical approval and informed consent were not required for the present research work.

  1. Asadova, A., Babayeva, S. and Abbasov, M. (2024). Molecular Characterization of [Vigna unguiculata (L.) Walp] collection based on ISSR markers. Agricultural Reviews. 45(2): 354-359. doi: 10.18805/ag.RF-241.

  2. Ba, F.S., Remy, S. Pasquet and Paul, G. (2004). Genetic diversity in cowpea [Vigna unguiculata (L.) Walp.] as revealed by RAPD markers. Genetic Resources and Crop Evolution. 51: 539-550.

  3. Diallo, S., Badiane, F.A., Kabkia, B.N.P.A., Diédhiou, I., Diouf, M. and Diouf, D. (2024). Genetic diversity and population structure of cowpea mutant collection using SSR and ISSR molecular markers. Scientific Reports. 14(1): 31833.

  4. Dwivedi, S.L., Gurtu S., Chandra S., Yuejin, W. and Nigam S.N. (2001). Assessment of genetic diversity among selected groundnut germplasm. I: RAPD analysis. Plant Breeding. 120: 345-349.

  5. Food and Agriculture Organization of the United Nations. (2023). Production: Crops and livestock products. In: FAOSTAT. Rome.

  6. Gajera, H.P., Domadiya, R.K., Patel, S.V. and Golakiya, B.A. (2014). Appraisal of RAPD and ISSR markers for genetic diversity analysis among cowpea (Vigna unguiculata L.) genotypes. Journal of Crop Science and Biotechnology. 17(2): 79-88.

  7. Govindaraj, M., Vetriventhan, M. and Srinivasan, M. (2015).  Importance of genetic diversity assessment in crop plants and its recent advances: An overview of its analytical perspectives. Genet. Res. Int. pp 1-14.

  8. Giachino, R.R.A. (2019). Investigation of the genetic variation of anise (Pimpinella anisum L.) using RAPD and ISSR markers. Genetic Resources and Crop Evolution. 67: 763-780. 

  9. Kaushik, A.K., Bisht, K., Yadav, S.K. and Srivastava, P., (2016). Impact of various insecticides on natural enemies in cowpea ecosystem. Journal of Plant Development Sciences. 8: 547-550.

  10. Khan, M.A., Han Y, Zhao, Y.F. and Korban, S.S. (2012). A high- throughput apple SNP genotyping platform using the golden gate TM assay. Genetics. 494: 196-201.

  11. Khare, A., Gupta, N., Srivastava, N. and Bhagyawant, S.S. (2024). Genetic polymorphism amongst moth bean [Vigna aconitifolia (Jacq.) Marechal] accessions evaluated using various molecular markers. Ecological Genetics and Genomics. 33: 100284.

  12. Nath, A., Maloo, S.R., Meena, B.L., Gangarani, D.A. and Tak, S. (2017). Assessment of genetic diversity using issr markers in green gram [Vigna radiata (L.) Wilczek]. Int. J. Curr. Microbiol. App. Sci. 6(5): 1150-1158.

  13. Nounagnon, M., Roko, G., Agbodjato, N.A., Dah-Nouvlessounon, D., Babalola, O.O. and Baba-Moussa, L. (2024). Cowpea (Vigna unguiculata). Potential Pulses: Genetic and Genomic Resources. pp 58-77.

  14. Moolendra, S.T., Mali, A. and Tandan, N. (2018). RAPD based genetic diversity assessment of cowpea. Plant Archives. 18(2): 2383-2388.

  15. Munyao, R.K. (2023). Assessment of Diversity among Cowpea Accessions from Semi-Arid Areas of Kenya (Doctoral Dissertation, JKUAT-COANRE).

  16. Ofem, N.P., Kamal, N.M., Pearson, S., Shatte, T., Jordan, D., Mace, E. and Ishii, T. (2025). Genetic diversity and structure of Asian cowpea germplasm. Scientific Reports. 15(1): 27909.

  17. Olabanji, I.O., Ajayi, O.S., Oluyemi, E.A., Olawuni, I.J., Adeniji, A.O., Olasupo, O.F. and Olusesi, I.M. (2018). Nutraceuticals in different varieties of cowpeas. American Journal of Food Science and Technology. 6(2): 68-75.

  18. Paiva E., Vilarinhos A.D., Debarros E.G., Sediyama C.S. and Moreira M.A. (1994). Use of the random amplified polymorphic DNA technique to characterize soybean [Glycine max (L.) Mer- rill] genotypes. Brazilian Journal of Genetics. 17: 287-290.

  19. Raza, A.F., Ashraf, S.S.M. and Khan, R.S.A. (2019). Evaluation of genetic diversity among exotic sorghum (Sorghum bicolor L. Moench) genotypes through molecular based analysis (RAPD-PCR). Journal of Crop Health. 71: 187- 196.

  20. Saxena, A. and Rukam, T.S. (2020). Assessment of genetic diversity in cowpea [Vigna unguiculata (L.) Walp.] through ISSR marker. Research Journal of Biotechnology. 13(3): 66-71.

  21. Sharma, P., Sharma, V. and Kumar, V. (2014). Genetic diversity analysis of cluster bean [Cyamopsis tetragonoloba (L.) Taub] genotypes using RAPD and ISSR markers. Journal of Agricultural Science and Technology. 16(2): 433-443.

  22. Subramanian, V., Gurtu S., Nageswara, R.R.C. and Nigam, S.N. (2000). Identification of DNA polymorphism in cultivated groundnut using random amplified polymorphic DNA (RAPD) assay. Genome. 43: 656-660.

  23. Vyas, D., Joshi, A., Rajamani, G., Jain, D. and Kaur, G. (2016). Genetic diversity analysis in different genotypes of black gram [Vigna mungo (L.) Hepper] using RAPD marker. Legume Research. 39(5): 690-697. doi: 10.18805/lr.v0iOF.9430.

Molecular Characterization of Cowpea Genotypes for Genetic Diversity using RAPD and ISSR Techniques

M
I
Ishu Kumar Khute1
R
Roshan Parihar2,*
S
Sunil Kumar Nag1
S
Shubha Banerjee3
N
Nirmodh Prabha4
1Department of Genetics and Plant Breeding, College of Agriculture, Indira Gandhi Krishi Vishwavidyalaya, Raipur-492 012,  Chhattisgarh, India.
2Zonal Agricultural Research Station, (Indira Gandhi Krishi Vishwavidyalaya, Raipur), Kumhrawand, Jagdalpur-494 005, Chhattisgarh, India.
3Department of Plant Molecular Biology and Biotechnology, College of Agriculture, Indira Gandhi Krishi Vishwavidyalaya, Raipur-492 012, Chhattisgarh, India.
4College of Agriculture and Research Station, (Indira Gandhi Krishi Vishwavidyalaya, Raipur), Kawardha, Raipur-492 012,  Chhattisgarh, India.
  • Submitted04-12-2025|

  • Accepted19-02-2026|

  • First Online 28-03-2026|

  • doi 10.18805/LR-5619

Background: Cowpea (Vigna unguiculata (L.) Walp) is an important pulse crop with wide nutritional and agronomic value. Genetic diversity is essential for its improvement and molecular markers such as RAPD and ISSR provide valuable information for assessing phylogenetic relationships and species distinctiveness.

Methods: In the present investigation, 78 cowpea germplasm lines, including three check varieties, were evaluated in a randomized complete block design with three replications. Out of these, 48 diverse genotypes were selected for molecular characterization using 22 primers (10 RAPD and 12 ISSR). DNA was extracted from young leaves and subjected to PCR amplification. Nine RAPD primers generated 1180 scoreable bands, while eight ISSR primers amplified 1215 bands across the 48 genotypes.

Result: The polymorphic information content (PIC) values for RAPD primers ranged from 0.27 (RAPD-1) to 0.41 (RAPD-4), with a mean of 0.35. ISSR primers showed PIC values between 0.23 (UBC 815) and 0.42 (UBC 835), with a mean of 0.35. Similarity coefficients among genotypes based on RAPD markers ranged from 0.58 to 1.00, while ISSR-based similarity values ranged from 0.55 to 0.97 with a mean of 0.729. Dendrogram analysis grouped the genotypes into two major clusters. Combined RAPD and ISSR data produced similarity coefficients between 0.58 and 0.94 and the integrated dendrogram revealed 78–83% similarity with individual marker systems. The results indicate a high degree of genetic variation among the cowpea accessions, which can be effectively utilized in pre-breeding programs to broaden the genetic base and facilitate selection of superior lines.

Cowpea (Vigna unguiculata (L.) Walp.) (2n = 22) is an important pulse crop belonging to the often cross-pollinated category. It plays a vital role in global food security and human nutrition due to its rich nutritional and nutraceutical properties (Olabanji et al., 2018). Cowpea is an annual herb with a strong tap root system and spreading growth habit. It is a multipurpose crop used as food for humans, feed and fodder for livestock and as green manure. The crop also provides livelihood opportunities for farmers and small-scale entrepreneurs (Moolendra et al., 2018). Cowpea is popularly known as southern pea or black-eyed pea.
       
Globally, cowpea was cultivated on about 15 million hectares in 2022, with more than 90 percent of the area located in Africa (FAOSTAT, 2023). World production was approximately 9.8 million tonnes, with Niger and Nigeria being the major producers, occupying nearly 6 and 4.8 million hectares, respectively (FAO, 2023). In India, cowpea is grown on about 1.5 million hectares with an annual production of 2.21 million tonnes and an average productivity of 683 kg ha{ ¹ (Kaushik et al., 2016). The crop is widely cultivated in states such as Maharashtra, Gujarat, Rajasthan, Punjab, Haryana, Madhya Pradesh, Kerala, Karnataka, Andhra Pradesh and Tamil Nadu. Chhattisgarh is one of the important cowpea-growing states, contributing about 0.81 per cent of national area and 0.25 per cent of production, with approximately 43.13 thousand hectares under cultivation and an average productivity of 314 kg ha-1 (Agriculture Department, Government of Chhattisgarh; FAOSTAT, 2023).
       
Genetic diversity is the major factor influencing population improvement and selection efficiency. Assessment of genetic variability in cowpea is essential for developing improved varieties with higher yield and stress tolerance. However, cultivated cowpea possesses a narrow genetic base due to its predominantly self-pollinating nature and domestication bottlenecks (Govindaraj et al., 2015). Therefore, evaluation of genetic diversity is crucial for effective conservation and utilization of genetic resources (Khan et al., 2012).
       
Molecular marker techniques such as RAPD (Random Amplified Polymorphic DNA) and ISSR (Inter Simple Sequence Repeat) are widely used for assessing genetic diversity and establishing phylogenetic relationships. RAPD markers are simple, rapid and cost-effective tools for detecting genetic variation among cultivars (Paiva et al., 1994; Subramanian et al., 2000; Dwivedi et al., 2001). ISSR markers target microsatellite regions and provide highly reproducible and polymorphic multilocus profiles. DNA polymorphism-based molecular techniques have become essential for characterizing germplasm and identifying genetic variation in breeding materials. The development of molecular tools that directly detect DNA-level variation has significantly improved the efficiency and accuracy of genetic diversity analysis in crops lacking extensive classical genetic information (Nath et al., 2017; Asadova et al., 2024, Ofem et al., 2025).
Genotype collection
 
78 germplasm lines in all (Table 1), including three check varieties, were assessed using a randomized block design with three replicates at Research cum Instructional Farm, Department of Genetics and Plant Breeding, College of Agriculture, Indira Gandhi Krishi Vishwavidalaya, Raipur, Chhattisgarh. Out of these only 48 selected diverse genotypes of cowpea were analyzed to assess the molecular diversity.

Table 1: List of 78 cowpea genotypes collected from different locations of the Chhattisgarh state.


 
DNA extraction and electrophoresis
 
From each of the 48 genotypes, young and healthy leaves were harvested independently and the modified mini prep CTAB procedure was used to extract the genomic DNA. Place the chopped leaf pieces in a 2.0 ml Eppendorf tube, add 700 μl of CTAB buffer and then leave it at 4°C for three to four hours. Next, Put the leaf in a water bath at 65°C for 20 minutes, grind it if necessary, add extra CTAB buffer, add 700 μl of chloroform: iso amyl alcohol (24:1), vortex the sample and centrifuge it for 20 minutes at 14,000 rpm. Next, move the top layer of clear material into a fresh 1.5 ml Eppendorf tube, add 700 μl of ice-cold isopropanol in and 175 μl of potassium acetate and incubate for two hours at 4°C. After 20 minutes of centrifuging at 14,000 rpm and discarding the supernatant, wash the pellet with 50 μl of 70% ethanol and centrifuge it once more for 10 minutes after the centrifugation air dry the pellet and add 100 μl of TE buffer to dissolve the pellet. The Nano Drop Spectrophotometer was used to quantify the DNA samples. 2.0 μl of sample was taken and placed on the Nano Drop followed by closing the lid, enter 1st blank (TE) whose range should be 0.0 to 0.1. After quantification, the final DNA concentration for PCR amplification was 50 ηg/μl after the DNA was diluted with TE or sterile water (nuclease free water). 1 μl of diluted template DNA of each genotype was dispensed at the bottom of PCR plate (AXYGEN) for both RAPD and ISSR primers. About 19 μl to each sample and the PCR was performed.
       
Initial denaturation at 93ºC for 3 min and again denaturation step is carried out for 1 min. followed by annealing for 1 min at 36ºC. The amplification was completed with 2 min final extension for 72ºC. The electrophoresis of the amplification products was performed on a 1.2% agarose gel that was produced in 1X TAE buffer and stained with ethidium bromide. After electrophoresis, the products are visualized under UV transilluminator and recorded under GEL-DOC system.
 
Molecular marker RAPD and ISSR assay
 
Nine RAPD primer and eight ISSR primers were used for PCR reaction (Table 2).

Table 2: Total number of bands and PIC value of RAPD and ISSR markers.

The PIC (Polymorphic information content) index indicated extent of polymorphic bands generated by a primer has been used extensively in many genetic diversity studies. It depends on the number of detectable alleles and distribution of their frequency and is equivalent to gene diversity. The abundant allelic variation is desirable for rapid genotyping in large populations with limited primers, since it favors the attainment of high PIC, especially if the markers are associated with gene of interest. The moderate values of PIC for the RAPD primers and ISSR primers could be attributed to the diverse nature of the cowpea genotypes and/or highly informative RAPD and ISSR markers used in this study. A total number of 22 primers (RAPD 10 primers, ISSR 12 primers) were used to estimate genetic variability between cowpea genotypes (Table 2). Out of 10 RAPD markers 9 RAPD marker and 08 out of 12 ISSR markers are amplified in all genotypes. A total 1180 bands were observed among the 48 cowpea genotypes based on RAPD analysis with 9 primers. There were 1215 scoreable bands in all from ISSR markers. The likelihood that polymorphism would occur between two randomly chosen genotypes at a given locus is determined by using the PIC value as a measure of variability at that locus. With an average of 0.35, the estimated PIC value varied from 0.27 for primer RAPD 1 to 0.41 for RAPD 4 on the other hand, the average PIC for ISSR primers was 0.35, with a range of 0.23 for primer UBC 815 to 0.42 for primer UBC 835. With an average of 131.11 bands per primer, the number of scoreable markers generated varied from 41 to 226.
       
The similarity coefficient analysis was performed unit UPGMA coefficient (NT SYSPC 2014, software) based on 9 RAPD markers ranged from 0.58 to 1.00 (Fig 1). The genotype RCC63 and RCC64 showed the highest similarity index (1.00) and CP1 and CP41 showed the lowest similarity index (0.58). The mean similarity index was 0.79. RAPD markers to generate a dendrogram of cowpea genotypes. The dendrogram separated cowpea genotypes and showing two clusters; the cluster I (24 genotypes) divided into the two sub cluster A (4 genotypes) and B (20 genotypes). Sub cluster A further divided two cluster a‘ and a‘‘ with cluster a‘ containing only one genotypes RCC47 and Cluster a‘‘ having 3 genotype RCC35, RCC51, RCC40. Sub cluster B included two cluster b‘ and b‘‘. Cluster b‘ have 3 genotype Vellyani-1, RCC5 and RCC6 whereas 17 genotypes containing in cluster b‘‘. Among the 48 genotypes, eight ISSR markers yielded 1215 scoreable bands. With an average of 151.88 bands per primer, the number of scoreable markers generated varied from 57 to 330 (Similar results also reported by Nounagnon et al., 2024 and Munyao 2023). Whereas for ISSR marker ISSR primer the similarity coefficients for the cowpea genotypes based on 1215 ISSR markers ranged between 0.55 and 0.97 (Fig 2). The mean similarity index was 0.729. The genetic similarity matrix was used to obtain dendrogram through UPGMA based cluster analysis showing two clusters; the cluster II (24 genotypes) further divided sub cluster C (1 genotype) and D (23 genotypes). Sub cluster C have only one genotype RCC66. Sub cluster D included two cluster d‘ and d‘‘. Cluster d‘ contained the 10 genotypes and the cluster d‘‘ included 13 genotype. Cluster I (24 genotypes) further divided into two sub cluster A (3 genotypes) and B (21 genotypes). The sub cluster A classified into two cluster a‘ and a‘‘. Cluster a‘ included only one genotype RCC51 and cluster a‘‘ have two genotype RCC47, RCC49. Sub cluster B also further divided into two cluster b‘ and b‘‘. The cluster b‘ containing 15 genotypes, whereas b‘‘ containing 6 genotypes MFC-09-3, BL1, BL2, RCC40, RCC47 and RCC40. Highest genetic distance between genotype MFC-8-14 and TSFC-12-15 and the least distance UPC-5286 with MFC-09-3, RCC7 with RCC8 and RCC61 with RCC68 (0.97) In general, the genetic similarity between the majorities of the genotypes analyzed was quite high. The similarity coefficient analysis was performed using UPGMA coefficient (NT SYSPC 2014, software) based on 9 RAPD markers ranged from 0.58 to 1.00 which indicates there was high variability among the cowpea cultivar under study. Similar results also reported by Raza et al. (2019) and Giachino (2019) for many studies, including population differentiation and fingerprinting. To determine the genetic similarity and dissimilarity between all the genotypes, the dendogram divided the cowpea genotypes into two clusters. Similarity coefficients for the cowpea genotypes based on 1215 bands of ISSR markers ranged between 0.55 and 0.97. Similar results also reported by Vyas et al., 2016, Diallo et al., 2024, Saxena and Rukam 2020). The mean similarity index was 0.729. Comparative analysis of RAPD and ISSR markers was done. RAPD and ISSR primers amplified an average of 131.11 and 151.88 fragments, respectively (Plate 1-8). The range of the similarity coefficient was 0.58 to 0.94. The combined cluster of ISSR and RAPD was 78–83% similar to each of the individual cluster analyses of ISSR or RAPD. Similar findings also reported by Moolendra et al. (2018); Ba et al., (2004); Gajera et al. (2014) and Sharma et al. (2014) and Khare et al(2024).

Fig 1: UPGMA based clustering of 48 cowpea genotypes using genetic distance matrix for RAPD markers.



Fig 2: UPGMA based clustering of 48 cowpea genotypes using genetic distance matrix for ISSR markers.



Plate 1: Polymorphic RAPD banding pattern for the locus RAPD 2, L=Molecular size ladder 100 bp (CP1 to CP24).



Plate 2: Polymorphic RAPD banding pattern for the locus RAPD 2, L=Molecular size ladder 100 bp. (CP25 to CP48).



Plate 3: Polymorphic RAPD banding pattern for the locus RAPD 9, L=Molecular size ladder 100 bp (CP1 to CP24).



Plate 4: Polymorphic RAPD banding pattern for the locus RAPD 9, L= Molecular size ladder 100 bp (CP25 to CP48).



Plate 5: Polymorphic ISSR banding pattern for the locus UBC 809, L=Molecular size ladder 100 bp (CP1 to CP24).



Plate 6: Polymorphic ISSR banding pattern for the locus UBC 809, L=Molecular size ladder 100 bp (CP25 to CP48).



Plate 7: Polymorphic ISSR banding pattern for the locus UBC 810, L=Molecular size ladder 100 bp (CP1 to CP24).



Plate 8: Polymorphic ISSR banding pattern for the locus UBC 810, L=Molecular size ladder 100 bp (CP25 to CP 48).

The likelihood that polymorphism would occur between two randomly chosen genotypes at a given locus is determined by using the PIC value as a measure of variability at that locus. Cowpea genotypes were divided into two clusters by the dendogram and the combined cluster of RAPD and ISSR was 78–83% similar to each of the individual cluster analyses (RAPD or ISSR). The estimated PIC value ranged from 0.27 for primer RAPD 1 to 0.41 for RAPD 4 with an average of 0.35. Whereas for ISSR primers PIC ranged from 0.23 for primer UBC 815 to 0.42 for UBC 835 with average of 0.35. The number of scorable markers produced per primer ranged from 41 to 226 with an average of 131.11 bands per primer. Eight ISSR markers produced a total of 1215 scorable bands among the 48 genotypes. The number of scorable markers produced per primer ranged from 57 to 330 with an average of 151.88 bands per primer. Based on the current findings, it is clear that the cowpea accessions exhibit a high degree of genetic diversity. Therefor it will be helpful in pre-breeding for generating variability for selection of superior lines.
       
Based on the overall results, both RAPD and ISSR markers were found to be effective in detecting genetic polymorphism among the cowpea genotypes. However, ISSR markers generated a higher average number of scorable bands and exhibited slightly higher polymorphism levels and similarity resolution compared to RAPD markers. The relatively higher PIC values and greater number of amplified fragments obtained with ISSR primers indicate that ISSR markers are comparatively more informative and reliable for assessing genetic diversity in cowpea. Identification of greater genetic distance and higher levels of polymorphism among genotypes can significantly facilitate effective hybridization. Utilization of genetically diverse parents in breeding programs enhances the chances of obtaining superior recombinants and broadens the genetic base, thereby contributing to the development of improved cowpea varieties.
The authors sincerely acknowledge the facilities and support provided by the Research cum Instructional Farm and Molecular biology Laboratory of Department of Genetics and Plant Breeding, College of Agriculture, Indira Gandhi Krishi Vishwavidyalaya, Raipur, Chhattisgarh, for successful conduct of the field experiments and for providing the laboratory facilities for carrying out RAPD and ISSR marker analysis. The authors express their heartfelt gratitude to the Head of the Department and all faculty members for their constant encouragement, guidance and valuable suggestions throughout the course of this research work. We also extend our sincere thanks to the technical and supporting staff for their cooperation and assistance in completion of the study.
 
Disclaimers
 
The views and conclusions expressed in this article are solely those of the authors and do not necessarily represent the views of other research 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.
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. This study did not involve any experiments on humans or animals. Therefore, ethical approval and informed consent were not required for the present research work.

  1. Asadova, A., Babayeva, S. and Abbasov, M. (2024). Molecular Characterization of [Vigna unguiculata (L.) Walp] collection based on ISSR markers. Agricultural Reviews. 45(2): 354-359. doi: 10.18805/ag.RF-241.

  2. Ba, F.S., Remy, S. Pasquet and Paul, G. (2004). Genetic diversity in cowpea [Vigna unguiculata (L.) Walp.] as revealed by RAPD markers. Genetic Resources and Crop Evolution. 51: 539-550.

  3. Diallo, S., Badiane, F.A., Kabkia, B.N.P.A., Diédhiou, I., Diouf, M. and Diouf, D. (2024). Genetic diversity and population structure of cowpea mutant collection using SSR and ISSR molecular markers. Scientific Reports. 14(1): 31833.

  4. Dwivedi, S.L., Gurtu S., Chandra S., Yuejin, W. and Nigam S.N. (2001). Assessment of genetic diversity among selected groundnut germplasm. I: RAPD analysis. Plant Breeding. 120: 345-349.

  5. Food and Agriculture Organization of the United Nations. (2023). Production: Crops and livestock products. In: FAOSTAT. Rome.

  6. Gajera, H.P., Domadiya, R.K., Patel, S.V. and Golakiya, B.A. (2014). Appraisal of RAPD and ISSR markers for genetic diversity analysis among cowpea (Vigna unguiculata L.) genotypes. Journal of Crop Science and Biotechnology. 17(2): 79-88.

  7. Govindaraj, M., Vetriventhan, M. and Srinivasan, M. (2015).  Importance of genetic diversity assessment in crop plants and its recent advances: An overview of its analytical perspectives. Genet. Res. Int. pp 1-14.

  8. Giachino, R.R.A. (2019). Investigation of the genetic variation of anise (Pimpinella anisum L.) using RAPD and ISSR markers. Genetic Resources and Crop Evolution. 67: 763-780. 

  9. Kaushik, A.K., Bisht, K., Yadav, S.K. and Srivastava, P., (2016). Impact of various insecticides on natural enemies in cowpea ecosystem. Journal of Plant Development Sciences. 8: 547-550.

  10. Khan, M.A., Han Y, Zhao, Y.F. and Korban, S.S. (2012). A high- throughput apple SNP genotyping platform using the golden gate TM assay. Genetics. 494: 196-201.

  11. Khare, A., Gupta, N., Srivastava, N. and Bhagyawant, S.S. (2024). Genetic polymorphism amongst moth bean [Vigna aconitifolia (Jacq.) Marechal] accessions evaluated using various molecular markers. Ecological Genetics and Genomics. 33: 100284.

  12. Nath, A., Maloo, S.R., Meena, B.L., Gangarani, D.A. and Tak, S. (2017). Assessment of genetic diversity using issr markers in green gram [Vigna radiata (L.) Wilczek]. Int. J. Curr. Microbiol. App. Sci. 6(5): 1150-1158.

  13. Nounagnon, M., Roko, G., Agbodjato, N.A., Dah-Nouvlessounon, D., Babalola, O.O. and Baba-Moussa, L. (2024). Cowpea (Vigna unguiculata). Potential Pulses: Genetic and Genomic Resources. pp 58-77.

  14. Moolendra, S.T., Mali, A. and Tandan, N. (2018). RAPD based genetic diversity assessment of cowpea. Plant Archives. 18(2): 2383-2388.

  15. Munyao, R.K. (2023). Assessment of Diversity among Cowpea Accessions from Semi-Arid Areas of Kenya (Doctoral Dissertation, JKUAT-COANRE).

  16. Ofem, N.P., Kamal, N.M., Pearson, S., Shatte, T., Jordan, D., Mace, E. and Ishii, T. (2025). Genetic diversity and structure of Asian cowpea germplasm. Scientific Reports. 15(1): 27909.

  17. Olabanji, I.O., Ajayi, O.S., Oluyemi, E.A., Olawuni, I.J., Adeniji, A.O., Olasupo, O.F. and Olusesi, I.M. (2018). Nutraceuticals in different varieties of cowpeas. American Journal of Food Science and Technology. 6(2): 68-75.

  18. Paiva E., Vilarinhos A.D., Debarros E.G., Sediyama C.S. and Moreira M.A. (1994). Use of the random amplified polymorphic DNA technique to characterize soybean [Glycine max (L.) Mer- rill] genotypes. Brazilian Journal of Genetics. 17: 287-290.

  19. Raza, A.F., Ashraf, S.S.M. and Khan, R.S.A. (2019). Evaluation of genetic diversity among exotic sorghum (Sorghum bicolor L. Moench) genotypes through molecular based analysis (RAPD-PCR). Journal of Crop Health. 71: 187- 196.

  20. Saxena, A. and Rukam, T.S. (2020). Assessment of genetic diversity in cowpea [Vigna unguiculata (L.) Walp.] through ISSR marker. Research Journal of Biotechnology. 13(3): 66-71.

  21. Sharma, P., Sharma, V. and Kumar, V. (2014). Genetic diversity analysis of cluster bean [Cyamopsis tetragonoloba (L.) Taub] genotypes using RAPD and ISSR markers. Journal of Agricultural Science and Technology. 16(2): 433-443.

  22. Subramanian, V., Gurtu S., Nageswara, R.R.C. and Nigam, S.N. (2000). Identification of DNA polymorphism in cultivated groundnut using random amplified polymorphic DNA (RAPD) assay. Genome. 43: 656-660.

  23. Vyas, D., Joshi, A., Rajamani, G., Jain, D. and Kaur, G. (2016). Genetic diversity analysis in different genotypes of black gram [Vigna mungo (L.) Hepper] using RAPD marker. Legume Research. 39(5): 690-697. doi: 10.18805/lr.v0iOF.9430.
In this Article
Published In
Legume Research

Editorial Board

View all (0)