Island demographics and trait associations in white-tailed deer.

When a population is isolated and composed of few individuals, genetic drift is the paramount evolutionary force and results in the loss of genetic diversity. Inbreeding might also occur, resulting in genomic regions that are identical by descent, manifesting as runs of homozygosity (ROHs) and the expression of recessive traits. Likewise, the genes underlying traits of interest can be revealed by comparing fixed SNPs and divergent haplotypes between affected and unaffected individuals. Populations of white-tailed deer (Odocoileus virginianus) on islands of Saint Pierre and Miquelon (SPM, France) have high incidences of leucism and malocclusions, both considered genetic defects; on the Florida Keys islands (USA) deer exhibit smaller body sizes, a polygenic trait. Here we aimed to reconstruct island demography and identify the genes associated with these traits in a pseudo case-control design. The two island populations showed reduced levels of genomic diversity and a build-up of deleterious mutations compared to mainland deer; there was also significant genome-wide divergence in Key deer. Key deer showed higher inbreeding levels, but not longer ROHs, consistent with long-term isolation. We identified multiple trait-related genes in ROHs including LAMTOR2 which has links to pigmentation changes, and NPVF which is linked to craniofacial abnormalities. Our mixed approach of linking ROHs, fixed SNPs and haplotypes matched a high number (~50) of a-priori body size candidate genes in Key deer. This suite of biomarkers and candidate genes should prove useful for population monitoring, noting all three phenotypes show patterns consistent with a complex trait and non-Mendelian inheritance.

Chromosome-level assembly of the Rangifer tarandus genome and validation of cervid and bovid evolution insights.

BACKGROUND: Genome assembly into chromosomes facilitates several analyses including cytogenetics, genomics and phylogenetics. Despite rapid development in bioinformatics, however, assembly beyond scaffolds remains challenging, especially in species without closely related well-assembled and available reference genomes. So far, four draft genomes of Rangifer tarandus (caribou or reindeer, a circumpolar distributed cervid species) have been published, but none with chromosome-level assembly. This emblematic northern species is of high interest in ecological studies and conservation since most populations are declining. RESULTS: We have designed specific probes based on Oligopaint FISH technology to upgrade the latest published reindeer and caribou chromosome-level genomes. Using this oligonucleotide-based method, we found six mis-assembled scaffolds and physically mapped 68 of the largest scaffolds representing 78% of the most recent R. tarandus genome assembly. Combining physical mapping and comparative genomics, it was possible to document chromosomal evolution among Cervidae and closely related bovids. CONCLUSIONS: Our results provide validation for the current chromosome-level genome assembly as well as resources to use chromosome banding in studies of Rangifer tarandus.

Genomic health is dependent on long-term population demographic history.

Current genetic methods of population assessment in conservation biology have been challenged by genome-scale analyses due to their quantitatively novel insights. These analyses include assessments of runs-of-homozygosity (ROH), genomic evolutionary rate profiling (GERP), and mutational load. Here, we aim to elucidate the relationships between these measures using three divergent ungulates: white-tailed deer, caribou, and mountain goat. The white-tailed deer is currently expanding, while caribou are in the midst of a significant decline. Mountain goats remain stable, having suffered a large historical bottleneck. We assessed genome-wide signatures of inbreeding using the inbreeding coefficient F and %ROH (FROH ) and identified evolutionarily constrained regions with GERP. Mutational load was estimated by identifying mutations in highly constrained elements (CEs) and sorting intolerant from tolerant (SIFT) mutations. Our results showed that F and FROH are higher in mountain goats than in caribou and white-tailed deer. Given the extended bottleneck and low Ne of the mountain goat, this supports the idea that the genome-wide effects of demographic change take time to accrue. Similarly, we found that mountain goats possess more highly constrained CEs and the lowest dN/dS values, both of which are indicative of greater purifying selection; this is also reflected by fewer mutations in CEs and deleterious mutations identified by SIFT. In contrast, white-tailed deer presented the highest mutational load with both metrics, in addition to dN/dS, while caribou were intermediate. Our results demonstrate that extended bottlenecks may lead to reduced diversity and increased FROH in ungulates, but not necessarily an increase in mutational load, probably due to the purging of deleterious alleles in small populations.

Birth date determines early calf survival in migratory caribou.

The decline of most caribou (Rangifer tarandus) populations underlines the need to understand the determinants of key demographic parameters. In migratory caribou, we have limited information on rates and drivers of pre-weaning mortality. We fitted 60 pregnant females of the Rivière-aux-Feuilles caribou herd with GPS camera collars to track the survival of calves from birth to weaning in 2016-2018. Over the three years, calf survival rate before weaning, i.e. to 01-Sep, approximately three months of age, was 0.63 (CI 0.50-0.77). Summer mortality risk was mainly influenced by calf birth date, with calves born earlier in the calving season having a lower mortality risk than those born later. Mortality also increased when calves experienced low or high temperature during calving. This study provides the first estimates of pre-weaning survival of migratory caribou calves in this herd, illustrating the value of new technologies to collect data otherwise difficult to obtain in widely distributed migratory populations. This approach can easily be extended to other large herbivores and predators. Our study brings new insights on how climate change may affect summer juvenile survival given the increased temperatures and faster changes in plant phenology expected in the future.

Fluctuating optimum and temporally variable selection on breeding date in birds and mammals.

Temporal variation in natural selection is predicted to strongly impact the evolution and demography of natural populations, with consequences for the rate of adaptation, evolution of plasticity, and extinction risk. Most of the theory underlying these predictions assumes a moving optimum phenotype, with predictions expressed in terms of the temporal variance and autocorrelation of this optimum. However, empirical studies seldom estimate patterns of fluctuations of an optimum phenotype, precluding further progress in connecting theory with observations. To bridge this gap, we assess the evidence for temporal variation in selection on breeding date by modeling a fitness function with a fluctuating optimum, across 39 populations of 21 wild animals, one of the largest compilations of long-term datasets with individual measurements of trait and fitness components. We find compelling evidence for fluctuations in the fitness function, causing temporal variation in the magnitude, but not the direction of selection. However, fluctuations of the optimum phenotype need not directly translate into variation in selection gradients, because their impact can be buffered by partial tracking of the optimum by the mean phenotype. Analyzing individuals that reproduce in consecutive years, we find that plastic changes track movements of the optimum phenotype across years, especially in bird species, reducing temporal variation in directional selection. This suggests that phenological plasticity has evolved to cope with fluctuations in the optimum, despite their currently modest contribution to variation in selection.

Design and validation of a 63K genome-wide SNP-genotyping platform for caribou/reindeer (Rangifer tarandus).

BACKGROUND: Development of large single nucleotide polymorphism (SNP) arrays can make genomic data promptly available for conservation problematic. Medium and high-density panels can be designed with sufficient coverage to offer a genome-wide perspective and the generated genotypes can be used to assess different genetic metrics related to population structure, relatedness, or inbreeding. SNP genotyping could also permit sexing samples with unknown associated metadata as it is often the case when using non-invasive sampling methods favored for endangered species. Genome sequencing of wild species provides the necessary information to design such SNP arrays. We report here the development of a SNP-array for endangered Rangifer tarandus using a multi-platform sequencing approach from animals found in diverse populations representing the entire circumpolar distribution of the species. RESULTS: From a very large comprehensive catalog of SNPs detected over the entire sample set (N = 894), a total of 63,336 SNPs were selected. SNP selection accounted for SNPs evenly distributed across the entire genome (~ every 50Kb) with known minor alleles across populations world-wide. In addition, a subset of SNPs was selected to represent rare and local alleles found in Eastern Canada which could be used for ecotype and population assignments - information urgently needed for conservation planning. In addition, heterozygosity from SNPs located in the X-chromosome and genotyping call-rate of SNPs located into the SRY gene of the Y-chromosome yielded an accurate and robust sexing assessment. All SNPs were validated using a high-throughput SNP-genotyping chip. CONCLUSION: This design is now integrated into the first genome-wide commercially available genotyping platform for Rangifer tarandus. This platform would pave the way to future genomic investigation of populations for this endangered species, including estimation of genetic diversity parameters, population assignments, as well as animal sexing from genetic SNP data for non-invasive samples.

Association of Environmental Factors with Seasonal Intensity of Erysipelothrix rhusiopathiae Seropositivity among Arctic Caribou.

Several caribou (Rangifer tarandus) populations have been declining concurrently with increases in infectious diseases in the Arctic. Erysipelothrix rhusiopathiae, a zoonotic bacterium, was first described in 2015 as a notable cause of illness and death among several Arctic wildlife species. We investigated epidemiologic and environmental factors associated with the seroprevalence of E. rhusiopathiae in the Arctic and found that seropositivity was highest during warmer months, peaking in September, and was highest among adult males. Summer seroprevalence increases tracked with the oestrid index from the previous year, icing and snowing events, and precipitation from the same year but decreased with growing degree days in the same year. Seroprevalence of E. rhusiopathiae varied more during the later years of the study. Our findings provide key insights into the influence of environmental factors on disease prevalence that can be instrumental for anticipating and mitigating diseases associated with climate change among Arctic wildlife and human populations.

Absence of founder effect and evidence for adaptive divergence in a recently introduced insular population of white-tailed deer (Odocoileus virginianus).

Islands are generally colonized by few individuals which could lead to a founder effect causing loss of genetic diversity and rapid divergence by strong genetic drift. Insular conditions can also induce new selective pressures on populations. Here, we investigated the extent of genetic differentiation within a white-tailed deer (Odocoileus virginianus) population introduced on an island and its differentiation with its source mainland population. In response to their novel environmental conditions, introduced deer changed phenotypically from mainland individuals, therefore we investigated the genetic bases of the morphological differentiation. The study was conducted on Anticosti Island (Québec, Canada) where 220 individuals were introduced 120 years ago, resulting in a population size over 160,000 individuals. We used genotyping-by-sequencing (GBS) to generate 8,518 filtered high-quality SNPs and compared patterns of genetic diversity and differentiation between the continental and Anticosti Island populations. Clustering analyses indicated a single panmictic island population and no sign of isolation by distance. Our results revealed a weak, albeit highly significant, genetic differentiation between the Anticosti Island population and its source population (mean FST  = 0.005), which allowed a population assignment success of 93%. Also, the high genetic diversity maintained in the introduced population supports the absence of a strong founder effect due to the large number of founders followed by rapid population growth. We further used a polygenic approach to assess the genetic bases of the divergent phenotypical traits between insular and continental populations. We found loci related to muscular function and lipid metabolism, which suggested that these could be involved in local adaptation on Anticosti Island. We discuss these results in a harvest management context.

Adult survival in migratory caribou is negatively associated with MHC functional diversity.

Genes of the major histocompatibility complex (MHC) are involved in acquired immunity in vertebrates. Only a few studies have investigated the fitness consequences of MHC gene diversity in wild populations. Here, we looked at the association between annual survival and body mass and MHC-DRB exon 2 (MHC-DRB) genetic diversity, obtained from high-throughput sequencing, in two declining migratory caribou (Rangifer tarandus) herds. To disentangle the potential direct and general effects of MHC-DRB genetic diversity, we compared different indices of diversity that were either based on DNA-sequence variation or on physicochemical divergence of the translated peptides, thereby covering a gradient of allelic-to-functional diversity. We found that (1) body mass was not related to MHC-DRB diversity or genotype, and (2) adult survival probability was negatively associated with point accepted mutation distance, a corrected distance that considers the likelihood of each amino acid substitution to be accepted by natural selection. In addition, we found no evidence of fluctuating selection over time on MHC-DRB diversity. We concluded that direct effects were involved in the negative relationship between MHC functional diversity and survival, although the mechanism underlying this result remains unclear. A possible explanation could be that individuals with higher MHC diversity suffer higher costs of immunity (immunopathology). Our results suggest that genetic diversity is not always beneficial even in genes that are likely to be strongly shaped by balancing selection.

Heritability Estimates of Antler and Body Traits in White-Tailed Deer (Odocoileus virginianus) From Genomic-Relatedness Matrices.

Estimating heritability (h2) is required to predict the response to selection and is useful in species that are managed or farmed using trait information. Estimating h2 in free-ranging populations is challenging due to the need for pedigrees; genomic-relatedness matrices (GRMs) circumvent this need and can be implemented in nearly any system where phenotypic and genome-wide single-nucleotide polymorphism (SNP) data are available. We estimated the heritability of 5 body and 3 antler traits in a free-ranging population of white-tailed deer (Odocoileus virginianus) on Anticosti Island, Quebec, Canada. We generated classic and robust GRMs from >10,000 SNPs: hind foot length, dressed body mass, and peroneus muscle mass had high h2 values of 0.62, 0.44, and 0.55, respectively. Heritability in male-only antler features ranged from 0.07 to 0.33. We explored the influence of filtering by minor allele frequency and data completion on h2: GRMs derived from fewer SNPs had reduced h2 estimates and the relatedness coefficients significantly deviated from those generated with more SNPs. As a corollary, we discussed limitations to the application of GRMs in the wild, notably how skewed GRMs, specifically many unrelated individuals, can increase variance around h2 estimates. This is the first study to estimate h2 on a free-ranging population of white-tailed deer and should be informative for breeding designs and management as these traits could respond to selection.

Maternal longevity and offspring sex in wild ungulates.

In species with sexual size dimorphism, the offspring of the larger sex usually have greater energy requirements and may lead to greater fitness costs for parents. The effects of offspring sex on maternal longevity, however, have only been tested in humans. Human studies produced mixed results and considerable debate mainly owing to the difficulty of distinguishing the effects of sexual dimorphism from sociocultural factors. To advance this debate, we examined how the relative number of sons influenced maternal longevity in four species of free-living ungulates (Soay sheep Ovis aries; bighorn sheep, Ovis canadensis; red deer, Cervus elaphus; mountain goat, Oreamnos americanus), with high male-biased sexual size dimorphism but without complicating sociocultural variables. We found no evidence for a higher cumulative cost of sons than of daughters on maternal longevity. For a given number of offspring, most females with many sons in all four populations lived longer than females with few sons. The higher cost of sons over daughters on maternal lifespan reported by some human studies may be the exception rather than the rule in long-lived iteroparous species.

No evidence of inbreeding depression in fast declining herds of migratory caribou.

Identifying inbreeding depression early in small and declining populations is essential for management and conservation decisions. Correlations between heterozygosity and fitness (HFCs) provide a way to identify inbreeding depression without prior knowledge of kinship among individuals. In Northern Quebec and Labrador, the size of two herds of migratory caribou (Rivière-George, RG and Rivière-aux-Feuilles, RAF) has declined by one to two orders of magnitude in the last three decades. This raises the question of a possible increase in inbreeding depression originating from, and possibly contributing to, the demographic decline in those populations. Here, we tested for the association of genomic inbreeding indices (estimated with 22,073 SNPs) with body mass and survival in 400 caribou sampled in RG and RAF herds between 1996 and 2016. We found no association of individual heterozygosity or inbreeding coefficient with body mass or annual survival. Furthermore, those genomic inbreeding indices remained stable over the period monitored. These results suggest that the rapid and intense demographic decline of the herds did not cause inbreeding depression in those populations. Although we found no evidence for HFCs, if demographic decline continues, it is possible that such inbreeding depression would be triggered.

Long-term studies of bighorn sheep and mountain goats reveal fitness costs of reproduction.

Fitness costs of reproduction are expected when resources are limited. Costs drive the evolution of life-history strategies and can affect population dynamics if females change their allocation of resources to reproduction. We studied fitness costs of reproduction in mountain ungulates in Alberta, Canada. We monitored two populations of bighorn sheep (Ovis canadensis) for 44 and 30 years, and one of mountain goats (Oreamnos americanus) for 30 years. Both species are highly iteroparous. Heterogeneity in individual reproductive potential makes fitness costs of reproduction difficult to detect and quantify without manipulations. In capital breeders, individual differences can be partly accounted for by considering body mass and other correlates of reproductive potential. Long-term monitoring can reveal costs that only manifest under stressful conditions such as disease or resource scarcity. Despite individual differences in reproductive potential, we detected fitness costs of reproduction in females. Costs, in terms of mass gain and survival, are almost entirely born by subsequent offspring, as mothers prioritize their own maintenance and survival. Costs are greater for primiparous females, decrease with increasing body mass and increase as resource availability declines, and sons are costlier than daughters. Costs may increase for senescent females that appear to reduce allocation to reproduction. In bighorn sheep, costs mostly involve reduced mass gain and lower survival of subsequent offspring. In mountain goats, costs include reductions in mass gain, subsequent fecundity and juvenile survival. In males, fitness costs derive mostly from attempts to reproduce rather than from siring success and likely depend upon individual competitiveness. In the absence of selective harvests, dominant males may enjoy high fitness and possibly lower costs compared to subordinates. The conservative reproductive tactic of mountain ungulate females likely explains why density dependence mostly involves later primiparity and lower recruitment, but rarely affects adult survival. Future research will seek to better account for heterogeneity in reproductive potential, assess cumulative reproductive costs and investigate the potential effects of fathers on maternal allocation tactics.

Les coûts de la reproduction émergent lorsque les ressources sont limitées et influencent l'évolution des stratégies d'histoire de vie. Si les femelles modifient l'allocation des ressources dans la reproduction pour mitiger ces coûts, ils peuvent aussi affecter la dynamique de population. Nous avons étudié les coûts de la reproduction chez des ongulés en Alberta, Canada. Nous avons suivi deux populations de mouflons d'Amérique (Ovis canadensis) pour 44 et 30 ans, et une de chèvre de montagne (Oreamnos americanus) pour 30 ans. Ces deux espèces sont itéropares. L'hétérogénéité dans le potentiel reproducteur des individus peut masquer les coûts en fitness liés à la reproduction. Ils sont donc plus souvent détectables à l'aide de manipulation expérimentale. Chez les espèces avec reproduction « sur capital ¼, il est possible de contrôler pour les différences individuelles en incluant des variables liées au potentiel reproducteur telles que la masse corporelle. Les suivis à long terme permettent de détecter des coûts qui se manifestent seulement lors de conditions environnementales défavorables. Malgré les différences en potentiel reproducteur, nous avons détecté des coûts chez les femelles. Ces coûts, tels qu'une réduction en croissance ou en survie, sont presque toujours subis par les jeunes nés lors d'épisodes de reproduction subséquents. En effet, les mères priorisent leur propre maintien et leur survie. Les coûts sont supérieurs pour les femelles primipares et dans les conditions défavorables. Ils sont plus faibles pour les femelles plus lourdes et il semble que les fils soient plus coûteux que les filles. Les coûts augmentent aussi chez les femelles sénescentes qui semblent réduire leur allocation dans la reproduction. Chez les mouflons, les coûts impliquent une réduction du gain en masse et une survie plus faible des jeunes nés lors de reproductions subséquentes. Chez les chèvres de montagne, ces coûts sont accompagnés d'une réduction de la fécondité future. Chez les mâles, les coûts sont surtout influencés par l'investissement en temps et en énergie dans le rut plutôt que par le succès d'accouplement. Ces coûts devraient donc surtout dépendre de l'habileté compétitive individuelle. En absence de chasse sélective, les mâles dominants devraient avoir un fort succès reproducteur avec des coûts potentiellement plus faibles comparativement aux subordonnées. La tactique conservative adoptée par les ongulés de montagne femelles peut expliquer pourquoi les effets densité-dépendant mènent à un retard dans l'âge à la primiparité et à un recrutement plus faible, mais affectent rarement la survie adulte. Nos recherches futures tenteront de mieux quantifier l'hétérogénéité individuelle, évalueront les coûts cumulatifs à la reproduction et investigueront l'effet potentiel des pères sur les tactiques d'allocation maternelles.

Reproductive plasticity of female white-tailed deer at high density and under harsh climatic conditions.

Life-history strategies of female ungulates usually depend on density-dependent and independent processes affecting body condition. Using a long-term data set on life-history traits of female white-tailed deer (2002-2014), we investigated the influence of population density and environmental factors on the reproductive effort of females. We also evaluated post-reproductive consequences on body condition using body mass, body fat, and body protein contents in the autumn following conception. We found that under high densities, females had a lower reproductive rate, which corresponds to a conservative reproduction strategy. However, females born at high density were more likely to reproduce and conceive larger litter size than females born at low density, a possible consequence of strong selective pressure in early life. Body condition was affected by reproduction; lactation had a large negative impact on body mass and body reserves, and conception, irrespectively of litter size, had a negative impact on body fat. Our long-term study demonstrates that plasticity in life-history strategies is a major determinant of reproductive potential for females living at high density and under harsh climates.

Faecal metabolites and hair cortisol as biological markers of HPA-axis activity in the Rocky mountain goat.

Monitoring glucocorticoids in faeces and hair is increasingly used in ecological studies and provides a powerful and minimally intrusive mean to identify physiological challenges faced by wild animals. Using a cortisol and a corticosterone immunoassays, we conducted an adrenocorticotropic (ACTH) challenge with five weekly repeated injections to validate the use of faecal glucocorticoid metabolites and hair cortisol concentration as biological markers of the HPA-axis activity in captive mountain goats (Oreamnos americanus). We also investigated the effect of endogenous (age, sex, reproductive status) and methodological (faecal sample collection date, freezing delay and hair type) variables on cortisol values using faecal and hair samples collected from marked wild mountain goats during a long-term study. The cortisol enzyme immunoassay was reliable for mountain goat faeces and hair, and was sensitive enough to detect a clear rise in glucocorticoid concentration following ACTH injections for both matrices. Age and sex had no detectable effect on faecal glucocorticoid metabolites, but hair cortisol concentration was higher in kids and yearlings than in older goats, and lower in adult males compared to adult females. Reproductive status had no detectable effect on both faecal and hair measurements. Faecal metabolite concentrations increased with sample collection date in late spring until mid-summer and decreased afterward until early fall. Guard hair had nearly twice as much cortisol per gram as undercoat hair. Prolonged delay to freezing reduced the concentration of faecal glucocorticoid metabolites, but degradation seemed limited when samples were exposed to wind and sun or when ambient temperature was low. We conclude that faeces and hair can be used as valid biomarkers of the HPA-axis activity in mountain goat provided that confounding variables are taken into account when interpreting measurements.

Influence of Early Reproductive Success on Longevity and Late Reproductive Success in an Alpine Ungulate.

The life-history theories of aging predict lifetime trade-offs between early reproductive allocation and late-life survival, reproduction, or both components of fitness. Recent studies in wild populations have found evidence for these early-late life trade-offs, but rarely have they been found across multiple traits while exploring the additional effects of variation in environmental conditions and individual quality. Benefiting from longitudinal data on adult female mountain goats (Oreamnos americanus), we investigated the influence of age at first reproduction (AFR) and early reproductive success (ERS) on longevity, late reproductive success, and senescence rates while accounting for the influence of natal environmental conditions and individual quality. Contrary to predictions, we did not find evidence for early-late life trade-offs. Instead, an earlier AFR and a greater ERS had positive but weak direct effects on late reproductive success. Natal population density, however, was the strongest determinant of all life-history traits, having a direct negative effect on female longevity, late reproductive success, AFR, and ERS. Although natal density reduced the probability of annual reproduction and annual survival during adulthood, higher allocation to reproduction in early life and poorer natal conditions did not lead to accelerated rates of senescence during adulthood. The results of this investigation provide an integrated picture of early-late life trade-offs, underscoring the importance of accounting for environmental conditions because of their potentially strong implications for population dynamics.

Experimental warming alters migratory caribou forage quality.

Global declines in caribou and reindeer (Rangifer) populations have drawn attention to the myriad of stressors that these Arctic and boreal forest herbivores currently face. Arctic warming has resulted in increased tundra shrub growth and therefore Rangifer forage quantity. However, its effects on forage quality have not yet been addressed although they may be critical to Rangifer body condition and fecundity. We investigated the impact of 8 yrs of summer warming on the quality of forage available to the Bathurst caribou herd using experimental greenhouses (n = 5) located in mesic birch hummock tundra in the central Canadian Low Arctic. Leaf forage quality and digestibility characteristics associated with nutrients (nitrogen and phosphorus), phenolics, and fiber were measured on the deciduous shrub Betula glandulosa (an important Rangifer diet component) at six time points through the growing season, and on five other very common vascular plant and lichen species in late summer. Experimental warming reduced B. glandulosa leaf nitrogen concentrations by ~10% in both late June and mid-July, but not afterwards. It also reduced late summer forage quality of the graminoid Eriophorum vaginatum by increasing phenolic concentrations 38%. Warming had mixed effects on forage quality of the lichen Cetraria cucullata in that it increased nutrient concentrations and tended to decrease fiber contents, but it also increased phenolics. Altogether, these warming-induced changes in forage quality over the growing season, and response differences among species, highlight the importance of Rangifer adaptability in diet selection. Furthermore, the early season reduction in B. glandulosa nitrogen content is a particular concern given the importance of this time for calf growth. Overall, our demonstration of the potential for significant warming impacts on forage quality at critical times for these animals underscores the importance of effective Rangifer range conservation to ensure sufficient appropriate habitat to support adaptability in forage selection in a rapidly changing environment.

Recognizing the importance of an all-inclusive approach to brown bear conservation now and into the future.

In their critique of our recent article in Oecologia (Pigeon et al. Oecologia 181:1101-1116, 2016a) investigating the influence of ambient temperature on brown bear habitat selection, Ordiz et al. (2017, current issue) argue that we downplay the role of human disturbance on bear behavior, and that we wrongly report on the findings of Ordiz et al. (Oecologia 166:59-67, 2011). We argue that our previous article in Oecologia (Pigeon et al. 2016a) by no means downplays the influence of human factors on bear behavior, and that we correctly stated that Ordiz et al. (2011) did not adequately consider the potential influence of temperature on their findings. Finally, we stress the relevance of considering all-inclusive approaches to the common goal of successful wildlife conservation.

Loss of connectivity among island-dwelling Peary caribou following sea ice decline.

Global warming threatens to reduce population connectivity for terrestrial wildlife through significant and rapid changes to sea ice. Using genetic fingerprinting, we contrasted extant connectivity in island-dwelling Peary caribou in northern Canada with continental-migratory caribou. We next examined if sea-ice contractions in the last decades modulated population connectivity and explored the possible impact of future climate change on long-term connectivity among island caribou. We found a strong correlation between genetic and geodesic distances for both continental and Peary caribou, even after accounting for the possible effect of sea surface. Sea ice has thus been an effective corridor for Peary caribou, promoting inter-island connectivity and population mixing. Using a time series of remote sensing sea-ice data, we show that landscape resistance in the Canadian Arctic Archipelago has increased by approximately 15% since 1979 and may further increase by 20-77% by 2086 under a high-emission scenario (RCP8.5). Under the persistent increase in greenhouse gas concentrations, reduced connectivity may isolate island-dwelling caribou with potentially significant consequences for population viability.

Staying cool in a changing landscape: the influence of maximum daily ambient temperature on grizzly bear habitat selection.

To fulfill their needs, animals are constantly making trade-offs among limiting factors. Although there is growing evidence about the impact of ambient temperature on habitat selection in mammals, the role of environmental conditions and thermoregulation on apex predators is poorly understood. Our objective was to investigate the influence of ambient temperature on habitat selection patterns of grizzly bears in the managed landscape of Alberta, Canada. Grizzly bear habitat selection followed a daily and seasonal pattern that was influenced by ambient temperature, with adult males showing stronger responses than females to warm temperatures. Cutblocks aged 0-20 years provided an abundance of forage but were on average 6 °C warmer than mature conifer stands and 21- to 40-year-old cutblocks. When ambient temperatures increased, the relative change (odds ratio) in the probability of selection for 0- to 20-year-old cutblocks decreased during the hottest part of the day and increased during cooler periods, especially for males. Concurrently, the probability of selection for 21- to 40-year-old cutblocks increased on warmer days. Following plant phenology, the odds of selecting 0- to 20-year-old cutblocks also increased from early to late summer while the odds of selecting 21- to 40-year-old cutblocks decreased. Our results demonstrate that ambient temperatures, and therefore thermal requirements, play a significant role in habitat selection patterns and behaviour of grizzly bears. In a changing climate, large mammals may increasingly need to adjust spatial and temporal selection patterns in response to thermal constraints.