3.1 Genetic Diversity and Epigenetic Genetic Diversity Analysis
The numbers of alleles (Na) observed by SSR ranged from 1.656 to 1.750, while the number of effective alleles (Ne) range from 1.319 to 1.530. The average values of Na and Ne were 1.703 and 1.406, respectively. In MSAP, Na ranges from 1.689 to 2.000, while the Ne ranges from 1.148 to 1.344, with the average values of 1.798 and 1.275, respectively. The average values of gene diversity (H) and the Shannon diversity index (I) of SSR and MSAP Nei were 0.236, 0.184 and 0.354, 0.303, respectively, with little variance (Table 3). In summary, the genetic diversity of population in HK was the highest compared with other populations, key fact to remember its epigenetic diversity was significantly lower than that of the other three populations. The genetic diversity of JH is significantly lower than that of the other three regions. In the genetic distances (Table 4), HK and FZ were the closest, while GZ and JH were the farthest. The epigenetics of JH and FZ were the most recent, while JH and HK were the furthest ones.
Table 3 Genetic and epigenetic diversity analysis on the four Wedelia trilobata populations from different invasive areas
Population
|
|
PIC
|
Na
|
Ne
|
H
|
I
|
FZ
|
genetic
|
75.00%
|
1.750
|
1.375
|
0.228
|
0.349
|
epigenetic
|
83.47%
|
1.835
|
1.344
|
0.224
|
0.357
|
HK
|
genetic
|
65.62%
|
1.656
|
1.319
|
0.193
|
0.297
|
epigenetic
|
60.33%
|
1.603
|
1.148
|
0.101
|
0.174
|
JH
|
genetic
|
68.75%
|
1.688
|
1.399
|
0.232
|
0.348
|
epigenetic
|
75.21%
|
1.752
|
1.327
|
0.207
|
0.326
|
GZ
|
genetic
|
71.88%
|
1.719
|
1.530
|
0.291
|
0.421
|
epigenetic
|
100.00%
|
2.000
|
1.279
|
0.206
|
0.353
|
SSR Mean
|
|
70.30%
|
1.703
|
1.406
|
0.236
|
0.354
|
SSR Species level
|
|
|
1.969
|
1.528
|
0.313
|
0.472
|
MSAP Mean
|
|
79.75%
|
1.798
|
1.275
|
0.184
|
0.303
|
MSAP Species level
|
|
|
2.000
|
1.457
|
0.281
|
0.440
|
Table 4 The distance of genetic and the distance of epigenetic genetic within populations
|
genetic distance
|
|
epigenetic distance
|
Population
|
FZ
|
HK
|
JH
|
GZ
|
Population
|
FZ
|
HK
|
JH
|
GZ
|
FZ
|
***
|
0.875
|
0.870
|
0.134
|
FZ
|
***
|
0.906
|
0.962
|
0.837
|
HK
|
0.134
|
***
|
0.788
|
0.142
|
HK
|
0.099
|
***
|
0.018
|
0.647
|
JH
|
0.139
|
0.238
|
***
|
0.078
|
JH
|
0.039
|
0.983
|
***
|
0.722
|
GZ
|
0.874
|
0.868
|
0.925
|
***
|
GZ
|
0.178
|
0.436
|
0.325
|
***
|
Note: Nei’s genetic/epigenetic distance (left bottom) and genetic/epigenetic identity (top right).
3.2 Phylogenetic Tree
It was depicted in the Fig.1 that most samples from HK and GZ regions were obviously clustered into independent branches, and samples from JH and FZ regions were partially intertwined. The cluster of MSAP molecular markers was more independent compared with SSR in each region. Presumably, this may be due to the fact that MSAP amplifies more specific sites than SSR.
3.3 PCoA and AMOVA Analysis
The main coordinate (PCoA) results showed that the JH and GZ populations were close based on SSR , but they were separated from the FZ population (Fig.2a). The PCoA were consistent with the results from UPGMAclustering. Based on MSAP, FZ, HK and JH populations were close, the second axis separated the three populations from GZ (Fig.2b).
3.4 Mantel Test
Mantel analysis of relationship between geographical distance and genetic or epigenetic distance showed that the R values were 0.0093 and 0.475, respectively, both with the 0.01 P values (Fig. 3a, 3b). The Mantel analysis of relationship between genetic and epigenetic distances showed that the R value were 0.2 (Fig. 3c), suggesting that a stronger correlation as a contrast to genetic distance.
3.5 AMOVA analysis
AMOVA analysis (Table 5) revealed that the epigenetic fixation index (Fst) of W. trilobata in different regions was 0.475, which was higher than genetic fixation index (0.277). It was suggested that the differentiation among populations was mainly attributed to epigenetic. Genetic and epigenetic gene flow (Nm) were 0.652 and 0.276, respectively, which were less than 1, indicating that gene flow among W. trilobata populations in different regions was limited by geographical distance.
Table 5 AMOVA analysis of Carex accessions
Source
|
|
df
|
SS
|
MS
|
Est. Var.
|
%
|
Fst
|
Nm
|
P value
|
SSR
|
Among Pops
|
3
|
141.233
|
47.078
|
1.529
|
18%
|
0.277
|
0.652
|
0.001
|
|
Among Indiv
|
56
|
67.800
|
1.211
|
0.000
|
0%
|
|
|
|
|
Within Indiv
|
60
|
406.000
|
6.767
|
6.767
|
82%
|
|
|
|
|
Total
|
119
|
615.033
|
|
8.296
|
100%
|
|
|
|
MSAP
|
Among Pops
|
3
|
3806.350
|
1268.783
|
40.222
|
47%
|
0.475
|
0.276
|
0.001
|
|
Among Indiv
|
56
|
3478.133
|
62.110
|
17.646
|
21%
|
|
|
|
|
Within Indiv
|
60
|
1609.000
|
26.817
|
26.817
|
32%
|
|
|
|
|
Total
|
119
|
8893.483
|
|
84.686
|
100%
|
|
|
|
3.6 Mode Structure of Epigenetic Regions
The total methylated rates in the four populations (Fig. 4) were 23.82 %, 18.25 %, 32.92 % and 27.55 %, respectively; the fully methylated rates were 16.25 %, 11.51 %, 25.75 % and 24.43 %, respectively; the hemi-methylation rates were 7.57 %, 6.75 %, 7.16 % and 3.12 %, respectively; moreover, the non-methylated rates are 76.18 %, 81.75 %, 67.08 % and 72.45 %, respectively.