Molecular Diversity Assessment Using Sequence Related Amplified Polymorphism (SRAP) Markers in Vicia faba L.

Sequence-related amplified polymorphism (SRAP) markers were used to assess the genetic diversity and relationship among 58 faba bean (Vicia faba L.) genotypes. Fourteen SRAP primer combinations amplified a total of 1036 differently sized well-resolved peaks (fragments), of which all were polymorphic with a 0.96 PIC value and discriminated all of the 58 faba bean genotypes. An average pairwise similarity of 21% was revealed among the genotypes ranging from 2% to 65%. At a similarity of 28%, UPGMA clustered the genotypes into three main groups comprising 78% of the genotypes. The local landraces and most of the Egyptian genotypes in addition to the Sudan genotypes were grouped in the first main cluster. The advanced breeding lines were scattered in the second and third main clusters with breeding lines from the ICARDA and genotypes introduced from Egypt. At a similarity of 47%, all the genotypes formed separated clusters with the exceptions of Hassawi 1 and Hassawi 2. Group analysis of the genotypes according to their geographic origin and type showed that the landraces were grouped according to their origin, while others were grouped according to their seed type. To our knowledge, this is the first application of SRAP markers for the assessment of genetic diversity in faba bean. Such information will be useful to determine optimal breeding strategies to allow continued progress in faba bean breeding.     

Genetic Analysis of 430 Chinese Cynodon dactylon Accessions Using Sequence-Related Amplified Polymorphism Markers

Although Cynodon dactylon (C. dactylon) is widely distributed in China, information on its genetic diversity within the germplasm pool is limited. The objective of this study was to reveal the genetic variation and relationships of 430 C. dactylon accessions collected from 22 Chinese provinces using sequence-related amplified polymorphism (SRAP) markers. Fifteen primer pairs were used to amplify specific C. dactylon genomic sequences. A total of 481 SRAP fragments were generated, with fragment sizes ranging from 260–1800 base pairs (bp). Genetic similarity coefficients (GSC) among the 430 accessions averaged 0.72 and ranged from 0.53–0.96. Cluster analysis conducted by two methods, namely the unweighted pair-group method with arithmetic averages (UPGMA) and principle coordinate analysis (PCoA), separated the accessions into eight distinct groups. Our findings verify that Chinese C. dactylon germplasms have rich genetic diversity, which is an excellent basis for C. dactylon breeding for new cultivars.     

Genetic Differentiation between Natural and Hatchery Stocks of Japanese Scallop (Mizuhopecten yessoensis) as Revealed by AFLP Analysis

Japanese scallop (Mizuhopecten yessoensis) is a cold-tolerant bivalve that was introduced to China for aquaculture in 1982. In this study, amplified fragment length polymorphism (AFLP) markers were used to investigate levels of genetic diversity within M. yessoensis cultured stocks and compare them with wild populations. Six pairs of primer combinations generated 368 loci among 332 individuals, in four cultured and three wild populations. High polymorphism at AFLP markers was found within both cultured and wild M. yessoensis populations. The percentage of polymorphic loci ranged from 61.04% to 72.08%, while the mean heterozygosity ranged from 0.2116 to 0.2596. Compared with wild populations, the four hatchery populations showed significant genetic changes, such as lower expected heterozygosity and percentage of polymorphic loci, and smaller frequency of private alleles, all indicative of a reduction in genetic diversity. Some genetic structures were associated with the geographical distribution of samples; with all samples from Dalian and Japan being closely related, while the population from Russia fell into a distinct clade in the phylogenetic analysis. The genetic information derived from this study indicated that intentional or accidental release of selected Japanese scallops into natural sea areas might result in disturbance of local gene pools and loss of genetic variability. We recommend monitoring the genetic variability of selected hatchery populations to enhance conservation of natural Japanese scallop resources.     

Relationships among the A Genomes of Triticum L. Species as Evidenced by SSR Markers, in Iran

The relationships among 55 wheat accessions (47 accessions collected from Iran and eight accessions provided by the Institute of Plant Biology of the University of Zurich, Switzerland) belonging to eight species carrying A genome (Triticum monococcum L., T. boeoticum Boiss., T. urartu Tumanian ex Gandilyan, T. durum Desf., T. turgidum L., T. dicoccum Schrank ex Schübler, T. dicoccoides (Körn. ex Asch. & Graebner) Schweinf. and T. aestivum L.) were evaluated using 31 A genome specific microsatellite markers. A high level of polymorphism was observed among the accessions studied (PIC = 0.77). The highest gene diversity was revealed among T. durum genotypes, while the lowest genetic variation was found in T. dicoccoides accessions. The analysis of molecular variance (AMOVA) showed a significant genetic variance (75.56%) among these accessions, representing a high intra-specific genetic diversity within Triticum taxa in Iran. However, such a variance was not observed among their ploidy levels. Based on the genetic similarity analysis, the accessions collected from Iran were divided into two main groups: diploids and polyploids. The genetic similarity among the diploid and polyploid species was 0.85 and 0.89 respectively. There were no significant differences in A genome diversity from different geographic regions. Based on the genetic diversity analyses, we consider there is value in a greater sampling of each species in Iran to discover useful genes for breeding purposes.     

The Possible Physical Barrier and Coastal Dispersal Strategy for Japanese Grenadier Anchovy,Coilia nasus in the East China Sea and Yellow Sea: Evidence from AFLP Markers

In order to ascertain the taxonomic status of the Ariake Sea population of Japanese grenadier anchovy, Coilia nasus, and assess the contemporary possible genetic barrier between the west and east coastal waters of the East China Sea, we used amplified fragment length polymorphism (AFLP) markers to detect the genetic structure of C. nasus, in the East China Sea and Yellow Sea. Eighty-one individuals of C. nasus were collected from five locations and 12 individuals of Coilia mystus were sampled from the Yangtze River Estuary. A total of 371 loci were detected by five primer combinations, 310 of which were polymorphic (83.56%). Analysis of molecular variation (AMOVA) and pairwise fixation index (F_ST) revealed significant genetic differentiation among five samples, indicating limited gene flow among populations. The dendrogram for populations by neighbor-joining (NJ) cluster analysis provided evidence of a clear relationship between genetic and geographic patterns, supporting significant genetic differentiation between China coastal populations and Ariake Sea populations. Compared to the genetic divergence between C. nasus and C. mystus, the level of genetic differentiation between China and the Ariake Sea populations of C. nasus is obvious below the species level, indicating isolated populations of C. nasus in the Ariake Sea. Isolation by distance analysis revealed that direct ocean distance with deep-water at the continental slope and high salinity between west and east coastal waters of the East China Sea served as major physical barrier to C. nasus, supporting the coastal dispersal pattern in this estuarine species, and rejecting offshore dispersal strategy.     

Use of Repetitive Sequences for Molecular and Cytogenetic Characterization of Avena Species from Portugal

Genomic diversity of Portuguese accessions of Avena species—diploid A. strigosa and hexaploids A. sativa and A. sterilis—was evaluated through molecular and cytological analysis of 45S rDNA, and other repetitive sequences previously studied in cereal species—rye subtelomeric sequence (pSc200) and cereal centromeric sequence (CCS1). Additionally, retrotransposons and microsatellites targeting methodologies—IRAP (inter-retrotransposon amplified polymorphism) and REMAP (retrotransposon-microsatellite amplified polymorphism)—were performed. A very high homology was detected for ribosomal internal transcribed sequences (ITS1 and ITS2) between the species analyzed, although nucleolar organizing regions (NOR) fluorescent in situ hybridization (FISH) analysis revealed distinct number of Nor loci between diploid and hexaploid species. Moreover, morphological diversity, evidenced by FISH signals with different sizes, was observed between distinct accessions within each species. pSc200 sequences were for the first time isolated from Avena species but proven to be highly similar in all genotypes analyzed. The use of primers designed for CCS1 unraveled a sequence homologous to the Ty3/gypsy retrotransposon Cereba, that was mapped to centromeric regions of diploid and hexaploid species, being however restricted to the more related A and D haplomes. Retrotransposon-based methodologies disclosed species- and accessions-specific bands essential for the accurate discrimination of all genotypes studied. Centromeric, IRAP and REMAP profiles therefore allowed accurate assessment of inter and intraspecific variability, demonstrating the potential of these molecular markers on future oat breeding programs.     

Application of Novel Polymorphic Microsatellite Loci Identified in the Korean Pacific Abalone (Haliotis diversicolor supertexta (Haliotidae)) in the Genetic Characterization of Wild and Released Populations

The small abalone, Haliotis diversicolor supertexta, of the family Haliotidae, is one of the most important species of marine shellfish in eastern Asia. Over the past few decades, this species has drastically declined in Korea. Thus, hatchery-bred seeds have been released into natural coastal areas to compensate for the reduced fishery resources. However, information on the genetic background of the small abalone is scarce. In this study, 20 polymorphic microsatellite DNA markers were identified using next-generation sequencing techniques and used to compare allelic variation between wild and released abalone populations in Korea. Using high-throughput genomic sequencing, a total of 1516 (2.26%; average length of 385 bp) reads containing simple sequence repeats were obtained from 86,011 raw reads. Among the 99 loci screened, 28 amplified successfully, and 20 were polymorphic. When comparing allelic variation between wild and released abalone populations, a total of 243 different alleles were observed, with 18.7 alleles per locus. High genetic diversity (mean heterozygosity = 0.81; mean allelic number = 15.5) was observed in both populations. A statistical analysis of the fixation index (F_ST) and analysis of molecular variance (AMOVA) indicated limited genetic differences between the two populations (F_ST = 0.002, p > 0.05). Although no significant reductions in the genetic diversity were found in the released population compared with the wild population (p > 0.05), the genetic diversity parameters revealed that the seeds released for stock abundance had a different genetic composition. These differences are likely a result of hatchery selection and inbreeding. Additionally, all the primer pair sets were effectively amplified in another congeneric species, H. diversicolor diversicolor, indicating that these primers are useful for both abalone species. These microsatellite loci may be valuable for future aquaculture and population genetic studies aimed at developing conservation and management plans for these two abalone species.     

Genetic Variation in Safflower (Carthamus tinctorious L.) for Seed Quality-Related Traits and Inter-Simple Sequence Repeat (ISSR) Markers

Safflower (Carthamus tinctorious L.) is an oilseed crop that is valued as a source of high quality vegetable oil. The genetic diversity of 16 safflower genotypes originated from different geographical regions of Iran and some with exotic origin were evaluated. Eight different seed quality-related traits including fatty acid composition of seed oil (stearic acid, palmitic acid, oleic acid and linoleic acid), the contents of, oil, protein, fiber and ash in its seeds, as well as 20 inter-simple sequence repeat (ISSR) polymorphic primers were used in this study. Analysis of variance showed significant variation in genotypes for the seed quality-related traits. Based on ISSR markers, a total of 204 bands were amplified and 149 bands (about 70%) of these were polymorphic. Cluster analysis based on either biochemical or molecular markers classified the genotypes into four groups, showing some similarities between molecular and biochemical markers for evaluated genotypes. A logical similarity between the genotype clusters based on molecular data with their geographical origins was observed.     

Genetic Diversity and Population Structure of Siberian apricot (Prunus sibirica L.) in China

The genetic diversity and population genetic structure of 252 accessions from 21 Prunus sibirica L. populations were investigated using 10 ISSR, SSR, and SRAP markers. The results suggest that the entire population has a relatively high level of genetic diversity, with populations HR and MY showing very high diversity. A low level of inter-population genetic differentiation and a high level of intra-population genetic differentiation was found, which is supported by a moderate level of gene flow, and largely attributable to the cross-pollination and self-incompatibility reproductive system. A STRUCTURE (model-based program) analysis revealed that the 21 populations can be divided into two main groups, mainly based on geographic differences and genetic exchanges. The entire wild Siberia apricot population in China could be divided into two subgroups, including 107 accessions in subgroup (SG) 1 and 147 accessions in SG 2. A Mantel test revealed a significant positive correlation between genetic and geographic distance matrices, and there was a very significant positive correlation among three marker datasets. Overall, we recommend a combination of conservation measures, with ex situ and in situ conservation that includes the construction of a core germplasm repository and the implement of in situ conservation for populations HR, MY, and ZY.     

Molecular Characterization of Wheat Stripe Rust Pathogen (Puccinia striiformis f. sp. tritici) Collections from Nine Countries

Stripe rust, caused by Puccinia striiformis f. sp. tritici (Pst), is one of the most important diseases of wheat worldwide. To understand the worldwide distribution of its molecular groups, as well as the diversity, differentiation, and migration of the Pst populations, 567 isolates collected from nine countries (China, Pakistan, Italy, Egypt, Ethiopia, Canada, Mexico, Ecuador, and the U.S.) in 2010–2018 were genotyped using 14 codominant simple sequence repeat markers. A total of 433, including 333 new multi-locus genotypes (MLGs), were identified, which were clustered into ten molecular groups (MGs). The MGs and country-wise populations differed in genetic diversity, heterozygosity, and correlation coefficient between the marker and virulence data. Many isolates from different countries, especially the isolates from Mexico, Ecuador, and the U.S., were found to be identical or closely related MLGs, and some of the MGs were present in all countries, indicating Pst migrations among different countries. The analysis of molecular variance revealed 78% variation among isolates, 12% variation among countries, and 10% variation within countries. Only low levels of differentiation were found by the pairwise comparisons of country populations. Of the 10 MGs, 5 were found to be involved in sexual and/or somatic recombination. Identical and closely related MLGs identified from different countries indicated international migrations. The study provides information on the distributions of various Pst genetic groups in different countries and evidence for the global migrations, which should be useful in understanding the pathogen evolution and in stressing the need for continual monitoring of the disease and pathogen populations at the global scale.     

A Molecular Genetic Linkage Map of Eucommia ulmoides and Quantitative Trait Loci (QTL) Analysis for Growth Traits

Eucommia ulmoides is an economically important tree species for both herbal medicine and organic chemical industry. Effort to breed varieties with improved yield and quality is limited by the lack of knowledge on the genetic basis of the traits. A genetic linkage map of E. ulmoides was constructed from a full-sib family using sequence-related amplified polymorphism, amplified fragment length polymorphism, inter-simple sequence repeat and simple sequence repeat markers. In total, 706 markers were mapped in 25 linkage groups covering 2133 cM. The genetic linkage map covered approximately 89% of the estimated E. ulmoides genome with an average of 3.1 cM between adjacent markers. The present genetic linkage map was used to identify quantitative trait loci (QTL) affecting growth-related traits. Eighteen QTLs were found to explain 12.4%–33.3% of the phenotypic variance. This genetic linkage map provides a tool for marker-assisted selection and for studies of genome in E. ulmoides.     

Genetic Variability and Population Structure of Salvia lachnostachys: Implications for Breeding and Conservation Programs

The genetic diversity and population structure of Salvia lachnostachys Benth were assessed. Inter Simple Sequence Repeat (ISSR) molecular markers were used to investigate the restricted distribution of S. lachnostachys in Parana State, Brazil. Leaves of 73 individuals representing three populations were collected. DNA was extracted and submitted to PCR-ISSR amplification with nine tested primers. Genetic diversity parameters were evaluated. Our analysis indicated 95.6% polymorphic loci (stress value 0.02) with a 0.79 average Simpson’s index. The Nei-Li distance dendrogram and principal component analysis largely recovered the geographical origin of each sample. Four major clusters were recognized representing each collected population. Nei’s gene diversity and Shannon’s information index were 0.25 and 0.40 respectively. As is typical for outcrossing herbs, the majority of genetic variation occurred at the population level (81.76%). A high gene flow (Nm = 2.48) was observed with a correspondingly low fixation index. These values were generally similar to previous studies on congeneric species. The results of principal coordinate analysis (PCA) and of arithmetic average (UPGMA) were consistent and all three populations appear distinct as in STRUCTURE analysis. In addition, this analysis indicated a majority intrapopulation genetic variation. Despite the human pressure on natural populations our study found high levels of genetic diversity for S. lachnostachys. This was the first molecular assessment for this endemic species with medicinal proprieties and the results can guide for subsequent bioprospection, breeding programs or conservation actions.     

Isolation and Characterization of Novel Microsatellite Markers in Pomegranate (Punica granatum L.)

Pomegranate (Punica granatum L.) has been cultivated from ancient times for its economic, ornamental and medicinal properties globally. Here, we report the isolation and characterization of 12 polymorphic microsatellite markers from a repeat-enriched genomic library of Punica granatum L. The genetic diversity of these loci was assessed in 60 genotypes of Punica granatum L. All loci were variable: the number of polymorphic alleles per locus ranged from two to five (average 2.9). The observed and expected heterozygosities ranged from 0.15 to 0.87 and 0.29 to 0.65, respectively. The polymorphic information content ranged from 0.26 to 0.61 (average: 0.43). To the best of our knowledge, this is the first time that polymorphic microsatellite markers have been reported for P. granatum L. These new markers should allow studies of the population structure and genetic diversity of pomegranate to be performed in the future.     

Assessment of Genetic Diversity of Bermudagrass (Cynodon dactylon) Using ISSR Markers

Bermudagrass (Cynodon spp.) is a major turfgrass for home lawns, public parks, golf courses and sport fields and is known to have originated in the Middle East. Morphological and physiological characteristics are not sufficient to differentiate some bermudagrass genotypes because the differences between them are often subtle and subjected to environmental influences. In this study, twenty seven bermudagrass accessions and introductions, mostly from different parts of Iran, were assayed by inter-simple sequence repeat (ISSR) markers to differentiate and explore their genetic relationships. Fourteen ISSR primers amplified 389 fragments of which 313 (80.5%) were polymorphic. The average polymorphism information content (PIC) was 0.328, which shows that the majority of primers are informative. Cluster analysis using the un-weighted paired group method with arithmetic average (UPGMA) method and Jaccard's similarity coefficient (r = 0.828) grouped the accessions into six main clusters according to some degree to geographical origin, their chromosome number and some morphological characteristics. It can be concluded that there exists a wide genetic base of bermudograss in Iran and that ISSR markers are effective in determining genetic diversity and relationships among them.     

Microsatellite Loci in the Gypsophyte Lepidium subulatum (Brassicaceae), and Transferability to Other Lepidieae

Polymorphic microsatellite markers were developed for the Ibero-North African, strict gypsophyte Lepidium subulatum to unravel the effects of habitat fragmentation in levels of genetic diversity, genetic structure and gene flow among its populations. Using 454 pyrosequencing 12 microsatellite loci including di- and tri-nucleotide repeats were characterized in L. subulatum. They amplified a total of 80 alleles (2–12 alleles per locus) in a sample of 35 individuals of L. subulatum, showing relatively high levels of genetic diversity, H_O = 0.645, H_E = 0.627. Cross-species transferability of all 12 loci was successful for the Iberian endemics Lepidium cardamines, Lepidium stylatum, and the widespread, Lepidium graminifolium and one species each of two related genera, Cardaria draba and Coronopus didymus. These microsatellite primers will be useful to investigate genetic diversity, population structure and to address conservation genetics in species of Lepidium.     

Isolation and Characterisation of 11 Polymorphic Microsatellite Markers in Papaver rhoeas L. (Corn Poppy), a Major Annual Plant Species from Cultivated Areas

Papaver rhoeas, an annual plant species in the Papaveraceae family, is part of the biodiversity of agricultural ecosystems and also a noxious agronomic weed. We developed microsatellite markers to study the genetic diversity of P. rhoeas, using an enriched microsatellite library coupled with 454 next-generation sequencing. A total of 13,825 sequences were obtained that yielded 1795 microsatellite loci. After discarding loci with less than six repeats of the microsatellite motif, automated primer design was successful for 598 loci. We tested 74 of these loci for amplification with a total of 97 primer pairs. Thirty loci passed our tests and were subsequently tested for polymorphism using 384 P. rhoeas plants originating from 12 populations from France. Of the 30 loci, 11 showed reliable polymorphism not affected by the presence of null alleles. The number of alleles and the expected heterozygosity ranged from 3 to 7.4 and from 0.27 to 0.73, respectively. A low but significant genetic differentiation among populations was observed (F_ST = 0.04; p < 0.001). The 11 validated polymorphic microsatellite markers developed in this work will be useful in studies of genetic diversity and population structure of P. rhoeas, assisting in designing management strategies for the control or the conservation of this species.     

Genome-Wide SNP Markers for Genotypic and Phenotypic Differentiation of Melon (Cucumis melo L.) Varieties Using Genotyping-by-Sequencing

Melon (Cucumis melo L.) is an economically important horticultural crop with abundant morphological and genetic variability. Complex genetic variations exist even among melon varieties and remain unclear to date. Therefore, unraveling the genetic variability among the three different melon varieties, muskmelon (C. melo subsp. melo), makuwa (C. melo L. var. makuwa), and cantaloupes (C. melo subsp. melo var. cantalupensis), could provide a basis for evolutionary research. In this study, we attempted a systematic approach with genotyping-by-sequencing (GBS)-derived single nucleotide polymorphisms (SNPs) to reveal the genetic structure and diversity, haplotype differences, and marker-based varieties differentiation. A total of 6406 GBS-derived SNPs were selected for the diversity analysis, in which the muskmelon varieties showed higher heterozygote SNPs. Linkage disequilibrium (LD) decay varied significantly among the three melon varieties, in which more rapid LD decay was observed in muskmelon (r² = 0.25) varieties. The Bayesian phylogenetic tree provided the intraspecific relationships among the three melon varieties that formed, as expected, individual clusters exhibiting the greatest genetic distance based on the posterior probability. The haplotype analysis also supported the phylogeny result by generating three major networks for 48 haplotypes. Further investigation for varieties discrimination allowed us to detect a total of 52 SNP markers that discriminated muskmelon from makuwa varieties, of which two SNPs were converted into cleaved amplified polymorphic sequence markers for practical use. In addition to these markers, the genome-wide association study identified two SNPs located in the genes on chromosome 6, which were significantly associated with the phenotypic traits of melon seed. This study demonstrated that a systematic approach using GBS-derived SNPs could serve to efficiently classify and manage the melon varieties in the genebank.     

Analysis of MTHFR and MTRR Gene Polymorphisms in Iranian Ventricular Septal Defect Subjects

Ventricular septal defect (VSD) is one of the most common types of congenital heart defects (CHD). There are vivid multifactorial causes for VSD in which both genetic and environmental risk factors are consequential in the development of CHD. Methionine synthase reductase (MTRR) and methylenetetrahydrofolate reductase (MTHFR) are two of the key regulatory enzymes involved in the metabolic pathway of homocysteine. Genes involved in homocysteine/folate metabolism may play an important role in CHDs. In this study; we determined the association of A66G and C524T polymorphisms of the MTRR gene and C677T polymorphism of the MTHFR gene in Iranian VSD subjects. A total of 123 children with VSDs and 125 healthy children were included in this study. Genomic DNA was extracted from the buccal cells of all the subjects. The restriction fragment length polymorphism polymerase chain reaction (PCR-RFLP) method was carried out to amplify the A66G and C524T polymorphism of MTRR and C677T polymorphism of MTHFR genes digested with Hinf1, Xho1 and Nde1 enzymes, respectively. The genotype frequencies of CC, CT and TT of MTRR gene among the studied cases were 43.1%, 40.7% and 16.3%, respectively, compared to 52.8%, 43.2% and 4.0%, respectively among the controls. For the MTRR A66G gene polymorphism, the genotypes frequencies of AA, AG and GG among the cases were 33.3%, 43.9% and 22.8%, respectively, while the frequencies were 49.6%, 42.4% and 8.0%, respectively, among control subjects. The frequencies for CC and CT genotypes of the MTHFR gene were 51.2% and 48.8%, respectively, in VSD patients compared to 56.8% and 43.2% respectively, in control subjects. Apart from MTHFR C677T polymorphism, significant differences were noticed (p < 0.05) in C524T and A66G polymorphisms of the MTRR gene between cases and control subjects.     

Structure and Principal Components Analyses Reveal an Intervarietal Fusion in Malaysian Mistletoe Fig (Ficus deltoidea Jack) Populations

Genetic structure and biodiversity of the medicinal plant Ficus deltoidea have rarely been scrutinized. To fill these lacunae, five varieties, consisting of 30 F. deltoidea accessions were collected across the country and studied on the basis of molecular and morphological data. Molecular analysis of the accessions was performed using nine Inter Simple Sequence Repeat (ISSR) markers, seven of which were detected as polymorphic markers. ISSR-based clustering generated four clusters supporting the geographical distribution of the accessions to some extent. The Jaccard’s similarity coefficient implied the existence of low diversity (0.50–0.75) in the studied population. STRUCTURE analysis showed a low differentiation among the sampling sites, while a moderate varietal differentiation was unveiled with two main populations of F. deltoidea. Our observations confirmed the occurrence of gene flow among the accessions; however, the highest degree of this genetic interference was related to the three accessions of FDDJ10, FDTT16 and FDKT25. These three accessions may be the genetic intervarietal fusion points of the plant’s population. Principal Components Analysis (PCA) relying on quantitative morphological characteristics resulted in two principal components with Eigenvalue >1 which made up 89.96% of the total variation. The cluster analysis performed by the eight quantitative characteristics led to grouping the accessions into four clusters with a Euclidean distance ranged between 0.06 and 1.10. Similarly, a four-cluster dendrogram was generated using qualitative traits. The qualitative characteristics were found to be more discriminating in the cluster and PCA analyses, while ISSRs were more informative on the evolution and genetic structure of the population.