Statistical analysis
Effects of workload, temperature and interactions on gonadal weight, body mass and gene expression levels were determined using type I two-way ANOVAs. Tukey HSD post-hoc pairwise comparisons were used to compare groups. Statistical significance was determined at α =0.05. Statistical results can be found in the supplementary information (Table S4). Analyses were performed using RStudio (version 1.2.1335) (R Core Team, 2013), and all figures were generated using the ggplot2 package (Wickham 2016).