Ancient hybridization patterns between bighorn and thinhorn sheep.

Whole-genome sequencing has advanced the study of species evolution, including the detection of genealogical discordant events such as ancient hybridization and incomplete lineage sorting (ILS). The evolutionary history of bighorn (Ovis canadensis) and thinhorn (Ovis dalli) sheep present an ideal system to investigate evolutionary discordance due to their recent and rapid radiation and putative secondary contact between bighorn and thinhorn sheep subspecies, specifically the dark pelage Stone sheep (O. dalli stonei) and predominately white Dall sheep (O. dalli dalli), during the last ice age. Here, we used multiple genomes of bighorn and thinhorn sheep, together with snow (O. nivicola) and the domestic sheep (O. aries) as outgroups, to assess their phylogenomic history, potential introgression patterns and their adaptive consequences. Among the Pachyceriforms (snow, bighorn and thinhorn sheep) a consistent monophyletic species tree was retrieved; however, many genealogical discordance patterns were observed. Alternative phylogenies frequently placed Stone and bighorn as sister clades. This relationship occurred more often and was less divergent than that between Dall and bighorn. We also observed many blocks containing introgression signal between Stone and bighorn genomes in which coat colour genes were present. Introgression signals observed between Dall and bighorn were more random and less frequent, and therefore probably due to ILS or intermediary secondary contact. These results strongly suggest that Stone sheep originated from a complex series of events, characterized by multiple, ancient periods of secondary contact with bighorn sheep.

Towards robust evolutionary inference with integral projection models.

Integral projection models (IPMs) are extremely flexible tools for ecological and evolutionary inference. IPMs track the distribution of phenotype in populations through time, using functions describing phenotype-dependent development, inheritance, survival and fecundity. For evolutionary inference, two important features of any model are the ability to (i) characterize relationships among traits (including values of the same traits across ages) within individuals, and (ii) characterize similarity between individuals and their descendants. In IPM analyses, the former depends on regressions of observed trait values at each age on values at the previous age (development functions), and the latter on regressions of offspring values at birth on parent values as adults (inheritance functions). We show analytically that development functions, characterized this way, will typically underestimate covariances of trait values across ages, due to compounding of regression to the mean across projection steps. Similarly, we show that inheritance, characterized this way, is inconsistent with a modern understanding of inheritance, and underestimates the degree to which relatives are phenotypically similar. Additionally, we show that the use of a constant biometric inheritance function, particularly with a constant intercept, is incompatible with evolution. Consequently, current implementations of IPMs will predict little or no phenotypic evolution, purely as artefacts of their construction. We present alternative approaches to constructing development and inheritance functions, based on a quantitative genetic approach, and show analytically and through an empirical example on a population of bighorn sheep how they can potentially recover patterns that are critical to evolutionary inference.

Demography, not inheritance, drives phenotypic change in hunted bighorn sheep.

Selective harvest, such as trophy hunting, can shift the distribution of a quantitative character such as body size. If the targeted character is heritable, then there will be an evolutionary response to selection, and where the trait is not, then any response will be plastic or demographic. Identifying the relative contributions of these different mechanisms is a major challenge in wildlife conservation. New mathematical approaches can provide insight not previously available. Here we develop a size- and age-based two-sex integral projection model based on individual-based data from a long-term study of hunted bighorn sheep (Ovis canadensis) at Ram Mountain, Canada. We simulate the effect of trophy hunting on body size and find that the inheritance of body mass is weak and that any perceived decline in body mass of the bighorn population is largely attributable to demographic change and environmental factors. To our knowledge, this work provides the first use of two-sex integral projection models to investigate the potential eco-evolutionary consequences of selective harvest.

Harnessing cross-species alignment to discover SNPs and generate a draft genome sequence of a bighorn sheep (Ovis canadensis).

BACKGROUND: Whole genome sequences (WGS) have proliferated as sequencing technology continues to improve and costs decline. While many WGS of model or domestic organisms have been produced, a growing number of non-model species are also being sequenced. In the absence of a reference, construction of a genome sequence necessitates de novo assembly which may be beyond the ability of many labs due to the large volumes of raw sequence data and extensive bioinformatics required. In contrast, the presence of a reference WGS allows for alignment which is more tractable than assembly. Recent work has highlighted that the reference need not come from the same species, potentially enabling a wide array of species WGS to be constructed using cross-species alignment. Here we report on the creation a draft WGS from a single bighorn sheep (Ovis canadensis) using alignment to the closely related domestic sheep (Ovis aries). RESULTS: Two sequencing libraries on SOLiD platforms yielded over 865 million reads, and combined alignment to the domestic sheep reference resulted in a nearly complete sequence (95% coverage of the reference) at an average of 12x read depth (104 SD). From this we discovered over 15 million variants and annotated them relative to the domestic sheep reference. We then conducted an enrichment analysis of those SNPs showing fixed differences between the reference and sequenced individual and found significant differences in a number of gene ontology (GO) terms, including those associated with reproduction, muscle properties, and bone deposition. CONCLUSION: Our results demonstrate that cross-species alignment enables the creation of novel WGS for non-model organisms. The bighorn sheep WGS will provide a resource for future resequencing studies or comparative genomics.

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.

Sex-based differences in the adaptive value of social behavior contrasted against morphology and environment.

The adaptive nature of sociality has long been a central question in ecology and evolution. However, the relative importance of social behavior for fitness, compared to morphology and environment, remains largely unknown. We assessed the importance of sociality for fitness (lamb production and survival) in a population of mark6d bighorn sheep (Ovis canadensis) over 16 years (n = 1022 sheep-years). We constructed social networks from observations (n = 38,350) of group membership (n = 3150 groups). We then tested whether consistent individual differences in social behavior (centrality) exist and evaluated their relative importance compared to factors known to affect fitness: mass, age, parental effects, and population density. Sheep exhibited consistent individual differences in social centrality. Controlling for maternal carryover effects and age, the positive effect of centrality in a social network on adult female lamb production and survival was equal or greater than the effect of body mass or population density. Social centrality had less effect on male survival and no effect on adult male lamb production or lamb survival. Through its effect on lamb production and survival, sociality in fission-fusion animal societies may ultimately influence population dynamics equally or more than morphological or environmental effects.

Estimating genome-wide heterozygosity: effects of demographic history and marker type.

Heterozygosity-fitness correlations (HFCs) are often used to link individual genetic variation to differences in fitness. However, most studies examining HFCs find weak or no correlations. Here, we derive broad theoretical predictions about how many loci are needed to adequately measure genomic heterozygosity assuming different levels of identity disequilibrium (ID), a proxy for inbreeding. We then evaluate the expected ability to detect HFCs using an empirical data set of 200 microsatellites and 412 single nucleotide polymorphisms (SNPs) genotyped in two populations of bighorn sheep (Ovis canadensis), with different demographic histories. In both populations, heterozygosity was significantly correlated across marker types, although the strength of the correlation was weaker in a native population compared with one founded via translocation and later supplemented with additional individuals. Despite being bi-allelic, SNPs had similar correlations to genome-wide heterozygosity as microsatellites in both populations. For both marker types, this association became stronger and less variable as more markers were considered. Both populations had significant levels of ID; however, estimates were an order of magnitude lower in the native population. As with heterozygosity, SNPs performed similarly to microsatellites, and precision and accuracy of the estimates of ID increased as more loci were considered. Although dependent on the demographic history of the population considered, these results illustrate that genome-wide heterozygosity, and therefore HFCs, are best measured by a large number of markers, a feat now more realistically accomplished with SNPs than microsatellites.

Sexually antagonistic association between paternal phenotype and offspring viability reinforces total selection on a sexually selected trait.

The evolution of conspicuous sexually selected traits, such as horns or antlers, has fascinated biologists for more than a century. Elaborate traits can only evolve if they substantially increase reproduction, because they probably incur survival costs to the bearer. Total selection on these traits, however, includes sexual selection on sires and viability selection on offspring and can be influenced by changes in each of these components. Non-random associations between paternal phenotype and offspring viability may thus affect total selection on sexually selected traits. Long-term data on wild bighorn sheep (Ovis canadensis) provide the first evidence in nature that association between paternal phenotype and lamb viability strengthens total selection on horn size of adult rams, a sexually selected trait. The association of paternal horn length and offspring viability was sexually antagonistic: long-horned males sired sons with high viability but daughters of low viability. These results shed new light on the evolutionary dynamics of an iconic sexually selected trait and have important implications for sustainable wildlife management.

Development of a high-density sub-species-specific targeted SNP assay for Rocky Mountain bighorn sheep (Ovis canadensis canadensis).

Due to their abundance and relative ease of genotyping, single nucleotide polymorphisms (SNPs) are a commonly used molecular marker for contemporary population genetic and genomic studies. A high-density and cost-effective way to type SNP loci is Allegro targeted genotyping (ATG), which is a form of targeted genotyping by sequencing developed and offered by Tecan genomics. One major drawback of this technology is the need for a reference genome and information on SNP loci when designing a SNP assay. However, for some non-model species genomic information from other closely related species can be used. Here we describe our process of developing an ATG assay to target 50,000 SNPs in Rocky Mountain bighorn sheep, using a reference genome from domestic sheep and SNP resources from prior bighorn sheep studies. We successfully developed a high accuracy, high-density, and relatively low-cost SNP assay for genotyping Rocky Mountain bighorn sheep that genotyped ~45,000 SNP loci. These loci were relatively evenly distributed throughout the genome. Furthermore, the assay produced genotypes at tens of thousands of SNP loci when tested on other mountain sheep species and subspecies.

Assessing shared respiratory pathogens between domestic (Ovis aries) and bighorn (Ovis canadensis) sheep; methods for multiplex PCR, amplicon sequencing, and bioinformatics to characterize respiratory flora.

Respiratory disease is responsible for dramatic population declines in bighorn sheep (Ovis canadensis), and respiratory pathogen diagnostics contribute to the management of bighorn populations. To create a comprehensive and consistent approach to bighorn sheep respiratory diagnostics, we created a culture-independent assay to detect and strain type Mannheimia haemolytica, Bibersteinia trehalosi, Pasteurella multocida, and Mycoplasma ovipneumoniae. The assay also detects and characterizes the Pasteurellaceae leukotoxin A gene, and broadly assesses the bacterial composition of each sample based on 16S rRNA sequences. The assay is based on a three-step approach: 1) Multiplex PCR to amplify targets including eight loci for each bacterial species, the Pasteurellaceae lktA gene, and the 16S rRNA gene 2) Library preparation, barcoding, and short-read Illumina sequencing to determine the genetic sequences of each target, and 3) Bioinformatics in the form of automated software to analyze genetic sequences. The assay was designed to assess shared pathogens between domestic and bighorn sheep, but could be useful for many applications in bighorn sheep respiratory disease research and management.

Population genomics as a new tool for wildlife management.

Admixture and introgression have varied effects on population viability and fitness. Admixture might be an important source of new alleles, particularly for small, geographically isolated populations. However, admixture might also cause outbreeding depression if populations are adapted to different ecological or climatic conditions. Because of the emerging use of translocation and admixture as a conservation and wildlife management strategy to reduce genetic load (termed genetic rescue), the possible effects of admixture have practical consequences (Bouzat et al. 2009; Hedrick & Fredrickson 2010). Importantly, genetic load and local adaptation are properties of individual loci and epistatic interactions among loci rather than properties of genomes. Likewise, the outcome and consequences of genetic rescue depend on the fitness effects of individual introduced alleles. In this issue of Molecular Ecology, Miller et al. (2012) use model-based, population genomic analyses to document locus-specific effects of a recent genetic rescue in the bighorn sheep population within the National Bison Range wildlife refuge (NBR; Montana, USA). They find a subset of introduced alleles associated with increased fitness in NBR bighorn sheep, some of which experienced accelerated introgression following their introduction. These loci mark regions of the genome that could constitute the genetic basis of the successful NBR bighorn sheep genetic rescue. Although population genomic analyses are frequently used to study local adaptation and selection (e.g. Hohenlohe et al. 2010; Lawniczak et al. 2010), this study constitutes a novel application of this analytical framework for wildlife management. Moreover, the detailed demographic data available for the NBR bighorn sheep population provide a rare and powerful source of information and allow more robust population genomic inference than is often possible.

Genomic consequences of genetic rescue in an insular population of bighorn sheep (Ovis canadensis).

Genetic rescue is a management intervention whereby a small population is supplemented with individuals from other populations in an attempt to reverse the effects of inbreeding and increased genetic load. One such rescue was recently documented in the population of bighorn sheep (Ovis canadensis) within the National Bison Range wildlife refuge (Montana, USA). Here, we examine the locus-specific effects of rescue in this population using a newly developed genome-wide set of 195 microsatellite loci and first-generation linkage map. We found that the rate of introgression varied among loci and that 111 loci, 57% of those examined, deviated from patterns of neutral inheritance. The most common deviation was an excess of homozygous genotypes relative to neutral expectations, indicative of directional selection. As in previous study of this rescue, individuals with more introduced alleles had higher reproductive success and longevity. In addition, we found 30 loci, distributed throughout the genome, which seem to have individual effects on these life history traits. Although the potential for outbreeding depression is a major concern when translocating individuals between populations, we found no evidence of such effects in this population.

QTL mapping for sexually dimorphic fitness-related traits in wild bighorn sheep.

Dissecting the genetic architecture of fitness-related traits in wild populations is key to understanding evolution and the mechanisms maintaining adaptive genetic variation. We took advantage of a recently developed genetic linkage map and phenotypic information from wild pedigreed individuals from Ram Mountain, Alberta, Canada, to study the genetic architecture of ecologically important traits (horn volume, length, base circumference and body mass) in bighorn sheep. In addition to estimating sex-specific and cross-sex quantitative genetic parameters, we tested for the presence of quantitative trait loci (QTLs), colocalization of QTLs between bighorn sheep and domestic sheep, and sex × QTL interactions. All traits showed significant additive genetic variance and genetic correlations tended to be positive. Linkage analysis based on 241 microsatellite loci typed in 310 pedigreed animals resulted in no significant and five suggestive QTLs (four for horn dimension on chromosomes 1, 18 and 23, and one for body mass on chromosome 26) using genome-wide significance thresholds (Logarithm of odds (LOD) >3.31 and >1.88, respectively). We also confirmed the presence of a horn dimension QTL in bighorn sheep at the only position known to contain a similar QTL in domestic sheep (on chromosome 10 near the horns locus; nominal P<0.01) and highlighted a number of regions potentially containing weight-related QTLs in both species. As expected for sexually dimorphic traits involved in male-male combat, loci with sex-specific effects were detected. This study lays the foundation for future work on adaptive genetic variation and the evolutionary dynamics of sexually dimorphic traits in bighorn sheep.

Short reads, circular genome: skimming solid sequence to construct the bighorn sheep mitochondrial genome.

As sequencing technology improves, an increasing number of projects aim to generate full genome sequence, even for nonmodel taxa. These projects may be feasibly conducted at lower read depths if the alignment can be aided by previously developed genomic resources from a closely related species. We investigated the feasibility of constructing a complete mitochondrial (mt) genome without preamplification or other targeting of the sequence. Here we present a full mt genome sequence (16,463 nucleotides) for the bighorn sheep (Ovis canadensis) generated though alignment of SOLiD short-read sequences to a reference genome. Average read depth was 1240, and each base was covered by at least 36 reads. We then conducted a phylogenomic analysis with 27 other bovid mitogenomes, which placed bighorn sheep firmly in the Ovis clade. These results show that it is possible to generate a complete mitogenome by skimming a low-coverage genomic sequencing library. This technique will become increasingly applicable as the number of taxa with some level of genome sequence rises.

Rapid decrease in horn size of bighorn sheep: environmental decline, inbreeding depression, or evolutionary response to trophy hunting?

There are numerous examples demonstrating that selection has greatly influenced phenotypes in wild-harvested species. Here, a significant reduction in horn size in trophy desert bighorn sheep rams over 30 years in a reintroduced population in Aravaipa Canyon, Arizona is documented. After examining the potential effects of a detrimental change in the environment, inbreeding depression, and hunter-caused evolutionary change, it appears that environmental deterioration, apparently from the effects of drought, may be a major cause of the decline in horn size. In particular, the reduction in ram horn size is positively associated with reduced winter lifetime rainfall over the 3 decades. Over the same period, the demographic indicator lamb-to-ewe ratio has also declined in the Aravaipa population. On the other hand, lamb-to-ewe ratio has not declined statewide in Arizona, and the population size in Aravaipa appears to be increasing, suggesting local- and trait-specific effects. Using a theoretical context, neither inbreeding depression nor hunter selection by themselves appear to the sole causes of the lower horn size. However, some combination of environmental factors, inbreeding depression, and hunter selection may have caused the decrease in observed horn size. It is not clear what management actions might be successful in countering the environmental effects on horn size, but supplemental feeding and cattle removal are suggested while translocation is suggested to counter the effects of inbreeding depression and reduced hunting and translocation are suggested to counter the effects of hunter selection.

Whole-genome sequence analysis unveils different origins of European and Asiatic mouflon and domestication-related genes in sheep.

The domestication and subsequent development of sheep are crucial events in the history of human civilization and the agricultural revolution. However, the impact of interspecific introgression on the genomic regions under domestication and subsequent selection remains unclear. Here, we analyze the whole genomes of domestic sheep and their wild relative species. We found introgression from wild sheep such as the snow sheep and its American relatives (bighorn and thinhorn sheep) into urial, Asiatic and European mouflons. We observed independent events of adaptive introgression from wild sheep into the Asiatic and European mouflons, as well as shared introgressed regions from both snow sheep and argali into Asiatic mouflon before or during the domestication process. We revealed European mouflons might arise through hybridization events between a now extinct sheep in Europe and feral domesticated sheep around 6000-5000 years BP. We also unveiled later introgressions from wild sheep to their sympatric domestic sheep after domestication. Several of the introgression events contain loci with candidate domestication genes (e.g., PAPPA2, NR6A1, SH3GL3, RFX3 and CAMK4), associated with morphological, immune, reproduction or production traits (wool/meat/milk). We also detected introgression events that introduced genes related to nervous response (NEURL1), neurogenesis (PRUNE2), hearing ability (USH2A), and placental viability (PAG11 and PAG3) into domestic sheep and their ancestral wild species from other wild species.

Genetic linkage map of a wild genome: genomic structure, recombination and sexual dimorphism in bighorn sheep.

BACKGROUND: The construction of genetic linkage maps in free-living populations is a promising tool for the study of evolution. However, such maps are rare because it is difficult to develop both wild pedigrees and corresponding sets of molecular markers that are sufficiently large. We took advantage of two long-term field studies of pedigreed individuals and genomic resources originally developed for domestic sheep (Ovis aries) to construct a linkage map for bighorn sheep, Ovis canadensis. We then assessed variability in genomic structure and recombination rates between bighorn sheep populations and sheep species. RESULTS: Bighorn sheep population-specific maps differed slightly in contiguity but were otherwise very similar in terms of genomic structure and recombination rates. The joint analysis of the two pedigrees resulted in a highly contiguous map composed of 247 microsatellite markers distributed along all 26 autosomes and the X chromosome. The map is estimated to cover about 84% of the bighorn sheep genome and contains 240 unique positions spanning a sex-averaged distance of 3051 cM with an average inter-marker distance of 14.3 cM. Marker synteny, order, sex-averaged interval lengths and sex-averaged total map lengths were all very similar between sheep species. However, in contrast to domestic sheep, but consistent with the usual pattern for a placental mammal, recombination rates in bighorn sheep were significantly greater in females than in males (~12% difference), resulting in an autosomal female map of 3166 cM and an autosomal male map of 2831 cM. Despite differing genome-wide patterns of heterochiasmy between the sheep species, sexual dimorphism in recombination rates was correlated between orthologous intervals. CONCLUSIONS: We have developed a first-generation bighorn sheep linkage map that will facilitate future studies of the genetic architecture of trait variation in this species. While domestication has been hypothesized to be responsible for the elevated mean recombination rate observed in domestic sheep, our results suggest that it is a characteristic of Ovis species. However, domestication may have played a role in altering patterns of heterochiasmy. Finally, we found that interval-specific patterns of sexual dimorphism were preserved among closely related Ovis species, possibly due to the conserved position of these intervals relative to the centromeres and telomeres. This study exemplifies how transferring genomic resources from domesticated species to close wild relative can benefit evolutionary ecologists while providing insights into the evolution of genomic structure and recombination rates of domesticated species.