PARASITE-MEDIATED SELECTION AGAINST INBRED SOAY SHEEP IN A FREE-LIVING ISLAND POPULATON

Parasites are thought to provide a selective force capable of promoting genetic variation in natural populations. One rarely considered pathway for this action is via parasite-mediated selection against inbreeding. If parasites impose a fitness cost on their host and the offspring of close relatives have greater susceptibility to parasites due to the increased homozygosity that results from inbreeding, then parasite-mediated mortality may select against inbred individuals. This hypothesis has not yet been tested within a natural vertebrate population. Here we show that relatively inbred Soay sheep (Ovis aries), as assessed by microsatellite heterozygosity, are more susceptible to parasitism by gastrointestinal nematodes, with interactions indicating greatest susceptibility among adult sheep at high population density. During periods of high overwinter mortality on the island of Hirta, St. Kilda, Scotland, highly parasitised individuals were less likely to survive. More inbred individuals were also less likely to survive, which is due to their increased susceptibility to parasitism, because survival was random with respect to inbreeding among sheep that were experimentally cleared of their gastrointestinal parasite burden by anthelminthic treatment. As a consequence of this selection, average microsatellite heterozygosity increases with age in St. Kildan Soay sheep. We suggest that parasite-mediated selection acts to maintain genetic variation in this small island population by removing less heterozygous individuals.     

Survival costs of reproduction are mediated by parasite infection in wild Soay sheep

A trade‐off between current and future fitness potentially explains variation in life‐history strategies. A proposed mechanism behind this is parasite‐mediated reproductive costs: individuals that allocate more resources to reproduction have fewer to allocate to defence against parasites, reducing future fitness. We examined how reproduction influenced faecal egg counts (FEC) of strongyle nematodes using data collected between 1989 and 2008 from a wild population of Soay sheep in the St. Kilda archipelago, Scotland (741 individuals). Increased reproduction was associated with increased FEC during the lambing season: females that gave birth, and particularly those that weaned a lamb, had higher FEC than females that failed to reproduce. Structural equation modelling revealed future reproductive costs: a positive effect of reproduction on spring FEC and a negative effect on summer body weight were negatively associated with overwinter survival. Overall, we provide evidence that parasite resistance and body weight are important mediators of survival costs of reproduction.     

POSITIVE GENETIC CORRELATION BETWEEN PARASITE RESISTANCE AND BODY SIZE IN A FREE-LIVING UNGULATE POPULATION

Abstract Parasite resistance and body size are subject to directional natural selection in a population of feral Soay sheep (Ovis aries) on the island of St. Kilda, Scotland. Classical evolutionary theory predicts that directional selection should erode additive genetic variation and favor the maintenance of alleles that have negative pleiotropic effects on other traits associated with fitness. Contrary to these predictions, in this study we show that there is considerable additive genetic variation for both parasite resistance, measured as fecal egg count (FEC), and body size, measured as weight and hindleg length, and that there are positive genetic correlations between parasite resistance and body size in both sexes. Body size traits had higher heritabilities than parasite resistance. This was not due to low levels of additive genetic variation for parasite resistance, but was a consequence of high levels of residual variance in FEC. Measured as coefficients of variation, levels of additive genetic variation for FEC were actually higher than for weight or hindleg length. High levels of additive genetic variation for parasite resistance may be maintained by a number of mechanisms including high mutational input, balancing selection, antagonistic pleiotropy, and host‐parasite coevolution. The positive genetic correlation between parasite resistance and body size, a trait also subject to sexual selection in males, suggests that parasite resistance and growth are not traded off in Soay sheep, but rather that genetically resistant individuals also experience superior growth.     

Quantitative trait loci (QTL) mapping of resistance to strongyles and coccidia in the free-living Soay sheep (Ovis aries)

A genome-wide scan was performed to detect quantitative trait loci (QTL) for resistance to gastrointestinal parasites and ectoparasitic keds segregating in the free-living Soay sheep population on St. Kilda (UK). The mapping panel consisted of a single pedigree of 882 individuals of which 588 were genotyped. The Soay linkage map used for the scans comprised 251 markers covering the whole genome at average spacing of 15cM. The traits here investigated were the strongyle faecal egg count (FEC), the coccidia faecal oocyst count (FOC) and a count of keds (Melophagus ovinus). QTL mapping was performed by means of variance component analysis so that the genetic parameters of the study traits were also estimated and compared with previous studies in Soay and domestic sheep. Strongyle FEC and coccidia FOC showed moderate heritability (h(2)=0.26 and 0.22, respectively) in lambs but low heritability in adults (h(2)<0.10). Ked count appeared to have very low h(2) in both lambs and adults. Genome scans were performed for the traits with moderate heritability and two genomic regions reached the level of suggestive linkage for coccidia FOC in lambs (logarithm of the odds=2.68 and 2.21 on chromosomes 3 and X, respectively). We believe this is the first study to report a QTL search for parasite resistance in a free-living animal population and therefore may represent a useful reference for similar studies aimed at understanding the genetics of host-parasite co-evolution in the wild.     

Metazoan-protozoan parasite co-infections and host body weight in St Kilda Soay sheep

For hundreds of years, the unmanaged Soay sheep population on St Kilda has survived despite enduring presumably deleterious co-infections of helminth, protozoan and arthropod parasites and intermittent periods of starvation. Important parasite taxa in young Soay sheep are strongyles (Trichostrongylus axei, Trichostrongylus vitrinus and Teladorsagia circumcincta), coccidia (11 Eimeria species) and keds (Melophagus ovinus) and in older animals, Teladorsagia circumcincta. In this research, associations between the intensity of different parasite taxa were investigated. Secondly, the intensities of different parasite taxa were tested for associations with variation in host weight, which is itself a determinant of over-winter survival in the host population. In lambs, the intensity of strongyle eggs was positively correlated with that of Nematodirus spp. eggs, while in yearlings and adults strongyle eggs and coccidia oocysts were positively correlated. In lambs and yearlings, of the parasite taxa tested, only strongyle eggs were significantly and negatively associated with host weight. However, in adult hosts, strongyles and coccidia were independently and negatively associated with host weight. These results are consistent with the idea that strongyles and coccidia are exerting independent selection on Soay sheep.     

Host isolation and patterns of genetic variability in three populations of Teladorsagia from sheep

We have used a mitochondrial marker to explore the population genetics of an economically important parasite of sheep, Teladorsagia. We examined diversity within and between parasites from three very different host populations, as well as within and between individual hosts. One of our study populations, the Soay sheep on Hirta, St Kilda, is unusually isolated with no sheep having been introduced to the island since 1932. Worm haplotypes from Hirta were compared with those from two other host populations. Remarkably, despite its historical isolation the Hirta population shows similar levels of within-population diversity to the other study populations. No divergence between the three Teladorsagia populations was found, consistent with gene flow between the populations. The high diversity within Teladorsagia populations provides compelling evidence that this variability is a general feature of parasitic nematode populations. Such diversity may be caused by high effective population size, coupled with an increased mutation rate for mtDNA, which has important implications for the spread of anthelmintic resistance in nematode populations.     

Genetic parameters for resistance to trichostrongylid infection in dairy sheep

In sheep, the traditional chemical control of gastrointestinal nematode (GIN) parasites with anthelmintics has led to the widespread development of anthelmintic resistance. The selection of sheep with enhanced resistance to GIN parasites has been suggested as an alternative strategy to develop sustainable control of parasite infections. Most of the estimations of the genetic parameters for sheep resistance to GIN parasites have been obtained from young animals belonging to meat- and/or wool-specialised breeds. We present here the estimated genetic parameters for four parasite resistance traits studied in a commercial population of adult Spanish Churra dairy ewes. These involved two faecal egg counts (FECs) (LFEC0 and LFEC1) and two serum indicator traits, the anti-Teladorsagia circumcincta fourth stage larvae IgA (IgA) and the pepsinogen (Peps) levels. In addition, this study has allowed us to identify the environmental factors influencing parasite resistance in naturally infected Spanish Churra sheep and to quantify the genetic component of this complex phenotype. The heritabilities estimated for the two FECs analysed (0.12 for LFEC0 and 0.09 for LFEC1) were lower than those obtained for the examined serum indicators (0.19 for IgA and 0.21 for Peps). The genetic correlations between the traits ranged from 0.43 (Peps-IgA) to 0.82 (LFEC0-LFEC1) and were higher than their phenotypic counterparts, which ranged between 0.07 and 0.10. The heritabilities estimated for the studied traits were lower than previously reported in lambs. This may be due to the differences in the immune mechanisms controlling the infection in young (antibody reactions) and adult (hypersensitivity reactions) animals/sheep. In summary, this study demonstrates the presence of heritable variation in parasite resistance indicator traits in the Churra population studied, which suggests that genetic improvement is feasible for this complex trait in this population. However, further studies in which the experimental variables are controlled as much as possible are needed to identify the best trait that could be measured routinely in adult sheep as an indicator of parasite resistance.     

Grazing decisions of Soay sheep,Ovis aries,on St Kilda: a consequence of parasite distribution?

Relationships between the height of grass swards and herbage intake have been established for a number of large herbivores and sward height plays a major role in determining the selection of herbivore diets. However, tall grass swards also represent a more stable damp microclimate for gastrointestinal parasite larvae and tall swards contain generally greater numbers of parasites than short swards. Herbivores may then trade-off the need to maximise nutrient intake through the selection of tall swards with the need to avoid parasite ingestion. Stratified sward sampling techniques were used to determine the distribution of nutrients and parasites in a heterogeneous sward structure on the island of Hirta, St Kilda, which is grazed by a population of feral Soay sheep. The disparity between the short gap vegetation and the tall tussock vegetation was greatest in spring, when gap vegetation was some three-fold shorter than tussock vegetation; this led to tussocks offering greater nutrient and energy intake rate when compared to gap vegetation in spring. Parasites were concentrated in tussocks in spring, thus creating trade-offs. Such trade-offs were not present in summer when parasites were more evenly distributed across the gap/tussock sward structure and the nutritional advantages associated with grazing tussocks were diminished. The diet selection of the resident population of sheep was determined by recording the number of bites taken from gap and tussock vegetation by randomly selected focal animals over repeated 5-min grazing periods. Overall, all animals avoided grazing tussocks, which were most strongly avoided in the spring, and older animals avoided grazing tussocks to a greater extent. Overall, females with lambs and males avoided tussocks to a similar degree and both avoided them less than barren females. Faecal egg counts (FEC) of female sheep were negatively correlated with tussock selectivity and vice versa for male sheep. The interaction between the grazing behaviour of each sex and FEC is discussed in relation to the immunocompetence and nutrient requirements of the different sexes. The maximum disparity between the costs and benefits of the trade-off studied occurred in late winter/early spring which also represents the time of greatest nutrient demand in the Soay sheep, since many are close to starvation and/or are in the advanced stages of pregnancy. Grazing decisions of individuals at this time, determining relative nutrient and parasite intake, may then be related to subsequent fitness and survival.     

Generalist nematodes dominate the nemabiome of roe deer in sympatry with sheep at a regional level

The growth of livestock farming and the recent expansion of wild ungulate populations in Europe favor opportunities for direct and/or indirect cross-transmission of pathogens. Comparatively few studies have investigated the epidemiology of gastro-intestinal nematode parasites, an ubiquitous and important community of parasites of ungulates, at the wildlife/livestock interface. In this study, we aimed to assess the influence of livestock proximity on the gastrointestinal nematode community of roe deer in a rural landscape located in southern France. Using ITS-2 rDNA nemabiome metabarcoding on fecal larvae, we analysed the gastrointestinal nematode communities of roe deer and sheep. In addition, we investigated Haemonchus contortus nad4 mtDNA diversity to specifically test parasite circulation among domestic and wild host populations. The dominant gastrointestinal nematode species found in both the roe deer and sheep were generalist species commonly found in small ruminant livestock (e.g. H. contortus), whereas the more specialised wild cervid nematode species (e.g. Ostertagia leptospicularis) were only present at low frequencies. This is in marked contrast with previous studies that found the nemabiomes of wild cervid populations to be dominated by cervid specialist nematode species. In addition, the lack of genetic structure of the nad4 mtDNA of H. contortus populations between host species suggests circulation of gastrointestinal nematodes between roe deer and sheep. The risk of contact with livestock only has a small influence on the nemabiome of roe deer, suggesting the parasite population of roe deer has been displaced by generalist livestock parasites due to many decades of sheep farming, not only for deer grazing close to pastures, but also at a larger regional scale. We also observed some seasonal variation in the nemabiome composition of roe deer. Overall, our results demonstrate significant exchange of gastrointestinal nematodes between domestic and wild ungulates, with generalist species spilling over from domestic ungulates dominating wild cervid parasite communities.     

The prediction of adaptive evolution: empirical application of the secondary theorem of selection and comparison to the breeder's equation

Adaptive evolution occurs when fitness covaries with genetic merit for a trait (or traits). The breeder's equation (BE), in both its univariate and multivariate forms, allows us to predict this process by combining estimates of selection on phenotype with estimates of genetic (co)variation. However, predictions are only valid if all factors causal for trait-fitness covariance are measured. Although this requirement will rarely (if ever) be met in practice, it can be avoided by applying Robertson's secondary theorem of selection (STS). The STS predicts evolution by directly estimating the genetic basis of trait-fitness covariation without any explicit model of selection. Here we apply the BE and STS to four morphological traits measured in Soay sheep (Ovis aries) from St. Kilda. Despite apparently positive selection on heritable size traits, sheep are not getting larger. However, although the BE predicts increasing size, the STS does not, which is a discrepancy that suggests unmeasured factors are upwardly biasing our estimates of selection on phenotype. We suggest this is likely to be a general issue, and that wider application of the STS could offer at least a partial resolution to the common discrepancy between naive expectations and observed trait dynamics in natural populations.     

A long-term genetic survey of an ungulate population reveals balancing selection acting on MHC through spatial and temporal fluctuations in selection

We explored a 13-year genetic survey of the major histocompatibility complex (MHC) and neutral loci of the Soay sheep population of St Kilda to test the existence and causes of balancing selection at the MHC. The sheep population experiences demographic fluctuations, partly driven by the nematode Teladorsagia circumcincta. The spatial differentiation detected at the MHC was comparable to that at neutral loci between 1988 and 1996, but significantly lower between 1996 and 2000. The rate of temporal genetic differentiation was higher at the MHC, but within the Eastern heft only. These comparisons of spatial and temporal divergence at MHC and non-MHC loci provide strong evidence of balancing selection at the MHC, acting through spatial and temporal heterogeneity in selection pressure. This heterogeneity could be due to fluctuations in the selection imposed by parasites, either directly, because the prevalence in T. circumcincta varies in space and time, or indirectly, because the fitness costs of infection may vary with resource availability.     

Maternal genetic effects set the potential for evolution in a free-living vertebrate population

Heritable maternal effects have important consequences for the evolutionary dynamics of phenotypic traits under selection, but have only rarely been tested for or quantified in evolutionary studies. Here we estimate maternal effects on early-life traits in a feral population of Soay sheep (Ovis aries) from St Kilda, Scotland. We then partition the maternal effects into genetic and environmental components to obtain the first direct estimates of maternal genetic effects in a free-living population, and furthermore test for covariance between direct and maternal genetic effects. Using an animal model approach, direct heritabilities (h2) were low but maternal genetic effects (m2) represented a relatively large proportion of the total phenotypic variance for each trait (birth weight m2=0.119, birth date m2=0.197, natal litter size m2=0.211). A negative correlation between direct and maternal genetic effects was estimated for each trait, but was only statistically significant for natal litter size (ram= -0.714). Total heritabilities (incorporating variance from heritable maternal effects and the direct-maternal genetic covariance) were significant for birth weight and birth date but not for natal litter size. Inadequately specified models greatly overestimated additive genetic variance and hence direct h2 (by a factor of up to 6.45 in the case of birth date). We conclude that failure to model heritable maternal variance can result in over- or under-estimation of the potential for traits to respond to selection, and advocate an increased effort to explicitly measure maternal genetic effects in evolutionary studies.     

Modeling linkage disequilibrium in natural populations: the example of the Soay sheep population of St. Kilda, Scotland

The use of linkage disequilibrium to localize the genes underlying quantitative traits has received considerable attention in the livestock genetics community over the past few years. This has resulted in the investigation of linkage disequilibrium structures of several domestic livestock populations to assess their potential use in fine-mapping efforts. However, the linkage disequilibrium structure of free-living populations has been less well investigated. As the direct evaluation of linkage disequilibrium can be both time consuming and expensive the use of simulations that include as many aspects of population history as possible is advocated as an alternative. A simulation of the linkage disequilibrium structure of the Soay sheep population of St. Kilda, Scotland, is provided as an example. The simulated population showed significant decline of linkage disequilibrium with genetic distance and low levels of background linkage disequilibrium, indicating that the Soay sheep population is a viable resource for linkage disequilibrium fine mapping of quantitative trait loci.     

Population cycles can maintain foraging polymorphism

The maintenance of genetic variability in morphological traits that affect fitness is poorly understood. We present a simple Mendelian model of genetic traits affecting foraging efficiency in grazing ungulates, based on a trade-off between rates of energy extraction at low versus high levels of plant abundance. The model suggests that variation in foraging efficiency could be maintained via lottery competition arising as a direct consequence of dynamically unstable interactions between consumers and their food resources. Lottery competition is a plausible mechanism for explaining wide variability in foraging efficiency, such as that documented in unstable Soay sheep populations on the St Kilda archipelago.     

Factors that affect fertility in a feral population of sheep

Feral livestock offer an excellent opportunity lo study factors affecting fertility as the physiology of their husbanded relatives is well known and social and environmental influences can be studied free of man's interference. This is so in the population of about 1250 sheep of the primitive Soay breed on the islands of St Kilda, Scotland. There is a high coincidence of oestrus amongst the ewes in mid November, and the breeding season is constant within a few days, from year to year. The breeding season in other British breeds of sheep demonstrates a strong correlation with latitude. On St Kilda, mortality rates in male sheep are higher than in females so that the ratio of rams to ewes at mating time is about 1:5. Competition amongst rams is intense. The optimal time for mating and the place of ram lambs in the social system are discussed.     

Modelling non-additive and nonlinear signals from climatic noise in ecological time series: Soay sheep as an example

Understanding how climate can interact with other factors in determining patterns of species abundance is a persistent challenge in ecology. Recent research has suggested that the dynamics exhibited by some populations may be a non-additive function of climate, with climate affecting population growth more strongly at high density than at low density. However, we lack methodologies to adequately explain patterns in population growth generated as a result of interactions between intrinsic factors and extrinsic climatic variation in non-linear systems. We present a novel method (the Functional Coefficient Threshold Auto-Regressive (FCTAR) method) that can identify interacting influences of climate and density on population dynamics from time-series data. We demonstrate its use on count data on the size of the Soay sheep population, which is known to exhibit dynamics generated by nonlinear and non-additive interactions between density and climate, living on Hirta in the St Kilda archipelago. The FCTAR method suggests that climate fluctuations can drive the Soay sheep population between different dynamical regimes--from stable population size through limit cycles and non-periodic fluctuations.     

The use of marker‐based relationship information to estimate the heritability of body weight in a natural population: a cautionary tale

A number of procedures have been developed that allow the genetic parameters of natural populations to be estimated using relationship information inferred from marker data rather than known pedigrees. Three published approaches are available; the regression, pair‐wise likelihood and Markov Chain Monte Carlo (MCMC) sib‐ship reconstruction methods. These were applied to body weight and molecular data collected from the Soay sheep population of St. Kilda, which has a previously determined pedigree. The regression and pair‐wise likelihood approaches do not specify an exact pedigree and yielded unreliable heritability estimates, that were sensitive to alteration of the fixed effects. The MCMC method, which specifies a pedigree prior to heritability estimation, yielded results closer to those determined using the known pedigree. In populations of low average relationship, such as the Soay sheep population, determination of a reliable pedigree is more useful than indirect approaches that do not specify a pedigree.     

Development of a linkage map and mapping of phenotypic polymorphisms in a free-living population of Soay sheep (Ovis aries)

An understanding of the determinants of trait variation and the selective forces acting on it in natural populations would give insights into the process of evolution. The combination of long-term studies of individuals living in the wild and better genomic resources for nonmodel organisms makes achieving this goal feasible. This article reports the development of a complete linkage map in a pedigree of free-living Soay sheep on St. Kilda and its application to mapping the loci responsible for three morphological polymorphisms for which the maintenance of variation demands explanation. The map was derived from 251 microsatellite and four allozyme markers and covers 3350 cM (approximately 90% of the sheep genome) at approximately 15-cM intervals. Marker order was consistent with the published sheep map with the exception of one region on chromosome 1 and one on chromosome 12. Coat color maps to chromosome 2 where a strong candidate gene, tyrosinase-related protein 1 (TYRP1), has also been mapped. Coat pattern maps to chromosome 13, close to the candidate locus Agouti. Horn type maps to chromosome 10, a location similar to that previously identified in domestic sheep. These findings represent an advance in the dissection of the genetic diversity in the wild and provide the foundation for QTL analyses in the study population.     

Natural Selection on Individual Variation in Tolerance of Gastrointestinal Nematode Infection

Hosts may mitigate the impact of parasites by two broad strategies: resistance, which limits parasite burden, and tolerance, which limits the fitness or health cost of increasing parasite burden. The degree and causes of variation in both resistance and tolerance are expected to influence host–parasite evolutionary and epidemiological dynamics and inform disease management, yet very little empirical work has addressed tolerance in wild vertebrates. Here, we applied random regression models to longitudinal data from an unmanaged population of Soay sheep to estimate individual tolerance, defined as the rate of decline in body weight with increasing burden of highly prevalent gastrointestinal nematode parasites. On average, individuals lost weight as parasite burden increased, but whereas some lost weight slowly as burden increased (exhibiting high tolerance), other individuals lost weight significantly more rapidly (exhibiting low tolerance). We then investigated associations between tolerance and fitness using selection gradients that accounted for selection on correlated traits, including body weight. We found evidence for positive phenotypic selection on tolerance: on average, individuals who lost weight more slowly with increasing parasite burden had higher lifetime breeding success. This variation did not have an additive genetic basis. These results reveal that selection on tolerance operates under natural conditions. They also support theoretical predictions for the erosion of additive genetic variance of traits under strong directional selection and fixation of genes conferring tolerance. Our findings provide the first evidence of selection on individual tolerance of infection in animals and suggest practical applications in animal and human disease management in the face of highly prevalent parasites.     

Predictors of early survival in Soay sheep: cohort-, maternal- and individual-level variation

A demographic understanding of population dynamics requires an appreciation of the processes influencing survival-a demographic rate influenced by parameters varying at the individual, maternal and cohort level. There have been few attempts to partition the variance in demography contributed by each of these parameter types. Here, we use data from a feral population of Soay sheep (Ovis aries), from the island of St Kilda, to explore the relative importance of these parameter types on early survival. We demonstrate that the importance of variation occurring at the level of the individual, and maternally, far outweighs that occurring at the cohort level. The most important variables within the individual and maternal levels were birth weight and maternal age class, respectively. This work underlines the importance of using individual based models in ecological demography and we, therefore, caution against studies that focus solely on population processes.