Wild cyclic voles maintain high neutral and MHC diversity without strong evidence for parasite-mediated selection

Neutral Forces and Balancing Selection Interplay to Shape the Major Histocompatibility Complex Spatial Patterns in the Striped Hamster in Inner Mongolia: Suggestive of Broad-Scale Local Adaptation

BACKGROUND

The major histocompatibility complex (MHC) plays a key role in the adaptive immune response to pathogens due to its extraordinary polymorphism. However, the spatial patterns of MHC variation in the striped hamster remain unclear, particularly regarding the relative contribution of the balancing selection in shaping MHC spatial variation and diversity compared to neutral forces.

METHODS

In this study, we investigated the immunogenic variation of the striped hamster in four wild populations in Inner Mongolia which experience a heterogeneous parasitic burden. Our goal was to identify local adaptation by comparing the genetic structure at the MHC with that at seven microsatellite loci, taking into account neutral processes.

RESULTS

We observed significant variation in parasite pressure among sites, with parasite burden showing a correlation with temperature and precipitation. Molecular analysis revealed a similar co-structure between MHC and microsatellite loci. We observed lower genetic differentiation at MHC loci compared to microsatellite loci, and no correlation was found between the two.

CONCLUSIONS

Overall, these results suggest a complex interplay between neutral evolutionary forces and balancing selection in shaping the spatial patterns of MHC variation. Local adaptation was not detected on a small scale but may be applicable on a larger scale.

Living with relatives offsets the harm caused by pathogens in natural populations

Living with relatives can be highly beneficial, enhancing reproduction and survival. High relatedness can, however, increase susceptibility to pathogens. Here, we examine whether the benefits of living with relatives offset the harm caused by pathogens, and if this depends on whether species typically live with kin. Using comparative meta-analysis of plants, animals, and a bacterium (nspecies = 56), we show that high within-group relatedness increases mortality when pathogens are present. In contrast, mortality decreased with relatedness when pathogens were rare, particularly in species that live with kin. Furthermore, across groups variation in mortality was lower when relatedness was high, but abundances of pathogens were more variable. The effects of within-group relatedness were only evident when pathogens were experimentally manipulated, suggesting that the harm caused by pathogens is masked by the benefits of living with relatives in nature. These results highlight the importance of kin selection for understanding disease spread in natural populations.

Macroimmunology: The drivers and consequences of spatial patterns in wildlife immune defence

The prevalence and intensity of parasites in wild hosts varies across space and is a key determinant of infection risk in humans, domestic animals and threatened wildlife. Because the immune system serves as the primary barrier to infection, replication and transmission following exposure, we here consider the environmental drivers of immunity. Spatial variation in parasite pressure, abiotic and biotic conditions, and anthropogenic factors can all shape immunity across spatial scales. Identifying the most important spatial drivers of immunity could help pre-empt infectious disease risks, especially in the context of how large-scale factors such as urbanization affect defence by changing environmental conditions. We provide a synthesis of how to apply macroecological approaches to the study of ecoimmunology (i.e. macroimmunology). We first review spatial factors that could generate spatial variation in defence, highlighting the need for large-scale studies that can differentiate competing environmental predictors of immunity and detailing contexts where this approach might be favoured over small-scale experimental studies. We next conduct a systematic review of the literature to assess the frequency of spatial studies and to classify them according to taxa, immune measures, spatial replication and extent, and statistical methods. We review 210 ecoimmunology studies sampling multiple host populations. We show that whereas spatial approaches are relatively common, spatial replication is generally low and unlikely to provide sufficient environmental variation or power to differentiate competing spatial hypotheses. We also highlight statistical biases in macroimmunology, in that few studies characterize and account for spatial dependence statistically, potentially affecting inferences for the relationships between environmental conditions and immune defence. We use these findings to describe tools from geostatistics and spatial modelling that can improve inference about the associations between environmental and immunological variation. In particular, we emphasize exploratory tools that can guide spatial sampling and highlight the need for greater use of mixed-effects models that account for spatial variability while also allowing researchers to account for both individual- and habitat-level covariates. We finally discuss future research priorities for macroimmunology, including focusing on latitudinal gradients, range expansions and urbanization as being especially amenable to large-scale spatial approaches. Methodologically, we highlight critical opportunities posed by assessing spatial variation in host tolerance, using metagenomics to quantify spatial variation in parasite pressure, coupling large-scale field studies with small-scale field experiments and longitudinal approaches, and applying statistical tools from macroecology and meta-analysis to identify generalizable spatial patterns. Such work will facilitate scaling ecoimmunology from individual- to habitat-level insights about the drivers of immune defence and help predict where environmental change may most alter infectious disease risk.

Pathogen-mediated selection in free-ranging elk populations infected by chronic wasting disease

Pathogens can exert a large influence on the evolution of hosts via selection for alleles or genotypes that moderate pathogen virulence. Inconsistent interactions between parasites and the host genome, such as those resulting from genetic linkages and environmental stochasticity, have largely prevented observation of this process in wildlife species. We examined the prion protein gene (PRNP) in North American elk (Cervus elaphus nelsoni) populations that have been infected with chronic wasting disease (CWD), a contagious, fatal prion disease, and compared allele frequency to populations with no history of exposure to CWD. The PRNP in elk is highly conserved and a single polymorphism at codon 132 can markedly extend CWD latency when the minor leucine allele (132L) is present. We determined population exposure to CWD, genotyped 1,018 elk from five populations, and developed a hierarchical Bayesian model to examine the relationship between CWD prevalence and PRNP 132L allele frequency. Populations infected with CWD for at least 30-50 y exhibited 132L allele frequencies that were on average twice as great (range = 0.23-0.29) as those from uninfected populations (range = 0.04-0.17). Despite numerous differences between the elk populations in this study, the consistency of increase in 132L allele frequency suggests pathogen-mediated selection has occurred due to CWD. Although prior modeling work predicted that selection will continue, the potential for fitness costs of the 132L allele or new prion protein strains to arise suggest that it is prudent to assume balancing selection may prevent fixation of the 132L allele in populations with CWD.

Selection on MHC in a Context of Historical Demographic Change in 2 Closely Distributed Species of Tuco-tucos (Ctenomys australis and C. talarum)

Selection necessarily acts within the same current and historical demographic framework as neutral evolutionary processes, and the outcome of the interplay between these forces may vary according to their relative strength. In this study, we compare the variation at a major histocompatibility complex (MHC) locus (DRB exon 2), typically subject to strong diversifying selection, and mitochondrial diversity (D-loop) across populations encompassing the entire distribution of 2 species of South American subterranean rodents: Ctenomys australis and C. talarum (tuco-tucos). Although these species are parapatric along most of their distribution, historically they have followed distinct demographic trajectories associated with sea level changes during the Quaternary. We surveyed 8 populations of C. australis and 15 of C. talarum, from which we analyzed 70 and 212 D-loop haplotypes and 91 and 346 DRB genotypes, respectively. Both species have gone through a recent demographic expansion; however, the signal of this process only encompasses the entire distribution of one of the species: C. australis. While balancing selection on MHC in C. talarum-enhanced DRB diversity at the local level compared to D-loop, although not promoting divergence among populations, in C. australis local diversifying selection may have driven higher population differentiation at DRB than at D-loop. Our findings reinforce the idea that the relative strength of selection acting on MHC genes varies spatially and temporally within and among species, even between species using the same macrohabitat and exposed to similar immune challenges.

Can MHC-assortative partner choice promote offspring diversity? A new combination of MHC-dependent behaviours among sexes in a highly successful invasive mammal

Sexual selection involving genetically disassortative mate choice is one of several evolutionary processes that can maintain or enhance population genetic variability. Examples of reproductive systems in which choosers (generally females) select mates depending on their major histocompatibility complex (MHC) genes have been reported for several vertebrate species. Notably, the role of MHC-dependent choice not in mating contexts, but in other kinds of social interactions such as in the establishment of complex social systems, has not yet drawn significant scientific interest and is virtually absent from the literature. We have investigated male and female MHC-dependent choice in an invasive population of North American raccoons (Procyon lotor) in Germany. Both male and female raccoons rely on olfaction for individual recognition. Males have an unusually complex social system in which older individuals choose unrelated younger ones to form stable male coalitions that defend territories and a monopoly over females. We have confirmed that females perform MHC-disassortative mate choice and that this behaviour fosters genetic diversity of offspring. We have also observed that males build coalitions by choosing male partners depending on their MHC, but in an assortative manner. This is the first observation of antagonistic MHC-dependent behaviours among sexes. We show that this is the only combination of MHC-dependent partner choice that leads to outbreeding. In the case of introduced raccoons, such behaviours can act together to promote the invasive potential of the species by increasing its adaptive genetic divergence.

Microevolution of bank voles (Myodes glareolus) at neutral and immune-related genes during multiannual dynamic cycles: Consequences for Puumala hantavirus epidemiology

Understanding how host dynamics, including variations of population size and dispersal, may affect the epidemiology of infectious diseases through ecological and evolutionary processes is an active research area. Here we focus on a bank vole (Myodes glareolus) metapopulation surveyed in Finland between 2005 and 2009. Bank vole is the reservoir of Puumala hantavirus (PUUV), the agent of nephropathia epidemica (NE, a mild form of hemorrhagic fever with renal symptom) in humans. M. glareolus populations experience multiannual density fluctuations that may influence the level of genetic diversity maintained in bank voles, PUUV prevalence and NE occurrence. We examine bank vole metapopulation genetics at presumably neutral markers and immune-related genes involved in susceptibility to PUUV (Tnf-promoter, Tlr4, Tlr7 and Mx2 gene) to investigate the links between population dynamics, microevolutionary processes and PUUV epidemiology. We show that genetic drift slightly and transiently affects neutral and adaptive genetic variability within the metapopulation. Gene flow seems to counterbalance its effects during the multiannual density fluctuations. The low abundance phase may therefore be too short to impact genetic variation in the host, and consequently viral genetic diversity. Environmental heterogeneity does not seem to affect vole gene flow, which might explain the absence of spatial structure previously detected in PUUV in this area. Besides, our results suggest the role of vole dispersal on PUUV circulation through sex-specific and density-dependent movements. We find little evidence of selection acting on immune-related genes within this metapopulation. Footprint of positive selection is detected at Tlr-4 gene in 2008 only. We observe marginally significant associations between Mx2 genotype and PUUV genogroups. These results show that neutral processes seem to be the main factors affecting the evolution of these immune-related genes at a contemporary scale, although the relative effects of neutral and adaptive forces could vary temporally with density fluctuations. Immune related gene polymorphism may in turn partly influence PUUV epidemiology in this metapopulation.

Development of novel associations between MHC alleles and susceptibility to parasitic infections in an isolated population of an endangered mammal

The role of pathogens in dynamics of endangered species is not fully understood, and the effect of infection often interacts with other processes affecting those species, such as fragmentation and isolation or loss of genetic variation. Small, isolated populations are prone to losing functional alleles due to demographic processes and genetic drift, which may diminish their ability to resist infection if immune genes are affected. Demographic processes may also alter the selective pressure exerted by a parasite, as they influence the rate of parasite transmission between individuals. In the present paper we studied changes in parasite infection levels and genetic variability in an isolated population of spotted suslik (Spermophillus suslicus). Over a three-year period (approx. three generations), when the population size remained relatively stable, we observed a considerable increase in parasite prevalence and infection intensity, followed by the development of novel associations between MHC DRB alleles and parasite burden. Contrary to expectations, the change in MHC allele frequency over time was not consistent with the effect of the allele - for instance, Spsu-DRB*07, associated with higher intensity of infection with a nematode Capillaria sp., increased in frequency from 11.8 to 20.2%. Yet, we found no signatures of selection in the studied loci. Our results show that an isolated, stable population may experience a sudden increase in parasitic infections, resulting in a development of novel associations between MHC alleles and parasite susceptibility/resistance, even though no signatures of selection can be found.

Evidence for selection maintaining MHC diversity in a rodent species despite strong density fluctuations

Strong spatiotemporal variation in population size often leads to reduced genetic diversity limiting the adaptive potential of individual populations. Key genes of adaptive variation are encoded by the immune genes of the major histocompatibility complex (MHC) playing an essential role in parasite resistance. How MHC variation persists in rodent populations that regularly experience population bottlenecks remains an important topic in evolutionary genetics. We analysed the consequences of strong population fluctuations on MHC class II DRB exon 2 diversity in two distant common vole (Microtus arvalis) populations in three consecutive years using a high-throughput sequencing approach. In 143 individuals, we detected 25 nucleotide alleles translating into 14 unique amino acid MHC alleles belonging to at least three loci. Thus, the overall allelic diversity and amino acid distance among the remaining MHC alleles, used as a surrogate for the range of pathogenic antigens that can be presented to T-cells, are still remarkably high. Both study populations did not show significant population differentiation between years, but significant differences were found between sites. We concluded that selection processes seem to be strong enough to maintain moderate levels of MHC diversity in our study populations outcompeting genetic drift, as the same MHC alleles were conserved between years. Differences in allele frequencies between populations might be the outcome of different local parasite pressures and/or genetic drift. Further understanding of how pathogens vary across space and time will be crucial to further elucidate the mechanisms maintaining MHC diversity in cyclic populations.

Next-generation genotyping of hypervariable loci in many individuals of a non-model species: technical and theoretical implications

BACKGROUND

Across species, diversity at the Major Histocompatibility Complex (MHC) is critical to disease resistance and population health; however, use of MHC diversity to quantify the genetic health of populations has been hampered by the extreme variation found in MHC genes. Next generation sequencing (NGS) technology generates sufficient data to genotype even the most diverse species, but workflows for distinguishing artifacts from alleles are still under development. We used NGS to evaluate the MHC diversity of over 300 captive and wild ring-tailed lemurs (Lemur catta: Primates: Mammalia). We modified a published workflow to address errors that arise from deep sequencing individuals and tested for evidence of selection at the most diverse MHC genes.

RESULTS

In addition to evaluating the accuracy of 454 Titanium and Ion Torrent PGM for genotyping large populations at hypervariable genes, we suggested modifications to improve current methods of allele calling. Using these modifications, we genotyped 302 out of 319 individuals, obtaining an average sequencing depth of over 1000 reads per amplicon. We identified 55 MHC-DRB alleles, 51 of which were previously undescribed, and provide the first sequences of five additional MHC genes: DOA, DOB, DPA, DQA, and DRA. The additional five MHC genes had one or two alleles each with little sequence variation; however, the 55 MHC-DRB alleles showed a high dN/dS ratio and trans-species polymorphism, indicating a history of positive selection. Because each individual possessed 1-7 MHC-DRB alleles, we suggest that ring-tailed lemurs have four, putatively functional, MHC-DRB copies.

CONCLUSIONS

In the future, accurate genotyping methods for NGS data will be critical to assessing genetic variation in non-model species. We recommend that future NGS studies increase the proportion of replicated samples, both within and across platforms, particularly for hypervariable genes like the MHC. Quantifying MHC diversity within non-model species is the first step to assessing the relationship of genetic diversity at functional loci to individual fitness and population viability. Owing to MHC-DRB diversity and copy number, ring-tailed lemurs may serve as an ideal model for estimating the interaction between genetic diversity, fitness, and environment, especially regarding endangered species.

Does the parasite-mediated selection drive the MHC class IIB diversity in wild populations of European chub (Squalius cephalus)?

The genes of major histocompatibility complex (MHC) provide an excellent opportunity to study host-parasite relationships because they are expected to evolve in response to parasites and variation in parasite communities. In this study, we investigated the potential role of parasite-mediated selection acting on MHC class IIB (DAB) genes in European chub (Squalius cephalus) natural populations. We found significant differences between populations in metazoan parasites, neutral and adaptive genetic diversities. The analyses based on pairwise data revealed that populations with dissimilar MHC allelic profiles were geographically distant populations with significantly different diversity in microsatellites and a dissimilar composition of parasite communities. The results from the generalized estimating equations method (GEE) on the level of individuals revealed that metazoan parasite load in European chub was influenced by the diversity of DAB alleles as well as by the diversity of neutral genetic markers and host traits reflecting condition and immunocompetence. The multivariate co-inertia analysis showed specific associations between DAB alleles and parasite species. DAB1-like alleles were more involved in associations with ectoparasites, while DAB3-like alleles were positively associated with endoparasites which could suggest potential differences between DAB genes caused by different selection pressure. Our study revealed that parasite-mediated selection is not the only variable affecting MHC diversity in European chub; however, we strongly support the role of neutral processes as the main driver of DAB diversity across populations. In addition, our study contributes to the understanding of the evolution of MHC genes in wild living fish.

Spatially variable coevolution between a haemosporidian parasite and the MHC of a widely distributed passerine

The environment shapes host-parasite interactions, but how environmental variation affects the diversity and composition of parasite-defense genes of hosts is unresolved. In vertebrates, the highly variable major histocompatibility complex (MHC) gene family plays an essential role in the adaptive immune system by recognizing pathogen infection and initiating the cellular immune response. Investigating MHC-parasite associations across heterogeneous landscapes may elucidate the role of spatially fluctuating selection in the maintenance of high levels of genetic variation at the MHC. We studied patterns of association between an avian haemosporidian blood parasite and the MHC of rufous-collared sparrows (Zonotrichia capensis) that inhabit environments with widely varying haemosporidian infection prevalence in the Peruvian Andes. MHC diversity peaked in populations with high infection prevalence, although intra-individual MHC diversity was not associated with infection status. MHC nucleotide and protein sequences associated with infection absence tended to be rare, consistent with negative frequency-dependent selection. We found an MHC variant associated with a ∽26% decrease in infection probability at middle elevations (1501-3100 m) where prevalence was highest. Several other variants were associated with a significant increase in infection probability in low haemosporidian prevalence environments, which can be interpreted as susceptibility or quantitative resistance. Our study highlights important challenges in understanding MHC evolution in natural systems, but may point to a role of negative frequency-dependent selection and fluctuating spatial selection in the evolution of Z. capensisMHC.