Wildlife translocation: the conservation implications of pathogen exposure and genetic heterozygosity
© Boyce et al; licensee BioMed Central Ltd. 2011
Received: 27 April 2010
Accepted: 1 February 2011
Published: 1 February 2011
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 (Ovis canadensis) 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 30 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 (survival) among translocated animals. While heterosis or genetic rescue effects may occur in F1 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.
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 (MHC) genes and non-MHC genes [1, 2]. 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 (1-2 generations) 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 (innate immunity) or re-exposure (adaptive immunity).
Bighorn sheep (Ovis canadensis) 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 [3, 4]. The history of population die-offs dates back to European settlement of the western United States >200 yrs ago , indicating that novel pathogens were likely introduced by contact with domestic sheep (Ovis aries). 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 (SAM) 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 (Puma concolor) 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 (KNWR), Arizona, and a less diverse captive herd in the Red Rock Wildlife Area (RRWA) that was originally derived from the native SAM population. In November 2002, 51 bighorn sheep were translocated into the SAM from the KNWR (n = 20) and the RRWA (n = 31), and 30 more bighorn were translocated from KNWR in November 2005. We examined genetic variation at 33 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 2007 to determine survivorship and cause specific mortality.
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 200 animals prior to a psoroptic scabies epizootic that began in the late 1970s. The population remained around 30 animals through the early 1990's and then declined to a single ewe that was temporarily brought into captivity in 1999 for less than two weeks. Nine animals sampled just prior to extirpation of the population had low mean heterozygosity at MHC (0.075) and microsatellite (0.359) loci . Immediately prior to the reintroduction effort that began in November 2002, the SAM population consisted of 4 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 81 bighorn sheep were captured, sampled, and translocated to the SAM in 2002 and 2005 from RRWA and KNWR. In 2007, 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 9 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 1972 with founders from the SAM, and has served as the source for >260 bighorn sheep translocated within New Mexico since 1979 . This population has been managed as a closed herd (no immigrants), and has low genetic diversity (0.36) relative to other free-ranging populations in the desert southwest (0.44-0.63) . Genetic and disease samples for this study were collected when bighorn sheep were captured for translocation to the SAM in 2002 (n = 31). The size of the RRWA population was <100 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 (0.60) 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 2002 (n = 21) and 2005 (n = 30). The population at KNWR has averaged about 700-800 animals since 1981, but dropped to <400 animals by 2006. Since 1979, >560 bighorn sheep have been translocated from KNWR.
Genomic DNA was extracted from blood and genetic variability was assessed by examining 33 microsatellite loci (OarFCB128, MAF209, OarFCB304, MAF33, MAF48, MAF65, OarFCB11, MAF36, OarFCB266, ETH152 (= D5S2), DRB3 (MHC-linked loci), BMC1009, CELJP23, BM203, OARCP026, TGLA94, FCB193, IRBP, CELB9, BM6506, CELJP15, BM4107, CSRD247, HSC, INRA023, INRA063, INRA105, MAF214, McM527, OarAE129, OarCP49, OarFCB20, SPS113). Many of these loci have been examined in previous publications [6, 9, 10], and we have used them to examine heterozygosity in different populations of bighorn sheep across their range in North America (unpublished). Briefly, using fluorescently labeled microsatellite primers, microsatellites present in genomic DNA were amplified by PCR (polymerase chain reaction). PCR products were separated by capillary electrophoresis using an Applied Biosystems 3730 DNA Analyzer (Applied Biosystems, Foster City, CA), 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 (HWE) and allele frequencies were examined with GenePop .
Multi-locus heterozygosity (MLH) values were calculated for each animal as the preferred metric of genetic variability. Chapman et al.  showed that different genetic metrics (MLH, SH, IR) are highly correlated and non-independent, and they advocated the use of the simplest metric, MLH, in future studies of heterozygote-fitness correlations (HFC). For our multivariate analyses, MLH values were categorized into two categories using a cutoff > 0.5 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 (BRSV), bluetongue virus (BTV), and parainfluenza-3 virus (PI3). 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:20 were considered positive for BRSV, and titers > 1:16 were considered positive for PI3.
A VHF radiocollar (Telonics Inc., Lotek Wireless Inc., Advanced Telemetry Systems, Inc., and Telemetry Solutions) 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 48 hrs of capture. Resight of individual animals by radio-signal was conducted by state and federal wildlife biologists through February 2007 to detect mortalities and determine specific causes of death. Mortalities were investigated immediately and a field investigation and necropsy typically performed within 72 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 1, 2007.
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 (or end of study period if survived) and death due to suspected pneumonia. Independent variables evaluated for their impact on these three outcomes included sex, age-class (< 3 yr, 3 to < 6 yr, ≥ 6 yr), source population (RRWA vs KNWR), year released (2002, 2005), exposure status to each pathogen (BRSV, BTV, PI3) at time of release, and genetic variability (MLH). 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 (as a continuous variable) 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 (n = 50). Sheep that died from predation or unknown causes were excluded from this analysis. Variables were selected by backward stepwise elimination (likelihood ratio test P < 0.1) 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 < 0.1. Variables significantly associated with survival in these univariate analyses were evaluated for their relationship to failure rate (days until death) 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 ≤ 0.1 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 11.1 software (STATACorp, 4905 Lakeway Drive, College Station, Texas 77845 USA).
Characteristics of Source Populations
Bighorn sheep translocated from RRWA were predominantly male (18/20), while 72% (34/47) of the sheep translocated from KNWR were female (Fisher exact P < 0.001). Age class was not significantly associated with source population (Fisher exact test P = 0.068), but twice as many sheep in the 3-6 yr old range were captured in KNWR (23/36) compared to RRWA (11/31). All sheep relocated to SAM in 2005 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 2002.
Multi-locus heterozygosity (MLH) and prevalence of antibodies to pathogens among bighorn sheep from the Kofa National Wildlife Refuge (KNWR) in Arizona, the Red Rock Wildlife Area (RRWA) in New Mexico, and the San Andres Mountains (SAM) in New Mexico, USA.
n = 30
n = 20
n = 31
n = 9
Pathogen Exposure b
Genetic Diversity c
Three of the 33 microsatellite loci tested were monomorphic and were excluded from further analyses (BM4107, INRA063, TGLA94). Complete genotypes were determined for all animals at each of the remaining 30 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 (Table 1), was significantly higher in sheep captured in KNWR (mean = 0.579, 95% CI 0.555 - 0.603) compared to sheep captured in RRWA (mean = 0.394, 95% CI = 0.340-0.449). 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 25 sheep. Pneumonia was suspected as a cause of mortality in 17 sheep, and mountain lion predation was identified as the cause of death in 8 sheep. Specific cause of death could not be determined for 9 sheep, but predation and pneumonia were excluded as their cause of death. Deaths due to lion predation occurred in 2002, 2003, and 2006, while deaths due to pneumonia occurred every year (2002-2007). Most deaths due to pneumonia (65%; 11/17) occurred in the fall from September thru November, and seven of these occurred in a cluster in the fall of 2006. Pneumonia impacted all age classes with 27% (4/15) of 1-3 year old sheep, 35% (9/26) of 3-6 year old sheep and 44% (4/9) of > 6 year old sheep dying of pneumonia. Mean MLH was not significantly different among the 33 sheep that survived (mean = 0.50, 95% CI 0.45 - 0.55) compared to the 17 sheep that died of pneumonia (mean = 0.44, 95% CI 0.37 - 0.51), but this analyses lacked power due to low sample size (two-sided t test P = 0.161, power = 0.31). 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 (log rank test P = 0.079) in the univariate analyses. More than half of the sheep released at 1 to < 3 years old survived the entire study period (1,509 days), while median survival time was 1,234 days for sheep 3 years to < 6 years old and only 475 days for sheep released when 6 years or older. Genetic diversity and pathogen exposure were not related to survival time, even after stratifying by source population.
Factors significantly related to time to death in the final Cox proportional hazards model of survival of bighorn sheep translocated to the San Andres Mountain, New Mexico, USA.
95% Conf. Interval
Age 3 to <6 yr a
0.60 - 3.56
Age ≥6 yr a
1.22 - 9.32
BRSV antibodies b
0.14 - 0.97
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 [10, 20, 21], 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 (mean MLH = 0.39; Table 1) because it had been managed for >20 years as an inbred population. In contrast, the KNWR population was much more diverse (mean MLH = 0.57-0.60; Table 1), 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 (MLH) 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 (SH, IR) 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 (microsatellites) 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 [6, 9, 22, 23]. 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 (reduced gene flow), 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 [24, 25]. We found no evidence of such an effect among translocated adult sheep (survival was not associated with source population). Although we did not detect fitness differences among translocated animals, heterosis or genetic rescue effects may occur in F1 and later generations as the two source populations interbreed [25, 26]. While we were able to sample nine offspring born in the SAM for genetic and pathogen analyses (Table 1), 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% (17/34) 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 [4, 20, 27, 28]. Several bacteria (Pasteurella and Manheimia spp.) and viruses (BRSV, PI3) 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 (RRWA only) as well as BRSV and PI3 (both RRWA and KNWR; Table 1). 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 > 6 yrs of age (Table 2). 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 > 600 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.
This project was supported by the U.S. Fish and Wildlife Service, New Mexico Department of Game and Fish (NMDGF), White Sands Missile Range, and the National Fish and Wildlife Foundation. We appreciate the support and cooperation of Elise Goldstein, Eric Rominger, and Darrel Weybright from NMDGF; Grace Lee from the University of California, Davis, as well as personnel from Arizona Department of Game and Fish, Arizona Desert Bighorn Sheep Society, and the New Mexico Chapter of Foundation for North American Wild Sheep. This study was performed in compliance with state and federal regulations.
The findings and conclusions in this article are those of the author(s) and do not necessarily represent the views of the U.S. Fish and Wildlife Service.
- Acevedo-Whitehouse K, Cunningham AA: Is MHC enough for understanding wildlife immunogenetics?. Trends Ecology Evolution. 2006, 21: 433-438. 10.1016/j.tree.2006.05.010.View ArticlePubMedGoogle Scholar
- Smith KF, Acevedo-Whitehouse K, Pederson AB: The role of infectious disease in biological conservation. Animal Conservation. 2009, 12: 1-12. 10.1111/j.1469-1795.2008.00228.x.View ArticleGoogle Scholar
- Buechner HK: The bighorn sheep of the United States: its past, present, and future. Wildlife Monographs. 1960, 4: 1-174.Google Scholar
- Cassirer EF, Sinclair ARE: Dynamics of pneumonia in a bighorn sheep metapopulation. Journal Wildlife Management. 2007, 71: 1080-1088. 10.2193/2006-002.View ArticleGoogle Scholar
- Monello RJ, Murray DL, Cassirer EF: Ecological correlates of pneumonia epizootics in bighorn sheep herds. Canadian Journal Zoology. 2001, 79: 1423-1432. 10.1139/cjz-79-8-1423.View ArticleGoogle Scholar
- Boyce WM, Hedrick PW, Muggli-Cockett NE, Kalinowski S, Penedo MCT, Ramey RR: Genetic variation of major histocompatibility complex and microsatellite loci: a comparison in bighorn sheep. Genetics. 1996, 145: 421-433.Google Scholar
- Boyce WM, Weisenberger M: The rise and fall of psoroptic scabies in bighorn sheep in the San Andres Mountains, New Mexico. Journal Wildlife Diseases. 2005, 41: 525-531.View ArticleGoogle Scholar
- Rominger EM, Whitlaw HA, Weybright DL, Dunn WC, Ballard WB: The influence of mountain lion predation on bighorn sheep translocations. Journal Wildlife Management. 2004, 68: 993-999. 10.2193/0022-541X(2004)068[0993:TIOMLP]2.0.CO;2.View ArticleGoogle Scholar
- Gutiérrez-Espeleta GA, Kalinowski ST, Boyce WM, Hedrick PW: Genetic variation and population structure in desert bighorn sheep: implications for conservation. Conservation Genetics. 2000, 1: 3-15.View ArticleGoogle Scholar
- Whittaker DG, Ostermann SD, Boyce WM: Genetic variability of reintroduced California bighorn sheep in Oregon. Journal Wildlife Management. 2004, 68: 850-859. 10.2193/0022-541X(2004)068[0850:GVORCB]2.0.CO;2.View ArticleGoogle Scholar
- Hughes SS: STRand Nucleic Acids Analysis Software. 2000, University of California, Davis, CA, USAGoogle Scholar
- Raymond M, Rousset F: GENEPOP 1.02: population genetics software for exact tests and ecumenicism. Journal Heredity. 1995, 86: 248-249.Google Scholar
- Chapman JR, Nakagawa S, Coltman DW, Slate J, Sheldon BC: A quantitative review of heterozygote-fitness correlations in animal populations. Molecular Ecology. 2009, 18: 2746-2765. 10.1111/j.1365-294X.2009.04247.x.View ArticlePubMedGoogle Scholar
- Hayes CL, Rubin ES, Jorgensen MC, Botta RA, Boyce WM: Mountain lion predation of bighorn sheep in the Peninsular Ranges, California. Journal Wildlife Management. 2000, 64: 954-959. 10.2307/3803204.View ArticleGoogle Scholar
- Mantel N, Haenszel W: Statistical aspects of the analysis of data from retrospective studies of disease. Journal National Cancer Institute. 1959, 22: 719-748.Google Scholar
- Breslow N: Covariance analysis of censored survival data. Biometrics. 1974, 30: 89-99. 10.2307/2529620.View ArticlePubMedGoogle Scholar
- Caro TM, Laurenson MK: Ecological and genetic factors in conservation: a cautionary tale. Science. 1994, 263: 485-486. 10.1126/science.8290956.View ArticlePubMedGoogle Scholar
- Kinnison MT, Hendry AP, Stockwell CA: Contemporary evolution meets conservation biology II: impediments to integration and application. Ecological Research. 2007, 22: 947-954. 10.1007/s11284-007-0416-6.View ArticleGoogle Scholar
- Pedersen AB, Jones KE, Nunn CL, Altizer S: Infectious disease and extinction risk in wild mammals. Conservation Biology. 2007, 21: 1269-1279. 10.1111/j.1523-1739.2007.00776.x.View ArticlePubMedGoogle Scholar
- Singer FJ, Papouchis C, Symonds K: Translocations as a tool for restoring bighorn sheep. Restoration Ecology. 2000, 8: 6-13. 10.1046/j.1526-100x.2000.80061.x.View ArticleGoogle Scholar
- United States Fish and Wildlife Service: Recovery plan for bighorn sheep in the Peninsular Ranges, California. 2000, U.S. Fish and Wildlife Service, Portland, OregonGoogle Scholar
- Gutierrez-Espeleta GA, Hedrick PW, Kalinowski S, Garrigan D, Boyce WM: Is the decline of desert bighorn sheep from infectious disease the result of low MHC variation?. Heredity. 2001, 86: 439-450. 10.1046/j.1365-2540.2001.00853.x.View ArticlePubMedGoogle Scholar
- Forbes SH, Hogg JT, Buchanan FC, Crawford AM, Allendorf FW: Microsatellite evolution in congeneric mammals: domestic and bighorn sheep. Molecular Biology Evolution. 1995, 12: 1106-1113.PubMedGoogle Scholar
- Hedrick PW, Gutierrez-Espeleta GA, Lee RN: Founder effect in an island population of bighorn sheep. Molecular Ecology. 2001, 10: 851-857. 10.1046/j.1365-294X.2001.01243.x.View ArticlePubMedGoogle Scholar
- Tallman DA, Luikart G, Waples RS: The alluring simplicity and complex reality of genetic rescue. Trends Ecology Evolution. 2004, 19: 489-496. 10.1016/j.tree.2004.07.003.View ArticleGoogle Scholar
- Hogg JT, Forbes SH, Steele BM, Luikart G: Genetic rescue of an insular population of large mammals. Proceedings Royal Society B. 2006, 273: 1491-1499. 10.1098/rspb.2006.3477.View ArticleGoogle Scholar
- Gross JE, Singer FJ, Moses ME: Effect of disease, dispersal and area on bighorn sheep restoration. Restoration Ecology. 2000, 8: 25-37. 10.1046/j.1526-100x.2000.80063.x.View ArticleGoogle Scholar
- George JL, Martin DJ, Lukus PM, Miller MW: Epidemic pasteurellosis in a bighorn sheep population coinciding with the appearance of a domestic sheep. Journal Wildlife Diseases. 2008, 44: 388-403.View ArticleGoogle Scholar
- Rudolph KM, Hunter DL, Rimler RB, Cassirer EF, Foreyt WJ, Delong WJ, Weiser GC, Ward ACS: Microorganisms associated with a pneumonic epizootic in Rocky Mountain bighorn sheep (Ovis canadensis canadensis). Journal Zoo and Wildlife Medicine. 2007, 38: 548-558. 10.1638/2006-0027R.1.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.