Oly methylation analysis, Nov 20, 2019
Revisiting Oly methylation data. We now have two lists of loci:
- 1) DMLs between two Olympia oyster populations, Hood Canal and South Sound, which were identified using MethylKit.
- 2) Loci where methylation status is associated with oyster shell length, filtered by a) loci have 10x coverage in all samples, and b) loci have 10x coverage in any sample.
Today I re-plotted heatmaps using MACAU loci, based on feedback from Steven & Katherine:
- Only use loci with 10x coverage
- Add heatmaps where samples are NOT ordered by cluster analysis, but instead by 1) tree from MethylKit, and 2) shell length. See my notebook, 06-analyzing-MACAU-results-rev1.html, and here’s one of the new heatmaps, with samples (columns) ordered by shell length, and a barplot of shell length below (red = Hood Canal oysters, green = South Sound oysters).
Then I used bedtools to see where DMLs and MACAU loci are located, see my notebook here:
07-DML-MACAU-annotation.ipynb.
Finally, I began annotating loci locations for DMLs and MACAU loci; see my notebook: 08-Annotations.html. Here’s a barplot showing which features overlap with the DMLs:
Written on November 20, 2019