Patterns, sources, and consequences of variation in age-specific vital rates: Insights from a long-term study of Weddell seals.

Variations in the reproductive and survival abilities of individuals within a population are ubiquitous in nature, key to individual fitness, and affect population dynamics, which leads to strong interest in understanding causes and consequences of vital-rate variation. For long-lived species, long-term studies of large samples of known-age individuals are ideal for evaluating vital-rate variation. A population of Weddell seals in Erebus Bay, Antarctica, has been studied each Austral spring since the 1960s. Since 1982, all newborns have been tagged each year and multiple capture-mark-recapture (CMR) surveys have been conducted annually. Over the past 20 years, a series of analyses have built on results of earlier research by taking advantage of steady improvements in the project's long-term CMR data and available analytical methods. Here, I summarize progress made on four major topics related to variation in age-specific vital rates for females: early-life survival and age at first reproduction, costs of reproduction, demographic buffering, and demographic senescence. Multistate modelling found that age at first reproduction varies widely (4-14 years of age) and identified contrasting influences of maternal age on survival and recruitment rates of offspring. Subsequent analyses of data for females after recruitment revealed costs of reproduction to both survival and future reproduction and provided strong evidence of demographic buffering. Recent results indicated that important levels of among-individual variation exist in vital rates and revealed contrasting patterns for senescence in reproduction and survival. Sources of variation in vital rates include age, reproductive state, year, and individual. The combination of luck and individual quality results in strong variation in individual fitness outcomes: ~80% of females born in the population produce no offspring, and the remaining 20% vary strongly in lifetime reproductive output (range: 1-23 pups). Further research is needed to identify the specific environmental conditions that lead to annual variation in vital rates and to better understand the origins of individual heterogeneity. Work is also needed to better quantify the relative roles of luck, maternal effects, and environmental conditions on variation in vital rates and to learn the importance of such variation to demographic performance of offspring and on overall population dynamics.

Evaluating the importance of individual heterogeneity in reproduction to Weddell seal population dynamics using integral projection models.

Identifying and accounting for unobserved individual heterogeneity in vital rates in demographic models is important for estimating population-level vital rates and identifying diverse life-history strategies, but much less is known about how this individual heterogeneity influences population dynamics. We aimed to understand how the distribution of individual heterogeneity in reproductive and survival rates influenced population dynamics using vital rates from a Weddell seal population by altering the distribution of individual heterogeneity in reproduction, which also altered the distribution of individual survival rates through the incorporation of our estimate of the correlation between the two rates and assessing resulting changes in population growth. We constructed an integral projection model (IPM) structured by age and reproductive state using estimates of vital rates for a long-lived mammal that has recently been shown to exhibit large individual heterogeneity in reproduction. Using output from the IPM, we evaluated how population dynamics changed with different underlying distributions of unobserved individual heterogeneity in reproduction. Results indicate that the changes to the underlying distribution of individual heterogeneity in reproduction cause very small changes in the population growth rate and other population metrics. The largest difference in the estimated population growth rate resulting from changes to the underlying distribution of individual heterogeneity was less than 1%. Our work highlights the differing importance of individual heterogeneity at the population level compared to the individual level. Although individual heterogeneity in reproduction may result in large differences in the lifetime fitness of individuals, changing the proportion of above- or below-average breeders in the population results in much smaller differences in annual population growth rate. For a long-lived mammal with stable and high adult-survival that gives birth to a single offspring, individual heterogeneity in reproduction has a limited effect on population dynamics. We posit that the limited effect of individual heterogeneity on population dynamics may be due to canalization of life-history traits.

Evidence of an Absence of Inbreeding Depression in a Wild Population of Weddell Seals (Leptonychotes weddellii).

Inbreeding depression can reduce the viability of wild populations. Detecting inbreeding depression in the wild is difficult; developing accurate estimates of inbreeding can be time and labor intensive. In this study, we used a two-step modeling procedure to incorporate uncertainty inherent in estimating individual inbreeding coefficients from multilocus genotypes into estimates of inbreeding depression in a population of Weddell seals (Leptonychotes weddellii). The two-step modeling procedure presented in this paper provides a method for estimating the magnitude of a known source of error, which is assumed absent in classic regression models, and incorporating this error into inferences about inbreeding depression. The method is essentially an errors-in-variables regression with non-normal errors in both the dependent and independent variables. These models, therefore, allow for a better evaluation of the uncertainty surrounding the biological importance of inbreeding depression in non-pedigreed wild populations. For this study we genotyped 154 adult female seals from the population in Erebus Bay, Antarctica, at 29 microsatellite loci, 12 of which are novel. We used a statistical evidence approach to inference rather than hypothesis testing because the discovery of both low and high levels of inbreeding are of scientific interest. We found evidence for an absence of inbreeding depression in lifetime reproductive success, adult survival, age at maturity, and the reproductive interval of female seals in this population.

Sources of variation in maternal allocation in a long-lived mammal.

Life history theory predicts allocation of energy to reproduction varies with maternal age, but additional maternal features may be important to the allocation of energy to reproduction. We aimed to characterize age-specific variation in maternal allocation and assess the relationship between maternal allocation and other static and dynamic maternal features. Mass measurements of 531 mothers and pups were used with Bayesian hierarchical models to explain the relationship between diverse maternal attributes and both the proportion of mass allocated by Weddell seal mothers, and the efficiency of mass transfer from mother to pup during lactation as well as the weaning mass of pups. Our results demonstrated that maternal mass was strongly and positively associated with the relative reserves allocated by a mother and a pup's weaning mass but that the efficiency of mass transfer declines with maternal parturition mass. Birthdate was positively associated with proportion mass allocation and pup weaning mass, but mass transfer efficiency was predicted to be highest at the mean birthdate. The relative allocation of maternal reserves declined with maternal age but the efficiency of mass transfer to pups increases, suggestive of selective disappearance of poor-quality mothers. These findings highlight the importance of considering multiple maternal features when assessing variation in maternal allocation.

Maladaptive nest-site selection by a sagebrush dependent species in a grazing-modified landscape.

Animals are expected to select habitats that maximize their fitness over evolutionary time scales. Yet in human-modified landscapes, habitat selection might not always lead to increased fitness because animals undervalue high-quality resources that appear less attractive than those of lower quality. In the American West, agriculture has modified landscapes, yet little is known about whether agricultural changes alter the reliability of the cues animals use to identify habitat quality; ultimately forming maladaptive breeding strategies where behavioral cues are mismatched with survival outcomes. Using the greater sage-grouse, a species highly dependent upon sagebrush landscapes, we (1) evaluated how females select nesting habitats based on sagebrush type, along with livestock grazing related linear and point features, and other biotic, abiotic characteristics, given hypothesized influences on hiding cover, microclimate and predator travel routes and perches, (2) compared habitat selection information with results for nest survival estimates to evaluate if selection appears to be adaptive or not, and (3) used our results to evaluate the most appropriate strategies for this species in a grazing-modified landscape. Nest-site selection for sagebrush type appears to be maladaptive: in the most-preferred sagebrush type, nest survival rate was one-fourth the rate realized by females nesting in the sagebrush type avoided. Nest survival was four times higher for nests placed away from (>100 m), rather than next to (1 m), the nearest fence, and survival was lower within sites with higher cow pie density (a proxy for previous grazing intensity). Live and dead grasses influenced selection and survival in opposing ways such that dead grass was selected for but resulted in reduced survival while live grass was avoided but resulted in increased survival. Results collectively provide the first empirical evidence that a specific type of sagebrush acts as an ecological trap while another sagebrush type is undervalued. These results also suggest that adding more fences to control livestock grazing systems will likely reduce sage-grouse nest survival.

Variation in the vital rates of an Antarctic marine predator: the role of individual heterogeneity.

Variation in life-history traits such as lifespan and lifetime reproductive output is thought to arise, in part, due to among-individual differences in the underlying probabilities of survival and reproduction. However, the stochastic nature of demographic processes can also generate considerable variation in fitness-related traits among otherwise-identical individuals. An improved understanding of life-history evolution and population dynamics therefore depends on evaluating the relative role of each of these processes. Here, we used a 33-yr data set with reproductive histories for 1,274 female Weddell seals from Erebus Bay, Antarctica, to assess the strength of evidence for among-individual heterogeneity in the probabilities of survival and reproduction, while accounting for multiple other sources of variation in vital rates. Our analysis used recent advances in Bayesian model selection techniques and diagnostics to directly compare model fit and predictive power between models that included individual effects on survival and reproduction to those that did not. We found strong evidence for costs of reproduction to both survival and future reproduction, with breeders having rates of survival and subsequent reproduction that were 3% and 6% lower than rates for non-breeders. We detected age-related changes in the rates of survival and reproduction, but the patterns differed for the two rates. Survival rates steadily declined from 0.92 at age 7 to 0.56 at the maximal age of 31 yr. In contrast, reproductive rates increased from 0.68 at age 7 to 0.79 at age 16 and then steadily declined to 0.37 for the oldest females. Models that included individual effects explained more variation in observed life histories and had better estimated predictive power than those that did not, indicating their importance in understanding sources of variation among individuals in life-history traits. We found that among-individual heterogeneity in survival was small relative to that for reproduction. Our study, which found patterns of variation in vital rates that are consistent with a series of predictions from life-history theory, is the first to provide a thorough assessment of variation in important vital rates for a long-lived, high-latitude marine mammal while taking full advantage of recent developments in model evaluation.

Patterns of age-related change in reproductive effort differ in the pre-natal and post-natal periods in a long-lived mammal.

Age-related changes in maternal reproductive allocation for long-lived species are a key prediction from life-history theory. Theoretical and empirical work suggests that allocation may increase with age due to constraint (increases with experience) or restraint (increases with age in the face of declining residual reproductive value), and may decrease among the oldest aged animals due to senescence in reproductive function. Here, we use a hierarchical modelling approach to investigate the age-related patterns of change in maternal reproductive effort in the Weddell seal, a long-lived marine mammal with a protracted period of maternal care during which mothers allocate a large proportion of body mass while feeding little. We find that maternal allocation increases with age for young mothers during both the pre-natal and post-natal periods. In contrast, older mothers demonstrate a senescent decline in pre-natal allocation but allocate more of their declining resources to their offspring during the post-natal period. We also find strong evidence for the importance of individual effects in reproductive allocation among mothers: some mothers consistently produce heavier (or lighter) pups than expected. Our results indicate that maternal allocation changes over a mother's reproductive life span and that age-specific differences differ in notable ways in pre-natal and post-natal periods.

Tight coupling of primary production and marine mammal reproduction in the Southern Ocean.

Polynyas are areas of open water surrounded by sea ice and are important sources of primary production in high-latitude marine ecosystems. The magnitude of annual primary production in polynyas is controlled by the amount of exposure to solar radiation and sensitivity to changes in sea-ice extent. The degree of coupling between primary production and production by upper trophic-level consumers in these environments is not well understood, which prevents reliable predictions about population trajectories for species at higher trophic levels under potential future climate scenarios. In this study, we find a strong, positive relationship between annual primary production in an Antarctic polynya and pup production by ice-dependent Weddell seals. The timing of the relationship suggests reproductive effort increases to take advantage of high primary production occurring in the months after the birth pulse. Though the proximate causal mechanism is unknown, our results indicate tight coupling between organisms at disparate trophic levels on a short timescale, deepen our understanding of marine ecosystem processes, and raise interesting questions about why such coupling exists and what implications it has for understanding high-latitude ecosystems.

Female Weddell seals show flexible strategies of colony attendance related to varying environmental conditions.

Many animal life cycles involve movements among different habitats to fulfill varying resource demands. There are inherent costs associated with such movements, and the decision to leave or stay at a given location ought to be motivated by the benefits associated with potential target habitats. Because movement patterns, especially those associated with reproduction, can have important implications for the success (survival, reproduction) of individual animals, and therefore a population's dynamics, it is important to identify and understand their sources of variation (environmental and individual). Here, using a mark-recapture, multistate modeling approach, we investigated a set of a priori hypotheses regarding sources and patterns of variation in breeding-colony attendance for Weddell seal (Leptonychotes weddellii) females on sabbatical from pup production. For such females, colony attendance might be motivated by predation avoidance and positive social interactions related to reproduction, but some costs, such as reduced foraging opportunities or aggressive interactions with conspecifics, might also exist. We expected these benefits and costs to vary with a female's condition and the environment. Results revealed that the probability of being absent from colonies was higher (1) in years when the extent of local sea ice was larger, (2) for the youngest and oldest individuals, and (3) for females with less reproductive experience. We also found substantial levels of residual individual heterogeneity in these rates. Based on our a priori predictions, we postulate that the decision to attend breeding colonies or not is directly influenced by an individual's physiological condition, as well as by the ice-covered distance to good foraging areas, availability of predator-free haul-out sites, and the level of negative interactions with conspecifics inside colonies. Our results support the idea that in iteroparous species, and colonial animals in particular, seasonal and temporary movements from/to reproductive sites represent flexible behavioral strategies that can play an important role in coping with environmental variability.

An evolutionary perspective on reproductive individual heterogeneity in a marine vertebrate.

Although the quantification of individual heterogeneity in wild populations' vital rates has recently attracted growing interest among ecologists, the investigation of its evolutionary consequences remains limited, mainly because of the difficulties in assessing fitness and heritability from field studies on free-ranging animals. In the presence of individual variability, evaluation of fitness consequences can notably be complicated by the existence of trade-offs among different vital rates. In this study, to further assess the evolutionary significance of previously quantified levels of individual heterogeneity in female Weddell seal (Leptonychotes weddellii Lesson) reproductive rates (Chambert et al. ), we investigated how several life-history characteristics of female offspring were related to their mother's reproductive rate, as well as to other maternal traits (age and experience) and environmental conditions at birth. The probability and age of first reproduction (recruitment) of female offspring was not related to their mother's reproductive rate, suggesting the absence of a maternal trade-off between the number and quality of offspring a female produces. Evidence of a positive, but relatively weak, relationship between the reproductive rates of a mother and her female offspring was found, suggesting some degree of heritability in this trait. Using a simulation approach based on these statistical findings, we showed that substantial differences in the number of grandchildren, produced through female progeny, can be expected among females with different reproductive rates. Despite the presence of substantial stochastic variability, due to environmental fluctuations and other unidentified mechanisms, and in the light of the fact that the metrics obtained do not provide a full measure of real fitness, our results do suggest that the individual reproductive variability found in female Weddell seals could potentially have important fitness consequences.

Time-lagged variation in pond density and primary productivity affects duck nest survival in the Prairie Pothole Region.

The Prairie Pothole Region (PPR) is the primary breeding region for most species of North American dabbling ducks (Anas spp.). Conservation of these species is guided in part by knowledge of relationships between nest survival probability and habitat features. Positive relationships between duck nest survival and amount and configuration of herbaceous perennial vegetation have been observed in previous studies, but these 2- to 4-year studies might not have adequately characterized the temporal effect of wet-dry episodes on nest survival. Over an eight-year period, we studied nest survival of five species of ducks in the PPR relative to spatial and temporal variation in pond density, primary productivity, and hydrologic status of wetlands, soil, and vegetation on 52 study sites selected to span a gradient of spatial variation in proportion of herbaceous perennial vegetation and in number of wetland basins. We observed the fate of 12 754 nests. Consistent with past studies, 90% of nests that failed to hatch were destroyed by predators. Nest survival probability was positively related to current-year pond density and primary productivity, negatively related to pond density and primary productivity during the previous two years, and positively related to the number of wetland basins on the study site. Predicted relationships between nest survival and proportion or configuration of herbaceous perennial vegetation in the surrounding landscape were not supported. For mallard (Anas platyrhynchos), median estimated nest survival probability ranged from 0.02 (SE = 0.01) to 0.22 (SE = 0.02). Estimated nest survival was greatest on sites with numerous wetland basins that had transitioned from dry, unproductive conditions to wet, productive conditions in the previous 1-2 years. Our results were consistent with time-lagged responses of food webs to resource pulses in a broad array of ecosystems. Our study highlighted the importance of wetland basins and wet-dry episodes to duck nest survival in the PPR. Current habitat conservation efforts focus on landscapes with numerous wetland basins and a high proportion of herbaceous perennial vegetation. Our results suggest that future conservation efforts should focus on preserving high-density wetland complexes across as large a geographic extent as possible even in cropland-dominated landscapes.

Variability in temporary emigration rates of individually marked female Weddell seals prior to first reproduction.

In many species, temporary emigration (TE) is a phenomenon, often indicative of life-history characteristics such as dormancy, skipped reproduction, or partial migration, whereby certain individuals in a population are temporarily unobservable at a particular site. TE may be a flexible condition-dependent strategy that allows individuals to mitigate effects of adverse conditions. Consequently, TE rates ought to be highly variable, but sources of variations are poorly understood for most species. We used data from known-aged female Weddell seals (Leptonychotes weddellii) tagged in Erebus Bay, Antarctica, to investigate sources of variation in TE rates prior to reproduction and to evaluate possible implications for age-specific probability of first reproduction. TE rates were near 1 the year after birth, decreased to an average of 0.15 (SE = 0.01) by age 8, and were similar thereafter. TE rates varied substantially from year-to-year and were lower for seals that attended reproductive colonies the previous year than for seals that did not attend (e.g., ψ(i,age 8)(UU) - ψ(i, age 8)(PU) = 0.22). Recruitment rates were marginally greater for seals that did attend than for seals that did not attend colonies the previous year. For Weddell seals specifically, our results suggest that (1) motivation to attend colonies varied temporally, (2) as seals grew older they had increased motivation to attend even before reproductive maturity, and (3) seals appear to follow various attendance strategies. More broadly, our results support the idea of TE as a variable, condition-dependent strategy, and highlight the utility of TE models for providing population and life-history insights for diverse taxa.

Environmental extremes versus ecological extremes: impact of a massive iceberg on the population dynamics of a high-level Antarctic marine predator.

Extreme events have been suggested to play a disproportionate role in shaping ecological processes, but our understanding of the types of environmental conditions that elicit extreme consequences in natural ecosystems is limited. Here, we investigated the impact of a massive iceberg on the dynamics of a population of Weddell seals. Reproductive rates of females were reduced, but survival appeared unaffected. We also found suggestive evidence for a prolonged shift towards higher variability in reproductive rates. The annual number of females attending colonies showed unusual swings during the iceberg period, a pattern that was apparently the consequence of changes in sea-ice conditions. In contrast to the dramatic effects that were recorded in nearby populations of emperor penguins, our results suggest that this unusual environmental event did not have an extreme impact on the population of seals in the short-term, as they managed to avoid survival costs and were able to rapidly re-achieve high levels of reproduction by the end of the perturbation. Nevertheless, population projections suggest that even this modest impact on reproductive rates could negatively affect the population in the long run if such events were to occur more frequently, as is predicted by models of climate change.

Evaluating the demographic buffering hypothesis with vital rates estimated for Weddell seals from 30 years of mark-recapture data.

1. Life-history theory predicts that those vital rates that make larger contributions to population growth rate ought to be more strongly buffered against environmental variability than are those that are less important. Despite the importance of the theory for predicting demographic responses to changes in the environment, it is not yet known how pervasive demographic buffering is in animal populations because the validity of most existing studies has been called into question because of methodological deficiencies. 2. We tested for demographic buffering in the southern-most breeding mammal population in the world using data collected from 5558 known-age female Weddell seals over 30 years. We first estimated all vital rates simultaneously with mark-recapture analysis and then estimated process variance and covariance in those rates using a hierarchical Bayesian approach. We next calculated the population growth rate's sensitivity to changes in each of the vital rates and tested for evidence of demographic buffering by comparing properly scaled values of sensitivity and process variance in vital rates. 3. We found evidence of positive process covariance between vital rates, which indicates that all vital rates are affected in the same direction by changes in annual environment. Despite the positive correlations, we found strong evidence that demographic buffering occurred through reductions in variation in the vital rates to which population growth rate was most sensitive. Process variation in vital rates was inversely related to sensitivity measures such that variation was greatest in breeding probabilities, intermediate for survival rates of young animals and lowest for survival rates of older animals. 4. Our work contributes to a small but growing set of studies that have used rigorous methods on long-term, detailed data to investigate demographic responses to environmental variation. The information from these studies improves our understanding of life-history evolution in stochastic environments and provides useful information for predicting population responses to future environmental change. Our results for an Antarctic apex predator also provide useful baselines from a marine ecosystem when its top- and middle-trophic levels were not substantially impacted by human activity.

Challenges and opportunities for comparative studies of survival rates: An example with male pinnipeds.

Survival rates are a central component of life-history strategies of large vertebrate species. However, comparative studies seldom investigate interspecific variation in survival rates with respect to other life-history traits, especially for males. The lack of such studies could be due to the challenges associated with obtaining reliable datasets, incorporating information on the 0-1 probability scale, or dealing with several types of measurement error in life-history traits, which can be a computationally intensive process that is often absent in comparative studies. We present a quantitative approach using a Bayesian phylogenetically controlled regression with the flexibility to incorporate uncertainty in estimated survival rates and quantitative life-history traits while considering genetic similarity among species and uncertainty in relatedness. As with any comparative analysis, our approach makes several assumptions regarding the generalizability and comparability of empirical data from separate studies. Our model is versatile in that it can be applied to any species group of interest and include any life-history traits as covariates. We used an unbiased simulation framework to provide "proof of concept" for our model and applied a slightly richer model to a real data example for pinnipeds. Pinnipeds are an excellent taxonomic group for comparative analysis, but survival rate data are scarce. Our work elucidates the challenges associated with addressing important questions related to broader ecological life-history patterns and how survival-reproduction trade-offs might shape evolutionary histories of extant taxa. Specifically, we underscore the importance of having high-quality estimates of age-specific survival rates and information on other life-history traits that reasonably characterize a species for accurately comparing across species.

Nutrient sources for offspring formation: diet-mother and mother-offspring isotope discrimination in domesticated gallinaceous birds.

Stable isotope techniques can be used to assess nutrient acquisition and allocation strategies used to produce offspring. Before stable isotope techniques can be employed, researchers need reliable isotope discrimination values. In this context, isotope discrimination compares the difference in the isotope ratio between the maternal-offspring tissue that occurs during nutrient transfer prior to egg laying. Currently, isotope discrimination values are unknown between the maternal blood constituents - that reflect different temporal scales of integration - and downy feathers of their offspring. In this study, we experimentally derive isotope discrimination relationships between maternal diet-blood constituents for egg laying, and between maternal blood constituents-down feathers of offspring in an experiment with 3 types of domesticated gallinaceous birds raised on known diets. Our experiment is the first to report isotope discrimination values for maternal blood constituents-down of offspring in avian taxa and provides a new sampling technique that is less invasive than previously available as collecting down does not require sampling viable eggs or individuals. Future researchers can use these results to assist in identifying nutrient sources used by adult birds to produce young.

An evaluation of density-dependent and density-independent influences on population growth rates in Weddell seals.

Much of the existing literature that evaluates the roles of density-dependent and density-independent factors on population dynamics has been called into question in recent years because measurement errors were not properly dealt with in analyses. Using state-space models to account for measurement errors, we evaluated a set of competing models for a 22-year time series of mark-resight estimates of abundance for a breeding population of female Weddell seals (Leptonychotes weddellii) studied in Erebus Bay, Antarctica. We tested for evidence of direct density dependence in growth rates and evaluated whether equilibrium population size was related to seasonal sea-ice extent and the Southern Oscillation Index (SOI). We found strong evidence of negative density dependence in annual growth rates for a population whose estimated size ranged from 438 to 623 females during the study. Based on Bayes factors, a density-dependence-only model was favored over models that also included environmental covariates. According to the favored model, the population had a stationary distribution with a mean of 497 females (SD = 60.5), an expected growth rate of 1.10 (95% credible interval = 1.08-1.15) when population size was 441 females, and a rate of 0.90 (95% credible interval = 0.87-.93) for a population of 553 females. A model including effects of SOI did receive some support and indicated a positive relationship between SOI and population size. However, effects of SOI were not large, and including the effect did not greatly reduce our estimate of process variation. We speculate that direct density dependence occurred because rates of adult survival, breeding, and temporary emigration were affected by limitations on per capita food resources and space for parturition and pup-rearing. To improve understanding of the relative roles of various demographic components and their associated vital rates to population growth rate, mark-recapture methods can be applied that incorporate both environmental covariates and the seal abundance estimates that were developed here. An improved understanding of why vital rates change with changing population abundance will only come as we develop a better understanding of the processes affecting marine food resources in the Southern Ocean.

Evaluating sample size to estimate genetic management metrics in the genomics era.

Inbreeding and relationship metrics among and within populations are useful measures for genetic management of wild populations, but accuracy and precision of estimates can be influenced by the number of individual genotypes analysed. Biologists are confronted with varied advice regarding the sample size necessary for reliable estimates when using genomic tools. We developed a simulation framework to identify the optimal sample size for three widely used metrics to enable quantification of expected variance and relative bias of estimates and a comparison of results among populations. We applied this approach to analyse empirical genomic data for 30 individuals from each of four different free-ranging Rocky Mountain bighorn sheep (Ovis canadensis canadensis) populations in Montana and Wyoming, USA, through cross-species application of an Ovine array and analysis of approximately 14,000 single nucleotide polymorphisms (SNPs) after filtering. We examined intra- and interpopulation relationships using kinship and identity by state metrics, as well as FST between populations. By evaluating our simulation results, we concluded that a sample size of 25 was adequate for assessing these metrics using the Ovine array to genotype Rocky Mountain bighorn sheep herds. However, we conclude that a universal sample size rule may not be able to sufficiently address the complexities that impact genomic kinship and inbreeding estimates. Thus, we recommend that a pilot study and sample size simulation using R code we developed that includes empirical genotypes from a subset of populations of interest would be an effective approach to ensure rigour in estimating genomic kinship and population differentiation.

Respiratory pathogens and their association with population performance in Montana and Wyoming bighorn sheep populations.

Respiratory disease caused by Mycoplasma ovipneumoniae and Pasteurellaceae poses a formidable challenge for bighorn sheep (Ovis canadensis) conservation. All-age epizootics can cause 10-90% mortality and are typically followed by multiple years of enzootic disease in lambs that hinders post-epizootic recovery of populations. The relative frequencies at which these epizootics are caused by the introduction of novel pathogens or expression of historic pathogens that have become resident in the populations is unknown. Our primary objectives were to determine how commonly the pathogens associated with respiratory disease are hosted by bighorn sheep populations and assess demographic characteristics of populations with respect to the presence of different pathogens. We sampled 22 bighorn sheep populations across Montana and Wyoming, USA for Mycoplasma ovipneumoniae and Pasteurellaceae and used data from management agencies to characterize the disease history and demographics of these populations. We tested for associations between lamb:ewe ratios and the presence of different respiratory pathogen species. All study populations hosted Pasteurellaceae and 17 (77%) hosted Mycoplasma ovipneumoniae. Average lamb:ewe ratios for individual populations where both Mycoplasma ovipneumoniae and Pasteurellaceae were detected ranged from 0.14 to 0.40. However, average lamb:ewe ratios were higher in populations where Mycoplasma ovipneumoniae was not detected (0.37, 95% CI: 0.27-0.51) than in populations where it was detected (0.25, 95% CI: 0.21-0.30). These findings suggest that respiratory pathogens are commonly hosted by bighorn sheep populations and often reduce recruitment rates; however ecological factors may interact with the pathogens to determine population-level effects. Elucidation of such factors could provide insights for management approaches that alleviate the effects of respiratory pathogens in bighorn sheep. Nevertheless, minimizing the introduction of novel pathogens from domestic sheep and goats remains imperative to bighorn sheep conservation.

Assessing respiratory pathogen communities in bighorn sheep populations: Sampling realities, challenges, and improvements.

Respiratory disease has been a persistent problem for the recovery of bighorn sheep (Ovis canadensis), but has uncertain etiology. The disease has been attributed to several bacterial pathogens including Mycoplasma ovipneumoniae and Pasteurellaceae pathogens belonging to the Mannheimia, Bibersteinia, and Pasteurella genera. We estimated detection probability for these pathogens using protocols with diagnostic tests offered by a fee-for-service laboratory and not offered by a fee-for-service laboratory. We conducted 2861 diagnostic tests on swab samples collected from 476 bighorn sheep captured across Montana and Wyoming to gain inferences regarding detection probability, pathogen prevalence, and the power of different sampling methodologies to detect pathogens in bighorn sheep populations. Estimated detection probability using fee-for-service protocols was less than 0.50 for all Pasteurellaceae and 0.73 for Mycoplasma ovipneumoniae. Non-fee-for-service Pasteurellaceae protocols had higher detection probabilities, but no single protocol increased detection probability of all Pasteurellaceae pathogens to greater than 0.50. At least one protocol resulted in an estimated detection probability of 0.80 for each pathogen except Mannheimia haemolytica, for which the highest detection probability was 0.45. In general, the power to detect Pasteurellaceae pathogens at low prevalence in populations was low unless many animals were sampled or replicate samples were collected per animal. Imperfect detection also resulted in low precision when estimating prevalence for any pathogen. Low and variable detection probabilities for respiratory pathogens using live-sampling protocols may lead to inaccurate conclusions regarding pathogen community dynamics and causes of bighorn sheep respiratory disease epizootics. We recommend that agencies collect multiples samples per animal for Pasteurellaceae detection, and one sample for Mycoplasma ovipneumoniae detection from at least 30 individuals to reliably detect both Pasteurellaceae and Mycoplasma ovipneumoniae at the population-level. Availability of PCR diagnostic tests to wildlife management agencies would improve the ability to reliably detect Pasteurellaceae in bighorn sheep populations.