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.

Selection of summer feeding sites and food resources by female migratory caribou (Rangifer tarandus) determined using camera collars.

Migratory caribou (Rangifer tarandus) is a socioeconomically and culturally key species for northern communities in the Arctic, and most of its populations are experiencing a sharp decline. Female migratory caribou depend on the availability of summer habitat resources to meet the needs associated with lactation and the accumulation of fat reserves to survive when resources are less abundant. Because of the large scales at which habitat and resource data are usually available, information on how female migratory caribou select habitat and resources at fine scales in the wild is lacking. To document selection of summer feeding sites, we equipped 60 female caribou with camera collars from 2016 to 2018. We collected a total of 65,150 10-sec videos between June 1st and September 1st for three years with contrasted spring phenology. We determined the selection at the feeding site scale (3rd scale of Johnson) and food item scale (4th scale of Johnson) using resource selection probability functions. Wetlands were highly selected as feeding sites in June and July while they were avoided in August. Shrublands were mostly selected in July and August. At the resources scale, lichen, birch, willow, and mushrooms were the most strongly selected resources. Our results provide precise and novel information on habitat selection at feeding sites and food resources selected by female caribou in the wild. This information will help understand foraging patterns and habitat selection behavior of female migratory caribou and will contribute to the management and conservation of its declining 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.

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.

Demographic history and conservation genomics of caribou (Rangifer tarandus) in Québec.

The loss of genetic diversity is a challenge many species are facing, with genomics being a potential tool to inform and prioritize decision-making. Most caribou (Rangifer tarandus) populations have experienced significant recent declines throughout Québec, Canada, and are considered of concern, threatened or endangered. Here, we calculated the ancestral and contemporary patterns of genomic diversity of five representative caribou populations and applied a comparative population genomics framework to assess the interplay between demographic events and genomic diversity. We first calculated a caribou specific mutation rate, µ, by extracting orthologous genes from related ungulates and estimating the rate of synonymous mutations. Whole genome re-sequencing was then completed on 67 caribou: from these data we calculated nucleotide diversity, θ π and estimated the coalescent or ancestral effective population size (N e), which ranged from 12,030 to 15,513. When compared to the census size, N C, the endangered Gaspésie Mountain caribou population had the highest ancestral N e:N C ratio which is consistent with recent work suggesting high ancestral N e:N C is of conservation concern. In contrast, values of contemporary N e, estimated from linkage-disequilibrium, ranged from 11 to 162, with Gaspésie having among the highest contemporary N e:N C ratio. Importantly, classic conservation genetics theory would predict this population to be of less concern based on this ratio. Interestingly, F varied only slightly between populations, and despite evidence of bottlenecks across the province, runs of homozygosity were not abundant in the genome. Tajima's D estimates mirrored the demographic models and current conservation status. Our study highlights how genomic patterns are nuanced and potentially misleading if viewed only through a contemporary lens; we argue a holistic conservation genomics view should integrate ancestral N e and Tajima's D into management decisions.

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.

Large herbivores trigger spatiotemporal changes in forest plant diversity.

Large herbivores can exert top-down control on terrestrial plant communities, but the magnitude, direction, and scale dependency of their impacts remain equivocal, especially in temperate and boreal forests, where multiple disturbances often interact. Using a unique, long-term, and replicated landscape experiment, we assessed the influence of a high density of white-tailed deer (Odocoileus virginianus) on the spatiotemporal dynamics of diversity, composition, and successional trajectories of understorey plant assemblages in recently logged boreal forests. This experiment provided a rare opportunity to test whether deer herbivory represents a direct filter on plant communities or if it mainly acts to suppress dominant plants, which, in turn, release other plant species from strong negative plant-plant interactions. These two hypotheses make different predictions about changes in community composition and alpha and beta diversity in different vegetation layers and at different spatial scales. Our results showed that deer had strong effects on plant community composition and successional trajectories, but the resulting impacts on plant alpha and beta diversity patterns were markedly scale dependent in both time and space. Responses of tree and non-tree vegetation layers were strongly asymmetric. Deer acted both as a direct filter and as a suppressor of dominant plant species during early forest succession, but the magnitude of both processes was specific to tree and non-tree vegetation layers. Although our data supported the ungulate-driven homogenization hypothesis, compositional shifts and changes of alpha diversity were poor predictors of beta diversity loss. Our findings underscore the importance of long-term studies in revealing nonlinear temporal community trends, and they challenge managers to prioritize particular community properties and scales of interest, given contrasting trends of composition and alpha and beta diversity across spatial scales.

CNVs with adaptive potential in Rangifer tarandus: genome architecture and new annotated assembly.

Rangifer tarandus has experienced recent drastic population size reductions throughout its circumpolar distribution and preserving the species implies genetic diversity conservation. To facilitate genomic studies of the species populations, we improved the genome assembly by combining long read and linked read and obtained a new highly accurate and contiguous genome assembly made of 13,994 scaffolds (L90 = 131 scaffolds). Using de novo transcriptome assembly of RNA-sequencing reads and similarity with annotated human gene sequences, 17,394 robust gene models were identified. As copy number variations (CNVs) likely play a role in adaptation, we additionally investigated these variations among 20 genomes representing three caribou ecotypes (migratory, boreal and mountain). A total of 1,698 large CNVs (length > 1 kb) showing a genome distribution including hotspots were identified. 43 large CNVs were particularly distinctive of the migratory and sedentary ecotypes and included genes annotated for functions likely related to the expected adaptations. This work includes the first publicly available annotation of the caribou genome and the first assembly allowing genome architecture analyses, including the likely adaptive CNVs reported here.

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.

Combining stable isotopes, morphological, and molecular analyses to reconstruct the diet of free-ranging consumers.

Accurate estimates of animal diet composition are essential to untangle complex interactions in food webs. Biomarkers and molecular tools are increasingly used to estimate diet, sometimes alongside traditional dietary tracing methods. Yet only a few empirical studies have compared the outcomes and potential gains of using a combination of these methods, especially using free-ranging animals with distinct foraging preferences.We used stable isotopes, morphological, and molecular analyses to investigate the diet of free-ranging consumers with two distinct diet types, that is, carnivore and omnivore. By combining the three analytical methods to assess the diet of consumers during the same period, we aimed to identify the limits of each method and to assess the potential benefits of their combined use to derive diet estimates.Our results showed that the different methods led to a consistent diet description for carnivores, which have a relatively simple diet mixture, but their outcomes somewhat differed for omnivore, which have a more complex diet. Still, the combined use of morphological and molecular analyses enhanced the diversity of food sources detected compared to the use of a single method independently of diet types. Precision of diet estimates derived from stable isotope analyses was improved by the addition of priors obtained from morphological and molecular diet analyses of the same population.Although we used free-ranging animals without a known diet, our empirical testing of three of the most widely used methods of diet determination highlights the limits of relying over a single approach, especially in systems with few or no a priori information about the foraging habits of consumers. The choice of an appropriate approach of diet description should be a key step when planning dietary studies of free-ranging populations. We recommend using more than one dietary determination methods especially for species with complex diet mixtures.

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.

Muskox status, recent variation, and uncertain future.

Muskoxen (Ovibos moschatus) are an integral component of Arctic biodiversity. Given low genetic diversity, their ability to respond to future and rapid Arctic change is unknown, although paleontological history demonstrates adaptability within limits. We discuss status and limitations of current monitoring, and summarize circumpolar status and recent variations, delineating all 55 endemic or translocated populations. Acknowledging uncertainties, global abundance is ca 170 000 muskoxen. Not all populations are thriving. Six populations are in decline, and as recently as the turn of the century, one of these was the largest population in the world, equaling ca 41% of today's total abundance. Climate, diseases, and anthropogenic changes are likely the principal drivers of muskox population change and result in multiple stressors that vary temporally and spatially. Impacts to muskoxen are precipitated by habitat loss/degradation, altered vegetation and species associations, pollution, and harvest. Which elements are relevant for a specific population will vary, as will their cumulative interactions. Our summaries highlight the importance of harmonizing existing data, intensifying long-term monitoring efforts including demographics and health assessments, standardizing and implementing monitoring protocols, and increasing stakeholder engagement/contributions.

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.

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.

Space, time and captivity: quantifying the factors influencing the fecal microbiome of an alpine ungulate.

The community of microorganisms in the gut is affected by host species, diet and environment and is linked to normal functioning of the host organism. Although the microbiome fluctuates in response to host demands and environmental changes, there are core groups of microorganisms that remain relatively constant throughout the hosts lifetime. Ruminants are mammals that rely on highly specialized digestive and metabolic modifications, including microbiome adaptations, to persist in extreme environments. Here, we assayed the fecal microbiome of four mountain goat (Oreamnos americanus) populations in western North America. We quantified fecal microbiome diversity and composition among groups in the wild and captivity, across populations and in a single group over time. There were no differences in community evenness or diversity across groups, although we observed a decreasing diversity trend across summer months. Pairwise sample estimates grouped the captive population distinctly from the wild populations, and moderately grouped the southern wild group distinctly from the two northern wild populations. We identified 33 genera modified by captivity, with major differences in key groups associated with cellulose degradation that likely reflect differences in diet. Our findings are consistent with other ruminant studies and provide baseline microbiome data in this enigmatic species, offering valuable insights into the health of wild alpine ungulates.