A genome-wide set of SNPs detects population substructure and long range linkage disequilibrium in wild sheep

The development of genomic resources for wild species is still in its infancy. However, cross-species utilization of technologies developed for their domestic counterparts has the potential to unlock the genomes of organisms that currently lack genomic resources. Here, we apply the OvineSNP50 BeadChip, developed for domestic sheep, to two related wild ungulate species: the bighorn sheep (Ovis canadensis) and the thinhorn sheep (Ovis dalli). Over 95% of the domestic sheep markers were successfully genotyped in a sample of fifty-two bighorn sheep while over 90% were genotyped in two thinhorn sheep. Pooling the results from both species identified 868 single-nucleotide polymorphisms (SNPs), 570 were detected in bighorn sheep, while 330 SNPs were identified in thinhorn sheep. The total panel of SNPs was able to discriminate between the two species, assign population of origin for bighorn sheep and detect known relationship classes within one population of bighorn sheep. Using an informative subset of these SNPs (n=308), we examined the extent of genome-wide linkage disequilibrium (LD) within one population of bighorn sheep and found that high levels of LD persist over 4 Mb.     

Re-sequencing regions of the ovine Y chromosome in domestic and wild sheep reveals novel paternal haplotypes

The male-specific region of the ovine Y chromosome (MSY) remains poorly characterized, yet sequence variants from this region have the potential to reveal the wild progenitor of domestic sheep or examples of domestic and wild paternal introgression. The 5' promoter region of the sex-determining gene SRY was re-sequenced using a subset of wild sheep including bighorn ( Ovis canadensis ), thinhorn ( Ovis dalli spp.), urial ( Ovis vignei ), argali ( Ovis ammon ), mouflon ( Ovis musimon ) and domestic sheep ( Ovis aries ). Seven novel SNPs ( oY 2– oY 8) were revealed; these were polymorphic between but not within species. Re-sequencing and fragment analysis was applied to the MSY microsatellite SRYM18 . It contains a complex compound repeat structure and sequencing of three novel size fragments revealed that a pentanucleotide element remained fixed, whilst a dinucleotide element displayed variability within species. Comparison of the sequence between species revealed that urial and argali sheep grouped more closely to the mouflon and domestic breeds than the pachyceriforms (bighorn and thinhorn). SNP and microsatellite data were combined to define six previously undetected haplotypes. Analysis revealed the mouflon as the only species to share a haplotype with domestic sheep, consistent with its status as a feral domesticate that has undergone male-mediated exchange with domestic animals. A comparison of the remaining wild species and domestic sheep revealed that O. aries is free from signatures of wild sheep introgression.     

Restoration of a bighorn sheep population impeded by Mycoplasma ovipneumoniae exposure

Bighorn sheep (Ovis canadensis) were once extirpated from the Black Hills region of South Dakota, U.S.A., mirroring declining populations throughout North America. Since the 1960s, several reintroductions have occurred in the Black Hills to reestablish populations, with varying success. We translocated 26 bighorn sheep from Alberta, Canada, to the Black Hills (February 2015) to restore bighorn sheep to their historic range. Due to prior examinations of cause‐specific survival, subsequent genetic diversity and disease prevalence analyses were required to evaluate success of the restoration effort. We measured a mean allelic diversity of 5.23 (SE = 0.44 [mean number of alleles]) and an observed heterozygosity of 0.71 (SE = 0.06; expected = 0.64 ± 0.05) in the translocated individuals. Translocated bighorn sheep tested negative for Mycoplasma ovipneumoniae at capture. An autogenous vaccine was administered prior to release in an attempt to safeguard the translocated bighorn sheep from infection with a strain known to be resident in adjacent bighorn sheep populations. However, the year following the translocation, a different strain of M. ovipneumoniae was associated with a pneumonia outbreak that resulted in 57.9% mortality. Our results suggest that allelic diversity and heterozygosity were sufficient for long‐term herd establishment, reducing the potential for founder effects. However, the overwhelming mortality associated with pneumonia, via the transfer of M. ovipneumoniae from an unknown source, limited the success or our reintroduction efforts. Successful attempts to restore bighorn sheep to their historic ranges must consider and mitigate potential routes for M. ovipneumoniae transmission pre‐ and post‐reintroduction.     

Bighorn Sheep Genetic Structure in Wyoming Reflects Geography and Management

Aligning wildlife management boundaries with accurate biological units promotes effective conservation and management practices that reflect ecological and evolutionary processes. Neutral genetic markers allow for quantitative delineation of population structure without a priori assumptions or biases. In the United States, bighorn sheep (Ovis canadensis) are a charismatic component of Wyoming's biodiversity and a species that provides important viewing and hunting opportunities. Bighorn sheep abundances are relatively stable throughout Wyoming, and the species is managed by administrative units identified using expert knowledge, distribution and movement data, and geographic and administrative boundaries. We used a panel of 38 variable microsatellite loci and 512 base pairs of mitochondrial DNA sequence to identify the genetic structure throughout the state and in translocation source herds, quantify the extent of genetic diversity within each genetic cluster, and estimate the degree of gene flow among herds using blood and tissue samples collected 1989–2017. We identified genetic structure of Rocky Mountain bighorn sheep in the major mountain ranges of Wyoming, with strong support for ≥5 genetic clusters using microsatellite loci. These genetic clusters generally aligned with current management units, whereas mitochondrial data showed a more complex mosaic that was not geographically patterned. Genetic variation estimated from both markers was high within each herd and comparable among herds. The assignment of individuals reflected a combination of geographic isolation and translocation, which has been extensive. Our results provide a state-wide assessment of genetic diversity and structure that will enhance management by understanding the outcomes of translocation, identifying the source of unknown individuals, and parameterizing disease ecology models. © 2020 The Wildlife Society.     

Desert bighorn sheep habitat selection,group size,and mountain lion predation risk

Bighorn sheep (Ovis canadensis) evolved for thousands of years in the presence of numerous predators, including mountain lions (Puma concolor). Bighorn sheep have presumably developed predator avoidance strategies; however, the effectiveness of these strategies in reducing risk of mountain lion predation is not well understood. These strategies are of increasing interest because mountain lion predation on bighorn sheep has been identified as a leading cause of mortality in some sheep populations. Therefore, we investigated how mountain lions affect both bighorn sheep habitat selection and risk of mortality in Arizona, USA. We used 2 approaches to investigate the predator-prey relationship between mountain lions and bighorn sheep. We fit 103 bighorn sheep (81 females and 22 males) with global positioning system radio-collars in 2 Arizona populations from 2013 to 2017, and used a negative binomial resource selection probability function to evaluate whether bighorn sheep selected for habitat features in accordance with presumed predator avoidance strategies, including terrain ruggedness, slope, topographic position, and horizontal obstruction, in 2 seasons (winter and summer). We then estimated how habitat features such as terrain ruggedness, slope, horizontal obstruction, and group size, influence the risk of mortality due to mountain lion predation using an Andersen-Gill proportional hazards model. Generally, both sexes selected areas with lower horizontal obstruction and intermediate ruggedness and slope, but selection patterns differed between seasons and sexes. The use of more rugged areas and steeper slopes decreased the risk of mortality due to mountain lion predation, consistent with presumed predator avoidance strategies. Increased group size decreased risk of bighorn sheep mortality due to mountain lion predation but this effect became marginal at approximately 10 individuals/group. We did not identify a relationship between horizontal obstruction and bighorn sheep mortality risk. Our findings can be used in habitat and population management decisions such as the prioritization of habitat restoration sites or selection of translocation sites. In addition, we suggest that augmentation of low-density bighorn sheep populations may reduce mountain lion predation risk by increasing group size, and that releasing large groups of bighorn sheep in population augmentation and reintroduction efforts may help to reduce mountain lion predation.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Development, characterization and transferability of microsatellite markers for the plant genus Phlox (Polemoniaceae)

In order to study diversification and microevolution in Phlox, we developed nine polymorphic microsatellite loci. In 20 individuals of Phlox pilosa from a single population, the average number of alleles per locus was 10.0 ± 5.1, and average observed and expected heterozygosities were 0.611 ± 0.234 and 0.769 ± 0.170, respectively. Most of these markers amplified successfully in 11 additional species of Phlox, representing a broad diversity of the genus, and some also amplified in more distantly related members of the Polemoniaceae. These microsatellite markers will be valuable for investigation of evolutionary processes in this important study system.     

Globally dispersed Y chromosomal haplotypes in wild and domestic sheep

To date, investigations of genetic diversity and the origins of domestication in sheep have utilised autosomal microsatellites and variation in the mitochondrial genome. We present the first analysis of both domestic and wild sheep using genetic markers residing on the ovine Y chromosome. Analysis of a single nucleotide polymorphism (oY1) in the SRY promoter region revealed that allele A-oY1 was present in all wild bighorn sheep (Ovis canadensis), two subspecies of thinhorn sheep (Ovis dalli), European Mouflon (Ovis musimon) and the Barbary (Ammontragis lervia). A-oY1 also had the highest frequency (71.4%) within 458 domestic sheep drawn from 65 breeds sampled from Africa, Asia, Australia, the Caribbean, Europe, the Middle East and Central Asia. Sequence analysis of a second locus, microsatellite SRYM18, revealed a compound repeat array displaying fixed differences, which identified bighorn and thinhorn sheep as distinct from the European Mouflon and domestic animals. Combined genotypic data identified 11 male-specific haplotypes that represented at least two separate lineages. Investigation of the geographical distribution of each haplotype revealed that one (H6) was both very common and widespread in the global sample of domestic breeds. The remaining haplotypes each displayed more restricted and informative distributions. For example, H5 was likely founded following the domestication of European breeds and was used to trace the recent transportation of animals to both the Caribbean and Australia. A high rate of Y chromosomal dispersal appears to have taken place during the development of domestic sheep as only 12.9% of the total observed variation was partitioned between major geographical regions.     

Infectious keratoconjunctivitis in bighorn sheep, Silver Bell Mountains, Arizona, USA

An infectious keratoconjunctivitis (IKC) epizootic in bighorn sheep (Ovis canadensis) occurred in the Silver Bell Mountains, Arizona, USA, from 1 December 2003 to 31 March 2004. We used standard culture methods and polymerase chain reaction (PCR) amplification of the 16S rRNA gene to test for the causative agents of IKC and other diseases reported to be associated with bighorn sheep populations. All bighorn sheep and domestic goat test results were negative except for Mycoplasma spp. and Branhamella spp. The culture and PCR results differed. Conjunctival swabs from four of 19 IKC-affected bighorn sheep tested by culture were positive for Mycoplasma spp., whereas 22 of 22 bighorn sheep samples tested by PCR were positive for Mycoplasma spp. None of 13 domestic goats tested positive by culture for Mycoplasma spp., whereas five of 16 tested positive by PCR. Three of 16 domestic goats and seven of 24 IKC-affected bighorn sheep tested positive for Branhamella spp. by culture. Bighorn sheep began showing clinical signs of IKC between 21 and 28 days following initial detection of domestic goats in bighorn sheep habitat. The IKC epizootic lasted 122 days, and individual bighorn sheep were blind for an average of 38.4 days. Given the clear potential for disease transmission to bighorn sheep, we recommend that land managers not allow the pasturing of domestic goats near occupied bighorn sheep habitat.     

Habitat use by desert bighorn sheep during lambing

Habitat used for predator escape may be a factor limiting restoration of desert bighorn sheep (Ovis canadensis mexicana) in small mountain ranges in New Mexico. Female bighorn sheep seek isolation in discrete areas for parturition. Although parturition sites are used only for <3 days they play an important role in neonate survival. We compared habitat characteristics at pre-parturition (n = 21), parturition (n = 38), random (n = 38), and post-parturition sites (n = 21). At each site we calculated distance to steep terrain, elevation, ruggedness, slope, and visibility. Parturition and post-parturition sites were higher in elevation and more rugged than sites used during pre-parturition. Post-parturition sites were closer to terrain with 100% slope than the pre-parturition or parturition sites. Post-parturition sites had higher visibility and steeper slopes than the pre-parturition sites. Parturition sites were steeper, higher in elevation, more rugged, and had lower visibility than the random sites. Fidelity to parturition sites was observed on 1 of 38 occasions. Because lambing sites are used for short periods they are poorly understood and often ignored by land managers. Additional research is recommended in other desert bighorn sheep populations to understand better and predict habitat use during the lambing period.     

Comparative aspects of biochemical polymorphism in the blood of Caprinae species and their hybrids

Biochemical variation at 14 blood loci was reviewed, and specific features compared experimentally in sheep Ovis aries, mouflon Ovis musimon, goat Capra hircus, aoudad Ammotragus lervia and in 2 stillborn aoudad × goat hybrids. Variation at 3 loci was also studied in dall sheep Ovis dalli, bighorn sheep Ovis canadensis and rocky mountain goat Oreamnos americanus. Haemoglobin C production in an anaemic Hb AB mouflon and in mouflon × sheep hybrids was examined. Mouflon differ from domestic sheep in that synthesis of both Hb β^AHb β^Bchains is switched off during Hb C production. The mouflon × sheep hybrids switched off one or both chains depending on whether they had inherited sheep or mouflon Hb β chain genes. In general aoudad showed a closer affinity to goats than to sheep.     

Microsatellite evolution in congeneric mammals: domestic and bighorn sheep

We compared genotypes at eight (AC)n microsatellite loci in domestic sheep (Ovis aries) and wild Rocky Mountain bighorn sheep (O. canadensis). The domestic sheep had greater genetic variation, higher allele-size variances, and larger allele sizes than the wild sheep. Accumulating evidence from higher taxonomic comparisons shows that these parameters are biased if microsatellite loci are selected in one taxon and used in another. Our results demonstrate similar biases between congeneric species. We compared standard measures of genetic variation, differentiation, and distance within and between species (H, D, FST) to newer measures based on allele-size variance (SW, SB, RST). The size-based distances better detected species-level divergence, but standard measures better distinguished allopatric populations. Empirical calibration of these measures at the subspecies level is needed to establish their useful ranges.      

Using Genetic Tools to Track Desert Bighorn Sheep Colonizations

ABSTRACT Understanding colonization is vital for managing fragmented populations. We employed mitochondrial DNA haplotypes and 14 microsatellite (nuclear DNA) markers to infer the origins of newly established populations of desert bighorn sheep (Ovis canadensis nelsoni) and to assess loss of genetic diversity during natural colonizations. We used haplotype distribution, F-statistics, Bayesian population clustering, and assignment tests to infer source populations for 3 recent colonies and identified a previously undetected colonization from multiple source populations. Allelic richness declined in 3 of 4 colonies in comparison to the primary source populations, but not as much as has been reported for translocated populations. Heterozygosity declined in only one colony. We also demonstrated that both native and translocated desert bighorn sheep have naturally recolonized empty habitats and suggest that colonization may partially offset population extinction in the region as long as connectivity is maintained. Genetic techniques and mitochondrial DNA haplotypes we described will allow managers to determine the origins of future colonizations by bighorn sheep in California, USA, and prioritize protection of linkages between known sources and colonies.     

Population response of reintroduced bighorn sheep after observed commingling with domestic sheep

Bighorn sheep (Ovis canadensis) often die from respiratory disease after commingling with domestic sheep. From 2000 to 2009, we observed commingling between domestic and reintroduced bighorn sheep in 3 populations in UT, USA. We investigated how commingling affected survival of radio-collared female bighorns that were released initially (founder) and those that were subsequently released (augmented). We predicted that the proportion of young surviving to their first winter and population growth would be lower after observed commingling with domestic sheep. We observed groups of bighorns year-round on 2,712 occasions and commingling between domestic sheep and bighorns in 6 instances. On Mount Timpanogos, survival rates were best modeled as constant for females (n?=?57) before and after observed commingling with domestic sheep. Survival rates of female bighorns, however, decreased significantly in Rock Canyon (n?=?21) and on Mount Nebo (n?=?22) for founder, but not augmented bighorns after observed commingling with domestic sheep. Also, the proportion of young surviving to their first winter was almost 3 times lower and population growth was reduced for bighorns after observed commingling with domestic sheep in Rock Canyon and on Mount Nebo. Commingling between domestic and bighorn sheep reduced population parameters in 2 of 3 bighorn populations we studied; however, on Mount Timpanogos, interactions between those 2 species were not fatal for radio-collared female bighorns. Wildlife biologists should manage for spatial separation of these 2 species and consider the location of hobby farms and trailing operations of domestic sheep near release sites for bighorns.     

Disease and secondary sexual traits: effects of pneumonia on horn size of bighorn sheep

Secondary sexual traits (e.g., horns and antlers) have ecological and evolutionary importance and are of management interest for game species. Yet, how these traits respond to emerging threats like infectious disease remains underexplored. Infectious pneumonia threatens bighorn sheep (Ovis canadensis) populations across North America and we hypothesized it may also reduce horn growth in male sheep. We assess the effect of pneumonia on horn size in male bighorn sheep using 12 herd datasets from across the western United States that had horn growth and disease data. Disease resulted in 12–35% reduction in increment (yearly) length and 3–13% reduction in total horn length in exposed individuals. The disease effect was prolonged when pathogens continued to circulate in sheep populations. Further, disease likely delays the age at which horns reach ¾-curl and prevents achievement of full-curl. This is further evidenced with 6 of the 12 herds experiencing an increase in average age at harvest following die-off events.     

Polymorphic microsatellite markers for the critically endangered Balearic shearwater, Puffinus mauretanicus

Ten novel polymorphic microsatellite loci were isolated and characterized from the Balearic shearwater (Puffinus mauretanicus), a critically endangered seabird. The developed loci revealed a relatively low number of alleles per locus, as well as low levels of polymorphism (H_O = 0.377 ± 0.241). One of the loci appeared to be W‐linked. All polymorphic loci were successfully amplified in its closely related species, the Yelkouan shearwater (Puffinus yelkouan). These microsatellite markers would be useful for assessing population structure in the Balearic shearwater and the possible hybridization process between both shearwaters species.     

Characterization of 29 polymorphic artiodactyl microsatellite markers for the mountain goat (Oreamnos americanus)

We report the results of a cross‐species amplification test of 156 bovine, ovine and cervid microsatellite markers in a wild population of mountain goats, Oreamnos americanus, inhabiting Caw Ridge, Alberta, Canada. Twenty‐nine markers were found to be low to moderately polymorphic with between two to nine alleles per locus. Observed heterozygosity ranged from 0.14 to 0.85 for a sample of 215 mountain goats. This set of markers will be used in parentage analyses to construct the pedigree of the long‐term studied population and to investigate the effects of individual genetic variability on life‐history traits.     

Survival of bighorn sheep (Ovis canadensis) commingled with domestic sheep (Ovis aries) in the absence of Mycoplasma ovipneumoniae

To test the hypothesis that Mycoplasma ovipneumoniae is an important agent of the bighorn sheep (Ovis canadensis) pneumonia that has previously inevitably followed experimental commingling with domestic sheep (Ovis aries), we commingled M. ovipneumoniae-free domestic and bighorn sheep (n=4 each). One bighorn sheep died with acute pneumonia 90 days after commingling, but the other three remained healthy for >100 days. This unprecedented survival rate is significantly different (P=0.002) from that of previous bighorn-domestic sheep contact studies but similar to (P>0.05) bighorn sheep survival following commingling with other ungulates. The absence of epizootic respiratory disease in this experiment supports the hypothesized role of M. ovipneumoniae as a key pathogen of epizootic pneumonia in bighorn sheep commingled with domestic sheep.     

Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

The confounding effects of population structure complicate efforts to identify regions of the genome under the influence of selection in natural populations. Here we test for evidence of selection in three genes involved in vertebrate immune function - the major histocompatibility complex (MHC), interferon gamma (IFNG) and natural resistance associated macrophage polymorphism (NRAMP) - in highly structured populations of wild thinhorn sheep (Ovis dalli). We examined patterns of variation at microsatellite loci linked to these gene regions and at the DNA sequence level. Simple Watterson's tests indicated balancing selection at all three gene regions. However, evidence for selection was confounded by population structure, as the Watterson's test statistics from linked markers were not outside of the range of values from unlinked and presumably neutral microsatellites. The translated coding sequences of thinhorn IFNG and NRAMP are fixed and identical to those of domestic sheep (Ovis aries). In contrast, the thinhorn MHC DRB locus shows significant evidence of overdominance through both an excess of nonsynonymous substitution and trans-species polymorphism. The failure to detect balancing selection at microsatellite loci linked to the MHC is likely the result of recombination between the markers and expressed gene regions.