Data analysis

The DNA profiles were scored visually from gel photographs. The clear and reproducible amplified bands were chosen for the analyses. Polymorphism information content (PIC) value was estimated using the following equation described by Anderson et al. (1993): PIC=1-∑j=1nPij2, where Pij is the frequency of the ith allele for marker j and the summation extends over n amplicon calculated for each locus. The discrimination power was calculated by dividing the number of polymorphic markers amplified for each primer by the total number of polymorphic bands obtained (Khierallah et al. 2011). The presence of a band was designated as “1” and the absence of a band was recorded as “0”. The data obtained by scoring the RAPD and ISSR profiles, both individually as well as collectively, were subjected to the calculation of a similarity matrix using Jaccard’s coefficients. Cluster analysis was performed to construct dendrograms with the unweighted pair-group method by arithmetic averages (UPGMA) from the similarity data matrices using PAST3.11 software (Hammer et al. 2001). The co-phenetic correlation coefficient was used to check the goodness of fit a cluster analysis to the associated similarity matrix. A bootstrap analysis of 1000 replicates was performed using PAST3.11 software to estimate structural stability of clusters. Mantel test (Mantel 1967) was also performed using the PAST3.11 software in order to investigate the relationship between the genetic and phytochemical distances of the cultivars.