Soay Sheep: The Back-story

The primitive race of Soay sheep from the St Kilda archipelago in northwest Scotland has played an important role in narratives of the history of domestic sheep. The Soays, apparently a ‘Bronze Age’ race of sheep, were probably confined to the precipitous isle of Soay as soon as ‘Iron Age’ sheep were introduced to Hirta, St Kilda's main island, owing to the competitive edge of the ferocious Soay rams over the new arrivals. In the 1880s, Pitt-Rivers, following his archaeozoological interests, was the first to keep Soays in his park, their epic journey from the edge of the Atlantic to southern England enabled by his acquaintance with their owner. In the early twentieth century, Soays featured in animal bone reports for archaeological sites, were kept in parks and involved in breeding experiments, particularly around Horsham, Sussex (where their owner lived), and in Edinburgh. The transfer of 107 Soays to Hirta in 1934 and 1935, after humans had evacuated St Kilda, was a remarkable feat, enabling the important long-running Soay Sheep Project. The historical exploitation of ‘feral’ Soay sheep by the islanders of St Kilda has significant cultural ramifications.

MHC class IIa haplotypes derived by high throughput SNP screening in an isolated sheep population

Investigating the current evolutionary processes acting on a highly polymorphic gene region, such as the major histocompatibility complex (MHC), requires extensive population data for both genotypes and phenotypes. The MHC consists of several tightly linked loci with both allelic and gene content variation, making it challenging to genotype. Eight class IIa haplotypes have previously been identified in the Soay sheep (Ovis aries) of St. Kilda using Sanger sequencing and cloning, but no single locus is representative of all haplotypes. Here, we exploit the closed nature of the island population of Soay sheep and its limited haplotypic variation to identify a panel of SNPs that enable imputation of MHC haplotypes. We compared MHC class IIa haplotypes determined by Sanger sequence-based genotyping of 135 individuals to their SNP profiles generated using the Ovine Infinium HD BeadChip. A panel of 11 SNPs could reliably determine MHC diplotypes, and two additional SNPs within the DQA1 gene enabled detection of a recombinant haplotype affecting only the SNPs downstream of the expressed genes. The panel of 13 SNPs was genotyped in 5951 Soay sheep, of which 5349 passed quality control. Using the Soay sheep pedigree, we were able to trace the origin and inheritance of the recombinant SNP haplotype. This SNP-based method has enabled the rapid generation of locus-specific MHC genotypes for large numbers of Soay sheep. This volume of high-quality genotypes in a well-characterized population of free-living sheep will be valuable for investigating the mechanisms maintaining diversity at the MHC.

Heritable variation in telomere length predicts mortality in Soay sheep

Telomere length (TL) is considered an important biomarker of whole-organism health and aging. Across humans and other vertebrates, short telomeres are associated with increased subsequent mortality risk, but the processes responsible for this correlation remain uncertain. A key unanswered question is whether TL–mortality associations arise due to positive effects of genes or early-life environment on both an individual’s average lifetime TL and their longevity, or due to more immediate effects of environmental stressors on within-individual TL loss and increased mortality risk. Addressing this question requires longitudinal TL and life history data across the entire lifetimes of many individuals, which are difficult to obtain for long-lived species like humans. Using longitudinal data and samples collected over nearly two decades, as part of a long-term study of wild Soay sheep, we dissected an observed positive association between TL and subsequent survival using multivariate quantitative genetic models. We found no evidence that telomere attrition was associated with increased mortality risk, suggesting that TL is not an important marker of biological aging or exposure to environmental stress in our study system. Instead, we find that among-individual differences in average TL are associated with increased lifespan. Our analyses suggest that this correlation between an individual’s average TL and lifespan has a genetic basis. This demonstrates that TL has the potential to evolve under natural conditions, and suggests an important role of genetics underlying the widespread observation that short telomeres predict mortality.

Seasonal dynamics and repeatability of gastrointestinal parasite burdens in wild Soay sheep

Seasonality is a ubiquitous feature of wildlife disease ecology, but is determined by a complex interplay of environmental, parasitological and host factors. Gastrointestinal parasites often exhibit strong seasonal dynamics in wild vertebrate populations due to, for example, environmental influences on free-living or vectored life stages, and variation in the physiological and immune status of hosts across their annual cycle. At the same time, wild populations are typically infected with multiple parasites. The seasonal dynamics of co-infecting parasites may differ depending on age and reproductive status, and associations among parasites may be driven by short-term within-individual changes or longer-term interactions that are consistent among hosts. Here, we used faecal samples and egg counts collected repeatedly from individually marked and monitored wild Soay sheep that were part of a long-term study to investigate seasonal dynamics of six gastrointestinal parasite groups (strongyle nematodes, coccidian protozoa, Capillaria, Strongyloides, Nematodirus, and Moniezia). Prevalence and abundance generally tended to be higher spring and summer, and burdens were higher in lambs than adults. Within the highly prevalent strongyle nematode group, we found differences in seasonality of egg counts depending on adult reproductive status. Reproductive ewes had increased counts in spring around the time of birth followed by a drop in abundance in summer, while barren ewes showed little evidence of seasonality. Males showed a sustained rise in egg counts through spring and summer, and sex differences were only strongly apparent in summer. In contrast, in similarly prevalent coccidia we found a peak in faecal oocyst counts in spring but no differences in seasonality among males, barren and pregnant ewes. Using multivariate mixed-effects models, we went on to show that both strongyle and coccidia counts are moderately repeatable across seasons among individuals. We further show that apparent positive correlation between strongyle and coccidia counts was driven by short-term within-individual changes in both parasite burdens rather than long-term among-individual covariation. Overall, our results demonstrate that seasonality varies across demographic and parasite groups and highlight the value of investigating fluctuating susceptibility and exposure over time for understanding epidemiology of a population.

Mutation load decreases with haplotype age in wild Soay sheep

Runs of homozygosity (ROH) are pervasive in diploid genomes and expose the effects of deleterious recessive mutations, but how exactly these regions contribute to variation in fitness remains unclear. Here, we combined empirical analyses and simulations to explore the deleterious effects of ROH with varying genetic map lengths in wild Soay sheep. Using a long-term dataset of 4,592 individuals genotyped at 417K SNPs, we found that inbreeding depression increases with ROH length. A 1% genomic increase in long ROH (>12.5cM) reduced the odds of first-year survival by 12%, compared to only 7% for medium ROH (1.56-12.5cM), while short ROH (<1.56cM) had no effect on survival. We show by forward genetic simulations that this is predicted: compared with shorter ROH, long ROH will have higher densities of deleterious alleles, with larger average effects on fitness and lower population frequencies. Taken together, our results are consistent with the idea that the mutation load decreases in older haplotypes underlying shorter ROH, where purifying selection has had more time to purge deleterious mutations. Finally, our study demonstrates that strong inbreeding depression can persist despite ongoing purging in a historically small population.