Schematic representation of “the physiological constraint hypothesis” and “the life-history hypothesis” and their predictions. Physiological constraints in the environment refers to variation over time of the energetic balance (mainly for heterothermic 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 or thermal conditions or vice versa. It predicts that the sex difference in dormancy phenology is explained by sex differences in physiological constraints, and that 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 the 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 beneficial envioronmental conditions for 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. Black, grey and dark blue horizontal arrows represent respectively time over the year, reproductive investment in males and reproductive investment in females. Red and orange triangles represent respectively immergence and emergence timing and black cross represent mating timing.

World map showing the study sites of the hibernating species studied. Green dots indicate the location of study sites. A red gradient represents the variation in minimum ambient temperature, a parameter used in this study.

Effects of body mass variation during mating on protandry. The minimum temperatures of the study sites are indicated by a color gradient with the warmest temperatures in red. The regression line in black indicates the effect of body mass variation during mating on protandry. Body mass variation during mating was represented as a percentage of body mass before mating.

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.

Effects of body mass variation through the end of mating on sex differences in immergence. Annual precipitation at the study sites are indicated by a color gradient with the highest precipitation in dark blue. The regression line in black indicates the effect of body mass variation through the end of mating on sex difference in immergence. Body mass variation through the end of mating was represented as a percentage of body mass at emergence. Males immerged earlier in species living in an environment with high precipitation (dark blue dots).

Effects of female specific reproductive effort on sex difference in immergence. The annual precipitation of the study sites are indicated by a color gradient with the highest precipitation in dark blue. The regression line in black indicates the effect of female specific reproductive effort on sex difference in immergence. Males immerged earlier in species living in an environment with high precipitation (dark blue dots).

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 S2 for references).

Consensus phylogenetic trees for the species under study: (a) model 1 (b) model 2 and (c) model 3. 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. Protandry was used as the dependent factor and Body mass change during mating, late mating, body mass male, strategy fatstoring or foodstoring, minimum temperature, precipitation and dimorphism of body mass at emergence were considered as independent factors. Protandry was calculated as follows: Female Julian date – male Julian date. 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 2 and 3 (same data base without Tachyglossus aculeatus). Sex difference in immergence was used as the dependent factor whereas strategy fatstoring or foodstoring, body mass change through the end of mating, body mass immergence, dimorphism at immergence maternal effort, maternal effort duration, female specific reproductive effort, minimum temperature and precipitation were considered as independent factors. The exact hibernation phenology data for Cricetus cricetus have been confirmed by the authors. Weaning body mass for Zapus princeps and Urocitellus mollis was calculated from body mass data at 52-60 days and 29 days, respectively, and corrected at 48 and 34 days, corresponding to weaning dates. For maternal effort duration of Urocitellus brunneus and Ictidomys parvidens, we used respectively the averages for the clade Marmotini obtained from Hayssen (2008). The body mass at weaning for Urocitellus brunneus correpsonds to the body mass at weaning of Urocitellus townsendii (phylogenetically close species) corrected for the difference in body mass of the females from Hayssen (2008) as body mass accounts for 76–84% of the variation in weaning mass in the clade Marmotini (Hayssen, 2008) Sex difference in immergence was calculated as follows: female Julian date – male Julian date.