Cold Acclimation Green Crab Experiment Part 35
Revising WA TTR data
Now that I have all my WA genotypes, I can revise my TTR analysis! All analysis was done in this R Markdown script.
Mixed effects models
After loading in the revised genotype data I noticed immediately that I FORGOT TO EXTRACT A CRAB. The row for crab 152 was hidden on my Google Sheet so I missed it. RIP.
I still looked at the new distribution of genotypes in the population. The population is not in Hardy-Weinberg equilibrium (p = 0.000036240950966106)
Figure 1. Genotype distribution of WA population
In any case, I used the genotype data for all but one crab to revise the stats. As expected, there was no impact of genotype on average TTR or failure to right.
Figures
My next step was to revise all my figures!
Figures 2-7. Various visualizations for crabs that failed to right
Figures 8-11. Various visualizations for Q10 differences for 28 crabs followed throughout the experiment.
Figures 12-13. Various visualizations for individual-level differences for crabs followed throughout the experiment and some extras as well.
Discrepancies between genotype sets
In making my figures, I noticed that previous versions had TT crabs failing to right, but in this version none of my TT crabs failed to right. This means that there are crabs I initially called TTs that were actually CTs. To identify these crabs, I went to a previous genotype commit. Crabs 154, 156, and 215 were called TTs even though the 12/10 gel had faint bands. When they were redone on 2/26 and 2/28, the CT band patterns were much clearer! Looking at the 2/28 gel, I also noticed that 257 was called a CT even though the gel looks like a TT with a false bottom band based on the spacing between the two bands. I will redo this PCR when I extract WA 152.
Going forward
- Write methods section in manuscript
- Add MA results to manuscript
- Add WA results to manuscript
- Add Catlin’s methods and results to manuscript
- Genotype Catlin’s samples
- Determine methods for comparing population responses
- Troubleshoot lipid assay protocol
- Conduct lipid assay for crabs of interest