Habitat characteristics and life history explain reproductive seasonality in lagomorphs

Lagomorphs show extensive seasonal variation in their reproduction. However, the factors causing this large variation have so far mostly been investigated intraspecifically and therefore provide only some exemplary comparisons of lagomorph reproductive seasonality. The present study applies both a categorical description (birth season categories 1–5) and a quantitative measure (birth season length in months) to summarize the degree of birth seasonality in the wild of 69 lagomorph species. Using a comparative approach, I tested the influence of 13 factors, comprising six habitat, five life history and two allometric variables on birth season length in lagomorphs. Leporids mainly show non-seasonal birthing patterns with high intraspecific variation. Their opportunistic breeding strategy with high reproductive output and their large distribution areas across wide latitude and elevation ranges might be the reasons for this finding. Ochotonids reproduce strictly seasonally, likely because they live at northern latitudes, are high-altitude specialists, and occur in limited distribution areas. The most important factors associated with variation in lagomorph birth seasonality are mid-latitude, mean annual temperature and precipitation of a species’ geographical range and life history adaptations including fewer but larger litters in seasonal habitats. Birth seasons become shorter with increasing latitude, colder temperatures, and less precipitation, corresponding to the decreasing length of optimal environmental conditions. Leporid species with shorter breeding seasons force maternal resources into few large litters to maximise reproductive output while circumstances are favourable. Since allometric variables were only weakly associated with reproductive seasonality, life history adaptations and habitat characteristics determine birth seasonality in Lagomorpha.


Introduction
Reproduction is energetically very expensive for most organisms. This is particularly true for placental mammals, where females bear the metabolic costs of pregnancy and lactation (McNab 2006;Speakman 2008;Heldstab et al. 2017). For mammals living in seasonal habitats with frequent periods of food scarcity, it is therefore essential that birth and lactation occur during the most favourable time of the year to increase their reproductive success (e.g., Lancaster and Lee 1965;Soto et al. 2004).
Lagomorphs show extensive seasonal variation in their reproduction. In some species, such as many pikas (Ochotonidae), all breeding activity takes place during a very restricted part of the year of only around one to two months (e.g., Chapman and Flux 1990;Smith et al. 2018). Other lagomorph species are capable of reproducing during a broader time window (e.g., Barkalow 1962;Newson 1964;Franken and Hik 2004;Rioja et al. 2011) or throughout the year but show seasonal birth peaks (e.g., Prakash and Taneja 1969;Stubbe et al. 1991;Vázquez et al. 2007). Lastly, in a fourth group of species, reproduction is entirely non-seasonal (Hayssen et al. 1993;Wilson et al. 2016).
Several ecological, physiological, and social factors have been proposed to explain this broad spectrum of reproductive seasonality in Lagomorpha (e.g., Barkalow 1962;Conaway et al. 1974;Swihart 1984;Boyd 1985;Ben Saad and Maurel 2002;Alves and Rocha 2003;Franken and Hik 2004;Portales et al. 2004;Rödel et al. 2005;Virgos et al. 2006;Tablado et al. 2009;Rioja et al. 2011;Schai-Braun et al. 2017). However, this knowledge derives almost exclusively from single-species studies and therefore provides only some exemplary comparisons of seasonality in reproduction between different lagomorph species. Furthermore, future studies enhancing our understanding of reproductive behaviour in the order Lagomorpha, with approximately one-fifth of species being threatened with extinction (Smith 2018), are also fundamental to captive-breeding programs and conservation management (Côté 2003). The importance of lagomorphs both as a pest of agricultural significance and as a key prey species and major human food resource additionally emphasizes the need to understand what controls their reproductive activity (Cowan and Bell 1986;Tablado et al. 2009;Delibes-Mateos et al. 2011;Lumpkin and Seidensticker 2011;Smith et al. 2019).
The first aim of this paper is therefore to collate and summarize data on birth seasonality from the wild of the majority of lagomorph species and compare them interspecifically. The second objective is to assess the influence of 13 factors on reproductive seasonality in lagomorphs that have been suggested in previous studies to play a role in birth seasonality: namely, six habitat characteristics (mid-latitude and latitude range of a species' geographical distribution (Sadleir 1969;Conaway et al. 1974;Bronson 1988), altitude (Conaway et al. 1974;Schai-Braun et al. 2017), temperature (Portales et al. 2004;Rödel et al. 2005;Vázquez et al. 2007;Tablado et al. 2009;Portales-Betancourt et al. 2012;Wells et al. 2016), precipitation (Wood 1980Franken and Hik 2004;Portales et al. 2004;Vázquez et al. 2007;Rioja et al. 2011;Portales-Betancourt et al. 2012), and home range size (Innes and Millar 1995;Spady et al. 2007)), two allometric variables (sexual size dimorphism (Valdespino 2007;Zerbe et al. 2012;Heldstab 2021) and body mass (Janson and Verdolin 2005;Rödel et al. 2005;Tökölyi et al. 2014)), and five life history traits (litter size, litters per year, gestation length, weaning age, and age at sexual maturity (Lord 1960;Conaway et al. 1974;Chapman 1984;Swihart 1984;Bronson 1989;Tablado et al. 2009;Tökölyi et al. 2014)). The influence of these 13 factors on birth seasonality will be tested across the whole order Lagomorpha and for each of the distinct four ecotypes (hares and jackrabbits, rabbits, rock-and talus-dwelling pikas, and burrowing pikas) separately. For a more detailed description of how these factors may influence birth seasonality, see Heldstab (2021) and supplementary material of Heldstab et al. (2018).
Lastly, the third aim is to compare reproductive seasonality of Lagomorpha with other mammalian orders. Relative to most other mammals, lagomorphs are characterised by extremely rapid foetal and neonatal growth rates (Case 1978;Eisenberg 1981;Swihart 1984) including post-partum oestrus (females can simultaneously suckle one litter while carrying another in the uterus) and superfetation in the genus Lepus (female hares are able to carry a second fertilized ovum in their uterus already containing a pregnancy) (Flux 1967). These characteristics allow lagomorphs to accomplish a reproductive cycle and maximize litter production within a very short time period but impose high energetic demands on females (Flux 1981;Keith 1981;Swihart 1984). Lagomorph species are therefore expected to breed only during optimal environmental conditions in seasonal habitats resulting in a pronounced birth season for species living at high latitudes and/or elevations. On the other hand, the fast reproductive cycle allows lagomorph species to breed opportunistically, e.g., responding to a sudden increase in food availability after rainfall, which may result in a non-seasonal birthing pattern with high intraspecific variation for those species living in regions with large variations in temperature and/or precipitation.

Materials and methods
To investigate which factors influence reproductive seasonality in the natural habitat of different lagomorphs, a broad dataset on birth seasonality covering 69 species was compiled from the literature. Birth seasonality was described both in a quantitative and a categorical way (Table 1). As a continuous measure, birth season length in months was used as in Heldstab (2021). For the categorical measure, each species was assigned to one of five birth seasonality categories based on its birthing pattern following (Zerbe et al. 2012;Heldstab et al. 2018Heldstab et al. , 2021Heldstab 2021): Category 1: narrow peak (less than approximately 60 days), no births for the rest of the year; category 2: expanded peak (> 60 days) or double peak, no births for the rest of the year; category 3: peak(s) accompanied by a small number of births throughout the year; category 4: births throughout the year, undulating pattern (preferred seasons possible). category 5: constant births throughout the year without a clear preferred season. If several sources of data were available for one species, preference was given to the study with the highest number of births, the longest study period or the birth season most often mentioned in the literature.
Data on all explanatory variables, namely, the habitat characteristics (mid-altitude, latitude range, mean altitude, mean annual temperature and precipitation of the species' geographical distribution, and home range size), allometric variables (adult body mass and sexual size dimorphism expressed as the ratio of male vs. female body mass), and life-history traits (gestation length, litters per year, litter size, weaning age, and age at sexual maturity) were compiled from the published literature. Although litter size is a
In the present study, preference was given to litter size values deriving from the number of young at parturition and litter size and litters per year values most commonly mentioned in the literature. Detailed sources of data for the whole dataset are given in Table 1. Statistical analyses were performed in JMP™ 13.0 (SAS Institute Inc 1989-2016) and in R 4.0.4 (R Core Team 2021). All analyses were controlled for phylogenetic non-independence using phylogenetic generalized least-squares regressions (PGLS: Martins and Hansen 1997) using the "caper" package (Orme et al. 2013) in R. PGLS is a modification of generalised least squares using the knowledge of phylogenetic relationships to produce an estimate of expected covariance in cross-species data (Symonds and Blomberg 2014). Closely related species are assumed to have more similar traits because of their shared ancestry and hence produce more similar residuals from the least-squares regression line. By taking into account the expected covariance structure of these residuals, modified slope and intercept estimates are generated that can account for interspecific autocorrelation due to phylogeny. Caper estimates PGLS model parameters in maximum likelihood (Orme et al. 2013) and the parameter lambda (λ), which quantifies the magnitude of the phylogenetic signal in the model residuals (Freckleton et al. 2002). The value of λ can vary between 0, indicating no phylogenetic signal, and 1, indicating that the observed pattern fits a Brownian motion model of trait evolution along the branches of the phylogeny such that similarity between species is directly proportional to relatedness.
The lagomorph phylogeny was extracted from the composite mammalian supertree including branch length estimations provided by Fritz et al. (2009) and is given in Fig. 1 To determine the most important explanatory variables of birth seasonality either measured in categories or number of Fig. 1 The phylogenetic tree used in the present study from Fritz et al. (2009), with an indication how a the birth season categories and b the birth season length in months are distributed across the tree, ranging from very short (a highly seasonal reproductive pattern, red) to very long (a non-seasonal reproductive pattern, dark blue) months, all thirteen explanatory variables mentioned above were tested in univariate PGLS analyses. Unfortunately, performing multivariate regressions or multi-model averaging was impossible because each explanatory variable had a different sample size (due to missing values for several understudied species) and there was problematic high multicollinearity between predictors (Rogerson 2001) (Online Resource 1). The critical level of significance (P = 0.050) was adjusted by Bonferroni correction for multiple testing to (P = 0.004). Model assumptions of normality and homogeneity of residuals were checked by visually inspecting histograms and qq-plots of the residuals as well as by plotting residuals against fitted values. The values of home range size, adult body mass, gestation length, weaning age, and age at sexual maturity were log-transformed before analysis to reduce the skew of their distribution. Additionally, to investigate the effect of habitat characteristics, allometric variables and life-history traits on birth seasonality in more detail, the four distinct ecotypes of lagomorphs (hares and jackrabbits, rabbits, rock-and talus-dwelling pikas, and burrowing pikas) were each analysed separately.

Results
A significant positive correlation between the five ordinal birth season categories and the continuous variable, the birth season length in months, was found (PGLS: p < 0.001, n = 69, λ < 0.001, β = 1.055, SE = 0.019, t = 14.000), indicating that the quantitative and the categorical approach yielded a comparable ranking of species. The ordinal birth seasonality categories showed a right-skewed distribution with the highest number of species in category 2, indicating that the majority of lagomorphs reproduce during an expanded but determined time window of more than 2 months with no births for the rest of the year (Fig. 2a). The birth season length in months ranged from 1 to 12 months and its distribution was also right skewed but showed a second peak at birth season length of 12 months with 14.5% of species reproducing year-round (Fig. 2b). Lagomorphs show a high degree of intraspecific variation in birth season length. As evident in Table 1, for many species, various birth season lengths have been described in the literature ranging in the most extreme cases of volcano rabbits (Romerolagus diazi) and cape hares (Lepus capensis) from 3 to 12 months and of mountain hares (Lepus timidus) from 1 to 9 months.
The phylogenetic patterns of reproductive seasonality are given in Fig. 1a and b. The most strictly seasonally reproducing group were the pikas (Ochotonidae) with an average birth season length of 3.23 months and the most seasonal species, the collared pika (Ochotona collaris) and the Ladak pika (Ochotona ladacensis) followed by hares with an average birth season length of 6.60 months, and the rabbits with an average birth season length of 7.73 months.
The most important factors influencing the degree of birth seasonality in lagomorphs were mid-latitude, mean annual temperature and precipitation of a species' geographical distribution, litter size, and litters per year (Tables 2 and 3). Birth seasonality either measured in months or categories was negatively correlated to the mid-latitude of a species' geographical distribution ( Fig. 3a and b) showing that reproductive periods became shorter with increasing latitude. Mean annual temperature ( Fig. 3c and d) and precipitation ( Fig. 3e and f) were positively correlated with reproductive seasonality, indicating that lagomorph species living in warmer habitats with more precipitation exhibit an expanded birth season. Litters per year was positively ( Fig. 3g and h) and litter size negatively ( Fig. 3i and j) correlated with birth season length, suggesting that species with a more pronounced birth seasonality have larger but fewer litters per year than species which have a more expanded birthing pattern. Adult body mass was positively linked to birth seasonality but only if measured in months ( Fig. 3k and l), where heavier species are less likely to be seasonal breeders. Latitude range, mean altitude of a species' geographical distribution, home range size, gestation length, weaning age, age at sexual maturity, and sexual size dimorphism had no influence on reproductive seasonality in lagomorphs (Tables 2 and 3).
Investigating each of the four distinct ecotypes separately, I found that for hares and jackrabbits, the most important factors influencing birth seasonality were litter size, litters per year and mean annual temperature (Online Resources 2 and 3). Mean annual temperature and litters per year were positively and litter size negatively correlated with birth seasonality either measured in the number of months or categories. Additionally, mean annual precipitation showed a positive and mid-latitude negative relationship with reproductive seasonality in hares and jackrabbits but only if birth seasonality was measured in categories (Online Resources 2 and 3). For the other three ecotypes (rabbits, rock-and talus-dwelling pikas, and burrowing pikas), none of the 13 tested factors had a significant influence on birth seasonality after Bonferroni correction for multiple testing (Online Resources 4-9). The only factor close to significance was litters per year in rabbits showing a positive trend on birth season length measured in months (Online Resources 4).

Discussion
Lagomorphs occur worldwide (except in Antarctica) across a wide range of habitats spanning altitudes from sea level to more than 5000 m and latitudes from the equator up to 80°N (Wilson et al. 2016;Smith et al. 2018). Not surprisingly, lagomorphs are important for maintaining the equilibrium  Ln (body mass in grams) of many different ecosystems given that they are primary consumers and key prey species for a large array of predators (e.g., Tablado et al. 2009;Vaughan et al. 2013;Badingqiuying et al. 2016). This emphasizes the need to understand what controls population dynamics in Lagomorpha including their reproductive behaviour. By applying both a categorical description (birth season categories 1-5) and a quantitative measure (birth season length in months), this study characterizes the degree of reproductive seasonality in the wild of 69 lagomorph species and tests the influence of 13 variables that were previously suggested to play a role in birth seasonality. I found that many lagomorph species exhibit a high degree of intraspecific variation in birth season length with the majority of species reproducing for an expanded but determined period of more than two months. I further show that the most important factors influencing the degree of birth seasonality were mid-latitude, mean annual temperature, and precipitation in a species' geographical range, litter size and litters per year. Although less strongly, body mass was also correlated with birth season length. These results signify that environmental conditions of the natural habitats, as well as life history adaptations, determine birth seasonality in the order Lagomorpha.

Characterization of birth seasonality in Lagomorpha
Consistent with previous studies (e.g., Conaway et al. 1974;Smith 1978;Wood 1980;Gibb et al. 1985;Franken and Hik 2004;Tablado et al. 2009), lagomorphs show a large intraspecific variation in birth seasonality. This is especially true for Leporidae (rabbits, jackrabbits and hares) with wide geographical ranges (Schai-Braun and Hackländer 2016) and some very successful invasive species that have been introduced in numerous countries (Flux and Fullagar 1992;Montuire 2001;Chapman and Flux 2008). Since populations of several leporid species are known to adjust the length of their birthing season to their geographical location (Conaway et al. 1974;Bronson 1989;Tablado et al. 2009), large within-species variation might be a consequence of their extensive distribution. In species with a more limited geographic range, intraspecific variation in birth seasonality might occur because lagomorphs are known to reproduce very opportunistically with high year-to-year variability, depending on various factors such as weather and snow cover (Flux 1970;Franken and Hik 2004;Portales et al. 2004;Rödel et al. 2005;Morrison and Hik 2007), food availability (Flux 1970;Millar 1972;Smith 1978) and predation risk influencing population density (Wood 1980;Boonstra et al. 1998).
In contrast to leporids, most pikas (Ochotonidae) reproduce strictly seasonally. Many species of this family live at northern latitudes Thapa et al. 2018), are high-altitude specialists (Smith 1974) and occur in limited distribution areas (Leach et al. 2015). These characteristics might explain their seasonal birthing pattern because the period of optimal environmental conditions for parturition and subsequent conception is short in northern alpine environments (Sparks and Menzel 2002), particularly for mammalian herbivores. The snow-free season is often shorter than two months and young must therefore be born late enough to avoid snow or low-quality food resources and early enough to allow growth and accumulation of body fat stores to survive the subsequent winter (O'Donoghue and Boutin 1995;Franken and Hik 2004). Thus, there is great selective pressure for early births, and short birthing periods for northern mountain-dwelling mammals such as many pika species (Smith 1978) (see below for the influence of latitude and altitude on birth seasonality).

Factors associated with seasonal reproduction in Lagomorpha
Habitat characteristics: latitude, latitude range, altitude, temperature, precipitation and home range size Overall, this study shows that mid-latitude, mean annual temperature and precipitation of a species' habitat have a significant influence on reproductive seasonality in lagomorphs. Previous intraspecific studies in lagomorphs similarly confirmed a strong relationship between latitude and birth seasonality (Barkalow 1962;Conaway et al. 1974;Flux 1981;Gibb et al. 1985;Bronson 1989;Chapman and Flux 2008;Tablado et al. 2009). Birth windows became narrower with increasing latitude of the natural habitat, because spring starts later and winter starts earlier in the year, resulting in a decline of the length of the favourable season at higher latitudes (Sparks and Menzel 2002). Since reproduction requires significant energy investments, the optimal timing of parturition to avoid hostile conditions when offspring are unlikely to survive exerts a strong selective pressure for short birthing periods at high latitudes (Prendergast et al. 2002). Latitudinal range, however, had no influence on reproductive seasonality in lagomorphs. In comparable studies across rodent, primate and Carnivora species, I found that species with a wider latitudinal range had longer birth seasons (Heldstab et al. 2018Heldstab 2021). However, the effect of latitudinal range was always much smaller than the effect of mid-latitude.
Another factor that is supposed to have similar effects as latitude on the reproductive strategy of lagomorphs is the mean altitude of species distribution. As for increasing latitude, the favourable season becomes shorter with increasing altitude (Mooney and Billings 1961;Sadleir 1969; Körner 1 3 2007). I therefore expected that the birth season would be more pronounced for species living in high-altitudinal habitats. However, unlike previous studies in the Eastern cottontail (Sylvilagus floridanus) and Alpine mountain hares (L. timidus varronis), which found an effect of altitude on the onset of reproduction (Conaway et al. 1974) or on litter size (Schai-Braun et al. 2017), I found no influence of altitude on reproductive seasonality across Lagomorpha.
As shown here and in previous single-species studies, habitats or years that were warmer (Hewson and Taylor 1975;Frylestam 1980;Myrberget 1983;Angerbjörn 1986;Tablado et al. 2009) and had a higher amount of precipitation (Vorhies and Taylor 1933;Bothma and Teer 1977;McKay and Verts 1978;Wilde 1979;Wood 1980) were associated with less seasonal reproduction. Low temperatures and low amounts of precipitation are associated with reduced food availability (e.g., reduced vegetation growth) (Boutin and Larsen 1993;Tablado et al. 2009) and with increased thermal stress (Gilbert et al. 1987;Rödel et al. 2004), resulting in poor body condition and a temporary cessation of breeding (Flux 1970;Andersson et al. 1981;Bell and Webb 1991;Rödel et al. 2005) consequently leading to a more seasonal reproduction.

Life history traits: gestation length, litter size, litters per year, weaning age and age at sexual maturity
Highly seasonal habitats provide shorter time windows of optimal conditions for breeding. Female mammals living in such environments are therefore expected to force maternal resources into a few large litters to maximise reproductive output while circumstances are favourable. Consistent with this prediction, I found that the majority of lagomorph species (with the possible exception of Ochotonidae, see below for the results of the four distinct ecotypes) with shorter birth seasons do have larger but fewer litters per year than species which a more expanded birthing pattern. A smaller comparative study across 22 lagomorph species also found that the number of litters decreases and the litter size augments with shorter breeding seasons (Swihart 1984). Studies using proxies for birth season length, such as latitude or temperature, similarly confirm the negative relationship between breeding season length and litter size intraspecifically in European hares (Lepus europaeus) (Hewson and Taylor 1975), Eastern cottontails (S. floridanus) (Conaway et al. 1974) and snowshoe hares (Lepus americanus) (Keith et al. 1966) and interspecifically in the genus Lepus (Flux 1981) and Sylvilagus (Lord 1960) and the family Leporidae (Virgos et al. 2006). Another study in European rabbits (Oryctolagus cuniculus) supports the present finding that the number of litters per year increased with the length of the breeding season (Tablado et al. 2009). In sum, there is ample evidence for Leporidae that litter size and the number of litters per year are adjusted according to their birth season length (see discussion of results for Ochotonidae below).
Other ways in which mammals gain selective advantages in seasonal habitats is by shortening gestation length and/ or accelerating foetal and neonatal growth to ensure that mating, birthing and raising offspring occur at a beneficial time of the year (Kiltie 1988;Zerbe et al. 2012;Heldstab et al. 2018Heldstab et al. , 2021Heldstab 2021). Across eight New World rabbits, Chapman (1984) showed that species living at northern latitudes with shorter breeding seasons had a shorter gestation length. However, a higher degree of reproductive seasonality was not linked to a shorter gestation length in the present study, correcting for phylogenetic non-independence of taxa and covering a much higher number of species. In support, a study across 17 leporids was also not able to find a relationship between the variation of monthly average temperatures and gestation times (Virgos et al. 2006). Relative to other eutherian mammals, lagomorphs exhibit very short gestation periods (Swihart 1984;Wilson et al. 2016). Furthermore, females of many lagomorph species experience postpartum oestrus so that they can simultaneously suckle one litter while carrying another in the uterus with only a small increase in the duration of pregnancy (Millar 1973(Millar , 1974Flux 1981;Keith 1981). In the genus Lepus, females are also known to superfetate, carrying a second fertilized ovum in their uterus during an already ongoing pregnancy (Flux 1967). These characteristics might be the reason for the absence of a correlation between gestation and birth season length in lagomorphs.
Concerning offspring growth in lagomorphs, Swihart (1984) found a correlation between the length of the breeding season and the age of first reproduction, but weaning age and age at sexual maturity had no influence on reproductive seasonality in the current study. These contrasting results suggest that the link between birth seasonality and ontogenetic development deserves future study in lagomorphs.

Allometric variables: adult body mass and sexual size dimorphism
Sexual size dimorphism had no influence on reproductive seasonality in lagomorphs. Another allometric variable, namely adult body mass, had no effect on birth seasonality measured in categories but correlated positively with birth season length measured in months. To my knowledge, there are no published studies investigating the relationship between breeding season length and either body mass or sexual size dimorphism in the order Lagomorpha. Analysing the breeding season length data of Swihart (1984) with a smaller sample size but a more precise measure of birth seasonality in days, the relationship between birth season length and both body mass (PGLS: p = 0.326,n = 20,λ < 0.001,β = 24.195,SE = 23.953, t = 1.010) as well as between birth season length and sexual size dimorphism (PGLS: p = 0.387,n = 19,λ < 0.001,SE = 369.580, were not significant either. In Swihart's study, environmental conditions also had a stronger influence than body mass on several other reproductive traits. Davis and Roth (2008) found little evidence for a correlation between sexual size dimorphism and various fecundity-related traits in comparisons among cottontail rabbits.

Factors associated with seasonal reproduction within the four distinct ecotypes
As shown for the whole order of Lagomorpha, habitats that were warmer, had a higher amount of precipitation and were located at lower latitudes were also associated with less seasonal reproduction in hares and jackrabbits. These findings are consistent with single-species studies in this ecotype, e.g., in European hares (L. europaeus) (Hewson and Taylor 1975;Frylestam 1980), mountain hares (L. timidus) (Myrberget 1983;Angerbjörn 1986), antelope jackrabbits (Lepus alleni) (Vorhies and Taylor 1933), and black-tailed jackrabbits (Lepus californicus) (Vorhies and Taylor 1933).
Again consistent with the results across the whole Lagomorpha dataset, I also found that hares and jackrabbits with shorter birth seasons do have larger but fewer litters per year than species which a more expanded birthing pattern. Similarly, in rabbits, there was a tendency for species with a more seasonal reproduction to have fewer litters per year. Previous studies in hares and jackrabbits likewise confirm that litter size is inversely related to the duration of the reproductive season (Keith et al. 1966;Hewson and Taylor 1975;Flux 1981). Earlier published work, for instance in mountain hares (L. timidus) (Flux 1970;Kauhala et al. 2005), Eastern cottontails (S. floridanus) (Conaway et al. 1974) and European rabbits (O. cuniculus) (Tablado et al. 2009), supports the positive relationship between the number of litters per year and the length of the breeding season found in the present study within rabbits and within hares and jackrabbits.
Within the two ecotypes burrowing pikas and rock-and talus-dwelling pikas, none of the 13 factors had a significant influence on birth seasonality. However, the datasets comprising in most of the analyses less than 10 species are very small and problematic for phylogenetic comparative analyses (Münkemüller et al. 2012). It is therefore difficult to draw any firm conclusion from these analyses of the two ecotypes of pikas. As for many lagomorph species too, additional data for pikas are needed to solve the lack of knowledge concerning their reproductive seasonality and life history to make phylogenetic comparisons such as the present one possible.

Comparison of reproductive seasonality in lagomorphs with other mammalian groups
Similarly to rodents (Heldstab 2021) but unlike Carnivora (Heldstab et al. 2018), ruminants (Zerbe et al. 2012) and primates , many lagomorph species also show a large degree of intraspecific variation in birth season length. A possible explanation for this high withinspecies variation is that rodents and lagomorphs have a faster reproductive cycle than the other three orders (Case 1978;Eisenberg 1981;Swihart 1984) including a post-partum oestrus (Dewsbury et al. 1979) and superfetation in the genus Lepus (Flux 1967). This allows them to breed much more opportunistically and adjust birth season length according to predation risk (Wood 1980;Boonstra et al. 1998), weather (Bronson 1988;Franken and Hik 2004) or food abundance (Leslie et al. 1952;Millar 1972;Smith 1978;Wube et al. 2008). What is additionally expected to contribute to the high intraspecific variation in birth season length in rodents and lagomorphs is their occurrence across a much larger geographical range (Schai-Braun and Hackländer 2016; Wilson et al. 2016;Smith et al. 2018) than, for instance, primates, which are mostly confined to the tropics .
Mid-latitude of the species' geographical range, as well as mean annual temperature and precipitation had strong effects on reproductive seasonality in lagomorphs. Previous studies in Rodentia (Heldstab 2021), Ruminantia (Zerbe et al. 2012), Carnivora (Heldstab et al. 2018), Primates (Di Bitetti and Janson 2000;Janson and Verdolin 2005;Heldstab et al. 2021) and Macroscelidea (Neal 1995) also found that species living at higher latitudes have shorter breeding seasons.
In the present study across lagomorphs, habitats that were warmer and had a higher amount of precipitation were associated with longer seasonal reproduction. Comparative studies in rodents (Tkadlec 2000;Heldstab 2021) and primates (Heldstab et al. 2018) also showed a positive effect of temperature on the length of the breeding season. Mean annual precipitation in a species' habitat did not, however, translate into a less seasonal birthing pattern in these two orders. The strictly herbivorous diet of the large majority of lagomorphs might be responsible for the stronger influence of precipitation on their birth seasonality. Diets of rodents and primates are much more variable including also fruits, seeds, nuts, meat, or insects in significant proportions not observed in lagomorphs (Wilman et al. 2014). The broader diet range in these two orders likely makes food availability and ultimately reproductive seasonality less dependent on precipitation; in contrast, the lagomorph diet is mostly forbs, herbs and grasses. In addition, since many lagomorphs, in particular, pikas, live at high altitudes , snow might also elicit a stronger influence of precipitation on reproductive seasonality in comparison to the other two orders inhabiting mostly low-altitudinal habitats.
Solitary species with large home ranges do not often encounter potential mates (Llewellyn and Enders 1954;Greenwald 1956). Females of such species are therefore expected to have long oestrous periods to alleviate the uncertainty of finding a male (Larivière and Ferguson 2003). This should lead to less seasonal birthing patterns. A study in ursids supports this prediction, although the observed effect of home range size was small (Spady et al. 2007). The present study, a study with a larger sample of Carnivora species (Heldstab et al. 2018), as well as comparative studies in rodents (Heldstab 2021) and primates  were all unable to detect an influence of home range size on seasonal reproduction.
Contrasting what has been reported in primates , carnivores (Tökölyi et al. 2014;Heldstab et al. 2018) and ruminants (Zerbe et al. 2012;Tökölyi et al. 2014) but as found in rodents (Heldstab 2021), a higher degree of reproductive seasonality was not linked to a shorter gestation length in lagomorphs in the present study. I suggest two reasons for the difference between rodents and lagomorphs in comparison to the other mammalian orders: their small body size and their post-partum oestrus allowing them to gestate again immediately following the birth of the young. Large animals in seasonal environments need to reduce the length of their reproductive cycle to squeeze it into a single year, so that they do not lose out on the reproductive opportunity in the next year (Kiltie 1988). Small mammals such as rodents and lagomorphs exhibit extremely rapid foetal and neonatal growth rates and thus already have very short gestation periods (Swihart 1984;Wilson et al. 2016). The post-partum oestrus shortens the time between lactation and the next pregnancy, enabling rodents and lagomorphs to fit the complete reproductive cycle in an even shorter time window. Although lagomorphs do not adjust gestation length, they still seem to adapt life history to seasonal habitats with short favourable time windows, as seasonal Leporid species do have fewer but larger litters. Selection in highly seasonal environments has also favoured fewer litters in rodents (Heldstab 2021), bats (Tökölyi et al. 2014), elephant shrews (Neal 1995), and Eulipotyphla (Symonds 1999), larger litter sizes in primates , earlier weaning age in rodents (Heldstab 2021) and seals (Gentry and Kooyman 2014), and older age at sexual maturity in Eulipotyphla (Symonds 1999) and Rodentia (Heldstab 2021).
As found in primates ) and carnivores (Heldstab et al. 2018), but unlike in ruminants (Zerbe et al. 2012) and rodents (Heldstab 2021), sexual size dimorphism in lagomorphs was not related to birth season length. Even in ruminants and rodents, sexual size dimorphism was also only marginally negatively associated with birth season length. The other allometric variable, body mass, also only had a small influence on lagomorph birth seasonality, as it was only correlated with birth season length measured in months but not in categories. Studies in primates (Plavcan et al. 2005;Heldstab et al. 2021), Carnivora (Heldstab et al. 2018), rodents (Heldstab 2021) and ruminants (Zerbe et al. 2012) demonstrate a similar weak impact of body mass on reproductive seasonality. In sum, habitat and life history characteristics seem to be much more important for lagomorph reproductive traits than allometric variables, as already suggested by Swihart (1984).
Acknowledgements I thank Marcus Clauss for constructive comments and gratefully acknowledge the many people and institutions that have contributed data to my compilation. I am also thankful to the two reviewers for their very useful comments on the manuscript.
Funding Open Access funding provided by Universität Zürich.
Data availability All data generated or analysed during this study are included in this published article and its supplementary information files.

Conflict of interest
The author declares that she does not have any conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.