Comparison of results for filtered epiRAD data without outliers identified by DAPC
Comparison of results for the following data after adding an additional filter (Filter2b) for removing outlier samples identified by DAPC and rerunning the full epiRAD pipeline.
Number of individuals left after this additional filter is 33.
- Filter4: Data obtained post BLAST filter, low read filter and low coverage filter
- Filter 5a: Data obtained post filter4 AND removing problematic loci identified by rad_haplotyper
- Filter 5b: Data obtained post filter4 AND keeping only those loci that remain after full SNP filtering pipeline
| Filter# | Number of loci | Methylation cutoff |
|---|---|---|
| 4 | 7795 | -2.848 |
| 5a | 6809 | -2.883959 |
| 5b | 5912 | -2.830382 |
Result1: MDS without dapc outliers
| Filter4 | Filter5a | Filter5b |
|---|---|---|
![]() |
![]() |
![]() |
Result2: Methylation heatmap clustered by rows and columns without dapc outliers
| Filter4 | Filter5a | Filter5b |
|---|---|---|
![]() |
![]() |
![]() |
Result3: %Methylation without dapc outliers
ANOVA for % methylation per site is not significant.
| Filter4 | Filter5a | Filter5b |
|---|---|---|
![]() |
![]() |
![]() |
Result5a: Pst
| Filter4 | Filter5a | Filter5b | |
|---|---|---|---|
| Mean Pst | 0.127 | 0.132 | 0.132 |
| Filter4 | Filter5a | Filter5b |
|---|---|---|
![]() |
![]() |
![]() |
Result5b: Pairwise Pst
| Table 1 Filter4 | Table 2 Filter5a | Table3 Filter5b |
|---|---|---|
| | | EOB | PBF | WOB | WOF | |-----|-----|-------|-------|-------| | EOB | x | 0.217 | 0.053 | 0.02 | | PBF | | x | 0.2 | 0.218 | | WOB | | | x | 0.072 | | WOF | | | | x | | | | EOB | PBF | WOB | WOF | |-----|-----|-------|-------|-------| | EOB | x | 0.22 | 0.055 | 0.018 | | PBF | | x | 0.199 | 0.219 | | WOB | | | x | 0.075 | | WOF | | | | x | | | | EOB | PBF | WOB | WOF | |-----|-----|-------|-------|-------| | EOB | x | 0.22 | 0.055 | 0.018 | | PBF | | x | 0.199 | 0.219 | | WOB | | | x | 0.075 | | WOF | | | | x | |
Result5c: Pooled Back Reef vs Pooled Fringe Reef Pst
| Filter4 | Filter5a | Filter5b | |
|---|---|---|---|
| Mean Pst | 0.053 | 0.053 | 0.053 |
find.clusters
find.clusters identifies one cluster for data post Filter4,5a,5b.
PERMANOVA
by site
| Filter4 | Filter5a | Filter5b | |
|---|---|---|---|
| R2 | 0.12809 | 0.12993 | 0.12993 |
| Pvalue | 0.011 | 0.011 | 0.007 |
by loc
| Filter4 | Filter5a | Filter5b | |
|---|---|---|---|
| R2 | 0.06397 | 0.06215 | 0.6215 |
| Pvalue | 0.001 | 0.002 | 0.003 |
Written on July 16, 2020











