The role of parasite-driven selection in shaping landscape genomic structure in red grouse (Lagopus lagopus scotica).

Landscape genomics promises to provide novel insights into how neutral and adaptive processes shape genome-wide variation within and among populations. However, there has been little emphasis on examining whether individual-based phenotype-genotype relationships derived from approaches such as genome-wide association (GWAS) manifest themselves as a population-level signature of selection in a landscape context. The two may prove irreconcilable as individual-level patterns become diluted by high levels of gene flow and complex phenotypic or environmental heterogeneity. We illustrate this issue with a case study that examines the role of the highly prevalent gastrointestinal nematode Trichostrongylus tenuis in shaping genomic signatures of selection in red grouse (Lagopus lagopus scotica). Individual-level GWAS involving 384 SNPs has previously identified five SNPs that explain variation in T. tenuis burden. Here, we examine whether these same SNPs display population-level relationships between T. tenuis burden and genetic structure across a small-scale landscape of 21 sites with heterogeneous parasite pressure. Moreover, we identify adaptive SNPs showing signatures of directional selection using F(ST) outlier analysis and relate population- and individual-level patterns of multilocus neutral and adaptive genetic structure to T. tenuis burden. The five candidate SNPs for parasite-driven selection were neither associated with T. tenuis burden on a population level, nor under directional selection. Similarly, there was no evidence of parasite-driven selection in SNPs identified as candidates for directional selection. We discuss these results in the context of red grouse ecology and highlight the broader consequences for the utility of landscape genomics approaches for identifying signatures of selection.

Genome-wide association and genome partitioning reveal novel genomic regions underlying variation in gastrointestinal nematode burden in a wild bird.

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).

Host migration impacts on the phylogeography of Lyme Borreliosis spirochaete species in Europe.

The geographic patterns of transmission opportunities of vector-borne zoonoses are determined by a complex interplay between the migration patterns of the host and the vector. Here we examine the impact of host migration on the spread of a tick-borne zoonotic disease, using Lyme Borreliosis (LB) spirochaetal species in Europe. We demonstrate that the migration of the LB species is dependent on and limited by the migration of their respective hosts. We note that populations of Borrelia spp. associated with birds (Borrelia garinii and B. valaisiana) show limited geographic structuring between countries compared with those associated with small mammals (Borrelia afzelii), and we argue that this can be explained by higher rates of migration in avian hosts. We also show the presence of B. afzelii strains in England and, through the use of the multi-locus sequence analysis scheme, reveal that the strains are highly structured. This pattern in English sites is very different from that observed at the continental sites, and we propose that these may be recent introductions.

The Heterogeneity, Distribution, and Environmental Associations of Borrelia burgdorferi Sensu Lato, the Agent of Lyme Borreliosis, in Scotland.

Lyme borreliosis is an emerging infectious human disease caused by the Borrelia burgdorferi sensu lato complex of bacteria with reported cases increasing in many areas of Europe and North America. To understand the drivers of disease risk and the distribution of symptoms, which may improve mitigation and diagnostics, here we characterize the genetics, distribution, and environmental associations of B. burgdorferi s.l. genospecies across Scotland. In Scotland, reported Lyme borreliosis cases have increased almost 10-fold since 2000 but the distribution of B. burgdorferi s.l. is so far unstudied. Using a large survey of over 2200 Ixodes ricinus tick samples collected from birds, mammals, and vegetation across 25 sites we identified four genospecies: Borrelia afzelii (48%), Borrelia garinii (36%), Borrelia valaisiana (8%), and B. burgdorferi sensu stricto (7%), and one mixed genospecies infection. Surprisingly, 90% of the sequence types were novel and, importantly, up to 14% of samples were mixed intra-genospecies co-infections, suggesting tick co-feeding, feeding on multiple hosts, or multiple infections in hosts. B. garinii (hosted by birds) was considerably more genetically diverse than B. afzelii (hosted by small mammals), as predicted since there are more species of birds than small mammals and birds can import strains from mainland Europe. Higher proportions of samples contained B. garinii and B. valaisiana in the west, while B. afzelii and B. garinii were significantly more associated with mixed/deciduous than with coniferous woodlands. This may relate to the abundance of transmission hosts in different regions and habitats. These data on the genetic heterogeneity within and between Borrelia genospecies are a first step to understand pathogen spread and could help explain the distribution of patient symptoms, which may aid local diagnosis. Understanding the environmental associations of the pathogens is critical for rational policy making for disease risk mitigation and land management.