Epigenetic modification of cytosine methylation states can be elicited by environmental stresses and may be a key process affecting phenotypic plasticity and adaptation. Parasites are potent stressors with profound physiological and ecological effects on their host, but there is little understanding in how parasites may influence host methylation states. Here, we estimate epigenetic diversity and differentiation among 21 populations of red grouse (Lagopus lagopus scotica) in north-east Scotland and test for association of gastrointestinal parasite load (caecal nematode Trichostrongylus tenuis) with hepatic genome-wide and locus-specific methylation states. Following methylation-sensitive AFLP (MSAP), 129 bands, representing 73 methylation-susceptible and 56 nonmethylated epiloci, were scored across 234 individuals. The populations differed significantly in genome-wide methylation levels and were also significantly epigenetically (FSC = 0.0227; P < 0.001) and genetically (FSC = 0.0058; P < 0.001) differentiated. Parasite load was not associated with either genome-wide methylation levels or epigenetic differentiation. Instead, we found eight disproportionately differentiated epilocus-specific methylation states (FST outliers) using bayescan software and significant positive and negative association of 35 methylation states with parasite load from bespoke generalized estimating equations (GEE), simple logistic regression (sam) and Bayesian environmental analysis (bayenv2). Following Sanger sequencing, genome mapping and geneontology (go) annotation, some of these epiloci were linked to genes involved in regulation of cell cycle, signalling, metabolism, immune system and notably rRNA methylation, histone acetylation and small RNAs. These findings demonstrate an epigenetic signature of parasite load in populations of a wild bird and suggest intriguing physiological effects of parasite-associated cytosine methylation.

Selective pressure from pathogens is considered a key selective force driving the evolution of components of the immune system. Since single components of the immune system may interact with many pathogens, and single pathogens may be recognized by multiple components of the immune system, gaining a better understanding of the mechanisms of parasite-driven selection requires the study of multiple genes and pathogens. Toll-like receptors (TLRs) are a large gene family that code for antigen-presenting components of the innate immune response. In the present paper we characterize polymorphism and signatures of selection in seven TLRs in free-living bank voles Myodes glareolus. We report the first evidence of balancing selection in several TLR genes, supported by positive values of Fu and Li’s D* in TLR2 and TLR5, and positive values of Tajima’s D in LRR regions within TLR1 and TLR2. We further found significant associations between amino-acid alleles of TLR1 and TLR5 and susceptibility to infection with the blood pathogen Bartonella. Interestingly, selection patterns in TLRs presenting virus-derived motifs (TLR7 and TLR9) differed considerably from those interacting with bacterial PAMPs. In contrast to the highly variable TLRs presenting bacterial motifs, TLR7 and TLR9 had low polymorphism and displayed signatures of directional selection. These findings suggest different functional responses across the TLR gene family and highlight the complexity of parasite-driven selection.

Identifying the genetic architecture underlying complex phenotypes is a notoriously difficult problem that often impedes progress in understanding adaptive eco-evolutionary processes in natural populations. Host-parasite interactions are fundamentally important drivers of evolutionary processes, but a lack of understanding of the genes involved in the host's response to chronic parasite insult makes it particularly difficult to understand the mechanisms of host life history trade-offs and the adaptive dynamics involved. Here, we examine the genetic basis of gastrointestinal nematode (Trichostrongylus tenuis) burden in 695 red grouse (Lagopus lagopus scotica) individuals genotyped at 384 genome-wide SNPs. We first use genome-wide association to identify individual SNPs associated with nematode burden. We then partition genome-wide heritability to identify chromosomes with greater heritability than expected from gene content, due to harbouring a multitude of additive SNPs with individually undetectable effects. We identified five SNPs on five chromosomes that accounted for differences of up to 556 worms per bird, but together explained at best 4.9% of the phenotypic variance. These SNPs were closely linked to genes representing a range of physiological processes including the immune system, protein degradation and energy metabolism. Genome partitioning indicated genome-wide heritability of up to 29% and three chromosomes with excess heritability of up to 4.3% (total 8.9%). These results implicate SNPs and novel genomic regions underlying nematode burden in this system and suggest that this phenotype is somewhere between being based on few large-effect genes (oligogenic) and based on a large number of genes with small individual but large combined effects (polygenic).

The extent to which genotypic variation at a priori identified candidate genes can explain variation in complex phenotypes is a major debate in evolutionary biology. Whereas some high-profile genes such as the MHC or MC1R clearly do account for variation in ecologically relevant characters, many complex phenotypes such as response to parasite infection may well be underpinned by a large number of genes, each of small and effectively undetectable effect. Here, we characterize a suite of novel candidate genes for variation in gastrointestinal nematode (Trichostrongylus tenuis) burden among red grouse (Lagopus lagopus scotica) individuals across a network of moors in north-east Scotland. We test for associations between parasite load and genotypic variation in twelve genes previously identified to be differentially expressed in experimentally infected red grouse or genetically differentiated among red grouse populations with overall different parasite loads. These genes are associated with a broad physiological response including immune system processes. Based on individual-level generalized linear models, genotypic variants in nine genes were significantly associated with parasite load, with effect sizes accounting for differences of 514-666 worms per bird. All but one of these variants were synonymous or untranslated, suggesting that these may be linked to protein-coding variants or affect regulatory processes. In contrast, population-level analyses revealed few and inconsistent associations with parasite load, and little evidence of signatures of natural selection. We discuss the broader significance of these contrasting results in the context of the utility of population genomics and landscape genomics approaches in detecting adaptive genomic signatures.

The red-billed chough (Pyrrhocorax pyrrhocorax) is of conservation concern in the British Isles and continental Europe, with historically declining populations and a highly fragmented distribution. We quantified the distribution of genetic variation within and among European populations to identify isolated populations that may need to be managed as demographically independent units, and assess whether individual populations are denuded of genetic diversity and so may show reduced viability. We genotyped 326 choughs from ten wild populations and 22 from one captive population at 16 nuclear microsatellite loci, and sequenced 34 individuals across three mitochondrial regions to quantify genetic structure, diversity and phylogeography. Microsatellite diversity was low (often \4 alleles per locus), but pairwise population differentiation was high (often D est [ 0.1), with a signature of isolationby-distance. Bayesian-inferred a posteriori genetic clusters coincided with a priori populations, supporting strong genetic structure. Microsatellites also allowed us to identify the probable origin of the captive choughs and one recently founded wild population. Mitochondrial DNA sequence diversity was low (p = 0.00103). Phylogeographic structure was consequently poorly resolved, but indicated that sampled continentalEuropean populations are ancestral to British Isles populations, which comprised a single clade. Our data suggest that British Isles chough populations are relatively isolated with infrequent gene flow and relatively genetically depauperate, potentially requiring genetic management. These findings should be integrated into conservation management policy to ensure longterm viability of chough populations.