Analyzing Catlin’s data

Let’s get right to it.

Temperature data

I analyzed the temperature data in this R Markdown script. HOBO data can be found here. All temperature calculations and statistical output can be found in this folder. Some notes because I’m too tired to do this any other way:

• WA HOBO loggers were put in during the heat ramp so there was no acclimation information. However, the tanks were kept on the counter without any crabs for a few days after the experiment so I can use that to extrapolate acclimation temperatures.
• ME had acclimation and ramp information.
• THERE WAS NO NH HEAT SHOCK! When I looked at the temperature data, there was maybe a brief heat ramp but it came down pretty quickly. I will use the NH TTR data from timepoint 0 and use it just for parasite and genotype analysis.

Figure 1. Temperature data

TTR data

I re-analyzed Catlin’s TTR data in this R Markdown script using similar mixed effects model methods to every other experiment. A couple of notes up top, then I’ll just get into results.

• The parasite data Catlin sent for NH is actually ME data! I emailed him to ask for the NH data, so I’ll revise that when I get to it
• Carolyn suggested I calculate a different parasite intensity for each parasite species by dividing parasite count by carapace width. This will really only come into play for the NH data since ME crabs were all infected by metacircarial cysts
• For ME, I included metacircarial parasite intensity in the model and tested its significance. I will do a similar thing with NH crabs once I have the correct parasite counts.
• All statistical output can be found in this file
• Average TTR by population at each timepoint can be found in this file
• Average TTR by timepoint and genotype for each population can be found in this file
• All exploratory figures can be found here
• All three populations were in Hardy-Weinberg equilibrium?!
• Only CC and CT crabs were tested in ME and NH. I’ll need to delineate between pre-experiment mortality and during experiment mortality in the results.

Alright! As expected, time, but not genotype, presence of the C allele, nor presence of the T allele significantly impacted righting response in WA crabs. Crab ID accounted for < 1% of variance. Time and genotype impacted righting response in ME crabs, with crab ID explaining 27.82% of the variance in the final model. Based on these results, I made a multipanel figure that showed righting response before and after heat stress in WA and ME. I also created a supplementary figure for the NH data.

Figure 2. Righting response figure

Figure 3. NH supplementary figure

There was a marginal impact of population on average TTR when I analyzed all the data together.

…..I think I just have to right now?!

Going forward

1. Revise NH parasite data analysis
2. Add Catlin’s methods and results to manuscript
3. Outline introduction
4. Outline discussion
5. Write introduction
6. Write discussion
7. Troubleshoot lipid assay protocol
8. Conduct lipid assay for crabs of interest