Using genomic prediction to detect microevolutionary change of a quantitative trait.

Detecting microevolutionary responses to natural selection by observing temporal changes in individual breeding values is challenging. The collection of suitable datasets can take many years and disentangling the contributions of the environment and genetics to phenotypic change is not trivial. Furthermore, pedigree-based methods of obtaining individual breeding values have known biases. Here, we apply a genomic prediction approach to estimate breeding values of adult weight in a 35-year dataset of Soay sheep (Ovis aries). Comparisons are made with a traditional pedigree-based approach. During the study period, adult body weight decreased, but the underlying genetic component of body weight increased, at a rate that is unlikely to be attributable to genetic drift. Thus cryptic microevolution of greater adult body weight has probably occurred. Genomic and pedigree-based approaches gave largely consistent results. Thus, using genomic prediction to study microevolution in wild populations can remove the requirement for pedigree data, potentially opening up new study systems for similar research.

Genetic architecture and lifetime dynamics of inbreeding depression in a wild mammal.

Inbreeding depression is ubiquitous, but we still know little about its genetic architecture and precise effects in wild populations. Here, we combine long-term life-history data with 417 K imputed SNP genotypes for 5952 wild Soay sheep to explore inbreeding depression on a key fitness component, annual survival. Inbreeding manifests in long runs of homozygosity (ROH), which make up nearly half of the genome in the most inbred individuals. The ROH landscape varies widely across the genome, with islands where up to 87% and deserts where only 4% of individuals have ROH. The fitness consequences of inbreeding are severe; a 10% increase in individual inbreeding FROH is associated with a 60% reduction in the odds of survival in lambs, though inbreeding depression decreases with age. Finally, a genome-wide association scan on ROH shows that many loci with small effects and five loci with larger effects contribute to inbreeding depression in survival.

Senescence in immunity against helminth parasites predicts adult mortality in a wild mammal.

Our understanding of the deterioration in immune function in old age-immunosenescence-derives principally from studies of modern human populations and laboratory animals. The generality and significance of this process for systems experiencing complex, natural infections and environmental challenges are unknown. Here, we show that late-life declines in an important immune marker of resistance to helminth parasites in wild Soay sheep predict overwinter mortality. We found senescence in circulating antibody levels against a highly prevalent nematode worm, which was associated with reduced adult survival probability, independent of changes in body weight. These findings establish a role for immunosenescence in the ecology and evolution of natural populations.

No evidence for parental age effects on offspring leukocyte telomere length in free-living Soay sheep.

In humans, the effect of paternal age at conception (PAC) on offspring leukocyte telomere length (LTL) is well established, with older fathers thought to pass on longer telomeres to their offspring in their sperm. Few studies have looked for PAC effects in other species, but it has been hypothesised that the effect will be exacerbated in polygamous species with higher levels of sperm competition and production. We test for maternal (MAC) and paternal age at conception effects on offspring LTL in Soay sheep, a primitive breed experiencing strong sperm competition. We use qPCR to measure relative telomere length in 389 blood samples (n = 318 individuals) collected from an unmanaged population of sheep on St Kilda, where individual age and parentage are known. We find no evidence that either MAC or PAC are associated with LTL in offspring across the age range, or when considering only young lambs (n = 164). This is the first study to test for parental age effects on offspring LTL in a wild mammal population, and the results contrast with the findings of numerous human studies that find a PAC effect, as well as predictions of a stronger PAC effect in polygamous species.

Robust extraction of quantitative structural information from high-variance histological images of livers from necropsied Soay sheep.

Quantitative information is essential to the empirical analysis of biological systems. In many such systems, spatial relations between anatomical structures is of interest, making imaging a valuable data acquisition tool. However, image data can be difficult to analyse quantitatively. Many image processing algorithms are highly sensitive to variations in the image, limiting their current application to fields where sample and image quality may be very high. Here, we develop robust image processing algorithms for extracting structural information from a dataset of high-variance histological images of inflamed liver tissue obtained during necropsies of wild Soay sheep. We demonstrate that features of the data can be measured in a fully automated manner, providing quantitative information which can be readily used in statistical analysis. We show that these methods provide measures that correlate well with a manual, expert operator-led analysis of the same images, that they provide advantages in terms of sampling a wider range of information and that information can be extracted far more quickly than in manual analysis.

Physiological proteins in resource-limited herbivores experiencing a population die-off.

Nutrient availability is predicted to interact with herbivore population densities. Competition for low quality food at high density may reduce summer food intake, and in turn winter survival. Conversely, low population density may favor physiological recovery through better access to better quality spring forage. Here, we take advantage of the long-term study of the Soay sheep population of St. Kilda (Scotland) to measure plasma protein markers and immunity in two consecutive summers with contrasting population densities. We show that, following a winter die-off resulting in a shift to low population density, albumin and total proteins increased, but only in adult sheep. The effect was not solely attributable to selective disappearance of malnourished sheep. Similarly, the concentration of antibodies was higher following the die-off, potentially indicating recovery of immune function. Overall, our results are consistent with the physiological recovery of surviving individuals after a harsh winter.

Inbreeding depression by environment interactions in a free-living mammal population.

Experimental studies often find that inbreeding depression is more severe in harsh environments, but the few studies of in situ wild populations available to date rarely find strong support for this effect. We investigated evidence for inbreeding depression by environment interactions in nine traits in the individually monitored Soay sheep population of St Kilda, using genomic inbreeding coefficients based on 37 037 single-nucleotide polymorphism loci, and population density as an axis of environmental variation. All traits showed variation with population density and all traits showed some evidence for depression because of either an individual's own inbreeding or maternal inbreeding. However, only six traits showed evidence for an interaction in the expected direction, and only two interactions were statistically significant. We identify three possible reasons why wild population studies may generally fail to find strong support for interactions between inbreeding depression and environmental variation compared with experimental studies. First, for species with biparental inbreeding only, the amount of observed inbreeding in natural populations is generally low compared with that used in experimental studies. Second, it is possible that experimental studies sometimes actually impose higher levels of stress than organisms experience in the wild. Third, some purging of the deleterious recessive alleles that underpin interaction effects may occur in the wild.

Accounting for female space sharing in St. Kilda Soay sheep (Ovis aries) results in little change in heritability estimates.

When estimating heritability in free-living populations, it is common practice to account for common environment effects, because of their potential to generate phenotypic covariance among relatives thereby biasing heritability estimates. In quantitative genetic studies of natural populations, however, philopatry, which results in relatives being clustered in space, is rarely accounted for. The two studies that have been carried out so far suggest absolute declines in heritability estimates of up to 43% when accounting for space sharing by relatives. However, due to methodological limitations these estimates may not be representative. We used data from the St. Kilda Soay sheep population to estimate heritabilities with and without accounting for space sharing for five traits for which there is evidence for additive genetic variance (birthweight, birth date, lamb August weight, and female post-mortem jaw and metacarpal length). We accounted for space sharing by related females by separately incorporating spatial autocorrelation, and a home range similarity matrix. Although these terms accounted for up to 18% of the variance in these traits, heritability estimates were only reduced by up to 7%. Our results suggest that the bias caused by not accounting for space sharing may be lower than previously thought. This suggests that philopatry does not inevitably lead to a large bias if space sharing by relatives is not accounted for. We hope our work stimulates researchers to model shared space when relatives in their study population share space, as doing so will enable us to better understand when bias may be of particular concern.

Exposure to viral and bacterial pathogens among Soay sheep (Ovis aries) of the St Kilda archipelago.

We assessed evidence of exposure to viruses and bacteria in an unmanaged and long-isolated population of Soay sheep (Ovis aries) inhabiting Hirta, in the St Kilda archipelago, 65 km west of Benbecula in the Outer Hebrides of Scotland. The sheep harbour many metazoan and protozoan parasites but their exposure to viral and bacterial pathogens is unknown. We tested for herpes viral DNA in leucocytes and found that 21 of 42 tested sheep were infected with ovine herpesvirus 2 (OHV-2). We also tested 750 plasma samples collected between 1997 and 2010 for evidence of exposure to seven other viral and bacterial agents common in domestic Scottish sheep. We found evidence of exposure to Leptospira spp., with overall seroprevalence of 6·5%. However, serological evidence indicated that the population had not been exposed to border disease, parainfluenza, maedi-visna, or orf viruses, nor to Chlamydia abortus. Some sheep tested positive for antibodies against Mycobacterium avium subsp. paratuberculosis (MAP) but, in the absence of retrospective faecal samples, the presence of this infection could not be confirmed. The roles of importation, the pathogen-host interaction, nematode co-infection and local transmission warrant future investigation, to elucidate the transmission ecology and fitness effects of the few viral and bacterial pathogens on Hirta.

Detecting genes for variation in parasite burden and immunological traits in a wild population: testing the candidate gene approach.

Identifying the genes underlying phenotypic variation in natural populations can provide novel insight into the evolutionary process. The candidate gene approach has been applied to studies of a number of traits in various species, in an attempt to elucidate their genetic basis. Here, we test the application of the candidate gene approach to identify the loci involved in variation in gastrointestinal parasite burden, a complex trait likely to be controlled by many loci, in a wild population of Soay sheep. A comprehensive literature review, Gene Ontology databases, and comparative genomics resources between cattle and sheep were used to generate a list of candidate genes. In a pilot study, these candidates, along with 50 random genes, were then sequenced in two pools of Soay sheep; one with low gastrointestinal nematode burden and the other high, using a NimbleGen sequence capture experiment. Further candidates were identified from single nucleotide polymorphisms (SNPs) that were highly differentiated between high- and low-resistance sheep breeds. A panel of 192 candidate and control SNPs were then typed in 960 individual Soay sheep to examine whether they individually explained variation in parasite burden, as measured as faecal egg count, as well as two immune measures (Teladorsagia circumcincta-specific antibodies and antinuclear antibodies). The cumulative effect of the candidate and control SNPs were estimated by fitting genetic relationship matrices (GRMs) as random effects in animal models of the three traits. No more significant SNPs were identified in the pilot sequencing experiment and association study than expected by chance. Furthermore, no significant difference was found between the proportions of candidate or control SNPs that were found to be significantly associated with parasite burden/immune measures. No significant effect of the candidate or control gene GRMs was found. There is thus little support for the candidate gene approach to the identification of loci explaining variation in parasitological and immunological traits in this population. However, a number of SNPs explained significant variation in multiple traits and significant correlations were found between the proportions of variance explained by individual SNPs across multiple traits. The significant SNPs identified in this study may still, therefore, merit further investigation.

Selection and microevolution of coat pattern are cryptic in a wild population of sheep.

Understanding the maintenance of genetic variation in natural populations is a core aim of evolutionary genetics. Insight can be gained by quantifying selection at the level of the genotype, as opposed to the phenotype. Here, we show that in a natural population of Soay sheep which is polymorphic for coat pattern, recessive genetic variants at the causal gene, agouti signalling protein (ASIP) are associated with reduced lifetime fitness. This was due primarily to a reduction in juvenile survival of uniformly coloured (self-type) sheep, which are homozygous recessive, and occurs despite significantly higher reproductive success in surviving self-type adults. Consistent with their relatively low fitness, we show that the frequency of self-type individuals has declined from 1985 to 2008. Remarkably though, the frequency of the underlying self-allele has increased, because the frequency of heterozygous individuals (who harbour the majority of all self-alleles) has increased. Indeed, the ratio of observed/expected heterozygous individuals has increased during the study, such that there is now a significant excess of heterozygotyes. By employing gene-dropping simulations, we show that microevolutionary trends in the frequency and excess of ASIP heterozygotes are too pronounced to be caused by genetic drift. Studying this polymorphism at the level of phenotype rather than underlying genotype would have failed to detect cryptic fitness differences. We would also have been unable to rule out genetic drift as an evolutionary force driving genetic change. This highlights the importance of resolving the underlying genetic basis of phenotypic variation in explaining evolutionary dynamics.

Natural selection on a measure of parasite resistance varies across ages and environmental conditions in a wild mammal.

Parasites detrimentally affect host fitness, leading to expectations of positive selection on host parasite resistance. However, as immunity is costly, host fitness may be maximized at low, but nonzero, parasite infection intensities. These hypotheses are rarely tested on natural variation in free-living populations. We investigated selection on a measure of host parasite resistance in a naturally regulated Soay sheep population using a longitudinal data set and found negative correlations between parasite infection intensity and annual fitness in lambs, male yearlings and adult females. However, having accounted for confounding effects of body weight, the effect was only significant in lambs. Associations between fitness and parasite resistance were environment-dependent, being strong during low-mortality winters, but negligible during harsher high-mortality winters. There was no evidence for stabilizing selection. Our findings reveal processes that may shape variation in parasite resistance in natural populations and illustrate the importance of accounting for correlated traits in selection analysis.

Mapping quantitative trait loci in a wild population using linkage and linkage disequilibrium analyses.

Historical information can be used, in addition to pedigree, traits and genotypes, to map quantitative trait locus (QTL) in general populations via maximum likelihood estimation of variance components. This analysis is known as linkage disequilibrium (LD) and linkage mapping, because it exploits both linkage in families and LD at the population level. The search for QTL in the wild population of Soay sheep on St. Kilda is a proof of principle. We analysed the data from a previous study and confirmed some of the QTLs reported. The most striking result was the confirmation of a QTL affecting birth weight that had been reported using association tests but not when using linkage-based analyses.

Trypanosoma melophagium from the sheep ked Melophagus ovinus on the island of St Kilda.

SUMMARYThe sheep ked has been largely eradicated in the UK but persists in the feral Soay sheep of St Kilda in the Outer Hebrides. Sheep keds transmit Trypanosoma melophagium, but parasitaemias are typically cryptic and this trypanosome has not been recorded in the St Kilda sheep. Trypanosomes were detected by PCR in preserved keds and were also found in gut smears from live keds; one infected gut was used to establish the trypanosome in vitro. Examination of the morphology of bloodstream forms from culture confirmed its identity as T. melophagium. Most keds were found to harbour the trypanosome, particularly those collected from lambs. DNA was extracted from preserved keds and from trypanosomes grown in vitro. Sequence analysis of the small subunit ribosomal RNA (SSU rRNA) gene and the spliced leader transcript showed the T. melophagium sequences to be very similar to those from T. theileri. A partial sequence of the ked SSU rRNA gene was also obtained. The close genetic relationship of T. melophagium and T. theileri suggests that T. melophagium represents a lineage of T. theileri that adapted to transmission by sheep keds and hence became a specific parasite of sheep.

Maternal effects and early-life performance are associated with parasite resistance across life in free-living Soay sheep.

Maternal effects occur when the maternal phenotype influences that of the offspring in addition to the effects of maternal genes, and may have a considerable influence on offspring parasite resistance. These effects, and the effects of early levels of reproduction and parasite resistance, may persist into later life and even influence ageing rates. Here we analyse a 20-year longitudinal data set collected on a free-living population of Soay sheep, to investigate the associations between a suite of maternal phenotypic traits and early-life performance on measures of parasite resistance across life. Our results show that maternal effects are important in determining offspring parasite resistance, since lambs born as twins and those born to the youngest and oldest mothers show higher parasite burdens. We show that the association between parasite resistance and natal litter size persists into adulthood. We also show that age-specific changes in parasite resistance in males are associated with natal litter size, and that age-specific changes in females are influenced by early-life levels of reproduction and parasite infection. These results add to the growing evidence that conditions experienced by individuals during development can have a profound influence on immediate and late-life performance and may even influence ageing.

The genetic basis of recessive self-colour pattern in a wild sheep population.

Bridging the genotype-phenotype gap for traits of ecological and evolutionary importance in natural populations can provide a novel insight into the origin and maintenance of variation. Here, we identify the gene and putative causal mutations underlying a recessive colour pattern phenotype ('self' or uniform colour) in a wild population of primitive Soay sheep. We targeted the agouti signalling protein (ASIP) gene, a positional candidate based on previous study that mapped the Coat pattern locus to a presumptive region on chromosome 13. We found evidence for three recessive mutations, including two functional changes in the coding sequence and a putative third cis-regulatory mutation that inactivates the promoter. These mutations define up to five haplotypes in Soays, which collectively explained the coat pattern in all but one member of a complex multi-generational pedigree containing 621 genotyped individuals. The functional mutations are in strong linkage disequilibrium in the study population, and are identical to those known to underlie the self phenotype in domestic sheep. This is indicative of a recent (and simultaneous) origin in Soay sheep, possibly as a consequence of past interbreeding with modern domestic breeds. This is only the second study in which ASIP has been linked to variation in pigmentation in a natural population. Knowledge of the genetic basis of self-colour pattern in Soay sheep, and the recognition that several mutations are segregating in the population, will aid future studies investigating the role of selection in the maintenance of the polymorphism.

Sex ratio and morphological polymorphism in an isolated, endemic Teladorsagia circumcincta population.

Teladorsagia circumcincta is a polygamous nematode that exhibits morphological polymorphism. Sex ratio is typically female biased and the male nematodes occur in association with the genetically similar, minor morphotypes Teladorsagia davtiani and Teladorsagia trifurcata. In experimental infections, sex ratio (proportion male) and the proportion of minor male morphs observed have been shown to be influenced by both host and nematode-related factors. As similar investigations from natural systems are rare, this study examined whether sex ratio and minor male morph frequency were associated with host age and sex and nematode infra-population size in the isolated Soay sheep population on St Kilda. Count data for Teladorsagia were analysed for sheep of all age classes and both sexes from the winters of three consecutive population crashes (1999, 2002 and 2005). Generally, the intensity of Teladorsagia nematodes increased with host age until the age of 2 years before decreasing. In 2005, abundance of nematodes was generally higher than in the previous crashes, nematode sex ratio was negatively associated with host age and tended to be negatively associated with nematode intensity. Within the male nematode subpopulation, T. circumcincta always predominated, followed by T. davtiani and then T. trifurcata, with little variation in the relative proportions between hosts. The presence of each minor morph was primarily associated with the intensity of male T. circumcincta and, in those hosts where all three male morphs were detected, intensity of each minor morph was most associated with intensity of Teladorsagia females. Therefore, in a year when the nematode was generally more abundant, sex ratio appeared to be influenced by both host and nematode-related factors, whereas in all years examined, the frequency of morphological polymorphism was primarily density dependent.

Trading offspring size for number in a variable environment: selection on reproductive investment in female Soay sheep.

1. Given energetic constraints, female reproductive strategy is expected to be shaped by a trade-off between offspring size and offspring number, the optimal resolution of which may vary with environmental conditions. 2. We tested the hypothesis that selection will favour the production of larger litters, even at some expense to offspring size, under good conditions (and vice versa in harsh environments) using data from a long-term study of an unmanaged population of Soay sheep on the islands of St Kilda, NW Scotland. 3. Both litter size (which is either 1 or 2) and offspring birth weight are under positive selection through female annual fitness, but the strength of selection varies systematically with environmental conditions. Age effects were also detected, with selection weakening as female age increases. 4. Consistent with theoretical predictions, the selective trade-off between litter and offspring size is shifted in favour of the latter under poor environmental conditions. Thus, direct selection on offspring birth weight increases under harsh environmental conditions, particularly for females producing twins. 5. However, singletons are only favoured when environmental conditions are very poor, and offspring weight is constrained to be low. Thus, the current low incidence of twinning (13.5% of litters produced since 1985) appears to be suboptimal with respect to the empirically estimated selection regime. Thus, litter size, a trait known to be heritable, may be expected to evolve upwards in this population. 6. Our study highlights the necessity of incorporating environmental heterogeneity and age structure into analyses of natural selection, and suggests that the common assumption of optimality used in models of life history may sometimes be problematic.

Re-establishment of nematode infra-community and host survivorship in wild Soay sheep following anthelmintic treatment.

Compared to mainland Britain, where there has been decades of anthelmintic use, the natural host-parasite relationship of the wild Soay sheep on the remote archipelago of St. Kilda has remained undisturbed. Small-scale anthelmintic bolus experiments on the island have previously shown that the removal of gastrointestinal nematodes can improve over-winter survival of young and male sheep in high host density years. This study, in which two-year-old sheep were treated, is the first to examine patterns of re-establishment of different nematode species in sheep following treatment and also investigates which species are likely to affect host survival. The experiment showed that, although all sheep were equally likely to die, host sex and weight influenced temporal survivorship in that females and heavier sheep tended to survive longer. Examination of the nematodes that had re-infected males six months after administration of an anthelmintic bolus showed that, compared to controls, the diversity of species was lower. Of the nine nematode species, there were fewer Trichostrongylus axei and Trichostrongylus vitrinus nematodes in treated sheep. In control sheep there were more T. axei and T. vitrinus nematodes in males than females. In addition there was an association between host over winter weight loss and the intensity of T. vitrinus. The fact that this species had higher numbers in males than females and was associated with over-winter weight loss implies it could have an important role in host mortality in high-density years.

Estimating the functional form for the density dependence from life history data.

Two contrasting approaches to the analysis of population dynamics are currently popular: demographic approaches where the associations between demographic rates and statistics summarizing the population dynamics are identified; and time series approaches where the associations between population dynamics, population density, and environmental covariates are investigated. In this paper, we develop an approach to combine these methods and apply it to detailed data from Soay sheep (Ovis aries). We examine how density dependence and climate contribute to fluctuations in population size via age- and sex-specific demographic rates, and how fluctuations in demographic structure influence population dynamics. Density dependence contributes most, followed by climatic variation, age structure fluctuations and interactions between density and climate. We then simplify the density-dependent, stochastic, age-structured demographic model and derive a new phenomenological time series which captures the dynamics better than previously selected functions. The simple method we develop has potential to provide substantial insight into the relative contributions of population and individual-level processes to the dynamics of populations in stochastic environments.