Genetic diversity and phylogenetic analyses are important tools for plant breeding and genetic research because they provide information that forms the basis for conservation and utilization of genetic resources in crop improvement. Moreover, after the creation of new genetic variants by mutation induction in crop plant, genetic diversity of the mutant lines needs to be evaluated using genetic based tools and advanced molecular techniques    , to facilitate their utility in cultivar development. Genetic diversity study is also important in crop breeding for the purpose of selecting genetically distant lines that are suitable for the production of heterotic hybrids  as well as for characterization and conservation of mutant lines.
Characterization of mutants with novel morpho-agronomic traits is fundamental to their utility in crop improvement program. Traditional characterization and selection technique depends on the application of morphological markers to reveal phenotypic variations in the plant population. However, the use of morphological markers is constrained by certain factors. It is influenced by the environments   , laborious and time consuming  , while some morphological markers have epistatic effects  . DNA molecular markers are genetic markers based on individual nucleotide sequence variation, which are the direct reflection of genetic polymorphisms at the DNA level  . PCR-based DNA markers are efficient tools used in plant breeding and genetics to estimate very precisely the genetic diversity  and differentiate among genotypes at species and sub-species level. Therefore, molecular marker-base characterization has been found very useful complement to phenotypic characterization of new mutants for the purpose of crop improvement and genetic studies. Kolade et al.  characterized some mutant lines of cowpea using random amplified polymorphic DNAs (RAPDs) and amplified fragment length polymorphism (AFLP) markers. Their study proved the discriminatory power of AFLP over RAPDs in characterizing the mutants.
Among the PCR-based DNA molecular markers, microsatellites or simple sequence repeats (SSR) are the most frequently used markers in the genetic diversity analysis of cowpea   -  . SSRs are DNA sequences with repeat lengths of a few base pairs. They provide an effective means for discriminating between genotypes   . Research findings had revealed that SSRs are more efficient than RFLPs, RAPDs and AFLPs in level of polymorphism and information content in genetic diversity assessment of legumes    . Microsatellites are one of the most useful molecular markers for genetic analyses due to their co-dominant inheritance, high variability and abundance, and ease of analysis    .
DNA sequence data have received a great deal of attention as a potential source of “phylogenetically informative” characters that are putatively less ambiguous than non-molecular characters  . Moreover, in genetic characterization of mutant lines, sequences analysis will reveal further details of the type of changes that had occurred at the nucleotides level and the precise location(s) where the mutations occurred in the plants genomes. Universal primers could be the best option when it is difficult to obtain monomorphic bands of samples’ DNA required for sequence analysis from other PCR-based primers. A universal primer designed from rbcL gene of chloroplast DNA (cpDNA), which encodes the large subunit of ribulose-1,5-bisphosphate carboxylase/oxygenase (RUBISCO), has been reported to be useful for the study of phylogenetic relationship of flowering plants at the species and generic level      . Moreover, universal primers for the amplification and sequencing of noncoding spacer between the atpB and rbcL genes of cpDNA has been reported  which are useful for phylogenetic assessment and species level of population studies. It has also been demonstrated that rbcL primers are useful for inter- and intraspecific evolutionary study of plants  .Therefore this study was conducted to reveal and characterize the cpDNA level of variations in selected cowpea mutants and their parents.
2. Materials and Methods
2.1. Plant Materials
The 32 cowpea genotypes used in this study are listed in Table 1. Twenty three out of the whole 32 plants were selected mutant lines induced by gamma irradiated cowpea seeds   , while only one line was a mutant induced by UVirradiated pollen of cowpea  .
2.2. Sample Preparation, Genomic DNA Extraction and Quantification
Seeds of cowpea lines were planted in plastic pots filled with top soil at the roof-top garden of CPEB, University of Ibadan. Young leaf samples of these lines were harvested at four weeks after planting. Harvested leaf samples were stored in freezer at −80˚C and later lyophilized for 72 hours to remove the moisture content from the leaf and present the leaf in a dried form for easy grinding. Steel balls were placed in the extraction tubes and samples were punched to bits into each extraction tube. The samples were ground in a 2000 geno-grinder machine for 2 minutes at 150 strokes per minute. Genomic DNA was extracted from the thirty two leaf samples using modified method of Cetyltrimethyl ammonium bromide (CTAB) and sodium dodecyl sulphate (SDS) extraction protocol in a mini prep format  . Extraction buffer was added to the ground sample,
Table 1. Cowpea genotypes used for molecular characterization.
U.I. = University of Ibadan, IITA = International Institute of Tropical Agriculture, SP = Spontaneous Mutation.
vortexed and placed on ice for 30 min. 20% SDS was added to each tube and were incubated in water bath at 65˚C for 10 min with continuous agitation. 5 M sodium chloride and CTAB buffer was added to each sample and incubated at 65˚C for 10 min. The samples were removed from the water bath and 400 µL (24:1) chloroform:isoamyl alcohol added. The sample was then centrifuged at 12,000 rpm for 10 min and the supernatant transferred into new 1.2 mL extraction tube and equal volume of isopropanol was added to it. This was stored at −20˚C for 1 hour for DNA precipitation. Samples were removed and centrifuged at 12,000 rpm for 10 min. The supernatant was decanted and the DNA pellet was washed with 70% cold ethanol and air dried. Pellet was re-suspended in 100 µl distilled water and 1.7 µl RNAse was added. DNA concentration and purity was determined using Nanodrop spectrophotometer at absorbance values of 260 nm and 280 nm. DNA quality was checked on 1 % agarose gel electrophoresis.
2.3. Simple Sequence Repeats (SSR) Genotyping
DNA sample from the 28 lines and 4 parental lines was genotyped using SSR markers. Sixteen VignaSSR markers were run on the 32 samples and only markers that showed polymorphic amplification among the lines were used (Table 2). Markers that were monomorphic among lines were not informative hence were not used.
PCR was carried out in a total volume of 15 µl containing 2 µl of 20 ng genomic DNA, 1.5 µl of 10X PCR buffer, 1 µl of 25 mM MgCl2, 1 µl of 2.5 mM dNTPs, 0.2 µl Taq polymerase (Inaqaba), 1 µl of tween 20, 1µl each of forward and reverse primer and 7.3 µl of sterile distilled water. Amplification conditions were: an initial denaturation step of 2 mins at 94˚C, followed by 28 cycles each consisting of a denaturation step of 1 min at 94˚C, annealing step of 1 min, and a final extension step of 72˚C for 5 mins. All amplification reactions were performed using PTC-200 Peltier thermal cycler (MJ Research Inc., Watertown, MA).
2.4. Electrophoresis of PCR Products
Mini acrylamide gel electrophoresis was used for better resolution following the standard procedures   . The plates were briefly treated by wiping both long and short plates with ethanol. Amplicons (amplified DNA fragments) were run on 6% acrylamide gel electrophoresis (Acrylamide, 10XTBE, 10% Ammonium persulfate and Temmed). The gel was allowed to polymerize and run with 0.5XTBE buffer (45 mMTris-acetate, 5 mM Boric acid, and 1 mM EDTA, pH
Table 2. List of SSR markers used in the study.
8.0) at 100 V for 1 hour. 100 bp ladder was used as a molecular size standard. Gels were visualized by staining with ethidium bromide solution (0.5 μg/ml) and banding patterns were photographed over UV light using UVP-computerized gel photo documentation system.
2.5. Sequencing Reaction and Analysis
The SSR markers used in the diversity study were all polymorphic and could not be used for sequencing. Therefore, the following pair of rbcL primer was used for the sequence:
H535―5’CTTTCCAAGGCCCGCCTCA3’ for forward sequence;
C705―5’CATCATCTTTGGTAAAATCAAGTCCA3’ for reverse sequence.
2.6. Cocktail Mix and Condition for the PCR
Genomic DNA was subjected to the following cocktail mix: 1.0 µL of 10× PCR buffer, 1.0 µL of 25 mM MgCl2, 0.5 µL of 5 pMol forward primer (H535), 0.5 µL of 5 pMol reverse primer (C705), 1.0 µL DMSO (Dimethyl sulfoxide), 0.8 µL of 2.5 mM DNTPs, 0.1 µL of 5 ug/µL Taq polymerase (Thermoscientific), 2.0 µL of 10 ng/µL of genomic DNA and 3.1 µL of water. PCR was carried out in a total volume of 10 µL containing 20 ng of genomic DNA. Amplification reactions were performed with Veriti 96 well thermal cycler (Applied Biosystems) using the following conditions: an initial denaturation step of 5 min at 94˚C, followed by 36 cycles each consisting of a denaturation step of 30 sec at 94˚C, annealing step of 30 sec at 56˚C, extension temperature of 72˚C for 45 sec, a final extension step of 72˚C for 7 min and a hold temperature of 10˚C. The PCR product (amplicon) was loaded on 1.5% Agarose to check the amplification, after which the amplicon was purified.
2.7. PCR Purification
The PCR product was purified by adding 2 volume (20 µL) of absolute ethanol. It was incubated at room temperature for 15 minutes, spun down at 10,000 rpm for 15 minutes and the supernatant was decanted. Then 100 µL of 70% ethanol was added, the mixture was votexed and spun down at 10,000 revolutions per minute (rpm) for 15 minutes. The supernatant was decanted and PCR product was air dried. 20 µL of DNAse and RNAse free H2O was added and the product was checked again on 1.5% agarose.
2.8. Sequencing Reaction
The reactions for 96-well reaction plates (microcentrifuge tubes) was prepared by adding the following reagents to a separate tube: 4.0 µL of 2.5X terminator ready reaction mix, 2.0 µL of 5X Bigdye Sequencing Buffer (BigDye Terminator v3.1), 3.2 pmol primer, 20 ng PCR product and deionized water was added to make a total volume of 20 µL. It was mixed well and spun down briefly. The sample was loaded on the Veriti 96 well thermal cycler (Applied Biosystems). The tubes were placed in a thermal cycler and set to the correct volume. Initial denaturation step of rapid thermal ramp to 96˚C for 1 min was performed. This was followed by 25 cycles each consisting of a denaturation step of 10 sec at 96˚C, annealing step of 5 sec at 50˚C, extension temperature of 60˚C for 4 min, a hold temperature of 4˚C and the contents of the tubes in a microcentrifuge were spin down.
2.9. Purification of Sequencing Product
To the PCR product in the 96-well reaction plate, 5 µL 125 mM EDTA was added. This was followed by the addition of 60 µL of 100% ethanol to each well to each well. The plate was sealed with plate septa, mixed by inverting 4 times and incubated at room temperature for 15 minutes. The plate was spin down at 3000 rpm for 45 minutes at 4˚C. The plate was inverted and it was spin up to 900 rpm. To each well, 60 µL of 70% ethanol was added and centrifuge at 3000 rpm for 15 min. The plate was inverted and it was spin up to 900 rpm for 1 minute and the samples were re-suspended in injection buffer i.e. HID Formamide (Applied Biosystems).
2.10. Data Scoring and Statistical Analysis
Microsatellite markers were scored as follows:
➢ 1 for present alleles.
➢ 0 for absent alleles.
➢ 9 for missing data.
The polymorphic information content (PIC); that is the level of polymorphism shown by each SSR marker for distinguishing cowpea lines, was determined following the procedures of Powell et al.  and Weir  :
where, He is the expected heterozygosity or PIC and Pi is the frequency of the ith allele.
Each SSR fragment was treated as binary matrix in which band presence was coded as 1 for present and 0 for absent. Based on the binary matrix, Euclidean dissimilarity index was computed. Subsequently, using neighbor joining clustering algorithm, a dendogram was generated with the unweighted pair group method using arithmetic average (UPGMA) algorithm of DARwin 5.0.158 software   .
The samples were loaded on 3130xl Genetic Analyzer in order to generate sequence data. DNA sequences were edited and analyzed using BioEdit and MEGA software. Phylogenetic analyses were conducted in MEGA5  .
3.1. Genetic Diversity Assessment
Four primers (Vm34, Vm37, Vm54 and Vm57) out of the sixteen SSR primers screened for the molecular study revealed polymorphic loci (Figure 1) while
Figure 1. The picture from 6% acryl amide gel showing polymorphic bands of amplified primers; (a) polymorphic bands from primer Vm34, (b) polymorphic bands from primer Vm57.
monomorphism were observed with the other twelve primers. The polymorphic markers were used to evaluate the genetic diversity of 32 cowpea lines. A total of 15 alleles were produced on these four loci with an average of 3.75 alleles per SSR locus (Table 3). The primer, Vm57 revealed the highest diversity index among all the lines and the polymorphic information content (PIC) varied from 0.33 to 0.63 with a mean of 0.51.
The dendrogram generated by UPGMA method showed the genetic relationship among cowpea lines (Figure 2). Seven distinct branches were revealed in the phylogenetic analysis of all the mutant lines including their parents. Comparing the clustering of the mutants with their parents, only IB-ER and IB-ER-2 were observed in the same cluster (Group I) together with their parent IB, while other mutant derivatives of IB clustered in group VI and VII. IT90K-284-2 clustered with its mutant lines IT90K-BS-4, IT90K-500EM, IT90K-284TRV, IT90K-BS-3 and IT90K-284FPL-2 in group I, while it’s other mutant derivatives clustered in group II and IV. IT86D-719 was observed in cluster IV, while its mutant derivatives clustered with group I, III, V and VI. IT89KD-374-57 was observed in cluster IV, while other mutant derived from it clustered with group I and II. These results show that the mutants were diverse in their genetic makeup from their parental lines.
3.2. Sequence Analysis
The sequence alignments of cowpea mutant lines and their parents with the rbcL reference sequence are presented in Figure 3. As a result of mutagenic treatment of cowpea lines with gamma rays, insertion of T and G was observed in the rbcL
Table 3. Summary statistics of SSR markers used in cowpea mutants’ characterization.
Figure 2. Dendrogram showing genetic diversity between the parents and mutant cowpea lines.
regions 2 and 3, respectively of IB-ER when compared with its parent IB, while G was substituted for A at region 4 of the sequence. A was deleted in IB-ER at position 129 of its sequence when compared with IB sequence. In IB-LT, insertion of AT was observed at positions 1 - 2, while base substitutions of T for G at region 3 and G for A at position 4 of its parent IB-CR sequence. In addition to these, deletion of 5 bases (AATTC) in IB-LT was observed at position 128 to 132 when compared with IB-CR. The sequence data also revealed the type of changes that occurred in the rbcL region of IB line as a result of spontaneous mutations that produced IB-BPC, IB-CR, IB-Y-1 and IB-Y-2. There was deletion of G and A at regions 4 and 129, respectively in IB-BPC. In IB-CR insertion of T and TTC
Figure 3. rbcL sequence alignment of cowpea mutants with their parents.
was observed at positions 3 and 130 - 133, respectively. Deletion of A was observed at position 129 in IB-Y-1, while in IB-Y-2 there was insertion of G at position 3 and deletion of A at position 129.
From the rbcL sequence data there were variations in the effects of mutagenic treatment of IT86D-719 with gamma rays on the resultant mutants in this study. In mutant line IT-719BN-1, G and A nucleotides were deleted at positions 3. There was a deletion of G at position 3, while ACC bases were substituted for CAA at region 4 - 6 in IT-719BN-2. In IT-719FPL and IT-719G100DW, C was substituted for A at region 6 of the sequences. At regions 3 - 4 of IT-719G200BT, GA was deleted. In IT-719G400MS, GA was deleted at region 3 - 4, substitutions of C for G, G for C, C for A and C for G were respectively, observed at regions 5, 13, 19 and 21 and insertions of A was revealed at regions 28. However, in IT-719Y, G was deleted at region 3, while A and C were respectively, substituted for G and A at regions 4 and 6.
The mutants derived from IT89KD-374-57 were different from their parent with respect to rbcL sequence. There was an insertion of C at position 5 of IT89KD-400HT. In IT89KD-400UF, the sequence data revealed insertions of GAT, T and G at positions 3 - 5, 12 and 28 respectively, while A, G, C, C and C were respectively, substituted for C, C, T, A and G at positions 6, 7, 15, 19 and 21. Insertion of four bases (GGAC) was observed at position 2 - 5 of IT86KD-NL.
The rbcL sequence data shows various changes that occurred in the mutants produced from gamma irradiation of IT90K-284-2 line. In IT90K-284TRV, there was insertion of A at positions 4. In IT90K-500EM, there were insertions of GA and T at positions 3 - 4 and 6 respectively. In IT90K-BS-1, insertion of TTA bases was observed at position 2 - 4, while there was deletion of C and insertion of A in IT90K-BS-2 at position 5. TG was respectively, inserted to region 3 - 4 in IT90K-BS-3, while there was only an insertion of G to region 4 in IT90K-BS-4. IT90K-UVFPL-REV was a four-primary leaf mutant derived from the UV treatment of pollen prior to pollination in IT90K-284-2. In this mutant, a base substitution of C for T was observed at region 20. When compared with other regions of the rbcL sequences, high mutation rates were observed between regions 2 - 6 of all samples studied with regions 3 - 4 being the highest mutable sites.
Analysis of the rbcL sequence of the mutants shows the presence of InDels and point mutations (base substitutions) as the two main classes of mutations induced in the plastid DNA of the mutants studied. Only one stable UV induced mutant plant was included in the sequence analysis, hence the type of mutation induced from the UV source could not be quantified by the sequence results. Also, the rates of different transitional substitutions are shown in bold and those of transversional substitutions are shown in italics (Table 4). The nucleotide frequencies were 26.95% (A), 34.43% (T), 24.09% (C), and 14.53% (G). A total of 61.38% AT rich region was identified, while GC rich region was found to be 38.62%. The analysis involved 30 nucleotide sequences. Codon positions included were 1st + 2nd +3rd + Noncoding. All ambiguous positions were removed for each sequence pair. There were a total of 134 positions in the final dataset. Observations on the positions towards the end of the sequence alignments were not accounted for in this report because they could not be justified.
3.3. Molecular Phylogenetic Analysis of the Mutants
The estimated genetic divergence among the sequences of cowpea were obtained using rbcL gene (Table S1). Phylogenetic reconstruction of the mutant lines including the parents was inferred using the Maximum Likelihood method based on the Kimura 2-parameter model  . The bootstrap consensus tree inferred from 1000 replicates is taken to represent the evolutionary history of the taxa analyzed  . Branches corresponding to partitions reproduced in less than 50% bootstrap replicates are collapsed. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown above the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach, and then selecting the topology with superior log likelihood value. The rate variation model allowed for some sites to be evolutionarily invariable ([+I], 0.0000% sites). Codon positions included were 1st +2nd +3rd +Noncoding. All positions with less than 95% site coverage were eliminated. That is, fewer than 5% alignment gaps, missing data, and ambiguous bases were allowed at any
Table 4. Maximum composite likelihood estimate of the pattern of nucleotide substitution.
NOTE: Each entry shows the probability of substitution (r) from one base (row) to another base (column). For simplicity, the sum of r values is made equal to 100.
position. The evolutionary history was inferred using the UPGMA method  . The optimal tree is shown in Figure 4. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown above the branches  . The phylogenetic distances were computed using the Kimura 2-parameter method  and are in the units of the number of base substitutions per site. The analysis involved 32 nucleotide sequences. Codon positions included were 1st +2nd +3rd +Noncoding. All ambiguous positions were removed for each sequence pair. There were a total of 1164 positions in the final dataset.
Four main cluster groups of cowpea lines were revealed from the phylogenetic tree (Figure 4). The number of genotypes grouped within cluster I, II, III and IV was 1, 23, 4 and 2, respectively (excluding 2 samples with bad sequence data). The grouping which was irrespective of the mutant origin indicated that similarity within the mutant populations was independent of the rbcL sequence data. The rbcL reference sequence retrieved from NCBI database was included in the computation and it was found within cluster II, while the GU140278_COI gene sequence from a different organism that is not related to cowpea was usedas an out-group which confirmed the validity of these results. The result indicates that IB is similar to its mutant derivatives except IB-CR which was differentiated into cluster I. IT86D-719 is different from all its mutant lines except IT-719G200BT, IT-719BN-1 and IT-719G400MS. IT89KD-374-57 is similar to all its mutant lines except IT89KD-400UF, while IT90K-284-2 is different from IT90K-500EM only.
Molecular characterization has gradually replaced phenotypic evaluation among plants and animals species since the evolution of gene theory. PCR-based DNA markers are very valuable tools to plant geneticists and breeders because they are useful for precise estimates of genetic diversity. The SSR markers used in this study were able to separate each of the parental lines from most of their mutant derivatives. Some DNA changes resulting from mutation events led to the variations detected by microsatellites between the mutants and their parents in this study. Hamzekhanlu et al.  had reported that polymorphic amplification
Figure 4. A Phylogeny of rbcL gene in cowpea mutants inferred by using the Maximum Likelihood method based on the Kimura 2-parameter model and the percentage of trees in which the associated taxa clustered together is shown above the branches.
products which represent one allele per locus can result from changes in either the sequence of the primer binding site, such as point mutations, or from changes altering the size or preventing successful amplification of a target DNA such as insertions, deletions, and inversions. Mutations resulting in polymorphisms are those occurring on primer binding sites, leading to an increase or decrease in the total number of primer binding sites, and consequently the number of amplified fragments  . Furthermore, some changes in the genome of the mutants were responsible for the differences observed in the clustering of most mutants from their parents in the dendrogram constructed from SSR markers. Although the clustering patterns of these mutants are not correlated with their phenotypic traits (morphology), further investigations may be needed to know if they could infer useful agronomic traits inherent in the mutants. The genetic diversity assessment presented in this study demonstrates that SSR marker could be used for molecular characterization of cowpea mutant lines. Selected cowpea mutants in this study were genetically diverse from their parents indicating that they are valuable genetic resources to be added to cowpea genepool.
The effects of any gene mutation on an organism will vary depending upon the functional region of the gene (promoter, intron, coding, or 3’ untranslated regions) where the mutation occurs, and whether the function of an essential protein has been altered  . The effect of mutations on the function of a gene will be determined by the type of changes that has occurred in the nucleotide sequence. Generally, mutations are classified based on the type and the extent of molecular changes in the nucleic acid affected by mutational event. Analysis of the rbcL sequence of the mutants shows the presence of insertions and deletions (InDels) and point mutations (base substitutions) as the two main classes of mutations induced in the plastid DNA of the mutants studied. Of all the mutations observed, InDels are more predominant mutations in this study. This is expected because most of the mutants studied were induced by gamma irradiation. Brown  listed some negative consequences of ionizing radiation overdoses such as deletions of DNA nucleotide sequences that may cause reading-frame shifts, inactive protein products, or faulty transcripts. These may have accounted for a wide range of morphological variations observed in among the mutants reported in Olasupo et al.  . The relatively higher rate of transitional substitution range (5.16 - 12.23) as compared to the rate of transversion range (4.83 - 11.45) among the mutants studied corroborates the report of Lee et al.  that transitions occur more frequently than transversions in mutational events. The higher percentage of AT rich region than GC rich region observed in this study indicates that there is a bias towards AT in the substitutional rate with regards to A, T, G and C. Yamane et al.  reported that the AT-richness of chloroplast DNA influences both the distribution of InDels and dynamics of nucleotide substitution. The highest rate of mutations observed in region 3 - 4 also indicates that the region is less conserved in cowpea rbcL gene.
The four groups obtained from the phylogenetic tree which is irrespective of the mutant origin indicated that similarity within the mutant populations is independent of the rbcL sequence data. The clustering together of parental lines and their respective mutant derivatives in the phylogenetic tree of rbcL gene of the cpDNA confirms each of the parents as the evolution origin of the mutant lines. Differentiation of IB-CR into cluster I and the mutant lines IT-719G400MS and IT89KD-400UF into cluster IV appears to correlate with deep green crinkled leaf morphology of IB-CR and deep green normal leaves of IT-719G400MS and IT89KD-400UF. These may suggest that the mutations had caused major changes in the rbcL gene of these cowpea lines. However, the sequence data could not separate the each of the parental lines from their mutant derivatives; this seems justified because each of the parents and their mutants are components within the single Vigna unguiculata entity. This observation indicates that the phenotypic traits of most of the cowpea mutant lines are controlled by nuclear genes that could not be differentiated by the rbcL marker. Nevertheless, rbcL primer has been used in this study to reveal genetic diversity which exists in the chloroplast genome of selected cowpea mutants and their parents. The findings also showed the changes that occurred at the nucleotide level to produce new mutant varieties of cowpea thus indicating new alleles being added to cowpea genepool.
Induced mutation is a valuable tool for creating new genetic variability to complement existing germplasm and broaden crop genetic base. Therefore, creating genetic diversity is needful for effective plant breeding programs. The present study proved that SSR markers are useful for the genetic diversity assessment and characterization of cowpea mutants. Most of the mutants studied were genetically diverse from their parents indicating that they are valuable genetic resources to be added to cowpea genepool. Although the sequence analysis of cpDNA could not reveal details of the mutations with respect to the mutants’ morphological characters, however it shows that rbcL markers may be used for phylogenetic and evolutionary studies of mutant lines. Therefore, there is the need for whole genome sequencing of these mutant lines in order to obtain detail information about the genes responsible for the mutants’ phenotypes and how they may enhance the genetic map and cultivar improvement of cowpea.
We thank Dr. Wole Oyatomi of the Genetic Resources Centre of International Institute of Tropical Agriculture for the SSR primers used in this study.
Table S1. Estimated genetic divergence among the sequences of cowpea obtained using rbcL gene.
 Bennici, A., Maria, A. and Giovanni, G.V. (2003) Genetic Stability and Uniformity of Foeniculum vulgare Mill. Regenerated Plants through Organogenesis and Somatic Embryogenesis. Plant Science, 161, 221-227.
 Hamzekhanlu, M., Izadi-Darbandi, A., Pirvali-Beiranvand, N., Teher-Halljian, M. and Majdabadi, A. (2011) Phenotypic and Molecular Analysis of M7 Generation of Soybean Mutant Lines through Random Amplified Polymorphic DNA (RAPD) Marker and Some Morphological Traits. African Journal of Agricultural Research, 6, 1779-1785.
 Shiran, B., Amirbakhtiar, N., Kiani, S., Mohammadi, S.H., Sayed-Tabatabaei, B.E. and Moradi, H. (2007) Molecular Characterization and Genetic Relationship among Almond Cultivars Assessed by RAPD and SSR Markers. Scientia Horticulturae, 111, 280-290.
 Khan, T., Reza, O.H., Khan, A., Hoque, S., Islam, S. and Khan, B. (2015) Genetic Diversity Analysis of Cowpea by RAPD Markers. International Journal of Innovation and Applied Studies, 10, 459-465.
 Meglic, V. and Staub, J.E. (1996) Inheritance and Linkage Relationships of Isozyme and Morphological Loci in Cucumber (Cucumis sativus L.). Theoretical and Applied Genetics, 92, 865-872.
 Udensi, O.U., Okon, E., Ikpeme, E.V., Onung, O.O. and Ogban, F.U, (2015) Assessing Genetic Diversity in Cowpea (Vigna unguiculata L. Walp.) Accessions Obtained from IITA, Nigeria Using Random Amplified Polymorphic DNA (RAPD) International Journal of Plant Breeding and Genetics, 10, 12-22.
 Brown, S.M., Hopkins, M.S., Mitchell, S.E., Senior, M.L., Wang, T.Y., Duncan, R.R., Gonzalez-Candelas, F. and Kresovich, S. (1996) Multiple Methods for the Identification of Polymorphic Simple Sequence Repeats (SSRs) in Sorghum [Sorghum bicolor (L.) Moench]. Theoretical and Applied Genetics, 93, 190-198.
 Kolade, O.A., Olowolafe, M.O. and Fawole, I. (2016) Characterization of Mutant Cowpea [Vigna unguiculata (L) Walp] Lines Using Random Amplified Polymorphic DNAs (RAPDs) and Amplified Fragment Length Polymorphism (AFLP) Markers. African Journal of Biotechnology, 15, 2530-2537.
 Ali, Z.B., Yao, K.N., Odeny, D.A., Kyalo, M., Skilton, R. and Eltahir, I.M. (2015) Assessing the Genetic Diversity of Cowpea [Vigna unguiculata (L.) Walp.] Accessions from Sudan Using Simple Sequence Repeat (SSR) Markers. African Journal of Plant Science, 9, 293-304.
 Asare, A.T., Gowda, B.S., Galyuon, K.A., Aboagye, L.L., Takrama, J.F. and Timko, M.P. (2010) Assessment of the Genetic Diversity in Cowpea (Vigna unguiculata L. Walp.) Germplasm from Ghana Using Simple Sequence Repeat Markers. Plant Genetic Resources, 8, 142-150.
 Badiane, F.A., Gowda, B.S., Cisse, N., Diouf, D., Sadio, O. and Timko, M.P. (2012) Genetic Relationship of Cowpea (Vigna unguiculata) Varieties from Senegal Based on SSR Markers. Genetics and Molecular Research, 11, 292-304.
 Lee, J.R., Back, H.J., Yoon, M.S., Park, S.K., Cho, Y.H. and Kim, C.Y. (2009) Analysis of Genetic Diversity of Cowpea Landraces from Korea Determined by Simple Sequence Repeats and Establishment of a Core Collection. Korean Journal of Breeding Science, 41, 369-376.
 Li, C.D., Fatokun, C.A., Ubi, B., Singh, B.B. and Scoles, G.J. (2001) Determining Genetic Similarities and Relationships among Cowpea Breeding Lines and Cultivars by Microsatellite Markers. Crop Science, 41, 189-197.
 Ogunkanmi, L.A., Ogundipe, O.T., Ng, N.Q. and Fatokun, C.A. (2008) Genetic Diversity in Wild Relatives of Cowpea (Vigna unguiculata) as Revealed by Simple Sequence Repeats (SSR) Markers. Journal of Food, Agriculture & Environment, 6, 263-268.
 Powell, W., Morgante, M., Andre, C., Hanafey, M., Vogel, J., Tingey, S. and Rafalski, A. (1996) Thecomparison of RFLP, RAPD, AFLP and SSR (Microsatellite) Markers for Germplasm Analysis. Molecular Breeding, 2, 225-238.
 Barkley, N.A., Dean, R.E., Pittman, R.N., Wang, M.L., Holbrook, C.C. and Pedrson, G.A. (2007) Genetic Diversity of Cultivated and Wild-Type Peanuts Evaluated with M13-Tailed SSR Markers and Sequencing. Genetic Research, 89, 93-106.
 Gioi, T.D., Boora, K.S. and Chaudhary, K. (2012) Identification and Characterization of SSR Markers Linked to Yellow Mosaic Resistance Gene(s) in Cowpea (Vigna unguiculata). International Journal of Plant Research, 2, 1-8.
 Gupta, S.K. and Gopalakrishna, T. (2010) Development of Unigene-Derived SSR Markers in Cowpea (Vigna unguiculata) and Their Transferability to Other Vigna Species. Genome, 53, 508-523.
 Kongjaimun, A., Kaga, A., Tomooka, N., Somta, P., Shimizu, T., Shu, Y., Isemura, T., Vaughan, D.A. and Srinives, P. (2012) An SSR-Based Linkage Map of Yardlong Bean (Vigna unguiculata (L.) Walp. Subsp. Unguiculata Sesquipedalis Group) and QTL Analysis of Pod Length. Genome, 55, 81-92.
 Chase, M.W., Soltis, D.E., Olmstead, R.G., Morgan, D., Les, D.H., Mishler, B.D., Duvall, M.R., Price, R.A., Hills, H.G., Qiu, Y.L., Kron, K.A., Rettig, J.H., Conti, E., Palmer, J.D., Manhart, J.R., Sytsma, K.J., Michaels, H.J., Kress, W.J., Karol, K.G., Clark, W.D., Hedren, M., Gaut, B.S., Jansen, R.K., Kim, K.J., Wimpee, C.F., Smith, J.F., Furnier, G.R., Strauss, S.H., Xiang, Q.Y., Plunkett, G.M., Soltis, P.S., Swensen, S.M., Williams, S.E., Gadek, P.A., Quinn, C.J., Eguiarte, L.E., Golenberg, E., Learn Jr,, G.H., Graham, S.W., Barrett, S.C., Dayanandan, S. and Albert, V.A. (1993) Phylogenetics of Seed Plants: An Analysis of Nucleotide Sequences from the Plastid Gene rbcL. Annals of the Missouri Botanical Garden, 80, 528-580.
 Palmer, J.D., Jansen, R.K., Michaels, H.J., Chase, M.W. and Manhart, J.R. (1988) Chloroplast DNA Variation and Plant Phylogeny. Annals of the Missouri Botanical Garden, 75, 1180-1206.
 Savolainen, D., Chase, M.W., Hoot, S.B., Morton, C.M., Soltis, D.E., Bayer, C., Fay, M.F. and de Bruijn, A.Y. (2000) Phylogenetics of Flowering Plants Based on Combined Analysis of Plastid atpB and rbcL Gene Sequences. Systematic Biology, 49, 306-316.
 Taberlet, P., Gielly, L., Pautou, G. and Boucet, J. (1991) Universal Primers for Amplification of Three Non-Coding Regions of Chlororplast DNA. Plant Molecular Biology, 17, 1105-1109.
 Chiang, T., Schaal, B.A. and Peng, C. (1998) Universal Primers for Amplification and Sequencing between the atpB and rbcL Genes of Chloroplast DNA. Botanical Bulletin of Academia Sinica, 39, 245-250.
 Olasupo, F.O., Ilori, C.O., Forster, B.P. and Bado, S. (2016) Mutagenic Effects of Gamma Radiation on Eight Accessions of Cowpea (Vigna unguiculata [L] Walp.).
American Journal of Plant Sciences, 7, 339-351.
 Olasupo, F.O., Ilori, C.O., Forster, B.P. and Bado, S. (2018) Selection for Novel Mutations Induced by Gamma Irradiation in Cowpea (Vigna unguiculata [L] Walp.).
International Journal of Plant Breeding and Genetics, 12, 1-12.
 Olasupo, F.O., Ilori, C.O. and Muyiwa, A.A. (2016) Radio-Sensitivity of Cowpea to Ultra-Violet Radiation by Pollen Treatment. Journal of Plant Breeding and Crop Science, 8, 339-351.
 Barril, P. and Nates, S. (2012) Introduction to Agarose and Polyacrylamide Gel Electrophoresis Matrices with Respect to Their Detection Sensitivities. In: Magdeldin, S., Ed., Gel Electrophoresis—Principles and Basics, InTech, London, 4-6.
 Perrier, X., Flori, A. and Bonnot, F. (2003) Methods for Data Analysis. In: Hamon, P., Seguin, M., Perrier, X. and Glazmann, J.C., Eds., Genetic Diversity of Cultivated Tropical Plants, Science Publisher Inc. and Cirad, Montpellier, 31-63.
 Perrier, X. and Jacquemoud-Collet, J.P. (2006) DARwin Software.
 Tamura, K., Peterson, D., Peterson, N., Stecher, G., Nei, M. and Kumar, S. (2011) MEGA5: Molecular Evolutionary Genetics Analysis Using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Molecular Biology and Evolution, 28, 2731-2739.
 Kimura, M. (1980) A Simple Method for Estimating Evolutionary Rate of Base Substitutions through Comparative Studies of Nucleotide Sequences. Journal of Molecular Evolution, 16, 111-120.
 Lee, S., Costanzo, S. and Jia, Y. (2012) The Structure and Regulation of Genes and Consequences of Genetic Mutations. In: Shu, Q.Y., Forster, B.P. and Nakagawa, H., Eds., Plant Mutation Breeding and Biotechnology. Plant Breeding and Genetics Section, Joint FAO/IAEA Division of Nuclear Techniques in Food and Agriculture International Atomic Energy Agency, Vienna, 31-45.
 Brown, N. (2013) Mutagenesis. Institution of Plant Breeding, Genetics and Genomics, University of Georgia, Athens.
 Yamane, K., Yano, K. and Kawahara, T. (2006) Pattern and Rate of Indel Evolution Inferred from Whole Chloroplast Intergenic Regions in Sugarcane, Maize and Rice. DNA Research, 13, 197-204.