Schematic representation of “the phenology and energy limitation hypothesis” and “the phenology and life-history balance hypothesis” and their predictions. Energetic limitation in the environment refers to variation over time of the energetic balance (mainly for endotherms) and the thermal window favorable to activity (mainly for ectotherms). Hypothesis H1 assumes that dormancy phenology occurs at the time of transition between favorable and unfavorable energetic conditions or vice versa. It predicts that the sex difference in dormancy phenology is explained by differences in energy limitation, and reproductive investment should be independent of this sex difference in phenology. In contrast, hypothesis H2 predicts that a phenology that would occur before or after this energetic transition may be associated with benefits to survival or reproduction. It predicts that the sex difference in dormancy phenology is associated with a sex difference in reproductive investment. This pattern is expected for species without paternal effort. But the concept can be applied to other types of mating strategies. The hibernation phenology presented for prediction (H1) are those expected for hibernating mammals. Note that the magnitude and order of sex differences in phenology is not an expected general trend, because it is assumed to vary between species according to energy demand (prediction H1) and reproductive investment (prediction H2). Nevertheless, the sex difference is assumed to be smaller with the H1 prediction because there is less sex difference in energy demand than in reproductive investment. Black, grey and dark blue horizontal arrows represent respectively time over the year, reproductive investment in males and reproductive investment in females.

Summary of full models tested and sample size. Crosses indicate variables included in the models. Stars indicate factors for which interactions were tested with log transformed relative testes mass (model 5) or Δ body mass during mating (model 6). The abbreviation “diff” stands for “difference”.

Regression results for the best models explaining preliminary assumption. The Z standardized model estimates and the phylogenetic effect are reciprocally estimated by β and γM. The factor “body mass gain before mating” corresponds to all the species that have a positive value for the factor “Δ body mass before mating”. Relative testes mass, Δ body mass before mating, Δ body mass during mating, Δ body mass through the end of mating was represented respectively as a percentage of body mass, body mass at emergence, body mass before mating and body mass at emergence.

Effects of relative testes mass (log-transformed) on protandry. The minimum temperatures of the study sites are indicated by a color gradient with the warmest temperatures in red. The regression lines in red, light grey and blue indicate respectively the effect of log-transformed relative testes mass on protandry when the annual minimum temperature is equal to the max, mean and min value among study sites.

Species with dimorphisms biased in favor of males or females and their body mass gain during the year. Body size dimorphism is calculated as male body size divided by female body size. See section “Sex differences in reproductive investment” for the determination of the body mass gain before mating. Superscript numbers indicate bibliographic references

Effects of Δ body mass before mating on protandry. The delay between female emergence and the beginning of the mating period is represented by a color gradient with the greatest delay in light blue. The regression line in red indicates the effect of Δ body mass before mating for the mean mating delay. Δ body mass before mating was represented as a percentage of body mass at emergence.

Effects of active time spent by males after mating on the sex difference in immergence date. The regression line in red indicates the effect of post-mating activity time for the mean maternal effort. The duration of the maternal effort is represented by a color gradient with the longest effort in blue. A negative value on the y-axis indicates that males immerge after females and a positive value indicates that males immerge before females.

Species with dimorphisms biased in favor of males or females and their body mass gain during the year. Body size dimorphism is calculated as male body size divided by female body size. See section “Sex differences in reproductive investment” for the determination of the body mass gain before mating (See supplementary materials S4 for references).

Consensus phylogenetic trees for the species under study: (a) model 1 and 2 (b) model 3 (c) model 4 (d) model 5 (e) model 6 (f) model 7 (g) model 8. Each consensus tree was built from 100 trees obtained from http://vertlife.org/phylosubsets/. Branch lengths were calculated using Grafen’s computations with the ‘ape’ package in R (see Materials and Methods).

Data on dependent and independent factors used in models 1 and 2. The body mass gain before mating was used as the dependent factor in model 1 and 2. The body mass change during mating and the relative testes mass (log-transformed) were considered as independent factors in models 1 and 2 respectively. The body mass and testes mass data were used to calculate relative testes mass. Stars indicate body mass extracted from graphs with the software Plot Digitizer at the time of the seasonal maximum in testes mass. For Cricetus cricetus and Tachyglossus aculeatus, the relative testes mass are directly available in the cited references.

Data on dependent and independent factors used in models 3 and 4. Active time after mating was used as the dependent factor in model 3 and 4. The body mass change during mating and the body mass change before and during mating were considered as independent factors in models 3 and 4 respectively. Immergence date for males Cricetus cricetus, used to calculate active time after mating, have been confirmed by the authors.

Data on dependent and independent factors used in models 5. Protandry was used as the dependent factor and relative testes mass (log-transformed), late mating, strategy fatstoring or foodstoring and minimum temperature were considered as independent factors. Protandry was calculated as follows: male Julian date – female Julian date. The body mass and testes mass data were used to calculate relative testes mass. Stars indicate body mass extracted from graphs with the software Plot Digitizer at the time of the seasonal maximum in testes mass. For Cricetus cricetus and Tachyglossus aculeatus, the relative testes mass are directly available in the cited references. The exact hibernation phenology data for Cricetus cricetus have been confirmed by the authors. See materials and methods for the acquisition of minimum temperature data.

Data on dependent and independent factors used in models 6. Protandry was used as the dependent factor and body mass change before mating, late mating, strategy fatstoring or foodstoring and minimum temperature were considered as independent factors. The exact hibernation phenology data for Cricetus cricetus have been confirmed by the authors. Protandry was calculated as follows: male Julian date – female Julian date. See materials and methods for the acquisition of minimum temperature data.

Data on dependent and independent factors used in models 7. Sex difference in immergence was used as the dependent factor whereas active time after mating and maternal effort were considered as independent factors. The exact hibernation phenology data for Cricetus cricetus have been confirmed by the authors. Sex difference in immergence was calculated as follows: male Julian date – female Julian date.

Regression results for the best models explaining sex difference in immergence without outlier. The Z standardized model estimates and the phylogenetic effect are reciprocally estimated by β and γM. The abbreviations “Diff” stands for “Difference”.

Regression results for the best models explaining variation in protandry and sex difference in immergence. The Z standardized model estimates and the phylogenetic effect are reciprocally estimated by β and γM. The abbreviations “diff” and “rel” “Min” stand respectively for “difference”, “relative” and “Minimum”. A negative value for the sex difference in immergence indicates that males immerge before females and a positive value indicates that females immerge before males. Relative testes mass, Δ body mass before mating, Δ body mass during mating, Δ body mass through the end of mating was represented respectively as a percentage of body mass, body mass at emergence, body mass before mating and body mass at emergence.