SCoT analysis
From 36 SCoT primers, 15 primers that can generate reproducible and high-resolution polymorphic fragments were selected for sample research (Table 2). A total of 365 credible bands were obtained, of which 355 were polymorphic, and the percentage of polymorphism was 97.2%. The bands amplified by each primer range from 19 (SCoT 11) to 32 (SCoT 7), with an average of 24.3 and the number of polymorphic bands ranged from 18 (SCoT 11) to 32 (SCoT 7), with an average of 23.7. The percentage of polymorphism ranged from 87.5% (SCoT 14) to 100% (SCoT 28), each primer has 18-32 polymorphic bands, and SCoT 7 produced the most polymorphic bands (PB = 32). PIC values ranged from 0.307 (SCoT 14) to 0.426 (SCoT 33), with an average of 0.379 (Table 2). The amplification results showed that SCoT markers could be used to detect the high polymorphism of Polygonatum germplasm.
Table 2. Sequence and polymorphism information of the 15 SCOT primers selected in this study
Primer code
|
Primer sequence (5′-3′)
|
No. of
amplified
bands
|
No. of
polymorphic
bands (PB)
|
Polymor-
phic (%)
|
Polymorphic information content (PIC)
|
SCoT1
|
CAACAATGGCTACCACCA
|
31
|
30
|
96.8
|
0.384
|
SCoT 2
SCoT 6
SCoT 7
SCoT 11
SCoT 13
SCoT 14
SCoT 15
SCoT 17
SCoT 20
SCoT 22
SCoT 28
SCoT 29
SCoT 31
SCoT 33
Mean
Total
|
CAACAATGGCTACCACCC
CAACAATGGCTACCACGC
CAACAATGGCTACCACGG
AAGCAATGGCTACCACCA
ACGACATGGCGACCATCG
ACGACATGGCGACCACGC
ACGACATGGCGACCGCGA
ACCATGGCTACCACCGAG
ACCATGGCTACCACCGCG
AACCATGGCTACCACCAC
CCATGGCTACCACCGCCA
CCATGGCTACCACCGGCC
CCATGGCTACCACCGCCT
CCATGGCTACCACCGCAG
-
-
|
22
23
32
19
25
24
27
21
22
21
20
25
28
25
24.3
365
|
22
22
32
18
24
21
27
21
21
20
20
25
27
25
23.7
355
|
100.0
95.7
100.0
94.7
96.0
87.5
100.0
100.0
95.5
95.2
100.0
100.0
96.4
100.0
97.2
-
|
0.413
0.386
0.380
0.363
0.352
0.307
0.409
0.341
0.383
0.341
0.410
0.406
0.377
0.426
0.379
-
|
Population genetic diversity
According to the provincial sources, 28 samples were divided into eleven populations, representing different provinces, table 3 shows their genetic diversity index. The values of Observed number of alleles (Na), Effective number of alleles (Ne), Nei's gene diversity (H) and Shannon's Information index (I) in Zhejiang Province were the highest, Na = 1.8630, Ne = 1.4561, H = 0.2796, I = 0.4275, and the standard deviations were 0.3443, 0.3129, 0.1601, 0.2194, respectively. The corresponding parameters of genetic diversity in Sichuan Province were the lowest, Na = 1.3315, Ne = 1.2344, H = 0.1373, I = 0.2005, and the standard deviations were 0.4714, 0.3333, 0.1953 and 0.2851, respectively. Among the 15 SCoT markers (Table 4), the spans of Ne, H, I values are 1.4057-1.7173, 0.2546-0.3947, 0.3959-0.5737, respectively and the standard deviations are 0.2466-0.3835, 0.1104-0.1834, 0.1302-0.2381 respectively.
According to Nei’s genetic identity and genetic distance, the highest genetic identity (0.9216) and the minimum genetic distance (0.0816) was recorded between Zhejiang province and Hunan province, while the lowest identity (0.5370) and the maximum genetic distance (0.6218) was recorded between Gansu province and Jiangxi province (Table 5). These results indicate that Polygonatum has high genetic diversity and SCoT markers can significantly distinguish the genetic differences of Polygonatum germplasm.
Table 3. Genetic diversity parameters for different provinces of Polygonatum
Province
|
Sample size
|
Na
|
Ne
|
H
|
I
|
Zhejiang
|
8
|
1.8630(0.3443)
|
1.4561(0.3129)
|
0.2796(0.1601)
|
0.4275(0.2194)
|
Hunan
Hubei
Yunnan
Guizhou
Sichuan
Gansu
Shanxi
Jiangxi
Fujian
Guangxi
|
4
4
2
3
2
1
1
1
1
1
|
1.6712(0.4704)
1.7315(0.4438)
1.3781(0.4856)
1.5973(0.4911)
1.3315(0.4714)
1.0000(0.0000)
1.0000(0.0000)
1.0000(0.0000)
1.0000(0.0000)
1.0000(0.0000)
|
1.4084(0.3634)
1.4552(0.3696)
1.2673(0.3434)
1.3822(0.3723)
1.2344(0.3333)
1.0000(0.0000)
1.0000(0.0000)
1.0000(0.0000)
1.0000(0.0000)
1.0000(0.0000)
|
0.2416(0.1915)
0.2260(0.1883)
0.1566(0.2011)
0.2244(0.1979)
0.1373(0.1953)
0.0000(0.0000)
0.0000(0.0000)
0.0000(0.0000)
0.0000(0.0000)
0.0000(0.0000)
|
0.3621(0.2732)
0.3974(0.2648)
0.2286(0.2936)
0.3340(0.2854)
0.2005(0.2851)
0.0000(0.0000)
0.0000(0.0000)
0.0000(0.0000)
0.0000(0.0000)
0.0000(0.0000)
|
Mean
|
-
|
1.5954
|
1.3673
|
0.2109
|
0.3250
|
Note: Standard deviations are in parentheses, Na Observed number of alleles, Ne Effective number of alleles, H Nei’s gene diversity, I Shannon’s Information index
Table 4. Analysis of genetic diversity of Polygonatum by SCoT markers
Maker
|
Ne
|
H
|
I
|
SCoT1
|
1.5237(0.2616)
|
0.3232(0.1255)
|
0.4936(0.1562)
|
SCoT2
SCoT6
SCoT7
SCoT11
SCoT13
SCoT14
SCoT15
SCoT17
SCoT20
SCoT22
SCoT28
SCoT29
SCoT31
SCoT33
Mean
|
1.6093(0.2975)
1.4946(0.2466)
1.4756(0.2683)
1.5820(0.3619)
1.5224(0.3298)
1.4057(0.3093)
1.5489(0.2575)
1.4914(0.3183)
1.5330(0.2892)
1.5803(0.3835)
1.7173(0.3115)
1.5684(0.2767)
1.5249(0.2955)
1.5971(0.2618)
1.4396
|
0.3562(0.1287)
0.3562(0.1701)
0.2992(0.1333)
0.3304(0.1752)
0.3111(0.1567)
0.2546(0.1638)
0.3357(0.1189)
0.3004(0.1450)
0.3238(0.1356)
0.3254(0.1834)
0.3947(0.1352)
0.3416(0.1250)
0.3195(0.1344)
0.3567(0.1104)
0.3286
|
0.5321(0.1535)
0.4811(0.1539)
0.4629(0.1705)
0.4893(0.2285)
0.4705(0.2055)
0.3959(0.2246)
0.5095(0.1460)
0.4639(0.1770)
0.4911(0.1749)
0.4813(0.2381)
0.5737(0.1617)
0.5156(0.1520)
0.4871(0.1691)
0.5353(0.1302)
0.4922
|
Note: Standard deviations are in parentheses, Na Observed number of alleles, Ne Effective number of alleles, H Nei’s gene diversity, I Shannon’s Information index
Table 5. Nei's genetic identity (above diagonal) and genetic distance (below diagonal) among eleven provinces of Polygonatum
POP ID
|
Zhejiang
|
Hunan
|
Hubei
|
Guizhou
|
Gansu
|
Shanxi
|
Jiangxi
|
Yunnan
|
Guangxi
|
Fujian
|
Sichuan
|
Zhejiang
|
****
|
0.9216
|
0.9185
|
0.9266
|
0.6892
|
0.7179
|
0.7809
|
0.8814
|
0.7765
|
0.7805
|
0.8947
|
Hunan
|
0.0816
|
****
|
0.8809
|
0.8679
|
0.6637
|
0.6797
|
0.7769
|
0.8501
|
0.7568
|
0.7668
|
0.8451
|
Hubei
|
0.0850
|
0.1268
|
****
|
0.8891
|
0.7064
|
0.7244
|
0.7365
|
0.8079
|
0.7524
|
0.7606
|
0.8381
|
Guizhou
|
0.0762
|
0.1416
|
0.1175
|
****
|
0.6881
|
0.6937
|
0.7295
|
0.8454
|
0.7334
|
0.7663
|
0.8644
|
Gansu
|
0.3722
|
0.4100
|
0.3475
|
0.3738
|
****
|
0.5973
|
0.5370
|
0.6538
|
0.5753
|
0.5726
|
0.6291
|
Shanxi
|
0.3315
|
0.3861
|
0.3224
|
0.3657
|
0.5154
|
****
|
0.6164
|
0.6364
|
0.6000
|
0.6137
|
0.6691
|
Jiangxi
|
0.2474
|
0.2524
|
0.3058
|
0.3154
|
0.6218
|
0.4838
|
****
|
0.7151
|
0.6822
|
0.6411
|
0.7126
|
Yunnan
|
0.1263
|
0.1624
|
0.2133
|
0.1679
|
0.4250
|
0.4520
|
0.3354
|
****
|
0.7238
|
0.7409
|
0.8369
|
Guangxi
|
0.2529
|
0.2786
|
0.2845
|
0.3101
|
0.5528
|
0.5108
|
0.3824
|
0.3232
|
****
|
0.7288
|
0.7303
|
Fujian
|
0.2478
|
0.2655
|
0.2736
|
0.2661
|
0.5576
|
0.4883
|
0.4446
|
0.2999
|
0.3164
|
****
|
0.7339
|
Sichuan
|
0.1112
|
0.1683
|
0.1767
|
0.1457
|
0.4635
|
0.4018
|
0.3389
|
0.178
|
0.3143
|
0.3093
|
****
|
Cluster analysis and Population structure
Cluster analysis of 28 germplasm materials was carried out by UPGMA program and a tree diagram was generated. On the whole, it is divided into two groups: Group I and Group II, when the similarity index is 0.606, these germplasm materials can be clustered into four subgroups (Fig. 2). Two-dimensional PCoA analysis of SCoT markers amplification products showed that Polygonatum germplasm was mainly divided into two groups, which was similar to the results shown by UPGMA dendrogram (Fig. 3). In addition, the correlation analysis also shows that the similarity matrix fits well and the similarity coefficient is 0.7335, indicating that the two analysis methods (PCoA and UPGMA) have obtained comparable results (Fig. 4).
STRUCTURE software ver.2.3.4 was used to estimate the population structure of 28 accessions based on fragments obtained from 15 SCoT primers. Structure Harvester was used to estimate the best number of populations using the ΔK method (Earl and vonHoldt 2012). The highest value for ΔK was obtained at K = 2 (ΔK = 137.01); indicating the entire germplasm could be divided into two groups (Fig. 5). The results indicated that STRUCTURE analysis, PCoA analysis and UPGMA analysis had similar results and the same germplasm was divided into two populations. In Group I, four germplasms belong to Zhejiang, two germplasms belong to Hunan, three germplasms belong to Hubei, and the rest belong to Guizhou, Shaanxi, Gansu. In Group II, Zhejiang contains four germplasms, Hunan, Yunnan, Guizhou, Sichuan contain two germplasms respectively and Hubei, Guangxi, Jiangxi, Fujian each have one.