From gestation to weaning: Combining robust design and multi‐event models unveils cost of lactation in a large herbivore

Reduced habitat quality increases intrinsic but not ecological costs of reproduction

Although the costs of reproduction are predicted to vary with the quality of the breeding habitat thereby affecting population dynamics and life‐history trade‐offs, empirical evidence for this pattern remains sparse and equivocal. Costs of reproduction can operate through immediate ecological mechanisms or through delayed intrinsic mechanisms. Ignoring these separate pathways might hinder the identification of costs and the understanding of their consequences. We experimentally investigated the survival costs of reproduction for adult little owls (Athene noctua) within a gradient of habitat quality. We supplemented food to nestlings, thereby relieving the parents’ effort for brood provisioning. We used radio‐tracking and Bayesian multistate modeling based on marked recapture and dead recovery to estimate survival rates of adult little owls across the year as a function of food supplementation and habitat characteristics. Food supplementation to nestlings during the breeding season increased parental survival not only during the breeding season but also during the rest of the year. Thus, the low survival of parents of unfed broods likely represents both, strong ecological and strong intrinsic costs of reproduction. However, while immediate ecological costs occurred also in high‐quality habitats, intrinsic costs carrying over to the post‐breeding period occurred only in low‐quality habitats. Our results suggest that immediate costs resulting from ecological mechanisms such as predation, are high also in territories of high habitat quality. Long‐term costs resulting from intrinsic trade‐offs, however, are only paid in low‐quality habitats. Consequently, differential effects of habitat quality on immediate ecological and delayed intrinsic mechanisms can mask the increase of costs of reproduction in low‐quality breeding habitats. Intrinsic costs may represent an underrated mechanism of habitat quality affecting adult survival rate thereby considerably accelerating population decline in degrading habitats. This study therefore highlights the need for a long‐term perspective to fully assess the costs of reproduction and the role of habitat quality in modifying these costs.

Positive effect of spring advance on the diet quality of an alpine herbivore

Changes in vegetation phenology related to global warming are having alarming effects on the life history traits of many herbivore species. Such changes are particularly critical in alpine ecosystems, where strong climate limitations on plant growth make seasonal synchronization imperative for the growth, reproduction and survival of herbivores. However, despite the pivotal role of resource-use strategies on the performances of such species, few studies have explicitly assessed the mechanistic impact of climate change on their diets. We aimed to fill this gap by studying the effect of spring onset on the dietary composition and quality of a medium-size alpine herbivore while considering density-dependent processes and age- and sex-specific differences in foraging behavior. Using an exceptional, long-term (24 years) direct individual-based dietary monitoring of a Pyrenean chamois population (Rupicapra pyrenaica pyrenaica), we showed that ongoing earlier onsets of spring are leading to an earlier access to high-quality forage and therefore a higher diet quality at a fixed date, without apparent changes in diet composition. We also showed that at high densities, intraspecific competition reduced diet quality by driving animals to feed more on woody plants and less on nutritious forbs and graminoids. By assessing the mechanistic effects of global warming on the dietary patterns of species at the center of trophic networks, this study is an essential step for predictive models aiming at understanding the ongoing ecosystem consequences of the global climatic crisis.

Maternal effects on reproduction in the precocial European hare (Lepus europaeus)

In female mammals, reproduction, and in particular lactation, is the energetically most exigent life-history phase. Reproduction is strongly controlled by body reserves and food availability, so females with better body condition or food supply are believed to have higher reproductive output. Additionally, the growth and mortality of young mammals depends on their postnatal development. Therefore, the degree of precociality affects energetic demands for both mothers and young. To study the reproductive performance of the precocial European hare (Lepus europaeus), we analysed relationships between six predictor variables describing maternal and environmental effects and nine response variables relating to reproduction from 217 captive females. We compared the data with those of precocial and altricial mammal species from an extensive literature search. For hares, we found: (1) Heavier females had heavier litters at birth. (2) In summer and spring, total litter mass was larger than in winter. (3) At the end of lactation, the litters of multiparous females were heavier than those of primiparous females. (4) Both older females and females giving birth for the first time had relatively high leveret mortality during lactation. Comparing our results with the literature for other mammals revealed that the body condition (i.e., body mass) of females before birth is predictive of reproductive parameters in both precocial and altricial species. In the precocial hare, female body condition is no longer predictive of reproductive parameters at the end of lactation, whereas in altricial species, female body condition remains predictive of reproduction (litter mass at the end of lactation, offspring mortality) until the end of lactation. We conclude that these effects are caused by precocial offspring feeding on solid food soon after birth and, thus, being less dependent on the mother's body condition during lactation than altricial offspring. In line with this, precociality might have evolved as a way of buffering offspring against maternal effects.

Anonymous fecal sampling and NIRS studies of diet quality: Problem or opportunity?

Investigating the drivers of diet quality is a key issue in wildlife ecology and conservation. Fecal near infrared reflectance spectroscopy (f-NIRS) is widely used to assess dietary quality since it allows for noninvasive, rapid, and low-cost analysis of nutrients. Samples for f-NIRS can be collected and analyzed with or without knowledge of animal identities. While anonymous sampling allows to reduce the costs of individual identification, as it neither requires physical captures nor DNA genotyping, it neglects the potential effects of individual variation. As a consequence, regression models fitted to investigate the drivers of dietary quality may suffer severe issues of pseudoreplication. I investigated the relationship between crude protein and ecological predictors at different time periods to assess the level of individual heterogeneity in diet quality of 22 marked chamois Rupicapra rupicapra monitored over 2 years. Models with and without individual grouping effect were fitted to simulate identifiable and anonymous fecal sampling, and model estimates were compared to evaluate the consequences of anonymizing data collection and analysis. The variance explained by the individual random effect and the value of diet repeatability varied with seasons and peaked in winter. Despite the occurrence of individual variation in dietary quality, ecological parameter estimates under identifiable or anonymous sampling were consistently similar. This study suggests that anonymous fecal sampling may provide robust estimates of the relationship between dietary quality and ecological correlates. However, since the level of individual heterogeneity in dietary quality may vary with species- or study-specific features, inconsequential pseudoreplication should not be assumed in other taxa. When individual differences are known to be inconsequential, anonymous sampling allows to optimize the trade-off between sampling intensity and representativeness. When pseudoreplication is consequential, however, no conclusive remedy exists to effectively resolve nonindependence.

Multiple nutritional currencies shape pregnancy in a large herbivore

Nutritional condition embodies environmental conditions experienced by animals with survival and reproductive consequences. Body fat is often associated with ungulate fecundity; however, other nutritional currencies may limit fecundity. Using data from 129 moose (Alces alces (Linnaeus, 1758)) monitored over 429 moose-years, we examined the limiting role of multiple nutritional currencies on pregnancy rates while concurrently assessing the influence of age and prior reproduction. Females tended to be pregnant in successive years, suggesting differences in individual or habitat quality. However, probability of pregnancy declined with survival of calves from prior litters, indicating a reproductive cost to rearing offspring. Pregnancy was positively associated with serum iron (Fe), body fat, body mass, and serum protein. The best model predicting pregnancy included serum Fe, body fat, and age class, with serum Fe being the strongest single predictor. Moose with Fe concentrations considered deficient in cattle (Bos taurus Linnaeus, 1758) had pregnancy rates of 33%–35%, whereas 89%–91% of individuals with sufficient Fe were pregnant. We subsequently evaluated hypotheses concerning factors potentially limiting Fe concentrations, including Fe-deficient diet, chronic infection, parasitism, and malnutrition. The best supported hypothesis was energy and protein malnutrition constrained Fe stores. We conclude that subclinical anemia due to nutritional constraints can limit or indicate limits in moose fecundity.

In the wrong place at the wrong time: Moose and deer movement patterns influence wildlife-vehicle collision risk

Mitigation strategies for wildlife-vehicle collisions require sufficient knowledge about why, where and when collisions occur in order to be an efficient tool to improve public safety. Collisions with cervids are known to be influenced by spatial factors such as topography and forest cover. However, temporal changes in animal and motorist behaviors are often overlooked although they can increase the odds of cervid-vehicle collisions. Consequently, we evaluated potential factors influencing the spatiotemporal distribution of 450 collisions with moose and white-tailed deer that occurred between 1990 and 2015 along the 100-km long highway in southeastern Québec, Canada. Both spatial and temporal factors efficiently explained moose-vehicle collisions but not collisions with white-tailed deer, suggesting that the latter occurred more randomly along the highway. The risk of moose-vehicle collisions was mainly modulated by topographic and habitat variables, as the interactions between slope and elevation and slope and distance to suitable moose habitats had a strong effect on collision risk. Road sinuosity and the proportion of mature coniferous stands around the collision site positively influenced deer-vehicle collisions. A temporal increase in collision numbers was noted in different biological periods during which movement rates are known to be higher (e.g. post-winter dispersal and rut). These results suggest that cervid movement is the main factor influencing collision risk and frequency. Our results indicate that mitigation strategies aimed at decreasing the probability of collision with cervids must be species-specific and should focus more closely on animal movement.

Joint estimation of survival and breeding probability in female dolphins and calves with uncertainty in state assignment

While the population growth rate in long-lived species is highly sensitive to adult survival, reproduction can also significantly drive population dynamics. Reproductive parameters can be challenging to estimate as breeders and nonbreeders may vary in resighting probability and reproductive status may be difficult to assess. We extended capture-recapture (CR) models previously fitted for data on other long-lived marine mammals to estimate demographic parameters while accounting for detection heterogeneity between individuals and state uncertainty regarding reproductive status. We applied this model to data on 106 adult female bottlenose dolphins observed over 13 years. The detection probability differed depending on breeding status. Concerning state uncertainty, offspring were not always sighted with their mother, and older calves were easier to detect than young-of-the-year (YOY), respectively, 0.79 (95% CI 0.59-0.90) and 0.58 (95% CI 0.46-0.68). This possibly led to inaccurate reproductive status assignment of females. Adult female survival probability was high (0.97 CI 95% 0.96-0.98) and did not differ according to breeding status. Young-of-the-year and 1-year-old calves had a significantly higher survival rate than 2-year-old (respectively, 0.66 CI 95% 0.50-0.78 and 0.45 CI 95% 0.29-0.61). This reduced survival is probably related to weaning, a period during which young are exposed to more risks since they lose protection and feeding from the mother. The probability of having a new YOY was high for breeding females that had raised a calf to the age of 3 or lost a 2-year-old calf (0.71, CI 95% 0.45-0.88). Yet, this probability was much lower for nonbreeding females and breeding females that had lost a YOY or a 1-year-old calf (0.33, 95% CI 0.26-0.42). The multievent CR framework we used is highly flexible and could be easily modified for other study questions or taxa (marine or terrestrial) aimed at modeling reproductive parameters.

Revisiting the functional response in habitat selection for large herbivores: a matter of spatial variation in resource distribution?

Most habitats are distributed heterogeneously in space, forcing animals to move according to both habitat characteristics and their needs for energy and safety. Animal space use should therefore vary according to habitat characteristics, a process known as the “functional response” in habitat selection. This response has often been tested vis-à-vis the proportion of a habitat category within areas available to individuals. Measuring sought-after resources in landscape where they are continuously distributed is a challenge and we posit here that both the mean availability of a resource and its spatial variation should be measured. Accordingly, we tested for a functional response in habitat selection according to these two descriptors of the resource available for a mountain herbivore. We hypothesized that selection should decrease with mean value of resources available and increase with its spatial variation. Based on GPS data from 50 chamois females and data on the actual foodscape (i.e., distribution of edible-only biomass in the landscape), we estimated individual selection ratio (during summer months) for biomass at the home range level, comparing edible biomass in individual home ranges and the mean and standard deviation of edible biomass in their available range. Chamois being a group-living species, available accessible ranges were shared by several individuals that formed socio-spatial groups (clusters) in the population. As expected, selection ratios increased with the standard deviation of edible resources in each cluster, but unlike our prediction, was unrelated to its mean. Selection of areas richer in resources hence did not fade away when more resources were available on average, a result that may be explained by the need for this capital breeder species to accumulate fat-reserve at a high rate during summer months. Low spatial variation could limit the selection of chamois, which highlights the importance of resource distribution in the process of habitat selection.

Detection of tick-borne pathogens in questing Ixodes ricinus in the French Pyrenees and first identification of Rickettsia monacensis in France

Ticks are important vectors of several human and animal pathogens. In this study, we estimated the prevalence of important tick-borne infections in questing ticks from an area in Southwestern France (Hautes-Pyrénées) inhabited by Pyrenean chamois (Rupicapra pyrenaica pyrenaica) experiencing high tick burden. We examined adult and nymph ticks collected by the flag dragging method from 8 to 15 sites in the Pic de Bazès during the years 2009, 2011, 2013 and 2015. PCR assays were conducted on selected ticks for the detection of Borrelia burgdorferi s.l., Babesia spp., Rickettsia spp., spotted fever group (SFG) Rickettsia and Anaplasma phagocytophilum. Randomly selected positive samples were submitted for sequence analysis. A total of 1971 questing ticks were collected including 95 males, 101 females and 1775 nymphs. All collected ticks were identified as Ixodes ricinus. Among them, 696 ticks were selected for pathogen detection and overall prevalence was 8.4% for B. burgdorferi s.l.; 0.4% for Babesia spp.; 6.1% for A. phagocytophilum; 17.6% for Rickettsia spp.; and 8.1% for SFG Rickettsia. Among the sequenced pathogens, we detected in this population of ticks the presence of Babesia sp. EU1 and Rickettsia helvetica, as well as Rickettsia monacensis for the first time in France. The detection of these pathogens in the Pic de Bazès highlights the potential infection risks for visitors to this area and the Pyrenean chamois population.

Making use of multiple surveys: Estimating breeding probability using a multievent-robust design capture-recapture model

Increased environmental stochasticity due to climate change will intensify temporal variance in the life-history traits, and especially breeding probabilities, of long-lived iteroparous species. These changes may decrease individual fitness and population viability and is therefore important to monitor. In wild animal populations with imperfect individual detection, breeding probabilities are best estimated using capture-recapture methods. However, in many vertebrate species (e.g., amphibians, turtles, seabirds), nonbreeders are unobservable because they are not tied to a territory or breeding location. Although unobservable states can be used to model temporary emigration of nonbreeders, there are disadvantages to having unobservable states in capture-recapture models. The best solution to deal with unobservable life-history states is therefore to eliminate them altogether. Here, we achieve this objective by fitting novel multievent-robust design models which utilize information obtained from multiple surveys conducted throughout the year. We use this approach to estimate annual breeding probabilities of capital breeding female elephant seals (Mirounga leonina). Conceptually, our approach parallels a multistate version of the Barker/robust design in that it combines robust design capture data collected during discrete breeding seasons with observations made at other times of the year. A substantial advantage of our approach is that the nonbreeder state became "observable" when multiple data sources were analyzed together. This allowed us to test for the existence of state-dependent survival (with some support found for lower survival in breeders compared to nonbreeders), and to estimate annual breeding transitions to and from the nonbreeder state with greater precision (where current breeders tended to have higher future breeding probabilities than nonbreeders). We used program E-SURGE (2.1.2) to fit the multievent-robust design models, with uncertainty in breeding state assignment (breeder, nonbreeder) being incorporated via a hidden Markov process. This flexible modeling approach can easily be adapted to suit sampling designs from numerous species which may be encountered during and outside of discrete breeding seasons.