BACKGROUND
innate and adaptive immune responses evolved in vertebrates as a first and secondary line of defense, respectively, against a diverse and changing array of pathogenic organisms.
the effectiveness of these immunologic responses, and hence the fitness of individuals, populations, and species, is driven by pathogen exposure history and the immunogenetic repertoire of major histocompatibility complex  genes and non-mhc genes  <cit> .
novel, highly virulent pathogens can overwhelm host immune responses not primed to their exposure, and such pathogens can be a strong selective force, reducing the distribution and abundance of a species over short timeframes  through effects on survival and reproductive success  <cit> .
over multiple generations, a history of ongoing pathogen exposure theoretically should select for more resistant immunogenotypes that limit fitness impacts by responding effectively upon initial exposure  or re-exposure .
bighorn sheep  are a useful model for examining this interplay between disease, demography, and genetics.
they are a polygynous, highly philopatric species found in small, fragmented populations in the mountainous regions of western north america  <cit> .
they are highly susceptible to infectious disease, and outbreaks of disease regularly cause high morbidity and mortality  <cit> .
the history of population die-offs dates back to european settlement of the western united states > <dig> yrs ago  <cit> , indicating that novel pathogens were likely introduced by contact with domestic sheep .
pneumonia epizootics appear to be driven by density-dependence, serving to constrain population size  <cit> , and presumably selecting for the most fit genotypes.
however, small populations of bighorn sheep also are prone to inbreeding and genetic drift, making it difficult to understand the relative importance of pathogen-mediated selection, drift, and inbreeding on genetic variability and fitness.
we approached this problem by testing whether genetic background and previous pathogen exposure influenced survival when animals from two different founder populations were simultaneously translocated into the san andres mountains  in new mexico, usa.
the sam once supported the largest population of native bighorn sheep in the state.
however, by the late 1990's, a combination of disease, mountain lion  predation, and drought had reduced this population to the point of extinction, and translocation from captive and/or free-ranging herds was necessary to reestablish a self-sustaining population in the sam.
the two founder populations chosen for reintroduction were a genetically diverse free-ranging herd in the kofa national wildlife refuge , arizona, and a less diverse captive herd in the red rock wildlife area  that was originally derived from the native sam population.
in november  <dig>   <dig> bighorn sheep were translocated into the sam from the knwr  and the rrwa , and  <dig> more bighorn were translocated from knwr in november  <dig>  we examined genetic variation at  <dig> microsatellite loci to compare genetic variability, and we conducted population health analyses at the time of capture to assess infectious disease exposure.
radiocollared sheep released into the sam were then monitored through february  <dig> to determine survivorship and cause specific mortality.
methods
study animals
the sam mountain range represents the largest amount of high quality bighorn sheep habitat in new mexico, and the native sam bighorn sheep population exceeded  <dig> animals prior to a psoroptic scabies epizootic that began in the late 1970s.
the population remained around  <dig> animals through the early 1990's and then declined to a single ewe that was temporarily brought into captivity in  <dig> for less than two weeks.
nine animals sampled just prior to extirpation of the population had low mean heterozygosity at mhc  and microsatellite  loci  <cit> .
immediately prior to the reintroduction effort that began in november  <dig>  the sam population consisted of  <dig> rams that had been translocated from the rrwa as part of a sentinel disease study, and the single native ewe that was released along with the sentinel rams  <cit> .
a total of  <dig> bighorn sheep were captured, sampled, and translocated to the sam in  <dig> and  <dig> from rrwa and knwr.
in  <dig>  nine offspring from these animals were captured and sampled for genetics and disease surveillance in the sam.
the survival analysis presented in this paper focuses only on animals translocated to the sam, and data from these  <dig> offspring born in the sam are presented for descriptive purposes only.
the rrwa is a 500-hectare enclosure for bighorn sheep located in central new mexico.
the rrwa population was initiated in  <dig> with founders from the sam, and has served as the source for > <dig> bighorn sheep translocated within new mexico since  <dig>  <cit> .
this population has been managed as a closed herd , and has low genetic diversity  relative to other free-ranging populations in the desert southwest   <cit> .
genetic and disease samples for this study were collected when bighorn sheep were captured for translocation to the sam in  <dig> .
the size of the rrwa population was < <dig> animals at the time of capture and translocation.
the knwr is located in southwestern arizona and contains a free-ranging population of bighorn sheep that has significantly higher heterozygosity  than the rrwa captive population  <cit> .
like the rrwa population, the knwr herd has been the source for transplants of hundreds of animals within arizona and to other western states.
genetic and disease samples for this study were collected at capture from animals translocated to the sam in  <dig>  and  <dig> .
the population at knwr has averaged about 700- <dig> animals since  <dig>  but dropped to < <dig> animals by  <dig>  since  <dig>  > <dig> bighorn sheep have been translocated from knwr.
genetic analysis
genomic dna was extracted from blood and genetic variability was assessed by examining  <dig> microsatellite loci , drb <dig> , bmc <dig>  celjp <dig>  bm <dig>  oarcp <dig>  tgla <dig>  fcb <dig>  irbp, celb <dig>  bm <dig>  celjp <dig>  bm <dig>  csrd <dig>  hsc, inra <dig>  inra <dig>  inra <dig>  maf <dig>  mcm <dig>  oarae <dig>  oarcp <dig>  oarfcb <dig>  sps113).
many of these loci have been examined in previous publications  <cit> , and we have used them to examine heterozygosity in different populations of bighorn sheep across their range in north america .
briefly, using fluorescently labeled microsatellite primers, microsatellites present in genomic dna were amplified by pcr .
pcr products were separated by capillary electrophoresis using an applied biosystems  <dig> dna analyzer , with a fluorescent-labeled base pair size standard in each lane.
image analysis and fragment size determination were carried out using strand software  <cit> .
deviations from hardy-weinberg equilibrium  and allele frequencies were examined with genepop  <cit> .
multi-locus heterozygosity  values were calculated for each animal as the preferred metric of genetic variability.
chapman et al.
<cit>  showed that different genetic metrics  are highly correlated and non-independent, and they advocated the use of the simplest metric, mlh, in future studies of heterozygote-fitness correlations .
for our multivariate analyses, mlh values were categorized into two categories using a cutoff >  <dig>  to indicate a generally low or high level of genetic diversity.
assessment of health and disease exposure status
antibodies are a marker of previous exposure and priming of the adaptive immune response to specific pathogens.
sera from blood collected at the time of capture were tested for the presence of antibodies to three viruses known to cause respiratory disease and pneumonia in bighorn sheep: bovine respiratory syncytial virus , bluetongue virus , and parainfluenza- <dig> virus .
assays were performed at the california animal health and food safety laboratory in davis, california, and results were classified as positive or negative.
results for btv were reported as positive or negative; titers ≥ 1: <dig> were considered positive for brsv, and titers > 1: <dig> were considered positive for pi <dig>
post-release monitoring
a vhf radiocollar  with mortality sensor was placed on animals to be translocated at the time of capture to facilitate monitoring after their release in the sam.
animals were transported to the san andres national wildlife refuge in the sam and released within  <dig> hrs of capture.
resight of individual animals by radio-signal was conducted by state and federal wildlife biologists through february  <dig> to detect mortalities and determine specific causes of death.
mortalities were investigated immediately and a field investigation and necropsy typically performed within  <dig> hrs of death.
mortalities were classified as lion predation based on criteria of hayes et al.
<cit> , or as pneumonia based on gross lesions present in the lungs.
all mortalities that could not be classified as predation or pneumonia based on field examination were categorized as unknown, but "other" than predation and pneumonia.
the number of days each animal survived was calculated as the difference between the date of release and the date of death, or the date sheep were last observed.
for individuals lost to follow up, the end date was the date individuals were last observed, and for sheep that were confirmed alive at the end of the study period, this date was february  <dig>   <dig>
statistical analysis
we used independent univariate and multivariate approaches to explore correlation within our dataset, and to identify factors most influencing pathogen exposure, number of days until death  and death due to suspected pneumonia.
independent variables evaluated for their impact on these three outcomes included sex, age-class , source population , year released , exposure status to each pathogen  at time of release, and genetic variability .
significant associations with outcomes variables, confounding, and effect modification were evaluated by stratified univariate analyses using exact categorical tests and the student's t test.
the relationship between mlh  and age and sex class were evaluated by the two-way anova in order to adjust for source population.
independent categorical variables were also examined for their relationship to loss to follow-up by the two-sided fisher exact test to determine whether variables were related to censorship, and therefore could not be included in the survival analyses.
binary logistic regression was used to evaluate whether independent variables were related to death due to pneumonia, which was the most common cause of death detected during the study period.
only sheep that died of suspected pneumonia or were confirmed alive at the end of the study were included in this analysis .
sheep that died from predation or unknown causes were excluded from this analysis.
variables were selected by backward stepwise elimination  and confidence intervals for the logistic model evaluating risk factors for pneumonia were calculated using conditional exact inference due to low sample size.
hosmer-lemeshow goodness of fit test was used to estimate overall fit of the final logistic model.
for the survival analyses, median survival time was calculated for all risk factors of interest measured at time of translocation.
the kaplan-meier product-limit estimator was used to estimate the survivor function, and the log-rank test  <cit>  was used to test the equality of survivor function for each independent variable with groups considered significantly different if p <  <dig> .
variables significantly associated with survival in these univariate analyses were evaluated for their relationship to failure rate  using the semi-parametric cox proportional hazards regression.
sex was significantly associated with censorship so this variable was excluded from survival analyses.
variables were selected for the cox proportional hazards model by manual backward elimination using a selection criterion of p ≤  <dig>  for terms to stay in the model, and all categories of any significant variable were retained in the model.
the breslow approximation method was used to handle tied failure times  <cit> .the proportional hazards assumption based on schoenfeld residuals was evaluated to determine if the relative risk for each variable of interest was the same in time for the duration of the study.
all statistical analyses were performed using stata se  <dig>  software .
RESULTS
characteristics of source populations
bighorn sheep translocated from rrwa were predominantly male , while 72%  of the sheep translocated from knwr were female .
age class was not significantly associated with source population , but twice as many sheep in the 3- <dig> yr old range were captured in knwr  compared to rrwa .
all sheep relocated to sam in  <dig> were from knwr, but these newly introduced sheep were demographically similar and did not differ with regard to pathogen exposure or genetic diversity from sheep introduced from knwr in  <dig>
serologic evidence of previous pathogen exposure varied substantially by source population for all pathogens, except pi <dig> .
sheep from knwr were more likely to be exposed to brsv  compared to sheep from rrwa .
none of the sheep captured in knwr had evidence of exposure to btv, while 12/ <dig> sheep from rrwa were seropositive to btv .
results from the univariate analysis stratified by source population indicated that age was a significant confounder related to btv exposure for sheep from rrwa.
in fact, older sheep  captured at rrwa were  <dig> times more likely to be exposed to btv than younger age classes .
significant associations with age class were not detected for the other pathogens after adjusting for the effect of source population.
2002
2002
2005
2007
a bighorn sheep from rrwa and knwr were sampled at the time of their translocation into the sam , and offspring from these animals were sampled in  <dig> in the sam.
b btv = bluetongue virus, pi <dig> = parainfluenza- <dig> virus, brsv = bovine respiratory syncitial virus
c based on analysis of  <dig> microsatellite loci
three of the  <dig> microsatellite loci tested were monomorphic and were excluded from further analyses .
complete genotypes were determined for all animals at each of the remaining  <dig> loci with the exception of a single animal at a single locus.
allele frequencies did not significantly differ from hwe expectations across loci, and null alleles were not detected.
the genetic diversity parameter, mlh , was significantly higher in sheep captured in knwr  compared to sheep captured in rrwa .
mean mlh did not differ by sex or age class once comparisons were adjusted by source population.
specific causes of mortality
cause of death could be determined by field post-mortem examinations for  <dig> sheep.
pneumonia was suspected as a cause of mortality in  <dig> sheep, and mountain lion predation was identified as the cause of death in  <dig> sheep.
specific cause of death could not be determined for  <dig> sheep, but predation and pneumonia were excluded as their cause of death.
deaths due to lion predation occurred in  <dig>   <dig>  and  <dig>  while deaths due to pneumonia occurred every year .
most deaths due to pneumonia  occurred in the fall from september thru november, and seven of these occurred in a cluster in the fall of  <dig>  pneumonia impacted all age classes with 27%  of 1- <dig> year old sheep, 35%  of 3- <dig> year old sheep and 44%  of >  <dig> year old sheep dying of pneumonia.
mean mlh was not significantly different among the  <dig> sheep that survived  compared to the  <dig> sheep that died of pneumonia , but this analyses lacked power due to low sample size .
similarly, risk factors were not significantly associated with death due to pneumonia in the multivariate analyses using exact conditional logistic regression.
risk factors associated with survival
year of release was not associated with overall survival in univariate or multivariate analyses and survival time did not differ significantly by source population.
however, age class at time of release was marginally associated with survivorship  in the univariate analyses.
more than half of the sheep released at  <dig> to <  <dig> years old survived the entire study period , while median survival time was  <dig>  days for sheep  <dig> years to <  <dig> years old and only  <dig> days for sheep released when  <dig> years or older.
genetic diversity and pathogen exposure were not related to survival time, even after stratifying by source population.
risk factors found to be significantly associated with time to death in the cox proportional hazards multivariate framework were the oldest age class and previous exposure to brsv .
based on the hazards ratios of the cox proportional hazards model, we found that sheep with serological evidence of previous exposure to brsv had approximately one-third the risk of death compared to sheep that had not been exposed to brsv prior to release.
age was a significant confounding factor in the survival analyses, and as expected, older age class was associated with increasing risk and a shorter time to death.
sheep released at ages ≥  <dig> years had  <dig>  times the hazard rate of sheep released at  <dig> to  <dig> years of age.
source population, mlh as a continuous and binary variable, year of release, and previous exposure to other pathogens were not associated with survival once we accounted for the influence of age and exposure to brsv on time to death.
associations with sex could not be evaluated in the modeling procedure because males were more likely to be lost to follow up.
evaluation of the proportional hazards assumption for the final model based on a test of schoenfeld residuals indicated that the relative risk for each variable of interest was the same in time for the duration of the study .
a youngest age class  designated as reference category
b brsv = bovine respiratory syncitial virus
discussion
a key challenge for conservation biologists is to determine the most appropriate demographic and genetic management strategies for wildlife populations threatened by disease  <cit> .
bighorn sheep in north america provide a useful model for examining this issue because they are highly susceptible to infectious diseases, and they are frequently translocated to re-establish or augment populations.
in this study, we took advantage of a management effort - the restoration of bighorn sheep to the sam - to test whether genetic background and previous pathogen exposure influenced the survival of translocated bighorn sheep.
variability and local adaptation are considered to be key genetic factors influencing the persistence of small populations of bighorn sheep  <cit> , and this strongly influenced the choice of source animals for the sam translocation.
since the rrwa herd was founded with animals from the sam, their translocation presumably maximized the retention of locally adapted gene complexes.
however, the rrwa herd had low genetic variability  because it had been managed for > <dig> years as an inbred population.
in contrast, the knwr population was much more diverse , ensuring that the newly established sam population would be more variable than one established only with animals from the rrwa captive population.
while the genetic management goal of achieving increased heterozygosity in the new sam population was accomplished, it was not associated with enhanced survivorship of the individual sheep that were translocated.
bighorn sheep from the knwr had significantly higher genetic diversity than those from the rrwa, but neither source population nor genetic background  influenced their survival in the sam.
our failure to detect a significant relationship between genetic diversity and adult survival is consistent with the chapman et al.
<cit>  meta-analysis of hfc studies.
in their comprehensive analysis, chapman et al.
<cit>  concluded that the effects of such correlations are very weak, and that other proposed measures of genetic variation  are no more powerful than mlh for detecting relationships.
we acknowledge that a larger sample size and more marker loci would have increased the power of our analysis, but any such effect is apparently very weak.
our failure to detect a fitness effect with neutral markers  is perhaps not surprising given that gutierrez et al.
<cit>  failed to find a strong association between mhc variation and disease resistance in bighorn sheep.
previous examination of genetic variability among bighorn sheep suggests that microsatellite diversity has been influenced primarily by neutral factors, and mhc diversity by balancing selection  <cit> .
these studies have shown that genetic distances between populations are roughly proportional to geographic distance, and that most genetic variability is apportioned within, rather than between, populations.
these results are best explained by isolation by distance , and provide little genetic evidence for local adaptation.
our results are consistent with this view - the rrwa animals derived from the original sam population did not have increased survivorship.
conversely, inbreeding within the rrwa population could have reduced fitness by increasing the accumulation and expression of deleterious recessive alleles  <cit> .
we found no evidence of such an effect among translocated adult sheep .
although we did not detect fitness differences among translocated animals, heterosis or genetic rescue effects may occur in f <dig> and later generations as the two source populations interbreed  <cit> .
while we were able to sample nine offspring born in the sam for genetic and pathogen analyses , we were unable to monitor their relative fitness and they were not included in any of our survival analyses.
pneumonia was the most frequent cause of death, accounting for 50%  of the documented mortalities.
this mirrors what has been reported for many other populations, illustrating the ongoing importance of pneumonia as a cause of morbidity and mortality across the range of bighorn sheep  <cit> .
several bacteria  and viruses  have been implicated as primary or secondary pathogens in enzootic and epizootic pneumonia, and we associate these organisms with pneumonic disease because they can be detected in sick and dead animals.
however, we know far less about causal relationships.
the multifactorial etiology of pneumonia, coupled with the complex interplay of host, pathogen, and environment, has made it very difficult to identify a single, specific cause of pneumonia.
instead of focusing on identifying pathogens after death, we looked at how previous pathogen exposure influenced the future survival of healthy bighorn sheep.
the presence of antibodies was used as a marker of pathogen exposure, indicating that viral transmission and infection had occurred at some point in the past.
animals translocated into the sam showed evidence of previous exposure to btv  as well as brsv and pi <dig> .
all animals were healthy at the time of translocation, yet previous brsv exposure was associated with increased survival regardless of age class.
this suggests that antibodies to this virus provided some level of protection against pneumonia after animals were released into the sam.
age class at the time of release was significantly related to survival, and younger animals lived longer than those >  <dig> yrs of age .
this intuitive result illustrates the demographic advantage of selecting young animals for reintroductions.
singer et al.
<cit>  found that translocation success was directly related to the number of animals - larger translocations resulted in populations that were more likely to persist over time.
however, pneumonia epizootics appear to be driven by density-dependence, and thus disease may constrain population growth and size  <cit> .
lethal removal of mountain lions was conducted in the sam throughout our study period  <cit> , likely influencing the relative proportion of mortality caused by predation versus pneumonia.
in the absence of controlled studies, we have no way of inferring whether lion control increased overall survivorship.
regardless, there were no significant relationships between mean mlh and death due to lion predation or pneumonia.
CONCLUSIONS
we conclude that increased heterozygosity did not increase survival following translocation, nor did it reduce individual susceptibility to dying specifically from pneumonia.
instead, we found that previous pathogen exposure was more important than genetic heterozygosity as a marker for predicting survival of translocated animals.
in their review of >  <dig> published hfc effects, chapman et al.
<cit>  reported a publication bias with non-significant effects being under-reported in the literature.
we acknowledge that confounding factors, such as age, and small sample sizes likely limited our ability to detect significant relationships between genetic variability, disease and fitness.
these limitations often plague investigators who conduct observational field studies in natural settings.
we encourage others to test the theoretical advantages of increased genetic variability using a variety of genetic markers in long term field studies that also evaluate ecological factors, such as disease, that are likely to influence survival.
every wildlife translocation is an experiment, but many translocations lack an experimental design.
whenever possible, translocations should be designed and evaluated to test hypotheses that will improve our understanding of how pathogen exposure and genetic variability influence fitness.
authors' contributions
wmb and mew designed and implemented the study.
mctp and ckj conducted genetic and statistical analyses, respectively.
all authors have approved the final manuscript.

 innate and adaptive immune responses evolved in vertebrates as a first and secondary line of defense, respectively, against a diverse and changing array of pathogenic organisms. the effectiveness of these immunologic responses, and hence the fitness of individuals, populations, and species, is driven by pathogen exposure history and the immunogenetic repertoire of major histocompatibility complex genes and non-mhc genes. novel, highly virulent pathogens can overwhelm host immune responses not primed to their exposure, and such pathogens can be a strong selective force, reducing the distribution and abundance of a species over short timeframes through effects on survival and reproductive success. over multiple generations, a history of ongoing pathogen exposure theoretically should select for more resistant immunogenotypes that limit fitness impacts by responding effectively upon initial exposure or re-exposure. bighorn sheep are a useful model for examining this interplay between disease, demography, and genetics. they are a polygynous, highly philopatric species found in small, fragmented populations in the mountainous regions of western north america. they are highly susceptible to infectious disease, and outbreaks of disease regularly cause high morbidity and mortality. the history of population die-offs dates back to european settlement of the western united states > yrs ago, indicating that novel pathogens were likely introduced by contact with domestic sheep. pneumonia epizootics appear to be driven by density-dependence, serving to constrain population size, and presumably selecting for the most fit genotypes. however, small populations of bighorn sheep also are prone to inbreeding and genetic drift, making it difficult to understand the relative importance of pathogen-mediated selection, drift, and inbreeding on genetic variability and fitness. we approached this problem by testing whether genetic background and previous pathogen exposure influenced survival when animals from two different founder populations were simultaneously translocated into the san andres mountains in new mexico, usa. the sam once supported the largest population of native bighorn sheep in the state. however, by the late 1990's, a combination of disease, mountain lion predation, and drought had reduced this population to the point of extinction, and translocation from captive and/or free-ranging herds was necessary to reestablish a self-sustaining population in the sam. the two founder populations chosen for reintroduction were a genetically diverse free-ranging herd in the kofa national wildlife refuge, arizona, and a less diverse captive herd in the red rock wildlife area that was originally derived from the native sam population. in november bighorn sheep were translocated into the sam from the knwr and the rrwa, and more bighorn were translocated from knwr in november we examined genetic variation at microsatellite loci to compare genetic variability, and we conducted population health analyses at the time of capture to assess infectious disease exposure. radiocollared sheep released into the sam were then monitored through february to determine survivorship and cause specific mortality. methods study animals the sam mountain range represents the largest amount of high quality bighorn sheep habitat in new mexico, and the native sam bighorn sheep population exceeded animals prior to a psoroptic scabies epizootic that began in the late 1970s. the population remained around animals through the early 1990's and then declined to a single ewe that was temporarily brought into captivity in for less than two weeks. nine animals sampled just prior to extirpation of the population had low mean heterozygosity at mhc and microsatellite loci. immediately prior to the reintroduction effort that began in november the sam population consisted of rams that had been translocated from the rrwa as part of a sentinel disease study, and the single native ewe that was released along with the sentinel rams. a total of bighorn sheep were captured, sampled, and translocated to the sam in and from rrwa and knwr. in nine offspring from these animals were captured and sampled for genetics and disease surveillance in the sam. the survival analysis presented in this paper focuses only on animals translocated to the sam, and data from these offspring born in the sam are presented for descriptive purposes only. the rrwa is a 500-hectare enclosure for bighorn sheep located in central new mexico. the rrwa population was initiated in with founders from the sam, and has served as the source for > bighorn sheep translocated within new mexico since. this population has been managed as a closed herd, and has low genetic diversity relative to other free-ranging populations in the desert southwest. genetic and disease samples for this study were collected when bighorn sheep were captured for translocation to the sam in. the size of the rrwa population was < animals at the time of capture and translocation. the knwr is located in southwestern arizona and contains a free-ranging population of bighorn sheep that has significantly higher heterozygosity than the rrwa captive population. like the rrwa population, the knwr herd has been the source for transplants of hundreds of animals within arizona and to other western states. genetic and disease samples for this study were collected at capture from animals translocated to the sam in and. the population at knwr has averaged about 700- animals since but dropped to < animals by since > bighorn sheep have been translocated from knwr. genetic analysis genomic dna was extracted from blood and genetic variability was assessed by examining microsatellite loci, drb, bmc celjp bm oarcp tgla fcb irbp, celb bm celjp bm csrd hsc, inra inra inra maf mcm oarae oarcp oarfcb sps113). many of these loci have been examined in previous publications, and we have used them to examine heterozygosity in different populations of bighorn sheep across their range in north america. briefly, using fluorescently labeled microsatellite primers, microsatellites present in genomic dna were amplified by pcr. pcr products were separated by capillary electrophoresis using an applied biosystems dna analyzer, with a fluorescent-labeled base pair size standard in each lane. image analysis and fragment size determination were carried out using strand software. deviations from hardy-weinberg equilibrium and allele frequencies were examined with genepop. multi-locus heterozygosity values were calculated for each animal as the preferred metric of genetic variability. chapman et al. showed that different genetic metrics are highly correlated and non-independent, and they advocated the use of the simplest metric, mlh, in future studies of heterozygote-fitness correlations. for our multivariate analyses, mlh values were categorized into two categories using a cutoff > to indicate a generally low or high level of genetic diversity. assessment of health and disease exposure status antibodies are a marker of previous exposure and priming of the adaptive immune response to specific pathogens. sera from blood collected at the time of capture were tested for the presence of antibodies to three viruses known to cause respiratory disease and pneumonia in bighorn sheep: bovine respiratory syncytial virus, bluetongue virus, and parainfluenza- virus. assays were performed at the california animal health and food safety laboratory in davis, california, and results were classified as positive or negative. results for btv were reported as positive or negative; titers ≥ 1: were considered positive for brsv, and titers > 1: were considered positive for pi post-release monitoring a vhf radiocollar with mortality sensor was placed on animals to be translocated at the time of capture to facilitate monitoring after their release in the sam. animals were transported to the san andres national wildlife refuge in the sam and released within hrs of capture. resight of individual animals by radio-signal was conducted by state and federal wildlife biologists through february to detect mortalities and determine specific causes of death. mortalities were investigated immediately and a field investigation and necropsy typically performed within hrs of death. mortalities were classified as lion predation based on criteria of hayes et al., or as pneumonia based on gross lesions present in the lungs. all mortalities that could not be classified as predation or pneumonia based on field examination were categorized as unknown, but "other" than predation and pneumonia. the number of days each animal survived was calculated as the difference between the date of release and the date of death, or the date sheep were last observed. for individuals lost to follow up, the end date was the date individuals were last observed, and for sheep that were confirmed alive at the end of the study period, this date was february statistical analysis we used independent univariate and multivariate approaches to explore correlation within our dataset, and to identify factors most influencing pathogen exposure, number of days until death and death due to suspected pneumonia. independent variables evaluated for their impact on these three outcomes included sex, age-class, source population, year released, exposure status to each pathogen at time of release, and genetic variability. significant associations with outcomes variables, confounding, and effect modification were evaluated by stratified univariate analyses using exact categorical tests and the student's t test. the relationship between mlh and age and sex class were evaluated by the two-way anova in order to adjust for source population. independent categorical variables were also examined for their relationship to loss to follow-up by the two-sided fisher exact test to determine whether variables were related to censorship, and therefore could not be included in the survival analyses. binary logistic regression was used to evaluate whether independent variables were related to death due to pneumonia, which was the most common cause of death detected during the study period. only sheep that died of suspected pneumonia or were confirmed alive at the end of the study were included in this analysis. sheep that died from predation or unknown causes were excluded from this analysis. variables were selected by backward stepwise elimination and confidence intervals for the logistic model evaluating risk factors for pneumonia were calculated using conditional exact inference due to low sample size. hosmer-lemeshow goodness of fit test was used to estimate overall fit of the final logistic model. for the survival analyses, median survival time was calculated for all risk factors of interest measured at time of translocation. the kaplan-meier product-limit estimator was used to estimate the survivor function, and the log-rank test was used to test the equality of survivor function for each independent variable with groups considered significantly different if p <. variables significantly associated with survival in these univariate analyses were evaluated for their relationship to failure rate using the semi-parametric cox proportional hazards regression. sex was significantly associated with censorship so this variable was excluded from survival analyses. variables were selected for the cox proportional hazards model by manual backward elimination using a selection criterion of p ≤ for terms to stay in the model, and all categories of any significant variable were retained in the model. the breslow approximation method was used to handle tied failure times.the proportional hazards assumption based on schoenfeld residuals was evaluated to determine if the relative risk for each variable of interest was the same in time for the duration of the study. all statistical analyses were performed using stata se software.  characteristics of source populations bighorn sheep translocated from rrwa were predominantly male, while 72% of the sheep translocated from knwr were female. age class was not significantly associated with source population, but twice as many sheep in the 3- yr old range were captured in knwr compared to rrwa. all sheep relocated to sam in were from knwr, but these newly introduced sheep were demographically similar and did not differ with regard to pathogen exposure or genetic diversity from sheep introduced from knwr in serologic evidence of previous pathogen exposure varied substantially by source population for all pathogens, except pi. sheep from knwr were more likely to be exposed to brsv compared to sheep from rrwa. none of the sheep captured in knwr had evidence of exposure to btv, while 12/ sheep from rrwa were seropositive to btv. results from the univariate analysis stratified by source population indicated that age was a significant confounder related to btv exposure for sheep from rrwa. in fact, older sheep captured at rrwa were times more likely to be exposed to btv than younger age classes. significant associations with age class were not detected for the other pathogens after adjusting for the effect of source population. 2002 2002 2005 2007 a bighorn sheep from rrwa and knwr were sampled at the time of their translocation into the sam, and offspring from these animals were sampled in in the sam. b btv = bluetongue virus, pi = parainfluenza- virus, brsv = bovine respiratory syncitial virus c based on analysis of microsatellite loci three of the microsatellite loci tested were monomorphic and were excluded from further analyses. complete genotypes were determined for all animals at each of the remaining loci with the exception of a single animal at a single locus. allele frequencies did not significantly differ from hwe expectations across loci, and null alleles were not detected. the genetic diversity parameter, mlh, was significantly higher in sheep captured in knwr compared to sheep captured in rrwa. mean mlh did not differ by sex or age class once comparisons were adjusted by source population. specific causes of mortality cause of death could be determined by field post-mortem examinations for sheep. pneumonia was suspected as a cause of mortality in sheep, and mountain lion predation was identified as the cause of death in sheep. specific cause of death could not be determined for sheep, but predation and pneumonia were excluded as their cause of death. deaths due to lion predation occurred in and while deaths due to pneumonia occurred every year. most deaths due to pneumonia occurred in the fall from september thru november, and seven of these occurred in a cluster in the fall of pneumonia impacted all age classes with 27% of 1- year old sheep, 35% of 3- year old sheep and 44% of > year old sheep dying of pneumonia. mean mlh was not significantly different among the sheep that survived compared to the sheep that died of pneumonia, but this analyses lacked power due to low sample size. similarly, risk factors were not significantly associated with death due to pneumonia in the multivariate analyses using exact conditional logistic regression. risk factors associated with survival year of release was not associated with overall survival in univariate or multivariate analyses and survival time did not differ significantly by source population. however, age class at time of release was marginally associated with survivorship in the univariate analyses. more than half of the sheep released at to < years old survived the entire study period, while median survival time was days for sheep years to < years old and only days for sheep released when years or older. genetic diversity and pathogen exposure were not related to survival time, even after stratifying by source population. risk factors found to be significantly associated with time to death in the cox proportional hazards multivariate framework were the oldest age class and previous exposure to brsv. based on the hazards ratios of the cox proportional hazards model, we found that sheep with serological evidence of previous exposure to brsv had approximately one-third the risk of death compared to sheep that had not been exposed to brsv prior to release. age was a significant confounding factor in the survival analyses, and as expected, older age class was associated with increasing risk and a shorter time to death. sheep released at ages ≥ years had times the hazard rate of sheep released at to years of age. source population, mlh as a continuous and binary variable, year of release, and previous exposure to other pathogens were not associated with survival once we accounted for the influence of age and exposure to brsv on time to death. associations with sex could not be evaluated in the modeling procedure because males were more likely to be lost to follow up. evaluation of the proportional hazards assumption for the final model based on a test of schoenfeld residuals indicated that the relative risk for each variable of interest was the same in time for the duration of the study. a youngest age class designated as reference category b brsv = bovine respiratory syncitial virus discussion a key challenge for conservation biologists is to determine the most appropriate demographic and genetic management strategies for wildlife populations threatened by disease. bighorn sheep in north america provide a useful model for examining this issue because they are highly susceptible to infectious diseases, and they are frequently translocated to re-establish or augment populations. in this study, we took advantage of a management effort - the restoration of bighorn sheep to the sam - to test whether genetic background and previous pathogen exposure influenced the survival of translocated bighorn sheep. variability and local adaptation are considered to be key genetic factors influencing the persistence of small populations of bighorn sheep, and this strongly influenced the choice of source animals for the sam translocation. since the rrwa herd was founded with animals from the sam, their translocation presumably maximized the retention of locally adapted gene complexes. however, the rrwa herd had low genetic variability because it had been managed for > years as an inbred population. in contrast, the knwr population was much more diverse, ensuring that the newly established sam population would be more variable than one established only with animals from the rrwa captive population. while the genetic management goal of achieving increased heterozygosity in the new sam population was accomplished, it was not associated with enhanced survivorship of the individual sheep that were translocated. bighorn sheep from the knwr had significantly higher genetic diversity than those from the rrwa, but neither source population nor genetic background influenced their survival in the sam. our failure to detect a significant relationship between genetic diversity and adult survival is consistent with the chapman et al. meta-analysis of hfc studies. in their comprehensive analysis, chapman et al. concluded that the effects of such correlations are very weak, and that other proposed measures of genetic variation are no more powerful than mlh for detecting relationships. we acknowledge that a larger sample size and more marker loci would have increased the power of our analysis, but any such effect is apparently very weak. our failure to detect a fitness effect with neutral markers is perhaps not surprising given that gutierrez et al. failed to find a strong association between mhc variation and disease resistance in bighorn sheep. previous examination of genetic variability among bighorn sheep suggests that microsatellite diversity has been influenced primarily by neutral factors, and mhc diversity by balancing selection. these studies have shown that genetic distances between populations are roughly proportional to geographic distance, and that most genetic variability is apportioned within, rather than between, populations. these results are best explained by isolation by distance, and provide little genetic evidence for local adaptation. our results are consistent with this view - the rrwa animals derived from the original sam population did not have increased survivorship. conversely, inbreeding within the rrwa population could have reduced fitness by increasing the accumulation and expression of deleterious recessive alleles. we found no evidence of such an effect among translocated adult sheep. although we did not detect fitness differences among translocated animals, heterosis or genetic rescue effects may occur in f and later generations as the two source populations interbreed. while we were able to sample nine offspring born in the sam for genetic and pathogen analyses, we were unable to monitor their relative fitness and they were not included in any of our survival analyses. pneumonia was the most frequent cause of death, accounting for 50% of the documented mortalities. this mirrors what has been reported for many other populations, illustrating the ongoing importance of pneumonia as a cause of morbidity and mortality across the range of bighorn sheep. several bacteria and viruses have been implicated as primary or secondary pathogens in enzootic and epizootic pneumonia, and we associate these organisms with pneumonic disease because they can be detected in sick and dead animals. however, we know far less about causal relationships. the multifactorial etiology of pneumonia, coupled with the complex interplay of host, pathogen, and environment, has made it very difficult to identify a single, specific cause of pneumonia. instead of focusing on identifying pathogens after death, we looked at how previous pathogen exposure influenced the future survival of healthy bighorn sheep. the presence of antibodies was used as a marker of pathogen exposure, indicating that viral transmission and infection had occurred at some point in the past. animals translocated into the sam showed evidence of previous exposure to btv as well as brsv and pi. all animals were healthy at the time of translocation, yet previous brsv exposure was associated with increased survival regardless of age class. this suggests that antibodies to this virus provided some level of protection against pneumonia after animals were released into the sam. age class at the time of release was significantly related to survival, and younger animals lived longer than those > yrs of age. this intuitive result illustrates the demographic advantage of selecting young animals for reintroductions. singer et al. found that translocation success was directly related to the number of animals - larger translocations resulted in populations that were more likely to persist over time. however, pneumonia epizootics appear to be driven by density-dependence, and thus disease may constrain population growth and size. lethal removal of mountain lions was conducted in the sam throughout our study period, likely influencing the relative proportion of mortality caused by predation versus pneumonia. in the absence of controlled studies, we have no way of inferring whether lion control increased overall survivorship. regardless, there were no significant relationships between mean mlh and death due to lion predation or pneumonia.  we conclude that increased heterozygosity did not increase survival following translocation, nor did it reduce individual susceptibility to dying specifically from pneumonia. instead, we found that previous pathogen exposure was more important than genetic heterozygosity as a marker for predicting survival of translocated animals. in their review of > published hfc effects, chapman et al. reported a publication bias with non-significant effects being under-reported in the literature. we acknowledge that confounding factors, such as age, and small sample sizes likely limited our ability to detect significant relationships between genetic variability, disease and fitness. these limitations often plague investigators who conduct observational field studies in natural settings. we encourage others to test the theoretical advantages of increased genetic variability using a variety of genetic markers in long term field studies that also evaluate ecological factors, such as disease, that are likely to influence survival. every wildlife translocation is an experiment, but many translocations lack an experimental design. whenever possible, translocations should be designed and evaluated to test hypotheses that will improve our understanding of how pathogen exposure and genetic variability influence fitness. authors' contributions wmb and mew designed and implemented the study. mctp and ckj conducted genetic and statistical analyses, respectively. all authors have approved the final manuscript. 
we explored this topic by examining whether genetic background and previous pathogen exposure influenced survival of translocated animals when captive-bred and free-ranging bighorn sheep were used to re-establish a population that had been extirpated in the san andres mountains in new mexico, usa. while heterosis or genetic rescue effects may occur in f and later generations as the two source populations interbreed, we conclude that previous pathogen exposure was a more important marker than genetic heterozygosity for predicting survival of translocated animals. every wildlife translocation is an experiment, and whenever possible, translocations should be designed and evaluated to test hypotheses that will further improve our understanding of how pathogen exposure and genetic variability influence fitness. a key challenge for conservation biologists is to determine the most appropriate demographic and genetic management strategies for wildlife populations threatened by disease. although the free-ranging source population had significantly higher multi-locus heterozygosity at microsatellite loci than the captive bred animals, neither source population nor genetic background significantly influenced survival or cause of death. the presence of antibodies to a respiratory virus known to cause pneumonia was associated with increased survival, but there was no correlation between genetic heterozygosity and the presence of antibodies to this virus.
 a key challenge for conservation biologists is to determine the most appropriate demographic and genetic management strategies for wildlife populations threatened by disease. we explored this topic by examining whether genetic background and previous pathogen exposure influenced survival of translocated animals when captive-bred and free-ranging bighorn sheep were used to re-establish a population that had been extirpated in the san andres mountains in new mexico, usa.  although the free-ranging source population had significantly higher multi-locus heterozygosity at microsatellite loci than the captive bred animals, neither source population nor genetic background significantly influenced survival or cause of death. the presence of antibodies to a respiratory virus known to cause pneumonia was associated with increased survival, but there was no correlation between genetic heterozygosity and the presence of antibodies to this virus.  although genetic theory predicts otherwise, increased heterozygosity was not associated with increased fitness among translocated animals. while heterosis or genetic rescue effects may occur in f and later generations as the two source populations interbreed, we conclude that previous pathogen exposure was a more important marker than genetic heterozygosity for predicting survival of translocated animals. every wildlife translocation is an experiment, and whenever possible, translocations should be designed and evaluated to test hypotheses that will further improve our understanding of how pathogen exposure and genetic variability influence fitness. 