Wild Population Of Soay Sheep Research Articles

Investigating pedigree- and SNP-associated components of heritability in a wild population of Soay sheep

Estimates of narrow sense heritability derived from genomic data that contain related individuals may be biased due to the within-family effects such as dominance, epistasis and common environmental factors. However, for many wild populations, removal of related individuals from the data would result in small sample sizes. In 2013, Zaitlen et al. proposed a method to estimate heritability in populations that include close relatives by simultaneously fitting an identity-by-state (IBS) genomic relatedness matrix (GRM) and an identity-by-descent (IBD) GRM. The IBD GRM is identical to the IBS GRM, except relatedness estimates below a specified threshold are set to 0. We applied this method to a sample of 8557 wild Soay sheep from St. Kilda, with genotypic information for 419,281 single nucleotide polymorphisms. We aimed to see how this method would partition heritability into population-level (IBS) and family-associated (IBD) variance for a range of genetic architectures, and so we focused on a mixture of polygenic and monogenic traits. We also implemented a variant of the model in which the IBD GRM was replaced by a GRM constructed from SNPs with low minor allele frequency to examine whether any additive genetic variance is captured by rare alleles. Whilst the inclusion of the IBD GRM did not significantly improve the fit of the model for the monogenic traits, it improved the fit for some of the polygenic traits, suggesting that dominance, epistasis and/or common environment not already captured by the non-genetic random effects fitted in our models may influence these traits.

The impact of SNP density on quantitative genetic analyses of body size traits in a wild population of Soay sheep.

Understanding the genetic architecture underpinning quantitative traits in wild populations is pivotal to understanding the processes behind trait evolution. The 'animal model' is a popular method for estimating quantitative genetic parameters such as heritability and genetic correlation and involves fitting an estimate of relatedness between individuals in the study population. Genotypes at genome-wide markers can be used to estimate relatedness; however, relatedness estimates vary with marker density, potentially affecting results. Increasing density of markers is also expected to increase the power to detect quantitative trait loci (QTL). In order to understand how the density of genetic markers affects the results of quantitative genetic analyses, we estimated heritability and performed genome-wide association studies (GWAS) on five body size traits in an unmanaged population of Soay sheep using two different SNP densities: a dataset of 37,037 genotyped SNPs and an imputed dataset of 417,373 SNPs. Heritability estimates did not differ between the two SNP densities, but the high-density imputed SNP dataset revealed four new SNP-trait associations that were not found with the lower density dataset, as well as confirming all previously-found QTL. We also demonstrated that fitting fixed and random effects in the same step as performing GWAS is a more powerful approach than pre-correcting for covariates in a separate model.

Functionally distinct T-helper cell phenotypes predict resistance to different types of parasites in a wild mammal

The adaptive immune system is critical to an effective response to infection in vertebrates, with T-helper (Th) cells pivotal in orchestrating these responses. In natural populations where co-infections are the norm, different Th responses are likely to play an important role in maintaining host health and fitness, a relationship which remains poorly understood in wild animals. In this study, we characterised variation in functionally distinct Th responses in a wild population of Soay sheep by enumerating cells expressing Th-subset specific transcription factors and quantifying Th-associated cytokines. We tested the prediction that raised Th1 and Th2 responses should predict reduced apicomplexan and helminth parasite burdens, respectively. All measures of Th-associated cytokine production increased with age, while Th17- and regulatory Th-associated cytokine production increased more rapidly with age in males than females. Independent of age, sex, and each other, IL-4 and Gata3 negatively predicted gastro-intestinal nematode faecal egg count, while IFN-γ negatively predicted coccidian faecal oocyst count. Our results provide important support from outside the laboratory that Th1 and Th2 responses predict resistance to different kinds of parasites, and illustrate how harnessing specific reagents and tools from laboratory immunology will illuminate our understanding of host-parasite interactions in the wild.

Vitamin D status is heritable and under environment-dependent selection in the wild.

Vitamin D has a well‐established role in skeletal health and is increasingly linked to chronic disease and mortality in humans and companion animals. Despite the clear significance of vitamin D for health and obvious implications for fitness under natural conditions, no longitudinal study has tested whether the circulating concentration of vitamin D is under natural selection in the wild. Here, we show that concentrations of dietary‐derived vitamin D2 and endogenously produced vitamin D3 metabolites are heritable and largely polygenic in a wild population of Soay sheep (Ovis aries). Vitamin D2 status was positively associated with female adult survival, and vitamin D3 status predicted female fecundity in particular, good environment years when sheep density and competition for resources was low. Our study provides evidence that vitamin D status has the potential to respond to selection, and also provides new insights into how vitamin D metabolism is associated with fitness in the wild.

Genomic prediction in the wild: A case study in Soay sheep.

Genomic prediction, the technique whereby an individual's genetic component of their phenotype is estimated from its genome, has revolutionised animal and plant breeding and medical genetics. However, despite being first introduced nearly two decades ago, it has hardly been adopted by the evolutionary genetics community studying wild organisms. Here, genomic prediction is performed on eight traits in a wild population of Soay sheep. The population has been the focus of a >30year evolutionary ecology study and there is already considerable understanding of the genetic architecture of the focal Mendelian and quantitative traits. We show that the accuracy of genomic prediction is high for all traits, but especially those with loci of large effect segregating. Five different methods are compared, and the two methods that can accommodate zero-effect and large-effect loci in the same model tend to perform best. If the accuracy of genomic prediction is similar in other wild populations, then there is a real opportunity for pedigree-free molecular quantitative genetics research to be enabled in many more wild populations; currently the literature is dominated by studies that have required decades of field data collection to generate sufficiently deep pedigrees. Finally, some of the potential applications of genomic prediction in wild populations are discussed.

Within-trio tests provide little support for post-copulatory selection on major histocompatibility complex haplotypes in a free-living population.

Sexual selection has been proposed as a force that could help maintain the diversity of major histocompatibility complex (MHC) genes in vertebrates. Potential selective mechanisms can be divided into pre-copulatory and post-copulatory, and in both cases, the evidence for occurrence is mixed, especially in natural populations. In this study, we used a large number of parent-offspring trios that were diplotyped for MHC class II genes in a wild population of Soay sheep (Ovis aries) to examine whether there was within-trio post-copulatory selection on MHC class II genes at both the haplotype and diplotype levels. We found there was transmission ratio distortion of one of the eight MHC class II haplotypes (E) which was transmitted less than expected by fathers, and transmission ratio distortion of another haplotype (A) which was transmitted more than expected by chance to male offspring. However, in both cases, these deviations were not significant after correction for multiple tests. In addition, we did not find any evidence of post-copulatory selection at the diplotype level. These results imply that, given known parents, there is no strong post-copulatory selection on MHC class II genes in this population.

The genetic architecture of helminth-specific immune responses in a wild population of Soay sheep (Ovis aries).

Much of our knowledge of the drivers of immune variation, and how these responses vary over time, comes from humans, domesticated livestock or laboratory organisms. While the genetic basis of variation in immune responses have been investigated in these systems, there is a poor understanding of how genetic variation influences immunity in natural, untreated populations living in complex environments. Here, we examine the genetic architecture of variation in immune traits in the Soay sheep of St Kilda, an unmanaged population of sheep infected with strongyle gastrointestinal nematodes. We assayed IgA, IgE and IgG antibodies against the prevalent nematode Teladorsagia circumcincta in the blood plasma of > 3,000 sheep collected over 26 years. Antibody levels were significantly heritable (h2 = 0.21 to 0.57) and highly stable over an individual’s lifespan. IgA levels were strongly associated with a region on chromosome 24 explaining 21.1% and 24.5% of heritable variation in lambs and adults, respectively. This region was adjacent to two candidate loci, Class II Major Histocompatibility Complex Transactivator (CIITA) and C-Type Lectin Domain Containing 16A (CLEC16A). Lamb IgA levels were also associated with the immunoglobulin heavy constant loci (IGH) complex, and adult IgE levels and lamb IgA and IgG levels were associated with the major histocompatibility complex (MHC). This study provides evidence of high heritability of a complex immunological trait under natural conditions and provides the first evidence from a genome-wide study that large effect genes located outside the MHC region exist for immune traits in the wild.

Conserved Genetic Architecture Underlying Individual Recombination Rate Variation in a Wild Population of Soay Sheep (Ovis aries).

Meiotic recombination breaks down linkage disequilibrium (LD) and forms new haplotypes, meaning that it is an important driver of diversity in eukaryotic genomes. Understanding the causes of variation in recombination rate is important in interpreting and predicting evolutionary phenomena and in understanding the potential of a population to respond to selection. However, despite attention in model systems, there remains little data on how recombination rate varies at the individual level in natural populations. Here we used extensive pedigree and high-density SNP information in a wild population of Soay sheep (Ovis aries) to investigate the genetic architecture of individual autosomal recombination rates. Individual rates were high relative to other mammal systems and were higher in males than in females (autosomal map lengths of 3748 and 2860 cM, respectively). The heritability of autosomal recombination rate was low but significant in both sexes (h2 = 0.16 and 0.12 in females and males, respectively). In females, 46.7% of the heritable variation was explained by a subtelomeric region on chromosome 6; a genome-wide association study showed the strongest associations at locus RNF212, with further associations observed at a nearby ∼374-kb region of complete LD containing three additional candidate loci, CPLX1, GAK, and PCGF3. A second region on chromosome 7 containing REC8 and RNF212B explained 26.2% of the heritable variation in recombination rate in both sexes. Comparative analyses with 40 other sheep breeds showed that haplotypes associated with recombination rates are both old and globally distributed. Both regions have been implicated in rate variation in mice, cattle, and humans, suggesting a common genetic architecture of recombination rate variation in mammals.

Plasma markers of oxidative stress are uncorrelated in a wild mammal.

Oxidative stress, which results from an imbalance between the production of potentially damaging reactive oxygen species versus antioxidant defenses and repair mechanisms, has been proposed as an important mediator of life-history trade-offs. A plethora of biomarkers associated with oxidative stress exist, but few ecological studies have examined the relationships among different markers in organisms experiencing natural conditions or tested whether those relationships are stable across different environments and demographic groups. It is therefore not clear to what extent studies of different markers can be compared, or whether studies that focus on a single marker can draw general conclusions regarding oxidative stress. We measured widely used markers of oxidative damage (protein carbonyls and malondialdehyde) and antioxidant defense (superoxide dismutase and total antioxidant capacity) from 706 plasma samples collected over a 4-year period in a wild population of Soay sheep on St Kilda. We quantified the correlation structure among these four markers across the entire sample set and also within separate years, age groups (lambs and adults), and sexes. We found some moderately strong correlations between some pairs of markers when data from all 4years were pooled. However, these correlations were caused by considerable among-year variation in mean marker values; correlation coefficients were small and not significantly different from zero after accounting for among-year variation. Furthermore, within each year, age, and sex subgroup, the pairwise correlation coefficients among the four markers were weak, nonsignificant, and distributed around zero. In addition, principal component analysis confirmed that the four markers represented four independent axes of variation. Our results suggest that plasma markers of oxidative stress may vary dramatically among years, presumably due to environmental conditions, and that this variation can induce population-level correlations among markers even in the absence of any correlations within contemporaneous subgroups. The absence of any consistent correlations within years or demographic subgroups implies that care must be taken when generalizing from observed relationships with oxidative stress markers, as each marker may reflect different and potentially uncoupled biochemical processes.

Asynchrony of senescence among phenotypic traits in a wild mammal population

The degree to which changes in lifespan are coupled to changes in senescence in different physiological systems and phenotypic traits is a central question in biogerontology. It is underpinned by deeper biological questions about whether or not senescence is a synchronised process, or whether levels of synchrony depend on species or environmental context. Understanding how natural selection shapes patterns of synchrony in senescence across physiological systems and phenotypic traits demands the longitudinal study of many phenotypes under natural conditions. Here, we examine the patterns of age-related variation in late adulthood in a wild population of Soay sheep (Ovis aries) that have been the subject of individual-based monitoring for thirty years. We examined twenty different phenotypic traits in both males and females, encompassing vital rates (survival and fecundity), maternal reproductive performance (offspring birth weight, birth date and survival), male rutting behaviour, home range measures, parasite burdens, and body mass. We initially quantified age-related variation in each trait having controlled for annual variation in the environment, among-individual variation and selective disappearance effects. We then standardised our age-specific trait means and tested whether age trajectories could be meaningfully grouped according to sex or the type of trait. Whilst most traits showed age-related declines in later life, we found striking levels of asynchrony both within and between the sexes. Of particular note, female fecundity and reproductive performance declined with age, but male annual reproductive success did not. We also discovered that whilst home range size and quality decline with age in females, home range size increases with age in males. Our findings highlight the complexity of phenotypic ageing under natural conditions and, along with emerging data from other wild populations and laboratory models, suggest that the long-standing hypothesis within evolutionary biology that fitness-related traits should senesce in a synchronous manner is seriously flawed.

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.

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.

The Impact of Environmental Heterogeneity on Genetic Architecture in a Wild Population of Soay Sheep

This work demonstrates that environmental conditions experienced by individuals can shape their development and affect the stability of genetic associations. The implication of this observation is that the environmental response may influence the evolution of traits in the wild. Here, we examined how the genetic architecture of a suite of sexually dimorphic traits changed as a function of environmental conditions in an unmanaged population of Soay sheep (Ovis aries) on the island of Hirta, St. Kilda, northwest Scotland. We examined the stability of phenotypic, genetic, and environmental (residual) covariance in males during the first year of life between horn length, body weight, and parasite load in environments of different quality. We then examined the same covariance structures across environments within and between the adult sexes. We found significant genotype-by-environment interactions for lamb male body weight and parasite load, leading to a change in the genetic correlation among environments. Horn length was genetically correlated with body weight in males but not females and the genetic correlation among traits within and between the sexes was dependent upon the environmental conditions experienced during adulthood. Genetic correlations were smaller in more favorable environmental conditions, suggesting that in good environments, loci are expressed that have sex-specific effects. The reduction in genetic correlation between the sexes may allow independent evolutionary trajectories for each sex. This study demonstrates that the genetic architecture of traits is not stable under temporally varying environments and highlights the fact that evolutionary processes may depend largely upon ecological conditions.

Environmental Heterogeneity Generates Fluctuating Selection on a Secondary Sexual Trait

In any population in which resources are limiting, the allocation of resources toward increased reproductive success may generate costs to survival [1-8]. The relationship between a sexually selected trait and fitness will therefore represent a balance between its relative associations with fecundity versus viability [3, 6, 7]. Because the risk of mortality in a population is likely to be heavily determined by ecological conditions, survival costs may vary as a function of the prevailing environment [7]. As a result, for populations experiencing heterogeneous ecological conditions, there may not be a single optimal level of allocation toward reproduction versus survival [9]. Here, we show that early viability and fecundity selection act in opposing directions on a secondary sexual trait and that their relative magnitude depends upon ecological conditions, generating fluctuating selection. In a wild population of Soay sheep (Ovis aries), phenotypic and genetic associations between male horn growth and lifetime reproductive success were positive under good environmental conditions (because of increased breeding success) and negative under poor environmental conditions (because of reduced survival). In an unpredictable environment, high allocation to early horn growth is a gamble that will only pay off if ensuing conditions are favorable. Such fluctuating selection may play an important role in preventing the erosion of genetic variance in secondary sexual traits.