ISSR and CDDP marker polymorphisms and discriminating features
The genetic diversity of the forty-eight Iris aucheri genotypes was assessed via ten ISSR primers and ten CDDP primers. The ISSR and CDDP primers generated 108 polymorphic bands (with an average of 10.8 bands per primer) and 134 polymorphic bands (with an average of 13.4 bands per primer), respectively (Table 3). For the ISSR primers, the total number of polymorphic bands (TPB) ranged from 7.00 (UBC-813, UBC-847, and UBC-849) to 15.00 (UBC-811), whereas the total number of polymorphic bands (TPB) varied from 7.00 (ERF1) to 21.00 (ABP1-1) for the CDDP primers.
Among the ISSR markers, UBC-813 and UBC-847 presented the highest number of observed alleles (Na) of 1.80, whereas UBC-834 presented the lowest value (1.49). The maximum values of the effective number of alleles (Ne), Shannon’s information index (I), expected heterozygosity or gene diversity (He), unbiased expected heterozygosity (uHe), and polymorphic information content (PIC) were recorded for UBC-823, with values of 1.57, 0.47, 0.32, 0.34, and 0.40, respectively. The lowest values of Ne, I, He, and uHe were found for UBC-812, with values of 1.30, 0.32, 0.20, and 0.21, respectively. The minimum PIC value was 0.22 for UBC-834. In the case of the CDDP markers, the highest values of Na, Ne I, He, and uHe were recorded with the Myb2 primer at 1.68, 1.47, 0.41, 0.28, and 0.29, respectively. The maximum PIC value was recorded at 0.37 with the KNOX3 primer. The lowest values of Na, Ne, I, He, and uHe were observed for the MADS-1 marker, with values of 1.26, 0.28, 0.17, and 0.18, respectively. The lowest PIC value was indicated by the Myb1 primer. The mean values of Na, Ne, I, He, uHe, and PIC for the ISSR markers were 1.71, 1.43, 0.39, 0.26, 0.27, and 0.32, whereas those of the CDDP markers were 1.53, 1.37, 0.34, 0.22, 0.23, and 0.26, respectively.
Cluster analysis using ISSR, CDDP, and ISSR+CDDP data.
The dendrogram generated via UPGMA demonstrated that the ISSR markers categorized all Iris aucheri genotypes into two primary clusters (Fig. 2A). The first cluster (in blue) contained six genotypes, and two subclusters were produced; each subcluster included three genotypes: G43, G44, and G45. G43, G44, and G45 produced the first subcluster, whereas the second subcluster comprised G42, G46, and G47. The second cluster (in red) had forty-two genotypes, and eight subclusters were generated: the first subcluster included only two genotypes (G36 and G41), and the second subcluster included five genotypes (G21, G34, G35, G38, and G39). The third subcluster included seven genotypes, whereas the fourth subcluster contained three genotypes. The fifth, sixth, seventh, and eighth subclusters comprised 6, 8, 3, and 8 genotypes, respectively. UPGMA was used to cluster the Iris aucheri genotypes on the basis of the CDDP data, resulting in two major clusters (Fig. 2B). The first cluster (in blue) included twenty-two genotypes, and four subclusters were generated, whereas the second cluster (in red) included twenty-six genotypes, with four subclusters. Two main clusters were generated from the combined ISSR and CDDP marker data (Fig. 2C). The first cluster (blue) included eight genotypes, and two subclusters were produced. The second cluster (red) comprised forty genotypes, and seven subclusters were generated.
Structural analysis of the Iris aucheri genotypes on the basis of marker datasets
Structural analysis was utilized to study the population stratification of the forty-eight Iris aucheri genotypes via ISSR markers. The number of genotype clusters was calculated in terms of the K value on the basis of genotypic data from the whole genome. The optimal K value was determined by graphing the number of clusters (K) versus K, with the largest peak occurring at K = 2 (Fig. 3A). The optimal K value suggested that two populations (Pop 1 in red and Pop 2 in green) presented the greatest possibility of population clustering (Fig. 3B). The criterion for membership probability was 0.80. Each subgroup of accessions was determined individually, with fractions less than 0.20 indicating admixture. The first population included forty-two I. aucheri genotypes, whereas the second population included only five I. aucheri genotypes: G16, G17, G18, G19, and G40. Only one genotype (G1) was identified as an admixed genotype between the two populations, which indicates that this genotype is not pure.
On the basis of the CDDP data, the optimal K value was determined via structure–harvester analysis, and the number of clusters (K) was plotted versus K, with a maximum peak occurring at K = 2. (Fig. 3C). The forty-eight genotypes of I. aucheri were divided into two genetic groupings (Fig. 3D). The first group (red) consisted of just four genotypes (G16, G17, G18, and G19), whereas the second group consisted of forty-one genotypes (green). All forty-five genotypes investigated were from pure populations, whereas the remaining three genotypes (G1, G40, and G43) were from an admixed population. The combined ISSR and CDDP genotyping data were used in STRUCTURE software to conduct population structure analysis on forty-eight I. aucheri genotypes via an admixture model. According to the Evano method, the optimal number of clusters is K=2 (Fig. 3E). A genotype was considered a pure member of a cluster if its chance of membership, in that cluster, was greater than 80%. Forty-six genotypes were classified as pure, whereas the remaining two genotypes (G1 and G40) were labeled admixed (Fig. 3F).
Analysis of molecular variance and diversity indices of Iris aucheri populations
Molecular variance analysis (AMOVA) was used to assess genetic differences within and among populations via ISSR, CDDP, and ISSR+CDDP marker data (Table 4). The significant PhiPT values for the ISSR markers, CDDP markers, and their combination were 0.10, with a p value of 0.001; 0.09, with a p value of 0.001; and 0.10, with a p value of 0.001, respectively. The ISSR, CDDP, and ISSR+CDDP marker data indicated that the variance among populations accounted for 10.41%, 9.36%, and 9.69%, respectively, of the total variation. However, the greatest variation occurred within the population, with 89.59, 90.64, and 90.31% for ISSR, CDDP, and their combined markers, respectively.
The diversity indices (Na, Ne, I, He, uHe, PPL, and NPB) of the I. aucheri populations are shown in Table 5. According to the ISSR data, Pop 5 had the greatest Na and PPL diversity indices of 1.81 and 88.89, respectively, whereas Pop 2 had the highest Ne, I, He, uHe, and NPB values of 1.47, 0.42, 0.28, 0.29, and 1.00, respectively. Pop 4 had the lowest values of these indices, with values of 1.42, 1.36, 0.33, 0.22, 0.23, 64.81, and 0.00. Among the CDDP markers, Pop 3 presented the greatest values of Na and PPL, with values of 1.69 and 80.60, respectively. Pop 4 presented the highest value of NPB (5.00 bands). Pop 2 presented the highest values for Ne, I, He, and UHe, i.e., 1.44, 0.39, 0.26, and 0.28, respectively. Pop 4 presented the lowest values of Na, Ne, I, He, uHe, and PPL, with values of 1.25, 1.29, 0.26, 0.17, 0.18, and 55.22, respectively, whereas Pop 5 presented the lowest value of NPB (0.00). On the basis of the combined ISSR and CDDP marker data, Pop 2 presented the greatest values of Ne, I, He, and uHe, with values of 1.45, 0.41, 0.27, and 0.28, respectively, whereas Pop 5, Pop 3, and Pop 4 presented the highest values of Na (1.69), PPL (82.57), and NBP (5.00), respectively. Pop 4 had the lowest Na, Ne, I, He, uHe, and PPL values, i.e., 1.32, 1.32, 0.29, 0.19, 0.20, and 59.34, respectively, whereas Pop 5 had the lowest NBP value (0.00). The mean values of Na, Ne, I, He, uHe, NPB, and PPL were 1.71, 1.43, 0.39, 0.26, 0.27, 82.00, and 0.40 for the ISSR markers; 1.52, 1.37, 0.34, 0.22, 0.23, 71.05, and 2.20 for the CDDP markers; and 1.61, 1.39, 0.36, 0.24, 0.25, 75.93, and 2.60 for the combined markers, respectively.