Genetic diversity and structure of Capsicum annuum as revealed by start codon targeted and directed amplified minisatellite DNA markers

Background Identification of high resolving DNA-based markers is of paramount importance to unlock the potential of genetic diversity and selection of unique accessions of Capsicum annuum L., within Cross River and Ebonyi States of Nigeria, for breeding and conservation. Therefore, we comparatively explored the effectiveness of start codon targeted (SCoT) and directed amplified minisatellite DNA (DAMD) markers for diversity analysis of the accessions. Fifteen accessions were collected for DNA extraction and amplifications with the markers. Results Dendrograms from SCoT and DAMD categorized the accessions into five and three genetic groups, respectively, while the principal component analysis identified five genetic clusters, each from the markers. The average values of allele, gene diversity and polymorphic information content detected with SCoT and DAMD demonstrate that the two markers were effective and efficient, especially, SCoT in genetic diversity study of the accessions of pepper. Number of polymorphic loci (NPL) and percentage polymorphic loci (PPL) from SCoT (NPL = 64, PPL = 80.00–95.73%) and DAMD (NPL = 56, PPL = 53.33–86.67%) were high, but higher in SCoT markers. Other effective genetic parameters (effective number of alleles, Nei’s genetic diversity and Shannon’s information indices) identified with the two marker systems elucidated the allelic richness, rich genetic diversity within the populations and informative nature of the markers, especially SCoT. The intraspecific genetic diversity, interspecific genetic diversity, and coefficient of differentiation obtained with SCoT and DAMD further exposed the genetic structure with more genetic divergence within than among the populations of the accessions. Estimate of gene flow from the SCoT markers was 3.8375 and 0.6.2042 for the DAMD markers. The estimate of gene flow values from the markers indicated extensiveness with SCoT (Nm = 3.8375) and extremely extensive with DAMD (Nm = 6.2042) among the populations. Conclusion This study shows that SCoT markers may be more useful and informative than DAMD in measuring genetic diversity and differentiation of the accessions of the genus Capsicum. Genetic parameters obtained with SCoT showed that the accessions from Cross River were more genetically diverse than the ones from Ebonyi State. Therefore, SCoT may be a preferred marker in evaluating genetic diversity for improvement and conservation of this spicy crop, C. capsicum.


Background
Pepper (Capsicum spp.), belonging to family Solanaceae and sub-family Solanoideae is grown worldwide for its vital uses as vegetables, spices, ornaments, medicines, lachrymatories, vitamins A and C [1][2][3][4][5][6][7]. It was discovered as early as 6000 years ago and it is regarded as the first spice utilized by human beings [3,8]. The ancestral sources are Central and Southern America hinged in Bolivia, while Brazil is the center for genetic diversity of this important spice [9][10][11][12]. Thirty-eight species of Capsicum have been reported out of which Capsicum annuum L., C. frutescens L., C. Chinense Jacq., C. baccatum L. (aji), and C. pubescens Ruiz and Pay were considered to have been domesticated and maintained through more than five independent occurrences [9,13,14]. Some of these species were considered to have undergone a process of an independent domestication in at least two zones of the New World including Mesoamerica (C. annum and C. frutescens) and South America (C. chinese, C. baccatum and C. pubescens) [9,15]. The cultivated species (C. annuum, C. frutenscens, C. chinese, C. pubescens and C. baccatum var. pendulum) are of great economic importance and they possess good attributes that can be maintained through domestication [8,16,17]. However, continuous selection during domestication brought about emergence of novel lines with phenotypic traits including larger, non-pungent fruit with greater shape variation and tremendous increases in fruit mass [18,19]. Breeding programs targeting crop improvement are accelerated by the availability of well-characterized genetically diverse germplasm of pepper accessions and this is obtainable through the application of informative molecular markers to complement the roles of a morphologybased identification process.
Many marker systems of different categories including amplified fragment length polymorphisms (AFLP), simple sequence repeat (SSR), random amplified polymorphic DNA (RAPD), isozymes, inter-simple sequence repeat (ISSR), restriction fragment length polymorphism (RFLP), gene-based markers of start codon targeted and expressed sequence tag-simple sequence repeat (SCoT and EST-SSR), and single nucleotide polymorphism (SNP) have been utilized in studying and characterizing genetic loci linked to diversity, relationships, variations, parental selection, cultivar identity, phenotypic traits, purity, population studies, and resistance to diseases in Capsicum species . The basic information derivable from genetic variability of these potential markers could provide useful information that will assist plant breeders to manage and improve germplasm [23]. Use of microsatellite markers that are polymorphic, multi-allelic, reproducible and widely distributed in the pepper accessions could facilitate the selection of traits of interest and potential breeding materials for introgression through molecular marker-assisted breeding and germplasm conservation [23,43,44]. Functional gene-based marker such as SCoT, has been applied in other crops (Elymus sibiricus) [45]; (Boehmeria nivea L. Gaudich) [46]; and (Vigna unguiculata L.) [22]. Also another potential marker, directed amplified minisatellite DNA (DAMD), has also been useful in other crops including Flavoparmelia caperata [5]; Musa acuminata colla [47]; and Citrus [48].
Start codon targeted markers that are derivable from transcribed regions of the genome can be potentially useful for various applications in plant genotyping since they are capable of exposing polymorphisms that might be directly related to gene functions. These molecular markers were initially described, based on the observation that the short conserved regions of plant genes are surrounded by the ATG translation start codon [49,50]. The utility of SCoT was advocated due to their inherent features such as reproducibility, possibility of getting codominance during amplification, accuracy in investigation of genetic relatedness and genetic diversity of plants [50][51][52][53][54]. Also, DAMD markers, believed to be tandemly repeated units of genomes are satellite marker systems that have become valuable and of great importance in identification of genotypes following discovery [55]. They exhibit many variations as inversional mutations that result in their orientation and distribution on both strands of the DNA sequences [56,57]. Till now, SCoT and DAMD marker techniques have not been comparatively applied in the assessment of genetic diversity of pepper. Therefore, the objective of this study is to compare the efficacies of SCoT and DAMD markers in assessing genetic diversity among the accessions of C. annum from Cross River and Ebonyi States of Nigeria.

Plant materials
Fifteen accessions of C. annuum were collected from Cross River and Ebonyi States of Nigeria (Additional file 1: Table S1). One accession was purposefully obtained from each of the fifteen locations across Local Government Areas (LGAs) in Cross River and Ebonyi States to assess the efficiency of the two DNA-based molecular markers in studying the genetic diversity of C. annuum.

DNA extraction
DNA extraction of the accessions was carried out using a modified Cetyltrimethylammonium bromide (CTAB) method [58]. Briefly, a young leaf of C. annuum was collected and weighed between 100 and 200 mg. It was thoroughly ground with 500 μL of the extraction buffer (EB: 2% CTAB, 100 mM Tris-HCl, 20 mM EDTA, 1.4 M NaCl and 0.2% β-mercaptoethanol) and later made up with additional 200 μL of the EB before transferring to a sterile 1.5 ml microcentrifuge tube for homogenous mixture. It was incubated at 65°C for 15 min followed by addition of equal volume of phenol, chloroform and isoamyl alcohol in the ratio of 25:24:1. It was mixed and centrifuged at 14,000 rpm (rpm) for 15 min after which 500 μL of the supernatant was transferred to a new sterile 1.5 ml microcentrifuge tube. It was precipitated by adding 400 μL of cold-isopropanol, followed by gentle inversion prior to incubation at -20°C for overnight. It was centrifuged at 12,000 rpm to sediment the DNA prior to washing with 600 μL of 70% ethanol. It was subsequently air-dried and eluted with 100 μL of sterile water for polymerase chain reaction (PCR) set up.  Table S2) and 15.80 μL of 500 ml diethylpyrocarbonate (DEPC)-treated water (Invitrogen Corporation, USA). The PCR cycling profile used for the reaction consisted of an initial step at 95°C for 5 min., 40 cycles of 94°C for 30s, 55-65°C for 35 s, 72°C for 1 min, and a 10-min final extension at 72°C. Eight (8) μL of the PCR products were electrophoresed in a 1.5% agarose gel containing 0.5 mg/ml ethidium bromide and photographed on Transilluminator UV light (Fotodyne Incorporated, Analyst Express, USA).

Data analysis
Data matrix of SCoT and DAMD marker profiles for fragments of similar molecular weights from each amplicon or individual were scored as 1 (presence of alleles) and 0 (absence of alleles). The data achieved from the scoring of reproducible bands from SCoT and DAMD amplicons after the second repeated experiments were used for dendrogram reconstruction using Unweighted Pair Group Method with Arithmetic Mean (UPGMA) and dissimilarity index in Jaccard's option [60]. The analysis was conducted using Numerical Taxonomy and Multivariate Analysis System (NTSYSpc) software version 2.02. Furthermore, the genetic diversity, allele frequency and polymorphic information content (PIC) were computed using PowerMarker (Version 3.25). DARwin version 5 was used for principal component analysis (PCA), while effective number of alleles (Ne), Nei's genetic diversity (H), and Shannon's information index (I) [61] and population structure analyses including total genetic diversity (intraspecific diversity, Ht), gene diversity within the population (interspecific diversity, Hs), coefficient of differentiation (G ST ) and estimate of gene flow (Nm) of the accessions were done using POPGENE software version 1.32.

Results
Genetic diversity of accessions of Capsicum annuum as identified by start codon targeted markers To discern the level of genetic diversity in 15 accessions of C. annuum studied, a total of 10 SCoT primers were used for the study. Out of the 10 primers used, five (5) of them produced reproducible bands that were used in the analyses as exemplified in the gel images of SCoT13 and SCoT20 markers ( Figs. 1 and 2). A dendrogram of the 15 accessions using UPGMA procedure clustered the accessions into five major groups at a genetic distance threshold of 0.666 (Fig. 3). Group I was subdivided into subgroups I and II denoted as SGI and SGII, respectively. The SGI clustered CrPe-1, while SGII had CrPe-3, CrPe-5, CrPe-6, EbPe-5, CrPe-8 and EbPe-7. The accession EbPe-7 was the most genetically isolated one in the group followed by CrPe-8. The accession (CrPe) in SGI was collected from Cross River, while those in SGII were from both Cross River and Ebonyi (EbPe) States. Group II had CrPe-7, EbPe-6 and EbPe-3 and they were from Cross River and Ebonyi States. Group III clustered together accessions EbPe-1 and EbPe-2 from Ebonyi. Group IV had only CrPe-2 accession that was collected from Cross River State. Group VI contained CrPe-4 and EbPe-4 accessions that were collected from Cross River and Ebonyi States, respectively. Principal component analysis of the generated amplicons resulted into five clusters (Additional file 3: Fig. 1). The clustering of each is a representative of unique accession or accessions of C. annuum.
The five SCoT markers amplified a total of 57 alleles ( Table 1). The amplified alleles from each primer ranged from 7 to 14, with a mean of 12.000. Polymorphic information content values ranged from 0.7735-0.9193, with a mean value of 0.8709. The SCoT markers including SCoT13, SCoT28, SCoT20, SCoT24 and SCoT16 were found to be polymorphic. The genetic diversity ranged from 0.7911-0.9244, with a mean of 0.8815, while major allele frequency spanned from 0.1333-0.4000, with a mean of 0.2000. The identified allelic counts and frequencies from each of the SCoT loci ranged from 1 to 6 and 0.0667-0.4000, respectively (Additional file 4: Table S3).
The genetic diversity in CrPe-2 was identified to be highest among the accessions, with effective number of alleles (Ne), Nei's genetic diversity (H), and Shannon's information index (I) values of 1.9996, 0.4999 and 0.6931, respectively (Table 2). On the other hand, the genetic diversity in CrPe-4 was the lowest, with Ne, H, and I values of 1.0270, 0.0263 and 0.0708, respectively. The genetic diversity indices of these accessions were ranked in an ascending order as CrPe-2 > (CrPe-8, EbPe-   15 accessions using UPGMA procedure clustered them into three major groups at a genetic distance threshold of 0.6400 (Fig. 6). Group I was subdivided into subgroups I and II denoted as SGI and SGII, respectively. The SGI consisted of EbPe-1 and EbPe-5 accessions from Ebonyi State, while SGII was further subdivided into a (EbPe-3, CrPe-6, CrPe-8 and CrPe-4), b (EbPe-4, CrPe-5 and EbPe-7) and c (CrPe-1 and CrPe-3). Group II had EbPe-6 and CrPe-2, while group III clustered accessions EbPe-2 and CrPe-7. Principal component analysis of yielded five clusters (Additional file 6: Fig. 2). Each cluster contained an accession or accessions of C. annuum based on their dissimilarity indices within their respective coordinates.
The five DAMD primers amplified a total of 46 alleles ( Table 4). The amplified alleles from each DAMD ranged from 6 to 11, with a mean value of 10.1667. Polymorphic information content values ranged from 0.6519-0.8887 with an average of 0.8268. The DAMD primers including DAMD17R, DAMD13, DAMD16, DAMD1F and DAMD25 were polymorphic. The genetic diversity ranged from 0.6933-0.8978 with a mean value of 0.8430, while the major allele frequency spanned between 0.1333 and 0.4667, with a mean of 0.2556. Detected allelic counts and frequencies from each of the DAMD loci were respectively 1-7 and 0.0667-0.4667 (Additional file 7: Table S5). The genetic diversity in EbPe-2 was the highest, with Ne, H, and I values of 1.9501, 0.4872 and 0.6803, respectively (Table 5). Contrarily, the genetic diversity in EbPe-4 was the lowest, with Ne, H, and I values of 1.0000, 0.0000 and 0.0000, respectively. The genetic diversity values of these accessions were ranked as EbPe-2 > EbPe-6 > CrPe-7 > CrPe-2 > (CrPe-1, EbPe-3) > EbPe-1 > CrPe-3 > EbPe-5 > (CrPe-6, CrPe-8) > EbPe-7 from high to low. The mean values of Ne, H and I obtained from the accessions collected within Cross River State were 1.4015, 0.2693 and 0.4333, respectively, while those from Ebonyi had 1.4520, 0.2776 and 0.4265, respectively. The overall mean values of Ne, H and I detected in the accessions from the two States using DAMD were respectively 1.4268, 0.4081 and 0.6466.
The genetic variation in pepper accessions from Cross River, assessed using DAMD markers revealed that the mean values of Ht, Hs, G ST and Nm were 0.2693, 0.2465, 0.0992 and 6.4983, respectively ( Table 6) Table S6).

Comparison of SCoT and DAMD data obtained from the accessions of Capsicum annuum
Results obtained from the two markers were compared (Additional file 9: Table S7). Based on the dendrograms, five and three major groups were obtained from SCoT and DAMD data, respectively, while five clusters based on PCA were detected in each of the markers. A total of 57 alleles were found using SCoT markers, while 46 alleles were linked to DAMD markers.   Also, Ht recorded in the SCoT marker data was 0.3936, while DAMD had 0.2732. The value of Hs in SCoT was 0.3482 and DAMD yielded 0.2528. In SCoT marker data matrix, the G ST was 0.1153 (12% of the total genetic divergence among the populations and 88% within the populations) and 0.0746 (7% of the total genetic divergence among the populations and 93% within the populations) in DAMD, while the Nm obtained with the SCoT marker data was 3.8375 and 0.6.2042 for the DAMD markers.

Discussion
The efficacy of molecular marker systems for characterization and diversity assessment of C. annuum is very important for breeding, improvement and conservation purposes. The estimation of genetic diversity among accessions is beneficial for the utilization and maintenance of genetic resources. Genetic estimation helps to widen the genetic base of the cultivars and prevents gene insertion or deletion [62]. Assessing degree of genetic indices in populations of crops within germplasm using easily affordable and high-proof reading molecular markers could be a useful technique for characterization, detection and selection of novel accessions in a plant breeding-based approach for germplasm conservation of viable genetic resources with potential traits [63,64]. It will also aid in the verification of identified morphological traits for adequate discrimination of   [65]. Many molecular markers have been used to explore the genetic ingenuity of planting materials and to identify significant phenotypic characters associated with them. Molecular markers of different types have been applied in the characterization, mapping and study of genetic diversity of pepper [23][24][25][26][27][28][29][30][31][32][33]. In this study, SCoT and DMAD markers were explored to comparatively assess the genetic diversity of pepper accessions from Cross River and Ebonyi States of Nigeria. This became necessary since it is widely accepted that comparing different markers to study genetic diversity can be more resourceful and informative in the classification and selection of the best resolving markers [66].
Our results revealed that both SCoT and DAMD markers have high degree of discriminating potential on the pepper accessions as shown by the measured genetic diversity parameters (number of alleles per primer, PIC values, percent polymorphism, and polymorphic loci). Each of the markers, SCoT and DAMD resolved the accessions into five and three major groups, respectively using dendrogram and five clusters each based on PCA. The accessions were differently grouped by the two different sets of markers and this is possible due to variation in the accessibility of different genomic loci by the markers. Two major groups were identified in a study involving 149 black pepper germplasms using EST-SSR markers [33], while five major groups were detected in 30 accessions using SCoT markers [31] and this is in agreement with our findings with SCoT markers. Further, two major groups were found in a study involving 11 Chili pepper using AFLP markers [67]. The discrepancy between dendrograms and clusters of the pepper accessions by the two marker approaches could be linked to the different nature of the markers with regards to their regions of coverage in the genome, level of polymorphism and the number of loci involved as earlier reported [68,69]. This further demonstrates the significance of total number of identifiable loci and degree of coverage in the whole genome to actualize the required genetic relatedness within accessions [68].
On the part of genetic diversity assessment, SCoT markers identified those from Cross River State to be more genetically isolated compared to those from Ebonyi as a result of higher genetic indices inherent in them.  On the other hand, DAMD markers classified accessions from Ebonyi to be more diverse than those from Cross River State. It has been reported that H and I are major parameters in genetic diversity studies since H estimates genetic diversity within and between accessions, while I is for genetic diversity within and between populations regardless of the number of accessions accessed [69,70]. Based on the higher values in the basic parameters for genetic diversity assessment achieved from the markers (SCoT: H = 0.5936, I = 0.8590; DAMD: H = 0.4081, I = 0.6466), SCoT markers were found to be more efficient than DAMD and this further suggests that accessions from Cross River State are more genetically diverse than those obtained from Ebonyi State. In the present study, H and I, which are independent of number of samples being considered [70], are high in the two markers but higher in SCoT. These identified genetic parameters are higher than the ones obtained from a previous report, where H, I and PIC generated from 30 accessions of pepper using six SCoT markers were 0.2550, 0.4030 and 0.2120, respectively [31]. This difference could be linked to various marker techniques assessing different positions of the genome as previously reported in another combined marker-technique involving ISSR and RAPD [71]. The total potential numbers of identified alleles (SCoT: 57 with mean of 12.0000; DAMD: 46 with mean of 10.1667) in this study demonstrated high allelic richness and genetic diversity within the accessions thereby  contributing meaningfully to the genetic diversity among the accessions. The number of alleles detected with the SCoT markers is highly comparable to 53 alleles recently reported using six SCoT and six ISSR markers in 30 accessions of pepper [31] but a bit higher than the one identified with DAMD markers in the present study. In addition, the number of alleles generated with SCoT is lower but its allelic mean value is almost in agreement with a previous report involving 151 alleles, with a mean value of 10.0600 in 39 female Jojoba (Simmondsia chinensis (Link) Schneider) genotypes amplified with 15 SCoT markers were identified [72]. The obtained PIC using SCoT and DAMD markers were very high but higher in SCoT (PIC range: 0.7735-0.9193; mean: 0.8709) than DAMD (PIC range: 0.6519-0.8887; mean: 0.8268). In other genetic studies, molecular markers with PIC values between 0.4800-0.8000, with a mean of 0.7200 or PIC values above 50% were highly informative and polymorphic ones [73][74][75]. Both markers demonstrated PIC values higher than the previously reported ones in other crops of female Jojoba genotypes analyzed with SCoT markers (PIC range: 0.22-0.4800; mean: 0.4000) [72], Musa acuminata with DAMD (PIC range: 0.2900-0.4200; mean: 0.3500) [47] and grapevine cultivars with DAMD (PIC range: 0.3800-0.4900; mean: 0.4400) [76]. The observed mean values of PIC in DAMD is less than that of SCoT which shows that the SCoT markers evaluated are more informative than the DAMD ones. These markers exhibited high percentage polymorphisms (80-95.73% for SCoT and 53.33-86.67% for DAMD markers) in the accessions studied as previously reported on 30 accessions of pepper using SCoT (75-100%) and ISSR (70-100%) [31]. In addition, 33-100% polymorphism was detected from 40 landrace chickpea (Cicer arietinum L.) genotypes with SCoT and DAMD markers [77]. The percentage polymorphism obtained is generally high for both marker types, with the SCoT markers revealing higher polymorphisms than DAMD, thereby indicating the efficacious nature of the SCoT marker systems used. High number of alleles from the two marker polymorphic loci (SCoT: 64; DAMD: 56) demonstrated that they contributed meaningfully to the genetic diversity among these accessions of pepper. This is in contrast with previous reports where 160 alleles were identified with EST-SSR markers in 148 accessions of black pepper [33]. Similarly, 75 and 53 high numbers of alleles of detected from 64 and 30 accessions of pepper were respectively high [28,31]. The identified genetic parameters from the two markers (SCoT: Ne = 1.6971, H = 0.5936 and I = 0.8590; DAMD: Ne = 1.4268, H = 0.4081 and I = 0.6466) are quite high but higher in SCoT, thereby classifying the accessions from Cross River as the most genetically dissimilar compared to those collected from Ebonyi. In the present study, the identified values of these vital genetic indices are higher than the ones obtained from earlier studies on 30 accessions of pepper using SCoT (H = 0.2550, I = 0.4030) and ISSR (H = 0.3060, I = 0.4590) [31]. The obtained population indices in these accessions using SCoT (Ht = 0.3936, Hs = 0.3482, G ST = 0.1153 and Nm = 3.8375) and DAMD (Ht = 0.2732, Hs = 0.2528, G ST = 0.0746 and Nm = 6.2042) were higher than the ones previously reported in black pepper analyzed with EST-SSR [33]., in pepper accessions with EST-SSR and SCoT [31]., and in wild and domesticated populations of ramie amplified with SCoT [46]. The two marker techniques explicitly detected wide genetic diversity within and between the accessions evaluated, with the SCoT markers showing higher potential in discriminating them.
According to Nei [78], G ST is classified as low when its value is < 0.0500, medium when it is 0.0500 < G ST < 0.1500 and high when G ST > 0.1500. The obtained values of G ST in the two marker systems (SCoT = 0.1153; DAMD: =0.0746) indicate that only about 7 and 12% of the total genetic divergence exist and this may be due to source locations of the accessions (external factors), while 93 and 88% were contributed by genotypic variations within the populations possibly due to internal factors. This finding is consistent with the previous reports [79,80], where G ST values of 0.1700 and 0.3200 were respectively identified but lower than the ones obtained from chili pepper using AFLP (G ST = 0.8600) [67], Commiphora wightii (RAPD: G ST = 0.5800) [81], Torreya jackii Chun (ISSR: G ST = 0.6300) [82] and Curculigo latifolia (SCoT: G ST = 0.4800) [83]. The noted differences could be attributable to the number of samples analyzed and the nature of the plant genome assessed. The overall genetic variation (Ht = 0.2732) obtained in the whole population with the DAMD marker was lower than the one of intra-population genetic diversity (Hs = 0.2528). On the part of SCoT marker, 0.3936 of Ht was higher than Hs generated as 0.3482. It has been known that the Nm is low when Nm < 1, moderate (Nm > 1) and extensive (Nm > 4) in nature [84]. The recorded Nm in this present study indicates extensiveness in SCoT (Nm = 3.8375) and extremely extensive in DAMD (Nm = 6.2042) gene flow among the populations. This also exhibits a great deal of variation and rich genetic diversity within the populations from different locations across Cross River and Ebonyi States of Nigeria and these unique accessions could be utilized in further hybridization processes in breeding program. This high genetic differentiation within populations/accessions is possibly linked to outcrossing pollination process and exchange of genetic materials as reported in Begonia species involving microsatellite markers [85]. In population genetics, the estimation of gene flow with respect to the number of individual genes involves migration of one population to the other and per generation. It induces population differentiation and it is possible to prevent substantial differentiation in a population as a result of genetic drift if the gene flow is greater than 1 (Nm > 1) [86,87]. Compared to the outcome of this present study, lower values of G ST and Nm have been reported in other crop species including Rhamnus persicifolia Moris and Lilium cernuum with SSR markers [88,89].

Conclusion
Our research revealed that both SCoT and DAMD markers demonstrated a high discriminatory role in the dissection of genetic diversity within the accessions of C. annuum from both Ebonyi and Cross River States. However, SCoT markers were found more discriminatory and informative than DAMD markers following the higher genetic diversity indices inherent in SCoT, thereby classifying the accessions from Cross River State as more genetically diverse than those from Ebonyi. The SCoT markers generally were more efficacious by demonstrating higher parameters of all the genetic diversity indices deployed than DAMD markers, except major allele frequency, allele count, allele frequency and gene flow, which may not be good indicators when accessing genetic diversity in C. annuum species. Selection from the accessions of the unique clusters and genetic indicators as revealed by the polymorphic SCoT markers could serve as stocks in a hybridization process to achieve superior hybrids capable of withstanding genetic barriers, environmental stressors and novel gene exchange in breeding programs of C. annuum.